Skip to main content

Reproduction of East-African bats may guide risk mitigation for coronavirus spillover

Abstract

Background

Bats provide important ecosystem services; however, current evidence supports that they host several zoonotic viruses, including species of the Coronaviridae family. If bats in close interaction with humans host and shed coronaviruses with zoonotic potential, such as the Severe Acute Respiratory Syndrome virus, spillover may occur. Therefore, strategies aiming to mitigate potential spillover and disease emergence, while supporting the conservation of bats and their important ecological roles are needed. Past research suggests that coronavirus shedding in bats varies seasonally following their reproductive cycle; however, shedding dynamics have been assessed in only a few species, which does not allow for generalization of findings across bat taxa and geographic regions.

Methods

To assess the generalizability of coronavirus shedding seasonality, we sampled hundreds of bats belonging to several species with different life history traits across East Africa at different times of the year. We assessed, via Bayesian modeling, the hypothesis that chiropterans, across species and spatial domains, experience seasonal trends in coronavirus shedding as a function of the reproductive cycle.

Results

We found that, beyond spatial, taxonomic, and life history differences, coronavirus shedding is more expected when pups are becoming independent from the dam and that juvenile bats are prone to shed these viruses.

Conclusions

These findings could guide policy aimed at the prevention of spillover in limited-resource settings, where longitudinal surveillance is not feasible, by identifying high-risk periods for coronavirus shedding. In these periods, contact with bats should be avoided (for example, by impeding or forbidding people access to caves). Our proposed strategy provides an alternative to culling – an ethically questionable practice that may result in higher pathogen levels – and supports the conservation of bats and the delivery of their key ecosystem services.

Background

The order Chiroptera is the second largest order of mammals with more than 1000 identified species [1]. The members of this order, bats, provide important ecosystem services (reviewed in [2, 3]). For example, insectivorous bats can reduce arthropod herbivory [4,5,6], increase agricultural yields [7], reduce the need for insecticides [8], and prevent large financial losses in agriculture [9,10,11]. Plant-visiting chiropterans provide pollination and seed-dispersing services (reviewed in [3]), certain nectivorous bats are pollinators of economically important plants [12], and frugivorous bats can be important for reforestation [13]. Finally, cave-roosting bats produce guano, the main energy source in many cave ecosystems [14, 15], and the mining of this product is an income source in poor communities [16].

However, current evidence supports that bats are a natural host of several disease-causing viruses across the globe, including zoonotic viruses, such as rabies virus (Rhabdoviridae, genus Lyssavirus); Hendra and Nipah viruses (Paramyxoviridae, genus Henipavirus); and Marburg and ebolaviruses (Filoviridae, genus Marburgvirus and Ebolavirus, respectively; [17, 18]). Bats are also hosts of several viruses of the Coronaviridae family [19,20,21,22]. Molecular evidence suggests that the Severe Acute Respiratory Syndrome Betacoronavirus (SARS-CoV betaCoV) and the Middle-East Respiratory Syndrome betaCoV (MERS-CoV) originated from bats [23,24,25,26,27,28,29,30,31]. Both viruses emerged in humans during the past two decades, specifically in China (2002) and Saudi Arabia (2012). The SARS-CoV pandemic included 8096 cases in 27 countries with a ~ 10% case fatality [32], while MERS-CoV has affected 2279 people in 27 countries with a case fatality of ~ 35% [33]. Incidental cases of MERS-CoV are still detected mainly in Saudi Arabia [32], and it is thought that camels are important for human infection [34,35,36].

Today it is known that: i) a high genetic diversity of coronaviruses (CoVs) is present in more than 100 bat species, including viruses related to SARS and MERS CoV [37]; ii) CoVs are prone to move and adapt to new host species [38]; iii) plausibly, all mammalian-adapted CoVs may have originated in bats [20, 21, 38], including a recently emerged highly fatal Alphacoronavirus in piglets [39] and the 229E human CoV [22, 40, 41]; and iv) CoVs found in bats can use human receptors for cell entry [21, 25, 42, 43]. These lines of evidence suggest that future spillover of coronaviruses humans is feasible.

Because CoVs are found in bat species that have adapted to be in close contact with humans, such as the straw-colored fruit bat (Eidolon helvum) and the Brazilian free-tailed bat (Tadarida brasiliensis [44, 45]), high contact “bat-human” interfaces currently exist around the world. If the bats in these interfaces shed CoVs with the ability to infect humans, then opportunities for spillover through direct exposure to feces [37] or the contamination of food are created, as these viruses can remain infectious in the environment for days [46]. Therefore, strategies aiming to mitigate human exposure to CoVs, and thus, the risk of spillover and disease emergence are needed, while supporting the conservation of bats and their important ecological roles.

Longitudinal sampling with specific species has shown that the proportion of bats shedding CoVs varies seasonally [47,48,49,50] and that fecal CoV-RNA loads can also be heterogeneous over time [51, 52]. If exposure through contact with bat feces is a main pathway for zoonotic CoV spillover to humans but shedding of these pathogens is not uniform over time, then mitigation strategies aiming to prevent bat CoV-shed exposure could be targeted temporally, directed especially at high-risk seasons. Such a strategy could guide policy in limited-resource settings where sampling bats for CoV testing is not feasible and it could support an ethically acceptable management to mitigate spillover risk. However, the few species and locations tested to date do not allow for identification of a potential seasonal shedding pattern to responsibly suggest temporal spillover risk management across species and geography. Therefore, assessment of the CoV dynamics in a broader range of bat species that show different life history traits, as well as in diverse geographic and ecological circumstances, could be extremely useful.

To this end, we evaluated the dynamics of CoV shedding in different bat species sampled in several locations in East Africa at different times of the year. This geographical region has been identified as a hotspot of pathogen emergence [53], where CoV host switching events seem to be higher compared to other areas [22], but, to our knowledge, no study on CoV dynamics in bats has been conducted. Specifically, we hypothesized that bat species exhibit seasonal trends in CoV shedding that are associated with the reproductive season. We assessed this hypothesis by fitting Bayesian statistical multivariable models to evaluate whether CoV shedding in bats is positively associated with the time period when pups are becoming independent from the dam. Beyond the inclusion of several species sampled in different countries at different times, we explicitly identified the reproductive events for each species at the time of sampling and also included other traits, such as the aggregation of individuals at the roost, that may be involved in CoV dynamics.

Methods

Sample collection

Samples (rectal swabs and fresh feces) were collected from bats captured in Uganda, Rwanda, and Tanzania (Fig. 1), between September 2011 and April 2014 with permission from local authorities and under the Institutional Animal Care and Use Committee at the University of California, Davis (protocol number: 16048).

Fig. 1
figure 1

Unique locations where samples from bats were obtained

Bats were captured in 36 unique locations between latitudes − 0.9739 and − 10.7506 (Fig. 1). These locations were selected because they represented potentially high-risk interfaces for contact between bats and humans, such as areas of land-use change, human-dwellings, ecotourism sites, markets, and places with potential for occupational exposure [54]. Locations in close proximity (Euclidean distance < 20 km) in which sampling was conducted within the same week were considered a single sampling event. The remaining sampling events that occurred in the same location or spatially close to others but conducted in different weeks were considered independent sampling events. As result, we collected samples from 30 unique sampling events.

All captures were conducted using mist nets set in the early morning or at dawn. Individuals were released after sample collection. Samples were handled as previously described [22, 55]. Each sample was immediately transferred to vials containing Viral Transport Media and NucliSens® Lysis Buffer (bioMérieux, Inc., Marcy-I’Étoile, France), which were maintained in liquid nitrogen until storage in a -80C freezer in each country.

Coronavirus detection

RNA was extracted from all samples, and cDNA was prepared as previously described [22, 55]. Two broadly reactive consensus PCR assays targeting non-overlapping fragments of the orf1ab were used to detect both known and novel CoVs [56, 57]. Amplified products of the expected size were cloned and sequenced as described in [22]. Sequences were analyzed and edited using Geneious Prime 2019.1.3 [58].

A sample was considered positive when at least one PCR assay yielded a sequence that matched coronaviruses in GenBank. Coronavirus sequences were classified as belonging to a specific taxa following previously described methodologies [22].

Bats age, sex, and species identification

