# Longitudinal Survey of Coronavirus Circulation and Diversity in Insectivorous Bat Colonies in Zimbabwe

Vimbiso Chidoti, Hélène de Nys, Valérie Pinarello, Getrude Mashura, Dorothée Missé, Laure Guerrini, Davies Pfukenyi, Julien Cappelle, Ngoni Chiweshe, Ahidjo Ayouba, et al.

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

## To cite this version:

Vimbiso Chidoti, Hélène de Nys, Valérie Pinarello, Getrude Mashura, Dorothée Missé, et al. Longitudinal Survey of Coronavirus Circulation and Diversity in Insectivorous Bat Colonies in Zimbabwe. Viruses, 2022, 14 (4), pp.781. 10.3390/v14040781. hal-03663465

HAL Id: hal-03663465

https://hal.inrae.fr/hal-03663465

Submitted on 10 May 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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
Article

# Longitudinal Survey of Coronavirus Circulation and Diversity in Insectivorous Bat Colonies in Zimbabwe

Vimbiso Chidoti<fnref n="1" />, Hélène De Nys 2,3 ID, Valérie Pinarello 1,2,3, Getrude Mashura 1, Dorothée Missé 4 ID,
Laure Guerrini 2,3, Davies Pfukenyi 1, Julien Cappelle 2,5 ID, Ngoni Chiweshe 3, Ahidjo Ayouba 6 ID, Gift Matope 1,
Martine Peeters 6, Elizabeth Gori 1 ID, Mathieu Bourgarel 2,3 ID and Florian Liégeois 1,4,* ID

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

check for
updates

**Citation:** Chidoti, V.; De Nys, H.; Pinarello, V.; Mashura, G.; Missé, D.; Guerrini, L.; Pfukenyi, D.; Cappelle, J.; Chiweshe, N.; Ayouba, A.; et al. Longitudinal Survey of Coronavirus Circulation and Diversity in Insectivorous Bat Colonies in Zimbabwe. *Viruses* **2022**, *14*, 781. https://doi.org/10.3390/v14040781

Academic Editor: Gilda Tachedjian

Received: 16 February 2022

Accepted: 3 April 2022

Published: 9 April 2022

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

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

Copyright: © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** **Background:** Studies have linked bats to outbreaks of viral diseases in human populations such as SARS-CoV-1 and MERS-CoV and the ongoing SARS-CoV-2 pandemic. Methods: We carried out a longitudinal survey from August 2020 to July 2021 at two sites in Zimbabwe with bat-human interactions: Magweto cave and Chirundu farm. A total of 1732 and 1866 individual bat fecal samples were collected, respectively. Coronaviruses and bat species were amplified using PCR systems. Results: Analysis of the coronavirus sequences revealed a high genetic diversity, and we identified different sub-viral groups in the *Alphacoronavirus* and *Betacoronavirus* genus. The established sub-viral groups fell within the described *Alphacoronavirus* sub-genera: *Decacovirus*, *Duvinacovirus*, *Rhinacovirus*, *Setracovirus* and *Minunacovirus* and for *Betacoronavirus* sub-genera: *Sarbecoviruses*, *Merbecovirus* and *Hibecovirus*. Our results showed an overall proportion for CoV positive PCR tests of 23.7% at Chirundu site and 16.5% and 38.9% at Magweto site for insectivorous bats and *Macronycteris gigas*, respectively. Conclusions: The higher risk of bat coronavirus exposure for humans was found in December to March in relation to higher viral shedding peaks of coronaviruses in the parturition, lactation and weaning months of the bat populations at both sites. We also highlight the need to further document viral infectious risk in human/domestic animal populations surrounding bat habitats in Zimbabwe.

**Keywords:** bat coronavirus (Bt CoVs); human-bat interaction; genetic diversity; reproductive phenology; Zimbabwe

1. Introduction

The *Coronaviridae* family is a monophyletic cluster in the order of *Nidovirales*. These are enveloped positive stranded RNA viruses of three classes of vertebrates: mammals (corona- and toroviruses), birds (coronaviruses) and fish (bafiniviruses) [1]. They are classified into four genera: *Alphacoronavirus* and *Betacoronavirus*, hosted by mammals; *Gammacoronavirus* and *Deltacoronavirus*, infecting avian species [2,3].

Coronaviruses are etiological agents of respiratory, enteric, hepatic and neurological diseases in animals and humans [4], with variable severity, from asymptomatic to severely ill individuals [5]. Surveillance of coronaviruses in wild animals has led to the discovery of a high diversity of viruses in bats and avian species, suggesting these animals as the natural reservoirs [6].

<footnote n="1">¹ Faculty of Veterinary Science, University of Zimbabwe, Harare P.O. Box MP 167, Zimbabwe; vimbisochidoti@gmail.com (V.C.); valerie.pinarello@cird.fr (V.P.); getrudemashura@gmail.com (G.M.); dmpfukenyi@vet.uz.ac.zw (D.P.); giftmatope@gmail.com (G.M.); gori.elizabeth@gmail.com (E.G.)</footnote>

<footnote n="2">² ASTRE, CIRAD, INRAE, University of Montpellier, 34980 Montpellier, France; helene.de_nys@cirad.fr (H.D.N.); laure.guerrini@cirad.fr (L.G.); julien.cappelle@cirad.fr (J.C.); mathieu.bourgarel@cirad.fr (M.B.)</footnote>

<footnote n="3">³ CIRAD, UMR ASTRE, Harare, Zimbabwe; chiweshengoni@gmail.com</footnote>

<footnote n="4">⁴ MIVEGEC, University of Montpellier, IRD, CNRS, 34394 Montpellier, France; dorothee.misse@ird.fr</footnote>

<footnote n="5">5 CIRAD, UMR ASTRE, 34398 Montpellier, France</footnote>

<footnote n="6">⁶ TransVIHMI, University of Montpellier, IRD, Inserm, 34394 Montpellier, France; ahidjo.ayouba@ird.fr (A.A.); martine.peeters@ird.fr (M.P.)</footnote>

<footnote>* Correspondence: florian.liegeois@ird.fr</footnote>
Seven coronaviruses have been identified in humans [3]. The first human coronaviruses were discovered in the early 1960s, nCoV-OC43 and HCoV-229E [7,8], and they are known to cause mild respiratory disease [9]. HCoV-NL63 and -HKU1 were discovered in 2004 and 2005, respectively [10], they also associated with mild respiratory disease [9]. Other high profile human disease outbreaks of the respiratory system caused by coronaviruses took place during the last two decades, namely severe acute respiratory syndrome type 1 and 2 (SARS-CoV-1 and SARS-CoV-2) and Middle East respiratory syndrome (MERS-CoV) [11–13].

Phylogenetic evidence shows that HCoV-NL63 and -229E originated from bat coronaviruses [7,8], whereas HCoV-OC43 and HCoV-HKU1 originated from rodents [8,14]. Bats are major and ancient reservoirs of several viral families [15], and recent molecular studies demonstrated bats as natural hosts of important zoonotic viruses such as *Filoviruses* and *Paramyxoviruses* among others [16].

Insectivorous bats have recently been shown to harbor the precursor of severe acute respiratory syndrome coronavirus (SARS-CoV-1) [17]. Studies have linked bats (*Rhinolophus* spp.) to SARS-CoV-1 outbreak in humans, transmitted through an intermediate host, i.e., civet cats [18]. Recent research on the outbreak of SARS-CoV-2 (COVID-19) suggest that bats are the original source of the coronavirus due to genetic similarity of the virus to other known bat coronaviruses [12,19,20]. Bat coronaviruses have also been identified in numerous insectivorous and frugivorous bats in Asia, Africa, America and Europe [21,22].

Bats are the second most diverse mammalian order after rodents [23], and the only order of mammals capable of flying. Bats constitute the order Chiroptera [23] and are classified into two groups: *Yinpterochiroptera*, megabats and a number of microbats, and *Yangochiroptera*, which mainly covers insectivorous microbats [24,25]. Bats possess characteristics that maximize their effectiveness as reservoir hosts. These include high species diversity, long life span, long flight distance and dispersal, interspecies dense roost aggregations, social behaviors, hibernation and torpor as well as unique innate immunology [26,27].

The circulation of viruses in bat populations is affected by seasonality and the prevailing environmental conditions [23,28]. Each of these directly affects reproduction periods, gregarious behavior, torpor, population density and rate of contacts between individuals in a population [27,28] consequently determining reproductive number of viruses and their transmission [28]. Hence, risks of viral spillover to other animal species and humans are strongly influenced by these factors. The effect of seasonality on viral circulation has been shown in several studies on coronaviruses in bats [29,30]. Temporal variations associated with the bat reproduction phenology have been shown in several coronavirus studies across the world in bats, with higher coronavirus prevalence detected in juvenile and sexually immature individuals [31,32].

Previous studies carried out in Zimbabwe included the characterization of *Dicistrovirus*, *Coronavirus* and *Paramyxovirus* in insectivorous bats (*Hipposideros caffer*) [33,34]. More than 60 bat species have been recorded in Zimbabwe and in rural areas, and people are in close contact with these animals [35], posing a risk of spillover of bat-associated viruses to humans. Temporal variations associated with the bat reproduction phenology have been shown in several coronavirus studies across the world in bats, with higher coronavirus prevalence detected in juvenile and sexually immature individuals [31,32]. The prevalence and genetic diversity of coronaviruses circulating in insectivorous bat colonies and their seasonal pattern are not well established. This study sought to establish the prevalence and the genetic diversity of coronaviruses in insectivorous bats according to bat reproduction phenology in two multi-species bat sites in Zimbabwe.

