Philippa Kaur

Chase, J. M. et al. The interaction between predation and competition: A review and synthesis. Ecol. Lett. 5, 302–315. https://doi.org/10.1046/j.1461-0248.2002.00315.x (2002).Article 

Google Scholar 
Droge, E., Creel, S., Becker, M. S. & M’Soka, J. Risky times and risky places interact to affect prey behaviour. Nat. Ecol. Evol. 1, 1123–1128. https://doi.org/10.1038/s41559-017-0220-9 (2017).Article 

Google Scholar 
Allesina, S. & Levine Jonathan, M. A competitive network theory of species diversity. Proc. Natl. Acad. Sci. U.S.A. 108, 5638–5642. https://doi.org/10.1073/pnas.1014428108 (2011).Article 
ADS 

Google Scholar 
Bairey, E., Kelsic, E. D. & Kishony, R. High-order species interactions shape ecosystem diversity. Nat. Commun. 7, 12285. https://doi.org/10.1038/ncomms12285 (2016).Article 
ADS 
CAS 

Google Scholar 
Letten, A. D. & Stouffer, D. B. The mechanistic basis for higher-order interactions and non-additivity in competitive communities. Ecol. Lett. 22, 423–436. https://doi.org/10.1111/ele.13211 (2019).Article 

Google Scholar 
Lotka, A. J. Elements of physical biology. Sci. Prog. Twent. Century (1919–1933) 21, 341–343 (1926).
Google Scholar 
Volterra, V. Variazioni e Fluttuazioni del Numero d’Individui in Specie Animali Conviventi. (Società Anonima Tipografica “Leonardo da Vinci”, 1926).Schmitz, O. J. Top predator control of plant biodiversity and productivity in an old-field ecosystem. Ecol. Lett. 6, 156–163. https://doi.org/10.1046/j.1461-0248.2003.00412.x (2003).Article 

Google Scholar 
Fey, K., Banks, P. B., Oksanen, L. & Korpimäki, E. Does removal of an alien predator from small islands in the Baltic Sea induce a trophic cascade?. Ecography 32, 546–552. https://doi.org/10.1111/j.1600-0587.2008.05637.x (2009).Article 

Google Scholar 
Terborgh John, W. Toward a trophic theory of species diversity. Proc. Natl. Acad. Sci. U.S.A. 112, 11415–11422. https://doi.org/10.1073/pnas.1501070112 (2015).Article 
ADS 
CAS 

Google Scholar 
Pringle, R. M. et al. Predator-induced collapse of niche structure and species coexistence. Nature 570, 58–64. https://doi.org/10.1038/s41586-019-1264-6 (2019).Article 
ADS 
CAS 

Google Scholar 
Sandom, C. et al. Mammal predator and prey species richness are strongly linked at macroscales. Ecology 94, 1112–1122. https://doi.org/10.1890/12-1342.1 (2013).Article 

Google Scholar 
Louette, G. & De Meester, L. Predation and priority effects in experimental zooplankton communities. Oikos 116, 419–426. https://doi.org/10.1111/j.2006.0030-1299.15381.x (2007).Article 

Google Scholar 
Johnston, N. K., Pu, Z. & Jiang, L. Predator identity influences metacommunity assembly. J. Anim. Ecol. 85, 1161–1170. https://doi.org/10.1111/1365-2656.12551 (2016).Article 

Google Scholar 
Karakoc, C., Radchuk, V., Harms, H. & Chatzinotas, A. Interactions between predation and disturbances shape prey communities. Sci. Rep. 8, 2968. https://doi.org/10.1038/s41598-018-21219-x (2018).Article 
ADS 
CAS 

Google Scholar 
Hubbell, S. P. The Unified Neutral Theory of Biodiversity and Biogeography (MPB-32) (Princeton University Press, 2011).Book 

Google Scholar 
MacArthur, R. H. & Wilson, E. O. The Theory of Island Biogeography (Princeton University Press, 2001).Book 

Google Scholar 
Daniel, J., Gleason, J. E., Cottenie, K. & Rooney, R. C. Stochastic and deterministic processes drive wetland community assembly across a gradient of environmental filtering. Oikos 128, 1158–1169. https://doi.org/10.1111/oik.05987 (2019).Article 

Google Scholar 
Lehner, B. & Döll, P. Development and validation of a global database of lakes, reservoirs and wetlands. J. Hydrol. 296, 1–22. https://doi.org/10.1016/j.jhydrol.2004.03.028 (2004).Article 
ADS 

Google Scholar 
Chase, J. M., Biro, E. G., Ryberg, W. A. & Smith, K. G. Predators temper the relative importance of stochastic processes in the assembly of prey metacommunities. Ecol. Lett. 12, 1210–1218. https://doi.org/10.1111/j.1461-0248.2009.01362.x (2009).Article 

Google Scholar 
Werner, E. E. & Peacor, S. D. A review of trait-mediated indirect interactions in ecological communities. Ecology 84, 1083–1100. https://doi.org/10.1890/0012-9658(2003)084[1083:AROTII]2.0.CO;2 (2003).Article 

Google Scholar 
Pearson, D. E., Ortega, Y. K., Eren, Ö. & Hierro, J. L. Community assembly theory as a framework for biological invasions. Trends Ecol. Evol. 33, 313–325. https://doi.org/10.1016/j.tree.2018.03.002 (2018).Article 

Google Scholar 
Duchesne, É. et al. Variable strength of predator-mediated effects on species occurrence in an arctic terrestrial vertebrate community. Ecography 44, 1236–1248. https://doi.org/10.1111/ecog.05760 (2021).Article 

Google Scholar 
Ryberg, W. A., Smith, K. G. & Chase, J. M. Predators alter the scaling of diversity in prey metacommunities. Oikos 121, 1995–2000. https://doi.org/10.1111/j.1600-0706.2012.19620.x (2012).Article 

Google Scholar 
Carrete Vega, G. & Wiens, J. J. Why are there so few fish in the sea?. Proc. R. Soc. B 279, 2323–2329. https://doi.org/10.1098/rspb.2012.0075 (2012).Article 

Google Scholar 
Barrett, M. et al. Living planet report 2018: Aiming higher. (2018).Reid, A. J. et al. Emerging threats and persistent conservation challenges for freshwater biodiversity. Biol. Rev. 94, 849–873. https://doi.org/10.1111/brv.12480 (2019).Article 

Google Scholar 
Di Marco, M. et al. Changing trends and persisting biases in three decades of conservation science. Glob. Ecol. Conserv. 10, 32–42. https://doi.org/10.1016/j.gecco.2017.01.008 (2017).Article 

Google Scholar 
Hammerschlag, N. et al. Ecosystem function and services of aquatic predators in the anthropocene. Trends Ecol. Evol. 34, 369–383. https://doi.org/10.1016/j.tree.2019.01.005 (2019).Article 

Google Scholar 
Wang, T. et al. Amur tigers and leopards returning to China: direct evidence and a landscape conservation plan. Landsc Ecol 31, 491–503. https://doi.org/10.1007/s10980-015-0278-1 (2016).Article 

Google Scholar 
Hong, S. et al. Stream health, topography, and land use influences on the distribution of the Eurasian otter Lutra lutra in the Nakdong River basin, South Korea. Ecol. Indic. 88, 241–249. https://doi.org/10.1016/j.ecolind.2018.01.004 (2018).Article 

Google Scholar 
Guter, A., Dolev, A., Saltz, D. & Kronfeld-Schor, N. Using videotaping to validate the use of spraints as an index of Eurasian otter (Lutra lutra) activity. Ecol. Indic. 8, 462–465. https://doi.org/10.1016/j.ecolind.2007.04.009 (2008).Article 

Google Scholar 
Sittenthaler, M., Bayerl, H., Unfer, G., Kuehn, R. & Parz-Gollner, R. Impact of fish stocking on Eurasian otter (Lutra lutra) densities: A case study on two salmonid streams. Mamm. Biol. 80, 106–113. https://doi.org/10.1016/j.mambio.2015.01.004 (2015).Article 

Google Scholar 
Zheng, B., Huang, H., Zhang, Y. & Dai, D. The Fishes of Tumen River (Jilin People’s Publishing House, 1980).
Google Scholar 
Fleishman, E., Murphy, D. D. & Brussard, P. F. A new method for selection of umbrella species for conservation planning. Ecol Appl 10, 569–579. https://doi.org/10.1890/1051-0761(2000)010[0569:ANMFSO]2.0.CO;2 (2000).Article 

Google Scholar 
Roberge, J.-M. & Angelstam, P. E. R. Usefulness of the umbrella species concept as a conservation tool. Conserv. Biol. 18, 76–85. https://doi.org/10.1111/j.1523-1739.2004.00450.x (2004).Article 

Google Scholar 
McGowan, J. et al. Conservation prioritization can resolve the flagship species conundrum. Nat. Commun. 11, 994. https://doi.org/10.1038/s41467-020-14554-z (2020).Article 
ADS 
CAS 

Google Scholar 
Katano, I., Doi, H., Eriksson, B. K. & Hillebrand, H. A cross-system meta-analysis reveals coupled predation effects on prey biomass and diversity. Oikos 124, 1427–1435. https://doi.org/10.1111/oik.02430 (2015).Article 

Google Scholar 
Leibold, M. A. A graphical model of keystone predators in food webs: Trophic regulation of abundance, incidence, and diversity patterns in communities. Am. Nat. 147, 784–812. https://doi.org/10.1086/285879 (1996).Article 

Google Scholar 
McPeek, M. A. The consequences of changing the top predator in a food web: A comparative experimental approach. Ecol. Monogr. 68, 1–23. https://doi.org/10.1890/0012-9615(1998)068[0001:TCOCTT]2.0.CO;2 (1998).Article 

Google Scholar 
Chase, J. M. & Leibold, M. A. Ecological Niches: Linking Classical and Contemporary Approaches (University of Chicago Press, 2003).Book 

Google Scholar 
Gravel, D., Canham, C. D., Beaudet, M. & Messier, C. Reconciling niche and neutrality: The continuum hypothesis. Ecol. Lett. 9, 399–409. https://doi.org/10.1111/j.1461-0248.2006.00884.x (2006).Article 

Google Scholar 
Yoshida, T., Jones, L. E., Ellner, S. P., Fussmann, G. F. & Hairston, N. G. Rapid evolution drives ecological dynamics in a predator–prey system. Nature 424, 303–306. https://doi.org/10.1038/nature01767 (2003).Article 
ADS 
CAS 

Google Scholar 
Yin, X., Wang, J., Yin, H. & Ruan, Y. Does inducible defense mitigate physiological stress responses of prey to predation risk?. Hydrobiologia 843, 173–181. https://doi.org/10.1007/s10750-019-04046-7 (2019).Article 

Google Scholar 
Chalcraft, D. R. & Resetarits, W. J. Jr. Predator identity and ecological impacts: Functional redundancy or functional diversity?. Ecology 84, 2407–2418. https://doi.org/10.1890/02-0550 (2003).Article 

Google Scholar 
Petchey, O. L. & Gaston, K. J. Functional diversity: Back to basics and looking forward. Ecol. Lett. 9, 741–758. https://doi.org/10.1111/j.1461-0248.2006.00924.x (2006).Article 

Google Scholar 
Burner, R. C. et al. Functional structure of European forest beetle communities is enhanced by rare species. Biol. Conserv. 267, 109491. https://doi.org/10.1016/j.biocon.2022.109491 (2022).Article 

Google Scholar  More

Institute of Carbon Neutrality, Sino-French Institute for Earth System Science, College of Urban and Environmental Sciences, Peking University, Beijing, ChinaQian Zhao, Yao Zhang & Shilong PiaoSchool of Urban Planning and Design, Shenzhen Graduate School, Peking University, Shenzhen, ChinaZaichun Zhu & Hui ZengKey Laboratory of Earth Surface System and Human—Earth Relations, Ministry of Natural Resources of China, Shenzhen Graduate School, Peking University, Shenzhen, ChinaZaichun Zhu & Hui ZengDepartment of Earth and Environment, Boston University, Boston, MA, USARanga B. MyneniCSIC, Global Ecology Unit CREAF-CSIC-UAB, Barcelona, Catalonia, SpainJosep PeñuelasCREAF, Barcelona, Catalonia, SpainJosep PeñuelasState Key Laboratory of Tibetan Plateau Earth System, Resources and Environment (TPESRE), Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing, ChinaShilong Piao More