Bats were categorized as adults and juveniles based on size, and morphological and behavioral characteristics were observed at sample collection. The sex of the bats was also recorded.

Identification of some bat species can be challenging in the field. For this reason, field-identified species were confirmed by DNA barcoding using the cytochrome b or cytochrome oxidase subunit 1 mitochondrial genes [59]. Obtained sequences were compared against sequences in the GenBank and Barcode of Life databases [60, 61]. When possible or necessary, sequences from both genes were used for species identification. A threshold of 97% nucleotide identity was used to confirm the species; sequences with 95–97% nucleotide identity were assigned a “conferre” (cf.) species status, and sequences below 95% nucleotide identity were either classified to the genus level or as unidentified. Sequences with > 97% nucleotide identity to more than one species for either gene, were classified to the genus level unless they clearly clustered with sequences from other animals in the same geographic area.

If barcoding results for all of the first ten bats tested per species were in agreement with the field identification, we assumed that the field identification for the remaining bats of that species in each country was correct. Otherwise, all of the remaining samples were barcoded to ensure correct speciation.

Bat life history

We recorded when sampled females were pregnant based on abdominal palpation, had attached pups (indicating recent parturition), and were lactating, as well as when juveniles were captured. Therefore, we were able to track pregnancy, lactation, and recent birth pulses. Moreover, we accessed the data in the PanTHERIA [62] and Amniote [63] databases, and we thoroughly reviewed the literature on the biology of the bats species we sampled for latitudes similar to our sampling locations. With the gathering of these information sources we established the timing of the birth pulses, lactation periods, and the weaning of pups for each species. For the details justifying the dates inferred for these three life history events for each species and the corresponding bibliographical references see Additional file 1.

Once the timing of these events was confirmed or inferred, we were able to establish 2 seasons that occur at least once during the year across all observed bat species: i) when juveniles are being weaned and female-pup contact decreases (“Recent weaning” [RW]) and ii) the rest of the year (hereafter “N-RW” for “Not recent weaning”). We chose to evaluate risk of CoV shedding for the first period because past longitudinal studies with microchiropterans in Germany and China found higher CoV-RNA loads approximately 1 month after parturition [51, 52]. Similarly, peaks were found 2 months after the formation of a maternity colony of Myotis macropus [49], which would match a post weaning period for that species. Here, we defined the end of the RW period as 1 month after the last pups were weaned. We assumed that 1 month would provide a reasonable time window for the colony to “clear” the CoV susceptibility status of this period and acquire the susceptibility corresponding to the season(s) when weaning does not occur (N-RW period), if differences actually exist.

Finally, once we determined these two seasons, we categorized each bat sample into one of them depending on the week of the year in which they were taken. Because some species had more than one litter per year, there could be more than one RW period during the year. It is worth noting that we were able to define these periods for those species in our dataset that have synchronized reproduction, whose biology was properly described, and whose taxonomy is generally accepted. When we could not assign a reproductive period to specific bats, this season was imputed (see Methods: Statistical Analysis).

Bat species traits

We characterized specific traits of each bat species studied based on previous scientific literature on pathogen dynamics in bats [51, 64,65,66,67,68,69,70,71]: colony size (small, medium, or large if a typical colony contains one to dozens, hundreds to thousands, or thousands to millions of individuals, respectively); roost type (“closed” if the species has been reported to use caves, mines, roofs, or other confined spaces; or “open” if the bats have been typically reported roosting in the foliage of trees); the aggregation of bats in clusters while roosting (no, yes); and the number of litters per year of the species at equatorial latitudes. References are provided in Additional file 2. We also considered data from PanTHERIA and Amniote [62, 63].

We could not include other species traits, such as the mating strategy (harem or another) and the segregation of females in maternity colonies, because available studies were incomplete or contradictory. We did not include factors, such as multi-species occupancy of the roost, because we did not observe all of the roosts. Further, we did not assess postpartum estrus, as within the study area, it is only known to occur in some Molossid bats [72], of which we only sampled a small number.

Statistical analysis

To statistically estimate the association between the RW and CoV detection we used a Bayesian inference to model the detection of CoVs as a Bernoulli process of the form:

$$ Co{V}_i\sim Bernoulli\left({p}_i\right)\ withi=1,\dots, I $$

where CoVi, the detection of CoV in rectal swabs (1 if detected, 0 otherwise) from the ith bat (sample), is assumed to follow a Bernoulli process parameterized by pi, the probability of CoV detection on the ith bat. This parameter was related to a set of candidate covariates as:

$$ logit\left({p}_i\right)={\alpha}_0+{\beta}_1{X}_{1i}+\gamma {Y}_i+{\rho}_1{C}_{1i}+\dots +{\rho}_l{C}_{li}+{S}_{j(i)}+S{p}_{k(i)}, $$
(1)

with Sj~Normal (0, σSand Spk~Normal (0, σSp).

Here X and Y are binary covariates representing the RW season and juvenile age category. We specifically included these two terms to separate the potential association of the season with CoV detection from the seasonal presence of juveniles. Because it was not feasible to allocate all species in the RW or N-RW seasons based on previous research, we assumed that these unknown reproductive seasons were “Missing at Random” and they were imputed as a function of the latitude at sampling, the day of the year of the sampling event, the number of litters per year of the corresponding species (one litter per year versus more than one litter per year), and the historical precipitation of the month at the sampling event location. The description of the imputation model is provided in Additional file 3.

The terms Sj and Spk represent the sampling event- and the species-specific intercepts, respectively, because we assumed that bats sampled in the same event and bats belonging to the same species were not independent with respect to CoV detection. No sampling event involved the same bat colony in successive RW and N-RW seasons, therefore, we assumed that CoV detection was not temporally correlated within sampling locations.

We constructed the model by adding other covariates one-at-a-time to this working model: the C1...l categorical variables; and they remained in the model if they were judged to confound the relationship between CoV detection and the reproductive seasons or between age and CoV detection (i.e., their inclusion caused meaningful changes in the Posterior Probability Distributions [PPDs] of the specific reproductive season or age coefficients). Finally, C1...l categorical variables could be retained as well if they were marginally associated with CoV detection (the corresponding coefficient PPD did not include zero). To assess the goodness-of-fit of the models, we evaluated the congruence between CoV detection in the data and in the posterior predictive distributions yielded by the models by: i) ages and seasons, ii) age, iii) per season, iv) per age and season, and v) per sampling event.

All models were constructed using “Stan” v. 2.17.0 [73] which was run from “R” v. 3.6.0 [74] through the package “RStan” v. 2.17.3 [75]. Weakly informative priors were assigned for all coefficients: Normal(0, 1.5) for the estimates of α0, β’s, γ, and ρ’s. The σS and σSp had a prior Half-Cauchy (0,5) following previous suggestions [76, 77]. The PPDs were estimated by sampling in parallel from 4 MCMC chains for 4000 iterations following 3000 iteration warm-up for a total of 4000 saved samples for each parameter PPD. Convergence was assessed by the Gelman-Rubin statistic [78] and graphically using trace plots. The code to fit the models is available at https://github.com/dmontecino/East-African-bats-and-CoV-shedding.

Results

Bat samples, age, sex and reproductive seasons

We sampled 753 individuals, all of them aged and successfully identified to belonging to 15 species. Nycteris thebaica, Pipistrellus hesperidus, and Rhinolophus clivosus were assigned the “conferre” status. Hipposiderids were assigned only to the genus level because barcoding did not provide certainty on species identification in line with previously recognized taxonomic difficulties [79,80,81]; however, the biology of the candidate hipposiderid species is similar (Fig. 2; [82, 83]), and we were able to used them for estimation purposes. We excluded Scotophilus viridis (n = 6) from the analysis because their reproductive traits are unknown, and this species also has taxonomic difficulties for species identification [84]. Therefore, these six individuals were removed, and 747 bats were included in the study.

Fig. 2
figure 2

Summary of the inferred reproductive periods of the bat species sampled. The red squares show the week of a typical year each bat species was sampled. The yellow, green, and blue polygons show the assigned extent of the birth pulse (yellow), lactation (green), and mating-pregnancy periods (blue). The grey polygons represent the inferred “Recent weaning” period(s) per bat species. In species with unknown lactation length (N. nana and T. persicus), the light green polygons represent the likely overextension of this period not including any other bat in the “Recent weaning” period. The question marks show the period we could not infer the corresponding reproductive activities