**2. Methods**

**2.1. Study Sites**

The study was conducted in Hurungwe district of Mashonaland West province in Zimbabwe, at two sites, Magweto cave and Chirundu farm, known for their colonies of
insectivorous bats. The two sites situated in communal lands represent habitats with close human-bat interactions (Figure 1).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** The study sites highlighted on the map show the sampling areas for insectivorous bat colonies in Hurungwe district.</figcaption></figure>

## 2.2. *Magweto Cave*

Magweto cave site has been subject to preliminary research on virus circulation in *Hipposideros caffer* bats [33]. The site is considered to be a sacred place by the local community. The cave itself is used for various activities such as cultural and religious ceremonies. People also collect guano (bat feces) and use it as fertilizer in their fields, while others sell it. Additionally, *Macronycteris gigas* bats are also present in this cave and are hunted for consumption.

## 2.3. *Chirundu Farm*

Chirundu farm site is surrounded by a banana plantation where a large group of people work daily. The building basement in which the bats roost is surrounded by a nursery of young trees. Due to irrigation of the plants, the tunnels are always flooded with water most part of the year. The bats fly around the property and feed on insects around the plantation.

## 2.4. *Sample and Data Collection*

On both sites, non-invasive sampling of bat feces was performed from August 2020 to July 2021 at monthly intervals. Two-square-meter plastic sheets were placed under the bat roosts overnight (5 sheets per site). Feces were collected in the mid-morning from each plastic sheet. Individual sampling was carried out by assuming that a space of 20 by 20 cm on the plastic sheet was representative of the defecation area of one individual: One dropping of fresh feces was selected and placed in 2 mL tubes with 0.5 mL of home-made RNA stabilization solution (http://www.protocol-online.org/prot/Protocols/RNAlater-3999.html; accessed on 1 May 2020). Samples were then stored at -80 °C before further laboratory analyses. During sampling, other observations were recorded such as females showing pregnancy, puppies attached to mothers and puppies in specific roost nurseries. This information was used to map a timeline of bat reproductive phenology. Periods of gestation, parturition, lactation, weaning and presence of 4-6 month old juveniles were determined based on observations (from captures and observation of the roosting bats) combined with the literature to confirm observations or to complete observation gaps (regarding periods or species) [35].
Additionally, two capture sessions were performed at each site In October 2020 and March 2021 to assess the bat diversity and reproductive status. Bats were captured using nest and hard traps. All animals were released after collection of dried blood spot, fecal and salivary swab samples as well as morphometric measurements [36].

## 2.5. RNA/DNA Extraction

Both RNA and DNA were extracted from 200 µL of fecal samples using 5X MagMax Pathogen RNA/DNA Kit (Thermo Fisher Scientific, Illkich, France) by using the fecal samples preserved in 0.5 mL of home-made preservative RNA solution. The feces were vortexed vigorously (30 Hz) using a Retsch MM400 Tissue lyser for 5 min to fully homogenize and mix well the fecal particles, followed by centrifugation at 16,000× g for ~3 min to fully separate the supernatant from the fecal pellet. A total of 130 µL of the supernatant was collected to isolate and purify the nucleic acids using a Mag Max extraction kit (Thermo Fisher Scientific, France) with the automatic KingFisher Duo Prime Purification System (Thermo Fisher Scientific, France) following the manufacturer's instructions. A final volume of 80 µL of eluted RNA/DNA was collected and stored at -80 °C.

## 2.6. Reverse Transcription (RT) Using Random Hexamers

A total of 5 µL of RNA/DNA sample template was reverse transcribed using 1 µL random hexamers, 0.5 µL Oligo dT primer, 0.4 µL of dNTPs (10 mM) (Thermo Fisher Scientific, France) and 5.5 µL molecular grade water and incubated for 5 min at 65 °C. Next, 4 µL of Buffer 5×, 2 µL of 0.1 M DTT (M-MLV Reverse Transcriptase, Thermo Fisher Scientific, France) and 1 µL of RNAse OUT were added to the previous mix. After incubation for 2 min at 37 °C, 1 µL of M-MLV reverse transcriptase (M-MLV Reverse Transcriptase, Invitrogen, Thermo Fisher Scientific, Illkrich, France) was added to the mixture followed by incubation at 25 °C for 10 min, 37 °C for 50 min and 70 °C for 15 min. The cDNA obtained was stored at -20 °C.

## 2.7. RNA-Dependent RNA-Polymerase (RdRp) Nested PCR on Coronaviruses Using Pan-Coronavirus Primers

The CoV partial *RdPd* gene was amplified by nested PCR. Both PCR rounds were realized in 50 µL reaction mixture containing 5 µL of cDNA template (first round) and 5 µL of first round amplified product (second round), 0.5 µL DNA polymerase, 5 µL Buffer 10X, 4 µL MgCl₂ (25 mM) (Firepol DNA Polymerase, Solis Bio Dyne, Tartu, Estonia), 0.4 µL dNTPs (25 mM) (Thermo Fisher Scientific, Illkrich, France), 31.1 µL of molecular grade water and 2 µL of pan-coronavirus primers (10 mM) (Thermo Fisher Scientific, Illkrich, France): Pan-CoV-F1: 5'-GGKTGGGAYTAYCCKAARTG-3' and Pan-CoV-R1: 3'TGYTGTSWRCARAAYTCRTG-5' for the first PCR round, Pan-CoVF2: 5'GGTTGGGACTATCCTAAGTGTGA-3', Pan-CoVR2: 3'-CCATCATCAGATAGAATCATCAT-5' for the second PCR round as previously described [37]. First round PCR conditions were as follows: 95 °C, 2 min denaturation followed by 40 cycles at 95 °C, 20 s denaturation, 50 °C, 30 s hybridization, 72 °C, 2 min elongation and a final extension of 72 °C for 5 min. Second round PCR conditions were as follows: 95 °C, 2 min denaturation, 40 cycles at 95 °C, 30 s denaturation, 52 °C, 30 s hybridization, 72 °C, 1 min elongation and a final extension of 72 °C for 5 min. For all PCRs, positive and negative controls were run in parallel. Visualization of positive PCR product was carried out with gel electrophoresis using ethidium bromide on a 1% gel. PCR positive products (440 bp) were directly purified (GeneClean Turbo Kit, MP Biomedicals, Illkirch, France) according to supplier's instructions and then sequenced using the Sanger sequencing method in both the 5' and 3' directions (LGC, Berlin, Germany).

## 2.8. Genetic Analyses

### 2.8.1. Viruses Sequence Identification and Phylogenetic Analysis

Generated Sanger viral sequence assemblage was performed using Geneious software package V. 2021.2.2 (Biomatters Ltd., Auckland, New Zealand). The generated contigs
were compared to a database of sequences online using Basic Local Alignment Search Tool (BLAST) (https://blast.ncbi.nlm.nih.gov/Blast.cgi) to identify the amplified sequence. Mega 7 [38] was used to align the CoV sequences obtained in this study with references obtained from GenBank. A maximum likelihood phylogenetic tree of aligned CoV sequences was implemented by using IQ-Tree [39]. The reliability of branching orders was tested using the bootstrap approach (1000 replicates), and the GTR + F+ I substitution model was determined as the best suited evolution model [39].

## 2.8.2. Bat Species Identification

Genotyping was carried out only on *Coronavirus* positive samples for bat species identification. Partial amplification of cytochrome B and 12S RNA mitochondrial genes followed by sequencing was used for identification of bat species or genus, respectively, for all the CoV positive samples from the two sites. Cytochrome B/12S RNA partial sequences were amplified by one step PCR in 50 µL reaction mix containing 5 µL of RNA template, 0.5 µL DreamTaq DNA polymerase, 5 µL DreamTaq Buffer 10× (DreamTaq DNA Polymerase, Thermo Fisher Scientific, US), 0.4 µL dNTPs (25 mM) (Thermo Fisher Scientific), 35.35 µL of molecular grade water and 2 µL of primers (10 mM) (Thermo Fisher Scientific, France): **Cytb-L14724**: 5'-CGAAGCTTGATATGAAAAACCATCGTTG-3', **Cytb-H15506**: 5'-AGTGGRTTRGCTGGTGTRTARTTGTC-3' [40]. For the **12S RNA, L1091**: (5'-AAAAAG-CTTCAAACTGGGATTAGATACCCACTAT-3') and **H1478**: (5'-TGACTGCAGAGGGTGACGGGCGGTGTGT-3') primers were used [41]. The PCR reaction was carried out under the following conditions: 95 °C, 3 min denaturation, 20 cycles at 94 °C, 20 s denaturation, 45 °C, 30 s hybridization, 72 °C, 1 min 30 s, followed by 20 cycles at 94 °C, 20 s denaturation, 50 °C, 30 s hybridization, 72 °C, 1 min 30 s elongation and a final extension of 72 °C for 5 min. Visualization of positive PCR product was performed with gel electrophoresis using ethidium bromide on a 1% gel. PCR positive products for cytochrome B (800 bp) and for 12S RNA (400 bp) were directly purified (GeneClean Turbo Kit, MP Biomedicals, Illkirch, France) according to supplier's instructions and then sequenced using Sanger sequencing method in both the 5' and 3' directions (LGC, Germany). Cytochrome B and 12S RNA sequences obtained were then compared to bat species sequences available in the GenBank database using BLAST tool (https://blast.ncbi.nlm.nih.gov/Blast.cgi).

## 2.8.3. Temporal Variations of Coronavirus Prevalence and Bat Reproduction Phenology

The prevalence was calculated at the community level for insectivorous bat species from the two sites. The *M. gigas* prevalence was calculated separately as fecal samples could easily be attributed to this species given their size and location in the cave. The reproductive phases of the bats were allocated according to field observations and also the literature review. The proportion of RNA CoV positive individuals was calculated by month and site with 95% confidence intervals (CIs) using the Wilson score test [42]. Descriptive analysis of prevalence of CoVs per month/per season at each site was performed using R software (R-studio) programming to construct box plots. The boxplots were plotted with the x-axis as the prevalence of CoVs and y-axis as the date and reproductive season when sampling occurred at each site. Similar analysis was also carried out to describe the effect of reproductive seasonality on the detection of the sub-viral groups described. Sub-viral groups with more than 30 sequences were analyzed; only Chirundu site analyses were performed as at Magweto there were insufficient data to allow an accurate conclusion to be drawn.

In addition, the influence of the different phases of the bat reproduction cycle (pregnancy, parturition/lactation, weaning, weaned juveniles of 4 to 6 months old) on the prevalence of coronaviruses at each site (Magweto and Chirundu) was tested by running two generalized linear mixed models (GLMMs). Parturition and lactation were merged because our sampling sessions did not allow the separation of the two periods. The response variable with a binomial distribution was the coronavirus positivity of the samples, and the explanatory variables with a fixed effect were the different phases of the reproduction
cycle (coded as 1 if it was during the corresponding reproduction phase and 0 if it was not). To account for clustered samples collected during the same sampling session at the same site, we included the session identification code as a random effect to control for repeated measures from the same trapping session.

### 2.8.4. GenBank Accession Numbers

The Coronavirus sequences have been deposited to GenBank under the following numbers: OM469940-OM470469.

The Cytochrome B sequences have been deposited to GenBank under the following numbers: OM487705-OM488020.

## 3. Results

### 3.1. *Sampling, Data Collection and Reproductive Phenology*

Due to travel restrictions linked to COVID-19 in Zimbabwe, no samples were collected in January, April or June 2021 at Chirundu farm and in Magweto none were collected from December 2020 to February 2021 or in May and July 2021. We analyzed 1866 fecal samples from different insectivorous bat species including 308 fecal samples from Macronycteris gigas species collected at Magweto cave site and a total of 1732 fecal samples from insectivorous bats collected at the Chirundu farm site (Table 1).

<table><caption>**Table 1.** Prevalence of CoVs and confidence intervals (CIs) per month at both sites in fecal samples from insectivorous bat communities.</caption><thead><tr><th>Site</th><th>Reproduction Cycle</th><th>Month Sampled</th><th>No of Samples Tested</th><th>No of Coronavirus Positives</th><th>Prevalence (%) + CI (95%)</th></tr></thead><tbody><tr><td rowspan="8">Chirundu</td><td>Non-gestation</td><td>August 2020</td><td>154</td><td>9</td><td>5.8 (3.1–10.7)</td></tr><tr><td>Pregnancy</td><td>October 2020</td><td>296</td><td>4</td><td>1.35 (0.5–3.4)</td></tr><tr><td>Parturition</td><td>November 2020</td><td>242</td><td>45</td><td>18 (14.2–23.9)</td></tr><tr><td>Lactation</td><td>December 2020</td><td>160</td><td>60</td><td>37 (30.4–45.2)</td></tr><tr><td>Weaning</td><td>Februar 2021</td><td>172</td><td>76</td><td>44.2 (37–51.7)</td></tr><tr><td>4–6 months juveniles</td><td>March 2021</td><td>240</td><td>84</td><td>35 (29.2–41.2)</td></tr><tr><td>4–6 months juveniles</td><td>May 2021</td><td>242</td><td>61</td><td>25.2 (20.2–31)</td></tr><tr><td>Non-gestation</td><td>July 2021</td><td>226</td><td>71</td><td>31 (25.7–37.7)</td></tr><tr><td colspan="2">Overall prevalence</td><td></td><td>1732</td><td>410</td><td>23.7 (21.73–25.73)</td></tr><tr><td rowspan="6">Magweto</td><td>Non-gestation</td><td>September 2020</td><td>348</td><td>7</td><td>2.0 (0.98–4.1)</td></tr><tr><td>Pregnancy</td><td>October 2020</td><td>257</td><td>27</td><td>10.5 (7.3–14.9)</td></tr><tr><td>Parturition</td><td>November 2020</td><td>228</td><td>82</td><td>35.9 (30–42.4)</td></tr><tr><td>4–6 months juveniles</td><td>March 2021</td><td>242</td><td>84</td><td>34.7 (28.9–40.9)</td></tr><tr><td>4–6 months juveniles</td><td>April 2021</td><td>241</td><td>40</td><td>16.6 (12.4–21.8)</td></tr><tr><td>Non-gestation</td><td>June 2021</td><td>242</td><td>17</td><td>7.02 (4.4–10.9)</td></tr><tr><td colspan="2">Overall prevalence</td><td></td><td>1558</td><td>257</td><td>16.5 (14.74–18.42)</td></tr><tr><td rowspan="6">Magweto M. gigas</td><td>Non-gestation</td><td>September 2020</td><td>2</td><td>2</td><td>100 (34.2–100)</td></tr><tr><td>Pregnancy</td><td>October 2020</td><td>124</td><td>70</td><td>56.4 (47.7–64.9)</td></tr><tr><td>Parturition</td><td>November 2020</td><td>74</td><td>26</td><td>35.1 (25.2–46.5)</td></tr><tr><td>4–6 months juveniles</td><td>March 2021</td><td>42</td><td>8</td><td>19.1 (9.98–33.3)</td></tr><tr><td>4–6 months juveniles</td><td>April 2021</td><td>39</td><td>9</td><td>23.1 (12.7–38.3)</td></tr><tr><td>Non-gestation</td><td>June 2021</td><td>27</td><td>5</td><td>18.5 (8.18–36.7)</td></tr><tr><td colspan="2">Overall prevalence</td><td></td><td>308</td><td>120</td><td>38.96 (33.68–44.51)</td></tr></tbody></table>

At both sites, the reproductive season, i.e., from gestation to weaning, was observed from September to February for the predominant insectivorous bat species. At both sites, the different insectivorous bat families observed were found to be synchronous regarding their reproductive cycles, and consensus reproduction periods were determined based on the literature and observations (Tables 1 and S1). The larger insectivorous *M. gigas* is a migratory bat species and was mainly present at Magweto cave from October 2020 to March 2021 during gestation, parturition and 4–6 month juvenile periods, and the reproductive cycle was also synchronous with that of the smaller insectivorous bats (Tables 1 and S1).
## 3.2. *Bat Species*

By using BLAST to compare the newly obtained 800 bp Cytochrome B sequences of the CoV positive samples to the CytB sequences available in GenBank we identified different bat species at the two sampling sites. We identified at least seven species; Hipposideros caffer and Miniopterus mossambicus were present at both sites. At Chirundu farm, we also identified Rhinolophus landeri and Nycteris thebaica, whereas at Magweto cave we also characterized Macronycteris gigas, Cleotis percivali and Rhinolophus simulator. The presence of these species was confirmed by capture sessions. We also characterized bats from the Rhinolophus genus, close to Rhinolophus eloquens from the captured animals. Results obtained by BLAST analysis were confirmed by maximum likelihood phylogenetic analysis (Figure S1).

## 3.3. *Prevalence and Seasonality of RNA Coronaviruses at the Bat Community Level*

The overall prevalence of coronaviruses at Chirundu farm site was 23.7% (95% CI: 22.08–26.13). The lowest prevalence was 1.35% (95% CI: 0.5–3.4) observed in samples collected in October 2020, and prevalence increased between December and March 2021, corresponding to the lactation, weaning and 4–6 month juvenile periods, with the highest prevalence of 44.19% (95% CI: 37.4–45.2) in February, corresponding to the weaning period (Table 1, Figure 2A). At Magweto cave site the sampling prevalence was calculated for the insectivorous bat species and the *Macronycteris gigas* separately. For *M. gigas* the overall prevalence of CoVs was 38.9% (95% CI: 33.68–44.51) with a peak of 56.4% (95% CI: 47.7–64.9) in October 2020 during the gestation period (Table 1, Figure 2B). For the bat species which comprised Hipposideros spp., Rhinolophus spp., Nycteris spp. and Miniopterus spp., the overall prevalence was 16.5% (95% CI:14.74–18.42) with the lowest prevalence (2%; 95% CI: 0.98–4.1)) detected in September 2020 and the highest prevalence detected in November 2020 (35.9%; 95% CI: 30–42.4), corresponding to the period of gestation and parturition, and March 2021 (34.7%; 95% CI: 28.9–40.9), which corresponds to the juvenile periods (Table 1, Figure 2C). Prevalence during the weaning period for this site remains unknown as no samples were collected during that time.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** *Cont.*</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 2.** Results of the longitudinal sampling at two sites for coronaviruses. (A) Estimation of the coronavirus prevalence (with 95% CI) at Chirundu site, (B,C) the same for bats and *Macronycteris gigas* (respectively) at Magweto cave site. The areas colored green, cyan, navy blue, magenta, red and pink show periods of gestation, parturition, lactation, weaning, 4–6 months old juvenile and non-gestation periods, respectively, observed in the dominant bat species at both sites.</figcaption></figure>

Results from the first GLMM showed higher coronavirus prevalence associated with the parturition/lactation (odds ratio (OR) = 5.57, 95% CI) = (1.16–29.56), *p*-value = 0.012) and weaned juveniles of 4 to 6 month old periods (OR = 4.22 (0.83–21.35, *p* = 0.037)) at the Magweto site. Results from the second GLMM showed higher coronavirus prevalence associated with the parturition/lactation period (OR = 6.04 (1.85–22.62, *p* = 0.0014)) and the weaning period (OR = 5.29 (1.16–26.56, *p* = 0.0016)) at the Chirundu site, while lower coronavirus prevalence was associated with the pregnancy period (OR = 0.165 (0.044–0.567, *p* = 0.0018)).
The three viral groups analyzed at Chirundu site, ASVG-01, ASVG-05 and BSVG-01, all showed varying prevalence of CoVs, with peaks in detection during different seasons (Figure S3). For ASVG-01, the highest prevalence was detected in the 4–6 month old juvenile periods, whereas the highest peaks for ASVG-05 were during the lactation and non-gestation periods (Figure S3). For the betacoronavirus sub-viral group BSVG-01, the highest prevalence corresponded to the weaning period (Figure S3).

