Ecology

Liu, W. et al. African origin of the malaria parasite Plasmodium vivax. Nat. Commun. 5, 3346 (2014).PubMed 

Google Scholar 
Liu, W. et al. Multigenomic delineation of Plasmodium species of the Laverania subgenus infecting wild-living chimpanzees and gorillas. Genome Biol. Evolution 8, 1929–1939 (2016).CAS 

Google Scholar 
Liu, W. et al. Single genome amplification and direct amplicon sequencing of Plasmodium spp. DNA from ape fecal specimens. Protocol Exchange 1–14 (2010).Liu, W. et al. Wild bonobos host geographically restricted malaria parasites including a putative new Laverania species. Nat. Commun. 8, 1635 (2017).PubMed 
PubMed Central 

Google Scholar 
Prugnolle, F. et al. African great apes are natural hosts of multiple related malaria species, including Plasmodium falciparum. Proc. Natl Acad. Sci. USA 107, 1458–1463 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Sharp, P. M., Plenderleith, L. J. & Hahn, B. H. Ape origins of human malaria. Annu. Rev. Microbiol. 74, 39–63 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Liu, W. et al. Origin of the human malaria parasite Plasmodium falciparum in gorillas. Nature 467, 420–425 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Otto, T. D. et al. Genomes of all known members of a Plasmodium subgenus reveal paths to virulent human malaria. Nat. Microbiol. 3, 687–697 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
Boundenga, L. et al. Diversity of malaria parasites in great apes in Gabon. Malar. J. 14, 1–8 (2015).CAS 

Google Scholar 
Délicat-Loembet, L. et al. No evidence for ape Plasmodium infections in humans in gabon. Plos One 10, e0126933 (2015).PubMed 
PubMed Central 

Google Scholar 
Sundararaman, S. A. et al. Plasmodium falciparum-like parasites infecting wild apes in southern Cameroon do not represent a recurrent source of human malaria. Proc. Natl Acad. Sci. USA 110, 7020–7025 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Junker, J. et al. Recent decline in suitable environmental conditions for African great apes. Diversity Distrib. 18, 1077–1091 (2012).
Google Scholar 
de Nys, H. M. et al. Age-related effects on malaria parasite infection in wild chimpanzees. Biol. Lett. 9, 20121160 (2013).PubMed 
PubMed Central 

Google Scholar 
de Nys, H. M. et al. Malaria parasite detection increases during pregnancy in wild chimpanzees. Malar. J. 13, 413 (2014).PubMed 
PubMed Central 

Google Scholar 
Kaiser, M. et al. Wild chimpanzees infected with 5 Plasmodium species. Emerg. Infect. Dis. 16, 1956–1959 (2010).PubMed 
PubMed Central 

Google Scholar 
Paupy, C. et al. Anopheles moucheti and Anopheles vinckei are candidate vectors of ape Plasmodium parasites, including Plasmodium praefalciparum in Gabon. PLoS ONE 8, e57294 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Makanga, B. et al. Ape malaria transmission and potential for ape-to-human transfers in Africa. Proc. Natl Acad. Sci. USA 113, 5329–5334 (2016).Loy, D. E. et al. Investigating zoonotic infection barriers to ape Plasmodium parasites using faecal DNA analysis. Int. J. Parasitol. 48, 531–542 (2018).Martin, M., Rayner, J., Gagneux, P., Barnwell, J. & Varki, A. Evolution of human–chimpanzee differences in malaria susceptibility: Relationship to human genetic loss of N-glycolylneuraminic acid. Proc. Natl Acad. Sci. USA 102, 12819–12824 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
Scully, E. J., Kanjee, U. & Duraisingh, M. T. Molecular interactions governing host-specificity of blood stage malaria parasites. Curr. Opin. Microbiol. 40, 21–31 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Sundararaman, S. A. et al. Genomes of cryptic chimpanzee Plasmodium species reveal key evolutionary events leading to human malaria. Nat. Commun. 7, 11078 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wanaguru, M., Liu, W., Hahn, B. H., Rayner, J. C. & Wright, G. J. RH5-Basigin interaction plays a major role in the host tropism of Plasmodium falciparum. Proc. Natl Acad. Sci. USA 110, 20735–20740 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ngoubangoye, B. et al. The host specificity of ape malaria parasites can be broken in confined environments. Int. J. Parasitol. 46, 737–744 (2016).PubMed 

Google Scholar 
Mapua, M. I. et al. A comparative molecular survey of malaria prevalence among Eastern chimpanzee populations in Issa Valley (Tanzania) and Kalinzu (Uganda). Malar. J. 15, 423 (2016).PubMed 
PubMed Central 

Google Scholar 
Wu, D. F. et al. Seasonal and inter-annual variation of malaria parasite detection in wild chimpanzees. Malar. J. 17, 1–5 (2018).CAS 

Google Scholar 
Craig, M., le Sueur, D. & Snow, B. A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol. Today 15, 105–111 (1999).CAS 
PubMed 

Google Scholar 
Mordecai, E. A. et al. Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol. Lett. 16, 22–30 (2013).PubMed 

Google Scholar 
Paaijmans, K. P. et al. Influence of climate on malaria transmission depends on daily temperature variation. Proc. Natl Acad. Sci. USA 107, 15135–15139 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Parham, P. E. & Michael, E. Modeling the effects of weather and climate change on malaria transmission. Environ. Health Perspect. 118, 620–626 (2010).PubMed 

Google Scholar 
LaPointe, D. A., Goff, M. L. & Atkinson, C. T. Thermal constraints to the sporogonic development and altitudinal distribution of avian malaria Plasmodium relictum in Hawai’i. J. Parasitol. 96, 318–324 (2010).PubMed 

Google Scholar 
Vanderberg, J. P. & Yoeli, M. Effects of temperature on sporogonic development of Plasmodium berghei. J. Parasitol. 52, 559–564 (1966).Macdonald, G. The Epidemiology and Control of Malaria (Oxford University Press, 1957).Ryan, S. J. et al. Mapping physiological suitability limits for malaria in Africa under climate change. Vector-Borne Zoonotic Dis. 15, 718–725 (2015).PubMed 
PubMed Central 

Google Scholar 
Gemperli, A. et al. Mapping malaria transmission in West and Central Africa. Tropical Med. Int. Health 11, 1032–1046 (2006).
Google Scholar 
Gething, P. W. et al. Modelling the global constraints of temperature on transmission of Plasmodium falciparum and P. vivax. Parasites Vectors 4, 92 (2011).PubMed 
PubMed Central 

Google Scholar 
Weiss, D. J. et al. Air temperature suitability for Plasmodium falciparum malaria transmission in Africa 2000–2012: a high-resolution spatiotemporal prediction. Malar. J. 13, 171 (2014).PubMed 
PubMed Central 

Google Scholar 
Lyons, C. L., Coetzee, M. & Chown, S. L. Stable and fluctuating temperature effects on the development rate and survival of two malaria vectors, Anopheles arabiensis and Anopheles funestus. Parasites Vectors 6, 104 (2013).PubMed 
PubMed Central 

Google Scholar 
Paaijmans, K. P., Wandago, M. O., Githeko, A. K. & Takken, W. Unexpected high losses of Anopheles gambiae larvae due to rainfall. PLoS One 2, e1146 (2007).PubMed 
PubMed Central 

Google Scholar 
Faust, C. & Dobson, A. P. Primate malarias: diversity, distribution and insights for zoonotic Plasmodium. One Health 1, 66–75 (2015).PubMed 
PubMed Central 

Google Scholar 
Tucker Lima, J. M., Vittor, A., Rifai, S. & Valle, D. Does deforestation promote or inhibit malaria transmission in the Amazon? A systematic literature review and critical appraisal of current evidence. Philos. Trans. R. Soc. Lond. Ser. B, Biol. Sci. 372, 20160125 (2017).
Google Scholar 
Borner, J. et al. Phylogeny of haemosporidian blood parasites revealed by a multi-gene approach. Mol. Phylogenetics Evolution 94, 221–231 (2016).CAS 

Google Scholar 
Emery Thompson, M., Muller, M. N., Machanda, Z. P., Otali, E. & Wrangham, R. W. The Kibale Chimpanzee Project: over thirty years of research, conservation, and change. Biol. Conserv. 252, 108857 (2020).
Google Scholar 
Langergraber, K. E., Mitani, J. C. & Vigilant, L. The limited impact of kinship on cooperation in wild chimpanzees. Proc. Natl Acad. Sci. USA 104, 7786–7790 (2007).CAS 
PubMed 
PubMed Central 

Google Scholar 
Arandjelovic, M. et al. Two-step multiplex polymerase chain reaction improves the speed and accuracy of genotyping using DNA from noninvasive and museum samples. Mol. Ecol. Resour. 9, 28–36 (2009).CAS 
PubMed 

Google Scholar 
Herbert, A. et al. Malaria-like symptoms associated with a natural Plasmodium reichenowi infection in a chimpanzee. Malar. J. 14, 220 (2015).PubMed 
PubMed Central 

Google Scholar 
Torres, J. R. Therapy of Infectious Diseases 597–613 (2003).Trampuz, A., Jereb, M., Muzlovic, I. & Prabhu, R. M. Clinical review: severe malaria. Crit. Care 7, 315 (2003).PubMed 
PubMed Central 

Google Scholar 
Akim, N. I. et al. Dynamics of P. falciparum gametocytemia in symptomatic patients in an area of intense perennial transmission in Tanzania. Am. J. Tropical Med. Hyg. 63, 199–203 (2000).CAS 

Google Scholar 
Mackinnon, M. J. & Read, A. F. Genetic relationships between parasite virulence and transmission in the rodent malaria Plasmodium chabaudi. Evolution 53, 689–703 (1999).PubMed 

Google Scholar 
Huelsenbeck, J. P. & Ronquist, F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17, 754–755 (2001).CAS 
PubMed 

Google Scholar 
Prugnolle, F. et al. African monkeys are infected by Plasmodium falciparum nonhuman primate-specific strains. Proc. Natl Acad. Sci. USA 108, 11948–11953 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ayouba, A. et al. Ubiquitous Hepatocystis infections, but no evidence of Plasmodium falciparum-like malaria parasites in wild greater spot-nosed monkeys (Cercopithecus nictitans). Int. J. Parasitol. 42, 709–713 (2012).PubMed 
PubMed Central 

Google Scholar 
Martinsen, E. S., Perkins, S. L. & Schall, J. J. A three-genome phylogeny of malaria parasites (Plasmodium and closely related genera): Evolution of life-history traits and host switches. Mol. Phylogenetics Evolution 47, 261–273 (2008).CAS 

Google Scholar 
Thurber, M. I. et al. Co-infection and cross-species transmission of divergent Hepatocystis lineages in a wild African primate community. Int. J. Parasitol. 43, 613–619 (2013).PubMed 
PubMed Central 

Google Scholar 
Baayen, R. H. Analyzing Linguistic Data: A Practical Introduction to Statistics (Cambridge University Press, 2008).Stanisic, D. I. et al. Acquisition of antibodies against Plasmodium falciparum merozoites and malaria immunity in young children and the influence of age, force of infection, and magnitude of response. Infect. Immun. 83, 646–660 (2015).PubMed 
PubMed Central 

Google Scholar 
Taylor, R. R., Allen, S. J., Greenwood, B. M. & Riley, E. M. IgG3 antibodies to Plasmodium falciparum merozoite surface protein 2 (MSP2): increasing prevalence with age and association with clinical immunity to malaria. Am. J. Tropical Med. Hyg. 58, 406–413 (1998).CAS 

Google Scholar 
World Malaria Report (World Health Organization, 2015).Shaman, J. Letter to the Editor: Caution needed when using gridded meteorological data products for analyses in Africa. Eur. Surveill. 19, 20930 (2014).
Google Scholar 
Tatem, A. J., Goetz, S. J. & Hay, S. I. Terra and Aqua: new data for epidemiology and public health. Int. J. Appl. Earth Observation Geoinf. 6, 33–46 (2004).
Google Scholar 
Adler, R. F. et al. The version-2 global precipitation climatology project (GPCP) monthly precipitation analysis (1979–present). J. Hydrometeorol. 4, 1147–1167 (2003).
Google Scholar 
Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850–853 (2013).CAS 
PubMed 

Google Scholar 
Carter, R. & Mendis, K. N. Evolutionary and historical aspects of the burden of malaria. Clin. Microbiol. Rev. 15, 564–594 (2002).PubMed 
PubMed Central 

Google Scholar 
Kwiatkowski, D. P. How malaria has affected the human genome and what human genetics can teach us about malaria. Am. J. Hum. Genet. 77, 171–192 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
Tarello, W. A fatal Plasmodium reichenowi infection in a chimpanzee? Rev. de. Med. Veterinaire 156, 503–505 (2005).
Google Scholar 
Taylor, D. W. et al. Parasitologic and immunologic studies of experimental Plasmodium falciparum infection in nonsplenectomized chimpanzees (Pan troglodytes). Am. J. Tropical Med. Hyg. 34, 36–44 (1985).CAS 

Google Scholar 
Krief, S., Martin, M., Grellier, P., Kasenene, J. & Sevenet, T. Novel antimalarial compounds isolated in a survey of self-medicative behavior of wild chimpanzees in Uganda. Antimicrobial Agents Chemother. 48, 3196–3199 (2004).CAS 

Google Scholar 
Cox-Singh, J. et al. Plasmodium knowlesi malaria in humans is widely distributed and potentially life threatening. Clin. Infect. Dis. 46, 165–171 (2008).CAS 
PubMed 

Google Scholar 
Singh, B. & Daneshvar, C. Human infections and detection of Plasmodium knowlesi. Clin. Microbiol. Rev. 26, 165–184 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Brasil, P. et al. Outbreak of human malaria caused by Plasmodium simium in the Atlantic Forest in Rio de Janeiro: a molecular epidemiological investigation. Lancet Global Health 5, e1038–e1046 (2017).Krief, S. et al. On the diversity of malaria parasites in African apes and the origin of Plasmodium falciparum from bonobos. PLoS Pathog. 6, e1000765 (2010).Pacheco, M. A., Cranfield, M., Cameron, K. & Escalante, A. A. Malarial parasite diversity in chimpanzees: the value of comparative approaches to ascertain the evolution of Plasmodium falciparum antigens. Malar. J. 12, 328 (2013).PubMed 
PubMed Central 