We were able to infer the reproductive season for all bats except for Lissonycteris angolensis, Rhinolophus cf. clivosus, and Pipistrellus cf. hesperidus (n = 117, 104 adults and 13 juveniles) that had limited available knowledge on biology and reproductive season. These reproductive seasons were imputed as explained above. This imputation process did not substantively affect the proportion of bats in each reproductive period or the crude CoV detection per age (Additional file 3: Figure S3.2). Moreover, we had partial data for Neoromicia nana and Triaenops persicus lactation period and we assigned one that is likely longer than it would be in reality based on the other species. Even under this overextended period, none of these bats (N. nana and T. persicus) could have being sampled during the RW season, so this knowledge gap did not risk a misclassification (Fig. 2).

In the end, 274 and 356 bats were allocated in the N-RW and RW periods, respectively. The distribution of bats across the reproductive periods per age and sampling event was heterogeneous as expected due to the opportunistically nature of sampling. Indeed, 233 and 273 adult bats were in the N-RW and RW periods, respectively, while 41 and 83 juvenile bats were sampled in these seasons, respectively.

Coronavirus detection and identification

In total, 30.79% bats were positive for CoVs (n = 230). Within the subset of adult bats, 26.22% were positive (n = 160) while 51.09% of the juvenile bats were positive (n = 70). The detection of CoV shedding was variable across seasons and bat species, as well as across sampling events (0, 100, 9.69, and 25.84% for the minimum, maximum, median, and mean detection, respectively). The CoVs found per species are shown in Additional file 4: Table S4.

Species traits

A summary of the roosting and reproductive traits of the bat species sampled is provided in Table 1. All bats except N. nana (n = 9) and E. helvum (n = 315) roosted in “closed” structures, such as caves, abandoned mines, and roofs. Within the group of bats using “open” structures, E. helvum was the only species with CoV positive individuals. Therefore, we did not use this covariate to assess a potential association with CoV shedding.

Table 1 Summary of traits by sampled chiropteran species

Statistical analysis

The models showed adequate sampling. The 4 Markov Chain Monte Carlo chains converged graphically, whilst all Gelman-Rubin statistics were < 1.004. The selected model had a number of effective samples for each coefficient of at least 1636. The data were properly fitted, as well (Additional file 5: Figure S5.1), although some predictions lacked precision. The PPDs of the fixed coefficients are shown in Additional file 5: Figure S5.2.

The selected model to assess periodic differences in CoV shedding included season and age, species-specific intercepts, and sampling event-specific intercepts. Beyond the species-specific terms, we included a binary categorical covariate equal to 1 for E. helvum and T. persicus and 0 otherwise. We incorporated this term because CoV detection in these species was remarkably higher than the other species. As expected, this fixed effect was correlated with the corresponding species-specific intercepts (the remaining correlations were all low); however, we decided to keep it to assess if the main findings hold even when accounting for the bat species with highest detection. The corresponding means, standard deviations, and 90% HPDI are shown in Table 2.

Table 2 Summary of the posterior probability distributions of the fixed-effects coefficients of the selected model

The coefficients’ PPDs from the selected model indicate an association between age and CoV shedding, with juveniles presenting 1.26–2.94 times higher odds to shed compared to adult bats (90% HPDI). The coefficients’ PPDs also point to an association between the reproductive season and CoV shedding as well, with an estimated odds 1.71–16.00 times higher to shed during the period when pups are being compared to other seasons (90% HPDI). The proportions of CoV shedders estimated by reversing the 90% HPDI of the logits were: 0.02–0.22, 0.09–0.59, 0.01–0.13, and 0.05–0.42 for juveniles during the “N-RW” and “RW” periods, and adults during the “N-RW” and “RW” periods, respectively (90% HPDI). These values refer to bats not belonging to the species E. helvum or T. persicus. Finally, the predicted CoV detections, based on the posterior predictive distributions, were 0.01–0.18, 0.06–0.54, 0.00–0.05, and 0.04–0.36, for these same groups (90% HPDI; Fig. 3, left). In practical terms, these last values imply that juveniles during the “RW” period are, on average, 3.34 times more likely to be detected shedding CoVs compared to juveniles in the “N-RW” period. Adults during the “RW” period are, on average, 3.93 times more likely to be detected shedding CoVs compared to adults in the “N-RW” period. In both seasons, juveniles are, on average, more likely to shed CoVs, than adults.

Fig. 3
figure 3

a Modeled proportion of chiropterans shedding coronavirus viral particles in two reproductive periods inferred for these bats by age class. b The estimated species-specific intercept coefficients. The lighter colors represent the range where the 90% of the estimated detection values are concentrated (the High Posterior Density Interval). The darker colors show the 51% High Posterior Density Interval

The selected model suggests a higher odds of CoV detection in E. helvum and T. persicus compared to other species. The species-specific intercept terms, once the E. helvum - T. persicus effect is included, suggest no further differences in terms of CoV shedding (Fig. 3, right); however, the estimates are not precise. The sampling event-specific random intercepts suggest that a few specific locations could show differential CoV shedding but that most of them do not explain further variation (Additional file 6).

Discussion

If coronavirus shedding by bats follows temporal patterns that are generalizable across species and locations, then mitigation strategies targeting the prevention of human exposure and potential spillover could be directed toward high-risk periods, through mechanisms that can also support bat-human coexistence and the provision of bat ecosystem services. Previous research has focused on viral identification in specific locations and in few species [47,48,49,50], resulting in a limited representation of viral dynamics in association with few ecological settings, biological traits, and reproductive strategies. Additionally, few studies of coronavirus shedding patterns have employed statistical models and, in consequence, the potential complex web of factors and causal relationships that may determine this process has not been fully explored. Here, we aimed to address these issues by statistically modeling coronavirus RNA detection in several bat species, captured at different times and locations in East Africa and involving different ecological contexts and life histories. Using data from several hundreds of bats, we found that, beyond spatial, taxonomic, and life history differences; i) the odds of coronavirus shedding is higher during the period when pups are being weaned (up to a month after the lactation period is over), and ii) juvenile bats have higher odds to shed these viruses. Moreover, the ratios of predicted detections per bat category (age and reproductive season) suggest that juvenile bats during the recent weaning period have relatively higher shedding compared to bats out of this period no matter their age. Caution must be taken with these ratios because we used a logit link and our data had high proportion of CoV shedding in specific groups. However, our results are consistent, and they are in agreement with previous research conducted in a restricted number of species and locations.

Similar seasonality of coronavirus shedding has been observed in Germany, Australia, Thailand, China, and Ghana (West Africa). In the specific species involved in these previous studies, higher coronavirus shedding and viral loads were detected weeks after the birth pulse [47,48,49,50,51,52, 85]. Further, and consistent with our results, detection of higher levels of coronavirus in juveniles has been reported in micro- as well megabats from Africa, Asia, Europe, and North and South America [22, 48, 50, 86,87,88,89].

It has been proposed that the increased detection of coronaviruses after the birth pulse is attributable to the waning of passively-received maternal antibodies in juveniles [51]. This idea has been frequently cited; however, we are not aware of any longitudinal age-specific coronavirus seroprevalence study in bats. Such studies are important to understand the drivers of pathogen persistence and spillover risk, and in consequence, to ethically manage and prevent bat pathogen exposure. Nevertheless, this kind of research is difficult to conduct due to logistical challenges, our questionable ability to obtain statistically-representative samples across age groups, cross-reactivity of serological assays, and the difficulties to differentiate serodynamics derived from closed-population processes from those caused by migratory movements. Although extrapolations for antibody dynamics across viruses and species are not simple [90], bat serodynamics for Hendra virus are congruent with the increased detection of coronaviruses after the birth pulse. Pups passively receive maternal Hendra virus antibodies which decline after the first month of age up to 6 months of age [91,92,93,94]. This decline would lead to a period in which young bats tend to be more susceptible to infection, become infected, and then shed virus. Consistently, coronavirus shedding peaked weeks after the birth pulse in a German and a Chinese species [51, 52] and immunologically naïve bats shed higher coronavirus loads [30]. Over time, as young bats clear Hendra virus infection, they become seropositive again [91,92,93,94]. Concordantly, capture-mark-recapture studies support the clearance of coronaviruses in infected bats [47, 95], which would become seropositive. However, young bats may not reach adult seroprevalence levels until they are older than a year, as occurs with Hendra virus [91, 92, 94]. Therefore, the population of juvenile bats would remain comparatively more susceptible to viral infection and shedding beyond the period immediately after weaning. Age-specific henipavirus seroprevalence in African E. helvum is in agreement with the serodynamics described for Hendra virus [90].