### 3.4. Genetic Diversity of Coronaviruses at Chirundu and Magweto Sites

A high genetic diversity of alpha and beta coronaviruses was observed in the insectivorous bat ecosystem at both sites. We amplified and sequenced a total of 532 *Coronaviruses* from both sites: 307 alphacoronaviruses (Magweto cave N = 103; Chirundu farm N = 204) and 225 betacoronaviruses (Magweto cave N = 101; Chirundu farm N = 124) (Table 2). A total of 279 out of the 532 CoV sequences were from *Hipposideros* species, 64 from *Macronycteris* spp., 89 from *Rhinolophus* spp., 17 from *Nycteris* spp. and 6 from *Miniopterus* spp. For the remaining 77 CoV sequences, the bat genus/species could not be determined mainly owing to bad sequence qualities (Table 2).

We observed different sub-viral groups with phylogenetic analyses (Figures 3 and 4) which we named from A-SVG-01 to A-SVG-08 for alphacoronaviruses and B-SVG-01 to B-SVG-08 for betacoronaviruses (Table 2, Figures 3 and 4).

### 3.5. Alphacoronaviruses

A-SVG-01 to A-SVG-03 groups belonged to the *Duvinacovirus* sub-genus together with the Human CoV 229E strain (Figure 3). Of the sequences comprising the A-SVG-01 group, 8 of 114 were identified at Magweto cave with four from *H. caffer* and four from *M. gigas*. The 99 remaining sequences were obtained from Chirundu farm with most of them amplified from *H. caffer* (N = 89) and the rest from *Rhinolophus* (N = 3) and *Nycteris* (N = 7) species (Table 2). Viruses belonging to the A-SVG-01 group were found at Chirundu site from August 2020 to May 2021 and at Magweto site, from March to June 2021 only (Figure S2). A-SVG-02 group was mainly composed of sequences identified in *M. gigas* (N = 41) sampled at Magweto. We also identified four other sequences obtained from Chirundu farm: two from *Rhinolophus* species and two from unclassified bat species (Figure 3, Table 2). These A-SVG-02 group sequences were collected in September and November 2020 at Magweto cave and from September to October 2020 at Chirundu site (Figure S2). The A-SVG-03 group was a subclade specific to *Hipposideros* spp. from Magweto site collected between September 2020 and March 2021. However, in this subgroup, two *Rhinolophus* and an unspecified bat species were also infected by this virus (Figures 3 and S2, Tables 2 and S1).