Google Scholar 
Etienne, L. et al. Noninvasive follow-up of simian immunodeficiency virus infection in wild-living nonhabituated western lowland gorillas in Cameroon. J. Virol. 86, 9760–9772 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Keele, B. F. et al. Chimpanzee reservoirs of pandemic and nonpandemic HIV-1. Science 313, 523–526 (2006).CAS 
PubMed 
PubMed Central 

Google Scholar 
Keele, B. F. et al. Increased mortality and AIDS-like immunopathology in wild chimpanzees infected with SIVcpz. Nature 460, 515–519 (2009).CAS 
PubMed 
PubMed Central 

Google Scholar 
Li, Y. et al. Eastern chimpanzees, but not bonobos, represent a simian immunodeficiency virus reservoir. J. Virol. 86, 10776–10791 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Neel, C. et al. Molecular epidemiology of simian immunodeficiency virus infection in wild-living gorillas. J. Virol. 84, 1464–1476 (2010).CAS 
PubMed 

Google Scholar 
Rudicell, R. S. et al. Impact of simian immunodeficiency virus infection on chimpanzee population dynamics. PLoS Pathog. 6, 1–17 (2010).
Google Scholar 
Darriba, D., Taboada, G. L., Doallo, R. & Posada, D. jModelTest 2: more models, new heuristics and parallel computing. Nat. Methods 9, 772 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Bates, D. & Maechler, M. Lme4: linear mixed-effects models using s4 classes. Cran R Project Website (2010). More

Morens, D. M. et al. The origin of COVID-19 and why it matters. Am. J. Trop. Med. Hyg. 103, 955–959 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pierson, T. C. & Diamond, M. S. The emergence of Zika virus and its new clinical syndromes. Nature 560, 573–581 (2018).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Gates, B. The next epidemic—Lessons from Ebola., https://doi.org/10.1056/NEJMp1502918 (2015).World Health Organization. Lassa fever research and development (R&D) roadmap. https://www.who.int/publications/m/item/lassa-fever-research-and-development-(r-d)-roadmap (2018).World Health Organization. Prioritizing diseases for research and development in emergency contexts. https://www.who.int/activities/prioritizing-diseases-for-research-and-development-in-emergency-contexts.Akpede, G. O. et al. Caseload and case fatality of Lassa fever in Nigeria, 2001–2018: A specialist center’s experience and its implications. Front. Public Health 7, https://doi.org/10.3389/fpubh.2019.00170 (2019).Eberhardt, K. A. et al. Ribavirin for the treatment of Lassa fever: A systematic review and meta-analysis. Int. J. Infect. Dis. 87, 15–20 (2019).CAS 
PubMed 
Article 

Google Scholar 
Lukashevich, I. S., Paessler, S. & de la Torre, J. C. Lassa virus diversity and feasibility for universal prophylactic vaccine. F1000Res 8, https://doi.org/10.12688/f1000research.16989.1 (2019).Purushotham, J., Lambe, T. & Gilbert, S. C. Vaccine platforms for the prevention of Lassa fever. Immunol. Lett. 215, 1–11 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mateo, M. et al. A single-shot Lassa vaccine induces long-term immunity and protects cynomolgus monkeys against heterologous strains. Sci. Transl. Med. 13, eabf6348 (2021).CAS 
PubMed 
Article 

Google Scholar 
McCormick, J. B. et al. Lassa Fever. N. Engl. J. Med. 314, 20–26 (1986).CAS 
PubMed 
Article 

Google Scholar 
Bell-Kareem, A. R. & Smither, A. R. Epidemiology of Lassa fever. in 1–23 (Springer, 2021). https://doi.org/10.1007/82_2021_234.Nigeria Centre for Disease Control. https://ncdc.gov.ng/diseases/sitreps/?cat=5&name=An%20update%20of%20Lassa%20fever%20outbreak%20in%20Nigeria.Manning, J. T., Forrester, N. & Paessler, S. Lassa virus isolates from Mali and the Ivory Coast represent an emerging fifth lineage. Front. Microbiol. 6, https://doi.org/10.3389/fmicb.2015.01037 (2015).Dzotsi, E. K. et al. The first cases of Lassa fever in Ghana. Ghana. Med. J. 46, 166–170 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Patassi, A. A. et al. Emergence of Lassa fever disease in northern Togo: Report of two cases in Oti District in 2016. Case Rep. Infect. Dis. 2017, 8242313 (2017).PubMed 
PubMed Central 

Google Scholar 
Yadouleton, A. et al. Lassa fever in Benin: Description of the 2014 and 2016 epidemics and genetic characterization of a new Lassa virus. Emerg. Microbes Infect. 9, 1761–1770 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
McCormick, J. B. & Fisher-Hoch, S. P. Lassa fever. Curr. Top. Microbiol. Immunol. 262, 75–109 (2002).CAS 
PubMed 

Google Scholar 
Monath, T. P., Newhouse, V. F., Kemp, G. E., Setzer, H. W. & Cacciapuoti, A. Lassa virus isolation from Mastomys natalensis rodents during an epidemic in Sierra Leone. Science 185, 263–265 (1974).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Stephenson, E. H., Larson, E. W. & Dominik, J. W. Effect of environmental factors on aerosol-induced Lassa virus infection. J. Med. Virol. 14, 295–303 (1984).CAS 
PubMed 
Article 

Google Scholar 
Wozniak, D. M. et al. Inoculation route-dependent Lassa virus dissemination and shedding dynamics in the natural reservoir – Mastomys natalensis. Emerg. Microbes Infect. 10, 2313–2325 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ter Meulen, J. et al. Hunting of peridomestic rodents and consumption of their meat as possible risk factors for rodent-to-human transmission of Lassa virus in the Republic of Guinea. Am. J. Trop. Med. Hyg. 55, 661–666 (1996).PubMed 
Article 

Google Scholar 
Downs, I. L. et al. Natural history of aerosol induced Lassa fever in non-human primates. Viruses 12, 593 (2020).CAS 
PubMed Central 
Article 

Google Scholar 
Lecompte, E. et al. Mastomys natalensis and Lassa Fever, West Africa. Emerg. Infect. Dis. 12, 1971–1974 (2006).PubMed 
PubMed Central 
Article 

Google Scholar 
Smither, A. R. & Bell-Kareem, A. R. Ecology of Lassa Virus. in 1–20 (Springer, 2021). https://doi.org/10.1007/82_2020_231.Ogbu, O., Ajuluchukwu, E. & Uneke, C. J. Lassa fever in West African sub-region: An overview. J. Vector Borne Dis. 44, 1–11 (2007).CAS 
PubMed 

Google Scholar 
Fichet-Calvet, E. et al. Fluctuation of abundance and Lassa virus prevalence in Mastomys natalensis in Guinea, West Africa. Vector Borne Zoonotic Dis. 7, 119–128 (2007).PubMed 
Article 

Google Scholar 
Fichet-Calvet, E., Becker-Ziaja, B., Koivogui, L. & Günther, S. Lassa serology in natural populations of rodents and horizontal transmission. Vector Borne Zoonotic Dis. 14, 665–674 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Lo Iacono, G. et al. Using modelling to disentangle the relative contributions of zoonotic and anthroponotic transmission: The case of Lassa fever. PLoS Negl. Trop. Dis. 9, e3398 (2015).Siddle, K. J. et al. Genomic analysis of Lassa virus during an increase in cases in Nigeria in 2018. N. Engl. J. Med. 379, 1745–1753 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kafetzopoulou, L. E. et al. Metagenomic sequencing at the epicenter of the Nigeria 2018 Lassa fever outbreak. Science 363, 74–77 (2019).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Andersen, K. G. et al. Clinical sequencing uncovers origins and evolution of Lassa virus. Cell 162, 738–750 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lalis, A. & Wirth, T. Mice and men: An evolutionary history of Lassa fever. in Biodiversity and Evolution (eds. Grandcolas, P. & Maurel, M.-C.) 189–212, https://doi.org/10.1016/B978-1-78548-277-9.50011-5 (Elsevier, 2018).Mylne, A. Q. N. et al. Mapping the zoonotic niche of Lassa fever in Africa. Trans. R. Soc. Trop. Med. Hyg. 109, 483–492 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Colangelo, P. et al. A mitochondrial phylogeographic scenario for the most widespread African rodent, Mastomys natalensis. Biol. J. Linn. Soc. 108, 901–916 (2013).Article 

Google Scholar 
Gryseels, S. et al. When viruses don’t go viral: The importance of host phylogeographic structure in the spatial spread of arenaviruses. PLoS Path 13, e1006073 (2017).Article 

Google Scholar 
Cuypers, L. N. et al. Three arenaviruses in three subspecific natal multimammate mouse taxa in Tanzania: Same host specificity, but different spatial genetic structure? Virus Evol. https://doi.org/10.1093/ve/veaa039 (2020).Vazeille, M., Gaborit, P., Mousson, L., Girod, R. & Failloux, A.-B. Competitive advantage of a dengue 4 virus when co-infecting the mosquito Aedes aegypti with a dengue 1 virus. BMC Infect. Dis. 16, 318 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Chan, K. F. et al. Investigating viral interference between influenza A virus and human respiratory syncytial virus in a ferret model of infection. J. Infect. Dis. 218, 406–417 (2018).CAS 
PubMed 
Article 

Google Scholar 
Meunier, D. Y., McCormick, J. B., Georges, A. J., Georges, M. C. & Gonzalez, J. P. Comparison of Lassa, Mobala, and Ippy virus reactions by immunofluorescence test. Lancet 1, 873–874 (1985).CAS 
PubMed 
Article 

Google Scholar 
Howard, C. R. Antigenic diversity among the Arenaviruses. in The Arenaviridae (ed. Salvato, M. S.) 37–49, https://doi.org/10.1007/978-1-4615-3028-2_3 (Springer US, 1993).Bhattacharyya, S., Gesteland, P. H., Korgenski, K., Bjørnstad, O. N. & Adler, F. R. Cross-immunity between strains explains the dynamical pattern of paramyxoviruses. Proc. Natl Acad. Sci. U. S. A. 112, 13396–13400 (2015).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Luis, A. D., Douglass, R. J., Mills, J. N. & Bjørnstad, O. N. Environmental fluctuations lead to predictability in Sin Nombre hantavirus outbreaks. Ecology 96, 1691–1701 (2015).Article 

Google Scholar 
Anderson, R. M., Jackson, H. C., May, R. M. & Smith, A. M. Population dynamics of fox rabies in Europe. Nature 289, 765–771 (1981).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Tian, H. et al. Anthropogenically driven environmental changes shift the ecological dynamics of hemorrhagic fever with renal syndrome. PLoS Pathog. 13, e1006198 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Redding, D. W., Moses, L. M., Cunningham, A. A., Wood, J. & Jones, K. E. Environmental-mechanistic modelling of the impact of global change on human zoonotic disease emergence: a case study of Lassa fever. Methods Ecol. Evol. 7, 646–655 (2016).Article 

Google Scholar 
Peterson, A. T., Moses, L. M. & Bausch, D. G. Mapping transmission risk of Lassa fever in West Africa: the importance of quality control, sampling bias, and error weighting. PLoS One 9, e100711 (2014).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Fichet-Calvet, E. & Rogers, D. J. Risk maps of Lassa fever in West Africa. PLoS. Negl. Trop. Dis. 3, e388 (2009).Basinski, A. J. et al. Bridging the gap: Using reservoir ecology and human serosurveys to estimate Lassa virus spillover in West Africa. PLoS Comput. Biol. 17, e1008811 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Iacono, G. L. et al. A unified framework for the infection dynamics of zoonotic spillover and spread. PLoS Negl. Trop. Dis. 10, e0004957 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Coumou, D. & Rahmstorf, S. A decade of weather extremes. Nat. Clim. Change 2, 491–496 (2012).ADS 
Article 

Google Scholar 
Coumou, D., Robinson, A. & Rahmstorf, S. Global increase in record-breaking monthly-mean temperatures. Clim. Change 118, 771–782 (2013).ADS 
Article 

Google Scholar 
Bathiany, S., Dakos, V., Scheffer, M. & Lenton, T. M. Climate models predict increasing temperature variability in poor countries. Sci. Adv. 4, eaar5809 (2018).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Arneth, A. Uncertain future for vegetation cover. Nature 524, 44–45 (2015).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Brandt, M. et al. Human population growth offsets climate-driven increase in woody vegetation in sub-Saharan Africa. Nat. Ecol. Evol. 1, 81 (2017).PubMed 
Article 

Google Scholar 
Herrmann, S. M., Brandt, M., Rasmussen, K. & Fensholt, R. Accelerating land cover change in West Africa over four decades as population pressure increased. Com. Earth Envir 1, 1–10 (2020).
Google Scholar 
Gibb, R., Moses, L. M., Redding, D. W. & Jones, K. E. Understanding the cryptic nature of Lassa fever in West Africa. Pathog. Glob. Health 111, 276–288 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Elith, J., Leathwick, J. R. & Hastie, T. A working guide to boosted regression trees. J. Anim. Ecol. 77, 802–813 (2008).CAS 
PubMed 
Article 

Google Scholar 
Frieler, K. et al. Assessing the impacts of 1.5 °C global warming—simulation protocol of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP2b). Geosci. Model Dev. 10, 4321–4345 (2017).ADS 
Article 

Google Scholar 
Soberón, J. & Nakamura, M. Niches and distributional areas: Concepts, methods, and assumptions. Proc. Natl Acad. Sci. U. S. A. 106, 19644–19650 (2009).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Moss, R. H. et al. The next generation of scenarios for climate change research and assessment. Nature 463, 747–756 (2010).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Lemey, P., Rambaut, A., Welch, J. J. & Suchard, M. A. Phylogeography takes a relaxed random walk in continuous space and time. Mol. Biol. Evol. 27, 1877–1885 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lukashevich, I. S. Generation of reassortants between African arenaviruses. Virology 188, 600–605 (1992).CAS 
PubMed 
Article 