Additionally, coronavirus transmission may be favored by high colony density created by the birth pulse, as previously proposed [85, 95], and then the seasonal influx of susceptible juveniles could accelerate viral spread across the entire colony, including adult bats. Indeed, adult Myotis macropus in an Australian colony showed a peak of coronavirus detection after the birth pulse [85]. The peak of coronavirus detection for two Hipposideros species and Nycteris cf. gambiensis sampled in Ghana occurred during the months that encompassed the birth pulse and nursing after accounting for the age of the sampled individuals [50]. Higher coronavirus infection has also been reported in lactating females [85, 88], which overlap with the period of pup weaning and decay of maternally-derived immunity; however, the opposite has also been found [48, 50, 85].

In practical terms, public health managers could anticipate high risk periods for coronavirus shedding to target interventions. Assuming that higher spillover risk is a function of higher viral shedding [67] and that all coronaviruses with zoonotic potential behave ecologically similarly to coronaviruses detected in this study, managers could target the prevention of human-bat direct (consumption) or indirect (bat droppings) contact specifically during the high-risk season: around and just after weaning, the timing of observable juveniles or individuals smaller than adults. For the species and interfaces defined herein, those management periods have now been determined (Fig. 2). For others, direct observation of bats at high-risk transmission interfaces could be used to identify time periods when non-adult sized bats are present. However, observation of dependent pups is not always easy [96]. Of course, for specific species, birth pulses and lactation seasons could also be used to more precisely establish high risk periods similar to the methods we used here, including a combination of direct observation, reports from previous literature, and consultation with knowledgeable bat biologists. Our proposed risk-driven strategy i) is evidence-based, as it builds upon coronavirus shedding patterns observed across several chiropteran species present around the world; ii) does not require the advanced laboratory capacity often lacking in resource-restricted settings where intense bat-human interfaces usually occur; iii) is a good alternative to the ideal but expensive and resource-intensive longitudinal surveys; and iv) it may prevent the exposure to viruses belonging to other taxa whose observed bat shedding dynamics resemble our findings for coronaviruses (e.g. paramyxoviruses [97]),

The lower coronavirus detections in African emballonurids (C. afra and T. mauritianus) and the higher detections in E. helvum, African hipposiderids (Hipposideros sp. and T. persicus), and R. cf. clivosus that we found are consistent with previous reports [22, 44, 50, 98] and should be considered by managers when providing risk-based spillover prevention strategies. Moreover, SARS-like coronaviruses in Africa have been found in hipposiderid, rhinolophid, and molossid bats [22, 57, 98], and MERS-related coronaviruses have been found in vespertillionid bats [22, 28, 29]. Therefore, it seems reasonable to prioritize the identification of birth pulses and lactation seasons, and thus determine high-risk periods of coronavirus shedding, for these bat families. Interestingly, E. helvum roost in tight clusters that can contain hundreds of individuals [99], similar to T. persicus. On the other hand, emballonurids, showing the lowest crude detection levels across families (represented by Coleura afra and Taphozous mauritanus here) tend not to cluster while roosting [100, 101]. We did find an association between coronavirus shedding and whether the species typically aggregate in clusters while roosting when the variable “E. helvumT. persicus” was not considered, but we chose a different model not including this term because we did not directly assess bat roosts and our categorization may oversimplify the continuum from mostly solitary roosting (e.g., Neromocia nana) to common tight aggregations of bats (e.g., Mops condylurus). Using this categorization could be misleading, as some species differentially cluster while roosting depending on temperature, colony size, colony type (e.g., maternity colony versus not a maternity colony), and season [72]. Future studies should consider the roosting habits of bats, as this trait could further support risk-based management to prevent or reduce human exposure.

The risk-driven strategy we propose provides a contact-reduction alternative that is ethically favorable compared to often-employed measures, such as culling or other reactive measures, that ensue when the public becomes aware of a health threat without a suggested practical option to reduce their risk for exposure. In addition to ethical concerns and being logistically difficult and expensive, culling has failed to reduce disease in wild populations and can result in even higher pathogen levels. For example, “badger culling can make no meaningful contribution to cattle tuberculosis control in Britain” [102]. A culling program to reduce Echinococcus multilocularis prevalence in red foxes (Vulpes vulpes) resulted in an increase of infection [103]. This strategy has also failed to control rabies in canids around the globe [104]. Similar results have been observed in bats. Culling failed to reduce rabies seroprevalence in Desmodus rotundus in Perú and could have increased the levels of exposure to the virus [105, 106]. In Argentina, the extermination of bats changed the direction of spread of rabies in livestock but did not prevent its advancement [107]. In Uganda, miners exterminated a colony of Rousettus aegyptiacus bats after an outbreak of Marburg virus in 2007 that involved 4 miners in close contact with these bats. Five years later, a new outbreak occurred in miners from the same mine. The second time, Marburg virus RNA was detected in a higher proportion in the R. aegyptiacus that recolonized the mine (13.3%, n = 400; [108]) compared to RNA detection before culling in this cave (5.1%, n = 611; [109]) and other caves in Uganda (2.5%, n = 1622; [67]) and Gabon (4.8%, n = 187; [110]) where culling has never been reported.

Culling can also cause demographic changes, leading to a higher proportion of juvenile individuals. This change may occur because of a disproportionate cull of older individuals; the potential increase in survival of pups at lower population densities, followed by higher recruitment of juvenile females into reproductive age [111]; the hypothetical increase of young dispersers immigrating from neighboring colonies into culled, less dense, and better resourced colonies [105, 112]; or by causing compensatory reproduction [113]. This last possibility may have not been studied in bats but seems unlikely due to their high conception rates and usual litter size of one. Examples of younger populations after culling have been reported in the red deer (Cervus elaphus), racoon (Procyon lotor), American mink (Mustela vison), and Australian brushtail possums (Trichosurus vulpecula), among others [114,115,116,117]. As our results and past research consistently show higher viral shedding and detection in young individuals, activities leading to a younger bat population are not advisable for viral spillover management. Similar results are expected when fruit bats are culled based on being categorized as “agricultural pests”; therefore, this kind of management may create higher risk of viral exposure to the human population.

Virological, ecological, and epidemiological research on bats over the last 15 years has helped to identify chiropterans as hosts of zoonotic viruses and to document that human-driven environmental change, human behavior, and human-to-human transmission are the key drivers for the creation of bat-human interfaces, spillover, and epidemics of emergent viruses, respectively [118,119,120]. In the context of the current biodiversity and bat conservation crisis [121, 122], we must not omit these facts when attempting to effectively, and responsibly frame and communicate disease risks associated with bats. Realistic, data-based risk communication is of paramount importance to avoid framing bats as a threat to humans and to support bat conservation given their important ecological roles [123, 124]. With this background, it seems a proper time for the scientific community studying “bat-associated” viruses to move the conversation from bat spillover risk assessments to the planning of pro-biodiversity and subsequently pro-ecosystem strategies aiming to mitigate spillover risk. Science is valued not only for the diagnosis of problems but because it finds solutions to them. Here, we have attempted to aid the progress of scientific and management dialogue by proposing, not only a management strategy to limit potential coronavirus spillover, but one that is context- and logistically-grounded and pro-conservation, promoting the delivery of the key ecosystem services provided by bats.

Conclusions

Data from hundreds of bats collected in East-Africa show that coronavirus shedding is expected to be more frequent when pups are becoming independent from the dam, independently of the age of the bats, their species, their location, and their life histories; however, the odds of shedding do differ by species. These results can guide temporal-based mitigation strategies to prevent bat-associated coronavirus exposure using non-lethal methods in limited-resource settings, where longitudinal surveillance is not feasible, by identifying high-risk periods for coronavirus shedding when contact with bats should be avoided.

Availability of data and materials

The datasets used and/or analyzed are available at https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.9917762.

Abbreviations

betaCoV:

Betacoronavirus

cf:

Conferre

