Selective deforestation and exposure of African wildlife to bat-borne viruses

Similar content being viewed by others

Introduction

Spillover of viruses from wildlife to humans is often thought to be preceded by viral transmission and amplification among wildlife. For example, human ebolavirus outbreaks in Africa follow sylvatic transmission cycles in non-human primates and ungulates, with humans likely becoming infected through contact with carcasses1,2,3. Similarly, epidemiological data and analyses of inferred viral genomic recombination suggest that approximately half of human-infecting coronaviruses underwent transmission from wildlife reservoirs to humans through intermediary hosts4,5. Despite the high social and economic costs of zoonoses6, the mechanisms underlying such antecedent virus transmission within animals remain poorly understood.

Budongo Forest Reserve, western Uganda, contains approximately 482 km2 of medium-altitude, semi-deciduous forest7 and is located in the Albertine Rift, a region of exceptional biodiversity and endemism8. Until approximately 2008, the swamp forests of Budongo contained Raphia farinifera, a palm that, when decaying, provided a high-quality source of essential dietary minerals to wildlife9. Between 2006 and 2012, tobacco farming increased markedly in the area due to rising international demand and incentives from tobacco companies with longstanding operations in Uganda10. As a result, local farmers nearly extirpated R. farinifera because of its usefulness for making strings on which to dry tobacco leaves9,11. Budongo's eastern chimpanzees (Pan troglodytes schweinfurthii) altered their feeding behavior in response to this loss of a primary source of dietary minerals, increasingly consuming alternative sources such as clay, termite mounds, and the decaying pith of other tree species9. In 2017, we observed a never-before documented behavior by several species of wildlife in Budongo, including chimpanzees: the consumption of bat guano.

Here we present the results of an investigation as to whether this behavior could be an adaptation to dietary mineral scarcity, analogous to what has been documented for similar behaviors at this site9. We also investigate whether guano consumption could be an ecological mechanism whereby wildlife such as chimpanzees might be exposed to bat-borne viruses. We document high frequencies of guano consumption by three species of wildlife in Budongo, high concentrations of essential dietary minerals in the guano, and diverse bat-borne viruses in the guano, including a novel betacoronavirus within the Hibecovirus subgenus. These results illustrate how remote upstream forces can induce unanticipated causal chains that alter wildlife ecology and behavior, one result of which may be to increase virus spillover risk.

Results

Field studies

On June 25, 2017, we first observed chimpanzees consuming bat guano from under a large, hollow tree (Mildbraediodendron excelsum) in which a colony of Noack's roundleaf bat (Hipposideros ruber) was roosting (Fig. 1). Using a trail camera, we captured video images of chimpanzees, black-and-white colobus (Colobus guereza occidentalis) and red duiker (Cephalophus natalensis) repeatedly consuming guano from beneath the tree (Fig. 1). Animals consumed the guano directly, and not incidentally (e.g. from consumption of adjacent contaminated clay and water), as evidenced by the clearly visible selection of the guano itself during all instances and the presence of excavations and characteristic hand prints in the guano after animals had left (Supplementary Fig. 1a). We recorded 92 separate instances of guano consumption by chimpanzees on 71 different days, with between 1 and 13 chimpanzees per instance. Chimpanzees removed and ate guano with their hands (Fig. 1b, Supplementary Videos 1 and 2), or they drank adjacent water using a leaf sponge (folded leaves used to collect water12; Supplementary Video 3). Cameras captured black-and-white colobus feeding on guano on 65 occasions during 56 different days, with between 1 and 9 individuals per instance. These primates ate guano directly (Fig. 1c; Supplementary Video 4). Cameras captured solitary red duikers on 682 occasions on 210 different days. Duikers either licked guano directly or drank adjacent water next to the pile (Fig. 1d; Supplementary Video 5). On one occasion, we observed a ~ 2 m human-modified pole, suggesting that local people had also visited this tree, perhaps to collect guano (Supplementary Fig. 1b).

Images of guano pile (a, arrow) and chimpanzees (b) black-and-white colobus (c), and red duiker (d) consuming bat guano were captured using trail cameras between July 5th and October 18th 2017 and between September 14th 2018 and April 28th 2019.

Dietary mineral analyses

Nutritional analyses revealed that the guano contained concentrations of magnesium, phosphorus and potassium higher than in any other recorded dietary source of minerals at Budongo (Table 1). The guano also contained concentrations of sodium approximately equal to that of decaying Cleistopholis patens wood, the primary alternative source of dietary sodium for chimpanzees subsequent to the loss of R. farinifera9 (Table 1). Concentrations of calcium, manganese and iron were within ranges of other sources at Budongo (Table 1).

Virus identification and characterization

