Abstract
The majority of human emerging infectious diseases are zoonotic, with viruses that originate in wild mammals of particular concern (for example, HIV, Ebola and SARS)1,2,3. Understanding patterns of viral diversity in wildlife and determinants of successful cross-species transmission, or spillover, are therefore key goals for pandemic surveillance programs4. However, few analytical tools exist to identify which host species are likely to harbour the next human virus, or which viruses can cross species boundaries5,6,7. Here we conduct a comprehensive analysis of mammalian host–virus relationships and show that both the total number of viruses that infect a given species and the proportion likely to be zoonotic are predictable. After controlling for research effort, the proportion of zoonotic viruses per species is predicted by phylogenetic relatedness to humans, host taxonomy and human population within a species range—which may reflect human–wildlife contact. We demonstrate that bats harbour a significantly higher proportion of zoonotic viruses than all other mammalian orders. We also identify the taxa and geographic regions with the largest estimated number of ‘missing viruses’ and ‘missing zoonoses’ and therefore of highest value for future surveillance. We then show that phylogenetic host breadth and other viral traits are significant predictors of zoonotic potential, providing a novel framework to assess if a newly discovered mammalian virus could infect people.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
$199.00 per year
only $3.90 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout




Similar content being viewed by others
References
Woolhouse, M. E. J. & Gowtage-Sequeria, S. Host range and emerging and reemerging pathogens. Emerg. Infect. Dis. 11, 1842–1847 (2005)
Jones, K. E. et al. Global trends in emerging infectious diseases. Nature 451, 990–993 (2008)
Taylor, L. H., Latham, S. M. & Woolhouse, M. E. J. Risk factors for human disease emergence. Phil. Trans. R. Soc. Lond. B 356, 983–989 (2001)
Morse, S. S. et al. Prediction and prevention of the next pandemic zoonosis. Lancet 380, 1956–1965 (2012)
Parrish, C. R. et al. Cross-species virus transmission and the emergence of new epidemic diseases. Microbiol. Mol. Biol. Rev. 72, 457–470 (2008)
Lipsitch, M. et al. Viral factors in influenza pandemic risk assessment. eLife 5, e18491 (2016)
Pulliam, J. R. C. & Dushoff, J. Ability to replicate in the cytoplasm predicts zoonotic transmission of livestock viruses. J. Infect. Dis. 199, 565–568 (2009)
Woolhouse, M. E., Haydon, D. T. & Antia, R. Emerging pathogens: the epidemiology and evolution of species jumps. Trends Ecol. Evol. 20, 238–244 (2005)
Cooper, N. et al. Phylogenetic host specificity and understanding parasite sharing in primates. Ecol. Lett. 15, 1370–1377 (2012)
Davies, T. J. & Pedersen, A. B. Phylogeny and geography predict pathogen community similarity in wild primates and humans. Proc. R. Soc. Lond. B 275, 1695–1701 (2008)
Gómez, J. M., Nunn, C. L. & Verdú, M. Centrality in primate-parasite networks reveals the potential for the transmission of emerging infectious diseases to humans. Proc. Natl Acad. Sci. USA 110, 7738–7741 (2013)
Luis, A. D. et al. A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special? Proc. R. Soc. Lond. B. 280, 20122753 (2013)
Brierley, L., Vonhof, M. J., Olival, K. J., Daszak, P. & Jones, K. E. Quantifying global drivers of zoonotic bat viruses: a process-based perspective. Am. Nat. 187, E53–E64 (2016)
Streicker, D. G. et al. Host phylogeny constrains cross-species emergence and establishment of rabies virus in bats. Science 329, 676–679 (2010)
Han, B. A., Schmidt, J. P., Bowden, S. E. & Drake, J. M. Rodent reservoirs of future zoonotic diseases. Proc. Natl Acad. Sci. USA 112, 7039–7044 (2015)
Fauquet, C ., Mayo, M. A ., Maniloff, J ., Desselberger, U. & Ball, L. A. Virus taxonomy: Eighth Report of the International Committee on Taxonomy of Viruses. (Elsevier Academic Press, 2005)
Dunn, R. R ., Davies, T. J ., Harris, N. C. & Gavin, M. C. Global drivers of human pathogen richness and prevalence. Proc. R. Soc. Lond. B 277, 2587–2595 (2010)
Levinson, J. et al. Targeting surveillance for zoonotic virus discovery. Emerg. Infect. Dis. 19, 743–747 (2013)
Morse, S. S. in Emerging Viruses (ed. Morse, S. S. ) 10–28 (Oxford University Press, 1993)
Zhou, P. et al. Contraction of the type I IFN locus and unusual constitutive expression of IFN-α in bats. Proc. Natl Acad. Sci. USA 113, 2696–2701 (2016)
Parker, I. M. et al. Phylogenetic structure and host abundance drive disease pressure in communities. Nature 520, 542–544 (2015)
Longdon, B., Brockhurst, M. A., Russell, C. A., Welch, J. J. & Jiggins, F. M. The evolution and genetics of virus host shifts. PLoS Pathog. 10, e1004395 (2014)
Ge, X.-Y. et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 503, 535–538 (2013)
Organtini, L. J., Allison, A. B., Lukk, T., Parrish, C. R. & Hafenstein, S. Global displacement of canine parvovirus by a host-adapted variant: A structural comparison between pandemic viruses with distinct host ranges. J. Virol. 89, 1909–1912 (2015)
Anthony, S. J. et al. A strategy to estimate unknown viral diversity in mammals. MBio 4, e00598–13 (2013)
Drexler, J. F. et al. Bats host major mammalian paramyxoviruses. Nat. Commun. 3, 796 (2012)
Geoghegan, J. L., Senior, A. M., Di Giallonardo, F. & Holmes, E. C. Virological factors that increase the transmissibility of emerging human viruses. Proc. Natl Acad. Sci. USA 113, 4170–4175 (2016)
Fritz, S. A., Bininda-Emonds, O. R. P. & Purvis, A. Geographical variation in predictors of mammalian extinction risk: big is bad, but only in the tropics. Ecol. Lett. 12, 538–549 (2009)
Nunn, C. L. & Altizer, S. M. The global mammal parasite database: An online resource for infectious disease records in wild primates. Evol. Anthropol. 14, 1–2 (2005)
Olival, K. J., Epstein, J. H., Wang, L. F., Field, H. E. & Daszak, P. in New Directions in Conservation Medicine: Applied Cases of Ecological Health (eds Aguirre, A. A., Ostfeld, R. S. & Daszak, P. ) Ch. 14, 195–212 (Oxford University Press, 2012)
Calisher, C. H., Childs, J. E., Field, H. E., Holmes, K. V. & Schountz, T. Bats: important reservoir hosts of emerging viruses. Clin. Microbiol. Rev. 19, 531–545 (2006)
Scherf, B. D. World Watch List for Domestic Animal Diversity. 3rd edn, (Food and Agriculture Organization of the United Nations, 2000)
IUCN. The IUCN Red List of Threatened Species. Version 2014.1, http://www.iucnredlist.org (2014)
Epstein, J. H. & Price, J. T. The significant but understudied impact of pathogen transmission from humans to animals. Mt. Sinai J. Med. 76, 448–455 (2009)
Messenger, A. M., Barnes, A. N. & Gray, G. C. Reverse zoonotic disease transmission (zooanthroponosis): a systematic review of seldom-documented human biological threats to animals. PLoS One 9, e89055 (2014)
Jones, K. E. et al. PanTHERIA: a species-level database of life history, ecology, and geography of extant and recently extinct mammals. Ecology 90, 2648 (2009)
de Magalhães, J. P. & Costa, J. A database of vertebrate longevity records and their relation to other life-history traits. J. Evol. Biol. 22, 1770–1774 (2009)
Cooper, N., Kamilar, J. M. & Nunn, C. L. Host longevity and parasite species richness in mammals. PLoS One 7, e42190 (2012)
Felsenstein, J. Phylogenies and the comparative method. Am. Nat. 125, 1–15 (1985)
Blomberg, S. P., Garland, T., Jr & Ives, A. R. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution 57, 717–745 (2003)
Grafen, A. The phylogenetic regression. Phil. Trans. R. Soc. Lond. B 326, 119–157 (1989)
Diniz-Filho, J. A. F. et al. On the selection of phylogenetic eigenvectors for ecological analyses. Ecography 35, 239–249 (2012)
Diniz-Filho, J. A. F., de Sant’Ana, C. E. R. & Bini, L. M. An eigenvector method for estimating phylogenetic inertia. Evolution 52, 1247–1262 (1998)
Bininda-Emonds, O. R. P. et al. The delayed rise of present-day mammals. Nature 446, 507–512 (2007)
Paradis, E ., Claude, J . & Strimmer, K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20, 289–290 (2004)
Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004)
Stamatakis, A., Hoover, P. & Rougemont, J. A rapid bootstrap algorithm for the RAxML Web servers. Syst. Biol. 57, 758–771 (2008)
Cuthill, J. H. & Charleston, M. A. A simple model explains the dynamics of preferential host switching among mammal RNA viruses. Evolution 67, 980–990 (2013)
Poulin, R., Krasnov, B. R. & Mouillot, D. Host specificity in phylogenetic and geographic space. Trends Parasitol. 27, 355–361 (2011)
Poulin, R. & Mouillot, D. Parasite specialization from a phylogenetic perspective: a new index of host specificity. Parasitology 126, 473–480 (2003)
R Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/ (2014)
Marra, G. & Wood, S. N. Practical variable selection for generalized additive models. Comput. Stat. Data Anal. 55, 2372–2387 (2011)
Pawitan, Y. In All Likelihood: Statistical Modelling and Inference Using Likelihood. (Oxford University Press, 2001)
Holt, B. G. et al. An update of Wallace’s zoogeographic regions of the world. Science 339, 74–78 (2013)
QGIS Geographic Information System. Open Source Geospatial Foundation Project http://www.qgis.org/ (2016)
Goldewijk, K. K., Beusen, A., van Drecht, G. & de Vos, M. The HYDE 3.1 spatially explicit database of human-induced global land-use change over the past 12,000 years. Glob. Ecol. Biogeogr. 20, 73–86 (2011)
raster: Geographic Data Analysis and Modeling version 2.3-40 https://cran.r-project.org/package=raster (2015)
sp: Classes and Methods for Spatial Data version 1.2-1 https://cran.r-project.org/package=sp (2015)
Acknowledgements
This research was supported by the United States Agency for International Development (USAID) Emerging Pandemic Threats PREDICT program; and NIH NIAID awards R01AI079231 and R01AI110964. The authors thank C. N. Basaraba, J. Baxter, L. Brierley, E. A. Hagan, J. Levinson, E. H. Loh, L. Mendiola, N. Wale and A. R. Willoughby for assistance with data collection, and B. M. Bolker, A. R. Ives, K. E. Jones, C. K. Johnson, A. M. Kilpatrick, J. A. K. Mazet and M. E. J. Woolhouse for comments.
Author information
Authors and Affiliations
Contributions
K.J.O., T.L.B. and P.D. designed the study and supervised the collection of data. N.R., P.R.H. and K.J.O. designed the statistical approach, wrote the code, and generated figures. K.J.O. performed phylogenetic analyses. C.Z.-T. performed spatial analyses. All authors were involved in writing the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Reviewer Information Nature thanks J. Dushoff and the other anonymous reviewer(s) for their contribution to the peer review of this work.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Figure 1 Conceptual model of zoonotic spillover, viral richness, and summary of models.
a, Conceptual model of zoonotic spillover showing primary risk factors examined, colour-coded according to generalized additive models used. b, Conceptual model of observed, predicted, and actual viral richness in mammals. c, GAMs used in our study to address specific components of a and b, colour-coded by model. Variables listed with ‘or’ under each GAM covaried and were provided as competing terms in model selection, and those in bold were included in the best-fit model using all host–virus associations. Significant variables from each best-fit GAM are noted with an asterisk. Zoonotic viral spillover first depends on the underlying total viral richness in mammal populations and the ecological, taxonomic, and life-history traits that govern this diversity (GAM 1). Second, host- and virus-specific factors may facilitate viral spillover. We examine the relative importance of host phylogenetic distance to humans, ecological opportunity for contact, or other species-specific life-history and taxonomic traits (GAM 2), and identify viral traits associated with a higher likelihood of an observed virus being zoonotic (GAM 3). We estimate the total and zoonotic viral richness per host species using GAMs 1 and 2, and calculate the missing viruses and missing zoonoses under a scenario of increased research effort (b, Methods). Owing to imperfect surveillance in both humans and wildlife and biases in viral detection, there may be uncertainty in the exact proportion of viruses that are zoonotic (b, light grey), and also between the actual, or true, viral richness (dotted lines) and the predicted maximum viral richness per host (dashed line).
Extended Data Figure 2 Heat map of observed total viral richness by mammalian order and viral family.
Dataset includes 754 mammalian species and 586 unique ICTV recognized viral species. Heat map aggregated by rows and columns to group taxa with similar levels of observed viral richness.
Extended Data Figure 3 Global distribution of viral and host species richness for all wild mammals.
a, Observed total viral richness (for n = 576 host spp.); b, predicted total viral richness given maximum research effort; c, missing viruses or predicted minus observed total viral richness; d, observed zoonotic viral richness (n = 584); e, predicted zoonotic viral richness given maximum research effort; f, missing zoonoses or predicted minus observed zoonotic viral richness (same as included in Fig. 3a); g, global mammal species richness (n = 5,290); h, mammal richness for species in our database (n = 753); i, mammal species with no described viruses in the literature. Warmer colours (larger values) in panels c and f highlight areas predicted to be of greatest value for discovering novel viruses or novel viral zoonoses, respectively, in mammals. Red/pink colours in panel i highlight areas with poor viral surveillance in mammal species to date. Hatched regions represent areas where model predictions deviate systematically for the collection of species in that grid cell (see Methods).
Extended Data Figure 4 Global distribution of viral and host species richness for wild carnivores (order Carnivora).
a, Observed total viral richness (for n = 55 host spp.); b, predicted total viral richness given maximum research effort; c, missing viruses or predicted minus observed total viral richness; d, observed zoonotic viral richness (n = 55); e, predicted zoonotic viral richness given maximum research effort; f, missing zoonoses or predicted minus observed zoonotic viral richness (same as included in Fig. 3b); g, global host species richness for Carnivora (n = 276); h, host species richness for Carnivora in our database (n = 79); i, species of the order Carnivora with no described viruses in the literature. Warmer colours (larger values) in c and f highlight areas predicted to be of greatest value for discovering novel viruses or novel viral zoonoses, respectively, in carnivores. Red/pink colours in panel i highlight areas with poor viral surveillance in carnivore species to date. Hatched regions represent areas where model predictions deviate systematically for the collection of species in that grid cell (see Methods).
Extended Data Figure 5 Global distribution of viral and host species richness for wild even-toed ungulates (order Cetartiodactyla).
a, Observed total viral richness (for n = 70 host spp.); b, predicted total viral richness given maximum research effort; c, missing viruses or predicted minus observed total viral richness; d, observed zoonotic viral richness (n = 70); e, predicted zoonotic viral richness given maximum research effort; f, missing zoonoses or predicted minus observed zoonotic viral richness (same as included in Fig. 3c); g, global host species richness for Cetartiodactyla (n = 229); h, host species richness for Cetartiodactyla in our database (n = 105); i, species of the order Cetartiodactyla with no described viruses in the literature. Warmer colours (larger values) in c and f highlight areas predicted to be of greatest value for discovering novel viruses or novel viral zoonoses, respectively, in even-toed ungulates. Red/pink colours in panel i highlight areas with poor viral surveillance in even-toed ungulates species to date. Hatched regions represent areas where model predictions deviate systematically for the collection of species in that grid cell (see Methods).
Extended Data Figure 6 Global distribution of viral and host species richness for bats (order Chiroptera).
a, Observed total viral richness (for n = 156 host spp.); b, predicted total viral richness given maximum research effort; c, missing viruses or predicted minus observed total viral richness; d, observed zoonotic viral richness (n = 157); e, predicted zoonotic viral richness given maximum research effort; f, missing zoonoses or predicted minus observed zoonotic viral richness (same as included in Fig. 3d); g, global host species richness for Chiroptera (n = 1117); h, host species richness for Chiroptera in our database (n = 192); i, species of the order Chiroptera with no described viruses in the literature. Warmer colours (larger values) in c and f highlight areas predicted to be of greatest value for discovering novel viruses or novel viral zoonoses, respectively, in bats. Red/pink colours in panel i highlight areas with poor viral surveillance in bat species to date. Hatched regions represent areas where model predictions deviate systematically for the collection of species in that grid cell (see Methods).
Extended Data Figure 7 Global distribution of viral and host species richness for primates (order Primates).
a, Observed total viral richness (for n = 71 host spp.); b, predicted total viral richness given maximum research effort; c, missing viruses or predicted minus observed total viral richness; d, observed zoonotic viral richness (n = 73); e, predicted zoonotic viral richness given maximum research effort; f, missing zoonoses or predicted minus observed zoonotic viral richness (same as included in Fig. 3e); g, global host species richness for Primates (n = 400); h, host species richness for Primates in our database (n = 98); i, primate species with no described viruses in the literature. Warmer colours (larger values) in c and f highlight areas predicted to be of greatest value for discovering novel viruses or novel viral zoonoses, respectively, in primates. Red/pink colours in panel i highlight areas with poor viral surveillance in primate species to date. Hatched regions represent areas where model predictions deviate systematically for the collection of species in that grid cell (see Methods).
Extended Data Figure 8 Global distribution of viral and host species richness for rodents (order Rodentia).
a, Observed total viral richness (for n = 178 host spp.); b, predicted total viral richness given maximum research effort; c, missing viruses or predicted minus observed total viral richness; d, observed zoonotic viral richness (n = 183); e, predicted zoonotic viral richness given maximum research effort; f, missing zoonoses or predicted minus observed zoonotic viral richness (same as included in Fig. 3f); g, global host species richness for Rodentia (n = 2206); h, host species richness for Rodentia in our database (n = 221); i, rodent species with no described viruses in the literature. Warmer colours (larger values) in c and f highlight areas predicted to be of greatest value for discovering novel viruses or novel viral zoonoses, respectively, in wild rodents. Red/pink colours in panel i highlight areas with poor viral surveillance in rodent species to date. Hatched regions represent areas where model predictions deviate systematically for the collection of species in that grid cell (see Methods).
Extended Data Figure 9 Order-level phylogenies showing residuals from zoonoses model.
a–e, Subtrees from cytochrome b maximum likelihood phylogeny for 558 mammal species (constrained to order-level topology of mammal supertree) for bats (a), carnivores (b), even-toed ungulates (c), rodents (d) and primates (e). Species included have at least one described virus association and available genetic data. Wildlife species names and terminal branches are colour-coded by the residuals (predicted minus observed) from the best-fit GAM to predict the number of zoonotic viruses using all data. Species with residual values between −1 and 1 (black) are accurately predicted within one virus. Warm colours represent species with positive residuals (orange >1 to 3; red >3). Cool colours represent species with negative residuals (green <−1 to −3; blue <−3). Marine mammals, domestic animals, and species with missing data and not included in the best-fit models are shown in grey.
Supplementary information
Supplementary Information
This file contains a Supplementary Discussion and additional references. (PDF 199 kb)
Supplementary Data
This file contains Supplementary Table 1. This file was added on 23 August 2017. (CSV 17 kb)
Supplementary Data
This file contains Supplementary Table 2. This file was added on 23 August 2017. (CSV 3 kb)
Supplementary Data
This file contains Supplementary Table 3. This file was added on 23 August 2017. (CSV 7 kb)
Supplementary Data
This file contains Supplementary Table 4. This file was added on 23 August 2017. (CSV 82 kb)
Rights and permissions
About this article
Cite this article
Olival, K., Hosseini, P., Zambrana-Torrelio, C. et al. Host and viral traits predict zoonotic spillover from mammals. Nature 546, 646–650 (2017). https://doi.org/10.1038/nature22975
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nature22975
This article is cited by
-
Receptor-binding proteins from animal viruses are broadly compatible with human cell entry factors
Nature Microbiology (2025)
-
Pathogens and planetary change
Nature Reviews Biodiversity (2025)
-
Host ecology and phylogeny shape the temporal dynamics of social bee viromes
Nature Communications (2025)
-
Bat genomes illuminate adaptations to viral tolerance and disease resistance
Nature (2025)
-
Identification and genetic characterization of five novel bat coronaviruses from Yunnan, China
BMC Veterinary Research (2024)