The A-SVG-04 subgroup belongs to the *Setravirus* sub-genus together with the human CoV NL63 strain. Only one sequence, identified at Chirundu from a *Miniopterus mossambicus* bat species in May 2021 comprised this group (Figures 3 and S2, Table 2).

A-SVG-05 and -06 subgroups belonged to the *Decacovirus* clade. The A-SVG-05 sub-group mainly comprised sequences isolated from *Rhinolophus* bats (N = 52) from September 2020 to July 2021 at Chirundu site. Furthermore, *Hipposideros* (N = 1), *Nycteris* (N = 1) and *Miniopterus* (N = 2) bat species were also identified as carriers of this virus at this site in the same period. At Magweto site, six A-SVG-05 sequences from *Hipposideros* species (N = 2) were identified in June 2021 and from *Macronycteris* spp. (N = 4) in April and June 2021. A-SVG-06 group sequences (N = 13) were only detected at Magweto site from *Rhinolophus* (N = 9), *Hipposideros* (N = 1), *Miniopterus* (N = 1) and from unspecified (N = 2) bat species. This viral group was present from September to November 2020 at this site (Figures 3 and S2, Table 2).
<table><caption>**Table 2.** Longitudinal detection of different viral groups in Magweto and Chirundu sites.</caption><thead><tr><th rowspan="2">Viral Group</th><th rowspan="2">Total Number of Sequences</th><th colspan="2">Number of Sequence per Site</th><th colspan="4">Number of Sequences per Bat Species</th><th colspan="2">Longitudinal Detection of the Different Viral Groups</th></tr><tr><th>Magweto</th><th>Chirundu</th><th>Hipposideros spp.</th><th>Macronycteris spp.</th><th>Rhinolophe spp.</th><th>Nycteris spp.</th><th>Miniopterus spp.</th><th>Unknown</th><th>Magweto</th><th>Chirundu</th></tr></thead><tbody><tr><td>A-SVG-01</td><td>114</td><td>8</td><td>106</td><td>89</td><td>4</td><td>3</td><td>7</td><td>0</td><td>11</td><td>March 2021 to June 2021</td><td>August 2020 to May 2021</td></tr><tr><td>A-SVG-02</td><td>45</td><td>42</td><td>3</td><td>0</td><td>41</td><td>2</td><td>0</td><td>0</td><td>2</td><td>September 2020 to November 2020</td><td>October 2020 to November 2020</td></tr><tr><td>A-SVG-03</td><td>32</td><td>32</td><td>0</td><td>29</td><td>0</td><td>2</td><td>0</td><td>0</td><td>1</td><td>September 2020 to March 2021</td><td>-</td></tr><tr><td>A-SVG-04</td><td>1</td><td>0</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td><td>-</td><td>May 2021</td></tr><tr><td>A-SVG-05</td><td>84</td><td>6</td><td>78</td><td>7</td><td>4</td><td>52</td><td>1</td><td>3</td><td>17</td><td>June 2021</td><td>September 2020 to July 2021</td></tr><tr><td>A-SVG-06</td><td>13</td><td>13</td><td>0</td><td>1</td><td>0</td><td>9</td><td>0</td><td>1</td><td>2</td><td>September 2020 to November 2020</td><td>-</td></tr><tr><td>A-SVG-07</td><td>2</td><td>2</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>October 2020</td><td>-</td></tr><tr><td>A-SVG-08</td><td>16</td><td>0</td><td>16</td><td>0</td><td>0</td><td>8</td><td>1</td><td>0</td><td>7</td><td>-</td><td>August-November-December 2020-July 2021</td></tr><tr><td>sub Total</td><td>307</td><td>103</td><td>204</td><td>126</td><td>49</td><td>76</td><td>9</td><td>6</td><td>41</td><td></td><td></td></tr><tr><th colspan="12" style="text-align: center;">Number of Sequence per site</th></tr><tr><th colspan="12" style="text-align: center;">Number of sequences per Bat species</th></tr><tr><td>B-SVG-01</td><td>105</td><td>5</td><td>100</td><td>76</td><td>2</td><td>4</td><td>1</td><td>0</td><td>22</td><td>March to April 2021</td><td>November 2020 to July 2021</td></tr><tr><td>B-SVG-02</td><td>21</td><td>21</td><td>0</td><td>19</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>November 2020 to April 2021</td><td>-</td></tr><tr><td>B-SVG-03</td><td>6</td><td>0</td><td>6</td><td>2</td><td>0</td><td>1</td><td>0</td><td>0</td><td>3</td><td>-</td><td>February 2021 May and June 2021</td></tr><tr><td>B-SVG-04</td><td>13</td><td>13</td><td>0</td><td>0</td><td>13</td><td>0</td><td>0</td><td>0</td><td>0</td><td>September 2020 to June 2021</td><td>-</td></tr><tr><td>B-SVG-05</td><td>51</td><td>51</td><td>0</td><td>48</td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td><td>November 2020 to July 2021</td><td>March 2021</td></tr><tr><td>B-SVG-06</td><td>8</td><td>8</td><td>0</td><td>0</td><td>0</td><td>5</td><td>1</td><td>0</td><td>2</td><td>October 2020 to March 2021</td><td>-</td></tr><tr><td>B-SVG-07</td><td>4</td><td>0</td><td>4</td><td>0</td><td>0</td><td>2</td><td>0</td><td>0</td><td>2</td><td>-</td><td>July 2021</td></tr><tr><td>B-SVG-08</td><td>17</td><td>3</td><td>14</td><td>8</td><td>0</td><td>1</td><td>6</td><td>0</td><td>2</td><td>April 2021</td><td>August 2020 to Mach 2021</td></tr><tr><td>sub Total</td><td>225</td><td>101</td><td>124</td><td>153</td><td>15</td><td>13</td><td>8</td><td>0</td><td>36</td><td></td><td></td></tr><tr><td>TOTAL</td><td>532</td><td>204</td><td>328</td><td>279</td><td>64</td><td>89</td><td>17</td><td>6</td><td>77</td><td></td><td></td></tr></tbody></table>
<figure><img src="image_8.png" /><figcaption>**Figure 3.** Phylogenetic tree of *AlphaCoV* partial *RdRp* gene. The sequences detected at Chirundu site are in red color and from Magweto site in blue color. Pink = more sequences detected at Chirundu than Magweto within the same subclade and purple = more sequences detected at Magweto site than Chirundu site within the subclade. The tree was built using the maximum likelihood method based on the GTR + G4 + I model. The robustness of nodes was assessed with 1000 bootstrap replicates. Bootstrap values >70 are denoted with asterisks, and those <70 are not shown.</figcaption></figure>