CoV:

Coronavirus

MERS-CoV:

Middle-East Respiratory Syndrome betaCoV

N-RW:

Not recent weaning

PPD:

Posterior probability distribution

RW:

Recent weaning

SARS-CoV:

Severe Acute Respiratory Syndrome betaCoV

References

  1. Mammal Diversity Database. 2019. www.mammaldiversity.org. American Society of Mammalogists. Accessed 20 Mar 2019.

  2. Ghanem SJ, Voigt CC. Increasing awareness of ecosystem services provided by bats. In: Brockmann HJ, Roper TJ, Naguib M, Mitani JC, Simmons LW, editors. Advances in the study of behavior. San Diego: Elsevier; 2012. p. 279–302.

    Google Scholar 

  3. Kunz TH, Braun de Torrez E, Bauer D, Lobova T, Fleming TH. Ecosystem services provided by bats. Ann N Y Acad Sci. 2011;1223:1–38.

    Article  PubMed  Google Scholar 

  4. Kalka MB, Smith AR, Kalko EKV. Bats limit arthropods and herbivory in a tropical forest. Science. 2008;320:71.

    Article  CAS  PubMed  Google Scholar 

  5. Williams-Guillén K, Perfecto I, Vandermeer J. Bats limit insects in a neotropical agroforestry system. Science. 2008;320:70.

    Article  PubMed  CAS  Google Scholar 

  6. Böhm SM, Wells K, Kalko EKV. Top-down control of herbivory by birds and bats in the canopy of temperate broad-leaved oaks (Quercus robur). PLoS One. 2011;6:e17857.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Maas B, Clough Y, Tscharntke T. Bats and birds increase crop yield in tropical agroforestry landscapes. Ecol Lett. 2013;16:1480–7.

    Article  PubMed  Google Scholar 

  8. Federico P, Hallam TG, McCracken GF, Purucker ST, Grant WE, Correa-Sandoval AN, et al. Brazilian free-tailed bats as insect pest regulators in transgenic and conventional cotton crops. Ecol Appl. 2008;18:826–37.

    Article  PubMed  Google Scholar 

  9. Cleveland CJ, Betke M, Federico P, Frank JD, Hallam TG, Horn J, et al. Economic value of the pest control service provided by Brazilian free-tailed bats in south-Central Texas. Front Ecol Environ. 2006;4:238–43.

    Article  Google Scholar 

  10. Boyles JG, Cryan PM, McCracken GF, Kunz TH. Conservation. Economic importance of bats in agriculture. Science. 2011;332:41–2.

    Article  PubMed  Google Scholar 

  11. Taylor PJ, Grass I, Alberts AJ, Joubert E, Tscharntke T. Economic value of bat predation services--a review and new estimates from macadamia orchards. Ecosyst Serv. 2018;30:372–81.

    Article  Google Scholar 

  12. Ducummon SL. Ecological and economic importance of bats. 2000. https://www.cbd.int/financial/values/g-ecobats.pdf. Accessed 7 Jul 2019.

    Google Scholar 

  13. Muscarella R, Fleming TH. The role of frugivorous bats in tropical forest succession. Biol Rev Camb Philos Soc. 2007 Nov;82:573–90.

    Article  PubMed  Google Scholar 

  14. Fenolio DB, Graening GO, Collier BA, Stout JF. Coprophagy in a cave-adapted salamander; the importance of bat guano examined through nutritional and stable isotope analyses. Proc Biol Sci. 2006;273:439–43.

    Article  CAS  PubMed  Google Scholar 

  15. Gnaspini P, Trajano E. Guano communities in tropical caves. In: Wilkins H, Culver DC, Humphreys WF, editors. Ecosystems of the world subterranean ecosystems. Amsterdam: Elsevier; 2000. p. 251–68.

    Google Scholar 

  16. Food and Agriculture Organization of the United Nations. In: Newman SH, Field H, Epstein J, de Jong C, editors. Investigating the role of bats in emerging zoonoses: balancing ecology, conservation and public health interest. Rome: FAO; 2011.

    Google Scholar 

  17. Hayman DTS. Bats as viral reservoirs. Annu Rev Virol. 2016;3:77–99.

    Article  CAS  PubMed  Google Scholar 

  18. Calisher CH, Childs JE, Field HE, Holmes KV, Schountz T. Bats: important reservoir hosts of emerging viruses. Clin Microbiol Rev. 2006;19:531–45.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Woo PCY, Lau SKP, Li KSM, Poon RWS, Wong BHL, Tsoi H-W, et al. Molecular diversity of coronaviruses in bats. Virology. 2006;351:180–7.

    Article  CAS  PubMed  Google Scholar 

  20. Vijaykrishna D, Smith GJD, Zhang JX, Peiris JSM, Chen H, Guan Y. Evolutionary insights into the ecology of coronaviruses. J Virol. 2007;81:4012–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Drexler JF, Corman VM, Drosten C. Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS. Antivir Res. 2014;101:45–56.

    Article  CAS  PubMed  Google Scholar 

  22. Anthony SJ, Johnson CK, Greig DJ, Kramer S, Che X, Wells H, et al. Global patterns in coronavirus diversity. Virus Evol. 2017;3:vex012.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Li W, Shi Z, Yu M, Ren W, Smith C, Epstein JH, et al. Bats are natural reservoirs of SARS-like coronaviruses. Science. 2005;310:676–9.

    Article  CAS  PubMed  Google Scholar 

  24. Anthony SJ, Gilardi K, Menachery VD, Goldstein T, Ssebide B, Mbabazi R, et al. Further evidence for bats as the evolutionary source of middle east respiratory syndrome coronavirus. MBio. 2017;8,e00373-17.

  25. Ge X-Y, Li J-L, Yang X-L, Chmura AA, Zhu G, Epstein JH, et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature. 2013;503:535–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Yuan J, Hon C-C, Li Y, Wang D, Xu G, Zhang H, et al. Intraspecies diversity of SARS-like coronaviruses in Rhinolophus sinicus and its implications for the origin of SARS coronaviruses in humans. J Gen Virol. 2010;91:1058–62.

    Article  CAS  PubMed  Google Scholar 

  27. Hu B, Ge X, Wang L-F, Shi Z. Bat origin of human coronaviruses. Virol J. 2015;12:221.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Ithete NL, Stoffberg S, Corman VM, Cottontail VM, Richards LR, Schoeman MC, et al. Close relative of human Middle East respiratory syndrome coronavirus in bat, South Africa. Emerg Infect Dis. 2013;19:1697–9.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Corman VM, Ithete NL, Richards LR, Corrie Schoeman M, Preiser W, Drosten C, et al. Rooting the phylogenetic tree of MERS-coronavirus by characterization of a conspecific virus from an African bat. J Virol. 2014;88:11297–303.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Lau SKP, Woo PCY, Li KSM, Huang Y, Tsoi H-W, Wong BHL, et al. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc Natl Acad Sci. 2005;102:14040–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Memish ZA, Mishra N, Olival KJ, Fagbo SF, Kapoor V, Epstein JH, et al. Middle East respiratory syndrome coronavirus in bats, Saudi Arabia. Emerg Infect Dis. 2013;19:1819–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. World Health Organization. Summary of probable SARS cases with onset of illness from 1 November 2002 to 31 July 2003. 2004. Emergencies preparedness, response. http://www.who.int/csr/sars/country/table2004_04_21/en/. Accessed 14 Aug 2018.

    Google Scholar 

  33. World Health Organization. Middle East respiratory syndrome coronavirus. https://www.who.int/emergencies/mers-cov/en/. Accessed 14 Oct 2018.

  34. Haagmans BL, Al Dhahiry SHS, Reusken CBEM, Raj VS, Galiano M, Myers R, et al. Middle East respiratory syndrome coronavirus in dromedary camels: an outbreak investigation. Lancet Infect Dis. 2014;14:140–5.

    Article  CAS  PubMed  Google Scholar 

  35. Azhar EI, Hashem AM, El-Kafrawy SA, Sohrab SS, Aburizaiza AS, Farraj SA, et al. Detection of the Middle East respiratory syndrome coronavirus genome in an air sample originating from a camel barn owned by an infected patient. MBio. 2014;5:e01450–14.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Azhar EI, El-Kafrawy SA, Farraj SA, Hassan AM, Al-Saeed MS, Hashem AM, et al. Evidence for camel-to-human transmission of MERS coronavirus. N Engl J Med. 2014;370:2499–505.

    Article  CAS  PubMed  Google Scholar 

  37. Chen L, Liu B, Yang J, Jin Q. DBatVir: the database of bat-associated viruses. Database. 2014;2014:bau021.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Woo PCY, Lau SKP, Huang Y, Yuen K-Y. Coronavirus diversity, phylogeny and interspecies jumping. Exp Biol Med. 2009;234:1117–27.

    Article  CAS  Google Scholar 

  39. Pan Y, Tian X, Qin P, Wang B, Zhao P, Yang Y-L, et al. Discovery of a novel swine enteric alphacoronavirus (SeACoV) in southern China. Vet Microbiol. 2017;211:15–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Corman VM, Baldwin HJ, Tateno AF, Zerbinati RM, Annan A, Owusu M, et al. Evidence for an ancestral association of human coronavirus 229E with bats. J Virol. 2015;89:11858–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Pfefferle S, Oppong S, Drexler JF, Gloza-Rausch F, Ipsen A, Seebens A, et al. Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats, Ghana. Emerg Infect Dis. 2009;15:1377–84.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Yang X-L, Hu B, Wang B, Wang M-N, Zhang Q, Zhang W, et al. Isolation and characterization of a novel bat coronavirus closely related to the direct progenitor of SARS coronavirus. J Virol. 2015;90:3253–6.

    Article  PubMed  CAS  Google Scholar 

  43. Luo C-M, Wang N, Yang X-L, Liu H-Z, Zhang W, Li B, et al. Discovery of novel bat coronaviruses in south China that use the same receptor as MERS coronavirus. J Virol. 2018;92:JVI.00116–8.

    Google Scholar 

  44. Tao Y, Shi M, Chommanard C, Queen K, Zhang J, Markotter W, et al. Surveillance of bat coronaviruses in Kenya identifies relatives of human coronaviruses NL63 and 229E and their recombination history. J Virol. 2017;91:JVI.01953–16.

    Google Scholar 

  45. de Sales Lima FE, Campos FS, Kunert Filho HC, Batista HB d CR, Júnior PC, Cibulski SP, et al. Detection of Alphacoronavirus in velvety free-tailed bats (Molossus molossus) and Brazilian free-tailed bats (Tadarida brasiliensis) from urban area of southern Brazil. Virus Genes. 2013;47:164–7.

    Article  CAS  Google Scholar 

  46. Geller C, Varbanov M, Duval RE. Human coronaviruses: insights into environmental resistance and its influence on the development of new antiseptic strategies. Viruses. 2012;4:3044–68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Lau SKP, Li KSM, Huang Y, Shek C-T, Tse H, Wang M, et al. Ecoepidemiology and complete genome comparison of different strains of severe acute respiratory syndrome-related Rhinolophus bat coronavirus in China reveal bats as a reservoir for acute, self-limiting infection that allows recombination events. J Virol. 2010;84:2808–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Wacharapluesadee S, Duengkae P, Chaiyes A, Kaewpom T, Rodpan A, Yingsakmongkon S, et al. Longitudinal study of age-specific pattern of coronavirus infection in Lyle’s flying fox (Pteropus lylei) in Thailand. Virol J. 2018;15:38.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Smith C. Persistent or long-term coronavirus infection in Australian bats. Microbiol Aust. 2017;38:8–11.

    Google Scholar 

  50. Baldwin HJ. Epidemiology and ecology of virus and host: bats and coronaviruses in Ghana, West Africa, Doctor of philosophy dissertation. New South Wales: Macquarie University; 2015.

    Google Scholar 

  51. Drexler JF, Corman VM, Wegner T, Tateno AF, Zerbinati RM, Gloza-Rausch F, et al. Amplification of emerging viruses in a bat colony. Emerg Infect Dis. 2011;17:449–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Wang M-N, Zhang W, Gao Y-T, Hu B, Ge X-Y, Yang X-L, et al. Longitudinal surveillance of SARS-like coronaviruses in bats by quantitative real-time PCR. Virol Sin. 2016;31:78–80.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, et al. Global trends in emerging infectious diseases. Nature. 2008;451:990–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Greger M. The human/animal interface: emergence and resurgence of zoonotic infectious diseases. Crit Rev Microbiol. 2007;33:243–99.

    Article  PubMed  Google Scholar 

  55. Anthony SJ, Epstein JH, Murray KA, Navarrete-Macias I, Zambrana-Torrelio CM, Solovyov A, et al. A Strategy to estimate unknown viral diversity in mammals. mBio. 2013;4:mbio.00598–13.

    Article  CAS  Google Scholar 

  56. Watanabe S, Masangkay JS, Nagata N, Morikawa S, Mizutani T, Fukushi S, et al. Bat coronaviruses and experimental infection of bats, the Philippines. Emerg Infect Dis. 2010;16:1217–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Quan P-L, Firth C, Street C, Henriquez JA, Petrosov A, Tashmukhamedova A, et al. Identification of a Severe Acute Respiratory Syndrome Coronavirus-like virus in a leaf-nosed bat in Nigeria. mBio. 2010;1:mBio.00208–10.

    Article  CAS  Google Scholar 

  58. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Townzen JS, Brower AVZ, Judd DD. Identification of mosquito bloodmeals using mitochondrial cytochrome oxidase subunit I and cytochrome b gene sequences. Med Vet Entomol. 2008;22:386–93.

    Article  CAS  PubMed  Google Scholar 

  60. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Ostell J, Pruitt KD, et al. GenBank. Nucleic Acids Res. 2018;46:D41–7.

    Article  CAS  PubMed  Google Scholar 

  61. Ratnasingham S, Hebert PDN. The barcode of life data system. Mol Ecol Notes. 2007;7:355–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Jones KE, Bielby J, Cardillo M, Fritz SA, O’Dell J, Orme CDL, et al. PanTHERIA: a species-level database of life history, ecology, and geography of extant and recently extinct mammals. Ecology. 2009;90:2648.

    Article  Google Scholar 

  63. Myhrvold NP, Baldridge E, Chan B, Sivam D, Freeman DL, Ernest SKM. An amniote life-history database to perform comparative analyses with birds, mammals, and reptiles. Ecology. 2015;96(11):3109.

    Article  Google Scholar 

  64. Luis AD, O’Shea TJ, Hayman DTS, Wood JLN, Cunningham AA, Gilbert AT, et al. Network analysis of host-virus communities in bats and rodents reveals determinants of cross-species transmission. Ecol Lett. 2015;18:1153–62.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Mühldorfer K, Speck S, Kurth A, Lesnik R, Freuling C, Müller T, et al. Diseases and causes of death in European bats: dynamics in disease susceptibility and infection rates. PLoS One. 2011;6:e29773.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Hayman DTS, Suu-Ire R, Breed AC, McEachern JA, Wang L, Wood JLN, et al. Evidence of henipavirus infection in west African fruit bats. PLoS One. 2008;3:e2739.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  67. Amman BR, Carroll SA, Reed ZD, Sealy TK, Balinandi S, Swanepoel R, et al. Seasonal pulses of Marburg virus circulation in juvenile Rousettus aegyptiacus bats coincide with periods of increased risk of human infection. PLoS Pathog. 2012;8:e1002877.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Kolodny O, Weinberg M, Reshef L, Harten L, Hefetz A, Gophna U, et al. Coordinated change at the colony level in fruit bat fur microbiomes through time. Nat Ecol Evol. 2018. https://0-doi-org.brum.beds.ac.uk/10.1038/s41559-018-0731-z.

    Article  Google Scholar 

  69. Willoughby AR, Phelps KL, PREDICT Consortium, Olival KJ. A Comparative analysis of viral richness and viral sharing in cave-roosting bats. Diversity. 2017;9:35.

    Article  Google Scholar 

  70. Hoyt JR, Langwig KE, White JP, Kaarakka HM, Redell JA, Kurta A, et al. Cryptic connections illuminate pathogen transmission within community networks. Nature. 2018;563:710–3.

    Article  CAS  PubMed  Google Scholar 

  71. van Schaik J, Kerth G. Host social organization and mating system shape parasite transmission opportunities in three European bat species. Parasitol Res. 2017;116:589–99.

    Article  PubMed  Google Scholar 

  72. Happold M, Happold D, editors. Mammals of Africa volume IV-hedgehogs, shrews and bats. London: Bloomsbury Publishing; 2013.

    Google Scholar 

  73. Carpenter B, Gelman A, Hoffman MD, Lee D, Goodrich B, Betancourt M, et al. Stan: A probabilistic programming language. J Stat Softw. 2017;76(1) Available from: https://www.osti.gov/biblio/1430202.

  74. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2019. https://www.R-project.org/.

  75. Stan Development Team. RStan: the R interface to Stan. 2018.

    Google Scholar 

  76. Gelman A. Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper). Bayesian Anal. 2006;1:515–34.

    Article  Google Scholar 

  77. Polson NG, Scott JG. On the half-Cauchy prior for a global scale parameter. Bayesian Anal. 2012;7:887–902.

    Article  Google Scholar 

  78. Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;7:457–72.

    Article  Google Scholar 

  79. Simmons NB. Order Chiroptera. In: Wilson DE, Reeder DM, editors. Mammal species of the world: a taxonomic and geographic reference; 2005. p. 312–529.

    Google Scholar 

  80. Vallo P, Benda P, Martínková N, Kauch P, Kalko EKV, Čeý J, et al. Morphologically uniform bats Hipposideros aff. ruber (Hipposideridae) exhibit high mitochondrial genetic diversity in southeastern Senegal. Acta Chiropt. 2011;13:79–88.

    Article  Google Scholar 

  81. Vallo P, Guillén-Servent A, Benda P, Pires DB, Koubek P. Variation of mitochondrial DNA in the Hipposideros caffer complex (Chiroptera: Hipposideridae) and its taxonomic implications. Acta Chiropt. 2008;10:193–206.

    Article  Google Scholar 

  82. Happold M. Hipposideros ruber Noack’s leaf-nosed bat. In: Happold M, David H, editors. Mammals of Africa volume IV: bats and shrews; 2013. p. 393–5.

    Google Scholar 

  83. Happold M. Hipposideros cafer Sundevall’s Roundleaf bat. In: Happold M, David H, editors. Mammals of Africa volume IV: bats and shrews; 2013. p. 375–8.

    Google Scholar 

  84. Schweiger BR. Elucidating the species limits and range boundaries of the African yellow house bats, genus Scotophilus. Senior thesis. Lake Forest college, Illinois, US; 2017.

  85. Smith CS. Australian bat coronaviruses, Philosophical doctor thesis. Queensland: University of Queensland; 2014.

    Google Scholar 

  86. Gloza-Rausch F, Ipsen A, Seebens A, Göttsche M, Panning M, Drexler JF, et al. Detection and prevalence patterns of group I coronaviruses in bats, northern Germany. Emerg Infect Dis. 2008;14:626–31.

    Article  PubMed  PubMed Central  Google Scholar 

  87. Rihtarič D, Hostnik P, Steyer A, Grom J, Toplak I. Identification of SARS-like coronaviruses in horseshoe bats (Rhinolophus hipposideros) in Slovenia. Arch Virol. 2010;155:507–14.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  88. Osborne C, Cryan PM, O’Shea TJ, Oko LM, Ndaluka C, Calisher CH, et al. Alphacoronaviruses in New World bats: prevalence, persistence, phylogeny, and potential for interaction with humans. PLoS One. 2011;6:e19156.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Annan A, Baldwin HJ, Corman VM, Klose SM, Owusu M, Nkrumah EE, et al. Human betacoronavirus 2c EMC/2012-related viruses in bats, Ghana and Europe. Emerg Infect Dis. 2013;19:456–9.

    Article  PubMed  PubMed Central  Google Scholar 

  90. Peel AJ, Baker KS, Hayman DTS, Broder CC, Cunningham AA, Fooks AR, et al. Support for viral persistence in bats from age-specific serology and models of maternal immunity. Sci Rep. 2018;8:3859.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  91. Plowright RK, Field HE, Smith C, Divljan A, Palmer C, Tabor G, et al. Reproduction and nutritional stress are risk factors for Hendra virus infection in little red flying foxes (Pteropus scapulatus). Proc Biol Sci. 2008;275:861–9.

    Article  PubMed  PubMed Central  Google Scholar 

  92. Epstein JH, Baker ML, Zambrana-Torrelio C, Middleton D, Barr JA, Dubovi E, et al. Duration of maternal antibodies against canine distemper virus and Hendra virus in pteropid bats. PLoS One. 2013;8:e67584.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Peel AJ, Baker KS, Crameri G, Barr JA, Hayman DTS, Wright E, et al. Henipavirus neutralising antibodies in an isolated island population of African fruit bats. PLoS One. 2012;7:e30346.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Breed AC, Breed MF, Meers J, Field HE. Evidence of endemic Hendra virus infection in flying-foxes (Pteropus conspicillatus)—implications for disease risk management. PLoS One. 2011;6:e28816.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Jeong J, Smith CS, Peel AJ, Plowright RK, Kerlin DH, McBroom J, et al. Persistent infections support maintenance of a coronavirus in a population of Australian bats (Myotis macropus). Epidemiol Infect. 2017;145:2053–61.

    Article  CAS  PubMed  Google Scholar 

  96. Peel AJ, Wood JLN, Baker KS, Breed AC, Carvalho AD, Fernández-Loras A, et al. How does africa’s most hunted bat vary across the continent? Population traits of the straw-coloured fruit bat (Eidolon helvum) and its interactions with humans. Acta Chiropt. 2017;19:77–92.

    Article  Google Scholar 

  97. Dietrich M, Wilkinson DA, Benlali A, Lagadec E, Ramasindrazana B, Dellagi K, et al. Leptospira and paramyxovirus infection dynamics in a bat maternity enlightens pathogen maintenance in wildlife. Environ Microbiol. 2015;17(11):4280–9.

    Article  PubMed  Google Scholar 

  98. Tong S, Conrardy C, Ruone S, Kuzmin IV, Guo X, Tao Y, et al. Detection of novel SARS-like and other coronaviruses in bats from Kenya. Emerg Infect Dis. 2009;15:482–5.

    Article  PubMed  PubMed Central  Google Scholar 

  99. DeFrees SL, Wilson DE. Eidolon helvum. Mamm Species. 1988;312:1–5.

    Article  Google Scholar 

  100. Dengis CA. Taphozous mauritianus. Mamm Species. 1996;522:1–5.

    Article  Google Scholar 

  101. Dunlop J. Coleura afra. Mamm Species. 1997;566:1–4.

    Article  Google Scholar 

  102. Independent Scientific Group on Cattle TB. Bovine TB: the scientific evidence. London: Defra Publications; 2007.

    Google Scholar 

  103. Comte S, Umhang G, Raton V, Raoul F, Giraudoux P, Combes B, et al. Echinococcus multilocularis management by fox culling: an inappropriate paradigm. Prev Vet Med. 2017;147:178–85.

    Article  PubMed  Google Scholar 

  104. Morters MK, Restif O, Hampson K, Cleaveland S, Wood JLN, Conlan AJK. Evidence-based control of canine rabies: a critical review of population density reduction. J Anim Ecol. 2013;82:6–14.

    Article  PubMed  Google Scholar 

  105. Streicker DG, Recuenco S, Valderrama W, Gomez Benavides J, Vargas I, Pacheco V, et al. Ecological and anthropogenic drivers of rabies exposure in vampire bats: implications for transmission and control. Proc Biol Sci. 2012;279:3384–92.

    Article  PubMed  PubMed Central  Google Scholar 

  106. Blackwood JC, Streicker DG, Altizer S, Rohani P. Resolving the roles of immunity, pathogenesis, and immigration for rabies persistence in vampire bats. Proc Natl Acad Sci. 2013;110:20837–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  107. Fornes A, Lord RD, Kuns ML, Larghi OP, Fuenzalida E, Lazara L. Control of bovine rabies through vampire bat control. J Wildl Dis. 1974;10:310–6.

    Article  CAS  PubMed  Google Scholar 

  108. Amman BR, Nyakarahuka L, McElroy AK, Dodd KA, Sealy TK, Schuh AJ, et al. Marburgvirus resurgence in Kitaka mine bat population after extermination attempts, Uganda. Emerg Infect Dis. 2014;20:1761–4.

    Article  PubMed  PubMed Central  Google Scholar 

  109. Towner JS, Amman BR, Sealy TK, Carroll SAR, Comer JA, Kemp A, et al. Isolation of genetically diverse Marburg viruses from Egyptian fruit bats. PLoS Pathog. 2009;5:e1000536.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  110. Maganga GD, Bourgarel M, Ella GE, Drexler JF, Gonzalez J-P, Drosten C, et al. Is Marburg virus enzootic in Gabon? J Infect Dis. 2011;204(Suppl 3):S800–3.

    Article  CAS  PubMed  Google Scholar 

  111. López-Roig M, Serra-Cobo J. Impact of human disturbance, density, and environmental conditions on the survival probabilities of pipistrelle bat (Pipistrellus pipistrellus). Popul Ecol. 2014;56:471–80.

    Article  Google Scholar 

  112. Field HE. Bats and emerging zoonoses: henipaviruses and SARS. Zoonoses Public Health. 2009;56:278–84.

    Article  CAS  PubMed  Google Scholar 

  113. Kirkpatrick JF, Turner JW Jr. Compensatory reproduction in feral horses. J Wildl Manag. 1991;4:649–52.

    Article  Google Scholar 

  114. Langvatn R, Loison A. Consequences of harvesting on age structure, sex ratio and population dynamics of red deer Cervus elaphus in Central Norway. Wildlife Biol. 1999;5:213.

    Article  Google Scholar 

  115. Beasley JC, Olson ZH, Beatty WS, Dharmarajan G, Rhodes OE Jr. Effects of culling on mesopredator population dynamics. PLoS One. 2013;8:e58982.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  116. Bonesi L, Harrington LA, Maran T, Sidorovich VE, Macdonald DW. Demography of three populations of American mink mustela vison in Europe. Mamm Rev. 2006;36:98–106.

    Article  Google Scholar 

  117. Cowan PE. Effects of intensive trapping on breeding and age structure of brushtail possums, Trichosurus vulpecula, on Kapiti Island, New Zealand. N Z J Zool. 1993;20:1–11.

    Article  Google Scholar 

  118. Lindahl JF, Grace D. The consequences of human actions on risks for infectious diseases: a review. Infect Ecol Epidemiol. 2015;5:30048.

    PubMed  Google Scholar 

  119. Morse SS. Factors in the emergence of infectious diseases. In: Price-Smith AT, editor. Plagues and politics: infectious disease and international policy. London: Palgrave Macmillan UK; 2001. p. 8–26.

    Chapter  Google Scholar 

  120. Murray KA, Daszak P. Human ecology in pathogenic landscapes: two hypotheses on how land use change drives viral emergence. Curr Opin Virol. 2013;3:79–83.

    Article  PubMed  PubMed Central  Google Scholar 

  121. Ceballos G, Ehrlich PR, Barnosky AD, García A, Pringle RM, Palmer TM. Accelerated modern human–induced species losses: entering the sixth mass extinction. Sci Adv. 2015;1:e1400253.

    Article  PubMed  PubMed Central  Google Scholar 

  122. Voigt CC, Kingston T. Bats in the Anthropocene. In: Voigt CC, Kingston T, editors. Bats in the Anthropocene: conservation of bats in a changing world. Cham: Springer International Publishing; 2016. p. 1–9.

    Chapter  Google Scholar 

  123. Buttke DE, Decker DJ, Wild MA. The role of one health in wildlife conservation: a challenge and opportunity. J Wildl Dis. 2015;51:1–8.

    Article  PubMed  Google Scholar 

  124. Terraube J, Fernández-Llamazares Á, Cabeza M. The role of protected areas in supporting human health: a call to broaden the assessment of conservation outcomes. Curr Opin Environ Sustain. 2017;25:50–8.

    Article  Google Scholar 