Google Scholar 
Vijaykrishna, D., Mukerji, R. & Smith, G. J. D. RNA virus reassortment: an evolutionary mechanism for host jumps and immune evasion. PLoS Path 11, e1004902 (2015).Article 

Google Scholar 
Whitmer, S. L. M. et al. New lineage of Lassa Virus, Togo, 2016. Emerg. Infect. Dis. 24, 599 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Ehichioya, D. U. et al. Phylogeography of Lassa virus in Nigeria. J. Virol. 93, e00929–19 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Dellicour, S., Rose, R., Faria, N. R., Lemey, P. & Pybus, O. G. SERAPHIM: studying environmental rasters and phylogenetically informed movements. Bioinformatics 32, 3204–3206 (2016).CAS 
PubMed 
Article 

Google Scholar 
Dellicour, S. et al. Using viral gene sequences to compare and explain the heterogeneous spatial dynamics of virus epidemics. Mol. Biol. Evol. 34, 2563–2571 (2017).CAS 
PubMed 
Article 

Google Scholar 
Dellicour, S. et al. Epidemiological hypothesis testing using a phylogeographic and phylodynamic framework. Nat. Commun. 11, 5620 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Dijkstra, E. W. A note on two problems in connexion with graphs. Numer. Math. 1, 269–271 (1959).MathSciNet 
MATH 
Article 

Google Scholar 
Strahler, A. N. Quantitative analysis of watershed geomorphology. Eos, Trans. Am. Geophys. Union 38, 913–920 (1957).Article 

Google Scholar 
Kass, R. E. & Raftery, A. E. Bayes factors. J. Am. Stat. Assoc. 90, 773–795 (1995).MathSciNet 
MATH 
Article 

Google Scholar 
Ehichioya, D. U. et al. Current molecular epidemiology of Lassa virus in Nigeria. J. Clin. Microbiol. 49, 1157 (2011).PubMed 
PubMed Central 
Article 

Google Scholar 
Oloniniyi, O. K. et al. Genetic characterization of Lassa virus strains isolated from 2012 to 2016 in southeastern Nigeria. PLoS Negl. Trop. Dis. 12, e0006971 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Olesen, J. E. et al. Uncertainties in projected impacts of climate change on European agriculture and terrestrial ecosystems based on scenarios from regional climate models. Clim. Change 81, 123–143 (2007).Article 

Google Scholar 
Simo Tchetgna, H. et al. Molecular characterization of a new highly divergent Mobala related arenavirus isolated from Praomys sp. rodents. Sci. Rep. 11, 10188 (2021).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Olayemi, A. et al. New hosts of the Lassa virus. Sci. Rep. 6, 25280 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zaidi, M. B. et al. Competitive suppression of dengue virus replication occurs in chikungunya and dengue co-infected Mexican infants. Parasit. Vectors 11, 378 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Olayemi, A. et al. Widespread arenavirus occurrence and seroprevalence in small mammals, Nigeria. Parasit. Vectors 11, 416 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Nigeria Centre for Disease Control. https://ncdc.gov.ng/diseases/sitreps/?cat=5&name=An%20update%20of%20Lassa%20fever%20outbreak%20in%20Nigeria.Norris, K. et al. Biodiversity in a forestagriculture mosaic: the changing face of west Africa rainforests. Biol. Conserv. 143, 2341–2350 (2010).Article 

Google Scholar 
Stocker, T. F. et al. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the fifth assessment report of the Intergovernmental Panel on Climate Change. (Cambridge University Press, Cambridge, United Kingdom and New York, 2013).Buba, M. I. et al. Mortality among confirmed Lassa fever cases during the 2015-2016 outbreak in Nigeria. Am. J. Public Health 108, 262–264 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Tobin, E. A. et al. Knowledge of secondary school children in Edo State on Lassa fever and its implications for prevention and control. West. Afr. J. Med. 34, 101–107 (2015).CAS 
PubMed 

Google Scholar 
Saez, A. M. et al. Rodent control to fight Lassa fever: Evaluation and lessons learned from a 4-year study in Upper Guinea. PLoS Negl. Trop. Dis. 12, e0006829 (2018).Article 

Google Scholar 
Ejembi, J. et al. Contact tracing in Lassa fever outbreak response, an effective strategy for control? Online J. Public Health Inf. 11, e378 (2019).
Google Scholar 
ECHO Flash List. https://erccportal.jrc.ec.europa.eu/ECHO-Flash/ECHO-Flash-List/yy/2018/mm/2.Pigott, D. M. et al. Local, national, and regional viral haemorrhagic fever pandemic potential in Africa: a multistage analysis. Lancet 390, 2662–2672 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Kraemer, M. U. G. et al. Past and future spread of the arbovirus vectors Aedes aegypti and Aedes albopictus. Nature Microbiology https://doi.org/10.1038/s41564-019-0376-y (2019).Elith, J. et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 29, 129–151 (2006).Article 

Google Scholar 
Dhingra, M. S. et al. Global mapping of highly pathogenic avian influenza H5N1 and H5Nx clade 2.3.4.4 viruses with spatial cross-validation. eLife 5, e19571 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 190, 231–259 (2006).Article 

Google Scholar 
Elith, J. et al. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 17, 43–57 (2011).Article 

Google Scholar 
Phillips, S. J. et al. Sample selection bias and presence-only distribution models: Implications for background and pseudo-absence data. Ecol. Appl. 19, 181–197 (2009).PubMed 
Article 

Google Scholar 
Valavi, R., Elith, J., Lahoz‐Monfort, J. J. & Guillera‐Arroita, G. blockCV: An r package for generating spatially or environmentally separated folds for k-fold cross-validation of species distribution models. Methods Ecol. Evol. 10, 225–232 (2019).Article 

Google Scholar 
Suchard, M. A. et al. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 4, vey016 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Randin, C. F. et al. Are niche-based species distribution models transferable in space? J. Biogeogr. 33, 1689–1703 (2006).Article 

Google Scholar 
Lange, S. Bias correction of surface downwelling longwave and shortwave radiation for the EWEMBI dataset. Earth Syst. Dyn. 9, 627–645 (2018).ADS 
Article 

Google Scholar 
Dunne, J. P. et al. GFDL’s ESM2 global coupled climate–carbon earth system models. Part I: physical formulation and baseline simulation characteristics. J. Clim. 25, 6646–6665 (2012).ADS 
Article 

Google Scholar 
Jones, C. D. et al. The HadGEM2-ES implementation of CMIP5 centennial simulations. Geosci. Model Dev. 4, 543–570 (2011).ADS 
Article 

Google Scholar 
Dufresne, J.-L. et al. Climate change projections using the IPSL-CM5 Earth System Model: from CMIP3 to CMIP5. Clim. Dyn. 40, 2123–2165 (2013).Article 

Google Scholar 
Watanabe, M. et al. Improved climate simulation by MIROC5: Mean states, variability, and climate sensitivity. J. Clim. 23, 6312–6335 (2010).ADS 
Article 

Google Scholar 
Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An overview of CMIP5 and the experiment design. Bull. Am. Meteor. Soc. 93, 485–498 (2012).ADS 
Article 

Google Scholar 
Hurtt, G. C. et al. Harmonization of global land-use change and management for the period 850-2100 (LUH2) for CMIP6. Geosci. Model Dev. 1–65 https://doi.org/10.5194/gmd-2019-360 (2020)Jones, B. & O’Neill, B. C. Spatially explicit global population scenarios consistent with the Shared Socioeconomic Pathways. Environ. Res. Lett. 11, 084003 (2016).ADS 
Article 

Google Scholar 
Katoh, K. & Standley, D. M. MAFFT Multiple Sequence Alignment Software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Larsson, A. AliView: a fast and lightweight alignment viewer and editor for large datasets. Bioinformatics 30, 3276–3278 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ayres, D. L. et al. BEAGLE 3: Improved performance, scaling, and usability for a high-performance computing library for statistical phylogenetics. Syst. Biol. https://doi.org/10.1093/sysbio/syz020 (2019).Tavaré, S. Some probabilistic and statistical problems in the analysis of DNA sequences. Lect. Math. Life Sci. 17, 57–86 (1986).MathSciNet 
MATH 

Google Scholar 
Rambaut, A., Drummond, A. J., Xie, D., Baele, G. & Suchard, M. A. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst. Biol. 67, 901–904 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Laenen, L. et al. Spatio-temporal analysis of Nova virus, a divergent hantavirus circulating in the European mole in Belgium. Mol. Ecol. 25, 5994–6008 (2016).CAS 
PubMed 
Article 

Google Scholar 
Dellicour, S. et al. Landscape genetic analyses of Cervus elaphus and Sus scrofa: comparative study and analytical developments. Heredity 123, 228–241 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Dellicour, S. et al. Phylodynamic assessment of intervention strategies for the West African Ebola virus outbreak. Nat. Commun. 9, 2222 (2018).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Dellicour, S. et al. Phylogeographic and phylodynamic approaches to epidemiological hypothesis testing. bioRxiv https://doi.org/10.1101/788059 (2020).Dellicour, S., Rose, R. & Pybus, O. G. Explaining the geographic spread of emerging epidemics: a framework for comparing viral phylogenies and environmental landscape data. BMC Bioinform 17, 1–12 (2016).Article 

Google Scholar 
McRae, B. H. Isolation by resistance. Evolution 60, 1551–1561 (2006).PubMed 
Article 

Google Scholar 
Jacquot, M., Nomikou, K., Palmarini, M., Mertens, P. & Biek, R. Bluetongue virus spread in Europe is a consequence of climatic, landscape and vertebrate host factors as revealed by phylogeographic inference. Proc. R. Soc. Lond. B 284, 20170919 (2017).
Google Scholar 
Gill, M. S. et al. Improving Bayesian population dynamics inference: A coalescent-based model for multiple loci. Mol. Biol. Evol. 30, 713–724 (2013).CAS 
PubMed 
Article 

Google Scholar 
Karcher, M. D., Palacios, J. A., Bedford, T., Suchard, M. A. & Minin, V. N. Quantifying and mitigating the effect of preferential sampling on phylodynamic inference. PLoS Comput. Biol. 12, e1004789 (2016).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

Liu LP. Characteristics of blue algal bloom in Dianchi Lake and analysis on its cause. Res Environ Sci. 1999;12:36–37.
Google Scholar 
Liu YM, Chen W, Li DH, Shen YW, Liu YD, Song LR. Analysis of paralytic shellfish toxins in Aphanizomenon DC-1 from Lake Dianchi, China. Environ Toxicol. 2006;21:289–95.CAS 
PubMed 
Article 

Google Scholar 
Dziallas C, Grossart HP. Temperature and biotic factors influence bacterial communities associated with the cyanobacterium Microcystis sp. Environ Microbiol. 2011;13:1632–41.PubMed 
Article 

Google Scholar 
Parveen B, Ravet V, Djediat C, Mary I, Quiblier C, Debroas D, Humbert JF. Bacterial communities associated with Microcystis colonies differ from free-living communities living in the same ecosystem. Environ Microbiol Rep. 2013;5:716–24.CAS 
PubMed 

Google Scholar 
Shi LM, Cai YF, Kong FX, Yu Y. Specific association between bacteria and buoyant Microcystis colonies compared with other bulk bacterial communities in the eutrophic Lake Taihu, China. Environ Microbiol Rep. 2012;4:669–78.CAS 
PubMed 

Google Scholar 
Kouzuma A, Watanabe K. Exploring the potential of algae/bacteria interactions. Curr Opin Biotech. 2015;33:125–9.CAS 
PubMed 
Article 

Google Scholar 
Cooper MB, Smith AG. Exploring mutualistic interactions between microalgae and bacteria in the omics age. Curr Opin Plant Biol. 2015;26:147–53.PubMed 
Article 

Google Scholar 
Yang L, Xiao L. Outburst, jeopardize and control of cyanobacterial bloom in lakes. Beijing: Science Press; 2011. p. 71–212.
Google Scholar 
de-Bashan LE, Antoun H, Bashan Y. Involvement of indole-3-acetic-acid produced by the growth-promoting bacterium Azospirillum spp. in promoting growth of Chlorella vulgaris. J Phycol. 2008;44:938–47.CAS 
PubMed 
Article 

Google Scholar 
Xiao Y, Wang L, Wang X, Chen M, Chen J, Tian BY, Zhang BH. Nocardioides lacusdianchii sp. nov., an attached bacterium of Microcystis aeruginosa. Antonie van Leeuwenhoek. 2022;115:141–53.PubMed 
Article 

Google Scholar 
Shirling EB, Gottlieb D. Methods for characterization of Streptomyces species. Int J Syst Bacteriol. 1966;16:313–40.Article 

Google Scholar 
Zhang BH, Chen W, Li HQ, Zhou EM, Hu WY, Duan YQ, Mohamad OA, Gao R, Li WJ. An antialgal compound produced by Streptomyces jiujiangensis JXJ 0074T. Appl Microbiol Biotechnol. 2015;99:7673–83.CAS 
PubMed 
Article 

Google Scholar 
Zhang BH, Salam N, Cheng J, Xiao M, Li HQ, Yang JY, Zha DM, Li WJ. Citricoccus lacusdiani sp. nov., an actinobacterium promoting Microcystis growth with limited soluble phosphorus. Antonie Van Leeuwenhoek. 2016;109:1457–65.CAS 
PubMed 
Article 

Google Scholar 
Zhang BH, Salam N, Cheng J, Li HQ, Yang JY, Zha DM, Guo QG, Li WJ. Microbacterium lacusdiani sp. nov., a phosphate–solubilizing novel actinobacterium isolated from mucilaginous sheath of Microcystis. J Antibiot. 2017;70:147–51.Article 

Google Scholar 
Smibert RM, Krieg NR. Phenotypic characterization. In: Gerhardt P, Murray RGE, Wood WA, Krieg NR, editors. Methods for general and molecular bacteriology. Washington, DC: American Society for Microbiology; 1994. p. 607–54.Dong XZ, Cai MY. Manual of systematic identification of common bacteria. Beijing: Science Press; 2001. p. p349–89.
Google Scholar 
Minnikin DE, Collins MD, Goodfellow M. Fatty acid and polar lipid composition in the classification of Cellulomonas, Oerskovia and related taxa. J Appl Bacteriol. 1979;47:87–95.CAS 
Article 

