Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Letter
  • Published:

Host and viral traits predict zoonotic spillover from mammals

An Erratum to this article was published on 23 August 2017

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

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Figure 1: Observed viral richness in mammals.
Figure 2: Host traits that predict total viral richness (top row) and proportion of zoonotic viruses (bottom row) per wild mammal species.
Figure 3: Global distribution of the predicted number of ‘missing zoonoses’ by order.
Figure 4: Traits that predict zoonotic potential of a virus.

Similar content being viewed by others

References

  1. Woolhouse, M. E. J. & Gowtage-Sequeria, S. Host range and emerging and reemerging pathogens. Emerg. Infect. Dis. 11, 1842–1847 (2005)

    Article  PubMed  PubMed Central  Google Scholar 

  2. Jones, K. E. et al. Global trends in emerging infectious diseases. Nature 451, 990–993 (2008)

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  3. 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)

    Article  CAS  Google Scholar 

  4. Morse, S. S. et al. Prediction and prevention of the next pandemic zoonosis. Lancet 380, 1956–1965 (2012)

    Article  PubMed  PubMed Central  Google Scholar 

  5. Parrish, C. R. et al. Cross-species virus transmission and the emergence of new epidemic diseases. Microbiol. Mol. Biol. Rev. 72, 457–470 (2008)

    Article  PubMed  PubMed Central  Google Scholar 

  6. Lipsitch, M. et al. Viral factors in influenza pandemic risk assessment. eLife 5, e18491 (2016)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. 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)

    Article  PubMed  Google Scholar 

  8. Woolhouse, M. E., Haydon, D. T. & Antia, R. Emerging pathogens: the epidemiology and evolution of species jumps. Trends Ecol. Evol. 20, 238–244 (2005)

    Article  PubMed  PubMed Central  Google Scholar 

  9. Cooper, N. et al. Phylogenetic host specificity and understanding parasite sharing in primates. Ecol. Lett. 15, 1370–1377 (2012)

    Article  PubMed  Google Scholar 

  10. 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)

    Article  Google Scholar 

  11. 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)

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  12. 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)

    Article  Google Scholar 

  13. 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)

    Article  PubMed  Google Scholar 

  14. Streicker, D. G. et al. Host phylogeny constrains cross-species emergence and establishment of rabies virus in bats. Science 329, 676–679 (2010)

    Article  CAS  ADS  PubMed  Google Scholar 

  15. 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)

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  16. 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)

  17. 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)

    Article  Google Scholar 

  18. Levinson, J. et al. Targeting surveillance for zoonotic virus discovery. Emerg. Infect. Dis. 19, 743–747 (2013)

    Article  PubMed  PubMed Central  Google Scholar 

  19. Morse, S. S. in Emerging Viruses (ed. Morse, S. S. ) 10–28 (Oxford University Press, 1993)

  20. 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)

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  21. Parker, I. M. et al. Phylogenetic structure and host abundance drive disease pressure in communities. Nature 520, 542–544 (2015)

    Article  CAS  ADS  PubMed  Google Scholar 

  22. 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)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Ge, X.-Y. et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 503, 535–538 (2013)

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  24. 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)

    Article  CAS  PubMed  Google Scholar 

  25. Anthony, S. J. et al. A strategy to estimate unknown viral diversity in mammals. MBio 4, e00598–13 (2013)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Drexler, J. F. et al. Bats host major mammalian paramyxoviruses. Nat. Commun. 3, 796 (2012)

    Article  ADS  CAS  PubMed  Google Scholar 

  27. 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)

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  28. 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)

    Article  PubMed  Google Scholar 

  29. 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)

    Article  Google Scholar 

  30. 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)

    Google Scholar 

  31. 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)

    Article  PubMed  PubMed Central  Google Scholar 

  32. Scherf, B. D. World Watch List for Domestic Animal Diversity. 3rd edn, (Food and Agriculture Organization of the United Nations, 2000)

  33. IUCN. The IUCN Red List of Threatened Species. Version 2014.1, http://www.iucnredlist.org (2014)

  34. 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)

    Article  PubMed  PubMed Central  Google Scholar 

  35. 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)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  36. 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)

    Article  Google Scholar 

  37. 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)

    Article  PubMed  Google Scholar 

  38. Cooper, N., Kamilar, J. M. & Nunn, C. L. Host longevity and parasite species richness in mammals. PLoS One 7, e42190 (2012)

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  39. Felsenstein, J. Phylogenies and the comparative method. Am. Nat. 125, 1–15 (1985)

    Article  Google Scholar 

  40. 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)

    Article  PubMed  Google Scholar 

  41. Grafen, A. The phylogenetic regression. Phil. Trans. R. Soc. Lond. B 326, 119–157 (1989)

    Article  CAS  ADS  Google Scholar 

  42. Diniz-Filho, J. A. F. et al. On the selection of phylogenetic eigenvectors for ecological analyses. Ecography 35, 239–249 (2012)

    Article  Google Scholar 

  43. 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)

    Article  PubMed  Google Scholar 

  44. Bininda-Emonds, O. R. P. et al. The delayed rise of present-day mammals. Nature 446, 507–512 (2007)

    Article  CAS  ADS  PubMed  Google Scholar 

  45. Paradis, E ., Claude, J . & Strimmer, K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20, 289–290 (2004)

    Article  CAS  PubMed  Google Scholar 

  46. Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Stamatakis, A., Hoover, P. & Rougemont, J. A rapid bootstrap algorithm for the RAxML Web servers. Syst. Biol. 57, 758–771 (2008)

    Article  PubMed  Google Scholar 

  48. 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)

    Article  PubMed  PubMed Central  Google Scholar 

  49. Poulin, R., Krasnov, B. R. & Mouillot, D. Host specificity in phylogenetic and geographic space. Trends Parasitol. 27, 355–361 (2011)

    Article  PubMed  Google Scholar 

  50. Poulin, R. & Mouillot, D. Parasite specialization from a phylogenetic perspective: a new index of host specificity. Parasitology 126, 473–480 (2003)

    Article  CAS  PubMed  Google Scholar 

  51. R Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org/ (2014)

  52. Marra, G. & Wood, S. N. Practical variable selection for generalized additive models. Comput. Stat. Data Anal. 55, 2372–2387 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  53. Pawitan, Y. In All Likelihood: Statistical Modelling and Inference Using Likelihood. (Oxford University Press, 2001)

  54. Holt, B. G. et al. An update of Wallace’s zoogeographic regions of the world. Science 339, 74–78 (2013)

    Article  CAS  ADS  PubMed  Google Scholar 

  55. QGIS Geographic Information System. Open Source Geospatial Foundation Project http://www.qgis.org/ (2016)

  56. 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)

    Article  Google Scholar 

  57. raster: Geographic Data Analysis and Modeling version 2.3-40 https://cran.r-project.org/package=raster (2015)

  58. sp: Classes and Methods for Spatial Data version 1.2-1 https://cran.r-project.org/package=sp (2015)

Download references

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

Authors

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

Correspondence to Kevin J. Olival or Peter Daszak.

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.

ae, 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.

Extended Data Table 1 Summary of best-fit GAMs for total and zoonotic viral richness per wild mammal species, and probability of a virus being zoonotic

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)

PowerPoint slides

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nature22975

This article is cited by

Search

Quick links

Nature Briefing Microbiology

Sign up for the Nature Briefing: Microbiology newsletter — what matters in microbiology research, free to your inbox weekly.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing: Microbiology