Metagenomic analyses of the bat guano revealed 27 novel putative eukaryotic viruses with 30.2-92.7% amino acid identity to known viruses of 12 families (Coronaviridae, Dicistroviridae, Hepeviridae, Iflaviridae, Nodaviridae, Parvoviridae, Picobirnaviridae, Picornaviridae, Permutotetraviridae, Polycipiviridae, Reoviridae and Totiviridae) and to 7 currently unclassified viruses (Supplementary Table 1). Individual guano samples analyzed contained an average of 14.5 viruses (standard deviation 3.6) that varied in prevalence from 9 to 100% among samples and in abundance over approximately 4 orders of magnitude, with arthropod-infecting viruses generally most prevalent and most abundant, consistent with the insectivorous diet of H. ruber (Supplementary Fig. 2). Sequences corresponding to a novel betacoronavirus (Coronaviridae: Betacoronavirus) were present in 6 samples (55%) (Supplementary Fig. 2). Due to the public health significance of betacoronaviruses, we intensively sequenced this virus, Buhirugu virus 1 (BHRGV-1), using the sample with the most abundant reads (sample 9 in Supplementary Fig. 2) and succeeded in obtaining 15,433 bases of the orf1a/b polyprotein gene and 3,181 bases of the spike protein gene (GenBank OP199247). Phylogenetic analyses (Fig. 2) show BHRGV-1 to be a novel member of the subgenus Hibecovirus, approximately equidistant from Bat Hp-betacoronavirus and Zaria bat coronavirus13,14. BHRGV-1 and the other hibecoviruses form a well-supported clade most closely related to viruses of the separate subgenus Sarbecovirus, which contains SARS-CoV and SARS-CoV-2 (Fig. 2)15.

Trees of BHRGV-1 (bold, red dot) and representative betacoronaviruses based on the viral RNA-dependent RNA polymerase (A) and spike protein (B) genes were constructed from codon-based nucleotide alignments (15,816 and 5334 positions, respectively). Taxon labels are abbreviations (See Table S1 for full descriptions) with viruses from humans (H) or from hosts of other mammalian orders (Ar = Artiodactyla; Ca = Carnivora; Ch = Chiroptera; Eu = Eulipotyphla; La = Lagomorpha; Pe = Perissodactyla; Ph = Pholidota; Ro = Rodentia) in parentheses. Circles on nodes indicate bootstrap values based on 1,000 replicates (black = 100%; gray = 75%-99%; white = 50%-74%); only values ≥ 50% are shown. Colored groupings correspond to Betacoronavirus subgenera (top to bottom): Sarbecovirus, Hibecovirus, Nobecovirus, Merbecovirus, and Embecovirus. Trees were midpoint rooted, yielding root placement the same as that of the International Committee on the Taxonomy of Viruses Betacoronavirus clade within the Coronaviridae15. Scale bar indicates nucleotide substitutions per site.

To investigate the potential host range of BHRGV-1, we conducted predicted protein structure and in silico docking analysis of the BHRGV-1 spike protein (Supplementary Fig. 3) and the angiotensin II (ACE2) receptors of humans and the other mammals observed consuming bat guano (Supplementary Fig. 4). In cases where the ACE2 nucleotide sequence of a particular animal was not available, we used sequences from a closely related species. The Ramachandran scores for BGHRV-1 S and the various species of ACE2 within the energetically favored region of the protein ranged from 95 to 98% (Supplementary Table 3). Docking analyses of the BHRGV-1 spike protein indicate non-permissive binding interactions between BHRGV-1 S and ACE2 receptors in all sequences analyzed (Supplementary Table 4), implying that ACE2 may not be the Hibecovirus receptor. To investigate other potential receptors, we repeated this analysis on Aminopeptidase-N, Dipeptidyl peptidase 4, and CEACAM116, with similar results indicating non-permissive binding interactions between BHRGV-1 S and each of these molecules (Supplementary Table 4, Supplementary Figs. 5-7 and Supplementary Tables 2-4). Finally, we examined whether BHRGV-1 has a predicted hemagglutinin-esterase region, which would implicate use of O-acetylated sialic acids as a receptor17, but we found no evidence of such a region by scanning for matches in the InterPro protein signature databases18,19.

Discussion

Minerals are essential for physiological functioning, growth, reproduction and immunity20. Minerals are also often limiting in the core diets of wild animals21. Some cave-dwelling invertebrates, fish and salamanders consume bat guano to obtain minerals in their nutrient-limited subterranean environments22. However, to our knowledge, guano ingestion by forest-dwelling mammals has not previously been reported. Bat guano also contains nutrients critical to plant growth, such as nitrogen, phosphate, and potassium, making guano an efficient and widely used fertilizer23. This may explain why people appear to have visited the same tree where we documented guano consumption by wildlife. We note that another betacoronavirus has been described in bat guano collected as fertilizer in Thailand24, and that harvesting bat guano for this purpose is a widespread but underappreciated practice that may increase pandemic risk25.

Our results suggest that guano consumption by Budongo wildlife may be a behavioral adaptation to mineral scarcity. This inference is supported by a decades-long body of evidence showing that wildlife in Budono have responded to the disappearance of R. farinifera by seeking alternative mineral sources9,26,27. The guano contained concentrations of potassium, magnesium, sodium and phosphorus equal to or in excess of concentrations in other dietary sources. Past studies have shown that consumption of alternative sources of minerals by Budongo chimpanzees began with the disappearance of R. farinifera9,26,27. Black-and-white colobus and duiker have not been as intensively studied in Budongo, so it is unknown whether guano consumption is also a new behavior for these animals. Black-and-white colobus frequently consume soil, clay, aquatic plants and even cement, demonstrating extreme dietary plasticity with respect to mineral acquisition28.