A-SVG-07 subgroup belonged to the *Minunacovirus* clade, with one sequence obtained from *Miniopterus* spp. and one from an unspecified species (Figure 3, Table 2).

A-SVG-08 subgroup belonged to the *Rhinacovirus* clade and comprised sequences (N = 16) obtained from samples collected in August, November and December 2020 at Magweto site and in July 2021 at Chirundu site. Eight were from *Rhinolophus*, one from *Nycteris* and seven from uncharacterized bat species (Figures 3 and S2, Table 2).
<figure><img src="image_9.png" /><figcaption>**Figure 4.** Phylogenetic tree of *BetaCoV* partial *RdRp* gene. The sequences detected at Chirundu site are in red color and at Magweto site in blue color. Pink = more sequences detected at Chirundu than Magweto site within the same subclade. The tree was built using the maximum likelihood method based on the GTR + G4 + I model. The robustness of nodes was assessed with 1000 bootstrap replicates. Bootstrap values >70 are denoted with asterisks, and those <70 are not shown.</figcaption></figure>

### 3.6. *Betacoronaviruses*

B-SVG-01 to B-SVG-05 groups belonged to the *Hibecovirus* clade. B-SVG-01 was the most important subclade with 105 sequences: five were from *Hipposideros* (N = 3) and *Macronycteris* (N = 2) bat species sampled at Magweto cave in March and April 2021. For the 100 remaining sequences, 73 were identified from *Hipposideros*, four from *Rhinolophus*, one from *Nycteris* and 22 from uncharacterized bat species collected from November 2020 to July 2021 (Figures 4 and S2, Table 2). B-SVG-02 and -04 were specific to *Hipposideros caffer* and *Macronycteris gigas* bat species from Magweto cave and were detected from November 2020 to April 2021 and from September to June 2021, respectively (Figures 4 and S2, Table 2). B-SVG-03 subclade comprised six sequences isolated at Chirundu site from *Hipposideros* (N = 2), *Rhinolophus* (N = 1) and from undetermined bat species (N = 3). This subclade was detected in February, May and June 2021 (Figures 4 and S2, Table 2). B-SVG-05 was composed of 52 sequences; 51 were detected at Magweto site and one at Chirundu site. Only *Hipposideros* bat species were infected by this subclade, and these viruses were present at Magweto site from November 2020 to July 2021 and only in March 2021 at Chirundu site (Figures 4 and S2, Table 2).
B-SVG-06 and -07 belong to the *Sabercovirus* clade together with SARS-CoV-1 and -2 as well as numerous SARS-related bat CoVs isolated from *Rhinolophus* bat species in Asia and Africa. The B-SVG-06 subclade was specific to Magweto site with five sequences characterized from *Rhinolophus* bat species, one from *Nycteris* bat species and two from uncharacterized bat species, and the presence of this subclade was detected from October 2020 to March 2021 (Figure 4 and Figure S2, Table 2). As for B-SVG-07 subclade, it was specific to Chirundu site and comprised only four sequences, (N = 2) from *Rhinolophus* bat species and (N = 2) from uncharacterized bat species collected in July 2021 (Figures 4 and S2, Table 2).

B-SVG-08 belonged to the Merbecovirus clade. This subclade was mainly composed of sequences isolated at Chirundu site from *Hipposideros* (N = 8), *Nycteris* (N = 6), *Rhinolophus* (N = 1) and two unspecified bat species sampled from August 2020 to March 2021 (Figures 4 and S2, Table 2). Three sequences were isolated at Magweto site from *Hipposideros* (N = 2) and *Rhinolophus* (N = 1) bat species collected in April 2021.

**4. Discussion**

The majority of bat coronavirus studies in Africa mainly focused on genetic diversity and prevalence of CoVs, mostly based on single sampling sessions [43]. However, longitudinal studies covering bat reproduction cycles including the gestation, parturition, lactation, weaning and juvenile periods are crucial to further understand whether viral diversity and viral shedding vary over time and to assess the high-risk season for CoV spillovers from bats to human and domestic animal populations in close contact with these bat populations [29,43,44]. Owing to logistical challenges, few longitudinal studies have been conducted on bat coronaviruses. Moreover, these studies highlighted the variability of viral shedding over time and particularly during the reproductive season [29,31,43,45,46]. We performed a longitudinal study in two multi-bat species roosting sites in Zimbabwe. We collected fecal samples monthly from August 2020 to July 2021. Unfortunately, due to the travel restrictions of the lockdown periods during the COVID-19 pandemic, only nine months at Chirundu farm site and six months at Magweto cave site were sampled instead of the expected 12 month follow up. Nonetheless, we collected more than 3000 fecal samples covering reproductive and non-reproductive periods of the bat colonies.

A high genetic diversity was observed among the newly identified coronaviruses including alphacoronaviruses from five sub-genera, *Duvinacovirus*, *Setracovirus*, *Decacovirus*, *Minunacovirus* and *Rhinacovirus*, as well as betacoronaviruses from three sub-genera, *Hibecovirus*, *Sarbecovirus* and *Merbecovirus*. All these sub-genera except for *Decacovirus*, *Minuncavirus* and *Rhinacovirus* have been previously described in Africa [43]. Our study enlarges the knowledge on CoV diversity in bats from Africa.

The *Hipposideros* bat species were the most prominent carriers of CoVs, representing 52.4% of our CoV sequences, followed by *Rhinolophus* bats (16.7%), *Macronycteris gigas* species (12%) and more anecdotally some *Nycteris* species (3.1%) and very few *Miniopterus* bats (1.1%). Bat species were only characterized for the positive CoV samples, and it is therefore difficult to link these results with bat species community constitution at both sites during the time of this study. Nonetheless, at Chirundu farm we visually observed a shift of bat populations in June and July. Indeed, *Hipposideros* bat species were predominant before June, whereas *Rhinolophus* species became predominant in June and July, a shift confirmed by cytochrome B analyses of all coronavirus positive samples. Thus, in June and July, most of our CoV positive samples were detected from *Rhinolophus bat* spp.

The *Duvinacovirus* sub-genus comprises the Hu-CoV 229E strain known to induce common cold in humans [32]. This virus originates from bats and more specifically from *Hipposideros* species [3]. In our study we confirmed *Hipposideros* species as the main carriers of CoV-229E-related strains in Zimbabwe. However, we also observed a *M. gigas* specific *Duvinacovirus*-related subclade. This viral strain was also detected in a few *Rhinolophus* and *Nycteris* bats, and presence in these species can be attributed to cross-species transmission from *Hipposideros* or *Macronycteris* bat species. However, cross-contamination among fecal
samples collected on the floor cannot be excluded. Our phylogenetic analyses also showed phylogeographic subclades according to the sampling site.

We also identified a CoV clade belonging to the *Setracovirus* sub-genus together with the Hu-CoV NL63 strain which also induces common cold in humans [32]. To date NL63-related bat CoV have been identified to originate from the *Triaenop* bat genus [43]. In our study we characterized a NL63-CoV-related strain from a *Miniopterus* bat species. However, only one sequence was recovered, suggesting a very rare circulation of this virus in our study sites. Moreover, it is possible that this virus was acquired by this bat through contact with a *Triaenops* bat species from another site. Nonetheless, only one *Triaenops* bat species has been recorded in this region, *Triaenops afer*, which seems to be marginally present in the Eastern region of Zimbabwe [47], thus suggesting that another bat species could also be a carrier of this CoV strain in Zimbabwe. Within the *Decacovirus* sub-genus, we observed two different subclades corresponding to the two sites with the majority of *Rhinolophus* bat species as carriers of this virus. Additionally, we also observed some cross-species transmissions of this virus from *Rhinolophus* bats to *Hipposideros*, *Nycteris* and *Miniopterus* bat species, although here too cross-contamination among fecal samples collected on the floor cannot be excluded. The last two alpha subgroup viruses, A-SVG-07 and -08, were detected at Magweto and Chirundu sites, respectively, with A-SVG-07 belonging to the sub-genus *Minunacovirus* with only two sequences characterized and A-SVG-08 mainly present in *Rhinolophus* bat species belonging to the *Rhinacovirus* sub-genus.