Google Scholar 
Tamaoka J, Katayama-Fujimura Y, Kuraishi H. Analysis of bacterial menaquinone mixtures by high performance liquid chromatography. J Appl Bacteriol. 1983;54:31–36.CAS 
Article 

Google Scholar 
Schleifer KH, Kandler O. Peptidoglycan types of bacterial cell walls and their taxonomic implications. Bacteriol Rev. 1972;36:407–77.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tang SK, Wang Y, Chen Y, Lou K, Cao LL, Xu LH, Li WJ. Zhihengliuella alba sp. nov., and emended description of the genus Zhihengliuella. Int J Syst Evol Microbiol. 2009;59:2025–32.CAS 
PubMed 
Article 

Google Scholar 
Yoon SH, Ha SM, Kwon S, Lim J, Kim Y, Seo H, Chun J. Introducing EzBiocloud: a taxonomically united database of 16S rRNA gene sequences and whole–genome assemblies. Int J Syst Evol Microbiol. 2017;67:1613–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Saitou N, Nei M. The neighbor–joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–42.CAS 
PubMed 

Google Scholar 
Fitch WM. Toward defining the course of evolution: minimum change for a specific tree topology. Syst Zool. 1971;20:406–16.Article 

Google Scholar 
Felsenstein J. Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981;17:368–76.CAS 
PubMed 
Article 

Google Scholar 
Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39:783–91.PubMed 
Article 

Google Scholar 
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Massouras A, Hens K, Gubelmann C, Uplekar S, Decouttere F, Rougemont J, Cole ST, Deplancke B. Primer-initiated sequence synthesis to detect and assemble structural variants. Nat Methods. 2010;7:485–6.CAS 
PubMed 
Article 

Google Scholar 
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.CAS 
PubMed 
Article 

Google Scholar 
Bland C, Ramsey TL, Sabree F, Lowe M, Brown K, Kyrpides NC, Hugenholtz P. CRISPR Recognition Tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinforma. 2007;8:209.Article 

Google Scholar 
Meier-Kolthoff JP, Auch AF, Klenk HP, Göker M. Genome sequence–based species delimitation with confidence intervals and improved distance functions. BMC Bioinforma. 2013;14:60.Article 

Google Scholar 
Xiao Y, Chen J, Chen M, Deng SJ, Xiong ZQ, Tian BY, Zhang BH. Mycolicibacterium lacusdiani sp. nov., an attached bacterium of Microcystis aeruginosa. Front Microbiol. 2022;13:861291.PubMed 
PubMed Central 
Article 

Google Scholar 
Vaz-Moreira I, Lopes AR, Faria C, Spröer C, Schumann P, Nunes OC, Manaia CM. Microbacterium invictum sp. nov., isolated from homemade compost. Int J Syst Evol Microbiol. 2009;59:2036–41.PubMed 
Article 

Google Scholar 
Ohta Y, Ito T, Mori K, Nishi S, Shimane Y, Mikuni K, Hatada Y. Microbacterium saccharophilum sp. nov., isolated from a sucrose-refining factory. Int J Syst Evol Microbiol. 2013;63:2765–9.CAS 
PubMed 
Article 

Google Scholar 
Kageyama A, Takahashi Y, Ōmura S. Microbacterium deminutum sp. nov., Microbacterium pumilum sp. nov. and Microbacterium aoyamense sp. nov. Int J Syst Evol Microbiol. 2006;56:2113–7.CAS 
PubMed 
Article 

Google Scholar 
Stackebrandt E, Ebers J. Taxonomic parameters revisited: tarnished gold standards. Microbiol Today. 2006;33:152–5.
Google Scholar 
Meier-Kolthoff JP, Auch AF, Klenk HP, Göker M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinforma. 2013;14:60.Article 

Google Scholar 
Kim M, Oh HS, Park SC, Chun J. Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int J Syst Evol Microbiol. 2014;64:346–51.CAS 
PubMed 
Article 

Google Scholar 
Chun J, Oren A, Ventosa A, Christensen H, Arahal DR, da Costa MS, Rooney AP, Yi H, Xu XW, De Meyer S, Trujillo ME. Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. Int J Syst Evol Microbiol. 2018;68:461–6.CAS 
PubMed 
Article 

Google Scholar 
Hoke AK, Reynoso G, Smith MR, Gardner MI, Lockwood DJ, Gilbert NE, Wilhelm SW, Becker IR, Brennan GJ, Crider KE, Farnan SR, Mendoza V, Poole AC, Zimmerman ZP, Utz LK, Wurch LL, Steffen MM. Genomic signatures of Lake Erie bacteria suggest interaction in the Microcystis phycosphere. PLoS ONE. 2021;16:e0257017.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zhang BH, Salam N, Cheng J, Li HQ, Yang JY, Zha DM, Zhang YQ, Ai MJ, Hozzein WN, Li WJ. Modestobacter lacusdianchii sp. nov., a phosphate-solubilizing actinobacterium with ability to promote Microcystis growth. PLoS ONE. 2016;11:e0161069.PubMed 
PubMed Central 
Article 

Google Scholar  More

Case studyThe experimental setup for the field studies that provided the bacterial population and weather data used here was previously described by Belias et al. [9]. Briefly, baby spinach and lettuce plants were spray-inoculated with E. coli and S. enterica (Salmonella) onto field plots established in Davis, CA (University of California, Plant Sciences Field Research Facility); Freeville, NY (Homer C. Thompson Research Farm, Cornell University); and Murcia, Spain (La Matanza Research Farm). The spinach and lettuce varieties were selected based on their suitability for baby leaf production: lettuce var. Tamarindo, and spinach var. Acadia F1 and Seaside F1. Four replicate trials at different times of the regional growing season were carried out per location. The plants were spray-inoculated with a 104 CFU/mL cocktail of rifampin-resistant strains of commensal E. coli and attenuated S. enterica serovar Typhimurium (Salmonella), and samples were collected for bacterial cell quantification by plate counts on selective and differential media at 0, 4, 8, 24, 48, 72 and 96 h post-inoculation. Concurrent with leaf sample collection, weather variables (temperature, relative humidity (RH), solar radiation intensity, and wind velocity) were recorded hourly for the respective field locations. The hourly dew point (DP) was calculated as a function of both the hourly temperature and RH.Model for persister formation on plantsMathematical modeling to characterize the switch rate from a non-persister bacterial cell (hereafter termed “normal cell”) to a persister cell in the phyllosphere under laboratory conditions was performed as described in our previously published study [24]. Briefly, persister cell fractions were quantified in culturable EcO157 populations after inoculation onto young lettuce plants cultivated in plant growth chambers. Persister cells recovered from the lettuce phyllosphere were identified using the antibiotic lysing method [23]. The greatest persister fraction in the EcO157 population on lettuce in our laboratory investigation above was observed during population decline on leaf surfaces of plants left to dry after inoculation. Using mathematical modeling, we calculated the switch rate from an EcO157 normal to persister cell on dry lettuce plants based on these data [24]. Importantly, our laboratory conditions mimicked inoculation conditions in which E. coli arrived via water on leaves, the surfaces of which progressively dried like under prevailing weather conditions in the field.Based on the main dynamic observed in the field study data [9] and building on our previous study [24], we assumed that the total enteric pathogen population is composed of (i) non-persister (normal) cells consisting of two sub-populations, characterized by fast (n1) (CFU/100g) and slow (n2) (CFU/100g) decay, and (ii) the persister population, leading to the following model from Munther et al. [24]:$$frac{{dn_1}}{{dt}} = – theta _{n_1}n_1 – alpha _dn_1 + beta _dleft( {1 – sigma } right)hat p,$$
(1a)
$$frac{{dn_2}}{{dt}} = – theta _{n_2}n_2 – alpha _dn_2 + beta _dsigma hat p,$$
(1b)
$$frac{{dhat p}}{{dt}} = – mu _{hat p}hat p – beta _dhat p + alpha _dleft( {n_1 + n_2} right),$$
(1c)
$$n_1left( 0 right) = n_{10},n_2left( 0 right) = n_{20},, hat pleft( 0 right) = widehat {p_0},$$
(1d)
where (theta _{n_i})(1/h) is the death rate of the normal cells (subscript i = 1 for fast and i = 2 for slow), (hat p) (CFU/100 g) represents the persister cell population at time t (h), (mu _{hat p}) (1/h) reflects the persister population inactivation rate, αd (1/h) is the switch rate from normal to persister state, βd (1/h) is the switch rate from persister to the normal state, and σ ∈ (0,1) is a constant, describing the fraction of persister cells switching back to the normal, slowly decaying state. Equation (1a) and (1b) reflect the assumption that times between switching states are exponentially distributed, using the expected values (frac{1}{{alpha _d}}) (h) and (frac{1}{{beta _d}}) (h) of the respective distributions.Lacking data for potential persister populations from the field trials, we assumed the persister population is a fraction 1  > k  > 0 of the tail population, as observed in Munther et al. [24]. Regarding the model above, this implies that (hat p approx kn_2) for (t ge t^ ast), where (t^ ast approx frac{1}{{theta _{n_1}}}) (the time scale of survival for the fast-decaying population (n1)). In accord with bi-phasic decay, for (t ge t^ ast), the main dynamics for slow decaying population (n2) is dictated by (- theta _{n_2}n_2) in Eq. (1b). This suggests that the effective switch rates from n2 to (hat p) and from (hat p) back to n2 balance, so that (beta _dsigma hat p approx alpha _dn_2) in Eq. (1b). Following these ideas, we simplified the model in Eq. (1a)–(1d) to:$$frac{{dn_1}}{{dt}} = – theta _{n_1}n_1 – alpha _dn_1,$$
(2a)
$$frac{{dn_2}}{{dt}} = – theta _{n_2}n_2,$$
(2b)
$$frac{{dhat p}}{{dt}} = – theta _{hat p}hat p + alpha _dn_1,$$
(2c)
$$n_1left( 0 right) = n_{10},n_2left( 0 right) = n_{20},, hat pleft( 0 right) = widehat {p_0},$$
(2d)
where we ignored (beta _dleft( {1 – sigma } right)hat p) in (1a) since the decay rate ((theta _{n_1})) dominates. Also, by setting (theta _{hat p} = mu _{hat p} + beta _d(1 – sigma )), and using (beta _dsigma hat p approx alpha _dn_2), we obtained Eq. (2c). Furthermore, because (hat p approx kn_2) for (t ge t^ ast), (theta _{hat p} approx) (theta _{n_2}).In particular, the assumption that (hat p approx kn_2) for (t ge t^ ast) characterizes the switch rate from normal to persister cells, αd, as (alpha _d approx kalpha), where α is a hypothetical switch rate assuming that the population is composed only of fast decaying normal cells (n1) and a hypothetical persister cell population (p). In this case, the hypothetical population p starts small at (widehat {p_0}), initially increases due to switching from population n1 and then slowly decays as the n1 population is effectively inactivated (i.e., the tail of the total population is comprised entirely of p). From this perspective we utilized the following equations:$$frac{{dn_1}}{{dt}} = – theta _{n_1}n_1 – alpha n_1,$$
(3a)
$$frac{{dp}}{{dt}} = alpha n_1 – theta _pp.$$
(3b)
$$n_1left( 0 right) = n_0,, pleft( 0 right) = widehat {p_0},$$
(3c)
For mathematical justification regarding the relationship (alpha _d approx kalpha), please see the appendix (Supplementary Information).The utility of the relationship (alpha _d approx kalpha), is twofold. First, we used model fitting (Eqs. (3a)–(3c)) to determine α from the respective field study data [9]. Note that using Eqs. (3a)–(3c), we actually fit for (theta _{n_1}), θp, and α using the field study data [9]. Please reference the “model fitting procedure” section as well as the appendix for details concerning the unique determination of the aforementioned parameters, i.e., the practical identifiability of these parameters, and justification regarding the legitimacy of measured tail populations relative to the respective field trial data [9]. Second, because we wanted to examine Spearman’s correlations (corr) between αd and various weather factors, given a particular weather factor (vec w) across trials (i = 1, ldots ,n), let k be the maximum persister fraction (of the tail) across these n trials, that is, for each i, we have (alpha _{d_i} approx k_ialpha _i), so (alpha _{d_i} lesssim kalpha _i). Thus kαi represents the maximum persister switch rate for each trial i, and since corr((kvec alpha ,vec w)) =corr((vec alpha ,vec w)), we conducted the correlation analysis with the fitted α values in lieu of the actual persister switch rate αd.The assumptions behind our approach are summarized below:

A.

The tails of pathogen populations surviving on plants in the field study [9] are comprised of some fraction k ∈(0,1) of persister cells since their decay rate is quite small and they remain culturable.

B.

Because (alpha _d approx kalpha), we hereafter utilize α from model (3a)–(3c) as the representative persister switch rate.

C.

Given that the experimental context [24] for modeling persister switching occurred during population decline, we only employed trials from Belias et al. [9] that exhibited bi-phasic decay. Namely, we did not include trials in which significant bacterial growth was observed at the time scale of successive data points (the time scale in the field study is on the order of 4–16 h for the 1st day and then 24 h thereafter.)

D.

The switch rate from normal to persister cell is on average a monotonic function of some measure of environmental stress.

Based on assumptions A–D above, we applied the model (3a)–(3c) to published pathogen population size and weather data from four replicate trials in Spain, two in California, and one in NY [9]. More specifically, we fit model (3a)–(3c) to the respective population data in order to:

1.

determine values for the maximum switch rate α relative to the produce/bacteria type at the field scale,

2.

describe the correlative relationship between α and weather factors in the respective field trials.