Download references

Acknowledgements

We appreciate the support of the American people through the United States Agency for International Development (USAID) Emerging Pandemic Threats PREDICT project (cooperative agreement number GHN-A-OO-09-00010-00). We thank the governments of Rwanda, Tanzania, and Uganda for permission to conduct this study, and the country PREDICT teams that led and conducted the field activities. We are mostly grateful of all the technicians involved in the processing, testing, and barcoding of the samples, especially Ruth Maganga (Sokoine University of Agriculture) and Brett Smith (One Health Institute, University of California – Davis). Finally, we appreciate the help of Dr. Evan Eskew for supporting with STAN coding and Dr. Wes Johnson for guidance with Bayesian inference.

Funding

This study was supported by the generous support of the American people through the United States Agency for International Development (cooperative agreement numbers GHN-A-OO-09-00010-00 and AID-OAA-A− 14-00102). The contents of this paper are the responsibility of the authors and do not necessarily reflect the views of the US Agency for International Development or the US government.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

DML = Designed the project; collected, analyzed, and interpreted data; and drafted, edited, and prepared the final manuscript and figures. TG = Guided laboratory analysis and generated data, especially barcoding and sample testing, and supported manuscript writing. KG, DW, EVW, RK = Organized data collection and supported manuscript writing. BS, JN, ZS, MC = Data collection and supported manuscript writing. The PREDICT Consortium = Provided support in the design and implementation of surveillance and interpretation of data. JKM = Designed and supervised the project, including sampling design and data collection; supported data analysis and interpretation; and drafted, edited, and prepared the final manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jonna A. K. Mazet.