The *Hibecovirus* sub-genus was the most important *Betacoronavirus* clade, representing 88% of all β-CoV sequences characterized in our study. Like the *Duvinacovirus* sub-genus, the majority of our *Hibecovirus* sequences were identified in *Hipposideros* bat species and a few *Rhinolophus* and *Nycteris* bats. The latter species might have been infected by cross-species transmission of the virus from *Hipposideros* bats. We also evidenced the presence of *Sarbecovirus* which is composed of SARS-CoV-01 and -02 strains and bat CoVs isolated from *Rhinolophus* bat species. This sub-genus was characterized at both sites. At Magweto site only eight samples were positive along with four from Chirundu, highlighting the low circulation of these viruses at both sites despite detection from October 2020 to March 2021 at Magweto site. Finally, our last β-CoV sub viral group belonged within the *Merbecovirus* sub-genus which also comprised the MERS-CoV virus. Both *Hipposideros* and *Nycteris* bat species were infected as well as one *Rhinolophus* species. Bat CoVs from this sub-genus were mainly associated with the Vespertilioniforme bat sub-order [43,48] from *Molosidae*, *Nycteridae* and *Emballonuridae* bat families [49–51]. In this study we also evidenced *Merbecovirus* in the Pteropodiforme bat sub-order and particularly in *Hipposideridae* and *Rhinolophidae* bat families which extended the range of bat families as carriers of *Merbecoviruses*.

Globally, we identified bat CoVs belonging to five out of the eight CoV sub-genera previously described for *Alphacoronavirus* genus and three out of the four sub-genera described for *Betacoronavirus* genus. We also observed what might indicate cross-species transmissions within the same site, but possibility of contamination cannot be excluded. Thus, the impact of these potential transmissions to new bat hosts has yet to be elucidated as well as whether these new hosts play a role in the dissemination of these viruses to other animal species or to humans. Additionally, we observed identical viruses in both sites suggesting that the two studied bat communities are somehow interconnected. Magweto and Chirundu sites are situated at 250 km apart, and numerous caves are present in-between, allowing bat population exchanges [52] alongside with their microorganisms over a large area.

CoV RNA was detected in 24% of the bat samples at Chirundu site and in 16.5% of the insectivorous bat samples and 39% of the *M. gigas* bat samples at Magweto site. However, at Magweto site this result was probably underestimated since we were not able to collect samples from December 2020 to March 2021, which corresponds to the lactation and weaning periods with numerous juveniles in the cave. In other longitudinal surveys, high RNA CoV detection was found to be associated with juvenile bats [29,31,46]. Pathogens respond to the influx of susceptible hosts (new recruits) through increased
transmission following the birth pulse [53–55]. At Chirundu site we observed a peak of RNA CoV prevalence during the lactation, weaning and juvenile periods with 37%, 44% and 35% of RNA CoV detection, respectively. We observed the same trend at Magweto site in the bats with an RNA CoV detection of 35.9% during the parturition period in November 2020 and of 34.7% during the juvenile period in March 2021. The transmission and prevalence of pathogens is facilitated by higher proportion of juveniles who are immunologically immature [31,55] explaining thus the high prevalence of CoVs during these seasons. Furthermore, the waning protection of maternal antibodies upon birth and weaning leads to increased prevalence in virus detection in bats [21]. Due to the similarity in dominant bat species between the two sites, it is likely that the Magweto site showed a similar trend in prevalence of CoVs during the unsampled and unobserved months corresponding to the lactation and weaning periods as observed at Chirundu site. Additionally, within the *M. gigas* bat species the peak of RNA CoV detection (56.4%) was in October 2020 during the gestation period which was also the period with the most important *M. gigas* population in the cave. It is important to bear in mind that this study was conducted over 1 year only, and that sampling over a long period of time would be needed to confirm the annual seasonal trends observed in this study. Moreover, the prevalence peaks of the analyzed viral subgroups which occur at different stages of the reproductive season show that the global seasonal trend observed at the CoV family level is not necessarily the case for the very viral subgroup. Further research should thus be performed on specific CoV groups which are of great importance to domestic animal and public health.

Coronavirus evolution and emergence in new hosts is driven by different factors such as recombination, horizontal gene transfer, gene duplication and alternative open reading frames which expand their functional and adaptative capacity for the current and new host. During the high viral shedding season, bats are carriers of quasi-species pools of viruses which contribute to their genetic diversity and increase their potential to jump and emerge in a new species [56]. Due to the great genetic diversity observed at both sites as well as the important RNA bat CoV prevalence during the reproductive season, our sites could be potential hotspots of new CoV spillover and potential zoonotic emergence.

The risk of CoV spillover to human population is higher at Magweto site than Chirundu site. At Chirundu, the roosting area is in the basement of an old, damaged building with very few direct interactions between the farm workers and the bat populations. Contrary to this, Magweto site could be more susceptible to a viral spillover event. Human populations use the cave for different activities. Throughout the year, bat guano is collected inside the cave and applied as fertilizers in neighboring fields and sold for commercial earnings, generating fecal dust which is inhaled by workers. The workers and buyers are thus directly exposed to bat dejections and aerosols. Moreover, this cave is considered to be sacred, with traditional and religious activities inside the cave, and finally some people hunt *M. gigas* bats for consumption. All these cause direct interactions between bats and humans. The interface between bats and human populations at Magweto cave is thus very tight, and increased surveillance of viral circulation and zoonotic transmissions should be recommended, especially during periods with peaks of viral shedding, between October and March.

**5. Conclusions**

Altogether, our biological results associated with the human activities at Magweto site and to a lesser extent at Chirundu site highlighted the need to further document the viral infectious risk in both human and domestic animal populations surrounding bat habitats in Zimbabwe. Particularly, the human population using caves is constantly exposed to viral shedding, and all factors required for viral spillover are met at least for the first and second stages of viral disease emergence [57,58]. This study falls within the One Health approach where the established results can be used to form a framework for further studies in human populations living in proximity to the study areas to further extrapolate the
extent, if any, of the bat coronavirus circulation in humans and domestic animals through a serological study.

Our study focused on coronaviruses; however, bats host a plethora of potential infectious viruses and numerous bat species known to harbor highly fatal viruses are present in Zimbabwe. Therefore, further studies are needed to better understand the bat viral ecology in Zimbabwe and to better assess the risk of viral disease emergence in this geographic area [43,48].

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14040781/s1, Table S1: The reproductive season of dominant bat species across months of the year. Figure S1: Phylogenetic tree of bat species identified by cytochrome B; Figure S2: Timeline of sub-viral group detection at Chirundu and Magweto sites; Figure S3: Results of the prevalence of CoVs in Chirundu site sub-viral groups.

**Author Contributions:** Conceptualization: V.C., H.D.N., M.B., M.P., A.A., L.G. and F.L.; Methodology: V.C., L.G., J.C., H.D.N., D.P., E.G., G.M. (Gift Matope) and F.L.; Software: V.C., J.C., H.D.N. and F.L.; Validation: V.C., H.D.N., J.C. and F.L.; Formal analysis: V.C., H.D.N., J.C. and F.L.; Investigation: V.C., V.P., N.C., G.M. (Getrude Mashura) and F.L.; Resources: D.P., G.M. (Gift Matope), M.P., M.B., H.D.N., L.G. and F.L.; Data curation: V.C., H.D.N., J.C. and F.L.; Writing—original draft preparation: V.C. and F.L.; Writing—review and editing: H.D.N., M.B., M.P., A.A., L.G., V.P., E.G., D.M., J.C., D.P., G.M. (Gift Matope) and F.L.; Supervision: H.D.N., D.P., G.M. (Gift Matope), M.B. and F.L.; Project administration: L.G., M.B., H.D.N. and F.L.; Funding acquisition: M.P., M.B., H.D.N., L.G. and F.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by grants of the Agence Nationale de Recherche sur le SIDA et les Maladies Infectieuses (ANRS-COV10), the ANR (the French National Research Agency) under the "Investissements d'avenir" programme with the reference ANR-16-IDEX-0006) and the French Ministry of Europe and Foreign Affairs (Fond de Solidarité pour les Projets Innovants, les sociétés civiles, la francophonie et le développement humain—CAZCOM Project) and the Agence Française pour le Développement (AFD-PACMAN-Project).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** We thank Billy Butete for his field assistance. The authors gratefully acknowledge and thank the International Atomic Energy Agency for making available the sequencing services. We thank the Research Council of Zimbabwe for approving this study (research registration certificate no. 03006) and the Hurungwe Rural District council for their assistance and facilitation. We thank the Animal Research Ethics Committee of Zimbabwe for their approval (ref number 002/2017). This work was conducted within the framework of the Research Platform "Production and Conservation in Partnership" (RP-PCP).