Identification of broad-host-range rhizoplane colonization genes by Tn-seqThis work was focused on SF2 harboring a typical multipartite genome of Sinorhizobium (chromosome, chromid, and symbiosis plasmid) [59]. To perform genome-wide survey of rhizoplane colonization genes of SF2 (Fig. 1), the input mutant library was inoculated on filter paper of plant culture dish, and output mutant libraries were collected from filter papers at 1 h post inoculation (F1h) and 7 days post inoculation (dpi; F7d), and from rhizoplane of cultivated soybean (CS7d), wild soybean (WS7d), rice (R7d), and maize (Z7d) at 7 dpi. To facilitate Tn-seq library construction, all output mutant libraries were subject to 32 h cultivation in the TY rich medium, with input libraries cultivated at the same condition as control (TY). Tn-seq revealed that transposon insertion density in three input and 21 output samples ranged from 57.03 to 86.99% (Table S3), which are above the threshold of 50% insertion density for a good Tn-seq dataset [49]. A reproducible rhizosphere effect was observed in three independent experiments (Fig. S1), i.e., rhizoplane samples (CS7d, WS7d, R7d, and Z7d) consistently formed distinct clusters compared to those of TY, F1h, and F7d. A considerable signature of three independent input libraries was also identified (Data S1, Data S2, and Fig. S1). These results highlight that stochastic variations among multiple independent input libraries should be considered before making conclusions on gene fitness, which has been largely overlooked in earlier studies based on just one input library [49].Based on gene fitness scores of rhizoplane samples (CS7d, WS7d, R7d and Z7d) compared to corresponding F1h datasets (Fig. S2A; Data S2), 93, 91, 127, and 206 genes were identified as rhizoplane colonization genes for test plants of cultivated soybean, wild soybean, maize, and rice, respectively, accounting for 1.4–3.1% of the SF2 genome (p values  More

van der Ent, R. J. & Tuinenburg, O. A. The residence time of water in the atmosphere revisited. Hydrol. Earth Syst. Sci. 21, 779–790 (2017).Article 
ADS 

Google Scholar 
Risi, C., Noone, D., Frankenberg, C. & Worden, J. Role of continental recycling in intraseasonal variations of continental moisture as deduced from model simulations and water vapor isotopic measurements. Water Resour. Res. 49, 4136–4156 (2013).Article 
ADS 
CAS 

Google Scholar 
Gupta, S. K. & Deshpande, R. D. Water for India in 2050: First-order assessment of available options. Curr. Sci. 86, 9 (2004).
Google Scholar 
Kuttippurath, J. et al. Observed rainfall changes in the past century (1901–2019) over the wettest place on Earth. Environ. Res. Lett. 16, 024018 (2021).Article 
ADS 

Google Scholar 
Rao, M. P. et al. Seven centuries of reconstructed Brahmaputra River discharge demonstrate underestimated high discharge and flood hazard frequency. Nat. Commun. 11, 6017 (2020).Article 
ADS 
CAS 

Google Scholar 
Bassi, N., Kumar, M. D., Sharma, A. & Pardha-Saradhi, P. Status of wetlands in India: A review of extent, ecosystem benefits, threats and management strategies. J. Hydrol. Regional Stud. 2, 1–19 (2014).Article 

Google Scholar 
Dikshit, K. R. & Dikshit, J. K. Natural vegetation: Forests and grasslands of North-East India. in North-East India: Land, People and Economy (eds. Dikshit, K. R. & Dikshit, J. K.) 213–255 (Springer, 2014). https://doi.org/10.1007/978-94-007-7055-3_9.Chakraborty, S. et al. Linkage between precipitation isotopes and biosphere-atmosphere interaction observed in northeast India. Npj Clim. Atmos. Sci. 5, 1–11 (2022).Article 

Google Scholar 
Pathak, A., Ghosh, S. & Kumar, P. Precipitation recycling in the Indian subcontinent during summer monsoon. J. Hydrometeorol. 15, 2050–2066 (2014).Article 
ADS 

Google Scholar 
Mahanta, R., Sarma, D. & Choudhury, A. Heavy rainfall occurrences in northeast India. Int. J. Climatol. 33, 1456–1469 (2013).Article 

Google Scholar 
Murata, F. et al. dominant synoptic disturbance in the extreme rainfall at cherrapunji, northeast India, based on 104 years of rainfall data (1902–2005). J. Clim. 30, 8237–8251 (2017).Article 
ADS 

Google Scholar 
Roy, S. C. & Chatterji, G. Origin of nor’westers. Nature 124, 481–481 (1929).Article 
ADS 

Google Scholar 
Dhar, O. N. & Nandargi, S. A study of floods in the Brahmaputra basin in India. Int. J. Climatol. 20, 771–781 (2000).Article 

Google Scholar 
Reager, J. T., Thomas, B. F. & Famiglietti, J. S. River basin flood potential inferred using GRACE gravity observations at several months lead time. Nat. Geosci. 7, 588–592 (2014).Article 
ADS 
CAS 

Google Scholar 
Gat, J. R. Isotope Hydrology: A Study of the Water Cycle (World Scientific, 2010).Book 

Google Scholar 
Salati, E., DallOlio, A., Matsui, E. & Gat, J. R. Recycling of water in the Amazon basin: An isotopic study. Water Resourc. Res. 15, 1250–1258 (1979).Article 
ADS 
CAS 

Google Scholar 
Victoria, R. L., Martinelli, L. A., Mortatti, J. & Richey, J. Mechanisms of water recycling in the Amazon basin: Isotopic insights. Ambio 20, 384–387 (1991).
Google Scholar 
Wright, J. S. et al. Rainforest-initiated wet season onset over the southern Amazon. PNAS 114, 8481–8486 (2017).Article 
ADS 
CAS 

Google Scholar 
Leite-Filho, A. T., Soares-Filho, B. S., Davis, J. L., Abrahão, G. M. & Börner, J. Deforestation reduces rainfall and agricultural revenues in the Brazilian Amazon. Nat. Commun. 12, 2591 (2021).Article 
ADS 
CAS 

Google Scholar 
Spracklen, D. V. & Garcia-Carreras, L. The impact of Amazonian deforestation on Amazon basin rainfall. Geophys. Res. Lett. 42, 9546–9552 (2015).Article 
ADS 

Google Scholar 
Lele, N. & Joshi, P. K. Analyzing deforestation rates, spatial forest cover changes and identifying critical areas of forest cover changes in North-East India during 1972–1999. Environ. Monit. Assess. 156, 159 (2008).Article 

Google Scholar 
Sudhakar Reddy, C. et al. Quantification and monitoring of deforestation in India over eight decades (1930–2013). Biodivers. Conserv. 25, 93–116 (2016).Article 

Google Scholar 
Kathayat, G. et al. Interannual oxygen isotope variability in Indian summer monsoon precipitation reflects changes in moisture sources. Commun. Earth Environ. 2, 1–10 (2021).Article 

Google Scholar 
Kathayat, G. et al. Protracted Indian monsoon droughts of the past millennium and their societal impacts. Proc. Natl. Acad. Sci. 119, e2207487119 (2022).Article 
CAS 

Google Scholar 
Breitenbach, S. F. M. et al. Strong influence of water vapor source dynamics on stable isotopes in precipitation observed in Southern Meghalaya, NE India. Earth Planet. Sci. Lett. 292, 212–220 (2010).Article 
ADS 
CAS 

Google Scholar 
Pradhan, R., Singh, N. & Singh, R. P. Onset of summer monsoon in Northeast India is preceded by enhanced transpiration. Sci. Rep. 9, 18646 (2019).Article 
ADS 
CAS 

Google Scholar 
Das, S., Sarkar, A., Das, M. K., Rahman, Md. M. & Islam, Md. N. Composite characteristics of Nor’westers based on observations and simulations. Atmos. Res. 158–159, 158–178 (2015).Article 

Google Scholar 
Vinay Kumar, P. & Venkateswara Naidu, C. Is pre-monsoon rainfall activity over India increasing in the recent era of global warming?. Pure Appl. Geophys. 177, 4423–4442 (2020).Article 
ADS 

Google Scholar 
Narayanan, P., Sarkar, S., Basistha, A. & Sachdeva, K. Trend analysis and forecast of pre-monsoon rainfall over India. Weather 71, 94–99 (2016).Article 
ADS 

Google Scholar 
Stevens, B. Atmospheric moist convection. Annu. Rev. Earth Planet. Sci. 33, 605–643 (2005).Article 
ADS 
MathSciNet 
CAS 
MATH 

Google Scholar 
Bershaw, J. Controls on deuterium excess across Asia. Geosciences 8, 257 (2018).Article 
ADS 

Google Scholar 
Clark, I. D. & Fritz, P. Environmental Isotopes in Hydrogeology (CRC Press, 1997). https://doi.org/10.1201/9781482242911.Book 

Google Scholar 
Cui, J., Tian, L., Biggs, T. W. & Wen, R. Deuterium-excess determination of evaporation to inflow ratios of an alpine lake: Implications for water balance and modeling. Hydrol. Process. 31, 1034–1046 (2017).Article 
ADS 

Google Scholar 
Laskar, A. H. et al. Stable isotopic characterization of Nor’westers of southern Assam, NE India. J. Clim. Chang. 1, 75–87 (2015).Article 

Google Scholar 
Tanoue, M. et al. Seasonal variation in isotopic composition and the origin of precipitation over Bangladesh. Prog. Earth Planet Sci. 5, 77 (2018).Article 
ADS 

Google Scholar 
Oza, H., Ganguly, A., Padhya, V. & Deshpande, R. D. Hydrometeorological processes and evaporation from falling rain in Indian sub-continent: Insights from stable isotopes and meteorological parameters. J. Hydrol. 591, 125601 (2020).Article 
CAS 

Google Scholar 
Chen, F. et al. Relationship between sub-cloud secondary evaporation and stable isotopes in precipitation of Lanzhou and surrounding area. Quatern. Int. 380–381, 68–74 (2015).Article 

Google Scholar 
Sinha, N. et al. Isotopic investigation of the moisture transport processes over the Bay of Bengal. J. Hydrol. X 2, 100021 (2019).Article 
CAS 

Google Scholar 
Rawson, H. M., Begg, J. E. & Woodward, R. G. The effect of atmospheric humidity on photosynthesis, transpiration and water use efficiency of leaves of several plant species. Planta 134, 5–10 (1977).Article 
CAS 

Google Scholar 
Chakraborty, S., Belekar, A. R., Datye, A. & Sinha, N. Isotopic study of intraseasonal variations of plant transpiration: An alternative means to characterise the dry phases of monsoon. Sci. Rep. 8, 8647 (2018).Article 
ADS 
CAS 

Google Scholar 
Sinha, N., Chakraborty, S. & Mohan, P. M. Modern rain-isotope data from Indian island and the mainland on the daily scale for the summer monsoon season. Data Brief 23, 103793 (2019).Article 

Google Scholar 
Grujic, D. et al. Formation of a Rain Shadow: O and H stable isotope records in authigenic clays from the Siwalik group in eastern Bhutan. Geochem. Geophys. Geosyst. 19, 3430–3447 (2018).Article 
CAS 

Google Scholar 
Lambs, L., Balakrishna, K., Brunet, F. & Probst, J. L. Oxygen and hydrogen isotopic composition of major Indian rivers: A first global assessment. Hydrol. Process. 19, 3345–3355 (2005).Article 
ADS 
CAS 

Google Scholar 
Ek, M. B. & Holtslag, A. A. M. Influence of soil moisture on boundary layer cloud development. J. Hydrometeorol. 5, 86–99 (2004).Article 
ADS 

Google Scholar 
Findell, K. L. & Eltahir, E. A. B. Atmospheric controls on soil moisture-boundary layer interactions. Part I: Framework development. J. Hydrometeorol. 4, 552–569 (2003).Article 
ADS 

Google Scholar 
Syroka, J. & Toumi, R. On the withdrawal of the Indian summer monsoon. Q. J. R. Meteorol. Soc. 130, 989–1008 (2004).Article 
ADS 

Google Scholar 
Bhatta, L. D. et al. Ecosystem service changes and livelihood impacts in the maguri-motapung wetlands of Assam. India. Land 5, 15 (2016).Article 

Google Scholar 
Choudhury, B. A., Saha, S. K., Konwar, M., Sujith, K. & Deshamukhya, A. Rapid drying of northeast India in the last three decades: Climate change or natural variability?. J. Geophys. Res. Atmos. 124, 227–237 (2019).Article 
ADS 

Google Scholar 
Das, D. Changing climate and its impacts on Assam, Northeast India. Bandung J. Glob. South 2, 26 (2016).Article 

Google Scholar 
Deka, R. L., Mahanta, C., Pathak, H., Nath, K. K. & Das, S. Trends and fluctuations of rainfall regime in the Brahmaputra and Barak basins of Assam, India. Theor. Appl. Climatol. 114, 61–71 (2013).Article 
ADS 

Google Scholar 
Maurya, A. S., Shah, M., Deshpande, R. D. & Gupta, S. K. Protocol for δ18O and δD analyses of water sample using Delta V plus IRMS in CF Mode with Gas Bench II for IWIN National Programme at PRL, Ahmedabad. in 11th ISMAS Triennial Conference of Indian Society for Mass Spectrometry vol. 314, 314–317 (Indian Society for Mass Spectrometry Hyderabad, 2009).Deshpande, R. D. & Gupta, S. K. Oxygen and hydrogen isotopes in hydrological cycle: new data from IWIN national programme. Proc. Indian Natl. Sci. Acad. 78, 321–331 (2012).CAS 

Google Scholar 
Deshpande, R. D. & Gupta, S. K. National programme on isotope fingerprinting of waters of India (IWIN). Glimpses of Geosciences Research in India, the Indian Report to IUGS, Indian National Science Academy (eds Singhvi, AK, Bhattacharya, A. & Guha, S.), 10–16 (2008).Oza, H., Padhya, V., Ganguly, A. & Deshpande, R. D. Investigating hydrometeorology of the Western Himalayas: Insights from stable isotopes of water and meteorological parameters. Atmos. Res. 268, 105997 (2022).Article 
CAS 

Google Scholar 
Stein, A. F. et al. NOAA’s HYSPLIT atmospheric transport and dispersion modeling system. Bull. Am. Meteor. Soc. 96, 2059–2077 (2015).Article 
ADS 

Google Scholar 
Sodemann, H., Schwierz, C. & Wernli, H. Interannual variability of Greenland winter precipitation sources: Lagrangian moisture diagnostic and North Atlantic Oscillation influence. J. Geophys. Res. Atmos. 113, D3 (2008).
Google Scholar 
Oza, H. et al. Hydrometeorological processes in semi-arid western India: insights from long term isotope record of daily precipitation. Clim. Dyn. 54, 2745–2757 (2020).Article 

Google Scholar 
Su, L., Yuan, Z., Fung, J. C. H. & Lau, A. K. H. A comparison of HYSPLIT backward trajectories generated from two GDAS datasets. Sci. Total Environ. 506–507, 527–537 (2015).Article 
ADS 

Google Scholar 
Ahmed, M., Seraj, R. & Islam, S. M. S. The k-means algorithm: A comprehensive survey and performance evaluation. Electronics 9, 1295 (2020).Article 

Google Scholar  More

Geographical and temporal settingPaleogeographic and biogeographic evidence suggests Madagascar has been an isolated insular unit since the splitting from Greater India around 88 million years ago (Mya)1,57. However, a worldwide mass extinction event is known to have taken place around the K-Pg boundary ~66 Mya affecting Madagascar’s biodiversity2, so using this event as an island age for the current biota may also be appropriate. While Madagascar had a rich Mesozoic vertebrate fauna, including many mammals58,59, molecular phylogenetic data suggest that no mammalian lineage that colonized before the K-T event has survived until the present, although the confidence intervals for the colonization time of the Malagasy Afrosoricida extend until before this age (Fig. 1). Furthermore, the existence of short-lived land-bridges connecting the island to continental Africa at various stages has been proposed60, although this hypothesis has been contested61. We accounted for the uncertainty in island age in our sensitivity analyses, using both 88 and 66 million years (Myr) as island age. We do not consider the possible existence of temporary land-bridges, as we estimate average rates (e.g. of colonization) throughout the entire history of the island (see “DAISIE models” section).Malagasy mammalian diversityWe compiled a comprehensive taxonomic/phylogenetic dataset of the entire assemblage of Malagasy non-marine mammals, including information on phylogenetic relationships, timing and causes of extinction and levels of threat (Supplementary Data S1). We first compiled a checklist of all mammalian species known to have been present on Madagascar in the late Holocene before human arrival and expansion, ~2500 years ago (Supplementary Data S1). The checklist includes all native species that are still extant today and all those that are known or presumed to have gone extinct in the late Holocene. Taxa known only from ancient early- or pre-Holocene fossils that are assumed to have gone extinct long before human presence were excluded, as well as all non-native species. We followed the taxonomy and nomenclature of the Mammal Diversity Database of the American Society of Mammalogists62,63, as of May 2022. Our checklist does differ from that database because of ongoing taxonomic and nomenclatural revisions based on recent molecular phylogenetic analyses, and recent discoveries or descriptions of new species on the island. All such cases and cases where our taxonomy or nomenclature differs from that used in the phylogeny of43 are explained in the column “Taxonomy note” in Supplementary Data S1. Molecular-based taxonomy and morphological taxonomy (such as that from the palaeontological record) can be incongruent because cryptic species “detected” in DNA-based analyses may not be identifiable based on fossil data. The number of species in the existing fossil record is therefore likely an underestimate—we address this in the section “Impact of increased knowledge on the ERT”.To compile the checklist, we used a variety of sources from both the neontological and palaeontological literature. The main sources on Malagasy extant and extinct mammalian species are4,8,25,29,52,57,58, but other published studies were used, particularly for bats, lemurs and cryptic and recently discovered species (Supplementary Data S1). We classified species as endemic or non-endemic to Madagascar, as information on endemicity status is one of the types of data that DAISIE uses to estimate rates of speciation. For non-endemic species, only represented in the dataset by certain species of bats, we noted their range outside of Madagascar in the column “Additional range note”. For the ERT analyses, we also compiled the IUCN Red List status for each species in 2010, 2015 and 2021. We used the digital archive “wayback machine” to obtain the 2010 and 2015 IUCN Red List data, as older versions of listings are not kept online by the IUCN8,41,42. We compiled the checklist using Excel v16.63.Extinct speciesIn the DAISIE analyses we treat species that went extinct due to non-anthropogenic causes (before or after human arrival) as if they were species not known to science. This is because the natural extinction rate that is estimated based on the colonization and branching times extracted from phylogenies (without these extinct species), already accounts for such missing species. These include species that went extinct before the late Holocene and may or may not be known from the fossil record, but also species that have gone naturally extinct after human arrival. In contrast, we consider that anthropogenic extinctions do not contribute to the natural extinction rate. Therefore, we treat species for which an anthropogenic cause of extinction is likely as if they had survived into the present and we include them in the phylogenies, following the approach of Valente et al.22. The rates of speciation, colonization and natural extinction that are estimated from such phylogenies are the natural average rates assuming that humans had no impact on the island—these would be the natural average background rates in the periods pre-dating human arrival.Our checklist includes all species that are hypothesized to have gone extinct in the late Holocene, hereafter termed “recently extinct species”. For these species, we compiled information on whether they have gone extinct before or after human arrival (columns “Last date 14C age BP” and “Extinction before/after humans” in Supplementary Data S1). In assessing whether species went extinct before or after human colonization and expansion, we employ the circa 2500 years BP date as time zero3,4,5,64,65. After permanent settlement, anthropogenic pressures on Madagascar’s biodiversity have intensified, with a visible increase in the past few decades4. Species which are considered in the literature to have gone extinct before human arrival were excluded. For some taxa, it is unclear whether they went extinct before or after human colonization, because the fossil record is insufficiently known. These were included in the list and the effect of their inclusion/exclusion was evaluated in sensitivity analyses (see section below). In addition, for all extinct species in our list, we compiled information on the hypothesized causes of extinction cited in the literature, classifying each extinction as anthropogenic, natural or uncertain. References for timing and causes of extinction are provided in Supplementary Data S1.To account for these uncertainties, we re-ran analyses for two datasets, assuming high and low human impact (see sensitivity analyses). For the high human impact scenario, we assumed all recent extinctions (less than 2500 years ago) to have an anthropogenic cause and therefore included them in the phylogenies and in the counts of pre-human species diversity. In this scenario we also include recently extinct species for which no 14C dates are available, but which have been hypothesized to have gone extinct in the last 2500 years, as well as species with a putative natural cause of extinction, because even natural recent extinctions may have had an indirect human influence. For the high human impact scenario, we thus assume all recent species loss is linked to potential human influence, and we include all those species in the phylogenies. In the low human impact scenario we assume a natural cause for all recent extinctions that have previously been hypothesized to have had a natural cause, for all recent extinctions whose cause is unknown, and for the cases for which it is unclear whether extinction pre-dates or post-dates human arrival (e.g. no radiocarbon dating available). For the low human impact scenario, we excluded all such species from the phylogenies and from the counts of pre-human diversity (it is assumed they are unknown). See column “Low human impact scenario” in Supplementary Data S1. For two of the species for which no date currently exists indicating whether they went extinction before or after humans (the lemurs Mesopropithecus dolichobrachion and Palaeopropithecus kelyus), an anthropogenic cause of extinction has been hypothesized despite the lack of a precise last occurrence date, and we thus considered them to be anthropogenic extinctions in both high and low human impact scenarios.Phylogenetic dataThe source of all our phylogenetic data—including divergence times of Malagasy lineages—is the phylogeny by Upham et al.43, the most comprehensive and complete mammalian phylogeny published to date, including 5911 species of mammals. From this tree, we extracted phylogenetic information with reference to Madagascar on the number of colonization events, the estimated dates of colonization (divergence times from the most closely related non-Malagasy relatives), number of species per monophyletic colonist lineage and the timing of within-island speciation events. We created a Madagascar-specific dataset consisting of a series of multiple subtrees drawn from the same Mammalia-wide dating framework, representing all colonization events for most known late Holocene native mammals on the island, including bats and recently extinct species. We visually inspected the trees using Figtree v1.4.4.Upham et al.43 used two approaches to calibrate their phylogeny: node dating and tip dating. Following the recommendations in that publication, we used the trees based on the node-dating approach, in which node-age priors were placed on the tree based on 17 mammalian fossils and one root constraint. Regarding molecular sampling, they produced two types of trees: DNA-only, with 4098 species sampled in the phylogeny based on molecular data; and completed trees, where they placed an additional 1813 species that were unsampled for DNA in the tree using taxonomic constraints (across multiple posterior trees). The DNA-only trees have the advantage that the topology is based on molecular data, and is likely more reliable, but the disadvantage that DNA sequences were not available for many Malagasy species and so these species needed to be added to the phylogeny for the DAISIE analyses. The completed trees have the advantage that they are near-complete, but the disadvantage that some Malagasy species—particularly several bats that are unsampled for DNA—were placed randomly within a given clade constraint, which may lead some trees in the posterior to have some incorrectly inferred colonizations. We ran analyses on data extracted from both types of trees (see sensitivity analyses).An alternative to using this phylogenetic dataset would be to extract data from separate individual trees from publications with phylogenies focusing on specific clades. There are many such studies, and indeed some of them include taxa that are not present in the Upham et al.43 tree—e.g. new recently described cryptic species that were only identified after molecular analyses, including, for example, the nesomyine rodent Eliurus tsingimbato66 and the mouse lemur species Microcebus jonahi67; or extinct species for which no molecular data exists, but which were included in phylogenetic dating analyses based on morphological data, such as members of the lemur genus Mesopropithecus46. However, we favoured using phylogenetic data from a single study to ensure divergence times are comparable (i.e. same models, assumptions and data), even though this is done at the expense of reduced species sampling. Although we use a single tree (or posterior distribution of single trees) for our dataset, DAISIE treats each Malagasy colonizing lineage as its own separate tree, so we deal with a “forest” of phylogenetic trees, each representing a single Malagasy lineage resulting from one colonization event. For example, the lineages that have radiated on the island have a tree that includes the stem age of the lineage (splitting from the closest sampled mainland relative) and all branching events within the radiation. Lineages with a single species on Madagascar (endemic or non-endemic) are essentially a tree with a single tip and with an age equal to the splitting of that species from its closest (sampled) continental relative.Alternative colonization scenariosThe number of colonization events of Malagasy mammals inferred from the phylogenetic data can vary depending on the placement of some missing taxa in the tree or because some clades have poor branch support and could be the result of one or more colonizations in different trees from the posterior. We considered two alternative colonization scenarios (CS), one where we favour fewer colonizations (CS1) and one where we favour more (CS2). The differences between the two scenarios are summarized in Table S1 (all colonizations shown in Supplementary Data S2, S3). We considered lemurs to be the result of a single colonization event in both scenarios. A recent study68 has suggested that Daubentonia is a separate colonization of Madagascar, but we assigned both the extant aye-aye and the extinct giant aye-aye to the single Lemuroidea (lemurs) clade because that is the only scenario supported by the mammal tree.Adding missing speciesA total of 34 out of 249 species in our Madagascar mammal checklist are not present in the mammal phylogenetic tree43. Most of these (23 species) are extinct species (Data S1). The other 11 species are recently described species (two bats, eight lemurs and one nesomyine rodent), these are indicated in column “Taxonomy note” in Supplementary Data S1. An additional 61 species are included in the completed trees, but not in the DNA-only trees, as no molecular data were available for these. We added the 34 species missing from the mammal tree to both DNA-only and completed trees, and the 61 species missing molecular data to the DNA-only trees. Instead of adding species directly to the posterior distribution of trees and then extracting information from the phylogenies, we assign those species to specific clades using the “missing species” option in DAISIE. This tool allows them to be placed anywhere within the Malagasy clade they are believed to belong to, without specifying a specific topological position within the clade—DAISIE does not use topological information for its estimates. For example, a species of lemur that was missing from the tree was added to the species count of the lemur clade. The information on the clade to which each missing species was added to (under either CS1 or CS2) is provided in Supplementary Data S1.Most recently extinct species are not included in the mammal tree because the original study was primarily focused on the extant mammalian taxa43. Three extinct species of lemur, Archaeolemur majori, Megaladapis edwardsi and Palaeopropithecus ingens are included in their tree based on molecular data obtained from subfossil material. One extinct species of lemur (A. edwardsi), one extinct species of carnivoran (Cryptoprocta spelea) and two extinct hippopotamus species (Hippopotamus madagascariensis and H. lemerlei) are included in their completed trees, i.e. not based on molecular data. We added the remaining extinct species to the phylogenies using the approach explained above (Supplementary Data S1). These were: 1 tenrec (assigned to the Malagasy Afrosoricida (CS1) or Tenrecidae (CS2)); the 2 bibymalagasy (assigned to Malagasy Afrosoricida (CS1) or to Bibymalagasia (CS2)); 1 hippopotamus (assigned to the single hippopotamus clade (CS1) or to one of the two hippopotamus clades (CS2); 1 euplerid carnivore (Cryptoprocta sp. nov., assigned to Eupleridae); 2 bats (1 Paratriaenops, assigned to the Paratriaenops clade; 1 Macronycteris assigned to Macronycteris (CS1) or as its own colonization (CS2)); 3 nesomyine rodents (assigned to Nesomyinae); and 13 lemurs (assigned to the Lemuroidea clade).In a few cases, all descendants from a colonization of Madagascar were missing from the mammal tree. These were added as a separate colonization, using the DAISIE_max_age option, which assumes that they could have colonized at any time since the given age and the present. These were: the two species of Bibymalagasia (CS2, using the stem age of Afrosoricida in the mammal tree as the maximum age of colonization); Chaerephon leucogaster (CS1 and CS2, using crown age of Molossidae family as maximum colonization time); Macronycteris cryptovalorona (CS1 and CS2, using crown age of Hipposideridae family as maximum colonization time); Macronycteris besaoka (CS2, using crown age of Hipposideridae family as maximum colonization time); Hippopotamus laloumena (CS2, using stem age of genus Hippopotamus as maximum colonization time); and Pipistrellus raceyi (absent from the DNA-only tree, we used the crown age of Vespertilionidae as the maximum colonization time). For Miniopterus, the phylogenetic resolution for this radiation is poor (including both Malagasy and non-Malagasy taxa), and we therefore used the crown age of the genus as a maximum colonization time of Madagascar (we chose the crown and not stem because Miniopterus of Madagascar do not diverge early in the genus).Colonization and branching timesFor endemic Malagasy clades (radiations (e.g. lemurs) or clades with a single endemic species, e.g. Pteropus rufus), we assumed the time of colonization of Madagascar coincides with the divergence time from its closest non-Malagasy lineage, i.e. the stem age of the clade. These ages are likely overestimates (e.g. if the tree is incompletely sampled, or if the closest continental ancestor has gone extinct, see ref. 18), but are a good approximation, and we repeated analyses over the posterior distribution of trees to account for age uncertainties. Non-endemic species are represented by a single tip in the mammal tree, and we therefore used the age of that tip as a maximum age of colonization, as the actual colonization time of the Madagascar population is most certainly younger than that age. DAISIE integrates through all possible ages between that maximum age and the present. The only exception are three non-endemic species belonging to the Chiroptera genus Miniopterus, which likely resulted from cladogenesis within Madagascar and became non-endemic by colonizing the Comoros. These are treated as part of the Miniopterus clade or clades (Table S1), and thus contribute to the estimates of cladogenesis on Madagascar (rather than being assigned their own colonist lineage). The branching times within Madagascar radiations were taken directly from the trees. When species within a radiation were missing from the phylogeny, they were included using the DAISIE missing species option (see section above), thus contributing to the estimates of cladogenesis rates for the given clade.We wrote an R script to extract colonization and branching times from the maximum clade credibility (MCC) and posterior trees from the Upham et al.43 mammal phylogeny, add missing species and assign maximum colonization times (if relevant), assuming a variety of scenarios (see “Main dataset and sensitivity analysis” section below). Once the data were extracted from the trees, the script creates DAISIE objects, i.e. datasets in DAISIE format that can be read by DAISIE functions. This script uses functions from phytools v1.2-0 R69, ape v5.6-270 and DAISIE47 v4.0.5 R packages. The R script describes all the steps taken to prepare the phylogenetic data for the DAISIE analyses. The script, the precise source trees that we used from the Upham et al. mammal phylogeny43, as well as all DAISIE objects for the main analyses and the sensitivity analyses are provided in an online repository (https://github.com/luislvalente/madagascar). Analyses were run in R v4.2.1 and RStudio v2022.02.3.DAISIEWe used the DAISIE R package47 to estimate rates of speciation, colonization and extinction (CES rates) of Madagascar mammals using maximum likelihood (ML) under a range of different models and to identify the preferred model given the phylogenetic data. The DAISIE likelihood inference approach is based on theory and methods developed for phylogenetic birth-death models71,72. It has been demonstrated that the shape of phylogenies of extant species contains information about natural extinction rates72. While these approaches have many known limitations73, we have shown in different studies that the DAISIE model is able to accurately estimate extinction rate from simulated datasets for which the extinction rate is known18,74. In addition, unlike most phylogenetic birth-death models, which are single-clade approaches and use only information from branching times, DAISIE has the advantage that it uses information from multiple independent clades and from both colonization and branching times, increasing its statistical power to estimate parameters.We fitted a set of 30 DAISIE models to the phylogenetic data, explained in Table S2. Models M1–M4 assume homogeneous CES rates for all Malagasy mammals, while for models M5–M30 we allow one or more of the CES parameters to vary between non-volant mammals and bats. The set of models include both diversity-dependent and diversity-independent models. Models can differ in the number of parameters: for example, M1 has five parameters (colonization, cladogenesis, anagenesis, K and extinction); M3 has four parameters (same as M1, except that anagenesis is fixed to zero); M5 has six parameters (the same five parameters as M1, plus a parameter for colonization rate which differs for bats); and M6 has five parameters (same as M5 but K is fixed to infinite, i.e. there is no diversity dependence).CES rate heterogeneityDAISIE estimates average CES rates for the island and assumes that these rates are constant through time (except for models that include diversity dependence, in which rates decline with increasing diversity). However, from the geology of the island and the fossil record, we can infer that rates have most likely not been constant. For example, periods of large-scale natural extinction may have taken place throughout the history of the island3,75,76. While there may have been important temporal rate changes, when estimating the future island evolutionary return time (the main purpose of our analyses), we seek to estimate the overall average natural background rates, which incorporate periods of both low and high rates (which will certainly also occur in the future). We fitted two models (M31—bats and non-volant mammals share same rates; and M32—bats and non-volant mammals have different rates) in which colonization rates can shift to a lower or higher rate46 at a certain point in time, but these models were not preferred (Table S3). Therefore, when we estimate the ERTs, we use the average rates for the island as a whole over its entire geological history in the absence of humans. Importantly, although the preferred models assume constant rates, the DAISIE model has been shown to perform very well for ancient continental islands (separated from the mainland very deep in geological time, such as Madagascar), in terms of accurately predicting the number of species, and the number of species and colonizations through time77. In addition, although rates may have been lower or higher at some periods, the average rates are nevertheless informative of the unique geographical setting of the island and the ecological characteristics of the target community—this is particularly valuable, for example, when comparing Malagasy mammals with ERTs from other systems, such as in Caribbean bats22 and New Zealand birds21, both in which rates have also most likely varied through time.There is evidence for rate variation among mammalian lineages78. We therefore chose to test for differential rates for two groups: non-volant mammals and bats. In the context of islands, it is likely that bats will have different rates of colonization due to their higher dispersal abilities, and they may also vary in other parameters. We used the two-type DAISIE model approach first applied to the birds of the Galápagos (Darwin’s finches vs other birds45. While there may also be differences in rates between specific non-volant and bat clades, we favour obtaining average ERTs across the whole fauna, rather than specific ERTs for each lineage. First, assigning unique rates to each lineage would lead to over parameterization, and estimating lineage-specific rates would not be reliable for some individual Malagasy clades that are the product of a single colonization and have few species (e.g. many bat lineages, hippopotamuses, euplerid carnivorans). Thus, we restricted the test of idiosyncrasies to the comparison between bats and non-volant species. Second, an advantage of our approach is that the rates we obtained are based entirely on the phylogenies of Malagasy species and therefore our rates are already very specific to the Malagasy context—whereas comparable methods use average rates worldwide and then extrapolate to the focal lineages24. Third, we are interested in whether total diversity will recover, not whether specific types of species will recover. A trait-dependent diversification model for insular communities that would allow us to obtain ERTs based on, for example, certain morphological traits that may promote diversification, does not currently exist.Main dataset and sensitivity analysesWe consider the main dataset (D1) to comprise: colonization scenario 1 (CS1) with high human impact, using the DNA-only mammal tree (MCC and posterior), island age of 88 Myr. The reason for this is that we consider the CS1 (fewer colonizations) and high human impact scenarios to be the most realistic given the level of isolation of the island and because the evidence for anthropogenic mammalian extinction on Madagascar is compelling and growing6. We also consider the DNA-only tree more appropriate, as species were sampled based on molecular data, and all missing species were included in clades using the DAISIE missing species option, i.e. placing them in a clade but without forcing a given topology within that clade.There are currently ~6500 species of mammals62, but this number was certainly different in the past and only a subset of these constitute the potential mainland pool for Madagascar, which would include African species and to a much lesser extent from the Indian subcontinent or other portions of Asia. For the main dataset we considered the number of species on the mainland pool (M) to be 1000 (approximately the current number of African terrestrial mammal species), but we re-ran analyses with 2000 and 5000 species. For the models where bats differ from non-volant mammals (M5–M30), the proportion of bat species in the mainland pool was set to 0.22, equivalent to the proportion of all mammal species that are bats today.To account for uncertainty in island age, mainland pool size, colonization scenarios, human impact, topology, dating (colonization and branching times), and tree sampling completeness, we ran a series of sensitivity analyses. We re-ran analyses for the MCC tree of the main dataset (that is: DNA-only tree, high impact, CS1), assuming an island age of 66 Myr, and varying pool sizes (for both island ages). Then, fixing the mainland pool to 1000 species and the island age to 88 Myr, we ran DAISIE analyses assuming colonization scenarios CS1 and CS2, high and low human impact. We also repeated analyses using the completed and DNA-only trees, using the corresponding MCC trees for each scenario. In total we ran 13 different scenarios (D1–D13) for the sensitivity analysis (Supplementary Data S4).We used the following approach for ML optimizations on the main dataset and the sensitivity analyses, using the DAISIE_ML function implemented in the DAISE R package. For the analyses on a single MCC tree (all 13 scenarios, including the main dataset), we fitted each of the 30 DAISIE models to each dataset 10 times, using different random sets of starting values for the likelihood optimization (30 × 10 = 300 ML optimizations per scenario, total 3900 ML optimizations). For each scenario, we selected the preferred model by comparing Bayesian information criterion (BIC) and Akaike information criterion (AIC) scores between models. For the main dataset, to examine if the same model is preferred across the posterior distribution of trees, we also ran analyses on the posterior, fitting each model 4 times to each of 100 datasets from the posterior (30 × 4 × 100 = 12,000 optimizations). To obtain confidence intervals for the preferred model of the main dataset, we ran analyses on 1000 trees from the posterior, with 2 random sets of starting values (2 × 1000 = 2000 ML optimizations, Table S4). We ran ERT analyses using the parameters of the preferred model for all 13 scenarios. For the main dataset D1 we also fitted two models with a temporal shift in colonization rate (M31 and M32, Table S3). All analyses were run on the Peregrine cluster of the University of Groningen.For the main dataset, we ran simulations of the best overall rate model using the DAISIE_sim function. Under the parameters of the model, we simulated 5000 islands for 88 million years. We then assessed the goodness of fit of the model to the data by comparing diversity metrics in the simulated datasets to those in the empirical data.In the sensitivity analyses, varying mainland pool size, island age, human impact, colonization scenario or phylogenetic dataset (DNA-only vs completed), had a limited impact in the preferred models or parameters values (Supplementary Data S4). Varying island age, human impact or colonization scenario generally led to only minor changes in parameter values. Varying mainland pool sizes affected the colonization rate, which decreases with mainland pool size because colonization rate is measured per mainland species. When using BIC as the criterion for model selection, M26 was the preferred model in 10 out of 13 scenarios, with M22 being the preferred model under one scenario (D8, DNA-only data, island age 88, M = 1000, CS2, high human impact), and M11 preferred under two scenarios (D12 and D13, completed trees, island age 88, M = 1000, C2, for both high and low human impact) (Supplementary Data S4). Using AIC, alternative models to M26 were preferred for two additional scenarios – M10 was preferred for D10 (completed trees, island age 88, M = 1000, CS1, high human impact) and M11 was preferred for D11 (completed trees, island age 88, M = 1000, CS1, low human impact). We consider the M26 model to be the preferred model overall, because we favour the DNA-only trees (for which M26 was consistently selected as the best model under both AIC and BIC) and the BIC criterion for model selection (shown to perform better than AIC when selecting between DAISIE models18). However, like M26, all three alternative models preferred in some of the sensitivity analyses (M10, M11 and M22) are two-rate models under which bats have a higher rate of colonization than non-volant mammals and differ from non-volant mammals in one or more parameters.Preferred modelThe preferred model is the M26 model. Under this model, the background rate of cladogenesis for Malagasy mammals is 0.33 (0.27–0.36) events per lineage Myr−1 and the rate of anagenesis is 1.47 (1.18–2.12) events per lineage Myr−1 (Table S3, S4). The model is diversity independent (for both bats and non-volant species), meaning that there is no carrying capacity per clade (K per clade is infinite). The rate of natural extinction for non-volant mammals is 0.29 (0.22–0.31) events per lineage Myr−1, and for bats it is 0.46 (0.40–0.50) species per lineage Myr−1. The rate of colonization for non-volant mammals is 0.00036 (0.00027–0.00038) events per mainland species Myr−1, equivalent to 0.28 colonizations per Myr (0.21–0.30). The rate of colonization for bats is much higher, at 0.034 (0.030–0.038) events per mainland lineage Myr−1, equivalent to 7.5 successful bat colonization events per Myr (6.6–8.4).Evolutionary return timesThe island evolutionary return time (ERT) metric estimates the time it would take for an insular community to reach a given species diversity level assuming a given model of macroevolution with certain rates of colonization, speciation and natural extinction22. To estimate ERT, we first counted the number of species that were present on Madagascar in the late Holocene (pre-human diversity) (Table S5). We also counted the number of mammal species estimated to remain extant on the island if currently threatened species go extinct. Threatened species we classified as those that fall under the IUCN categories Vulnerable (VU), Endangered (EN) or Critically Endangered (CR). We estimated the ERT of Malagasy mammals for the following scenarios: (1) the return from current diversity to pre-human diversity, and (2) the return from the diversity that will remain if threatened species (VU + EN + CR) go extinct back to current diversity.There were some differences in ERTs between datasets (D1–D13). Most differences were subtle and some appear counterintuitive (Table S6). For example, a lower number of colonization events in the data (CS1 vs CS2) or a larger mainland pool both lead to a lower colonization rate, which evidently should increase the ERT. However, changes in the mainland pool and the number of colonization events can also lead to very small changes in the diversification rate (cladogenesis minus extinction), which has a much higher impact on the ERT than colonization rate. The largest differences in ERT were between the high and low human impact scenarios, because not only do estimated rates differ, but also the start and target diversities.Change in ERT between 2010 and 2021To compare how ERT for threatened species has been changing through time as human impact and working knowledge increases, we repeated analysis (2) using the IUCN threat statuses from 2010, 2015 and 2021. As the threat status of some species increases, they are uplisted by the IUCN, e.g. from Near Threatened (NT) to Vulnerable (VU) (e.g. the Madagascar rousette (Rousettus madagascariensis)), i.e. becoming threatened under IUCN classification of threat (VU, EN, CR). Although there have been many changes in status between the categories that we considered as threatened (e.g. from VU to EN, or from EN to CR), these are not considered in our analyses, as we only consider changes from non-threatened to threatened. Therefore, our analyses can be considered conservative; hence, the increase in ERT from 2010 to 2021 is probably higher than what we estimate. It is important to keep in mind that for different reasons, which include more intensively studied Malagasy mammalian groups and biopolitics, the manner certain data are weighed during assessments and resulting statutes are not necessarily the same across mammal groups.A total of 72 species were uplisted by the IUCN from a non-threat to a threat category between 2010 and 2021 (Supplementary Data S1). Of these, 57 moved from Not Evaluated (NE) or Data Deficient (DD) categories to one of the threat categories (VU, EN, CR), so may simply represent an increase in knowledge rather than a real increase in threat status. In addition, IUCN categories are not comparable across assessments for species that have undergone taxonomic revisions that may have altered their threat status—e.g. splitting of one species into two allopatric species will lead to a range size reduction. We therefore also repeated analyses only for those species that changed from an evaluated non-threat category (LC, NT) to a threat category (VU, EN, CR) and which have not undergone taxonomic changes between 2010 and 2021 or for which a taxonomic change was not the cause the up-listing. For all species that changed from a non-threat to a threat category between 2010 and 2021 (15 species, Supplementary Data S5) we consulted the literature to find out whether changes in taxonomy took place for that taxon, and whether those changes influenced the up-listing. We identified 10 species (Supplementary Data S5) for which there was no taxonomic change between 2010 and 2021 or for which a taxonomic change did not lead to an up-listing between those years. We then calculated ERTs for a scenario where we assume that only those 10 species were uplisted between 2010 and 2021 (Table S8, scenario C).Species diversity lost versus ERTThe number of species lost through anthropogenic extinctions and the ERT are two alternative ways of looking at the impact of humans on island biota. To assess whether the number of species lost is a good proxy for ERT, we ran simulations to measure how ERT varies with the number of extinct species to be recovered. For example, we compared how the time to return to pre-human diversity varies with starting diversity, e.g. assuming an increasing number of species have gone extinct. We ran simulations using the parameters of the best overall model for the main dataset. We first created 10,000 random start diversities, sampling between 0 species and a target species diversity, assuming variable numbers of species have gone extinct. For each of the starting diversities, we randomly specified a proportion of endemic and non-endemic species. Simulations were run in R.The shape of the island species diversity curve (how the total number of species on an island varies through time) under different DAISIE models can vary at different stages. For example, in an equilibrium model, diversity increases rapidly at early stages and low diversities, but at later stages it plateaus and increases slowly. In non-equilibrium models, diversity increases can for example be low at early stages and faster later. This will have implications for how ERT relates to the number of species that need to be recovered. We therefore separated the results into (a) returning to pre-human diversity (capturing a later stage of the diversity curve), (b) return to contemporary diversity (capturing an intermediate stage) and (c) return to half of the contemporary diversity (capturing early stage of the diversity curve). We did this separately for non-volant mammals and bats.Impact of increased knowledge on the ERTNew species discoveries and increasingly complete IUCN Red List assessments are likely to affect ERT estimates in the future. The discovery of new extant species may lead to an increase in ERT because undiscovered species are more likely to already be threatened (e.g. due to small range and population sizes). Taxonomic revisions may lead to species splits, resulting in additional threatened species. The known fossil record may also include cryptic species that cannot be identified using molecular methods if DNA is not available. The discovery of more taxa that have gone extinct since humans arrived will also likely increase the ERT (return to pre-human diversity). However, if new species are discovered, the rates of colonization and speciation estimated in DAISIE will also increase, and therefore ERTs may not rise dramatically.To assess how future species discoveries may affect our results, we performed analyses where we assume 30 new mammal species (15 bats species and 15 non-volant mammal species) will be discovered in the next 10 years on Madagascar. This is likely an overestimate. We simulated datasets by adding these bat and non-volant species at random locations to the main phylogenetic dataset D1, and repeated this procedure 1000 times. We fitted the preferred DAISIE model to these 1000 datasets and estimated the ERT for each of them. We then assumed that the newly discovered species were (a) all threatened, (b) half of them threatened and half already extinct or (c) all already extinct (since human arrival). The results of these analyses are summarized in Fig. S4. We found that the ERTs for non-volant mammals do not change substantially, increasing slightly or even declining under some scenarios. This is because although the number of species to recover increases, the estimated DAISIE rates also increase. On the other hand, under some scenarios, an increase in the number of extinct or extant bat species leads to large increases in ERT, as it takes longer on average to recover bat species according to the preferred DAISIE model.IUCN Red List assessments should become more comprehensive in the future; currently 8% of recognized Malagasy mammal species are Not Evaluated or assessed as Data Deficient, corresponding to 18 species (10 bats and 8 non-volant species). Thus, we also estimated how the completion of IUCN assessments may affect our results. If all species yet to be assessed by the IUCN were evaluated as threatened in the future, the ERT to return to contemporary diversity would rise to 26.2 (20.8–32) Myr for non-volant mammals (~13% increase) and 6.6 (5.5–7.8) Myr for bats (more than double) (Table S8, scenario A). The increase is proportionally higher for bats because it takes longer to recover bat species and because there are more unevaluated bat species.Reporting summaryFurther information on research design is available in the Nature Portfolio Reporting Summary linked to this article. More

Magnaval, J.-F., Glickman, L. T., Dorchies, P. & Morassin, B. Highlights of human toxocariasis. Korean J. Parasitol. 39, 1 (2001).Article 

Google Scholar 
Parise, M. E., Hotez, P. J. & Slutsker, L. Neglected parasitic infections in the United States: Needs and opportunities. Am. J. Trop. Med. Hyg. 90, 783–785 (2014).Article 

Google Scholar 
Holland, C. V. Knowledge gaps in the epidemiology of Toxocara: The enigma remains. Parasitology 144, 81–94 (2017).Article 

Google Scholar 
Richards, D. T. & Lewis, J. W. Fecundity and egg output by Toxocara canis in the red fox Vulpes vulpes. J. Helminthol. 75, 157–164 (2001).
Google Scholar 
Glickman, L. T. & Schantz, P. M. Epidemiology and pathogenesis of zoonotic toxocariasis. Epidemiol. Rev. 3, 230–250 (1981).Article 

Google Scholar 
Keegan, J. D. & Holland, C. V. A comparison of Toxocara canis embryonation under controlled conditions in soil and hair. J. Helminthol. 87, 78–84 (2013).Article 

Google Scholar 
Overgaauw, P. A. M. & Nederland, V. Aspects of toxocara epidemiology: Toxocarosis in dogs and cats. Crit. Rev. Microbiol. 23, 233–251 (1997).Article 

Google Scholar 
Strube, C., Heuer, L. & Janecek, E. Toxocara spp. infections in paratenic hosts. Vet. Parasitol. 193, 375–89 (2013).Parsons, J. C. Ascarid infections of cats and dogs. Vet. Clin. N. Am. Small Anim. Pract. 17, 1307–1339 (1987).Article 

Google Scholar 
Brunaská, M., Dubinský, P. & Reiterová, K. Toxocara canis: Ultrastructural aspects of larval moulting in the maturing eggs. Int. J. Parasitol. 25, 683–690 (1995).Article 

Google Scholar 
Dubinsky, P., Havasiova-Reiterova, K., Petko, B., Hovorka, I. & Tomasovicova, O. Role of small mammals in the epidemiology of toxocariasis. Parasitology 110(Pt 2), 187–193 (1995).Ma, G. et al. Human toxocariasis. Lancet Infect. Dis. 18, e14–e24 (2018).Article 

Google Scholar 
Despommier, D. Toxocariasis: clinical aspects, epidemiology, medical ecology, and molecular aspects. Clin. Microbiol. Rev. 16, 265–272 (2003).Article 

Google Scholar 
Hotez, P. J. & Wilkins, P. P. Toxocariasis: America’s Most Common Neglected Infection of Poverty and a Helminthiasis of Global Importance?. PLoS Negl. Trop. Dis. 3, e400 (2009).Article 

Google Scholar 
Fan, C.-K., Holland, C. V., Loxton, K. & Barghouth, U. Cerebral toxocariasis: Silent progression to neurodegenerative disorders?. Clin. Microbiol. Rev. 28, 663–686 (2015).Article 

Google Scholar 
Nathwani, D., Laing, R. B. S. & Currie, P. F. Covert toxocariasis—A cause of recurrent abdominal pain in childhood. Br. J. Clin. Pract. (1992).Rostami, A. et al. Seroprevalence estimates for toxocariasis in people worldwide: A systematic review and meta-analysis. PLoS Negl. Trop. Dis. 13, e0007809 (2019).Article 

Google Scholar 
Borecka, A. & Kłapeć, T. Epidemiology of human toxocariasis in Poland—A review of cases 1978–2009. Ann. Agric. Environ. Med. 22, 28–31 (2015).Article 

Google Scholar 
Krasnov, B. R., Mouillot, D., Shenbrot, G. I., Khokhlova, I. S. & Poulin, R. Geographical variation in host specificity of fleas (Siphonaptera) parasitic on small mammals: The influence of phylogeny and local environmental conditions. Ecography 27, 787–797 (2004).Article 

Google Scholar 
Andreassen, H. P. et al. Population cycles and outbreaks of small rodents: Ten essential questions we still need to solve. Oecologia 195, 601–622 (2021).Article 
ADS 

Google Scholar 
Ylönen, H. Vole cycles and antipredatory behaviour. Trends Ecol. Evol. https://doi.org/10.1016/0169-5347(94)90125-2 (1994).Article 

Google Scholar 
Hanski, I., Hansson, L. & Henttonen, H. Specialist predators, generalist predators, and the microtine rodent cycle. J. Anim. Ecol. https://doi.org/10.2307/5465 (1991).Article 

Google Scholar 
Martinez-Bakker, M. & Helm, B. The influence of biological rhythms on host–parasite interactions. Trends Ecol. Evol. 30, 314–326 (2015).Article 

Google Scholar 
Bajer, A., Pawelczyk, A., Behnke, J. M., Gilbert, F. S. & Sinski, E. Factors affecting the component community structure of haemoparasites in bank voles (Clethrionomys glareolus) from the Mazury Lake District region of Poland. Parasitology 122(Pt 1), 43–54 (2001).Article 

Google Scholar 
Behnke, J. M., Lewis, J. W., Zain, S. N. & Gilbert, F. S. Helminth infections in Apodemus sylvaticus in southern England: interactive effects of host age, sex and year on the prevalence and abundance of infections. J. Helminthol. 73, 31–44 (1999).Article 

Google Scholar 
Ferrari, N., Cattadori, I. M., Nespereira, J., Rizzoli, A. & Hudson, P. J. The role of host sex in parasite dynamics: Field experiments on the yellow-necked mouse Apodemus flavicollis. Ecol. Lett. 7, 88–94 (2003).Article 

Google Scholar 
Grzybek, M. et al. Long-term spatiotemporal stability and dynamic changes in helminth infracommunities of bank voles (Myodes glareolus) in NE Poland. Parasitology 142, 1722–1743 (2015).Article 

Google Scholar 
Reiterová, K. et al. Small rodents—permanent reservoirs of toxocarosis in different habitats of Slovakia. Helminthologia 50, (2013).Habig, B., Doellman, M. M., Woods, K., Olansen, J. & Archie, E. A. Social status and parasitism in male and female vertebrates: A meta-analysis. Sci. Rep. 8, 3629 (2018).Article 
ADS 

Google Scholar 
Izhar, R. & Ben-Ami, F. Host age modulates parasite infectivity, virulence and reproduction. J. Anim. Ecol. 84, 1018–1028 (2015).Article 

Google Scholar 
Migalska, M. et al. Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics. Mol. Ecol. 31, 3400–3415 (2022).Article 

Google Scholar 
Grzybek, M. et al. Zoonotic viruses in three species of voles from Poland. Animals 10, 1820 (2020).Article 

Google Scholar 
Bajer, A. et al. Rodents as intermediate hosts of cestode parasites of mammalian carnivores and birds of prey in Poland, with the first data on the life-cycle of Mesocestoides melesi. Parasit. Vectors 13, 95 (2020).Article 

Google Scholar 
Rabalski, L. et al. Zoonotic spillover of SARS-CoV-2: Mink-adapted virus in humans. bioRxiv (2021). https://doi.org/10.1101/2021.03.05.433713.Rabalski, L. et al. Severe acute respiratory syndrome coronavirus 2 in Farmed Mink (Neovison vison) Poland. Emerg. Infect. Dis. 27, 2333–2339 (2021).Article 

Google Scholar 
Grzybek, M. et al. Seroprevalence of Trichinella spp. infection in bank voles (Myodes glareolus)—A long term study. Int. J. Parasitol. Parasites Wildl. 9, 144–148 (2019).Binder, F. et al. Heterogeneous Puumala orthohantavirus situation in endemic regions in Germany in summer 2019. Transbound Emerg. Dis. 67, 502–509 (2020).Article 

Google Scholar 
Tołkacz, K. et al. Prevalence, genetic identity and vertical transmission of Babesia microti in three naturally infected species of vole, Microtus spp. (Cricetidae). Parasit. Vectors 10, 1–12 (2017).Tołkacz, K. et al. Bartonella infections in three species of Microtus: Prevalence and genetic diversity, vertical transmission and the effect of concurrent Babesia microti infection on its success. Parasit. Vectors 11, 491 (2018).Article 

Google Scholar 
Behnke, J. M. et al. Variation in the helminth community structure in bank voles (Clethrionomys glareolus) from three comparable localities in the Mazury Lake District region of Poland. Parasitology 123, 401–414 (2001).Article 

Google Scholar 
Behnke, J. M. et al. Temporal and between-site variation in helminth communities of bank voles (Myodes glareolus) from N.E. Poland. 2. The infracommunity level. Parasitology 135, 999–1018 (2008).Behnke, J. M. et al. Temporal and between-site variation in helminth communities of bank voles (Myodes glareolus) from N.E. Poland. 1. Regional fauna and component community levels. Parasitology 135, 985–997 (2008).Tołkacz, K. et al. Prevalence, genetic identity and vertical transmission of Babesia microti in three naturally infected species of vole, Microtus spp. (Cricetidae). Parasit. Vectors 10, 66 (2017).Morris, P. A review of mammalian age determination methods. Mamm. Rev. 2, 69–104 (1972).Antolová, D. et al. Small mammals: Paratenic hosts for species of Toxocara in eastern Slovakia. J. Helminthol. 87, 52–58 (2013).Article 

Google Scholar 
Reiterová, K. et al. Small rodents—permanent reservoirs of toxocarosis in different habitats of Slovakia. Helminthologia 50, 20–26 (2013).Article 

Google Scholar 
Reperant, L. A., Hegglin, D., Tanner, I., Fischer, C. & Deplazes, P. Rodents as shared indicators for zoonotic parasites of carnivores in urban environments. Parasitology 136, 329–337 (2009).Article 

Google Scholar 
Savigny, D. H. In vitro maintenance of Toxocara canis larvae and a simple method for the production of Toxocara ES antigen for use in serodiagnostic tests for visceral larva migrans. J. Parasitol. 61, 781–782 (1975).Article 

Google Scholar 
Cuéllar, C., Fenoy, S. & Guillén, J. L. Cross-reactions of sera from Toxascaris leonina and Ascaris suum infected mice with Toxocara canis, Toxascaris leonina and Ascaris suum antigens. Int. J. Parasitol. 25, 731–739 (1995).Article 

Google Scholar 
Naguleswaran, A., Hemphill, A., Rajapakse, R. P. V. J. & Sager, H. Elaboration of a crude antigen ELISA for serodiagnosis of caprine neosporosis: Validation of the test by detection of Neospora caninum-specific antibodies in goats from Sri Lanka. Vet. Parasitol. 126, 257–262 (2004).Article 

Google Scholar 
Sokal, R. R. & Rohlf, F. J. Statistical Tables (Freeman, 1995).MATH 

Google Scholar 
Behnke, J. M. et al. Variation in the helminth community structure in bank voles (Clethrionomys glareolus) from three comparable localities in the mazury lake istrict region of Poland. Parasitology 123, 401–414 (2001).Article 

Google Scholar 
Grzybek, M., Bajer, A., Behnke-Borowczyk, J., Al-Sarraf, M. & Behnke, J. M. Female host sex-biased parasitism with the rodent stomach nematode Mastophorus muris in wild bank voles (Myodes glareolus). Parasitol. Res. 114, 523–533 (2014).Article 

Google Scholar 
Grzybek, M. et al. Seroprevalence of TBEV in bank voles from Poland-a long-term approach. Emerg. Microbes Infect. 7, 145 (2018).Article 

Google Scholar 
Pullan, R. L., Smith, J. L., Jasrasaria, R. & Brooker, S. J. Global numbers of infection and disease burden of soil transmitted helminth infections in 2010. Parasit Vectors 7, 37 (2014).Article 

Google Scholar 
GBD 2017 DALYs and HALE Collaborators. Global, regional, and national disability-adjusted life-years (DALYs) for 359 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet 392, 1859–1922 (2018).Reperant, L. A., Hegglin, D., Tanner, I., Fisher, C. & Deplazes, P. Rodents as shared indicators for zoonotic parasites of carnivores in urban environments. Parasitology 136, 329–337 (2009).Article 

Google Scholar 
Antolová, D., Reiterová, K., Miterpáková, M., Stanko, M. & Dubinský, P. Circulation of Toxocara spp. in suburban and rural ecosystems in the Slovak Republic. Vet. Parasitol 126, 317–324 (2004).Reiterová, K. et al. Small rodents—permanent reservoirs of toxocarosis in different habitats of Slovakia. Helminthologia (Poland) 50, 20–26 (2013).Article 

Google Scholar 
Dvorožňáková, E., Kołodziej-Sobocińska, M., Hurníková, Z., Víchová, B. & Zub, K. Prevalence of zoonotic pathogens in wild rodents living in the Białowieża Primeval Forest, Poland. Ann. Parasitol. 62 (2016).Dubinský, P., Havasiova-Reiterova, K., Peťko, B., Hovorka, I. & Tomašovičová, O. Role of small mammals in the epidemiology of toxocariasis. Parasitology 110, 187–193 (1995).Article 

Google Scholar 
Antolová, D., Reiterová, K., Miterpáková, M., Stanko, M. & Dubinský, P. Circulation of Toxocara spp. in suburban and rural ecosystems in the Slovak Republic. Vet. Parasitol. 126, 317–24 (2004).Hildebrand, J., Zalesny, G., Okulewicz, A. & Baszkiewicz, K. Preliminary studies on the zoonotic importance of rodents as a reservoir of toxocariasis from recreation grounds in Wroclaw (Poland). Helminthologia 46, 80–84 (2009).Article 

Google Scholar 
Dvorožňáková, E., Kołodziej-Sobocińska, M., Hurníková, Z., Víchová, B. & Zub, K. Prevalence of zoonotic pathogens in wild rodents living in the Białowieża Primeval Forest Poland. Ann. Parasitol. 62, 183 (2016).
Google Scholar 
Azam, D., Ukpai, O. M., Said, A., Abd-Allah, G. A. & Morgan, E. R. Temperature and the development and survival of infective Toxocara canis larvae. Parasitol. Res. 110, 649–656 (2012).Article 

Google Scholar 
Kloch, A., Bednarska, M. & Bajer, A. Intestinal macro- and microparasites of wolves (Canis lupus L.) from north-eastern Poland recovered by coprological study. Ann. Agric. Environ. Med. 12, 237–45 (2005).Mierzejewska, E. J. et al. The efficiency of live-capture traps for the study of Red Fox (Vulpes vulpes) cubs: A three-year study in Poland. Animals 10, 374 (2020).Article 

Google Scholar 
Karamon, J. et al. Intestinal helminths of raccoon dogs (Nyctereutes procyonoides) and red foxes (Vulpes vulpes) from the Augustów Primeval Forest (north-eastern Poland). J. Vet. Res. (Poland) 60, 273–277 (2016).Article 

Google Scholar 
Cisek, A., Ramisz, A., Balicka-Ramisz, A., Pilarczyk, B. & Laurans, L. The prevalence of Toxocara canis (Werner, 1782) in dogs and red foxes in north-west Poland. Wiad Parazytol 50, 641–646 (2004).
Google Scholar 
Jarošová, J., Antolová, D., Lukáč, B. & Maďari, A. A Survey of Intestinal Helminths of dogs in Slovakia with an emphasis on zoonotic species. Animals 11, 3000 (2021).Article 

Google Scholar 
Kidawa, D. & Kowalczyk, R. The effects of sex, age, season and habitat on diet of the red fox Vulpes vulpes in northeastern Poland. Acta Theriol. (Warsz) 56, 209–218 (2011).Article 

Google Scholar 
Grzybek, M. et al. Seroprevalence of Trichinella spp. infection in bank voles (Myodes glareolus)—A long term study. Int. J. Parasitol. Parasites Wildl. (2019). https://doi.org/10.1016/j.ijppaw.2019.03.005.Grzybek, M. et al. Seroprevalence of Toxoplasma gondii among Sylvatic Rodents in Poland. Animals 11, 1048 (2021).Article 

Google Scholar 
Maciag, L., Morgan, E. R. & Holland, C. Toxocara: Time to let cati ‘out of the bag’. Trends Parasitol. 38, 280–289 (2022).Article 

Google Scholar 
Foreman-Worsley, R., Finka, L. R., Ward, S. J. & Farnworth, M. J. Indoors or outdoors? an international exploration of owner demographics and decision making associated with lifestyle of pet cats. Animals 11, 253 (2021).Article 

Google Scholar 
Liberg, O. Food habits and prey impact by feral and house-based domestic cats in a rural area in southern Sweden. J. Mammal. 65, 424–432 (1984).Article 

Google Scholar 
Krücken, J. et al. Small rodents as paratenic or intermediate hosts of carnivore parasites in Berlin Germany. PLoS ONE 12, e0172829 (2017).Article 

Google Scholar 
Dunsmore, J. D., Thompson, R. C. A. & Bates, I. A. The accumulation of Toxocara canis larvae in the brains of mice. Int. J. Parasitol. 13, 517–521 (1983).Article 

Google Scholar 
Nijsse, R., Mughini-Gras, L., Wagenaar, J. A., Franssen, F. & Ploeger, H. W. Environmental contamination with Toxocara eggs: a quantitative approach to estimate the relative contributions of dogs, cats and foxes, and to assess the efficacy of advised interventions in dogs. Parasit. Vectors 8, 397 (2015).Article 

Google Scholar 
Dunsmore, J. D., Thompson, R. C. A. & Bates, I. A. Prevalence and survival of Toxocara canis eggs in the urban environment of Perth Australia. Vet. Parasitol. 16, 303–311 (1984).Article 

Google Scholar 
Duscher, G. G., Leschnik, M., Fuehrer, H.-P. & Joachim, A. Wildlife reservoirs for vector-borne canine, feline and zoonotic infections in Austria. Int. J. Parasitol. Parasites Wildl. 4, 88–96 (2015).Article 

Google Scholar 
Ghai, R. R. et al. A generalizable one health framework for the control of zoonotic diseases. Sci. Rep. 12, 8588 (2022).Article 
ADS 

Google Scholar 
Grzybek, M. et al. Zoonotic virus seroprevalence among bank voles, Poland, 2002–2010. Emerg. Infect. Dis. 25, 1607–1609 (2019).Article 

Google Scholar  More

One Health Institute, School of Veterinary Medicine, University of California, Davis, Davis, CA, 95616, USAPranav S. Pandit, Tracey Goldstein, Megan M. Doyle, Nicole R. Gardner, Brian Bird, Woutrina Smith, David Wolking, Kirsten Gilardi, Corina Monagin, Terra Kelly, Marcela M. Uhart, Lucy Keatts, Jonna A. K. Mazet & Christine K. JohnsonCenter for Infection and Immunity, Columbia University, New York, NY, 10032, USASimon J. AnthonyEcoHealth Alliance, 520 Eighth Avenue, New York, NY, 10018, USAKevin J. Olival, Jonathan H. Epstein, Catherine Machalaba, Melinda K. Rostal, Patrick Dawson, Emily Hagan, Ava Sullivan, Hongying Li, Aleksei A. Chmura, Alice Latinne, Ariful Islam, James Desmond, Tom Hughes, William B. Karesh & Peter DaszakLabyrinth Global Health, Inc., 546 15th Ave NE, St Petersburg, FL, 33704, USAChristian Lange, Tammie O’Rourke & Karen SaylorsWildlife Conservation Society, Health Program, Bronx, NY, USASarah Olson, A. Patricia Mendoza, Cátia Dejuste de Paula, Amanda Fine & Cátia Dejuste de PaulaWildlife Conservation Society (WCS), Peru Program, Lima, PeruA. Patricia Mendoza & Alberto PerezGlobal Health Program, Smithsonian’s National Zoological Park and Conservation Biology Institute, Washington, DC, USADawn Zimmerman, Marc Valitutto & Ohnmar AungMosaic/Global Viral Cameroon, Yaoundé, CameroonMatthew LeBreton, Moctar Mouiche & Suzan MurrayMetabiota Inc, Nanaimo, VC, CanadaDavid McIver & Soubanh SilithammavongInstitut Pasteur du Cambodge, 5 Monivong Blvd, PO Box 983, Phnom Penh, 12201, CambodiaVeasna DuongWuhan Institute of Virology, Chinese Academy of Sciences, Wuhan, ChinaZhengli ShiKinshasa School of Public Health, University of Kinshasa, Kinshasa, Democratic Republic of the CongoPrime MulembakaniMetabiota Inc., Kinshasa, Democratic Republic of the CongoCharles KumakambaEgypt National Research Centre, 12311, Dokki, Giza, EgyptMohamed AliAklilu Lemma Institute of Pathobiology, Addis Ababa University, Addis Ababa, EthiopiaNigatu KebedeMetabiota Cameroon Ltd, Yaoundé, Centre Region Avenue Mvog-Fouda Ada, Av 1.085, Carrefour Intendance, Yaoundé, BP 15939, CameroonUbald TamoufeMilitary Veterinarian (Rtd.), P.O. Box CT2585, Accra, GhanaSamuel Bel-NonoCentre de Recherche en Virologie (VRV) Projet Fievres Hemoraquiques en Guinée, BP 5680, Nongo/Contéya-Commune de Ratoma, GuineaAlpha CamaraPrimate Research Center, Bogor Agricultural University, Bogor, 16151, IndonesiaJoko PamungkasFaculty of Veterinary Medicine, Bogor Agricultural University, Darmaga Campus, Bogor, 16680, IndonesiaJoko PamungkasDepartment Environment and Health, Institut Pasteur de Côte d’Ivoire, PO BOX 490, Abidjan 01, Ivory CoastKalpy J. CoulibalyDepartment of Basic Medical Veterinary Sciences, College of Veterinary Medicine, Jordan University of Science and Technology, Ar-Ramtha, JordanEhab Abu-BashaMolecular Biology Laboratory, Institute of Primate Research, Nairobi, KenyaJoseph KamauDepartment of Biochemistry, University of Nairobi, Nairobi, KenyaJoseph KamauConservation Medicine, Sungai Buloh, Selangor, MalaysiaTom HughesWildlife Conservation Society (WCS), Mongolia Program, Ulaanbaatar, MongoliaEnkhtuvshin ShiilegdambaCenter for Molecular Dynamics Nepal (CMDN), Thapathali -11, Kathmandu, NepalDibesh KarmacharyaRegional Headquarters, Mountain Gorilla Veterinary Project, Musanze, RwandaJulius Nziza & Benard SsebideUniversité Cheikh Anta Diop, BP 5005, Dakar, SénégalDaouda NdiayeMetabiota, Inc. Sierra Leone, Freetown, Sierra LeoneAiah GbakimaDepartment of Veterinary Medicine and Public Health, College of Veterinary Medicine and Biomedical Sciences, Sokoine University of Agriculture, Morogoro, TanzaniaZikankuba sajaliThai Red Cross Emerging Infectious Diseases Clinical Center, King Chulalongkorn Memorial Hospital, Bangkok, ThailandSupaporn WacharapluesadeeWildlife Conservation Society (WCS), Bolivia Program, La Paz, BoliviaErika Alandia RoblesFacultad de Medicina Veterinaria y Zootecnia, Universidad Nacional Autónoma de México, México City, 04510, MexicoGerardo SuzánCentro de Biodiversidad y Genética, Universidad Mayor de San Simón, Cochabamba, BoliviaLuis F. AguirreLaboratório de Epidemiologia e Geoprocessamento (EpiGeo), Instituto de Medicina Veterinária (IMV) Universidade Federal do Pará (UFPA), BR-316 Km 31, Castanhal, Pará, 69746-360, BrazilMonica R. SolorioDepartment of Microbiology, Sanjay Gandhi Post Graduate Institute of Medical Sciences, Lucknow, Uttar Pradesh, IndiaTapan N. DholeWildlife Conservation Society (WCS), Vietnam Program, Hanoi, VietnamNguyen T. T. NgaMelbourne Veterinary School, Faculty of Veterinary and Agricultural Sciences, University of Melbourne, Werribee, VIC, 3030, AustraliaPeta L. HitchensNyati Health Consulting, 2175 Dodds Road, Nanaimo, BC, V9X0A4, CanadaDamien O. Joly More

Durden, L. Lice (Phthiraptera). In Medical and Veterinary Entomology 3rd edn (eds Mullen, G. & Durden, L.) 79–106 (Academic Press, 2019).Chapter 

Google Scholar 
Stork, N. E. & Lyal, C. H. C. Extinction or ‘co-extinction’ rates?. Nature 366, 307. https://doi.org/10.1038/366307a0 (1993).Article 
ADS 

Google Scholar 
Koh, L. P. et al. Species coextinctions and the biodiversity crisis. Science 305, 1632–1634. https://doi.org/10.1126/science.1101101 (2004).Article 
ADS 
CAS 

Google Scholar 
Gerlach J (2014) Haematopinus oliveri. The IUCN Red List of Threatened Species 2014: e.T9621A21423551. https://doi.org/10.2305/IUCN.UK.2014-1.RLTS.T9621A21423551.en.Mingozzi, T. & Stève, R. Analysis of a historical extirpation of the bearded vulture Gypaetus barbatus (L.) in the Western Alps (France-Italy): former distribution and causes of extirpation. Biol. Conserv. 79, 155–171. https://doi.org/10.1016/S0006-3207(96)00110-3 (1997).Article 

Google Scholar 
Schaub, M., Zink, R., Beissmann, H., Sarrazin, F. & Arlettaz, R. When to end releases in reintroduction programmes: demographic rates and population viability analysis of bearded vultures in the Alps. J. Appl. Ecol. 46, 92–100. https://doi.org/10.1111/j.1365-2664.2008.01585.x (2009).Article 

Google Scholar 
BirdLife International. Gypaetus barbatus (amended version of 2017 assessment). The IUCN Red List of Threatened Species 2017: e.T22695174A118590506. https://doi.org/10.2305/IUCN.UK.2017-3.RLTS.T22695174A118590506.en, accessed 07 Apr 2021 (2017).Price, R. D., Hellenthal, R. A., Palma, R. L., Johnson, K. P., & Clayton, D. H. The chewing lice: World checklist and biological overview. Illinois Natural History Survey Special Publication 24. Illinois (2003).Clay, T. Revisions of mallophaga genera. Degeeriella from the Falconiformes. Bull. Br. Mus. (Nat. Hist.) 7, 123–207 (1958).
Google Scholar 
Martín Mateo, M. P. Fauna Ibérica, Vol. 32. Phthiraptera, Ischnocera. Museo Nacional de Ciencias Naturales (CSIC), Madrid (2009).Hoberg, E. P., Brooks, D. R. & Siegel-Causey, D. Host-parasite co-speciation: history, principles, and prospects. In Host-Parasite Evolution: General Principles and Avian Models (eds Clayton, D. H. & Moore, J.) 212–235 (Oxford University Press, 1997).
Google Scholar 
Johnson, K. P., Weckstein, J. D., Witt, C. C., Faucett, R. C. & Moyle, R. G. The perils of using host relationships in parasite taxonomy: phylogeny of the Degeeriella complex. Mol. Phylogenet. Evol. 23, 150–157. https://doi.org/10.1016/S1055-7903(02)00014-3 (2002).Article 
CAS 

Google Scholar 
Catanach, T. A. & Johnson, K. P. Independent origins of the feather lice (Insecta: Degeeriella) of raptors. Biol. J. Linn. Soc. 114, 837–847. https://doi.org/10.1111/bij.12453 (2015).Article 

Google Scholar 
Pérez, J. M., Ruiz-Martínez, I. & Cooper, J. E. Occurrence of chewing lice on Spanish raptors. Ardeola 43, 129–138 (1996).
Google Scholar 
Ash, J. S. A study of the mallophagan of birds with particular reference to their ecology. Ibis 102, 93–110. https://doi.org/10.1111/j.1474-919X.1960.tb05095.x (1960).Article 

Google Scholar 
Askew, R. R. Parasitic Insects (Heinemann Educational, 1971).
Google Scholar 
Marshall, A. G. The Ecology of Parasitic Insects (Academic Press, 1981).
Google Scholar 
Bartlow, A. W., Villa, S. M., Thompson, M. W. & Bush, S. E. Walk or ride? Phoretic behaviour of amblyceran and ischnoceran lice. Int. J. Parasitol. 46, 221–227. https://doi.org/10.1016/j.ijpara.2016.01.003 (2016).Article 

Google Scholar 
Leonardi, M. S., Crespo, E. A., Raga, J. A. & Fernández, M. Scanning electron microscopy of Antarctophthirus microchir (Phthiraptera: Anoplura: Echinophthiriidae): Studying morphological adaptations to aquatic life. Micron 43, 929–936. https://doi.org/10.1016/j.micron.2012.03.009 (2012).Article 

Google Scholar 
Ortega Insaurralde, I., Minoli, S., Toloza, A. C., Picollo, M. I. & Barrozo, R. B. The sensory machinery of the head louse Pediculus humanus capitis: from the antennae to the brain. Front. Physiol. 10, 434. https://doi.org/10.3389/fphys.2019.00434 (2019).Article 

Google Scholar 
Ortega Insaurralde, I., Picollo, M. I. & Barrozo, R. B. Sensory features of the human louse antenna: New contributions and comparisons between ecotypes. Med. Vet. Entomol. 35, 219–224. https://doi.org/10.1111/mve.12485 (2021).Article 
CAS 

Google Scholar 
Page, R. D. M., Lee, P. L. M., Becher, S. A., Griffiths, R. & Clayton, D. H. A different tempo of mitochondrial DNA evolution in birds and their parasitic lice. Mol. Phylogenet. Evol. 9, 276–293. https://doi.org/10.1006/mpev.1997.0458 (1998).Article 
CAS 

Google Scholar 
Cruickshank, R. H. et al. Phylogenetic analysis of partial sequences of elongation factor 1α identifies major groups of lice (Insecta: Phthiraptera). Mol. Phylogenet. Evol. 19, 202–215. https://doi.org/10.1006/mpev.2001.0928 (2001).Article 
CAS 

Google Scholar 
Murrell, A. & Barker, S. C. Multiple origins of parasitism in lice: phylogenetic analysis of SSU rDNA indicates that the Phthiraptera and Psocoptera are not monophyletic. Parasitol. Res. 97, 274–280. https://doi.org/10.1007/s00436-005-1413-8 (2005).Article 

Google Scholar 
Whiteman, N. K., Kimball, R. T. & Parker, P. G. Co-phylogeography and comparative population genetics of the threatened Galápagos hawk and three ectoparasite species: ecology shapes population histories within parasite communities. Mol. Ecol. 16, 4759–4773. https://doi.org/10.1111/j.1365-294X.2007.03512.x (2007).Article 
CAS 

Google Scholar 
Palma, R. L. Slide-mounting of Lice: a detailed description of the Canada Balsam technique. N. Z. Entomol. 6, 432–436. https://doi.org/10.1080/00779962.1978.9722313 (1978).Article 

Google Scholar 
Soler-Cruz, M. D. & Martín-Mateo, M. P. Scanning electron microscopy of legs of two species of sucking lice (Anoplura: Phthiraptera). Micron 40, 401–408. https://doi.org/10.1016/j.micron.2008.10.001 (2009).Article 
CAS 

Google Scholar 
Hafner, M. S. et al. Disparate rates of molecular evolution in cospeciating hosts and parasites. Science 265, 1087–1090. https://doi.org/10.1126/science.8066445 (1994).Article 
ADS 
CAS 

Google Scholar 
Simon, C. et al. Evolution, weighting and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann. Entomol. Soc. Am. 87, 651–701. https://doi.org/10.1093/aesa/87.6.651 (1994).Article 
CAS 

Google Scholar 
Danforth, B. N. & Ji, S. Elongation factor-1α occurs as two copies in bees: Implications for phylogenetic analysis of EF-1α sequences in insects. Mol. Biol. Evol. 15, 225–235. https://doi.org/10.1093/oxfordjournals.molbev.a025920 (1998).Article 
CAS 

Google Scholar 
Smith, V. S., Page, R. D. M. & Johnson, K. P. Data incongruence and the problem of avian louse phylogeny. Zool. Scr. 33, 239–259. https://doi.org/10.1111/j.0300-3256.2004.00149.x (2004).Article 

Google Scholar 
Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. J. Mol. Biol. 215, 403–410. https://doi.org/10.1016/S0022-2836(05)80360-2 (1990).Article 
CAS 

Google Scholar 
Darriba, D., Taboada, G. L., Doallo, R. & Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 9, 772. https://doi.org/10.1038/nmeth.2109 (2012).Article 
CAS 

Google Scholar 
Huelsenbeck, J. P. & Ronquist, F. MRBAYES: Bayesian inference of phylogeny. Bioinformatics 17, 754–755. https://doi.org/10.1093/bioinformatics/17.8.754 (2001).Article 
CAS 

Google Scholar 
Zwickl, D. J. Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. Ph.D. Thesis Dissertation, The University of Texas at Austin, Texas (2006).Rambaut, A. FigTree v1.4.2. Institute of Evolutionary Biology (University of Edinburgh, 2014).
Google Scholar 
Brown, C. J. Plumages and measurements of the Bearded Vulture in Southern Africa. Ostrich 60, 165–171 (1989).Article 

Google Scholar 
Chatterjee, P., Payra, A. & Sen, S. Insecta: Phthiraptera. In Faunal Diversity of Indian Himalaya (eds Chandra, K. et al.) 297–304 (Zoological Survey of India, 2018).
Google Scholar 
Liébanas, G. et al. The morphology of Colpocephalum pectinatum (Phthiraptera: Amblycera: Menoponidae) under scanning electron microscopy. Arthropod Struct. Dev. 64, 101085. https://doi.org/10.1016/j.asd.2021.101085 (2021).Article 

Google Scholar 
Pérez, J. M. Sobre algunos aspectos de la parasitación por malógafos en aves de presa. Ph.D. thesis Dissertation, Granada University (1990).Arya, G., Ahmad, A., Bansal, N., Saxena, R. & Saxena, A. K. Nature of placodean sensilla of four ischnoceran Phthiraptera. Entomon 35, 199–202 (2010).
Google Scholar 
Khan, V., Bansal, N., Arya, G., Ahmad, A. & Saxena, A. K. Contribution to the morphology of Degeeriella regalis (Insecta, Phthiraptera, Ischnocera). J. Entomol. Res. 35, 93–96 (2011).
Google Scholar 
Agarwal, G. P. et al. Bio-ecology of the louse, Upupicola upupae, infesting the Common Hoopoe, Upupa epops. J. Insect. Sci. 11, 46. https://doi.org/10.1673/031.011.4601 (2011).Article 
CAS 

Google Scholar 
Singh, P., Gupta, N., Khan, G., Kumar, S. & Ahmad, A. Diagnostic characters of three nymphal instars and morphological features of adult Collard-dove louse Columbicola bacillus (Phthiraptera: Insecta). J. Appl. Nat. Sci. 11, 7–11. https://doi.org/10.31018/jans.v11i1.1855 (2019).Article 

Google Scholar 
Clayton, D. H. & Johnson, K. P. Linking coevolutionary history to ecological process: Doves and lice. Evolution 57, 2335–2341. https://doi.org/10.1111/j.0014-3820.2003.tb00245.x (2003).Article 

Google Scholar 
Barker, S. C. Lice, cospeciation and parasitism. Int J Parasitol 26, 219–222 (1996).Article 
CAS 

Google Scholar 
Page, R. D. M., Clayton, D. H. & Paterson, A. A. Lice and cospeciation: A response to Barker. Int. J. Parasitol. 26, 213–218. https://doi.org/10.1016/0020-7519(95)00115-8 (1996).Article 
CAS 

Google Scholar 
Paterson, A. M. & Gray, R. D. Host-parasite cospeciation, host-switching and missing the boat. In Host-Parasite Evolution: General Principles and Avian Models (eds Clayton, D. H. & Moore, J.) 236–250 (Oxford University Press, 1997).
Google Scholar 
Paterson, A. M., Palma, R. L. & Gray, R. D. How frequently do avian lice miss the boat? Implications for coevolutionary studies. Syst. Biol. 48, 214–223. https://doi.org/10.1080/106351599260544 (1999).Article 

Google Scholar 
Frey, H. & Walter, W. The reintroduction of the bearded vulture Gypaetus barbatus into the Alps. In Raptors in the Modern World (eds Meyburg, B. U. & Chancellor, R. D.) 341–344 (WWGBP, 1989).
Google Scholar 
Pérez, J. M., Sánchez, I. & Palma, R. L. The dilemma of conserving parasites: the case of Felicola (Lorisicola) isidoroi (Phthiraptera: Trichodectidae) and its host, the endangered Iberian lynx (Lynx pardinus). Insect. Conserv. Divers. 6, 680–686. https://doi.org/10.1111/icad.12021 (2013).Article 

Google Scholar  More