Guano consumption also appears to expose wildlife to bat-associated viruses, to the extent that the sequences we obtained represent infectious viruses. BHRGV-1 is a member of the subgenus Hibecovirus, which contains viruses that primarily infect bats of the genus Hipposideros but have also been documented in the bat genera Macronycteris, Nycteris, and Rhinolophus29, and Hibecovirus is sister taxon to the subgenus Sarbecovirus, which contains SARS-CoV and SARS-CoV-2. Predicted protein structure analysis of the BHRGV-1 S protein indicates the highest structural similarity to the SARS-CoV S protein (Supplementary Fig. 3 and Supplementary Table 2). However, the binding affinity of the BHRGV-1 S protein for all ACE2 proteins examined (including those of H. armiger, a close relative of H. ruber) is low (Supplementary Fig. 4 and Supplementary Tables 3 and 4), indicating that ACE2 is likely not a viable receptor for host cell entry in these mammals. We obtained similar results for the putative alternative receptors Aminopeptidase-N, Dipeptidyl peptidase 4, and CEACAM116 (Supplementary Figs. 5-7 and Supplementary Tables 3 and 4), and we found no evidence that BHRGV-1 has a predicted hemagglutinin-esterase region that might bind O-acetylated sialic acids, as has been shown for other coronaviruses17. If none of these molecules are, in fact, receptors for BHRGV-1 and other hibecoviruses, this would merit further study, especially for predicting the host range and zoonotic potential of coronaviruses in subgenera other than Sarbecovirus.

Coronavirus infections of wildlife have not been documented in Budongo to date. Examining local wildlife for evidence of BHRGV-1 or similar viruses in feces could yield additional information about the breadth of species that might have been exposed to bat-borne coronaviruses. However, multiple outbreaks of respiratory disease in the chimpanzees of Budongo have been observed, with the causes remaining undiagnosed. Respiratory disease outbreaks in other chimpanzee populations in Uganda have resulted from cross-species transmission of viruses from humans30,31, and human betacoronavirus OC43 can infect wild chimpanzees and cause clinical disease32.

Intriguingly, chimpanzees, black-and-white colobus, and red duiker have all been implicated in ebolavirus outbreaks in Central and West Africa1,33,34. The natural history of the ebolaviruses is poorly understood, but multi-host models of sylvatic ebolavirus transmission posit that outbreaks occur when primates and ungulates become infected by bats and serve as amplifying hosts2,3,35. Similarly, many bat-borne coronaviruses have emerged in humans after transmission through intermediary hosts4,5. Mechanisms of virus transmission from bats to other wildlife in nature remain poorly understood, although consumption of fruit contaminated by bats35,36 and contact with viruses shed into the environment37,38 have been hypothesized. Our data suggest another plausible ecological mechanism for exposure of wildlife to bat-associated viruses: consumption of bat guano as a source of dietary minerals. Tropical forest plants and soils are mineral-poor39. Depletion of primary sources of minerals such as R. farinifera could create conditions that favor guano consumption as a "fallback" mineral source.

Infectious disease emergence is often attributed to drivers such as land conversion, hunting, urbanization, climate change, and agricultural intensification40, but the ecological mechanisms whereby these drivers lead to cross-species pathogen transmission remain imprecisely understood. Our results provide an illustration of how these mechanisms might follow elaborate causal chains. In Budongo, international demand for tobacco caused local selective deforestation and loss of a primary source of dietary minerals, which led to fallback consumption of guano by wildlife and exposure of wildlife to bat-associated viruses, including a congeneric of the pandemic SARS coronaviruses. Mathematical tools for representing causal chains and causal networks are becoming widespread in epidemiology41 and might prove useful for assessing how environmental and social drivers ultimately lead to zoonotic transmission. This understanding, in turn, could lead to improved precision in the application of tools for preventing pandemics. For example, compared to the costs of a pandemic, the costs of offering local farmers substitutes to R. farinifera for making strings to dry tobacco leaves would likely be trivial6,11. In general, we suggest that understanding causal chains and identifying their "breakable" links holds promise for illuminating disease ecology and improving zoonoses prevention.

Methods

Field studies

The study took place in the Budongo Forest Reserve, Uganda. We first observed chimpanzees in the habituated Waibira community42,43 feeding on bat guano in a hollow tree on the 25th of June, 2017, even though chimpanzees had regularly been observed since 2011. On the 5th of July 2017, we installed a trail camera (Bushnell Trophy Cam, model 119,774) with the following default settings: 10 s video length and interval, auto sensor level, low night vision shutter, and 24 h camera mode. We mounted the camera to a tree 6 m from the tree atop the guano at a height of 1 m. Image capture occurred between the 6th of July and the 18th of October 2017 and again between the 14th of September 2018 and the 28th of April 2019. We analyzed the resulting 14,567 10 s video recordings (40.46 camera-hours in total) for the presence and number of animals feeding on the guano. We defined an "instance" as a set of sequential video recordings of animals of a given species feeding on the guano, separated from the previous instance by at least 30 min during which no individuals of the same species were recorded.