Model fitting procedureIn model (3a)–(3c) above, we supposed dp/dtt = 0  > 0, i.e., we assumed that bacteria experience stress from the change in conditions from culture growth and inoculum suspension preparation to those on the plant surface and therefore, that persister formation increases in the phyllosphere immediately following inoculation. The report that EcO157 persister formation increases as early as 1 h after inoculation into leaf wash water [23], which could be considered as a proxy for the average oligotrophic environment that bacterial cells experience after spray inoculation onto leaves or through irrigation in the field, supports this assumption. To avoid identifiability issues between the initial persister population (widehat {p_0}) and α regarding the model fits above, we assumed that (widehat {p_0})= 1 ((widehat {p_0}) = 0 gives the same results). Thus, the initial persister population at inoculation is at its lowest, an assumption supported by Munther et al. [24], who observed an average fraction of EcO157 persisters of 0.0043% in the inoculum population. This imparts the largest possible switch rate, α, onto the population, corresponding to the largest and hence most conservative food safety risk.Let yk (CFU/100 g of produce) be the average bacteria population measurement at time tk (h) and let Pk,X (CFU/100 g of produce) represent the model prediction (total population) at time tk relative to the parameter vector (X = [ {theta _{n_1} , theta_p , alpha } ]^T). Following Eqs. (3a) and (3b), this means that ({{{{{{{mathrm{P}}}}}}}}_{k,X} = n_1left( {t_k,X} right) + p(t_k,X)). Since the population data spans multiple orders of magnitude, we calculated the residuals as (e_{k,X} = log _{10}y_k – log _{10}P_{k,X}). To determine the optimal model fit (see the appendix for details regarding a priori bounds on parameter ranges), we utilized the fminsearch function in MATLAB (MATLAB 2020b, The MathWorks, Inc., Natick, Massachusetts, United States) to determine the parameter vector X that minimizes the 2-norm of the following function F:$$| | Fleft( X right) | |_2 = left( {mathop {sum }limits_k e_{k,X}^2} right)^{frac{1}{2}}$$Correlation analysisTo provide a statistical foundation from which to relate the switch rate α and measured weather factors, we utilized Spearman and partial Spearman correlation. First, we calculated the Spearman correlation coefficients between α and each of the respective factors: 8-h average of temperature, RH, solar radiation, wind speed post-inoculation, and then we calculated the partial Spearman correlation coefficients for each respective weather factor, while controlling for the other three factors and simultaneously controlling for produce type (using lettuce =1 and spinach =0) (For details regarding why 8-h weather variables were used, see the “model fitting” subsection of the results.) The correlation coefficients were determined using the corr and partialcorr functions in MATLAB 2020b (The MathWorks, Inc., Natick, MA, USA). Considering the significant association of Salmonella α with RH and temperature, we also examined the correlation between α and dew point. Figure 1 presents a logical flow of the statistical analysis. Partial correlations with a P value of less than 0.05 were deemed significant. If the 8-h average of a weather factor exhibited a significant correlation with the switch rate, the 8-h minimum and range of the weather factor were also tested.Fig. 1: Logical flow diagram for statistical analysis.Factors in Step 1: UV (average ultraviolet radiation intensity), RH (average air relative humidity), Wind (average wind speed), and Temp (average air temperature). All weather data used in the statistical analysis were obtained over 8 h post-inoculation of E. coli and Salmonella onto lettuce and spinach leaves in the field.Full size image More

Isbell, F. et al. High plant diversity is needed to maintain ecosystem services. Nature 477, 199–202 (2011).CAS 
PubMed 
Article 

Google Scholar 
van der Sande, M. T. et al. Biodiversity in species, traits, and structure determines carbon stocks and uptake in tropical forests. Biotropica 49, 593–603 (2017).Article 

Google Scholar 
Grace, O. M. et al. Plant power: opportunities and challenges for meeting sustainable energy needs from the plant and fungal kingdoms. Plants People Planet 2, 446–462 (2020).Article 

Google Scholar 
Howes, M. J. R. et al. Molecules from nature: reconciling biodiversity conservation and global healthcare imperatives for sustainable use of medicinal plants and fungi. Plants People Planet 2, 463–481 (2020).Article 

Google Scholar 
Ulian, T. et al. Unlocking plant resources to support food security and promote sustainable agriculture. Plants People Planet 2, 421–445 (2020).Article 

Google Scholar 
Brondizio, E., Diaz, S., Settele, J. & Ngo, H. T. (eds) Global assessment report on biodiversity and ecosystem services of the Intergovernmental Science-Policy Platform on biodiversity and ecosystem services. Zenodo https://doi.org/10.5281/zenodo.3831673 (2019).Bennun, L. et al. The value of the IUCN Red List for business decision-making. Conserv. Lett. 11, e12353 (2018).Betts, J. et al. A framework for evaluating the impact of the IUCN Red List of threatened species. Conserv. Biol. 34, 632–643 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Maira, L. et al. Achieving international species conservation targets: closing the gap between top-down and bottom-up approaches. Conserv. Soc. 19, 25–33 (2021).Article 

Google Scholar 
IUCN Red List version 2022-2: Table 1a (IUCN, 2022); https://www.iucnredlist.org/resources/summary-statistics#Figure2Rivers, M. The global tree assessment—red listing the world’s trees. BGjournal 14, 16–19 (2017).
Google Scholar 
Nic Lughadha, E. et al. Extinction risk and threats to plants and fungi. Plants People Planet 2, 389–408 (2020).Article 

Google Scholar 
Silva, S. V. et al. Global estimation and mapping of the conservation status of tree species using artificial intelligence. Front. Plant Sci. 13, 839792 (2022).ThreatSearch Online Database (Botanic Gardens Conservation International, accessed 12 October 2021); https://tools.bgci.org/threat_search.phpBachman, S. P., Nic Lughadha, E. M. & Rivers, M. C. Quantifying progress toward a conservation assessment for all plants. Conserv. Biol. 32, 516–524 (2018).PubMed 
Article 

Google Scholar 
Rondinini, C., Di Marco, M., Visconti, P., Butchart, S. H. M. & Boitani, L. Update or outdate: long-term viability of the IUCN Red List. Conserv. Lett. 7, 126–130 (2014).Article 

Google Scholar 
Cazalis, V. et al. Bridging the research–implementation gap in IUCN Red List assessments. Trends Ecol. Evol. 37, 359–370 (2022).PubMed 
Article 

Google Scholar 
Dauby, G. et al. ConR: an R package to assist large-scale multispecies preliminary conservation assessments using distribution data. Ecol. Evol. 7, 11292–11303 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Stévart, T. et al. A third of the tropical African flora is potentially threatened with extinction. Sci. Adv. 5, eaax9444 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Bland, L. M., Collen, B., Orme, C. D. L. & Bielby, J. Predicting the conservation status of data-deficient species. Conserv. Biol. 29, 250–259 (2015).PubMed 
Article 

Google Scholar 
Darrah, S. E., Bland, L. M., Bachman, S. P., Clubbe, C. P. & Trias-Blasi, A. Using coarse-scale species distribution data to predict extinction risk in plants. Divers. Distrib. 23, 435–447 (2017).Article 

Google Scholar 
Pelletier, T. A., Carstens, B. C., Tank, D. C., Sullivan, J. & Espíndola, A. Predicting plant conservation priorities on a global scale. Proc. Natl Acad. Sci. USA 115, 13027–13032 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zizka, A., Silvestro, D., Vitt, P. & Knight, T. M. Automated conservation assessment of the orchid family with deep learning. Conserv. Biol. 35, 897–908 (2021).PubMed 
Article 

Google Scholar 
Walker, B. E., Leão, T. C. C., Bachman, S. P., Bolam, F. C. & Nic Lughadha, E. Caution needed when predicting species threat status for conservation prioritization on a global scale. Front. Plant Sci. 11, 520 (2020).Lughadha, E. N. et al. The use and misuse of herbarium specimens in evaluating plant extinction risks. Philos. Trans. R. Soc. B 374, 20170402 (2019).Article 

Google Scholar 
Walker, B. E., Leão, T. C. C., Bachman, S. P., Lucas, E. & Nic Lughadha, E. M. Evidence-based guidelines for developing automated assessment methods. Preprint at https://ecoevorxiv.org/zxq6s/ (2021).Isaac, N. J. B., Turvey, S. T., Collen, B., Waterman, C. & Baillie, J. E. M. Mammals on the EDGE: conservation priorities based on threat and phylogeny. PLoS ONE 2, e296 (2007).PubMed 
PubMed Central 
Article 

Google Scholar 
Grenié, M., Denelle, P., Tucker, C. M., Munoz, F. & Violle, C. funrar: an R package to characterize functional rarity. Divers. Distrib. 23, 1365–1371 (2017).Article 

Google Scholar 
Lindegren, M., Holt, B. G., MacKenzie, B. R. & Rahbek, C. A global mismatch in the protection of multiple marine biodiversity components and ecosystem services. Sci. Rep. 8, 4099 (2018).Pollock, L. J. et al. Protecting biodiversity (in all its complexity): new models and methods. Trends Ecol. Evol. 35, 1119–1128 (2020).PubMed 
Article 

Google Scholar 
Arnan, X., Cerdá, X. & Retana, J. Relationships among taxonomic, functional, and phylogenetic ant diversity across the biogeographic regions of Europe. Ecography 40, 448–457 (2017).Article 

Google Scholar 
Wong, J. S. Y. et al. Comparing patterns of taxonomic, functional and phylogenetic diversity in reef coral communities. Coral Reefs 37, 737–750 (2018).Article 

Google Scholar 
Devictor, V. et al. Spatial mismatch and congruence between taxonomic, phylogenetic and functional diversity: the need for integrative conservation strategies in a changing world. Ecol. Lett. 13, 1030–1040 (2010).PubMed 

Google Scholar 
Brum, F. T. et al. Global priorities for conservation across multiple dimensions of mammalian diversity. Proc. Natl Acad. Sci. USA 114, 7641–7646 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pollock, L. J., Thuiller, W. & Jetz, W. Large conservation gains possible for global biodiversity facets. Nature 546, 141–144 (2017).CAS 
PubMed 
Article 

Google Scholar 
Strassburg, B. B. N. et al. Global priority areas for ecosystem restoration. Nature 586, 724–729 (2020).CAS 
PubMed 
Article 

Google Scholar 
Cámara-Leret, R. et al. Fundamental species traits explain provisioning services of tropical American palms. Nat. Plants 3, 16220 (2017).Saslis-Lagoudakis, C. H. et al. Phylogenies reveal predictive power of traditional,medicinein bioprospecting. Proc. Natl Acad. Sci. USA 109, 15835–15840 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
van Kleunen, M. et al. Economic use of plants is key to their naturalization success. Nat. Commun. 11, 3201 (2020).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Molina-Venegas, R., Rodríguez, M., Pardo-de-Santayana, M., Ronquillo, C. & Mabberley, D. J. Maximum levels of global phylogenetic diversity efficiently capture plant services for humankind. Nat. Ecol. Evol. 5, 583–588 (2021).PubMed 
Article 

Google Scholar 
Molina-Venegas, R. Conserving evolutionarily distinct species is critical to safeguard human well-being. Sci. Rep. 11, 24187 (2021).Zaman, W. et al. Predicting potential medicinal plants with phylogenetic topology: inspiration from the research of traditional Chinese medicine. J. Ethnopharmacol. 281, 114515 (2021).CAS 
PubMed 
Article 

Google Scholar 
Cámara-Leret, R. et al. Climate change threatens New Guinea’s biocultural heritage. Sci. Adv. 5, eaaz1455 (2019).Lima, V. P. et al. Climate change threatens native potential agroforestry plant species in Brazil. Sci. Rep. 12, 2267 (2022).Johnson, D. V. Tropical Palms 2010 Revision Non-Wood Forest Products 10 (FAO, 2010).Johnson, D. V. & Sunderland, T. C. H. Rattan Glossary and Compendium Glossary with Emphasis on Africa Non-Wood Forest Products 16 (FAO, 2004).Ter Steege, H. et al. Hyperdominance in the Amazonian tree flora. Science 342, 1243092 (2013).PubMed 
Article 
CAS 

Google Scholar 
Zona, S. & Henderson, A. A review of animal-mediated seed dispersal of palms. Selbyana 11, 6–21 (1989).
Google Scholar 
Kissling, W. D. et al. PalmTraits 1.0, a species-level functional trait database of palms worldwide. Sci. Data 6, 178 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Tomlinson, P. B. The uniqueness of palms. Bot. J. Linn. Soc. 151, 5–14 (2006).Article 

Google Scholar 
Díaz, S. et al. The global spectrum of plant form and function. Nature 529, 167–171 (2016).PubMed 
Article 
CAS 

Google Scholar 
Muscarella, R. et al. The global abundance of tree palms. Glob. Ecol. Biogeogr. 29, 1495–1514 (2020).Article 

Google Scholar 
Dransfield, J. et al. Genera Palmarum: The Evolution and Classification of Palms (Kew Publishing, 2008).Diazgranados, M. et al. World Checklist of Useful Plant Species (Royal Botanic Gardens, Kew, 2020).Couvreur, T. L. P. & Baker, W. J. Tropical rain forest evolution: palms as a model group. BMC Biol. 11, 2–5 (2013).Article 

Google Scholar 
Faurby, S., Eiserhardt, W. L., Baker, W. J. & Svenning, J. Molecular phylogenetics and evolution: an all-evidence species-level supertree for the palms (Arecaceae). Mol. Phylogenet. Evol. 100, 57–69 (2016).PubMed 
Article 

Google Scholar 
The IUCN Red List of Threatened Species Version 2021-2 (IUCN, accessed 12 October 2021); https://www.iucnredlist.orgBaker, W. J. & Dransfield, J. Beyond genera Palmarum: progress and prospects in palm systematics. Bot. J. Linn. Soc. 182, 207–233 (2016).Article 

Google Scholar 
Henderson, A. A revision of Calamus (Arecaceae, Calamoideae, Calameae, Calaminae). Phytotaxa https://doi.org/10.11646/phytotaxa.445.1.1 (2020).Rakotoarinivo, M., Dransfield, J., Bachman, S. P., Moat, J. & Baker, W. J. Comprehensive red list assessment reveals exceptionally high extinction risk to Madagascar palms. PLoS ONE 9, e103684 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Cosiaux, A. et al. Low extinction risk for an important plant resource: conservation assessments of continental African palms (Arecaceae/Palmae). Biol. Conserv. 221, 323–333 (2018).Article 