**Conflicts of Interest:** We declare that we have no conflict of interest.

## References

1. ICTV Positive Sense RNA Viruses: Coronaviridae. VIRUS Taxon. 2020. Available online: https://talk.ictvonline.org/ (accessed on 1 December 2021).

2. Al-Qahtani, A.A. Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2): Emergence, history, basic and clinical aspects. *Saudi J. Biol. Sci.* **2020**, *27*, 2531–2538. [CrossRef] [PubMed]

3. Corman, V.M.; Muth, D.; Niemeyer, D.; Drosten, C. Hosts and Sources of Endemic Human Coronaviruses. *Adv. Virus Res.* **2018**, *100*, 163–188. [PubMed]

4. Li, K.; Lin, X.; Wang, W.; Shi, M.; Guo, W.; Zhang, X.; Xing, J.; He, J.; Wang, K.; Li, M.; et al. Isolation and characterization of a novel arenavirus harbored by Rodents and Shrews in Zhejiang province, China. *Virology* **2015**, *476*, 37–42. [CrossRef] [PubMed]

5. Monchatre-Leroy, E.; Boué, F.; Boucher, J.M.; Renault, C.; Moutou, F.; Gouilh, M.A.; Umhang, G. Identification of alpha and beta coronavirus in wildlife species in france: Bats, rodents, rabbits, and hedgehogs. *Viruses* **2017**, *9*, 364. [CrossRef]

6. Wertheim, J.O.; Chu, D.K.W.; Peiris, J.S.M.; Kosakovsky Pond, S.L.; Poon, L.L.M. A Case for the Ancient Origin of Coronaviruses. *J. Virol.* **2013**, *87*, 7039–7045. [CrossRef]
7. Macnaughton, M.R.; Thomas, B.J.; Davies, H.A.; Patterson, S. Infectivity of human coronavirus strain 229E. *J. Clin. Microbiol.* **1980**, *12*, 462–468. [CrossRef]

8. Van Der Hoek, L.; Pyrc, K.; Jebbink, M.F.; Vermeulen-Oost, W.; Berkhout, R.J.M.; Wolthers, K.C.; Wertheim-Van Dillen, P.M.E.; Kaandorp, J.; Spaargaren, J.; Berkhout, B. Identification of a new human coronavirus. *Nat. Med.* **2004**, *10*, 368–373. [CrossRef]

9. Razanajatovo, N.H.; Nomenjanahary, L.A.; Wilkinson, D.A.; Razafimanahaka, J.H.; Goodman, S.M.; Jenkins, R.K.; Jones, J.P.G.; Heraud, J.-M. Detection of new genetic variants of Betacoronaviruses in Endemic Frugivorous Bats of Madagascar. *Virol. J.* **2015**, *12*, 42. [CrossRef]

10. Woo, P.C.Y.; Lau, S.K.P.; Li, K.S.M.; Poon, R.W.S.; Wong, B.H.L.; Tsoi, H.; Yip, B.C.K.; Huang, Y.; Chan, K.; Yuen, K. Molecular diversity of coronaviruses in bats. *Virology* **2006**, *351*, 180–187. [CrossRef]

11. Hu, B.; Ge, X.; Wang, L.; Shi, Z. Bat origin of human coronaviruses: Emerging and re-emerging pathogens in humans and animals Susanna Lau Positive-strand RNA viruses. *Virol. J.* **2015**, *12*, 209. [CrossRef]

12. Zhou, P.; Yang, X.; Wang, X.; Hu, B.; Zhang, L.; Zhang, W.; Guo, H.; Jiang, R.; Liu, M.; Chen, Y.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. *Nature* **2020**, *579*, 270–273. [CrossRef] [PubMed]

13. Vlasova, A.; Saif, L.J. Bovine Coronavirus and the Associated Diseases. *Front. Vet. Sci.* **2021**, *8*, 1–14. [CrossRef]

14. Woo, P.C.Y.; Lau, S.K.P.; Huang, Y.; Yuen, K.Y. Coronavirus Diversity, Phylogeny and Interspecies Jumping. *Exp. Biol. Med.* **2009**, *234*, 1117–1127. [CrossRef] [PubMed]

15. Luis, A.D.; O'Shea, T.J.; Hayman, D.T.S.; Wood, J.L.N.; Cunningham, A.A.; Gilbert, A.T.; Mills, J.N.; Webb, C.T. Network analysis of host—Virus communities in bats and rodents reveals determinants of cross-species transmission. *Ecol. Lett.* **2015**, *18*, 1153–1162. [CrossRef] [PubMed]

16. Drexler, J.F.; Corman, V.M.; Müller, M.A.; Maganga, G.D.; Vallo, P.; Binger, T.; Gloza-Rausch, F.; Cottontail, V.M.; Rasche, A.; Yordanov, S.; et al. Bats host major mammalian paramyxoviruses. *Nat. Commun.* **2012**, *3*, 796. [CrossRef] [PubMed]

17. Peiris, J.S.M.; Guan, Y.; Yuen, K.Y. Severe acute respiratory syndrome. *Nat. Med.* **2004**, *10*, S88–S97. [CrossRef] [PubMed]

18. El-Sayed, A.; Kamel, M. Coronaviruses in humans and animals: The role of bats in viral evolution. *Environ. Sci. Pollut. Res. Int.* **2021**, *28*, 19589–19600. [CrossRef]

19. Li, X.; Zai, J.; Zhao, Q.; Nie, Q.; Li, Y.; Foley, B.T.; Chaillon, A. Evolutionary history, potential intermediate animal host, and cross-species analyses of SARS-CoV-2. *J. Med. Virol.* **2020**, *92*, 602–611. [CrossRef]

20. Delaune, D.; Hul, V.; Karlsson, E.A.; Hassanin, A.; Ou, T.P.; Baidaliuk, A.; Gámbaro, F.; Prot, M.; Tu, V.T.; Chea, S.; et al. A novel SARS-CoV-2 related coronavirus in bats from Cambodia. *Nat. Commun.* **2021**, *12*, 6563. [CrossRef]

21. Coertse, J.; De Vries, L.; Geldenhuys, M.; Mortlock, M. Bat-borne viruses in Africa: A critical review. *J. Zool.* **2020**, *311*, 77–98.

22. Ar Gouilh, M.; Puechmaille, S.J.; Diancourt, L.; Vandenbogaert, M.; Serra-Cobo, J.; Lopez Roig, M.; Brown, P.; Moutou, F.; Caro, V.; Vabret, A.; et al. SARS-CoV related Betacoronavirus and diverse Alphacoronavirus members found in western old-world. *Virology* **2018**, *517*, 88–97. [CrossRef] [PubMed]

23. Letko, M.; Seifert, S.N.; Olival, K.J.; Plowright, R.K.; Munster, V.J. Bat-borne virus diversity, spillover and emergence. *Nat. Rev. Microbiol.* **2020**, *18*, 461–471. [CrossRef] [PubMed]

24. Woo, P.C.Y.; Lau, S.K.P. Viruses and Bats. *Viruses* **2019**, *11*, 884. [CrossRef] [PubMed]

25. Latinne, A.; Hu, B.; Olival, K.J.; Zhu, G.; Zhang, L.; Li, H.; Chmura, A.A.; Field, H.E.; Zambrana-Torrelio, C.; Epstein, J.H.; et al. Origin and cross-species transmission of bat coronaviruses in China. *Nat. Commun.* **2020**, *11*, 1–15. [CrossRef]

26. Calisher, C.H.; Childs, J.E.; Field, H.E.; Holmes, K.V.; Schountz, T. Bats: Important Reservoir Hosts of Emerging Viruses. *Clin. Microbiol. Rev.* **2006**, *19*, 531–545. [CrossRef]

27. Role, T.H.E.; Bats, O.F.; Zoonoses, I.N.E. *Investigating the Role of Bats in Emerging Zoonoses: Balancing Ecology, Conservation and Public Health Interests*; Newman, S., Field, H., de Jong, C., Epstein, J.H., Eds.; FAO Animal; FAO: Rome, Italy, 2011; ISBN 1810-1119.

28. Alfelt, A.; Devaux, C.; Serra-Cobo, J.; Frutos, R. Bats, Bat-Borne Viruses, and Environmental Changes. In *Bats; IntechOpen: London, UK*, 2018; pp. 113–132.

29. Wacharapluesadee, S.; Duengkae, P.; Chaiyes, A.; Kaewpom, T.; Rodpan, A.; Yingsakmongkon, S.; Petcharat, S.; Phengsakul, P.; Maneeorn, P.; Hemachudha, T. Longitudinal study of age-specific pattern of coronavirus infection in Lyle's flying fox (Pteropus lylei) in Thailand. *Virol. J.* **2018**, *15*, 1–10. [CrossRef]

30. Anthony, S.J.; Johnson, C.K.; Greig, D.J.; Kramer, S.; Che, X.; Wells, H.; Hicks, A.L.; Joly, D.O.; Wolfe, N.D.; Daszak, P.; et al. Global patterns in coronavirus diversity. *Virus Evol.* **2017**, *3*, 1–15. [CrossRef]

