More stories

  • in

    Tree mode of death and mortality risk factors across Amazon forests

    School of Geography, Earth and Enviornmental Sciences, University of Birmingham, Birmingham, UK
    Adriane Esquivel-Muelbert & Thomas A. M. Pugh

    School of Geography, University of Leeds, Leeds, UK
    Adriane Esquivel-Muelbert, Oliver L. Phillips, Roel J. W. Brienen, Martin J. P. Sullivan, Timothy R. Baker, Emanuel Gloor, Aurora Levesley, Simon L. Lewis, Karina Liana Lisboa Melgaço Ladvocat, Gabriela Lopez-Gonzalez, Nadir Pallqui Camacho, Julie Peacock, Georgia Pickavance & David Galbraith

    Birmingham Institute of Forest Research, University of Birmingham, Birmingham, UK
    Adriane Esquivel-Muelbert & Thomas A. M. Pugh

    School of Geography, Earth and Environmental Sciences, University of Plymouth, Plymouth, UK
    Sophie Fauset

    Department of Natural Sciences, Manchester Metropolitan University, Manchester, UK
    Martin J. P. Sullivan

    International Master Program of Agriculture, National Chung Hsing University, Taichung, Taiwan
    Kuo-Jung Chao

    Geography, College of Life and Environmental Sciences, University of Exeter, Exeter, UK
    Ted R. Feldpausch

    Instituto Nacional de Pesquisas da Amazônia, Manaus, Brazil
    Niro Higuchi, Adriano José Nogueira Lima & Carlos Quesada

    School of Mathematics, University of Leeds, Leeds, UK
    Jeanne Houwing-Duistermaat & Haiyan Liu

    Faculty of Natural Sciences, Department of Life, Imperial College London Sciences, London, UK
    Jon Lloyd

    Environmental Change Institute, School of Geography and the Environment, University of Oxford, Oxford, UK
    Yadvinder Malhi & Simone Matias de Almeida Reis

    UNEMAT – Universidade do Estado de Mato Grosso PPG-Ecologia e Conservação, Campus de Nova Xavantina, Nova Xavantina, MT, Brazil
    Beatriz Marimon, Ben Hur Marimon Junior, Paulo Morandi, Edmar Almeida de Oliveira & Simone Matias de Almeida Reis

    Jardín Botánico de Missouri, Oxapampa, Peru
    Abel Monteagudo-Mendoza, Victor Chama Moscoso, Luis Valenzuela Gamarra & Rodolfo Vasquez Martinez

    Forest Ecology and Forest Management Group, Wageningen University and Research, Wageningen, Netherlands
    Lourens Poorter, Frans Bongers, Marielos Peña-Claros & Pieter Zuidema

    Centro de Ciências Biológicas e da Natureza, Universidade Federal do Acre, Rio Branco, AC, Brazil
    Marcos Silveira

    Instituto de Investigaciones para el Desarrollo Forestal (INDEFOR), Universidad de Los Andes, Mérida, Venezuela
    Emilio Vilanova Torre & Julio Serrano

    University of California, Berkeley, CA, USA
    Emilio Vilanova Torre

    Escuela de Ciencias Agropecuarias y Ambientales, Universidad Nacional Abierta y a Distancia, Boyacá, Colombia
    Esteban Alvarez Dávila

    Fundación ConVida, Medellín, Colombia
    Esteban Alvarez Dávila

    Instituto de Investigaciones de la Amazonia Peruana, Iquitos, Peru
    Jhon del Aguila Pasquel, Gerardo A. Aymard C., Nallaret Davila Cardozo & Eurídice Honorio Coronado

    Instituto de Biodiversidade e Florestas, Universidade Federal do Oeste do Pará, Santarém, Brazil
    Everton Almeida

    Center for Tropical Conservation, Nicholas School of the Environment, University in Durham, Durham, NC, USA
    Patricia Alvarez Loayza

    Projeto Dinâmica Biológica de Fragmentos, Instituto Nacional de Pesquisas da Amazônia Florestais, Manaus, AM, Brazil
    Ana Andrade & José Luís Camargo

    National Institute for Space Research (INPE), São José dos Campos, SP, Brazil
    Luiz E. O. C. Aragão

    Museo de Historia Natural Noel Kempff Mercado, Universidad Autónoma Gabriel Rene Moreno, Santa Cruz de la Sierra, Bolivia
    Alejandro Araujo-Murakami & Marisol Toledo

    Wageningen Environmental Research, Wageningen University and Research, Wageningen, Netherlands
    Eric Arets

    Dirección de la Carrera de Biología, Universidad Autónoma Gabriel René Moreno, Santa Cruz de la Sierra, Bolivia
    Luzmila Arroyo

    INRAE, UMR EcoFoG, CNRS, Cirad, AgroParisTech, Université des Antilles, Université de Guyane, Kourou, France
    Michel Baisie, Damien Bonal, Benoit Burban, Aurélie Dourdain, Maxime Rejou-Machain & Clement Stahl

    Department of Biological Sciences, International Center for Tropical Botany, Florida International University, Miami, FL, USA
    Christopher Baraloto

    Centro de Energia Nuclear na Agricultura, Universidade de São Paulo, Piracicaba, Brazil
    Plínio Barbosa Camargo

    Universidade Federal do Acre, Campus Floresta, Cruzeiro do Sul, Brazil
    Jorcely Barroso

    UR Forest & Societies, CIRAD, Montpellier, France
    Lilian Blanc

    Department of Biology, Utrecht, Netherlands
    René Boot

    Woods Hole Research Center, Falmouth, MA, USA
    Foster Brown

    Laboratório de Botânica e Ecologia Vegetal, Universidade Federal do Acre, Rio Branco, AC, Brazil
    Wendeson Castro

    Laboratoire Evolution et Diversite Biologique, CNRS, Toulouse, France
    Jerome Chave

    Inventory and Monitoring Program, National Park Service, Fort Collins, CO, USA
    James Comiskey

    Proyecto Castaña, Madre de Dios, Peru
    Fernando Cornejo Valverde

    Instituto de Geociências, Faculdade de Meteorologia, Universidade Federal do Para, Belém, Brazil
    Antonio Lola da Costa

    Department of Anthropology and Primate Molecular Ecology and Evolution Laboratory, University of Texas, Austin, TX, USA
    Anthony Di Fiore

    National Museum of Natural History, Smithsonian Institute, Washington, DC, USA
    Terry Erwin

    Universidad Nacional Jorge Basadre de Grohmann, Tacna, Peru
    Gerardo Flores Llampazo

    Museu Paraense Emílio Goeldi, Belém, Brazil
    Ima Célia Guimarães Vieira & Rafael Salomão

    Instituto Venezolano de Investigaciones Científicas (IVIC), Caracas, Venezuela
    Rafael Herrera

    IIAMA, Universitat Politécnica de València, València, Spain
    Rafael Herrera

    Universidad Nacional de San Antonio Abad del Cusco, Cusco, Peru
    Isau Huamantupa-Chuquimaco

    Instituto Amazónico de Investigaciones Imani, Universidad Nacional de Colombia Sede Amazonia, Leticia, Colombia
    Eliana Jimenez-Rojas

    Agteca, Santa Cruz, Bolivia
    Timothy Killeen

    College of Science and Engineering, James Cook University, Cairns, QLD, Australia
    Susan Laurance & William Laurance

    Department of Geography, University College London, London, UK
    Simon L. Lewis

    Environmental Science and Policy, George Mason University, Fairfax, VA, USA
    Thomas Lovejoy

    Research School of Biology, Australian National University, Canberra, ACT, Australia
    Patrick Meir

    School of Geosciences, University of Edinburgh, Edinburgh, UK
    Patrick Meir

    Escuela de Ciencias Forestales, Unidad Académica del Trópico, Universidad Mayor de San Simón, Cochabamba, Bolivia
    Casimiro Mendoza

    Facultad de Ingeniería Ambiental, Universidad Estatal Amazónica, Puyo, Ecuador
    David Neill

    Universidad Nacional de San Antonio Abad del Cusco, Cusco, Perú
    Percy Nuñez Vargas, Nadir Pallqui Camacho & Javier Silva Espejo

    Universidad Autónoma del Beni José Ballivián, Trinidad, Bolivia
    Guido Pardo & Vincent Vos

    Universidad Regional Amazónica Ikiam, Ikiam, Ecuador
    Maria Cristina Peñuela-Mora

    Broward County Parks Recreation, Oakland Park, FL, USA
    John Pipoly

    Keller Science Action Center, Field Museum, Chicago, IL, USA
    Nigel Pitman

    Instituto de Ciencias Naturales, Universidad Nacional de Colombia, Bogotá, Colombia
    Adriana Prieto & Agustín Rudas

    Institute of Research for Forestry Development (INDEFOR), Universidad de los Andes, Mérida, Venezuela
    Hirma Ramirez-Angulo

    Socioecosistemas y Cambio Climatico, Fundacion Con Vida, Medellín, Colombia
    Zorayda Restrepo Correa

    Centro de Conservacion, Investigacion y Manejo de Areas Naturales, CIMA Cordillera Azul, Lima, Peru
    Lily Rodriguez Bayona

    Universidade Federal Rural da Amazônia, Belém, Brazil
    Rafael Salomão & Natalino Silva

    Departamento de Biología, Universidad de La Serena, La Serena, Chile
    Javier Silva Espejo

    Guyana Forestry Commission, Georgetown, Guyana
    James Singh

    Federal University of Alagoas, Maceió, Brazil
    Juliana Stropp

    Institute for Conservation Research, Escondido, CA, USA
    Varun Swamy

    Institute for Transport Studies, University of Leeds, Leeds, UK
    Joey Talbot

    Biodiversity Dynamics, Naturalis Biodiversity Center, Leiden, The Netherlands
    Hans ter Steege

    Systems Ecology, Free University, De Boelelaan 1087, Amsterdam, Netherlands
    Hans ter Steege

    Department of Biology, University of Florida, Gainesville, FL, USA
    John Terborgh

    Iwokrama International Centre for Rainforest Conservation and Development, Georgetown, Guyana
    Raquel Thomas

    Universidad de los Andes, Mérida, Venezuela
    Armando Torres-Lezama

    School of Geography, University of Nottingham, Nottingham, UK
    Geertje van der Heijden

    Van Hall Larenstein University of Applied Sciences, Leeuwarden, Netherlands
    Peter van der Meer

    Van der Hoult Forestry Consulting, Leeuwarden, The Netherlands
    Peter van der Hout

    Núcleo de Estudos e Pesquisas Ambientais – Universidade Estadual de Campinas, Campinas, Brazil
    Simone Aparecida Vieira

    Herbario del Sur de Bolivia, Universidad de San Francisco Xavier de Chuquisaca, Sucre, Bolivia
    Jeanneth Villalobos Cayo

    Tropenbos International, Wageningen, Netherlands
    Roderick Zagt

    A.E.-M. and D.G. designed the study with contributions from O.L.P., R.J.W.B., S.F. and M.J.P.S. A.E.-M. carried out the analyses with inputs from D.G., O.L.P., R.J.W.B., S.F., M.J.P.S., J.H.-.D. and H.L. A.E.-M. wrote a first draft with contributions from D.G., M.J.P.S., T.A.M.P., S.F. and O.L.P. O.L.P., R.J.W.B., S.F., M.J.P.S., T.R.B., K.-J.C., T.R.F., N.H., Y.M., B.M., B.H.M.J., A.M.-M., L.P., M.S., E.V.T., E.A.D., J.d.A.P., E.A., P.A.L., A.A., L.E.O.CA., A.A.-M., E.Arets, L.A., G.A.A.C., M.B., C.B., P.B.C., J.B., L.B., D.B., F.B., R.J.W.B., F.Brown, B.B., J.L.C., W.C., V.C.M., J.C., J.Comiskey, F.C.V., A.L.d.C., N.D.C., A.D.F., A.D., T.E., G.F.L., I.C.G.V., R.H., E.H.C., I.H.-C., E.J.-R., T.K., S.L., W.L., S.L.L., T.L., P.M., C.M., P.Morandi, D.N., A.J.N.L., P.N.V., E.A.d.O., N.P.C., G.Prado, J.Pipoly, M.P.-C., M.C.P.-M., N.P., A.P., C.Q., H.R.-A., S.M.d.A.R., M.R.-M., Z.R.C., L.R.B., A.R., R.S., J.S., J.S.E., N.S., J.Singh, C.S., J.Stroop, V.S., J.T., H.t.S., J.T., R.T., M.T., A.T.-L., L.V.G., G.v.d.H., P.v.d.M., P.v.d.H., R.V.M., S.A.V., J.V.C., V.V., R.Z. and P.Z. led field expeditions for data collection. O.L.P., J.L. and Y.M. conceived the RAINFOR forest plot network; D.G., E.G. and T.R.B. contributed to its development. O.L.P., R.J.W.B., T.R.F., T.R.B., A.M.‐M., L.E.O.C.A., E.A.D., B.M., B.H.M.J., N.H., E.V.T., J.C., E.G. and Y.M. coordinated data collection with the help of many co‐authors. O.L.P., T.R.B., S.L.L. and G.L.-G. conceived ForestPlots.net, and M.J.P.S., A.L., J.Peacock, G.P., K.L.L.M.L., D.G. and E.G. helped to develop it. All authors read and approved the manuscript (with important insights provided by O.L.P., L.P., H.t.S., T.E., W.C., S.M.d.A.R., E.G., E.A.d.O., P.M., M.J.P.S., D.B., G.v.d.H. and P.Z.). More

  • in

    Microevolution in our megadont relative

    1.
    Broom, R. Nature 142, 377–379 (1938).
    Article  Google Scholar 
    2.
    Grine, F. E. (ed.) Evolutionary History of the “Robust” Australopithecines (Aldine de Gruyter, 1988).

    3.
    Herries, A. I. R. & Adams, J. W. J. Hum. Evol. 65, 676–681 (2013).
    Article  Google Scholar 

    4.
    Pickering, R. et al. Nature 565, 226–229 (2019).
    CAS  Google Scholar 

    5.
    Herries, A. I. R. et al. Science 368, eaaw7293 (2020).
    CAS  Article  Google Scholar 

    6.
    Martin, J. M. et al. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-020-01319-6 (2020).

    7.
    Wood, B. & Grabowski, M. in Macroevolution (eds Serrelli, E. & Gontier, N.) 345–376 (Springer, 2015).

    8.
    Vrba, E. S. in Paleoclimate and Evolution with Emphasis on Human Origins (eds Vrba, E. S. et al.) 24–45 (Yale Univ. Press, 1995).

    9.
    Potts, R. & Faith, J. T. J. Hum. Evol. 87, 5–20 (2015).
    Article  Google Scholar  More

  • in

    Spatial patterns of pathogen prevalence in questing Ixodes ricinus nymphs in southern Scandinavia, 2016

    1.
    Pfäffle, M., Littwin, N., Muders, S. V. & Petney, T. N. The ecology of tick-borne diseases. Int. J. Parasitol. 43, 1059–1077 (2013).
    Article  Google Scholar 
    2.
    Han, B. A. & Yang, L. Predicting novel tick vectors of zoonotic disease. in ICML Workshop on #Data4Good: Machine Learning in Social Good Applications 71–75 (2016).

    3.
    de la Fuente, J., Estrada-Pena, A., Venzal, J. M., Kocan, K. M. & Sonenshine, D. E. Overview: ticks as vectors of pathogens that cause disease in humans and animals. Front. Biosci. 13, 6938–6946 (2008).
    Article  Google Scholar 

    4.
    Michelet, L. et al. High-throughput screening of tick-borne pathogens in Europe. Front. Cell. Infect. Microbiol. 4, 103 (2014).
    Article  Google Scholar 

    5.
    Estrada-Peña, A. & de la Fuente, J. The ecology of ticks and epidemiology of tick-borne viral diseases. Antiviral Res. 108, 104–128 (2014).
    Article  CAS  Google Scholar 

    6.
    Paul, R. E. L. et al. Environmental factors influencing tick densities over seven years in a French suburban forest. Parasit. Vectors 9, 309 (2016).
    Article  CAS  Google Scholar 

    7.
    Randolph, S. E. Tick-borne disease systems emerge from the shadows: the beauty lies in molecular detail, the message in epidemiology. Parasitology 136, 1403 (2009).
    CAS  Article  Google Scholar 

    8.
    Jore, S. et al. Multi-source analysis reveals latitudinal and altitudinal shifts in range of Ixodes ricinus at its northern distribution limit. Parasit. Vectors 4, 1–11 (2011).
    Article  Google Scholar 

    9.
    Bernstein, L. et al. Climate Change 2007: Synthesis Report. Contribution of working groups I, II and III to the fourth assessment report of the Intergovernmental Panel on Climate Change. Intergovernmental Panel on Climate Change [Core Writing Team, Pachauri, R.K and Reisinger, A. (eds.)] https://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr_full_report.pdf (2007).

    10.
    Kovats, R. S., Campbell-Lendrum, D. H., McMichael, A. J., Woodward, A. & Cox, J. S. Early effects of climate change: do they include changes in vector-borne disease?. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 356, 1057–1068 (2001).
    CAS  PubMed  PubMed Central  Google Scholar 

    11.
    Gage, K. L., Burkot, T. R., Eisen, R. J. & Hayes, E. B. Climate and vectorborne diseases. Am. J. Prev. Med. 35, 436–450 (2008).
    Article  Google Scholar 

    12.
    Medlock, J. M. et al. Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit. Vectors 6, 1–11 (2013).
    Article  Google Scholar 

    13.
    Andreassen, A. et al. Prevalence of tick borne encephalitis virus in tick nymphs in relation to climatic factors on the southern coast of Norway. Parasit. Vectors 5, 1–12 (2012).
    Article  Google Scholar 

    14.
    Soleng, A. et al. Distribution of Ixodes ricinus ticks and prevalence of tick-borne encephalitis virus among questing ticks in the Arctic Circle region of northern Norway. Ticks Tick. Borne. Dis. 9, 97–103 (2018).
    CAS  Article  Google Scholar 

    15.
    Kjelland, V. et al. Tick-borne encephalitis virus, Borrelia burgdorferi sensu lato, Borrelia miyamotoi, Anaplasma phagocytophilum and Candidatus Neoehrlichia mikurensis in Ixodes ricinus ticks collected from recreational islands in southern Norway. Ticks Tick. Borne. Dis. 9, 1098–1102 (2018).
    Article  Google Scholar 

    16.
    Paulsen, K. M. et al. Prevalence of tick-borne encephalitis virus in Ixodes ricinus ticks from three islands in north-western Norway. APMIS 123, 759–764 (2015).
    Article  Google Scholar 

    17.
    Kjær, L. J. et al. A large-scale screening for the taiga tick, Ixodes persulcatus, and the meadow tick, Dermacentor reticulatus, in southern Scandinavia, 2016. Parasit. Vectors 12, 338 (2019).
    Article  Google Scholar 

    18.
    Oechslin, C. P. et al. Prevalence of tick-borne pathogens in questing Ixodes ricinus ticks in urban and suburban areas of Switzerland. Parasit. Vectors 10, 558 (2017).
    Article  CAS  Google Scholar 

    19.
    Becker, N. S. et al. Recurrent evolution of host and vector association in bacteria of the Borrelia burgdorferi sensu lato species complex. BMC Genom. 17, 734 (2016).
    Article  CAS  Google Scholar 

    20.
    Bowman, A. S. & Nuttall, P. A. Ticks: Biology, Disease and Control (Cambridge University Press, Cambridge, 2004).
    Google Scholar 

    21.
    Hasle, G. et al. Transport of ticks by migratory passerine birds to Norway. J. Parasitol. 95, 1342–1351 (2009).
    PubMed  Google Scholar 

    22.
    Klitgaard, K. et al. Screening for multiple tick-borne pathogens in Ixodes ricinus ticks from birds in Denmark during spring and autumn migration seasons. Ticks Tick. Borne. Dis. 10, 546–552 (2019).
    PubMed  Google Scholar 

    23.
    Skarphédinsson, S. et al. Detection and identification of Anaplasma phagocytophilum, Borrelia burgdorferi, and Rickettsia helvetica in Danish Ixodes ricinus ticks. APMIS 115, 225–230 (2007).
    PubMed  Google Scholar 

    24.
    Fraenkel, C.-J., Garpmo, U. & Berglund, J. Determination of novel Borrelia genospecies in Swedish Ixodes ricinus ticks. J. Clin. Microbiol. 40, 3308–3312 (2002).
    CAS  PubMed  PubMed Central  Google Scholar 

    25.
    Wilhelmsson, P. et al. Prevalence, diversity, and load of Borrelia species in ticks that have fed on humans in regions of Sweden and Åland Islands, Finland with different Lyme borreliosis incidences. PLoS ONE 8, e81433 (2013).
    ADS  PubMed  PubMed Central  Google Scholar 

    26.
    Vennestrøm, J., Egholm, H. & Jensen, P. M. Occurrence of multiple infections with different Borrelia burgdorferi genospecies in Danish Ixodes ricinus nymphs. Parasitol. Int. 57, 32–37 (2008).
    PubMed  Google Scholar 

    27.
    Kjelland, V., Stuen, S., Skarpaas, T. & Slettan, A. Prevalence and genotypes of Borrelia burgdorferi sensu lato infection in Ixodes ricinus ticks in southern Norway. Scand. J. Infect. Dis. 42, 579–585 (2010).
    CAS  PubMed  Google Scholar 

    28.
    Klitgaard, K., Kjær, L. J., Isbrand, A., Hansen, M. F. & Bødker, R. Multiple infections in questing nymphs and adult female Ixodes ricinus ticks collected in a recreational forest in Denmark. Ticks Tick. Borne. Dis. 10, 1060–1065 (2019).
    PubMed  Google Scholar 

    29.
    Maraspin, V., Ruzic-Sabljic, E. & Strle, F. Lyme borreliosis and Borrelia spielmanii. Emerg. Infect. Dis. 12, 1177–1177 (2006).
    PubMed  PubMed Central  Google Scholar 

    30.
    Rudenko, N., Golovchenko, M., Grubhoffer, L. & Oliver, J. H. Updates on Borrelia burgdorferi sensu lato complex with respect to public health. Ticks Tick. Borne. Dis. 2, 123–128 (2011).
    Article  Google Scholar 

    31.
    Fertner, M. E., Mølbak, L., Pihl, T. P. B., Fomsgaard, A. & Bødker, R. First detection of tick-borne “Candidatus Neoehrlichia mikurensis” in Denmark 2011. Eurosurveillance 17, 20096 (2012).
    Google Scholar 

    32.
    Quarsten, H. et al. Candidatus Neoehrlichia mikurensis and Borrelia burgdorferi sensu lato detected in the blood of Norwegian patients with erythema migrans. Ticks Tick. Borne. Dis. 8, 715–720 (2017).
    CAS  Article  Google Scholar 

    33.
    Stuen, S., Granquist, E. G. & Silaghi, C. Anaplasma phagocytophilum—a widespread multi-host pathogen with highly adaptive strategies. Front. Cell. Infect. Microbiol. 3, 31 (2013).
    Article  CAS  Google Scholar 

    34.
    Fomsgaard, A. et al. Tick-borne encephalitis virus, Zealand, Denmark, 2011. Emerg. Infect. Dis. 19, 1171–1173 (2013).
    Article  Google Scholar 

    35.
    Jensen, P. M. et al. Transmission differentials for multiple pathogens as inferred from their prevalence in larva, nymph and adult of Ixodes ricinus (Acari: Ixodidae). Exp. Appl. Acarol. 71, 171–182 (2017).
    Article  Google Scholar 

    36.
    Lundkvist, Å., Wallensten, A., Vene, S. & Hjertqvist, M. Tick-borne encephalitis increasing in Sweden, 2011. Eurosurveillance 16, 19981 (2011).
    Article  Google Scholar 

    37.
    Svensson, J., Hunfeld, K.-P. & Persson, K. E. M. High seroprevalence of Babesia antibodies among Borrelia burgdorferi-infected humans in Sweden. Ticks Tick. Borne. Dis. 10, 186–190 (2019).
    Article  Google Scholar 

    38.
    Mørch, K., Holmaas, G., Frolander, P. S. & Kristoffersen, E. K. Severe human Babesia divergens infection in Norway. Int. J. Infect. Dis. 33, 37–38 (2015).
    Article  Google Scholar 

    39.
    Uhnoo, I. et al. First documented case of human babesiosis in Sweden. Scand. J. Infect. Dis. 24, 541–547 (2009).
    Article  Google Scholar 

    40.
    Dumler, J. S., Barat, N. C., Barat, C. E. & Bakken, J. S. Human granulocytic anaplasmosis and macrophage activation. Clin. Infect. Dis. 45, 199–204 (2007).
    CAS  Article  Google Scholar 

    41.
    Nilsson, K., Elfving, K. & Påhlson, C. Rickettsia helvetica in patient with meningitis, Sweden, 2006. Emerg. Infect. Dis. 16, 490–492 (2010).
    CAS  Article  Google Scholar 

    42.
    Frivik, J. O., Noraas, S., Grankvist, A., Wennerås, C. & Quarsten, H. En mann i 60-årene fra Sørlandet med intermitterende feber (In Norwegian). Tidsskr. Den Nor. legeforening 137, (2017).

    43.
    Grankvist, A. et al. Infections with the tick-borne bacterium ‘Candidatus Neoehrlichia mikurensis’ mimic noninfectious conditions in patients with B cell malignancies or autoimmune diseases. Clin. Infect. Dis. 58, 1716–1722 (2014).
    CAS  Article  Google Scholar 

    44.
    Welinder-Olsson, C., Kjellin, E., Vaht, K., Jacobsson, S. & Wenneras, C. First case of human ‘Candidatus Neoehrlichia mikurensis’ infection in a febrile patient with chronic lymphocytic leukemia. J. Clin. Microbiol. 48, 1956–1959 (2010).
    Article  Google Scholar 

    45.
    Rizzoli, A. et al. Ixodes ricinus and its transmitted pathogens in urban and peri-urban areas in Europe: new hazards and relevance for public health. Front. Public Health 2, 251 (2014).
    Article  Google Scholar 

    46.
    Michelitsch, A., Wernike, K., Klaus, C., Dobler, G. & Beer, M. Exploring the reservoir hosts of tick-borne encephalitis virus. Viruses vol. 11 (2019).

    47.
    Keesing, F. et al. Reservoir competence of vertebrate hosts for Anaplasma phagocytophilum. Emerg. Infect. Dis. 18, 2013–2016 (2012).
    Article  Google Scholar 

    48.
    Zhan, L. et al. Anaplasma phagocytophilum in livestock and small rodents. Vet. Microbiol. 144, 405–408 (2010).
    ADS  Article  Google Scholar 

    49.
    Portillo, A., Santibáñez, P., Palomar, A. M., Santibáñez, S. & Oteo, J. A. Candidatus Neoehrlichia mikurensis, Europe. New Microbes New Infect. 22, 30–36 (2018).
    CAS  Article  Google Scholar 

    50.
    Jenkins, A. et al. Detection of Candidatus Neoehrlichia mikurensis in Norway up to the northern limit of Ixodes ricinus distribution using a novel real time PCR test targeting the groEL gene. BMC Microbiol. 19, 199 (2019).
    PubMed  PubMed Central  Google Scholar 

    51.
    Obiegala, A. & Silaghi, C. Candidatus Neoehrlichia mikurensis—recent insights and future perspectives on clinical cases, vectors, and reservoirs in Europe. Curr. Clin. Microbiol. Rep. 5, 1–9 (2018).
    Google Scholar 

    52.
    Yabsley, M. J. & Shock, B. C. Natural history of zoonotic Babesia: role of wildlife reservoirs. Int. J. Parasitol. Parasites Wildl. 2, 18–31 (2013).
    PubMed  Google Scholar 

    53.
    Sprong, H. et al. Ixodes ricinus ticks are reservoir hosts for Rickettsia helvetica and potentially carry flea-borne Rickettsia species. Parasit. Vectors 2, 41 (2009).
    PubMed  PubMed Central  Google Scholar 

    54.
    Jaenson, T. G. T. et al. Risk indicators for the tick Ixodes ricinus and Borrelia burgdorferi sensu lato in Sweden. Med. Vet. Entomol. 23, 226–237 (2009).
    CAS  PubMed  Google Scholar 

    55.
    Hudson, P. J. et al. Tick-borne encephalitis virus in northern Italy: molecular analysis, relationships with density and seasonal dynamics of Ixodes ricinus. Med. Vet. Entomol. 15, 304–313 (2001).
    MathSciNet  CAS  Article  Google Scholar 

    56.
    Nazzi, F. et al. Ticks and Lyme borreliosis in an alpine area in northeast Italy. Med. Vet. Entomol. 24, 220–226 (2010).
    CAS  PubMed  Google Scholar 

    57.
    Hubalek, Z., Halouzka, J. & Juricova, Z. Longitudinal surveillance of the tick Ixodes ricinus for Borreliae. Med. Vet. Entomol. 17, 46–51 (2003).
    CAS  Article  Google Scholar 

    58.
    Lindström, A. & Jaenson, T. G. T. Distribution of the common tick, Ixodes ricinus (Acari: Ixodidae), in different vegetation types in southern Sweden. J. Med. Entomol. 40, 375–378 (2003).
    Article  Google Scholar 

    59.
    Mejlon, H. A. & Jaenson, T. G. T. Jaenson (1993) Seasonal prevalence of Borrelia burgdorferi in Ixodes ricinus in different vegetation types in Sweden. Scand. J. Infect. Dis. 25, 449–456 (2009).
    Article  Google Scholar 

    60.
    Tack, W. et al. Local habitat and landscape affect Ixodes ricinus tick abundances in forests on poor, sandy soils. For. Ecol. Manag. 265, 30–36 (2012).
    Google Scholar 

    61.
    Walhström, L. K. & Kjellander, P. Ideal free distribution and natal dispersal in female roe deer. Oecologia 103, 302–308 (1995).
    ADS  PubMed  Google Scholar 

    62.
    Zeman, P. Objective assessment of risk maps of tick-borne encephalitis and Lyme borreliosis based on spatial patterns of located cases. Int. J. Epidemiol. 26, 1121–1129 (1997).
    CAS  PubMed  Google Scholar 

    63.
    Jat, M. K. & Mala, S. Application of GIS and space-time scan statistic for vector born disease clustering. In ICEGOV ’17 Proceedings of the 10th International Conference on Theory and Practice of Electronic Governance (2017) https://doi.org/10.1145/3047273.3047361.

    64.
    Hönig, V. et al. Model of risk of exposure to Lyme borreliosis and tick-borne encephalitis virus-infected ticks in the border area of the Czech Republic (South Bohemia) and Germany (Lower Bavaria and Upper Palatinate). Int. J. Environ. Res. Public Health 16, 1173 (2019).
    PubMed Central  Google Scholar 

    65.
    Randolph, S. E. & Rogers, D. J. Fragile transmission cycles of tick-borne encephalitis virus may be disrupted by predicted climate change. Proc. Biol. Sci. 267, 1741–1744 (2000).
    CAS  PubMed  PubMed Central  Google Scholar 

    66.
    Stefanoff, P. et al. A Predictive model has identified tick-borne encephalitis high-risk areas in regions where no cases were reported previously, Poland, 1999–2012. Int. J. Environ. Res. Public Health 15, 677 (2018).
    PubMed Central  Google Scholar 

    67.
    Kjær, L. J. et al. Predicting and mapping human risk of exposure to Ixodes ricinus nymphs using climatic and environmental data, Denmark, Norway and Sweden, 2016. Eurosurveillance 24, 1800101 (2019).
    PubMed Central  Google Scholar 

    68.
    Kjær, L. J. et al. Predicting the spatial abundance of Ixodes ricinus ticks in southern Scandinavia using environmental and climatic data. Sci. Rep. 9, 18144 (2019).
    ADS  Google Scholar 

    69.
    Kjær, L. J. et al. Spatial data of Ixodes ricinus instar abundance and nymph pathogen prevalence, Scandinavia, 2016–2017. Collection https://doi.org/10.6084/m9.figshare.c.4938270.v1 (2020).
    Article  Google Scholar 

    70.
    Kjær, L. J. et al. Spatial data of Ixodes ricinus instar abundance and nymph pathogen prevalence, Scandinavia, 2016–2017. Sci. Data 7, 1–7 (2020).
    Google Scholar 

    71.
    Scharlemann, J. P. W. et al. Global data for ecology and epidemiology: a novel algorithm for temporal Fourier processing MODIS data. PLoS ONE 3, e1408 (2008).
    ADS  Article  Google Scholar 

    72.
    Corine Land Cover 2006 raster data. European Environment Agency https://www.eea.europa.eu/data-and-maps/data/clc-2006-raster (2010).

    73.
    Klitgaard, K., Chriél, M., Isbrand, A., Jensen, T. K. & Bødker, R. Identification of Dermacentor reticulatus ticks carrying Rickettsia raoultii on migrating jackal, Denmark. Emerg. Infect. Dis. 23, 2072–2074 (2017).
    Article  Google Scholar 

    74.
    Moutailler, S. et al. Co-infection of ticks: the rule rather than the exception. PLoS Negl. Trop. Dis. 10, e0004539 (2016).
    Article  CAS  Google Scholar 

    75.
    Reye, A. L. et al. Prevalence of tick-borne pathogens in Ixodes ricinus and Dermacentor reticulatus ticks from different geographical locations in Belarus. PLoS ONE 8, e54476 (2013).
    ADS  CAS  Article  Google Scholar 

    76.
    R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing https://www.r-project.org (2018).

    77.
    Cowling, D. W., Gardner, I. A. & Johnson, W. O. Comparison of methods for estimation of individual-level prevalence based on pooled samples. Prev. Vet. Med. 39, 211–225 (1999).
    CAS  Article  Google Scholar 

    78.
    ESRI. ArcGIS Desktop: Release 10.6.1. Redlands, CA: Environmental Systems Research Institute. (2017).

    79.
    Kulldorff M. and Information Management Services, I. SaTScanTM v9.6: Software for the spatial and space-time scan statistics www.satscan.org, 2018.

    80.
    Kleinman, K. rsatscan: Tools, classes, and methods for interfacing with SaTScan stand-alone software. (2015).

    81.
    Kulldorff, M. A spatial scan statistic. Communications in Statistics – Theory and Methods vol. 26 https://www.tandfonline.com/doi/abs/10.1080/03610929708831995 (1997).

    82.
    Han, J. et al. Using Gini coefficient to determining optimal cluster reporting sizes for spatial scan statistics. Int. J. Health Geogr. 15, 27 (2016).
    Article  Google Scholar 

    83.
    Kuhn., M., Contributions from Jed Wing, Steve Weston, Andre Williams, Chris Keefer, Allan Engelhardt, T., Cooper, Zachary Mayer, Brenton Kenkel, the R Core Team, Michael Benesty, Reynald Lescarbeau, A. Z. & Luca Scrucca, Yuan Tang, C. C. and T. H. caret: Classification and regression training. R package version 6.0-81. https://CRAN.R-project.org/package=caret. (2018).

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

    85.
    Vapnik, V., Golowich, S. E. & Smola, A. Support vector method for function approximation, regression estimation, and signal processing. in Advances in Neural Information Processing Systems 9 (eds. Mozer, M., Jordan, M. & Petsche, T.) 281–287 (MIT Press., 1997).

    86.
    Sanz, H., Valim, C., Vegas, E., Oller, J. M. & Reverter, F. SVM-RFE: selection and visualization of the most relevant features through non-linear kernels. BMC Bioinform. 19, 432 (2018).
    Article  Google Scholar 

    87.
    Ghojogh, B., Ca, B., Crowley, M. & Ca, M. The theory behind overfitting, cross validation, rRegularization, bagging, and boosting: tTutorial. https://arxiv.org/abs/1905.12787 [stat.ML] 1–23 (2019).

    88.
    Skarphédinsson, S., Jensen, P. M. & Kristiansen, K. Survey of tickborne infections in Denmark. Emerg. Infect. Dis. 11, 1055–1061 (2005).
    Article  Google Scholar 

    89.
    Quarsten, H., Skarpaas, T., Fajs, L., Noraas, S. & Kjelland, V. Tick-borne bacteria in Ixodes ricinus collected in southern Norway evaluated by a commercial kit and established real-time PCR protocols. Ticks Tick. Borne. Dis. 6, 538–544 (2015).
    CAS  Article  Google Scholar 

    90.
    Wilhelmsson, P. et al. Prevalence and diversity of Borrelia species in ticks that have bitten humans in Sweden. J. Clin. Microbiol. 48, 4169–4176 (2010).
    CAS  PubMed  PubMed Central  Google Scholar 

    91.
    Strnad, M., Hönig, V., Růžek, D., Grubhoffer, L. & Rego, R. O. M. Europe-wWide meta-analysis of Borrelia burgdorferi sensu lato prevalence in questing Ixodes ricinus ticks. Appl. Environ. Microbiol. 83, 3838 (2017).
    Google Scholar 

    92.
    Mysterud, A. et al. Tick abundance, pathogen prevalence, and disease incidence in two contrasting regions at the northern distribution range of Europe. Parasit. Vectors 11, 309 (2018).
    PubMed  PubMed Central  Google Scholar 

    93.
    Severinsson, K., Jaenson, T. G., Pettersson, J., Falk, K. & Nilsson, K. Detection and prevalence of Anaplasma phagocytophilum and Rickettsia helvetica in Ixodes ricinus ticks in seven study areas in Sweden. Parasit. Vectors 3, 66 (2010).
    PubMed  PubMed Central  Google Scholar 

    94.
    Karlsson, M. E. & Andersson, M. O. Babesia species in questing Ixodes ricinus, Sweden. Ticks Tick. Borne. Dis. 7, 10–12 (2016).
    PubMed  Google Scholar 

    95.
    Øines, Ø., Radzijevskaja, J., Paulauskas, A. & Rosef, O. Prevalence and diversity of Babesia spp. in questing Ixodes ricinus ticks from Norway. Parasit. Vectors 5, 156 (2012).
    PubMed  PubMed Central  Google Scholar 

    96.
    Andersson, M., Bartkova, S., Lindestad, O. & Råberg, L. Co-Infection with ‘Candidatus Neoehrlichia mikurensis’ and Borrelia afzelii in Ixodes ricinus Ticks in Southern Sweden. Vector-Borne Zoonotic Dis. 13, 438–442 (2013).
    PubMed  Google Scholar 

    97.
    Pedersen, B. N. et al. Distribution of Neoehrlichia mikurensis in Ixodes ricinus ticks along the coast of Norway: the western seaboard is a low-prevalence region. Zoonoses Public Health https://doi.org/10.1111/zph.12662 (2019).
    Article  PubMed  Google Scholar 

    98.
    Kantsø, B., Bo Svendsen, C., Moestrup Jensen, P., Vennestrøm, J. & Krogfelt, K. A. Seasonal and habitat variation in the prevalence of Rickettsia helvetica in Ixodes ricinus ticks from Denmark. Ticks Tick. Borne. Dis. 1, 101–103 (2010).
    PubMed  Google Scholar 

    99.
    Solano-Gallego, L., Sainz, Á., Roura, X., Estrada-Peña, A. & Miró, G. A review of canine babesiosis: the European perspective. Parasit. Vectors 9, 336 (2016).
    Article  CAS  Google Scholar 

    100.
    Randolph, S. E. The shifting landscape of tick-borne zoonoses: tick-borne encephalitis and Lyme borreliosis in Europe. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 356, 1045–1056 (2001).
    CAS  PubMed  PubMed Central  Google Scholar 

    101.
    Sumilo, D. et al. Tick-borne encephalitis in the Baltic States : Identifying risk factors in space and time. Int. J. Med. Microbiol. 296(Suppl), 76–79 (2006).
    PubMed  Google Scholar 

    102.
    Sumilo, D. et al. Socio-economic factors in the differential upsurge of tick-borne encephalitis in Central and Eastern Europe. Rev. Med. Virol. 18, 81–95 (2008).
    PubMed  Google Scholar 

    103.
    Randolph, S. E., Green, R. M., Peacey, M. F. & Rogers, D. J. Seasonal synchrony : the key to tick-borne encephalitis foci identified by satellite data. Parasitology 121, 15–23 (2000).
    PubMed  Google Scholar 

    104.
    Halos, L. et al. Ecological factors characterizing the prevalence of bacterial tick-borne pathogens in Ixodes ricinus ticks in pastures and woodlands. Appl. Environ. Microbiol. 76, 4413–4420 (2010).
    CAS  PubMed  PubMed Central  Google Scholar 

    105.
    Sjörs, H. Nordisk växtgeografi (in Swedish) (Bonniers, Scandinavian University Books, 1967).
    Google Scholar  More

  • in

    A first look at the metabolic rate of Greenland sharks (Somniosus microcephalus) in the Canadian Arctic

    1.
    Cavicchioli, R., Amils, R., Wagner, D. & McGenity, T. Life and applications of extremophiles. Environ. Microbiol. 13, 1903–1907 (2011).
    Article  Google Scholar 
    2.
    Riesch, R., Tobler, M. & Plath, M. Extremophile Fishes (Springer, New York, 2015).
    Google Scholar 

    3.
    Wharton, D. A. Life at the Limits: Organisms in Extreme Environments (Cambridge University Press, Cambridge, 2007).
    Google Scholar 

    4.
    Lear, K. O. et al. Divergent field metabolic rates highlight the challenges of increasing temperatures and energy limitation in aquatic ectotherms. Oecologia 193, 311–323 (2020).
    ADS  Article  Google Scholar 

    5.
    Elliott, K. H. et al. High flight costs, but low dive costs, in auks support the biomechanical hypothesis for flightlessness in penguins. Proc. Natl. Acad. Sci. 110, 9380–9384 (2013).
    ADS  CAS  Article  PubMed  PubMed Central  Google Scholar 

    6.
    Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Toward a metabolic theory of ecology. Ecology 85, 1771–1789 (2004).
    Article  Google Scholar 

    7.
    Clarke, A. & Johnston, N. M. Scaling of metabolic rate with body mass and temperature in teleost fish. J. Anim. Ecol. 68, 893–905 (1999).
    Article  Google Scholar 

    8.
    Schulte, P. M. The effects of temperature on aerobic metabolism: towards a mechanistic understanding of the responses of ectotherms to a changing environment. J. Exp. Biol. 218, 1856–1866 (2015).
    Article  Google Scholar 

    9.
    Kleiber, M. Body size and metabolism. ENE 1, 315–353 (1932).
    Google Scholar 

    10.
    Glazier, D. S. A unifying explanation for diverse metabolic scaling in animals and plants. Biol. Rev. 85, 111–138 (2010).
    Article  PubMed  Google Scholar 

    11.
    Jerde, C. L. et al. Strong evidence for an intraspecific metabolic scaling coefficient near 0.89 in fish. Front. Physiol. 10, 1166 (2019).
    Article  PubMed  PubMed Central  Google Scholar 

    12.
    van der Meer, J. Metabolic theories in ecology. Trends Ecol. Evol. 21, 136–140 (2006).
    Article  PubMed  Google Scholar 

    13.
    Luongo, S. M. & Lowe, C. G. Seasonally acclimated metabolic Q10 of the California horn shark, Heterodontus francisci. J. Exp. Mar. Bio. Ecol. 503, 129–135 (2018).
    Article  Google Scholar 

    14.
    White, C. R., Alton, L. A. & Frappell, P. B. Metabolic cold adaptation in fishes occurs at the level of whole animal, mitochondria and enzyme. Proc. R. Soc. B Biol. Sci. 279, 1740–1747 (2011).
    Article  CAS  Google Scholar 

    15.
    Krogh, A. The Quantitative Relation Between Temperature and Standard Metabolism in Animals (Internationale Zeitschrift fuÈr Physikalisch-Chemische Biologie, New York, 1914).
    Google Scholar 

    16.
    Messamah, B., Kellermann, V., Malte, H., Loeschcke, V. & Overgaard, J. Metabolic cold adaptation contributes little to the interspecific variation in metabolic rates of 65 species of Drosophilidae. J. Insect Physiol. 98, 309–316 (2017).
    CAS  Article  PubMed  Google Scholar 

    17.
    Holeton, G. F. Metabolic cold adaptation of polar fish: fact or artefact?. Physiol. Zool. 47, 137–152 (1974).
    Article  Google Scholar 

    18.
    Steffensen, J. F. Metabolic cold adaptation of polar fish based on measurements of aerobic oxygen consumption: fact or artefact? Artefact!. Comp. Biochem. Physiol. A. 132, 789–795 (2002).
    Article  Google Scholar 

    19.
    Peck, L. S. A cold limit to adaptation in the sea. Trends Ecol. Evol. 31, 13–26 (2016).
    Article  PubMed  Google Scholar 

    20.
    Chabot, D., Steffensen, J. F. & Farrell, A. P. The determination of standard metabolic rate in fishes. J. Fish Biol. 88, 81–121 (2016).
    CAS  Article  PubMed  Google Scholar 

    21.
    Lawson, C. L. et al. Powering ocean giants : the energetics of shark and ray megafauna. Trends Ecol. Evol. 34, 1–13 (2019).
    MathSciNet  Article  Google Scholar 

    22.
    Lowe, C. Metabolic rates of juvenile scalloped hammerhead sharks (Sphyrna lewini). Mar. Biol. 139, 447–453 (2001).
    Article  Google Scholar 

    23.
    Payne, N. L. et al. A new method for resolving uncertainty of energy requirements in large water breathers: the ‘mega-flume’ seagoing swim-tunnel respirometer. Methods Ecol. Evol. 6, 668–677 (2015).
    Article  Google Scholar 

    24.
    Byrnes, E. E., Lear, K. O., Morgan, D. L. & Gleiss, A. C. Respirometer in a box: development and use of a portable field respirometer for estimating oxygen consumption of large-bodied fishes. J. Fish Biol. 96, 1045–1050 (2020).
    CAS  Article  PubMed  Google Scholar 

    25.
    MacNeil, M. A. et al. Biology of the greenland shark Somniosus microcephalus. J. Fish Biol. 80, 991–1018 (2012).
    CAS  Article  PubMed  Google Scholar 

    26.
    Edwards, J. E. et al. Advancing research for the management of long-lived species: a case study on the Greenland shark. Front. Mar. Sci. 6, 12 (2019).
    Article  Google Scholar 

    27.
    Augustine, S., Lika, K. & Kooijman, S. A. L. M. Comment on the ecophysiology of the Greenland shark, Somniosus microcephalus. Polar Biol. 40, 2429–2433 (2017).
    Article  Google Scholar 

    28.
    Nielsen, J. et al. Eye lens radiocarbon reveals centuries of longevity in the Greenland shark (Somniosus microcephalus). Science 353, 702–704 (2016).
    ADS  CAS  Article  PubMed  Google Scholar 

    29.
    Watanabe, Y. Y., Lydersen, C., Fisk, A. T. & Kovacs, K. M. The slowest fish: Swim speed and tail-beat frequency of Greenland sharks. J. Exp. Mar. Biol. Ecol. 426–427, 5–11 (2012).
    Article  Google Scholar 

    30.
    Hussey, N. E. et al. Rescaling the trophic structure of marine food webs. Ecol. Lett. 17, 239–250 (2014).
    Article  Google Scholar 

    31.
    Devine, B. M., Wheeland, L. J. & Fisher, J. A. D. First estimates of Greenland shark (Somniosus microcephalus) local abundances in Arctic waters. Sci. Rep. 8, 1–10 (2018).
    CAS  Article  Google Scholar 

    32.
    Wilson, E. E. & Wolkovich, E. M. Scavenging: how carnivores and carrion structure communities. Trends Ecol. Evol. 26, 129–135 (2011).
    Article  Google Scholar 

    33.
    Lear, K. O. et al. Correlations of metabolic rate and body acceleration in three species of coastal sharks under contrasting temperature regimes. J. Exp. Biol. 220, 397–407 (2017).
    Article  Google Scholar 

    34.
    Killen, S. S., Atkinson, D. & Glazier, D. S. The intraspecific scaling of metabolic rate with body mass in fishes depends on lifestyle and temperature. Ecol. Lett. 13, 184–193 (2010).
    Article  Google Scholar 

    35.
    Lear, K. O., Whitney, N. M., Brewster, L. R. & Gleiss, A. C. Treading water: respirometer choice may hamper comparative studies of energetics in fishes. Mar. Freshw. Res. 70, 437–448 (2018).
    Article  Google Scholar 

    36.
    Whitney, N. M., Lear, K. O., Gaskins, L. C. & Gleiss, A. C. The effects of temperature and swimming speed on the metabolic rate of the nurse shark (Ginglymostoma cirratum, Bonaterre). J. Exp. Mar. Bio. Ecol. 477, 40–46 (2016).
    Article  Google Scholar 

    37.
    Sims, D. W. The effect of body size on the standard metabolic rate of the lesser spotted dogfish. J. Fish Biol. 48, 542–544 (1996).
    Article  Google Scholar 

    38.
    Semmens, J. M., Payne, N. L., Huveneers, C., Sims, D. W. & Bruce, B. D. Feeding requirements of white sharks may be higher than originally thought. Sci. Rep. 3, 10–13 (2013).
    Article  CAS  Google Scholar 

    39.
    Giacomin, M., Schulte, P. M. & Wood, C. M. Differential effects of temperature on oxygen consumption and branchial fluxes of urea, ammonia, and water in the dogfish shark (Squalus acanthias suckleyi). Physiol. Biochem. Zool. 90, 627–637 (2017).
    Article  PubMed  PubMed Central  Google Scholar 

    40.
    Lowe, C. G. Bioenergetics of free-ranging juvenile scalloped hammerhead sharks (Sphyrna lewini) in Kāne’ohe Bay, Ō’ahu, HI. J. Exp. Mar. Biol. Ecol. 278, 141–156 (2002).
    Article  Google Scholar 

    41.
    Ezcurra, J. M., Lowe, C. G., Mollet, H. F., Ferry, L. A. & O’Sullivan, J. B. Oxygen consumption rate of young-of-the-year white sharks, Carcharodon carcharias during transport to the Monterey Bay Aquarium. Glob. Perspect. Biol. Life Hist. 1, 17–26 (2012).
    Article  Google Scholar 

    42.
    Barnett, A. et al. The utility of bioenergetics modelling in quantifying predation rates of marine apex predators: ecological and fisheries implications. Sci. Rep. 7, 12982 (2017).
    ADS  CAS  Article  PubMed  PubMed Central  Google Scholar 

    43.
    Watanabe, Y. Y., Payne, N. L., Semmens, J. M., Fox, A. & Huveneers, C. Swimming strategies and energetics of endothermic white sharks during foraging. J. Exp. Biol. 222, 4 (2019).
    Article  Google Scholar 

    44.
    Secor, S. M. Specific dynamic action: a review of the postprandial metabolic response. J. Comp. Physiol. B 179, 1–56 (2009).
    ADS  Article  PubMed  Google Scholar 

    45.
    Auer, S. K., Dick, C. A., Metcalfe, N. B. & Reznick, D. N. Metabolic rate evolves rapidly and in parallel with the pace of life history. Nat. Commun. 9, 8–13 (2018).
    ADS  Article  CAS  Google Scholar 

    46.
    Drazen, J. C. & Seibel, B. A. Depth-related trends in metabolism of benthic and benthopelagic deep-sea fishes. Limnol. Oceanogr. 52, 2306–2316 (2007).
    ADS  CAS  Article  Google Scholar 

    47.
    Brett, J. R. & Groves, T. D. D. Physiological energetics. Fish Physiol. 8, 280–352 (1979).
    Google Scholar 

    48.
    Widdows, J. Application of calorimetric methods in ecological studies. Therm. Energy. Stud. Cell. Biol. Syst. 1, 182–215 (1987).
    Article  Google Scholar 

    49.
    Armstrong, J. B. & Schindler, D. E. Excess digestive capacity in predators reflects a life of feast and famine. Nature 476, 84–87 (2011).
    CAS  Article  PubMed  Google Scholar 

    50.
    Stirling, I. & McEwan, E. Caloric value of whole ringed seals (Phoca hispida) in relation to Polar Bear (Ursus maritimus) ecology and hunting behavior. Can. J. Zool. 53, 1021–1027 (1975).
    CAS  Article  PubMed  Google Scholar 

    51.
    Furey, N. B., Hinch, S. G., Mesa, M. G. & Beauchamp, D. A. Piscivorous fish exhibit temperature-influenced binge feeding during an annual prey pulse. J. Anim. Ecol. 85, 1307–1317 (2016).
    Article  PubMed  Google Scholar 

    52.
    Svendsen, M. B. S., Bushnell, P. G. & Steffensen, J. F. Design and setup of intermittent-flow respirometry system for aquatic organisms. J. Fish Biol. 88, 26–50 (2016).
    CAS  Article  PubMed  Google Scholar 

    53.
    Clark, T. D., Sandblom, E. & Jutfelt, F. Aerobic scope measurements of fishes in an era of climate change: respirometry, relevance and recommendations. J. Exp. Biol. 216, 2771–2782 (2013).
    Article  PubMed  Google Scholar 

    54.
    Leclerc, L.-M.E. et al. A missing piece in the Arctic food web puzzle? Stomach contents of Greenland sharks sampled in Svalbard, Norway. Polar Biol. 35, 1197–1208 (2012).
    Article  Google Scholar  More

  • in

    Increased mosquito abundance and species richness in Connecticut, United States 2001–2019

    Summary statistics
    To date, The Connecticut Agricultural Experiment Station (CAES) has collected and tested 4,602,240 female mosquitoes comprised of 47 species in 8 genera. Approximately 98% of these collections were obtained from 92 trapping sites in 73 towns throughout the state, while the remainder of collections were from an additional 365 supplemental sites sampled between 1996 and 2007. Eighty-eight percent of collections come from CDC Light Traps, CDC Gravid Traps and Biogents BG Sentinel Traps (beginning in 2012). There have been several other collection methods used throughout the years that account for 11.6% of the mosquitoes collected (S. Table 1). Overall, there was considerable variation in mosquito abundance, surveillance effort, species richness/evenness, and the proportion of single species detections across CT (Fig. 1). One clear trend was that surveillance effort was greatest in CT’s human population centers (predominately CT’s southwestern and central counties) where WNV is commonly detected and along the CT-Rhode Island border where EEEV is most commonly detected (Fig. 1A). Another noticeable visual trend was that species evenness tends to be higher in the eastern portion of CT (Fig. 1B).
    Figure 1

    Maps of total mosquito abundance (log10 transformed) (A), total number of trap nights (A), average annual mosquito species richness (B), average annual mosquito species evenness (B), and average annual prevalence of single species detections (C) across 87 mosquito surveillance sites throughout Connecticut, U.S. sampled with ground level CDC CO2-baited light traps from 2001 to 2019. (A) Point sizes represent abundance while colors represent trap-nights; (B) point sizes represent species richness while colors represent species evenness; (C) point sizes represent prevalence of single species detections. (A–C) Solid black lines represent county political boundaries. The figure was created in R V 3.6.3 using the following packages: ggplot2 and maps.

    Full size image

    Objective 1: annual collections of mosquito populations among sites
    Our first objective was to identify spatial and temporal linear and nonlinear trends in mosquito abundance among sites. We also examined coarse-scale correlations between statewide (i.e., annual) and site-wide abundance and weather and land classification variables. All regression results and tables are provided as supporting information in Supporting Information: Regression Tables.
    Mosquito abundance
    Temporal regressions
    After accounting for trapping effort, regression parameters estimating the relationship between site-level mosquito abundance and year of collection were positive using generalized linear mixed effects models (GLMMs) (“Year”—Estimate 0.03, t-value 9.11) and generalized additive mixed effects models (GAMMs) (“Year”—Est. 0.77, t-value 2.7, p = 0.007), suggesting that site-level mosquito abundance has increased in CT since 2001 (Fig. 2A,B): this trend resulted in a predicted 60% increase in annual abundance from 2001 to 2019. While these regressions identified possible increasing trends in site-level abundance, they provided an overall poor-fit to the data: AIC scores from fixed effect GLMMs were higher than random effects-only models (ΔAIC 415.1). This poor model fit may be in part driven by directly modeling Year as a fixed continuous effect; Year as a random categorical effect may better capture variation in mosquito collections30. Despite large differences in AIC scores between fixed and random effects-only models, we detected a pattern of increasing intercept values when examining “Year” as a random effect (S. Fig. 1), providing further evidence of an increasing temporal trend in site-level mosquito abundance.
    Figure 2

    Average annual mosquito abundance (A), number of trap nights (B), mosquito species richness (C), mosquito species evenness (D), the annual correlation between mosquito species richness and evenness (E), and the prevalence of single mosquito species detections (F) across 87 mosquito surveillance sites throughout Connecticut, U.S. sampled with ground level CDC CO2-baited light traps from 2001 – 2019. For (A)–(D) and (F), points represent the average across all sites, solid lines represent the standard error of the average, and dashed lines are added to aid interpreting each plot as a time series. For (E), points represent the average across all sites while solid lines represent the 95% CI of the correlation point estimate. The figure was created in R V 3.6.3 using base functions.

    Full size image

    Spatial regressions
    After accounting for trapping effort, regression parameters estimating the relationship between site-level mosquito abundance and latitude/longitude were positive using a GLMM (“Latitude (centered)”—Est. 0.49, t-value 5.48; Longitude (centered)”—Est. 0.20, t-value 4.78), indicating that mosquito abundance tends to increase on a south to north and west to east gradient (which reflects the overall transition in land cover from developed to forested in CT). The best fitting fixed effect GAMM included Longitude by Latitude smoothing terms, which also predicted positive relationships between abundance and site coordinates (Smoothing term 1: Est. 0.24, p = 0.06; Smoothing term 2: Est. 0.05, p = 0.67). GAMM predictions of site-level mosquito abundance were considerably more complex than GLMM predictions, yet still supported the overall trend of increasing abundance from south to north and west to east (S. Fig. 2). Overall, the fixed effect GLMMs provided an extremely poor fit to the data compared to random effects-only GLMMs (Latitude—ΔAIC 1092.7; Longitude—ΔAIC 1099.8). These poor model fits may be in part driven by directly modeling coordinate (i.e., site) as a fixed continuous effect: GAMM predictions that account for nonlinear relationships between abundance and spatial location may provide a more appropriate fit to the data while site as a categorical random effect in the GLMMs may better capture variation in mosquito collections30.
    Weather correlations
    When comparing statewide annual mosquito abundance to weather variables, we found no correlations between summer temperatures, spring temperatures or precipitation. This was despite detecting a slight annual increase in temperatures across all three seasons examined (average daily temperature GLMM Est., Season/Summer: 0.05 °C, Prior Spring: 0.02 °C, Prior Winter: 0.07 °C) and a slight annual decline in within season and prior spring precipitation (total precipitation GLMM Est., Season/Summer: − 4.23 mm, Prior Spring: − 3.38 mm; Prior Winter: 2.22 mm) in CT since 2001. However, we did find a positive correlation between total summer precipitation and annual statewide mosquito abundance (r = 0.50, CI 0.07–0.78).
    Land cover correlations
    When comparing total site-wide abundance to land cover classifications, we found positive correlations between percent land cover categorized as barren (r = 0.22, CI 0.01–0.41), forested wetland (r = 0.34, 0.14–0.52), and non-forested wetland (r = 0.21, 0.004–0.41). We also found a negative association in total site-level abundance and percent land cover categorized as grass (r = − 0.35, − 0.52 to − 0.15).
    Species richness
    Temporal regressions
    After accounting for trapping effort, regression parameters estimating the relationship between site-level species richness and year of collection were positive using both GLMMs (“Year (centered)”—Est. 0.10, t-value 9.46) and GAMMs (“Year”—Est. 1.78, t-value 1.93, p = 0.05) (Fig. 2C): this trend resulted in a predicted 10% increase in site-level species richness from 2001 to 2019. Overall, fixed effects GLMMs of species richness provided an overall poor fit to the data when compared to a random effects-only model (ΔAIC 319.37). However, we did observe a pattern of increasing intercept values when examining “Year” as a random effect (S. Fig. 3), further indicating that mosquito species richness has annually increased across sites in CT since 2001.
    Spatial regressions
    Similar to models of site-level mosquito abundance, GLMMs of species richness by coordinate predicted positive relationships (Latitude (centered): Est. 0.63, t-value = 2.11; Longitude (centered): Est. 1.26, t-value = 9.34), indicating the species richness tends to increase along a south to north, west to east gradient. The best fitting GAMM included Longitude by Latitude smoothing terms, which also predicted positive relationships between species richness and site coordinate (Smoothing term 1: Est. 1.45, p = 0.0001; Smoothing term 2: Est. 0.70, p = 0.05). The GAMM further predicted a complex relationship of species richness among sites, yet overall predicted richness was lowest in the southwest/central portions of CT (areas of greatest development) and highest along coastal/eastern portions of CT (areas of non-forested and forested wetlands) (S. Fig. 4). The fixed effect GLMMs provided very poor fits to the data compared with random effects-only models (Latitude: ΔAIC 953.01; Longitude: ΔAIC 871.93; see the above results for Site-level collections: spatial regressions for possible reasons for these poor fits).
    Weather correlations
    We found no correlations of note between mosquito species richness and seasonal temperatures and precipitation.
    Land cover correlations
    Positive correlations of note for site-level species richness included: coniferous forest (r = 0.25, 0.04–0.43), deciduous forest (r = 0.56, 0.40–0.69), and forested wetland (r = 0.43, 0.23–0.58). Negative correlations included: barren (r = − 0.30, − 0.48 to − 0.10), developed (r = − 0.66, − 0.77 to − 0.53), grass (r = − 0.24, − 0.43 to − 0.03), and open water (r = − 0.31, − 0.49 to − 0.11).
    Species evenness
    Temporal regressions
    Trends in species evenness were negative using both GLMMs (“Year”—Est. − 0.01, t-value − 7.86) and GAMMs (“Year (centered)”—Est. − 0.04, t-value − 5.58, p = 0.000) (Fig. 2D): this trend resulted in a predicted 12% decrease in site-level species evenness from 2001 to 2019. Similar to fixed effects GLMMs of species richness, fixed effects GLMMs of species evenness were less informative than a random effects-only model (ΔAIC 66.5). Declining intercept values were evident when evaluating “Year” as a random effect (S. Fig. 5), further supporting an overall annual decline in species evenness estimates among sites.
    Spatial regressions
    Similar to spatial models of species richness, GLMMs predicted positive relationships between species evenness and coordinate (Latitude (centered): Est. 0.36, t-value = 7.63; Longitude (centered): Est. 0.18, t-value = 8.54); the best fitting GAMM, which included Longitude by Latitude smoothing terms, also predicted positive relationships (Smoothing term 1: Est. 0.12, p = 0.01; Smoothing term 2: Est. 0.16, p = 0.004). GAMM predictions of site-level species evenness were equally complex to predictions of abundance and richness, and predicted evenness to be highest in southcentral and eastern CT (S. Fig. 6). Fixed effect GLMMs provided very poor fits to the data compared with random effects-only models (Latitude: ΔAIC 502.6; Longitude: ΔAIC 488.4; see the above results for Site-level collections: spatial regressions for possible reasons for these poor fits).
    Weather correlations
    We did find a negative correlation between statewide prior spring minimum temperatures and mosquito species evenness (r = − 0.49, − 0.77 to − 0.04).
    Land cover correlations
    Positive correlations of note for species evenness included: deciduous forest (r = 0.46, 0.28–0.61) and forested wetland (r = 0.22, 0.01–0.41). Negative correlations included: barren (r = − 0.37, − 0.54 to − 0.18), developed (r = − 0.45, − 0.60 to − 0.26), and open water (r = − 0.32, − 0.50 to − 0.12).
    Correlations between abundance, richness, and evenness
    The relationships between abundance, richness, and evenness varied depending on the scale examined. Across all years of data at the site-level, the correlation between abundance and richness was positive (r = 0.53, 0.36–0.67), the correlation between abundance and evenness as negative (r = − 0.35, − 0.52 to − 0.15), and there was no correlation of note between richness and evenness. Across all sites at the year-level, there were no correlations of note between abundance, richness, and evenness. Annual statewide correlations between richness and evenness (RRE) were positive for all years yet there was no noticeable annual trend in these correlations (Fig. 2E). Spatially, the average site-level RRE was 0.15 (± 0.03 SE). Furthermore, the magnitude and direction of RRE tended to increase on a south to north gradient (r = 0.31, 0.11–0.49), yet there was no apparent relationship in RRE along a west to east gradient (S. Fig. 7). We did detect a positive correlation between RRE and average maximum spring temperatures (r = 0.46, 0.01–0.76) as well as a positive correlation between RRE and percent land cover classified as coniferous forest (r = 0.23, 0.02–0.42).
    Single detection events
    Single detection events were defined as the prevalence of single species detections at a site (i.e., number of species with a single pool divided by species richness). Changes in single species detections could indirectly indicate range expansion among species (i.e., the prevalence of single detections decreases with time) and/or areas of unique mosquito diversity (i.e., the prevalence of single detections changes across space).
    Temporal regressions
    We detected no overall pattern of increasing/decreasing annual prevalence of single-species detections among sites (GLMM, “Year”—Est. − 0.13, t-value = − 1.12, p = 0.22; GAMM, “Year”—Est. 0.02, t value = − 0.31, p = 0.75) (Fig. 2F). These models were considered equivalent to a random effects-only GLMM (ΔAIC  More

  • in

    Ecological and social constraints combine to promote evolution of non-breeding strategies in clownfish

    1.
    Emlen, S. T. In Behavioural Ecology: An Evolutionary Approach (eds Krebs, J. R. & Davies, N. B.) 301–337 (Blackwell, Oxford, 1991).
    2.
    Hamilton, W. D. The evolution of social behavior. J. Theor. Biol. 7, 1–52 (1964).
    CAS  Article  Google Scholar 

    3.
    West-Eberhard, M. J. The evolution of social behavior by kin selection. Q. Rev. Biol. 50, 1–33 (1975).
    Google Scholar 

    4.
    Emlen, S. T. & Wrege, P. H. A test of alternate hypotheses for helping behavior in white-fronted bee-eaters of Kenya. Behav. Ecol. Sociobiol. 25, 303–319 (1989).
    Article  Google Scholar 

    5.
    Reeve, H. K., Westneat, D. F., Noon, W. A., Sherman, P. W. & Aquadro, C. F. DNA fingerprinting’ reveals high levels of inbreeding in colonies of the eusocial naked mole- rat. Proc. Natn. Acad. Sci. 87, 2496–2500 (1990).
    CAS  Article  Google Scholar 

    6.
    Brouwere, L., Heg, D. & Taborsky, M. Experimental evidence for helper effects in a cooperatively breeding fish. Behav. Ecol. 16, 667–673 (2005).
    Article  Google Scholar 

    7.
    Wiley, R. H. & Rabenold, K. N. The evolution of cooperative breeding by delayed reciprocity and queuing for favorable social positions. Evolution 38, 97–107 (1984).
    Article  Google Scholar 

    8.
    Kokko, H. & Johnstone, R. A. Social queuing in animal societies: a dynamic model of reproductive skew. Proc. R. Soc. Lond. B 266, 571–578 (1999).
    Article  Google Scholar 

    9.
    Woolfenden, G. E. & Fitzpatrick, J. V. The inheritance of territory in group-breeding birds. Bioscience 28, 104–108 (1978).
    Article  Google Scholar 

    10.
    Creel, S. R. & Waser, P. M. Inclusive fitness and reproductive strategies in dwarf mongooses. Behav. Ecol. 5, 339–348 (1994).
    Article  Google Scholar 

    11.
    Balshine-Earn, S., Neat, F., Reid, H. & Taborsky, M. Paying to stay or paying to breed? Field evidence for directbenefits of helping in a cooperatively breeding fish. Behav. Ecol. 9, 432–438 (1998).
    Article  Google Scholar 

    12.
    Emlen, S. T. The evolution of helping. I. an ecological constraints model. Am. Nat. 119, 29–53 (1982).
    Article  Google Scholar 

    13.
    Hatchwell, B. J. & Komoder, J. Ecological constraints, life history traits and the evolution of cooperative breeding. Anim. Behav. 59, 1079–1086 (2000).
    CAS  Article  Google Scholar 

    14.
    Komdeur, J. Importance of habitat saturation and territory quality for evolution of cooperative breeding in the Seychelles warbler. Nature 358, 493–495 (1992).
    Article  Google Scholar 

    15.
    Faulkes, C. G. et al. Ecological constraints drive social evolution in the Africa mole rats. Proc. R. Soc. B 264, 1619–1627 (1997).
    CAS  Article  Google Scholar 

    16.
    Bergmüller, R., Heg, D. & Taborsky, M. Helpers in a cooperatively breeding cichlid stay and pay or disperse and breed, depending on ecological constraints. Proc. R. Soc. B 272, 325–331 (2005).
    Article  Google Scholar 

    17.
    Koenig, W. D. & Pitelka, F. A. Relatedness and inbreeding avoidance: Counterploys in the communally nesting acorn woodpecker. Science 206, 1103–1105 (1979).
    CAS  Article  Google Scholar 

    18.
    Cant, M. A., Hodge, S. J., Bell, M. B. V., Gilchrist, J. S. & Nichols, H. J. Reproductive control via eviction (but not the threat of eviction) in banded mongooses. Proc. R. Soc. B 277, 2219–2226 (2010).
    Article  Google Scholar 

    19.
    Dey, C. J., Tan, J. Q. Y., O’Connor, C. M., Reddon, A. R. & Caldwell, R. J. Dominance network structure across reproductive contexts in the cooperatively breeding cichlid fish Neolamprologus pulcher. Curr. Zool. 61, 45–54 (2015).
    Article  Google Scholar 

    20.
    Cant, M. A. The role of threats in animal cooperation. Proc. R. Soc. B 278, 170–178 (2011).
    Article  Google Scholar 

    21.
    Sherman, P. W., Lacey, E. A., Reeve, H. K. & Keller, L. Forum: The eusociality continuum. Behav. Ecol. 6, 102–108 (1995).
    Article  Google Scholar 

    22.
    Hing, M. L., Klanten, O. S., Dowton, M. & Wong, M. Y. L. The right tools for the job: cooperative breeding theory and an evaluation of the methodological approaches to understanding the evolution and maintenance of sociality. Front. Ecol. Evol. 5, 100 (2017).
    Article  Google Scholar 

    23.
    Fricke, H. & Fricke, S. Monogamy and sex change by aggressive dominance in coral reef fish. Nature 266, 830–832 (1977).
    CAS  Article  Google Scholar 

    24.
    Buston, P. M. Size and growth modification in clownfish. Nature 424, 145–146 (2003).
    CAS  Article  Google Scholar 

    25.
    Buston, P. M. & Balshine, S. Cooperating in the face of uncertainty: a consistent framework for understanding the evolution of cooperation. Behav. Process. 76, 152–159 (2007).
    Article  Google Scholar 

    26.
    Kokko, H., Johnstone, R. A. & Wright, J. The evolution of parental and alloparental effort in cooperatively breeding groups: when should helpers pay to stay? Behav. Ecol. 13, 291–300 (2002).
    Article  Google Scholar 

    27.
    Buston, P. M., Bogdanowicz, S. M., Wong, A. & Harrison, R. G. Are clownfish groups composed of relatives? Analysis of microsatellite DNA variation in Amphiprion percula. Mol. Ecol. 16, 3671–3678 (2007).
    CAS  Article  Google Scholar 

    28.
    Buston, P. M. Does the presence of non-breeders enhance the fitness of breeders? An experimental analysis in the clown anemonefish Amphiprion percula. Behav. Ecol. Sociobiol. 57, 23–31 (2004).
    Article  Google Scholar 

    29.
    Buston, P. M. Territory inheritance in the clown anemonefish. Proc. R. Soc. B (Suppl.) 271, S252–S254 (2004).
    Google Scholar 

    30.
    Mariscal, R. N. The nature of the symbiosis between Indo-Pacific anemone fish and seaanemones. Mar. Biol. 6, 58–65 (1970).
    Article  Google Scholar 

    31.
    Verwey, J. Coral reef studies. I. The symbiosis between damselfishesand sea anemones in Batavia Bay. Treubia 12, 305–366 (1930).
    Google Scholar 

    32.
    Fautin, D. G. The anemonefish symbiosis: what is known and what is not. Symbiosis 10, 23–46 (1991).
    Google Scholar 

    33.
    Elliott, J. K., Elliott, J. M. & Mariscal, R. N. Host selection, location and association behaviors of anemonefishes in field settlement experiments. Mar. Biol. 122, 377–389 (1995).
    Article  Google Scholar 

    34.
    Elliott, J. & Mariscal, R. Coexistence of nine anemonefish species: differential host and habitat utilization, size and recruitment. Mar. Biol. 138, 23–36 (2001).
    Article  Google Scholar 

    35.
    Buston, P. M. Forcible eviction and prevention of recruitment in the clown anemonefish. Behav. Ecol. 14, 576–582 (2003).
    Article  Google Scholar 

    36.
    Almany, G. et al. Larval fish dispersal in a coral-reef seascape. Nat. Ecol. Evol. 1, 0148 (2017).
    Article  Google Scholar 

    37.
    Buston, P. M. & Garcia, M. B. An extraordinary life span estimate for the clown anemonefish Amphiprion percula. J. Fish. Biol. 70, 1710–1719 (2007).
    Article  Google Scholar 

    38.
    Buston, P. M. Mortality is associated with social rank in the clown anemonefish (Amphiprion percula). Mar. Biol. 143, 811–815 (2003).
    Article  Google Scholar 

    39.
    Eibl-Eibesfeldt, I. Beobachtungen und Versuche an Anemonenfischen (Amphiprion) der Maldiven und der Nicobaren. Z. Tierpsychol. 17, 1–10 (1960).
    Article  Google Scholar 

    40.
    Moyer, J. T. & Nakazono, A. Protandrous hermaphroditism in six species of the anemonefish genus Amphiprion in Japan. Jpn J. Ichthyol. 25, 101–106 (1978).
    Google Scholar 

    41.
    Buston, P. M. & Cant, M. A. A new perspective on size hierarchies in nature: patterns, causes and consequences. Oecologia 149, 362–372 (2006).
    Article  Google Scholar 

    42.
    Wong. M.Y. Ecological constraints and benefits of philopatry promote group-living in a social but non-cooperatively breeding fish. Proc. R. Soc. B 277, 353–358 (2010).

    43.
    Kokko, H., Johnstone, R. A. & Clutton-Brock, T. H. The evolution of cooperative breeding through group augmentationProc. R. Soc. Lond. B. 268, 187–196 (2001).
    CAS  Article  Google Scholar 

    44.
    Bourke, A. F. G. Principles of Social Evolution (Oxford University Press, 2011).

    45.
    Nonacs, P. Go high or go low? Adaptive evolution of high and low relatedness societies in social hymenoptera. Front. Ecol. Evol. https://doi.org/10.3389/fevo.2017.00087 (2017)

    46.
    Grinsted, L. & Field, J. Market forces influence helping behaviour in cooperatively breeding paper wasps. Nat. Commun. 8, 13750 (2017).
    CAS  Article  Google Scholar 

    47.
    Wong, Y. L. M., Buston, P. M., Munday, P. L. & Jones, G. P. The threat of punishment enforces peaceful cooperation and stabilizes queues in a coral-reef fish. Proc. Biol. Sci. USA 274, 1093–1099 (2007).
    Google Scholar 

    48.
    Rueger, T. et al. Reproductve control via the threat of eviction in the clown anemonefish. Proc. R. Soc. B. 285, 20181295 (2018).
    Article  Google Scholar 

    49.
    Barbasch, T. et al. Substantial plasticity of reproduction and parental care in response to local resource availability. Oikos https://doi.org/10.1111/oik.07674 (2020).

    50.
    Salles, O. C. et al. Strong habitat and weak genetic effects shape the lifetime reproductive success in a wild clownfish population. Ecol. Lett. 23, 265–273 (2020).
    Article  Google Scholar 

    51.
    Holbrook, S. J. & Schmitt, R. J. Growth, reproduction and survival of a tropical sea anemone (Actiniaria): benefits of hosting anemone fish. Coral Reefs 24, 67–73 (2005).
    Article  Google Scholar 

    52.
    Cleveland, A., Verde, E. A. & Lee, R. W. Nutritional exchange in a tropical tripartite symbiosis: direct evidence for the transfer of nutrients from anemonefish to host anemone and zooxanthellae. Mar. Biol. 158, 589–602 (2011).
    Article  Google Scholar 

    53.
    Porat, D. & Chadwick-Furman, N. E. Effects of anemonefish on giant sea anemones: expansion behavior, growth, and survival. Hydrobiologia 530, 513–520 (2004).
    Google Scholar 

    54.
    Réale, D., Reader, S. M., Sol, D., McDougall, P. T. & Dingemanse, N. J. Integrating animal temperament within ecology and evolution. Biol. Rev. Camb. Philos. Soc. 82, 291–318 (2007).
    Article  Google Scholar 

    55.
    Schmiege, P. F. P., D’Aloia, C. C. & Buston P. M. Anemonefish personalities influence the strength of mutualistic interactions with host sea anemones. Mar Biol. https://doi.org/10.1007/s00227-016-3053-1 (2017).

    56.
    Barbasch, T. A. & Buston, P. M. Plasticity and personality of parental care in the clown anemonefish. Anim. Behav. 136, 65–73 (2018).
    Article  Google Scholar 

    57.
    Wong, M. Y. L. et al. Brief communication: consistent behavioural traits and behavioural syndromes in pairs of the false clown anemonefish Amphiprion ocellaris. J. Fish. Biol. 83, 207–213 (2013).
    CAS  Article  Google Scholar 

    58.
    Muthoo, A. A non-technical introduction to bargaining theory. World Econ. 1, 145–166 (2000).
    Google Scholar 

    59.
    Cant, M. A. & Johnstone, R. A. How threats influence the evolutionary resolution of within-group conflict. Am. Nat. 173, 759-771 (2009).

    60.
    Buston, P. M. & Zink, A. G. Reproductive skew and the evolution of conflict resolution: a synthesis of transactional and tug-of-war models. Behav. Ecol. 20, 672–684 (2009).
    Article  Google Scholar 

    61.
    Dixson, D. L. et al. Coral reef fish smell leaves to find island homes. Proc. R. Soc. Lond. B Biol. Sci. 275, 2831–2839 (2008).
    Google Scholar 

    62.
    Dixson, D. L. et al. Experimental evaluation of imprinting and the role innate preference plays in habitat selection in a coral reef fish. Oecologia 174, 99–107 (2014).
    Article  Google Scholar  More

  • in

    Climate reverses directionality in the richness–abundance relationship across the World’s main forest biomes

    1.
    Wright, A. J., Wardle, D. A., Callaway, R. & Gaxiola, A. The overlooked role of facilitation in biodiversity experiments. Trends Ecol. Evol. 32, 383–390 (2017).
    Article  PubMed  Google Scholar 
    2.
    Pacala, S. W. & Levin, S. A. Biologically generated spatial pattern and the coexistence of competing species. In Spatial Ecology: The Role of Space in Population Dynamics and Interspecific Interactions (Eds Tilman, D. and Kareiva, P.) (Princeton University Press, Princeton, NJ, 1997).

    3.
    Liang, J. et al. Positive biodiversity–productivity relationship predominant in global forests. Science 354, aaf8957 (2016).

    4.
    Poorter, L. et al. Biodiversity and climate determine the functioning of Neotropical forests. Glob. Ecol. Biogeogr. 26, 1423–1434 (2017).
    Article  Google Scholar 

    5.
    Ratcliffe, S. et al. Modes of functional biodiversity control on tree productivity across the European continent. Glob. Ecol. Biogeogr. 25, 251–262 (2016).
    Article  Google Scholar 

    6.
    Scherer-Lorenzen, M. & Schulze, E. D. Forest Diversity and Function: Temperate and Boreal Systems, Vol. 176 (Springer Science & Business Media, 2005).

    7.
    Paquette, A. & Messier, C. The effect of biodiversity on tree productivity: from temperate to boreal forests. Glob. Ecol. Biogeogr. 20, 170–180 (2011).
    Article  Google Scholar 

    8.
    Liu, X. et al. Tree species richness increases ecosystem carbon storage in subtropical forests. Proc. R. Soc. B. 285, 20181240 (2018).
    Article  PubMed  Google Scholar 

    9.
    Holzwarth, F., Rüger, N. & Wirth, C. Taking a closer look: disentangling effects of functional diversity on ecosystem functions with a trait-based model across hierarchy and time. R. Soc. Open Sci. 2, 140541 (2015).
    ADS  Article  PubMed  PubMed Central  Google Scholar 

    10.
    Bongers, F. J. et al. Growth–trait relationships in subtropical forest are stronger at higher diversity. J. Ecol. 108, 256–266 (2020).
    Article  Google Scholar 

    11.
    Adair, E. C., Hooper, D. U., Paquette, A. & Hungate, B. A. Ecosystem context illuminates conflicting roles of plant diversity in carbon storage. Ecol. Lett. 21, 1604–1619 (2018).
    Article  Google Scholar 

    12.
    Aponte, C. et al. Structural diversity underpins carbon storage in Australian temperate forests. Glob. Ecol. Biogeogr. 29, 789–802 (2020).
    Article  Google Scholar 

    13.
    Waide, R. B. et al. The relationship between productivity and species richness. Annu. Rev. Ecol. Evol. Syst. 30, 257–300 (1999).
    Article  Google Scholar 

    14.
    Dormann, C. F., Schneider, H. & Gorges, J. Neither global nor consistent: a technical comment on the tree diversity-productivity analysis of Liang et al. (2016). BioRxiv, 524363 (2019).

    15.
    Srivastava, D. S. & Lawton, J. H. Why more productive sites have more species: an experimental test of theory using tree-hole communities. Am. Nat. 152, 510–529 (1998).
    CAS  Article  PubMed  Google Scholar 

    16.
    Gaston, K. J. Global patterns in biodiversity. Nature 405, 220 (2000).
    CAS  Article  PubMed  Google Scholar 

    17.
    McGlynn, T. P., Weiser, M. D. & Dunn, R. R. More individuals but fewer species: testing the ‘more individuals hypothesis’ in a diverse tropical fauna. Biol. Lett. 6, 490–493 (2010).
    Article  PubMed  PubMed Central  Google Scholar 

    18.
    Yee, D. A. & Juliano, S. A. Abundance matters: a field experiment testing the more individuals hypothesis for richness–productivity relationships. Oecologia 153, 153–162 (2007).
    ADS  Article  PubMed  PubMed Central  Google Scholar 

    19.
    Storch, D., Bohdalková, E. & Okie, J. The more‐individuals hypothesis revisited: the role of community abundance in species richness regulation and the productivity–diversity relationship. Ecol. Lett. 21, 920–937 (2018).
    Article  Google Scholar 

    20.
    Ratcliffe, S. et al. Biodiversity and ecosystem functioning relations in European forests depend on environmental context. Ecol. Lett. 20, 1414–1426 (2017).
    Article  PubMed  Google Scholar 

    21.
    Zhai, P., Pörtner, H. O., & Roberts, D. (Eds). Summary for policymakers. In: global warming of 1.5 °C. An IPCC special report on the impacts of global warming of 1.5 °C above pre-industrial levels and related global greenhouse gas emission pathways. The Context of Strengthening the Global Response to the Threat of Climate Change, Sustainable Development, and Efforts to Eradicate Poverty, 32 (2018).

    22.
    Shipley, B. Testing causal explanations in organismal biology: causation, correlation and structural equation modelling. Oikos 86, 374–382 (1999).
    Article  Google Scholar 

    23.
    Yoda, K., Kira, T., Ogawa, H. & Hozami, K. Self-thinning in overcrowded pure stands under cultivated and natural conditions. J. Biol. Osaka City Univ. 14, 107–129 (1963).
    Google Scholar 

    24.
    Westoby, M. The self-thinning rule. Adv. Ecol. Res. 14, 167–225 (1984).
    Article  Google Scholar 

    25.
    Jump, A. S. et al. Structural overshoot of tree growth with climate variability and the global spectrum of drought‐induced forest dieback. Glob. Change Biol. 23, 3742–3757 (2017).
    ADS  Article  Google Scholar 

    26.
    Lefcheck, J. S. piecewiseSEM: piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579 (2016).
    Article  Google Scholar 

    27.
    Schluter, D. Experimental evidence that competition promotes divergence in adaptive radiation. Science 266, 798–801 (1994).
    ADS  CAS  Article  PubMed  Google Scholar 

    28.
    Francis, A. P. & Currie, D. J. A globally consistent richness–climate relationship for angiosperms. Am. Nat. 161, 523–536 (2003).
    Article  PubMed  Google Scholar 

    29.
    Grossiord, C. et al. Tree diversity does not always improve resistance of forest ecosystems to drought. Proc. Natl Acad. Sci. USA 111, 14812–14815 (2014).
    ADS  CAS  Article  PubMed  Google Scholar 

    30.
    Lavergne, S., Mouquet, N., Thuiller, W. & Ronce, O. Biodiversity and climate change: integrating evolutionary and ecological responses of species and communities. Annu. Rev. Ecol. Evol. Syst. 41, 321–350 (2010).
    Article  Google Scholar 

    31.
    Schemske, D. W., Mittelbach, G. G., Cornell, H. V., Sobel, J. M. & Roy, K. Is there a latitudinal gradient in the importance of biotic interactions? Annu. Rev. Ecol. Evol. Syst. 40, 245–269 (2009).
    Article  Google Scholar 

    32.
    Meehl, G. A. et al. Global climate projections. In Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (eds Solomon, S. et al.) (Cambridge University Press, Cambridge, United Kingdom and New York, NY, 2007).

    33.
    Lieth, H. Modelling the Primary Productivity of the Earth. Nature and Resources. UNESCO, VIII, 2:5-10 (1972).

    34.
    Lieth, H. Modeling the primary productivity of the world. In Primary Productivity of the Biosphere (eds Lieth, H. & Whittaker, R. H.) 237–264 (Springer-Verlag, New York, NY, 1975).

    35.
    Grieser, J., Gommes, R., Cofield, S. & Bernardi, M. In World Maps of Climatological Net Primary Production of Biomass, NPP (Food and Agriculture Organization of the United Nations, GEONETWORK, FAO, Rome, 2006).

    36.
    Revadekar, J. V. et al. Impact of altitude and latitude on changes in temperature extremes over South Asia during 1971–2000. Int. J. Climatol. 33, 199–209 (2013).
    Article  Google Scholar 

    37.
    Ali, S. et al. Effect of altitude on forest soil properties at Northern Karakoram. Eurasia. Soil Sci. 52, 1159–1169 (2019).
    ADS  Article  Google Scholar 

    38.
    Zhu, Z. X., Nizamani, M. M., Sahu, S. K., Kunasingam, A. & Wang, H. F. Tree abundance, richness, and phylogenetic diversity along an elevation gradient in the tropical forest of Diaoluo Mountain in Hainan, China. Acta Oecol. 101, 103481 (2019).
    Article  Google Scholar 

    39.
    Malizia, A. et al. Elevation and latitude drives structure and tree species composition in Andean forests: results from a large-scale plot network. PLoS ONE 15, e0231553 (2020).
    CAS  Article  PubMed  PubMed Central  Google Scholar 

    40.
    Homeier, J., Breckle, S. W., Günter, S., Rollenbeck, R. T. & Leuschner, C. Tree diversity, forest structure and productivity along altitudinal and topographical gradients in a species‐rich Ecuadorian montane rain forest. Biotropica 42, 140–148 (2010).
    Article  Google Scholar 

    41.
    Benavides, R. et al. Recruitment patterns of four tree species along elevation gradients in Mediterranean mountains: not only climate matters. For. Ecol. Manag. 360, 287–296 (2016).
    Article  Google Scholar 

    42.
    R Core Team R. A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna, 2013).

    43.
    Pinheiro, J., Bates, D., DebRoy, S. & Sarkar, D. R Core Team nlme: Linear and Nonlinear Mixed Effects Models. R Package Version 3.1-148. https://CRAN (2020). More

  • in

    Historical and projected future range sizes of the world’s mammals, birds, and amphibians

    Global land-use data
    For the historical time period 1700–2016, we used reconstructions of global cropland, pasture, and urban areas from the HYDE 3.2 dataset49 (available from https://doi.org/10.17026/dans-25g-gez3). Whilst HYDE 3.2 provides land-use data as far back as 10,000 BCE, we began our analysis in the year 1700, prior to which global land-use data are subject to increased uncertainty49,50. A total of 47 maps, including lower and upper uncertainty bounds, are available at 10-year intervals between 1700 and 2000, and at 1-year intervals between 2000 and 2016. These data were upscaled from their original spatial resolution of 0.083° to a 0.5° grid by summing up the cropland, pasture, or urban areas of all 0.083° grid cells contained in a given 0.5° cell.
    For the period 2020–2100, we used 0.5°-resolution 10-year time-step projections of global cropland, pasture, and urban areas from the AIM model51 (available from https://doi.org/10.7910/DVN/4NVGWA), covering Representative Concentration Pathways (RCPs) 2.6, 4.5, 6.0 and 8.5, and Shared Socio-economic Pathways (SSPs) 1–5. The dataset contains all possible combinations of these emission and socio-economic trajectories with the exception of RCP 2.6/SSP 3, and RCP 8.5/SSPs 1–4. The data were harmonised with the HYDE 3.2 data by adding the differences between HYDE 3.2 and AIM cropland, pasture and urban area maps in the year 2010 to the AIM future land use projections. We refer to refs. 27,28,29,52 for details of the emission and socio-economic pathways, and to ref. 28 for a comparison between the AIM model and other integrated assessment models.
    Global biome data
    We used the BIOME4 vegetation model53 (available from https://pmip2.lsce.ipsl.fr/synth/biome4.shtml) to simulate the distribution of global potential natural biomes between the years 1700 and 2000, and between 2020 and 2100 for each of the four climate-change scenarios considered here (RCPs 2.6, 4.5, 6.0, 8.5), at a spatial resolution of 0.5°. Inputs required by BIOME4 include global mean atmospheric CO2 concentration, and gridded monthly means of temperature, precipitation, and percent sunshine. Past and RCP-specific future CO2 levels were obtained from refs. 54 and 55, respectively. The climatic data were generated as follows. For the period 1700–1900, we used annual simulations from the HadCM3 climate model56 (available from https://esgf-node.llnl.gov/search/cmip5/; Experiments ‘past1000’ and ‘historical’, Ensemble ‘r1i1p1’). For the period 1901–2016, we used 0.5° resolution annual observational data57 (available from https://doi.org/10.5285/10d3e3640f004c578403419aac167d82). For the period 2020–2100, and for each RCP (2.6, 4.5, 6.0, 8.5), we used annual simulations from the HadGEM2-ES climate model58, the MIROC5 climate model59 and the CSIRO-Mk3.6.0 climate model60 (available from https://esgf-node.llnl.gov/search/cmip5/; for each climate model and each RCP, we used averages from Ensembles ‘r1i1p1’, ‘r2i1p1’, ‘r3i1p1’, ‘r4i1p1’). We downscaled and bias-corrected both the pre-1901 HadCM3 simulations and the future HadGEM2-ES, MIROC5, and CSIRO-Mk3.6.0 simulations using the delta method61. This method is based on applying the difference between simulated and observed climate at times at which both are available (here we used the 1900–1930 period for the historical data, and the year 2006 for the future data) to the simulated climate at points in time at which only simulated data exist (i.e., pre-1901 and post-2016) in order to correct systematic biases in the climate model61,62. The delta method also serves to spatially downscale the simulated climate to the 0.5° resolution of the observational data.
    For the computation of the global biome distribution at a point in time t, we used as inputs for BIOME4 the atmospheric CO2 concentration and gridded monthly climate values averaged across the time interval [t – 30 years, t]. Biome simulations were performed at 10-year intervals for both the historical and the future period. The complete time series of global biome simulations are available as Supplementary Movies 1–13.
    Estimation of species’ habitat ranges
    We estimated the geographic habitat ranges of an individual bird, mammal, and amphibian species through time following the general methodology in ref. 23. Our approach combines the following data:
    I.
    Spatial polygon data of species-specific extents of occurrence of all known birds63 (available from http://datazone.birdlife.org/species/requestdis), mammals, and amphibians64 (available from https://www.iucnredlist.org/).

    II.
    Species-specific biome requirements63,64 (data also available from the above websites).

    III.
    Maps of global potential natural biome distributions corresponding to the relevant climatic conditions through time (i.e., reconstructions for the past, and RCP-specific projections for the future).

    IV.
    Maps of global cropland, pasture, and urban areas through time (i.e., reconstructions for the past, and RCP- and SSP-specific projections for the future).

    The data I–IV were used to estimate the habitat range of individual species at a given point in time as illustrated in Fig. 4 and detailed in the following. In a first step, we used species-specific extents of occurrence (data I), which represent the outermost geographic limits of species’ observed, inferred or projected occurrences1. These spatial envelopes do not account for the distribution of natural or artificial land cover within that area, and therefore generally extend substantially beyond a species’ actual area of occupancy65,66. We first remapped extents of occurrence from their original spatial polygon format to a 0.083° resolution grid using the ‘rasterise’ function of the ‘raster’ package in R, which maps spatial polygons to those raster grid cells whose centres are contained within the polygons. For each species, we then determined the proportion of 0.083° cells contained in each 0.5° grid cell that represents the species’ extent of occurrence. This provides an estimate of the proportion of each 0.5° grid cell that is contained in the species’ extent of occurrence. Compared to the rasterising extent of occurrence directly to a 0.5° grid, this approach provides for more accurate estimates of species’ ranges and reduces the number of species that are not included in our analysis because their extents of occurrence do not overlap with any grid cell centre.
    Fig. 4: Method of estimating potential natural and actual range for the example of the bat-eared fox (Otocyon megalotis) in the year 1900.

    Here, for visualisation purposes, cropland, pasture, and urban areas were aggregated into one map; in reality, our method checks each of them separately against species’ artificial habitat preferences.

    Full size image

    In a second step, we refined the derived species-specific maps of the proportion of 0.5° grid cells contained in species’ extents of occurrence by combining them with species-specific biome requirements and maps of global biome distributions. Species-specific biome requirements (data II) include one or more habitat categories (cf. Supplementary Table 1), in which each species is known to occur. A species was estimated as being present in a grid cell contained in its previously derived extent of occurrence under the potential natural biome at a given point in time if the species’ list of habitat categories contained the local (i.e., grid cell-specific) potential natural biome at the relevant time (data III; see above). This required matching IUCN habitat categories (https://www.iucnredlist.org/resources/habitat-classification-scheme) with the biome categories of the Biome4 vegetation model, which was done as shown in Supplementary Table 1. In this way, we subset extents of occurrences by only retaining grid cells where the natural biome type is included in a species’ list of suitable habitat categories. The result of this step represents a species’ estimated potential natural habitat range (i.e., in the hypothetical absence of anthropogenic land use) at a given point in time.
    In a third step, we estimated actual habitat ranges by including maps of global land use through time. Each species’ actual habitat range at a given time was derived by removing any unsuitable anthropogenic land from the previously estimated potential natural range. Historical and projected future land use maps (data IV; see above) provide the fraction of each grid cell that is occupied by cropland, pasture or urban areas. These data were combined with information on which of these three artificial land cover types, if any, species can occur in, which is also included in the list of species’ biome requirements (data II). This allowed us, for each grid cell contained in a species’ potential natural range at a given time, to estimate the proportion of the grid cell that contained suitable habitat. A species’ actual habitat range size was then obtained as the sum of the areas of the remaining suitable habitat from all relevant grid cells.
    We applied the above method at each point in time for which global land use data is available (see above). In this way, we obtained potential natural ranges and actual ranges for 47 points in time between 1700 and 2016—using the baseline as well as lower and upper uncertainty bounds of the HYDE 3.2 land-use reconstructions—, and for nine points in time between 2020 and 2100—using the 16 combinations of future climatic and socio-economic pathways (see above), each of which, in turn, was considered based on climate data from three alternative models. Thus, we considered a total of 141 historical and 432 future scenarios.
    Since the global distribution of natural biomes varies over time as the result of (naturally or anthropogenically) changing climatic conditions, the sizes of potential natural habitat ranges are time-dependent. This motivates to consider range changes in relation to the potential natural ranges estimated at a particular reference time, for which we chose the year t0 = 1850, representing a modern pre-industrial baseline. Denoting the potential natural range and the actual range of a species i at a time t by (A_i^{{mathrm{potential}}}(t)) and (A_i^{{mathrm{actual}}}(t)), respectively, the range change associated with species i at time t as the result of the distribution of biomes and land use at that time was calculated at as

    $${Delta}A_ileft( t right) = 100{mathrm{% }} cdot left( {frac{{A_i^{{mathrm{actual}}}(t)}}{{A_i^{{mathrm{potential}}}(t_0)}} – 1} right).$$
    (1)

    Species whose potential natural habitat range size in the reference year t0 = 1850 (i.e., the range size estimated in the absence of anthropogenic land use and based on the global distribution of biomes in 1850) is zero, (A_i^{{mathrm{potential}}}left( {t_0} right) = 0), were not included in the analysis as, in this case, changes in range size are not defined. Based on the set (left{ {{Delta}A_ileft( t right)} right}_{i = 1,2, ldots }) of the individual range changes of all species through time, we calculated range change percentiles at each point in time (Fig. 1a), and determined the proportion of species that have experienced the loss of a given percentage of their baseline range (Fig. 1b). Similarly as in Eq. (1), we also computed the range change attributed only to climate-change-induced biome changes, (100% cdot left( {A_i^{{mathrm{potential}}}(t)/A_i^{{mathrm{potential}}}(t_0) – 1} right)) (Supplementary Fig. 1).
    Analyses were conducted using Matlab R2019a67 and R v3.6.368.
    Method discussion
    Whilst the available climate data for a given point in time only allows us to assign one primary natural biome type to each 0.5° grid cell, microclimates within cells may, in reality, result in the presence of different biomes in parts of a cell that are not represented in our data. By design of the approach used here, grid cells containing a non-primary biome that is suitable for a species, whilst the estimated primary biome is not, do not contribute to our estimation of the species’ habitat range. Conversely, grid cells containing a non-primary biome that is not suitable for a species, whilst the primary biome is suitable, would be included in their entirety in the species’ estimated range. This may lead us to underestimate the range sizes of species typically occurring in non-primary biomes in areas in which the estimated primary biomes are not suitable for the species, and to overestimate the range sizes of species typically occurring in the estimated primary biome in areas where other biomes also occur that are not suitable. Higher-resolution biome data could, in principle, reduce inaccuracies; however, generating such data in a reliable manner is not trivial. We are not aware of indications that this aspect of the approach would either systematically increase or decrease our overall estimates for range size changes across species in Fig. 1a.
    Our estimation of species’ habitat range sizes does not take into account habitat connectivity within or across grid cells. In principle, this can result in disconnected patches being included in a species’ estimated range, despite in reality being too small to represent potentially suitable habitat. However, neither species-specific data on the minimum size that spatially connected areas must not fall below before becoming non-viable nor reliable very-high-resolution land use and biome data, both of which would be needed to fully accommodate this issue, are currently available.
    Although species’ extents of occurrence are based not only on known, but also inferred and projected occurrences, the data remain very likely biased as the result of range contractions that occurred before the beginning of the systematic collection and mapping of species’ distributions, and that cannot be fully reconstructed. Whilst this may lead us to underestimate the absolute range sizes of species, it does not necessarily imply that we either systematically underestimate or overestimate the percentage change of species’ ranges through time.
    We chose the 0.5° resolution for our analysis as both the 1901–2016 observational climate data (and therefore also the pre-1901 and future climate data, which were downscaled using the observational data) and the projections of future land use are only available at this resolution. Attempts to further downscale these data would likely involve significant additional uncertainties. We are not aware of indications that an increase in the resolution of the analysis (if indeed the necessary datasets were available) would result in a systematic increase or decrease of either the absolute range sizes or the percentage change of range sizes relative to the baseline sizes, estimated here, at any point in time.
    Species-specific extents of occurrence and habitat preferences have been argued to be subject to uncertainty69; however, uncertainty estimates (quantitative or otherwise) are not provided with the data. In our main analysis, we therefore used the available data at face value. However, to verify that our results are not overly impacted by specific species, we performed the following bootstrapping analysis. Based on the set of species-specific range changes of all 16,919 species, estimated for the year 2016, we randomly sampled 16,919 values from this set with replacement a total of 104 times. For each of these 104 sets of range change estimates, we calculated 10%–90% percentiles analogous to Fig. 1a. For each percentile, we then calculated the mean and standard deviation of the computed 104 values. The result, shown in Supplementary Fig. 5, demonstrates that the uncertainties of our estimates with respect to specific species are very small, indicating that our results are robust with respect to potential uncertainties in the species data.
    Estimates of temporal delays in biome shifts in response to climatic changes70 are currently not available with the global coverage that would allow us to further refine our approach of assuming that biomes at a given point in time are determined by the climatic conditions in the preceding 30 years. This also applies to data on the dispersal speeds of plant functional types, and their effect on potential delays in colonisations of previously climatically unsuitable areas33; current studies on this topic are too spatially scarce to inform our approach. In our main analysis, we therefore followed the assumption commonly made in global vegetation models of no seed dispersal limitations71. However, to explore the impact of this assumption, we also repeated our analysis based on the extreme scenario of biomes not shifting at all between the present (year 2016) and 2100. The estimated range size changes (Supplementary Fig. 6) are quantitatively similar to the results of our main analysis (Fig. 3), consistent with our assessment of the overall stronger impact of land use compared to climate-driven biome changes. Qualitatively, i.e., in terms of how different RCP/SSP scenarios rank relative to each other, results are equivalent to those of our main analysis.
    As noted in the Introduction, our estimates of future habitat ranges represent upper estimates of species’ actual geographic distributions. In particular, our main analysis does not account for species’ ability to migrate to areas that will become suitable habitat at a future point in time but are not at present. However, our framework allows us to examine the effect of excluding such areas from the estimated habitat range. We repeated our analysis of future changes in habitat range sizes, but considered a grid cell as part of a species’ range only if the local biomes estimated for both the relevant point in the future and for the present (year 2016) were included in the species’ list of biome requirements. In other words, grid cells outside of species’ current potential natural habitat ranges were not counted towards their future range sizes, assuming that species are not able to migrate at all. This represents an extreme scenario that will underestimate most species’ mobility (e.g., over half of the species considered here can fly) and their ability to track biome shifts. Since the habitat range derived for a species in this manner is a subset of the one estimated in our main analysis, projected range losses based on this approach are, by design, higher (Supplementary Fig. 7). Qualitatively, results are equivalent to those in Fig. 3 in terms of how different RCP/SSP scenarios rank relative to each other.
    As the empirical data on species’ habitat preferences only provide categorical biome requirements, not continuous climatic envelopes, the method used here does not account for range changes due to changes in climatic conditions that are too small to manifest as biome changes. However, estimating precise climatic envelopes of species can be subject to considerable uncertainty and be highly sensitive to the way in which they are estimated (see below). By construction of the method used here, species’ ranges over time vary within the extents of occurrence provided with the empirical data, and do not exceed those. Justification for this assumption is provided by the fact that potential natural ranges (and, much more, actual ranges) are generally well-contained within extents of occurrence, with the former accounting for an average of 64% of the area of the latter in the reference year 1850, thus providing ample space for range shifts and expansions within the boundaries. Additional evidence that the restriction of habitat ranges to the extents of occurrence does not prevent significant range expansions can be seen in the sizeable number of species that have already experienced such range expansions (Fig. 1a and Supplementary Fig. 1) or are predicted to do so in future scenarios of strong global warming (Supplementary Fig. 1 and Supplementary Fig. 3a).
    Climate niche models estimate statistical relationships between climatic conditions and species’ spatial distributions, and apply these to climate projections in order to estimate future distribution patterns72. By design, they have great potential for mapping species’ distributions under a high degree of complexity in terms of possible predictor variables and their interactions, which has made the approach very useful in scenarios where the number of species, the geographic region and/or the temporal scale considered is relatively small so that statistical challenges are well-manageable73,74,75. In an analysis involving a large number of species, points in time, and different climatic and land-use scenarios considered here, the challenges commonly faced by climate nice models, specifically in terms of ensuring the robustness of the underlying statistical model and the estimated parameters, and avoiding unwanted artefacts in the extrapolation behaviour76,77,78,79,80,81, would be very difficult to manage. By operating directly and transparently on the empirical data of species’ extents of occurrence and biome requirements, and not being reliant on any particular statistical model or parameterisation, the approach used here provides the robustness needed at this scale of data23,82.
    Reporting summary
    Further information on research design is available in the Nature Research Reporting Summary linked to this article. More