We collected guano samples from this hollow tree monthly from 13th September 2018 to 29th April 2019 (Supplementary Fig. 2) and divided them for mineral content analysis and virus identification. We oven-dried samples for mineral content analysis (approximately 50 g) and stored and shipped them to the USA at ambient temperature. We placed samples for molecular analysis (approximately 0.9 ml) in 1.8 ml sterile cryovials, mixed them thoroughly with an equal volume of RNAlater Stabilization Solution (Thermo Fisher, Waltham, MA, USA), and stored them cold in the field and in liquid nitrogen prior to and during shipment to the USA.

All animal use was strictly non-invasive and observational. The study protocol was reviewed and approved by the Uganda Wildlife Authority and the Uganda National Council for Science and Technology, and was in compliance with the guidelines of the Animal Welfare and Ethical Review Body of the University of Stirling and all applicable regulations governing the protection of animals and research. We have complied with all relevant ethical regulations for animal use. The species of animals were Pan troglodytes schweinfurthii, Hipposideros ruber, Colobus guereza occidentalis, and Cephalophus natalensis, and all were wild-type and of undetermined sex and age.

Dietary mineral analyses

Prior to analysis, we inactivated samples with ultraviolet radiation and oven-dried them at 105 °C for 48 h. We then digested samples using the MARS 6 Microwave Digestion System (CEM Corporation, Matthews, NC, USA) and analyzed them for Ca, P, Mg, K, Na, Fe, Zn, Cu, Mn, and Mo on an iCAP 6300 inductively coupled plasma radial spectrometer (Thermo Fisher, Waltham, MA, USA)44.

Metagenomic analyses

We processed guano samples for virus identification using metagenomic methods31,45. Briefly, we added 200 μl of guano+RNAlater to 800 μl of Hanks' balanced salt solution and homogenized them in PowerBead Tubes (Qiagen, Hilden, Germany) containing 2.38 mm metal beads. We then treated the homogenate with nucleases to reduce unencapsidated nucleic acids46. We used the QIAmp MinElute Virus Spin Kit (Qiagen, Hilden, Germany) to isolate total nucleic acids, and we and converted RNA to double-stranded cDNA using the SuperScript double-stranded cDNA Synthesis Kit (Invitrogen, Carlsbad, CA, USA). We cleaned cDNA using Agencourt AmpureXP beads (Beckman Coulter, Brea, CA, USA) and synthesized DNA libraries using the Nextera XT DNA sample preparation kit (Illumina, San Diego, CA, USA). We sequenced libraries on an Illumina MiSeq instrument using 600 cycle v3 MiSeq Reagent Kits. We trimmed resulting sequences at a Phred quality score <30, discarded reads <50 bp, and removed sequences matching host genomes and known contaminants. We thereby sequenced samples to a mean depth after quality and length trimming and host genome subtraction of 1.9 M reads (standard error 0.2 M reads), ranging from 1.4 to 3.5 M reads per sample.

Data processing, bioinformatics, statistics, and reproducibility

We subjected sequence reads to de novo assembly using SPAdes 3.13.047, discarded resulting contiguous sequences (contigs) <500 nucleotides, and used cd-hit48 to remove redundant contigs (90% similarity threshold). We compared remaining contigs to custom databases of representative virus protein sequences and to the NCBI non-redundant protein sequence database using blastx49. We ran and analyzed blank samples in parallel to ensure that cross contamination had not occurred.

To investigate BHRGV-1 in greater detail, we queried the initial 15,433 bp contig representing this virus against the National Center for Biotechnology Information (NCBI) nucleotide database using blastn in BLAST + 50,51. Based on this analysis, we chose two reference sequences for downstream comparisons: bat Hp-betacoronavirus/Zhejiang2013 (NCBI accession ID NC_025217.1) with 74.5% nucleotide sequence identity to BHRGV-1 and Zaria bat coronavirus strain ZBCoV (Genbank: HQ166910.1) with 76.4% nucleotide sequence identity to BHRGV-1. We then mapped sequence reads from the sample containing the highest concentration of this virus (sample 9, Fig. S2) against each reference separately using CLC Genomics Workbench (Qiagen, Hilden, Germany), specifying a minimum length fraction of 0.5 and a minimum similarity fraction of 0.8. We extended mapped regions in CLC Genomics Workbench using iterative mapping of previously unmapped reads to consensus sequences extracted from each prior iteration. We collected reads thus identified, assembled them de novo, and aligned the resulting contigs to both references to create a genome scaffold (each reference was useful for assembling different regions of the BHRGV-1 genome). We then mapped contigs and reads that did not form contigs against the genome scaffold to create a final draft 32,594 bp sequence of the BHRGV-1 ORF1ab polyprotein and spike protein (S) genes. We constructed phylogenetic trees of the BHRGV-1 ORF1ab polyprotein and spike protein (S) genes using PhyML52 with smart model selection53 (GTR + I models selected in both cases) and 1000 bootstrap replicates to assess statistical confidence in clades.