Ethics declarations

Ethics approval and consent to participate

Samples were collected from bats with permission from national authorities and under the Institutional Animal Care and Use Committee at the University of California, Davis (protocol number: 16048). Bats were released after sample collection.

Consent for publication

No human subjects were involved in this study.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1:

Summary of the inferred start date of the birth pulse, the end of the lactation period, and the start date of the mating period per microchiropteran species.

Additional file 2:

Bibliographic references for the traits of the bat species included in the study.

Additional file 3:

Imputation of the reproductive season to those bats whose biology is insufficiently known. Model methods, results, and literature cited. Table S3. Summary of the coefficients’ posterior probability distributions of the selected model for the imputation of the reproductive seasons of bats missing this data. Figure S3.1. The discrepancy between the observed (inferred) period and the period predicted by the selected imputation model considering 5 thousand posterior predictions of the observed period. Values of zero indicate matching, (observed period - predicted period). Figure S3.2. Left: the proportion of bats in each period after the imputation for those bats with the non-inferred season. Each line connects the proportion of individuals per reproductive season in each of the 5,000 Markov Chain Monte Carlo iteration. Right: the distribution of the crude coronavirus detection per reproductive season across the 5,000 Markov Chain Monte Carlo sampling iterations after imputing the periods when un-inferred. The black boxplots show the distribution of the coronavirus detection per period, while the light and dark colored boxes above and below show the interquartile detection in non-adults and adult bats, respectively, per period.