Google Scholar 
Johnson, D. & UICN/SSC Palm Specialist Group (eds) Palms, Their Conservation and Sustained Utilization—Status Survey and Conservation Action Plan (Union Internationale pour la Conservation de la Nature et de ses Ressources, 1996).Bachman, S., Walker, B. E., Barrios, S., Copeland, A. & Moat, J. Rapid least concern: towards automating red list assessments. Biodivers. Data J. 8, e47018 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Enquist, B. J. et al. The commonness of rarity: global and future distribution of rarity across land plants. Sci. Adv. https://doi.org/10.1126/sciadv.aaz0414 (2019).Vieilledent, G. et al. Combining global tree cover loss data with historical national forest cover maps to look at six decades of deforestation and forest fragmentation in Madagascar. Biol. Conserv. 222, 189–197 (2018).Article 

Google Scholar 
Gaveau, D. L. A. et al. Rise and fall of forest loss and industrial plantations in Borneo (2000–2017). Conserv. Lett. 12, e12622 (2019).Gamoga, G., Turia, R., Abe, H., Haraguchi, M. & Iuda, O. The forest extent in 2015 and the drivers of forest change between 2000 and 2015 in Papua New Guinea: deforestation and forest degradation in Papua New Guinea. Case Stud. Environ. 5, 1442018 (2021).Cámara-Leret, R. & Bascompte, J. Language extinction triggers the loss of unique medicinal knowledge. Proc. Natl Acad. Sci. USA 118, e2103683118 (2021).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Henderson, A., Fischer, B., Scariot, A., Whitaker Pacheco, M. A. & Pardini, R. Flowering phenology of a palm community in a central Amazon forest. Brittonia 52, 149–159 (2000).Article 

Google Scholar 
Olivares, I. & Galeano, G. Leaf and inflorescence production of the wine palm (Attalea butyracea) in the dry Magdalena river valley, Colombia. Caldasia 35, 37–48 (2013).
Google Scholar 
Voeks, R. A. Disturbance pharmacopoeias: medicine and myth from the humid tropics. Ann. Assoc. Am. Geogr. 94, 868–888 (2004).
Google Scholar 
Pironon, S. et al. Potential adaptive strategies for 29 sub-Saharan crops under future climate change. Nat. Clim. Change 9, 758–763 (2019).Article 

Google Scholar 
Govaerts, R., Dransfield, J., Zona, S. & Henderson, A. World Checklist of Arecaceae (Royal Botanic Gardens, Kew, accessed 1 March 2018); http://wcsp.science.kew.org/Chamberlain, S. et al. rgbif: Interface to the Global Biodiversity Information Facility API. R package version 3.6.0 (2021).Zizka, A. et al. CoordinateCleaner: standardized cleaning of occurrence records from biological collection databases. Methods Ecol. Evol. 10, 744–751 (2019).Article 

Google Scholar 
Plants of the World Online (Royal Botanic Gardens, Kew, accessed 1 March 2018); http://www.plantsoftheworldonline.org/South, A. rworldmap v.1.3-6: Mapping global data (2016).Bivand, R. et al. maptools v.0.9-2: Tools for handling spatial objects (2017).Arel-Bundock, V., Enevoldsen, N. & Yetman, C. countrycode: an R package to convert country names and country codes. J. Open Source Softw. 3, 848 (2018).Article 

Google Scholar 
Becker, R. A., Wilks, A. R., Brownrigg, R., Minka, T. P. & Deckmyn, A. maps v.3.3.0: Draw geographical maps (2018).Pebesma, E. et al. sp v.1.2-7: Classes and methods for spatial data (2018).Wickham, H. et al. Welcome to the Tidyverse. J. Open Source Softw. 4, 1686 (2019).Article 

Google Scholar 
Wickham, H., Hester, J. & Chang, W. devtools v.1.13.5: Tools to make developing R packages easier (2018).World Geographic Scheme for Recording Plant Distributions Standard (TDWG, 2001); http://www.tdwg.org/standards/109Brummitt, R. K. World Geographical Scheme for Recording Plant Distributions (Hunt Institute for Botanical Documentation, 2001).Olson, D. M. et al. Terrestrial ecoregions of the world: a new map of life on Earth. Bioscience 51, 933–938 (2001).Article 

Google Scholar 
Moat, J. & Bachman, S. P. rCAT v.0.1.6: Conservation assessment tools (2017).Dinerstein, E. et al. An ecoregion-based approach to protecting half the terrestrial realm. Bioscience 67, 534–545 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Plants of the World Online (Royal Botanic Gardens, Kew, accessed 10 June 2020); http://www.plantsoftheworldonline.org/Csárdi, G. & FitzJohn, R. progress v.1.2.2: Terminal progress bars (2019).Microsoft Corporation & Weston, S. doParallel: Foreach parallel adaptor for the ‘parallel’ package. R package version 1.0.16 (2020).Microsoft Corporation & Weston, S. foreach: Provides foreach looping construct. R package version 1.5.0 (2020).Ooms, J., Lang, D. T. & Hilaiel, L. jsonlite v.1.7.2: A simple and robust JSON parser and generator for R (2020).Wickham, H. httr v.1.4.2: Tools for working with URLs and HTTP (2020).Global Human Footprint (Geographic), v2 (1995 – 2004) (SEDAC, accessed 14 May 2018); https://doi.org/10.7927/H4M61H5FFick, S. E. & Hijmans, R. J. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).Article 

Google Scholar 
Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850–853 (2013).CAS 
PubMed 
Article 

Google Scholar 
Wickham, H. plyr v.1.8.6: Tools for splitting, applying and combining data (2021).Wickham, H. & RStudio. tidyr v.1.1.4: Tidy messy data (2021).Wickham, H., François, R., Henry, L. & Müller, K. dplyr v.1.0.7: A grammar of data manipulation (2021).Bivand, R. et al. rgdal v.1.5-8: Bindings for the ‘geospatial’ data abstraction library (2020).Greenberg, J. A. & Mattiuzzi, M. gdalUtils v.2.0.3.2: Wrappers for the Geospatial data Abstraction Library (GDAL) utilities (2020).Hijmans, R. J. et al. raster v.3.1-5: Geographic data analysis and modeling (2020).The IUCN Red List of Threatened Species (IUCN, accessed 22 March 2018); https://www.iucnredlist.org/ThreatSearch Online Database (Botanic Gardens Conservation International, accessed 1 March 2018); https://tools.bgci.org/threat_search.phpChamberlain, S., ROpenSci & Salmon, M. rredlist: ‘IUCN’ Red List client (2020).Wickham, H. stringr v.1.4.0: Simple, consistent wrappers for common string operations (2019).Gagolewski, M. & Tartanus, B. stringi v.1.7.5: Character string processing facilities (2021).Kuhn, M. caret: Classification and regression training. R package version 6.0-86 (2020).Torgo, L. Data Mining with R, Learning with Case Studies (Chapman and Hall/CRC, 2010).Chawla, N. V., Bowyer, K. W., Hall, L. O. & Kegelmeyer, P. SMOTE: synthetic minority over-sampling technique. J. Artif. Intell. Res. 16, 321–357 (2020).Article 

Google Scholar 
Stokely, M. HistogramTools: Utility functions for R histograms. R package version 0.3.2 (2015).Sarkar, D. et al. lattice v.0.20-40: Trellis graphics for R (2020).Wickham, H. ggplot2 Elegant Graphics for Data Analysis (Springer, 2016).Auguie, B. & Antonov, A. gridExtra v.2.3: Miscellaneous functions for ‘grid’ graphics (2017).Pruim, R., Kaplan, D. T. & Horton, N. J. mosaic v.1.6.0: Project MOSAIC statistics and mathematics teaching utilities (2020).Meyer, D. & Buchta, C. proxy v.0.4-23: Distance and similarity measures (2019).Wickham, H. & Seidel, D. scales v.1.1: Scale functions for visualization (2019).Branco, P., Ribeiro, R. & Torgo, L. UBL v.0.0.6: An implementation of re-sampling approaches to utility-based learning for both classification and regression tasks (2017).Liaw, A. & Wiener, M. Classification and regression by randomForest. R News 2, 18–22 (2002).
Google Scholar 
Cohen, J. A coefficient of agreement for nominal scales. Educ. Psychol. Meas. 20, 37–46 (1960).Article 

Google Scholar 
Ripley, B. & Venables, W. nnet v.7.3-13: Feed-forward neural networks and multinomial log-linear models (2020).Warnes, G. R. et al. gdata v.2.18.0: Various R programming tools for data manipulation (2017).Wright, M. N., Wager, S. & Probst, P. ranger v.0.12.1: A fast implementation of random forests (2020).Arya, S., Mount, D., Kemp, S. E. & Jefferis, G. RANN v.2.6.1: Fast nearest neighbour search (wraps ANN Library) using L2 metric (2019).Meyer, D. et al. e1071 v.1.7-3: Misc Functions of the Department of Statistics, Probability Theory Group (formerly: E1071), TU Wien (2019).Lundberg, S. M. & Lee, S.-I. A unified approach to interpreting model predictions. Adv. Neural Inf. Process. Syst. 30, 4765–4774 (2017).
Google Scholar 
Greenwell, B. fastshap v.0.0.7: Fast approximate Shapley values (2021).Greenwell, B. vip v.0.3.2: Variable importance plots (2020).Donoghoe, M. W. glm2 v.1.2.1: Fitting generalized linear models (2018).Wickham, H. reshape2 v.1.4.4: Flexibly reshape data: a reboot of the reshape package (2020).Robin, X. et al. pROC v.1.18.0: Display and analyze ROC curves (2020).Warnes, G. R. et al. gplots v.3.0.3: Various R programming tools for plotting data (2019).Müller, K. & Bryan, J. here v.1.0.1: A simpler way to find your files (2017).Wickham, H., Hester, J., Francois, R., Jylänki, J. & Jørgensen, M. readr v.1.3.1: Read rectangular text data (2018).Wickham, H. et al. readxl v.1.3.1: Read Excel files (2019).Henry, L. & Wickham, H. purrr v.0.3.4: Functional programming tools (2020).Lin Pedersen, T. ggforce v.0.3.1: Accelerating ‘ggplot2’ (2019).Lin Pedersen, T. patchwork v.1.0.0: The composer of plots (2019).Hester, J. glue v.1.3.1: Interpreted string literals (2019).Ooms, J. & McNamara, J. writexl v.1.2: Export data frames to Excel ‘xlsx’ format (2019).Horikoshi, M. et al. ggfortify v.0.4.8: Data visualization tools for statistical analysis results (2019).Liaw, A. randomForest v.4.6-14: Breiman and Cutler’s random forests for classification and regression (2018).Kassambara, A. ggpubr v.0.2.5: ‘ggplot2’ based publication ready plots (2020).Gruca, M., Blach-Overgaard, A. & Balslev, H. African palm ethno-medicine. J. Ethnopharmacol. 165, 227–237 (2015).PubMed 
Article 

Google Scholar 
Cámara–Leret, R. & Dennehy, Z. Indigenous knowledge of New Guinea’s useful plants: a review. Econ. Bot. 73, 405–415 (2019).Article 

Google Scholar 
Macía, M. J. et al. Palm uses in Northwestern South America: a quantitative review. Bot. Rev. 77, 462–570 (2011).Article 

Google Scholar 
Orme, D. et al. caper: Comparative analyses of phylogenetics and evolution in R. R package version 1.0.1 https://cran.r-project.org/package=caper (2018).Kowarik, A. & Templ, M. Imputation with the R package VIM. J. Stat. Softw. 74, 1–16 (2016).Alfons, A. & Templ, M. Estimation of social exclusion indicators from complex surveys: the R package laeken. J. Stat. Softw. 54, 1–25 (2013).Article 

Google Scholar 
Milliken, W., Walker, B. E., Howes, M. J. R., Forest, F. & Nic Lughadha, E. Plants used traditionally as antimalarials in Latin America: mining the tree of life for potential new medicines. J. Ethnopharmacol. 279, 114221 (2021).PubMed 
Article 

Google Scholar 
Fritz, S. A. & Purvis, A. Selectivity in mammalian extinction risk and threat types: a new measure of phylogenetic signal strength in binary traits. Conserv. Biol. 24, 1042–1051 (2010).PubMed 
Article 

Google Scholar 
Suchard, M. A. et al. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 4, vey016 (2018).Paradis, E. & Schliep, K. Ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528 (2019).CAS 
PubMed 
Article 

Google Scholar 
Govaerts, R., Nic Lughadha, E., Black, N., Turner, R. & Paton, A. The World Checklist of Vascular Plants, a continuously updated resource for exploring global plant diversity. Sci. Data 8, 215 (2021).Yu, G. ggplotify v.0.0.4: Convert plot to ‘grob’ or ‘ggplot’ object (2019).Yu, G. aplot v.0.0.3: Decorate a ‘ggplot’ with associated information (2020).Slowikowski, K. et al. ggrepel v.0.8.1: Automatically position non-overlapping text labels with ‘ggplot2’ (2019).Schloerke, B. et al. GGally v.1.4.0: Extension to ‘ggplot2’ (2018).Rubis, B. et al. hrbrthemes v.0.6.0: Additional themes, theme components and utilities for ‘ggplot2’ (2019).Henry, L., Wickham, H. & Chang, W. ggstance v.0.3.3: Horizontal ‘ggplot2’ components (2019).Yu, G., Smith, D. K., Zhu, H., Guan, Y. & Lam, T. T. Y. Ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol. 8, 28–36 (2017).Article 

Google Scholar 
Brown, C. hash v.2.2.6.1: Full feature implementation of hash/associated arrays/dictionaries (2019).Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag, 2016).R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2020).RStudio Team. RStudio: Integrated Development for R (RStudio, 2021).Bellot, S. et al. Workflow and code used to perform palm extinction risk and regional palm use resilience analyses. Zenodo https://doi.org/10.5281/zenodo.6678122 (2022). More

Bellot, S. et al. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-022-01858-0 (2022).Article 

Google Scholar 
Eiserhardt, W. L. et al. Ann. Bot. 108, 1391–1416 (2011).Article 

Google Scholar 
Muscarella, R. et al. Glob. Ecol. Biogeogr. 29, 1495–1514 (2020).Article 

Google Scholar 
Cámara-Leret, R. et al. Nat. Plants 3, 16220 (2017).Article 