31. Cappelle, J.; Furey, N.; Hoem, T.; Ou, T.P.; Lim, T.; Hul, V.; Heng, O.; Chevalier, V.; Dussart, P.; Duong, V. Longitudinal monitoring in Cambodia suggests higher circulation of alpha and betacoronaviruses in juvenile and immature bats of three species. *Sci. Rep.* **2021**, *11*, 1–11. [CrossRef]

32. Banerjee, A.; Kulcsar, K.; Misra, V.; Frieman, M.; Mossman, K. Bats and coronaviruses. *Viruses* **2019**, *11*, 41. [CrossRef]

33. Bourgarel, M.; Fukenyi, D.M.; Boué, V.; Talignani, L.; Chiweshe, N.; Diop, F.; Caron, A.; Matope, G.; Missé, D.; Liégeois, F. Circulation of Alphacoronavirus, Betacoronavirus and Paramyxovirus in Hipposideros bat species in Zimbabwe. *Infect. Genet. Evol.* **2018**, *58*, 253–257. [CrossRef]

34. Bourgarel, M.; Noël, V.; Fukenyi, D.; Michaux, J.; André, A.; Becquart, P.; Cerqueira, F.; Barrachina, C.; Boué, V.; Talignani, L.; et al. Next-generation sequencing on insectivorous bat guano: An accurate tool to identify arthropod viruses of potential agricultural concern. *Viruses* **2019**, *11*, 1102. [CrossRef] [PubMed]
35. Monadjem, A.; Taylor, P.J.; Cotterill, F.P.D.; Schoeman, M. *Bats of Southern and Central Africa: A Biogeographic and Taxonomic Synthesis*; Wits University Press: Johannesburg, South Africa, 2010.

36. De Nys, H.M.; Mbala Kingebeni, P.; Keita, A.K.; Butel, C.; Thaurignac, G.; Villabona-Arenas, C.J.; Lemarcis, T.; Geraerts, M.; Vidal, N.; Esteban, A.; et al. Survey of ebola viruses in frugivorous and insectivorous bats in Guinea, Cameroon, and the democratic republic of the Congo, 2015–2017. *Emerg. Infect. Dis.* **2018**, *24*, 2228–2240. [CrossRef]

37. Chu, D.K.W.; Leung, C.Y.H.; Gilbert, M.; Joyner, P.H.; Ng, E.M.; Tse, T.M.; Guan, Y.; Peiris, J.S.M.; Poon, L.L.M. Avian coronavirus in wild aquatic birds. *J. Virol.* **2011**, *85*, 12815–12820. [CrossRef]

38. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. *Mol. Biol. Evol.* **2016**, *33*, 1870–1874. [CrossRef] [PubMed]

39. Trifinopoulos, J.; Nguyen, L.T.; von Haeseler, A.; Minh, B.Q. W-IQ-TREE: A fast online phylogenetic tool for maximum likelihood analysis. *Nucleic Acids Res.* **2016**, *44*, W232–W235. [CrossRef] [PubMed]

40. Yoder, A.D.; Yang, Z. Divergence dates for Malagasy lemurs estimated from multiple gene loci: Geological and evolutionary context. *Mol. Ecol.* **2004**, *13*, 757–773. [CrossRef] [PubMed]

41. Kocher, T.D.; Thomas, W.K.; Meyer, A.; Edwards, S.V.; Paabo, S.; Villablanca, F.X.; Wilson, A.C. Dynamics of mitochondrial DNA evolution in animals: Amplification and sequencing with conserved primers. *Proc. Natl. Acad. Sci. USA* **1989**, *86*, 6196–6200. [CrossRef]

42. Agresti, A.; Coull, B.A. Approximate Is Better than "Exact" for Interval Estimation of Binomial Proportions. *Am. Stat.* **1998**, *52*, 119–126.

43. Geldenhuys, M.; Mortlock, M.; Epstein, J.H.; Paweńska, J.T.; Weyer, J.; Markotter, W. Overview of Bat and Wildlife Coronavirus Surveillance in Africa: A Framework for Global Investigations. *Viruses* **2021**, *13*, 936. [CrossRef]

44. Seltmann, A.; Corman, V.M.; Rasche, A.; Drosten, C. Seasonal Fluctuations of Astrovirus, But Not Coronavirus Shedding in Bats Inhabiting Human-Modified Tropical Forests. *EcoHealth* **2017**, *14*, 272–284. [CrossRef]

45. Drexler, J.F.; Corman, V.M.; Wegner, T.; Tateno, A.F.; Zerbinati, R.M.; Gloza-Rausch, F.; Seebens, A.; Müller, M.A.; Drosten, C. Amplification of emerging viruses in a bat colony. *Emerg. Infect. Dis.* **2011**, *17*, 449–456. [CrossRef] [PubMed]

46. Montecino-latorre, D.; Goldstein, T.; Gilardi, K.; Wolking, D.; Wormer, E.V.; Kazwala, R.; Ssebide, B.; Nziza, J.; Sijali, Z.; Cranfield, M.; et al. Reproduction of East-African bats may guide risk mitigation for coronavirus spillover. *One Health Outlook* **2020**, *8*, 1–13. [CrossRef] [PubMed]

47. Van Cakenberghe, V.; Seamark, E. *African Chiroptera Report 2020*; AfricanBats NPC: Pretoria, South Africa, 2020.

48. Wong, A.C.P.; Li, X.; Lau, S.K.P.; Woo, P.C.Y. Global epidemiology of bat coronaviruses. *Viruses* **2019**, *11*, 174. [CrossRef] [PubMed]

49. Anthony, S.J.; Ojeda-Flores, R.; Rico-Chávez, O.; Navarrete-Macias, I.; Zambrana-Torrelio, C.M.; Rostal, M.K.; Epstein, J.H.; Tipps, T.; Liang, E.; Sanchez-Leon, M.; et al. Coronaviruses in bats from Mexico. *J. Gen. Virol.* **2013**, *94*, 1028–1038. [CrossRef]

50. Annan, A.; Baldwin, H.J.; Corman, V.M.; Klose, S.M.; Owusu, M.; Nkrumah, E.E.; Badu, E.K.; Anti, P.; Agbenyega, O.; Meyer, B.; et al. Human betacoronavirus 2c EMC/2012-related viruses in bats, Ghana and Europe. *Emerg. Infect. Dis.* **2013**, *19*, 456–459. [CrossRef]

51. Memish, Z.A.; Mishra, N.; Olival, K.J.; Fagbo, S.F.; Kapoor, V.; Epstein, J.H.; AlHakeem, R.; Al Asmari, M.; Islam, A.; Kapoor, A.; et al. Middle East respiratory syndrome coronavirus in Bats, Saudi Arabia. *Emerg. Infect. Dis.* **2013**, *19*, 1819–1823. [CrossRef]

52. Monadjem, A.; Soisook, P.; Thong, V.; Kingston, T. Family Hipposideridae (Old World Leaf-nosed Bats). *Handb. Mamm. World* **2019**, *9*, 210–227.

53. Plowright, R.K.; Eby, P.; Hudson, P.J.; Smith, I.L.; Westcott, D.; Bryden, W.L.; Middleton, D.; Reid, P.A.; Mcfarlane, R.A.; Martin, G.; et al. Ecological dynamics of emerging bat virus spillover. *Proc. R. Soc. B Biol. Sci.* **2014**, *282*, 20142124. [CrossRef]

54. Hayman, D.T.S.; Bowen, R.A.; Cryan, P.M.; Mccracken, G.F.; O'Shea, T.J.; Peel, A.J.; Gilbert, A.; Webb, C.T.; Wood, J.L.N. Ecology of Zoonotic Infectious Diseases in Bats: Current Knowledge and Future Directions. *Zoonoses Public Health* **2013**, *60*, 2–21. [CrossRef]

55. Bergner, L.M.; Orton, R.J.; Benavides, J.A.; Becker, D.J.; Tello, C.; Biek, R.; Streicker, D.G. Demographic and environmental drivers of metagenomic viral diversity in vampire bats. *Mol. Ecol.* **2020**, *29*, 26–39. [CrossRef]

56. Menachery, V.D.; Yount, B.L.; Debbink, K.; Agnihothram, S.; Gralinski, L.E.; Plante, J.A.; Graham, R.L.; Scobey, T.; Ge, X.Y.; Donaldson, E.F.; et al. A SARS-like cluster of circulating bat coronaviruses shows potential for human emergence. *Nat. Med.* **2015**, *21*, 1508–1513. [CrossRef] [PubMed]

57. Plowright, R.K.; Parrish, C.R.; Mccallum, H.; Hudson, P.J.; Ko, A.I.; Graham, A.L.; Lloyd-smith, J.O. Pathways to zoonotic spillover. *Nat. Rev. Microbiol.* **2017**, *15*, 502–510. [CrossRef] [PubMed]

58. Parrish, C.R.; Holmes, E.C.; Morens, D.M.; Park, E.; Burke, D.S.; Calisher, C.H.; Laughlin, C.A.; Saif, L.J.; Daszak, P. Cross-Species Virus Transmission and the Emergence of New Epidemic Diseases. *Microbiol. Mol. Biol. Rev.* **2008**, *72*, 457–470. [CrossRef] [PubMed]