We modeled 3D-structures of the receptor binding domain (RBD) of BHRGV-1 S and the putative receptor proteins ACE2, Aminopeptidase-N, Dipeptidyl peptidase 4, and CEACAM116 from selected animal species using Modeler 10.254. In cases for which a particular species lacked an available representative sequence, we chose the closest phylogenetic relative of that species for which a representative sequence was available: Hipposideros armiger (great roundleaf bat) in place of Hipposideros ruber (Noack's roundleaf bat); Capra hircus (goat) in place of Cephalophus natalensis (red duiker); and Colobus angolensis palliates (Angolan black-and-white colobus) in place of Colobus guereza occidentalis (black-and-white colobus) (Table S2). We modeled the BHRGV-1 S protein using the well characterized SARS-CoV spike fusion protein (PDB ID: 2BEZ) as a homologous structural template. Similarly, we used the human ACE2 (PDB ID: 1R42), Aminopeptidase-N (PDB ID: 5LHD), Dipeptidyl peptidase 4 (PDB ID: 2QT9), and CEACAM1 (PDB ID: 4QXW) proteins as homologous structural templates in the other species analyzed (Table S2). We evaluated the quality of the resulting models using the GA341 score55, DOPE (Discrete Optimized Protein Energy) method scores56, and the SWISS-MODEL structure assessment server57. We then refined structures with the lowest DOPE scores using molecular dynamics (MD) simulations and further analyzed them for quality using Ramachandran Plot and MolProbity in SWISS-MODEL (49) (Table S3).

We used HDOCK58 to model putative receptor/BGHRV-1 binding complexes. HDOCK maps receptor and ligand protein molecules onto grids, then "docks" two molecules using a hierarchical approach based on fast Fourier transformation. To minimize bias, we applied the template-free docking method with structures generated by Modeler 10.2. We optimized the final docked protein complexes using the AMBER99SB-ILDN force field in GROMACS59. Specifically, docked complexes were immersed in a truncated octahedron box of TIP3P water molecules. The solvated box was further neutralized with Na + or Cl − counter ions using the tleap program. We used Particle Mesh Ewald (PME) to calculate long-range electrostatic interactions, with a cut-off distance for long-range van der Waals (VDW) energy term of 12.0 Å, and the system minimized without restraints. We applied 2500 cycles of steepest descent minimization followed by 5000 cycles of conjugate gradient minimization. We initiated MD simulations by heating each system in the NVT ensemble from 0 to 300 K for 50 ps using a Langevin thermostat with a coupling coefficient of 1.0/ps and a force constant of 2.0 kcal/mol·Å2 on the complex. We ran the MD simulation for 100 ns at a constant temperature of 300 K in the NPT ensemble with periodic boundary conditions for each system. During the MD procedure, we applied the SHAKE algorithm was to all covalent bonds involving hydrogen atoms, with a time step of 2 fs. We calculated free energies of binding for all simulated docked structures using the molecular mechanics Poisson Boltzmann surface area (MM-PBSA) tool in GROMACS 202260 (Table S4). To examine BGHRV-1 for a predicted hemagglutinin-esterase region, we searched for matches in the InterPro protein signature databases using InterProScan 5.65-97.018,19.

Data availability

All raw sequence reads were deposited in the NIH National Center for Biotechnology Information (NCBI) Sequence Read Achieve under BioProject PRJNA1087330 (accession numbers SAMN40440184-SAMN40440195). All assembled virus genome sequences were deposited in NCBI GenBank under accession numbers OP199247 and OP834146-OP834171.

References

1. Leroy, E. M. et al. Multiple Ebola virus transmission events and rapid decline of central African wildlife. Science 303, 387-390 (2004).

   Article  CAS  PubMed  Google Scholar 

2. Caron, A. et al. Ebola virus maintenance: if not (only) bats, what else? Viruses 10, https://doi.org/10.3390/v10100549 (2018).

3. Gonzalez, J. P., Souris, M., Sylla, M., Veas, F. & Vincent, T. In Emerging Challenges in Filovirus Infections (ed. Okware, S. I.) 9-28 (IntechOpen, 2019).

4. Latif, A. A. & Mukaratirwa, S. Zoonotic origins and animal hosts of coronaviruses causing human disease pandemics: A review. Onderstepoort J. Vet. Res 87, e1-e9 (2020).

   Article  PubMed  Google Scholar 

5. Ye, Z. W. et al. Zoonotic origins of human coronaviruses. Int J. Biol. Sci. 16, 1686-1697 (2020).

   Article  CAS  PubMed  PubMed Central  Google Scholar 

6. Bernstein, A. S. et al. The costs and benefits of primary prevention of zoonotic pandemics. Sci. Adv. 8, eabl4183 (2022).

   Article  CAS  PubMed  PubMed Central  Google Scholar 

7. Eggling, W. J. Observations on the ecology of the Budongo Rain Forest, Uganda. J. Ecol. 34, 20-87 (1969).

   Article  Google Scholar 

8. Plumptre, A. J. et al. The biodiversity of the Albertine Rift. Biol. Conserv 134, 178-194 (2007).

   Article  Google Scholar 

9. Reynolds, V., Lloyd, A. W. & English, C. J. Adaptation by Budongo Forest chimpanzees (Pan troglodytes schweinfurthii) to loss of a primary source of dietary sodium. Afr. Primates 7, 156-162 (2012).

   Google Scholar 

10. Sejjaaka, S. in International Businesses and the Challenges of Poverty in the Developing World: Case Studies on Global Responsibilities and Practices (eds Bird, F. & Herman, S. W.) 111-123 (Palgrave Macmillan UK, 2004).

11. Kinights, J. Raphia Conservation Project. (Budongo Conservation Field Station, 2009).

12. Sousa, C., Biro, D. & Matsuzawa, T. Leaf-tool use for drinking water by wild chimpanzees (Pan troglodytes): acquisition patterns and handedness. Anim. Cogn. 12, S115-S125 (2009).

    Article  PubMed  Google Scholar 

13. Quan, P. L. et al. Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria. mBio 1, https://doi.org/10.1128/mBio.00208-10 (2010).

14. Wu, Z. et al. Deciphering the bat virome catalog to better understand the ecological diversity of bat viruses and the bat origin of emerging infectious diseases. ISME J. 10, 609-620 (2016).

    Article  PubMed  Google Scholar 

15. Woo, P. C. Y. et al. ICTV virus taxonomy profile: Coronaviridae 2023. J. Gen. Virol. 104, https://doi.org/10.1099/jgv.0.001843 (2023).

16. Millet, J. K., Jaimes, J. A. & Whittaker, G. R. Molecular diversity of coronavirus host cell entry receptors. FEMS Microbiol. Rev. 45, https://doi.org/10.1093/femsre/fuaa057 (2020).

17. Robson, B. Bioinformatics studies on a function of the SARS-CoV-2 spike glycoprotein as the binding of host sialic acid glycans. Comput. Biol. Med 122, 103849 (2020).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

18. Jones, P. et al. InterProScan 5: genome-scale protein function classification. Bioinforma. (Oxf., Engl.) 30, 1236-1240 (2014).

    CAS  Google Scholar 

19. Blum, M. et al. The InterPro protein families and domains database: 20 years on. Nucleic Acids Res 49, D344-D354 (2021).

    Article  CAS  PubMed  Google Scholar 

20. McDowell, L. R. Minerals in animal and human nutrition. (Elsevier, 2003).

21. Kaspari, M. The seventh macronutrient: how sodium shortfall ramifies through populations, food webs and ecosystems. Ecol. Lett. 23, 1153-1168 (2020).

    Article  PubMed  Google Scholar 

22. Fenolio, D. B., Graening, G. O., Collier, B. A. & Stout, J. F. Coprophagy in a cave-adapted salamander; the importance of bat guano examined through nutritional and stable isotope analyses. Proc. Biol. Sci. 273, 439-443 (2006).

    CAS  PubMed  Google Scholar 

23. Sakoui, S. et al. The life hidden inside caves: ecological and economic importance of bat guano. Int. J. Ecol. 2020, 9872532 (2020).

    Article  Google Scholar 

24. Wacharapluesadee, S. et al. Group C betacoronavirus in bat guano fertilizer, Thailand. Emerg. Infect. Dis. 19, 1349-1351 (2013).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

25. Osofsky, S. A., Lieberman, S., Walzer, C., Lee, H. L. & Neme, L. A. An immediate way to lower pandemic risk: (not) seizing the low-hanging fruit (bat). Lancet Planet Health 7, e518-e526 (2023).

    Article  PubMed  Google Scholar 

26. Reynolds, V., Pascual-Garrido, A., Lloyd, A. W., Lyons, P. & Hobaiter, C. Possible mineral contributions to the diet and health of wild chimpanzees in three East African forests. Am. J. Primatol. 81, e22978 (2019).

    Article  PubMed  Google Scholar 

27. Reynolds, V. et al. Mineral acquisition from clay by Budongo Forest chimpanzees. PLoS One 10, e0134075 (2015).

    Article  PubMed  PubMed Central  Google Scholar 

28. Pebsworth, P. A., Huffman, M. A., Lambert, J. E. & Young, S. L. Geophagy among nonhuman primates: A systematic review of current knowledge and suggestions for future directions. Am. J. Phys. Anthropol. 168, 164-194 (2019).

    Article  PubMed  Google Scholar 

29. Chidoti, V. et al. Longitudinal survey of coronavirus circulation and diversity in insectivorous bat colonies in Zimbabwe. Viruses 14, https://doi.org/10.3390/v14040781 (2022).

30. Scully, E. J. et al. Lethal respiratory disease associated with human rhinovirus C in wild chimpanzees, Uganda, 2013. Emerg. Infect. Dis. 24, 267-274 (2018).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

31. Negrey, J. et al. Simultaneous outbreaks of respiratory disease in wild chimpanzees caused by distinct viruses of human origin. Emerg. Microbes Infect. 8, 139-149 (2019).

    Article  PubMed  PubMed Central  Google Scholar 

32. Patrono, L. V. et al. Human coronavirus OC43 outbreak in wild chimpanzees, Cote d Ivoire, 2016. Emerg. microbes Infect. 7, 118 (2018).

    Article  PubMed  PubMed Central  Google Scholar 

33. Formenty, P. et al. Ebola virus outbreak among wild chimpanzees living in a rain forest of Cote d'Ivoire. J. Infect. Dis. 179, S120-S126 (1999).

    Article  PubMed  Google Scholar 

34. Lahm, S. A., Kombila, M., Swanepoel, R. & Barnes, R. F. Morbidity and mortality of wild animals in relation to outbreaks of Ebola haemorrhagic fever in Gabon, 1994-2003. Trans. R. Soc. Trop. Med Hyg. 101, 64-78 (2007).

    Article  PubMed  Google Scholar 

35. Olival, K. J. & Hayman, D. T. Filoviruses in bats: current knowledge and future directions. Viruses 6, 1759-1788 (2014).

    Article  PubMed  PubMed Central  Google Scholar 

36. Becquart, P. et al. High prevalence of both humoral and cellular immunity to Zaire ebolavirus among rural populations in Gabon. PLoS One 5, e9126 (2010).

    Article  PubMed  PubMed Central  Google Scholar 

37. Forbes, K. M., Sironen, T. & Plyusnin, A. Hantavirus maintenance and transmission in reservoir host populations. Curr. Opin. Virol. 28, 1-6 (2018).

    Article  PubMed  Google Scholar 

38. Becker, D. J., Eby, P., Madden, W., Peel, A. J. & Plowright, R. K. Ecological conditions predict the intensity of Hendra virus excretion over space and time from bat reservoir hosts. Ecol. Lett. 26, 23-36 (2023).

    Article  PubMed  Google Scholar 

39. Camenzind, T., Hättenschwiler, S., Treseder, K. K., Lehmann, A. & Rillig, M. C. Nutrient limitation of soil microbial processes in tropical forests. Ecol. Monogr. 88, 4-21 (2018).

    Article  Google Scholar 

40. Plowright, R. K. et al. Land use-induced spillover: a call to action to safeguard environmental, animal, and human health. Lancet Planet Health 5, e237-e245 (2021).

    Article  PubMed  PubMed Central  Google Scholar 

41. Ackley, S. F., Lessler, J. & Glymour, M. M. Dynamical modeling as a tool for inferring causation. Am. J. Epidemiol. 191, 1-6 (2022).

    Article  PubMed  Google Scholar 

42. Reynolds, V. The Chimpanzees of the Budongo Forest: Ecology, Behaviour and Conservation. (Oxford University Press, 2005).

43. Samuni, L., Mundry, R., Terkel, J., Zuberbühler, K. & Hobaiter, C. Socially learned habituation to human observers in wild chimpanzees. Anim. Cognition 17, 997-1005 (2014).

    Article  Google Scholar 

44. Irwin, M. T., Raharison, J. L., Chapman, C. A., Junge, R. E. & Rothman, J. M. Minerals in the foods and diet of diademed sifakas: are they nutritional challenges? Am. J. Primatol. 79, 1-14 (2017).

    Article  CAS  PubMed  Google Scholar 

45. Negrey, J. D. et al. Demography, life-history trade-offs, and the gastrointestinal virome of wild chimpanzees. Philos. Trans. R. Soc. Lond. B Biol. Sci. 375, 20190613 (2020).

    Article  PubMed  PubMed Central  Google Scholar 

46. Allander, T., Emerson, S. U., Engle, R. E., Purcell, R. H. & Bukh, J. A virus discovery method incorporating DNase treatment and its application to the identification of two bovine parvovirus species. Proc. Natl. Acad. Sci. USA 98, 11609-11614 (2001).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

47. Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput Biol. 19, 455-477 (2012).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

48. Li, W. & Godzik, A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22, 1658-1659 (2006).

    Article  CAS  PubMed  Google Scholar 

49. States, D. J. & Gish, W. Combined use of sequence similarity and codon bias for coding region identification. J. Comput Biol. 1, 39-50 (1994).

    Article  CAS  PubMed  Google Scholar 

50. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403-410 (1990).

    Article  CAS  PubMed  Google Scholar 

51. Camacho, C. et al. BLAST+: architecture and applications. BMC Bioinforma. 10, 421 (2009).

    Article  Google Scholar 

52. Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 59, 307-321 (2010).

    Article  CAS  PubMed  Google Scholar 

53. Lefort, V., Longueville, J. E. & Gascuel, O. SMS: Smart model selection in PhyML. Mol. Biol. Evol. 34, 2422-2424 (2017).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

54. Sali, A. & Blundell, T. L. Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 234, 779-815 (1993).

    Article  CAS  PubMed  Google Scholar 

55. John, B. & Sali, A. Comparative protein structure modeling by iterative alignment, model building and model assessment. Nucleic Acids Res 31, 3982-3992 (2003).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

56. Shen, M. Y. & Sali, A. Statistical potential for assessment and prediction of protein structures. Protein Sci. 15, 2507-2524 (2006).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

57. Benkert, P., Biasini, M. & Schwede, T. Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics 27, 343-350 (2011).

    Article  CAS  PubMed  Google Scholar 

58. Yan, Y., Tao, H., He, J. & Huang, S. Y. The HDOCK server for integrated protein-protein docking. Nat. Protoc. 15, 1829-1852 (2020).

    Article  CAS  PubMed  Google Scholar 

59. Lundborg, M. & Lindahl, E. Automatic GROMACS topology generation and comparisons of force fields for solvation free energy calculations. J. Phys. Chem. B 119, 810-823 (2015).

    Article  CAS  PubMed  Google Scholar 

60. Kumari, R., Kumar, R., Open Source Drug Discovery, C. & Lynn, A. g_mmpbsa-a GROMACS tool for high-throughput MM-PBSA calculations. J. Chem. Inf. Model 54, 1951-1962 (2014).

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We are grateful to the Uganda Wildlife Authority and the Uganda National Council for Science and Technology for kindly granting permission to conduct field studies. We also thank the staff and management of Budongo Conservation Field Station for logistic support, and in particular Jacob Ariyo, Vicent Kiiza, Stephen Mugisha and Charles Rabu for assisting with field data collection. This work was supported in part through the European Research Council, grant agreement number 679787 ((to C. C.) and 802179) (to C.H.), the Royal Zoological Society of Scotland (to V.R.), the Armed Forces Health Surveillance Division (AFHSD), Global Emerging Infections Surveillance (GEIS) Branch, ProMIS ID P0167_22_NM (to K.A.B-L.), Navy WUN A1417 (to K.A.B-L.), and the University of Wisconsin-Madison John D. MacArthur Research Chair (to T.L.G.). The views expressed in this article are those of the authors and do not necessarily reflect the official policy or position of the Department of Defense, Department of the Navy, nor the U.S. Government. Several of the authors are U.S. Government employees. This work was prepared as part of their official duties. Title 17 U.S.C. § 105 provides that 'Copyright protection under this title is not available for any work of the United States Government.' Title 17 U.S.C. §101 defines a U.S. Government work as a work prepared by a military service member or employee of the U.S. Government as part of that person's official duties.

Author information

Author notes

1. These authors contributed equally: Pawel Fedurek, Caroline Asiimwe.

Authors and Affiliations

1. Division of Psychology, Faculty of Natural Sciences, University of Stirling, Stirling, FK9 4LA, UK

   Pawel Fedurek

2. Budongo Conservation Field Station, PO Box 362, Masindi, Uganda

   Pawel Fedurek, Caroline Asiimwe, Walter J. Akankwasa, Vernon Reynolds, Catherine Hobaiter, Geoffrey Muhanguzi & Klaus Zuberbühler

3. Biological Defense Research Directorate, Naval Medical Research Command, Fort Detrick, MD, 21702, USA

   Gregory K. Rice, Regina Z. Cer, Andrew J. Bennett & Kimberly A. Bishop-Lilly

4. Leidos, 1750 Presidents St, Reston, VA, 20190, USA

   Gregory K. Rice & Andrew J. Bennett

5. School of Anthropology, University of Oxford, 51/53 Banbury Road, Oxford, OX2 6PE, UK

   Vernon Reynolds

6. School of Psychology and Neuroscience, University of St Andrews; St Mary's Quad, South Street, St Andrews, KY16 9JP, UK

   Catherine Hobaiter & Klaus Zuberbühler

7. Department of Zoology, Entomology & Fisheries Sciences, Makerere University, PO Box 7062, Kampala, Uganda

   Robert Kityo

8. Institute of Biology, University of Neuchâtel, Rue Emile-Argand 11, CH-2000, Neuchâtel, Switzerland

   Klaus Zuberbühler

9. Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, 04103, Leipzig, Germany

   Catherine Crockford

10. Institut des Sciences Cognitives, 67 Bd Pinel, 69500, Bron, France

    Catherine Crockford

11. Department of Anthropology, Hunter College of the City University of New York, 695 Park Avenue, New York, NY, 10065, USA

    Jessica M. Rothman

12. School of Veterinary Medicine, Department of Pathobiological Sciences, University of Wisconsin-Madison, 1656 Linden Drive, Madison, WI, USA

    Tony L. Goldberg

Contributions

P.F., C.A., and T.L.G. contributed to the study conception and design. P.F., C.A., W.J.A., V.R., C.H., R. K., G.M., K.Z. and C.C. contributed to sample and data collection in the field. P.F., J.M.R., G.K.R., R.Z.C. A.J.B., K.B.L. and T.L.G. conducted laboratory analyses. G.K.R., R.Z.C., A.J.B., K.B.L. and T.L.G. performed statistical analyses and interpretations. P.F. and T.L.G. wrote the initial manuscript, and all authors read, edited and approved the final manuscript.

Corresponding author

Correspondence to Tony L. Goldberg.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval

The study protocol was reviewed and approved by the Uganda Wildlife Authority and the Uganda National Council for Science and Technology and was in compliance with the guidelines of the Animal Welfare and Ethical Review Body of the University of Stirling and all applicable regulations governing the protection of animals and research. We have complied with all relevant ethical regulations for animal use.

Peer review

Peer review information

Communications Biology thanks Daniel Becker, Katharina C. Wollenberg Valero and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: Tobias Goris. A peer review file is available.

Additional information

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

Supplementary information

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Fedurek, P., Asiimwe, C., Rice, G.K. et al. Selective deforestation and exposure of African wildlife to bat-borne viruses. Commun Biol 7, 470 (2024). https://doi.org/10.1038/s42003-024-06139-z

Download citation

• Received: 23 May 2023

• Accepted: 02 April 2024

• Published: 22 April 2024

• DOI: https://doi.org/10.1038/s42003-024-06139-z

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.