Google Scholar 
The IUCN Red List of Threatened Species (IUCN, 2018).BGCI ThreatSearch Online Database (BGCI, 2018).Carlos-Júnior, L. A. et al. Divers. Distrib. 25, 743–757 (2019).Article 

Google Scholar 
Pollock, L. J., Thuiller, W. & Jetz, W. Nature 546, 141–144 (2017).CAS 
Article 

Google Scholar 
Rice, J. et al. The IPBES Regional Assessment Report on Biodiversity and Ecosystem Services for the Americas (IPBES, 2018).Coelho de Souza, F. et al. Nat. Ecol. Evol. 3, 1754–1761 (2019).Article 

Google Scholar  More

Literature values of R
To and Q10 of leaf respirationData of RTo were read from texts, tables, and figures in all available literature (18 species; Supplementary Tables 1, 2) when measured more than once within a period of darkness in lab- and field studies where measurement temperature, To, was kept constant. The RTo-initial was defined as the initial measurement of RTo for each study/species, and further values of RTo at later points within the same night of the same study were read as well.Apparent- and inherent temperature sensitivities (Q10, Equation 1; Fig. 2b) were obtained from all available literature (ten species; Supplementary Table 2) where in the same study/species, both nocturnal values of Q10,app and of Q10,inh were obtained in response to long-term natural T-changes in the environment during the night (hours) and nocturnal values were obtained in response to short-term artificial T-changes (max 30 min), respectively.Measurements of R
To and Q10 of leaf respirationIn the field (United Kingdom, Denmark, Panama, Colombia and Brazil), RTo (µmol CO2 m−2 s−1) in 16 species (Supplementary Tables 1, 3) was measured through nocturnal periods at constant To (controlled either by block-T or leaf-T) with infra-red gas analysers (Li-Cor-6400(XT) or Li-Cor-6800, Lincoln, Nebraska, USA). Mature, attached leaves positioned in the sunlight throughout the day were chosen. Target [CO2] in the leaf cuvette was set to ambient, ranging from 390 to 410 ppm, depending on when measurements were made, and target RH = 65 ± 10%, with a flow rate of 300 µmol s–1. The RTo-initial was defined as RTo at first measurement after darkness 30 min after sunset (to conservatively avoid light-enhanced dark respiration, LEDR50,51. Leak tests were conducted prior to measurements52. The temporal resolution of measurements varied between every three minutes to once per hour for the different species. Data were subsequently binned in hourly bins.Measurements to derive Q10,inh and Q10,app were conducted in two species in a T-controlled growth cabinet and in six species in the field (Supplementary Table 2), where Q10,inh was measured in response to 10–30 min of artificial changes in T and Q10,app was calculated from measurements of RT in response to T of the environment (growth cabinet or field) at the beginning of the night and again at the end of the night (hours apart).Tree level measurements in whole-tree chambersThe night-time respiratory efflux of the entire above-ground portion (crown and bole) in large growing trees of Eucalyptus tereticornis was measured in whole-tree chambers (WTCs) in Richmond, New South Wales (Australia, (33°36ʹ40ʺS, 150°44ʹ26.5ʺE). The WTCs are large cylindrical structures topped with a cone that enclose a single tree rooted in soil (3.25 m in diameter, 9 m in height, volume of ~53 m3) and under natural sunlight, air temperature and humidity conditions. An automated system measured the net exchange of CO2 between the canopy and the atmosphere within each chamber at 15-min resolution. During the night, we used the direct measurements of CO2 evolution (measured with an infra-red gas analyser; Licor 7000, Li-Cor, Inc., Lincoln, NE)53,54 as a measure of respiration.Due to the high noise-to-signal ratio in the CO2-exchange measurements from this system when analysing the high-resolution temporal variation through each night, we chose to only analyse temporal variation in tree-RT for the nights when tree-RT-initial were amongst the top 10% of CO2-exchange signals for the entire data set. The resulting data spanned 62 nights and included hourly average measurements from three replicate chambers.Data analysis of R
To
Measurements of nocturnal leaf respiration under constant temperature conditions (RTo) were divided by the initial rate of respiration (RTo-initial) at the onset of each night. Hourly means of RTo/RTo-initial were calculated for each leaf replicate to remove measurement noise and reduce bias due to the measurement of some species at more frequent intervals throughout the night. For species with multiple leaf replicates, these hourly means of RTo/RTo-initial were then combined to create hourly averages of RTo/RTo-initial at the species level. For each species, these values were plotted as a function of time to demonstrate how RTo/RTo-initial decreases with time since the onset of darkness, from sunset until sunrise (Supplementary Fig. 1). For each species, hourly means of RTo/RTo-initial plotted as a function of time were linearised by log-transforming data and the slope of the relationship determined. To test whether the slopes of the lines differed significantly within plant functional groups (woody, non-woody), species originating from the same biome (temperate, tropical) or species measured under the same conditions (lab, field), the slopes of the lines for all species from a given functional group, biome or measurement condition were tested pairwise against each other using the slope, standard error and sample size (number of points on the x-axis) for each line and applying a 0.05 cut-off for p values after Bonferroni correction for multiple testing. 11 out of 701 comparisons came out as being significantly different, which is why within-group slope differences were considered to be overall non-significant for this analysis. t-tests were used to test whether the slopes differed between plant functional groups (tree, non-woody), species originating from different biomes (temperate, tropical) and species measured under different environmental conditions (lab, field). In these tests, the degrees of freedom varied according to the different sample sizes. Since RTo/RTo-initial plotted as a function of time always starts at 1, the intercepts do not differ between species. t-tests were performed on linearised power functions by log-transforming data in order to test potential differences between lab and field, origin of species, between woody and non-woody species and between temperate and tropical biomes. Since these functions were statistically indistinguishable in each pairing, all measurements of nocturnal leaf respiration under constant temperature conditions (n = 967 nights, 31 species) were collated into a single plot. The data were binned hourly since some studies had very few measurements on half-hourly steps. A power function was fitted with a weighting of each hourly binned value using 1/(standard error of the mean). The power function was chosen as it, better than the exponential- or linear function, can capture both sudden steep- as well as slower decrease in RTo/RTo-initial in different species. The 95% confidence interval of the power function, following the new model equation, overlaps with all the 95% confidence intervals of the hourly binned values (Fig. 1a). All data analysis, including statistical analysis and figures were performed using Python version 3.9.4.Evaluation of new equationWe performed four sets of simulations (S1-S4) using different representations of leaf and plant respiration as outlined in Supplementary Table 4. Evaluation of Equation 4 (S2; Equation 3 from Fig. 1a merged with Equation 1) in comparison with Equation 1 (S1) and Equation 5 (S4) in comparison with Equation 2 (S3), respectively, for predictions of nocturnal variation in response to natural variation in temperature, was conducted by use of independent sets of leaf level data and tree scale data. The effect of including variable nocturnal RTo is estimated as the difference between S1 and S2 and between S3 and S4, respectively.The first data set used for the evaluation consists of nine broad-leaf species for which spot measurements of leaf respiration under ambient conditions were taken at sunset and before sunrise in the field (Fig. 1b and Supplementary Fig. 2a). Of these nine species, three species (Fig. 1c) were further measured throughout the night at ambient conditions. Further, whole-tree measurements measured throughout the night at ambient conditions (Supplementary Fig. 3a–d) were also used for evaluation. Finally, comparisons of Q10,inh with Q10,app in another ten species were used to test if RTo appeared constant as assumed in Equation 1 (Supplementary Tables 2, 3 and Fig. 2b).To validate the suitability of Equation 4 and Equation 5 over equations with full temporal control, modelled respiration values were compared against observed measurements for three species at the leaf level (Supplementary Fig. 2b–d) and for Eucalyptus tereticornis at the whole-tree level using three chamber replicates and during 62 nights using hourly measurements (Supplementary Fig. 3a–d). Linear fits were applied, using ordinary least squares regressions, to plots of normalised respiration (({R}_{T}/{R}_{{T}_{0}})) predicted by the four models against the observed values. The first measurements of the night were excluded from the fits, as these were necessarily equal to unity. The standardised residuals (S) in Supplementary Figs. 2c, 3b are calculated using the equation ({S}_{i}=({R}_{{{{{{{rm{modelled}}}}}}}_{i}}/{R}_{{{{{{{rm{Modelled}}}}}}}_{0}}-{R}_{{T}_{i}}/{R}_{{T}_{0}})/sqrt{(mathop{sum }nolimits_{i}^{N}{({R}_{{{{{{{rm{modelled}}}}}}}_{i}}/{R}_{{{{{{{rm{Modelled}}}}}}}_{0}}-{R}_{{T}_{i}}/{R}_{{T}_{0}})}^{2})/{df}}), for the residual of the ith measurement, where the sum is over all measurements, df is the number of degrees of freedom, and Rmodelled are the respiration values modelled by the four equations in Supplementary Table 4.Evaluation is done by comparing observed and simulated RT/RT, initial. We evaluate the nocturnal evolution of RT/RT, initial and use (i) one-to-one line figures that include fitted regression line, R2, p value and RMSE, (ii) Taylor diagrams and (iii) use plots of standardised residuals against temperature and hours since darkness for a qualitative assessment of the simulations, to identify whether there are any model biases at specific times or temperatures. Model evaluation, statistical analysis and figures were done using python version 3.9.4.Global scale modelling of plant R and NPP
We applied the novel formulation derived in this study (Equation 4 and Equation 5) to quantify the impact of incorporating variable RTo on simulated plant R and NPP globally using the JULES land surface model32,33 following simulations outlined in Supplementary Table 4.Plant respiration in JULES and simulations for this study: The original leaf respiration representation in JULES follows either eqn 1 ({{R}_{T}={R}_{{T}_{0}}{Q}}_{10}^{(T-{T}_{0})/10}) with Q10 = 2 and To = 25 oC or Equation 1 with an additional denominator ({{R}_{T}={R}_{{T}_{0}}{Q}}_{10}^{(T-{T}_{0})/10}/leftlfloor left(1+{e}^{0.3(T-{T}_{{upp}})}right)times left(1+{e}^{0.3({T}_{{low}}-T)}right)rightrfloor) (Equation 6). For the purpose of this application, we have used Equation 1 to represent leaf respiration in standard JULES simulations. The remaining components of maintenance respiration in JULES, i.e. fine root and wood are represented as a function of leaf to root and leaf to wood nitrogen ratios and leaf respiration rates following RT (β + (Nr + Ns)/Nl) (Equation 6) with RT as leaf respiration, Nr, Ns and Nl as root, stem and leaf Nitrogen content respectively and β as a soil water factor (Equation 42 in ref. 32). This implies that any variation in leaf respiration is passed to root and wood respiration as well30,31,35. Growth respiration is estimated as a fraction (25%) of the difference between GPP and maintenance respiration (Rm) expressed as Rg = 0.25 (GPP-Rm).JULES version 5.2 was modified to simulate leaf and plant respiration using the various descriptions (Equations 1–5) outlined in the modelling protocol in Supplementary Table 4. JULES uses standard astronomical equations to calculate the times of sunrise and sunset on a given day at each grid point. We used the model leaf temperature and RT at the timestep at or immediately preceding sunset to represent Tsunset, and RT,sunset and at every timestep through the night, the time since sunset (h) was updated. We performed global simulations for the period 2000–2018 with JULES, using the global physical configuration GL8, which is an update from GL755. We used WFDEI meteorological forcing data56 available at 0.5-degree spatial resolution and 3-h temporal resolution, and disaggregated and run in JULES with a 15 min timestep. Simulations were performed using nine plant functional types (PFTs)33. To isolate the effects of the new formulation on simulated Rp and NPP from possible impacts on leaf area index (LAI) or vegetation dynamics, we prescribed vegetation phenology via seasonal LAI fields and vegetation fractional cover based on the European Space Agency’s Land Cover Climate Change Initiative (ESA LC_CCI) global vegetation distribution57, processed to the JULES nine PFTs and re-gridded to the WFDEI grid. Annual variable fields of CO2 concentrations are based on annual mean observations from Mauna Loa58. JULES was spun up using the three cycles of the 2000–2018 meteorological forcing data to equilibrate the soil moisture stores. The mean annual output of Rp and NPP over the study period (2000–2018) is computed for all simulations and the effect of the new formulation is presented as the difference between the temporal mean of simulations with and without nocturnal variation in whole plant RTo for NPP and vice versa for Rp and as percentage respect to simulations without nocturnal variation in RTo. Results are presented for grid cells where grid level NPP is >50 g m−2 yr −1 in the standard simulations to avoid excessively large % effects at very low NPP. Output from JULES was analysed and plotted using python version 2.7.16.PermitsNo permit was required in Denmark as measurements were taken in private land (of author) and public land and measurements were non-destructive. Data were collected under the Panama Department of the Environment (current name MiAmbiente) research permit under the name of Dr Kaoru Kitajima. Permit number: SE/P-16-12. Data in Brazil were collected under the minister of Environment (Ministério do Meio Ambiente—MMA), Instituto Chico Mendes de Conservação da Biodiversidade—ICMBio, Sistema de Autorização e Informação em Biodiversidade—SISBIO permit number 47080-3. No permit was required in Colombia as measurements were taken on private land, no plant samples were collected, and trees were part of an existing experiment for which one of the co-authors is the lead. No access permits were required in the UK as they were conducted on the campus of own university plus in their own private garden.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

White, T. D. et al. Pleistocene Homo sapiens from Middle Awash, Ethiopia. Nature 423, 742–747 (2003).Article 

Google Scholar 
McDougall, I., Brown, F. H. & Fleagle, J. G. Stratigraphic placement and age of modern humans from Kibish, Ethiopia. Nature 433, 733–736 (2005).Article 

Google Scholar 
Potts et al. Increased ecological resource variability during a critical transition in hominin evolution. Sci. Adv. 6, eabc8975 (2020).Article 

Google Scholar 
Faith, J. T. et al. Rethinking the ecological drivers of hominin evolution. Trends Ecol. Evol. 36, 797–807 (2021).Article 

Google Scholar 
Mounier, A. & Lahr, M. M. Deciphering African late Middle Pleistocene hominin diversity and the origin of our species. Nat. Commun. 10, 10–13 (2019).Article 

Google Scholar 
Cohen, A. et al. The Hominin Sites and Paleolakes Drilling Project: inferring the environmental context of human evolution from eastern African rift lake deposits. Sci. Drill. 21, 1–16 (2016).Article 

Google Scholar 
Foerster, V. et al. Climatic change recorded in the sediments of the Chew Bahir basin, southern Ethiopia, during the last 45,000 years. Quat. Int. 274, 25–37 (2012).Article 

Google Scholar 
Fischer, M. L. et al. Determining the pace and magnitude of lake level changes in southern Ethiopia over the last 20,000 years using lake balance modeling and SEBAL. Front. Earth Sci. https://doi.org/10.3389/feart.2020.00197 (2020).Roberts, H. M. et al. Using multiple chronometers to establish a long, directly-dated lacustrine record: constraining >600,000 years of environmental change at Chew Bahir, Ethiopia. Quat. Sci. Rev. 266, 107025 (2021).Article 

Google Scholar 
Galway-Witham, J., Cole, J. & Stringer, C. Aspects of human physical and behavioural evolution during the last 1 million years. J. Quat. Sci. https://doi.org/10.1002/jqs.3137 (2019).Bergström, A., Stringer, C., Hajdinjak, M., Scerri, E. M. L. & Skoglund, P. Origins of modern human ancestry. Nature 590, 229–237 (2021).Article 

Google Scholar 
Scerri, E. M. L. et al. Did our species evolve in subdivided populations across Africa, and why does it matter? Trends Ecol. Evol. 33, 582–594 (2018).Article 

Google Scholar 
Schaebitz, F. et al. Hydroclimate changes in eastern Africa over the past 200,000 years may have influenced early human dispersal. Commun. Earth Environ. 2, 1–10 (2021).Article 

Google Scholar 
Lahr, M. M. & Foley, R. A. Towards a theory of modern human origins: geography, demography, and diversity in recent human evolution. Yearb. Phys. Anthropol. 41, 137–176 (1998).Article 

Google Scholar 
Duesing, W. et al. Multiband wavelet age modeling for a 293 m (600 kyr) sediment core from Chew Bahir basin, southern Ethiopian Rift. Front. Earth Sci. 9, 1–15 (2021).Article 

Google Scholar 
Foerster, V. et al. Towards an understanding of climate proxy formation in the Chew Bahir basin, Southern Ethiopian Rift. Palaeogeogr. Palaeoclimatol. Palaeoecol. 501, 111–123 (2018).Article 

Google Scholar 
Kaboth-Bahr, S. et al. Paleo-ENSO influence on African environments and early modern humans. Proc. Natl Acad. Sci. USA 118, e2018277118 (2021).Article 

Google Scholar 
Duesing, W. et al. Changes in the cyclicity and variability of the eastern African paleoclimate over the last 620 kyrs. Quat. Sci. Rev. 273, 107219 (2021).Article 

Google Scholar 
Herbert, T., Cleaveland Peterson, L., Lawrence, K. T. & Liu, Z. Tropical ocean temperatures over the past 3.5 million years. Science 328, 1530–1534 (2010).Article 

Google Scholar 
Trauth, M. H. et al. High- and low-latitude forcing of Plio–Pleistocene East African climate and human evolution. J. Hum. Evol. 53, 475–486 (2007).Article 

Google Scholar 
Bereiter, B. et al. Revision of the EPICA Dome C CO2 record from 800 to 600 kyr before present. Geophys. Res. Lett. 42, 542–549 (2015).Article 

Google Scholar 
Clemens, S., Prell, W., Murray, D., Shimmield, G. & Weedon, G. Forcing mechanisms of the Indian Ocean monsoon. Nature 353, 720–725 (1991).Article 

Google Scholar 
Owen, R. B. et al. Progressive aridification in East Africa over the last half million years and implications for human evolution. Proc. Natl Acad. Sci. USA 115, 11174–11179 (2018).Article 

Google Scholar 
Grant, K. M. et al. A 3 million year index for North African humidity/aridity and the implication of potential pan-African humid periods. Quat. Sci. Rev. 171, 100–118 (2017).Article 

Google Scholar 
Trauth, M. H. et al. Recurring types of variability and transitions in the ~620 kyr record of climate change from the Chew Bahir basin, southern Ethiopia. Quat. Sci. Rev. 266, 106777 (2021).Article 

Google Scholar 
Wagner, B. et al. Mediterranean winter rainfall in phase with African monsoons during the past 1.36 million years. Nature 573, 256–260 (2019).Article 

Google Scholar 
Gunz, P. et al. Early modern human diversity suggests subdivided population structure and a complex out-of-Africa scenario. Proc. Natl Acad. Sci. USA 106, 6094–6098 (2009).Article 

Google Scholar 
Pearson, O. M., Royer, D. F., Grine, F. E. & Fleagle, J. G. A description of the Omo I postcranial skeleton, including newly discovered fossils. J. Hum. Evol. 55, 421–437 (2008).Article 

Google Scholar 
Sahle, Y., Morgan, L. E., Braun, D. R., Atnafu, B. & Hutchings, W. K. Chronological and behavioural contexts of the earliest Middle Stone Age in the Gademotta Formation, Main Ethiopian Rift. Quat. Int. 331, 6–19 (2014).Article 

Google Scholar 
McBrearty, S. & Tryon, C. in Transitions Before the Transition: Evolution and Stability in the Middle Paleolithic and Middle Stone Age (eds Hovers, E. & Kuhn, S. L.) Ch. 14 (Springer, 2006).Clark, J. D., Beyene, Y. & WoldeGabriel, G. Stratigraphic, chronological and behavioural contexts of Pleistocene Homo sapiens from Middle Awash, Ethiopia. Nature 423, 747–752 (2003).Article 

Google Scholar 
de la Torre, I., Mora, R., Arroyo, A. & Benito-Calvo, A. Acheulean technological behaviour in the Middle Pleistocene landscape of Mieso (East-Central Ethiopia). J. Hum. Evol. 76, 1–25 (2014).Article 

Google Scholar 
Vidal, C. et al. Age of the oldest known Homo sapiens from eastern Africa. Nature 601, 579–583 (2022).Article 

Google Scholar 
Sahle, Y. et al. in Modern Human Origins and Dispersal (eds. Sahle, Y. et al.) 73–104 (Kerns Verlag, 2019).Hewitt, G. M. Genetic consequences of climatic oscillations in the Quaternary. Phil. Trans. R. Soc. B 359, 183–195 (2004).Article 

Google Scholar 
Brandt, S., Hildebrand, E., Vogelsang, R., Wolfhagen, J. & Wong, H. A new MIS 3 radiocarbon chronology for Mochena Borago Rockshelter, SW Ethiopia: implications for the interpretation of Late Pleistocene chronostratigraphy and human behavior. J. Archaeol. Sci. Rep. 11, 352–369 (2017).
Google Scholar 
Creanza, N., Kolodny, O. & Feldman, M. W. Greater than the sum of its parts? Modelling population contact and interaction of cultural repertoires. J. R. Soc. Interface 14, 20170171 (2017).Article 

Google Scholar 
Derex, M., Perreault, C. & Boyd, R. Divide and conquer: intermediate levels of population fragmentation maximize cultural accumulation. Phil. Trans. R. Soc. B 373, 20170062 (2018).Article 

Google Scholar 
Derex, M. & Mesoudi, A. Cumulative cultural evolution within evolving population structures. Trends Cogn. Sci. 24, 654–667 (2020).Article 

Google Scholar 
Grove, M. Environmental complexity, life history, and encephalisation in human evolution. Biol. Phil. 32, 395–420 (2017).Article 

Google Scholar 
Haidle, M. N. et al. The Nature of Culture: an eight-grade model for the evolution and expansion of cultural capacities in hominins and other animals. J. Anthropol. Sci. 93, 43–70 (2015).
Google Scholar 
Brooks, A. S. et al. Long-distance stone transport and pigment use in the earliest Middle Stone Age. Science 360, 90–94 (2018).Article 

Google Scholar 
Murren, C. J., Julliard, R., Schichting, C. D., & Clobert, J. in Dispersal (eds Clobert, J. et al.) 261–273 (Oxford Univ. Press, 2001).Grove, M. et al. Climatic variability, plasticity, and dispersal: a case study from Lake Tana, Ethiopia. J. Hum. Evol. 87, 32–47 (2015).Article 

Google Scholar 
Hershkovitz, I. et al. The earliest modern humans outside Africa. Science 359, 456–459 (2018).Article 

Google Scholar 
Tryon, C. A. The Middle/Later Stone Age transition and cultural dynamics of Late Pleistocene East Africa. Evol. Anthropol. 28, 267–282 (2019).Article 

Google Scholar 
Grove, M. & Blinkhorn, J. Neural networks differentiate between Middle and Later Stone Age lithic assemblages in eastern Africa. PLoS ONE 15, e0237528–27 (2020).Article 

Google Scholar 
Grove, M. & Blinkhorn, J. Testing the integrity of the Middle and Later Stone Age cultural taxonomic division in eastern Africa. J. Paleolit. Archaeol. 4, 14 (2021).Article 

Google Scholar 
Assefa, Z., Lam, Y. M. & Mienis, H. K. Symbolic use of terrestrial gastropod opercula during the Middle Stone Age at Porc-Epic Cave, Ethiopia. Curr. Anthropol. 49, 746–756 (2008).Article 

Google Scholar 
Assefa, Z. et al. Engraved ostrich eggshell from the Middle Stone Age contexts of Goda Buticha, Ethiopia. J. Archaeol. Sci. Rep. 17, 723–729 (2018).
Google Scholar 
Laskar, J. et al. A long term numerical solution for the insolation quantities of Earth. Astron. Astrophys. 428, 261–285 (2004).Article 

Google Scholar 
Campisano, C. et al. The Hominin Sites and Paleolakes Drilling Project: high-resolution paleoclimate records from the East African Rift System and their implications for understanding the environmental context of hominin evolution. PaleoAnthropology https://doi.org/10.1130/abs/2017am-295426 (2017).Davidson, A. The Omo River Project: Reconnaissance Geology and Geochemistry of Parts of Ilubabor, Kefa, Gemu Gofa and Sidamo (Ethiopian Institute of Geological Surverys, 1983).Noren, A. HSPDP-CHB_public. OSF https://doi.org/10.17605/OSF.IO/M8QU5 (2022).Gebregiorgis, D. et al. Modern sedimentation and authigenic mineral formation in the Chew Bahir basin, southern Ethiopia: implications for interpretation of late Quaternary paleoclimate records. Front. Earth Sci. 9, 244 (2021).Article 

Google Scholar 
Folk, R. L. & Ward, W. C. Brazos River bar [Texas]: a study in the significance of grain size parameters. J. Sediment. Res. 27, 3–26 (1957).Article 

Google Scholar 
Blott, S. J. & Pye, K. GRADISTAT: a grain size distribution and statistics package for the analysis of unconsolidated sediments. Earth Surf. Process. Landf. 26, 1237–1248 (2001).Article 

Google Scholar 
Opitz, S. et al. Holocene lake stages and thermokarst dynamics in a discontinuous permafrost affected region, north-eastern Tibetan Plateau. J. Asian Earth Sci. 76, 85–94 (2013).Article 

Google Scholar 
Opitz, S. et al. Spatio-temporal pattern of detrital clay-mineral supply to a lake system on the north-eastern Tibetan Plateau, and its relationship to late Quaternary paleoenvironmental changes. Catena 137, 203–218 (2016).Article 

Google Scholar 
Viehberg, F. A. et al. Environmental change during MIS 4 and MIS 3 opened corridors in the Horn of Africa for Homo sapiens expansion. Quat. Sci. Rev. 202, 139–153 (2018).Article 

Google Scholar 
Koerting, F. et al. Drill core mineral analysis by means of the hyper spectral imaging spectrometer HySpex, XRD and ASD in proximity of the Mytina Maar, Czech Republic. Int. Arch. Photogram. Remote Sens. Spatial Inf. XL-1/W5, 417–424 (2015).Giosan, L., Flood, R. D. & Aller, R. C. Paleoceanographic significance of sediment color on western North Atlantic drifts: I. Origin of color. Mar. Geol. 189, 25–41 (2002).Article 

Google Scholar 
Weltje, G. J. et al. in Micro-XRF Studies of Sediment Cores (eds Croudace, I. W. & Rothwell, R. G.) 507–534 (Springer, 2015).Löwemark, L. et al. Normalizing XRF-scanner data: a cautionary note on the interpretation of high resolution records from organic-rich lakes. J. Asian Earth Sci. 40, 1250–1256 (2011).Article 

Google Scholar 
Lyle, M. et al. Data report: raw and normalized elemental data along the Site U1338 splice from X-ray fluorescence scanning. Proc. Integr. Ocean Drill. Program 320/321, 1–19 (2012).
Google Scholar 
Schlolaut, G. et al. An extended and revised Lake Suigetsu varve chronology from ∼50 to ∼10 ka bp based on detailed sediment micro-facies analyses. Quat. Sci. Rev. 200, 351–366 (2018).Article 

Google Scholar 
Surdam, R. C. & Eugster, H. P. Mineral reactions in the sedimentary deposits of the Lake Magadi region, Kenya. Geol. Soc. Am. Bull. 87, 1739–1752 (1976).Article 

Google Scholar 
Davies, S. J., Lamb, H. F. & Roberts, S. J. in Micro-XRF Studies of Sediment Cores (eds Croudace, I. W. & Rothwell, R. G.) 189–226 (Springer, 2015).Elbert, E. et al. Late Holocene air temperature variability reconstructed from sediments of Laguna Escondida, Patagonia, Chile (45° 30’ S). Palaeogeogr. Palaeoclimatol. Palaeoecol. 369, 482–492 (2013).Article 

Google Scholar 
Trauth, M. H., Larrasoaña, J. C. & Mudelsee, M. Trends, rhythms and events in Plio–Pleistocene African climate. Quat. Sci. Rev. 28, 399–411 (2009).Article 

Google Scholar 
Trauth, M. H. et al. Abrupt or gradual? Change point analysis of the Late Pleistocene–Holocene climate record from Chew Bahir, southern Ethiopia. Quat. Res. 90, 321–330 (2018).Article 

Google Scholar  More