Additional file 4: Table S4.

Summary of the Alpha- and Betacoronaviruses (alphaCoV and betaCoV, respectively) found in the microbats tested.

Additional file 5:

Results of the model to assess the association between CoV shedding and the “Recent weaning” season in Eidolon helvum and microbats. Figure S5.1. Posterior predictive distributions of coronavirus positive bats (histograms) and the observed coronavirus positive bats (vertical lines). The histograms show the distribution of four thousand predictions of detection in the sampled bats. The lines show the observed detection. A) Coronavirus detection across all bats. B) Coronavirus detection across age categories: adults (light blue) and non-adults (yellow). C) Coronavirus detection across the reproductive seasons: “Not recent weaning” (light blue) and “Recent weaning” (yellow). D) Coronavirus detection across the reproductive seasons in the non-adult bats: “Not recent weaning” (light blue) and “Recent weaning” (yellow). E) Coronavirus detection across the life history seasons in the adult individuals: “Not recent weaning” (light blue) and “Recent weaning” (yellow). Figure S5.2. The density of the coefficients’ posterior probability distributions of the selected model. RW refers to the “Recent weaning” period (versus “Not recent weaning” season).

Additional file 6:

Sampling events -specific random intercepts values.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Montecino-Latorre, D., Goldstein, T., Gilardi, K. et al. Reproduction of East-African bats may guide risk mitigation for coronavirus spillover. One Health Outlook 2, 2 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s42522-019-0008-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s42522-019-0008-8

Keywords