Ecology

Ollerton, J. Pollinators & Pollination: Nature and Society (Pelagic Publishing, 2021).
Google Scholar 
Delaplane, K. S. Crop pollination By Bees: Evolution, Ecology, Conservation, and Management (CABI, 2021).
Google Scholar 
Aizen, M. A., Garibaldi, L. A., Cunningham, S. A. & Klein, A. M. Long-term global trends in crop yield and production reveal no current pollination shortage but increasing pollinator dependency. Curr. Biol. 18, 1572–1575 (2008).CAS 
PubMed 

Google Scholar 
Koh, I. et al. Modeling the status, trends, and impacts of wild bee abundance in the United States. Proc. Natl. Acad. Sci. 113, 140–145 (2016).ADS 
CAS 
PubMed 

Google Scholar 
Jordan, A., Patch, H. M., Grozinger, C. M. & Khanna, V. Economic dependence and vulnerability of United States agricultural sector on insect-mediated pollination service. Environ. Sci. Technol. 55, 2243–2253 (2021).ADS 
CAS 
PubMed 

Google Scholar 
Biddinger, D. J. & Rajotte, E. G. Integrated pest and pollinator management: Adding a new dimension to an accepted paradigm. Curr. Opin. Insect Sci. 10, 204–209 (2015).PubMed 

Google Scholar 
Egan, P. A., Dicks, L. V., Hokkanen, H. M. T. & Stenberg, J. A. Delivering integrated pest and pollinator management (IPPM). Trends Plant Sci. 25, 577–589 (2020).CAS 
PubMed 

Google Scholar 
Flöhr, A., Stenberg, J. A. & Egan, P. A. The Joint Economic Impact Level (jEIL): A Decision Metric for Integrated Pest and Pollinator Management. In Integrative Biological Control 17–38 (Springer, 2020).
Google Scholar 
Krupke, C. H., Hunt, G. J., Eitzer, B. D., Andino, G. & Given, K. Multiple routes of pesticide exposure for honey bees living near agricultural fields. PLoS ONE 7, e29268 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Traynor, K. S. et al. In-hive pesticide exposome: Assessing risks to migratory honey bees from in-hive pesticide contamination in the Eastern United States. Sci. Rep. 6, 1–16 (2016).
Google Scholar 
Mullin, C. A. et al. High levels of miticides and agrochemicals in North American apiaries: Implications for honey bee health. PLoS ONE 5, e9754 (2010).ADS 
PubMed 
PubMed Central 

Google Scholar 
Ravoet, J., Reybroeck, W. & de Graaf, D. C. Pesticides for apicultural and/or agricultural application found in belgian honey bee wax combs. Bull. Environ. Contam. Toxicol. 94, 543–548 (2015).CAS 
PubMed 
PubMed Central 

Google Scholar 
Liang, H. C., Bilon, N. & Hay, M. T. Analytical methods for pesticide residues. Water Environ. Res. 86, 2132–2155 (2014).CAS 

Google Scholar 
Calatayud-Vernich, P., Calatayud, F., Simó, E. & Picó, Y. Efficiency of QuEChERS approach for determining 52 pesticide residues in honey and honey bees. MethodsX 3, 452–458 (2016).PubMed 
PubMed Central 

Google Scholar 
Fernandez, M., Pico, Y. & Manes, J. Analytical methods for pesticide residue determination in bee products. J. Food Prot. 65, 1502–1511 (2002).CAS 
PubMed 

Google Scholar 
Vázquez, P. P., Lozano, A., Uclés, S., Ramos, M. M. G. & Fernández-Alba, A. R. A sensitive and efficient method for routine pesticide multiresidue analysis in bee pollen samples using gas and liquid chromatography coupled to tandem mass spectrometry. J. Chromatogr. A 1426, 161–173 (2015).PubMed 

Google Scholar 
Stoner, K. A., Cowles, R. S., Nurse, A. & Eitzer, B. D. Tracking pesticide residues to a plant genus using palynology in pollen trapped from honey bees (Hymenoptera: Apidae) at ornamental plant nurseries. Environ. Entomol. 48, 351–362 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Colwell, M. J., Williams, G. R., Evans, R. C. & Shutler, D. Honey bee-collected pollen in agro-ecosystems reveals diet diversity, diet quality, and pesticide exposure. Ecol. Evol. 7, 7243–7253 (2017).PubMed 
PubMed Central 

Google Scholar 
Favaro, R. et al. Botanical origin of pesticide residues in pollen loads collected by honeybees during and after apple bloom. Front. Physiol. 10, 1069 (2019).PubMed 
PubMed Central 

Google Scholar 
Tosi, S., Costa, C., Vesco, U., Quaglia, G. & Guido, G. A 3-year survey of Italian honey bee-collected pollen reveals widespread contamination by agricultural pesticides. Sci. Total Environ. 615, 208–218 (2018).ADS 
CAS 
PubMed 

Google Scholar 
Chaimanee, V., Chantawannakul, P., Khongphinitbunjong, K., Kamyo, T. & Pettis, J. S. Comparative pesticide exposure to Apis mellifera via honey bee-collected pollen in agricultural and non-agricultural areas of Northern Thailand. J. Apic. Res. 58, 720–729 (2019).
Google Scholar 
Friedle, C., Wallner, K., Rosenkranz, P., Martens, D. & Vetter, W. Pesticide residues in daily bee pollen samples (April–July) from an intensive agricultural region in Southern Germany. Environ. Sci. Pollut. Res. https://doi.org/10.1007/s11356-020-12318-2 (2021).Article 

Google Scholar 
Urbanowicz, C. et al. Low maize pollen collection and low pesticide risk to honey bees in heterogeneous agricultural landscapes. Apidologie 50, 379–390 (2019).
Google Scholar 
Stoner, K. A. & Eitzer, B. D. Using a hazard quotient to evaluate pesticide residues detected in pollen trapped from honey bees (Apis mellifera) in Connecticut. PLoS ONE 8, e77550 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
McArt, S. H., Fersch, A. A., Milano, N. J., Truitt, L. L. & Böröczky, K. High pesticide risk to honey bees despite low focal crop pollen collection during pollination of a mass blooming crop. Sci. Rep. 7, 1–10 (2017).
Google Scholar 
Calatayud-Vernich, P., Calatayud, F., Simó, E., Pascual Aguilar, J. A. & Picó, Y. A two-year monitoring of pesticide hazard in-hive: High honey bee mortality rates during insecticide poisoning episodes in apiaries located near agricultural settings. Chemosphere 232, 471–480 (2019).ADS 
CAS 
PubMed 

Google Scholar 
García-Valcárcel, A. I., Martínez-Ferrer, M. T., Campos-Rivela, J. M. & Hernando Guil, M. D. Analysis of pesticide residues in honeybee (Apis mellifera L.) and in corbicular pollen: Exposure in citrus orchard with an integrated pest management system. Talanta 204, 153–162 (2019).PubMed 

Google Scholar 
Fulton, C. A. et al. An assessment of pesticide exposures and land use of honey bees in Virginia. Chemosphere 222, 489–493 (2019).ADS 
CAS 
PubMed 

Google Scholar 
Calatayud-Vernich, P., Calatayud, F., Simó, E. & Picó, Y. Pesticide residues in honey bees, pollen and beeswax: Assessing beehive exposure. Environ. Pollut. 241, 106–114 (2018).CAS 
PubMed 

Google Scholar 
Ostiguy, N. et al. Honey bee exposure to pesticides: A four-year nationwide study. Insects 10, 13 (2019).PubMed Central 

Google Scholar 
Seeley, T. D. The honey bee colony as a superorganism. Am. Sci. 77, 546–553 (1989).ADS 

Google Scholar 
Thompson, H. M. & Maus, C. The relevance of sublethal effects in honey bee testing for pesticide risk assessment. Pest Manag. Sci. 63, 1058–1061 (2007).CAS 
PubMed 

Google Scholar 
Sponsler, D. B. & Johnson, R. M. Mechanistic modeling of pesticide exposure: The missing keystone of honey bee toxicology. Environ. Toxicol. Chem. 36, 871–881 (2017).CAS 
PubMed 

Google Scholar 
Gradish, A. E. et al. Comparison of pesticide exposure in honey bees (Hymenoptera: Apidae) and Bumble Bees (Hymenoptera: Apidae): implications for risk assessments. Environ. Entomol. 48, 12–21 (2019).PubMed 

Google Scholar 
Tosi, S. & Nieh, J. C. Lethal and sublethal synergistic effects of a new systemic pesticide, flupyradifurone (Sivanto®), on honeybees. Proc. R. Soc. B 286, 20190433 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Iverson, A., Hale, C., Richardson, L., Miller, O. & McArt, S. Synergistic effects of three sterol biosynthesis inhibiting fungicides on the toxicity of a pyrethroid and neonicotinoid insecticide to bumble bees. Apidologie 50, 733 (2019).CAS 

Google Scholar 
Siviter, H. et al. Agrochemicals interact synergistically to increase bee mortality. Nature 596, 389–392 (2021).ADS 
CAS 
PubMed 

Google Scholar 
Graham, K. K. et al. Identities, concentrations, and sources of pesticide exposure in pollen collected by managed bees during blueberry pollination. Sci. Rep. 11, 1–13 (2021).
Google Scholar 
EFSA. Guidance on the risk assessment of plant protection products on bees (Apis mellifera, Bombus spp. and solitary bees). EFSA J. https://doi.org/10.2903/j.efsa.2013.3295 (2013).Article 

Google Scholar 
EPA. Guidance for assessing pesticide risks to bees. (2014).USDA APHIS. Wax Sampling Protocol for the National Honey Bee Disease Survey. (2018).European Committee for Standardization. Foods of plant origin – Multimethod for the determination of pesticide residues using GC- and LC-based analysis following acetonitrile extraction/partitioning and clean-up by dispersive SPE – Modular QuEChERS-method. (2018).Couvillon, M. J. et al. Honey bee foraging distance depends on month and forage type. Apidologie 46, 61–70 (2015).
Google Scholar 
Knight, M. E. et al. An interspecific comparison of foraging range and nest density of four bumblebee (Bombus) species. Mol. Ecol. 14, 1811–1820 (2005).CAS 
PubMed 

Google Scholar 
McArt, S. H., Urbanowicz, C., Mccoshum, S., Irwin, R. E. & Adler, L. S. Landscape predictors of pathogen prevalence and range contractions in US bumblebees. Proc. R. Soc. B 284, 20172181 (2017).PubMed 
PubMed Central 

Google Scholar 
USDA NASS. USDA National Agricultural Statistics Service Cropland Data Layer. (2018).R Core Team. R: A Language and Environment for Statistical Computing. (2019).GraphPad Software. GraphPad Prism. (2017).Hartig, F. DHARMa: Residual Diagnostics for Hierarchical (Multi-Level/Mixed) Regression Models. (2019).Barton, K. MuMIn: Multi-Model Inference. (2019).Fox, J. & Weisburg, S. An {R} Companion to Applied Regression. (2011).Lefcheck, J. S. piecewiseSEM: Piecewise structural equation modeling in R for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579 (2015).
Google Scholar 
Hothorn, T., Bretz, F. & Westfall, P. Simultaneous inference in general parametric models. Biom. J. 50, 346–363 (2008).MathSciNet 
PubMed 
MATH 

Google Scholar 
Fox, J. RcmdrMisc: R Commander Miscellaneous Functions. (2020).Bhattacharya, M., Primack, R. B. & Gerwein, J. Are roads and railroads barriers to bumblebee movement in a temperate suburban conservation area?. Biol. Conserv. 109, 37–45 (2003).
Google Scholar 
Fragoso, F. P. & Brunet, J. Patch fidelity of honey bees and bumble bees differs and is affected by spatial configuration. In Entomological Society of America Annual Meeting, Plant-Insect Ecosystems (2021).Javorek, S. K., Mackenzie, K. E. & Vander Kloet, S. P. Comparative (Hymenoptera: Apoidea) on Lowbush Blueberry (Ericaceae: Vaccinium angustifolium). Ann. Entomol. Soc. Am. 95, 345–351 (2002).
Google Scholar 
Sandrock, C. et al. Impact of chronic neonicotinoid exposure on honeybee colony performance and queen supersedure. PLoS ONE 9, e103592 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Wood, S. C., Kozii, I. V., Koziy, R. V., Epp, T. & Simko, E. Comparative chronic toxicity of three neonicotinoids on New Zealand packaged honey bees. PLoS ONE 13, e0190517 (2018).PubMed 
PubMed Central 

Google Scholar 
Degrandi-Hoffman, G., Graham, H., Ahumada, F., Smart, M. & Ziolkowski, N. The economics of honey bee (Hymenoptera: Apidae) management and overwintering strategies for colonies used to pollinate almonds. J. Econ. Entomol. https://doi.org/10.1093/jee/toz213 (2019).Article 
PubMed 

Google Scholar 
Biddinger, D. J. et al. Comparative toxicities and synergism of apple orchard pesticides to Apis mellifera (L.) and Osmia cornifrons (Radoszkowski). PLoS ONE 8, e72587 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Johnson, R. M., Dahlgren, L., Siegfried, B. D. & Ellis, M. D. Acaricide, fungicide and drug Interactions in Honey Bees (Apis mellifera). PLoS ONE 8, e54092 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Jaffe, B. D., Lois, A. N. & Guédot, C. Effect of fungicide on pollen foraging by honeybees (Hymenoptera: Apidae) in cranberry differs by fungicide type. J. Econ. Entomol. 112, 499–503 (2019).CAS 
PubMed 

Google Scholar 
Sanchez-Bayo, F. & Goka, K. Pesticide residues and bees: A risk assessment. PLoS ONE 9, e94482 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Wang, Y., Zhu, Y. C. & Li, W. Interaction patterns and combined toxic effects of acetamiprid in combination with seven pesticides on honey bee (Apis mellifera L.). Ecotoxicol. Environ. Saf. 190, 110100 (2020).CAS 
PubMed 

Google Scholar 
Wang, Y., Zhu, Y. C. & Li, W. Comparative examination on synergistic toxicities of chlorpyrifos, acephate, or tetraconazole mixed with pyrethroid insecticides to honey bees (Apis mellifera L.). Environ. Sci. Pollut. Res. 27, 6971–6980 (2019).
Google Scholar 
Sgolastra, F. et al. Synergistic mortality between a neonicotinoid insecticide and an ergosterol-biosynthesis-inhibiting fungicide in three bee species. Pest Manag. Sci. 73, 1236–1243 (2017).CAS 
PubMed 

Google Scholar 
Azpiazu, C. et al. Chronic oral exposure to field-realistic pesticide combinations via pollen and nectar: effects on feeding and thermal performance in a solitary bee. Sci. Reports 9, 1–11 (2019).CAS 

Google Scholar 
Becher, M. A., Hildenbrandt, H., Hemelrijk, C. K. & Moritz, R. F. A. Brood temperature, task division and colony survival in honeybees: A model. Ecol. Modell. 221, 769–776 (2010).
Google Scholar 
Zhu, W., Schmehl, D. R., Mullin, C. A. & Frazier, J. L. Four common pesticides, their mixtures and a formulation solvent in the hive environment have high oral toxicity to honey bee larvae. PLoS ONE 9, e77547 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Dai, P. et al. Chronic toxicity of clothianidin, imidacloprid, chlorpyrifos, and dimethoate to Apis mellifera L. larvae reared in vitro. Pest Manag. Sci. 75, 29–36 (2019).CAS 
PubMed 

Google Scholar 
Campbell, J. B. et al. The fungicide Pristine® inhibits mitochondrial function in vitro but not flight metabolic rates in honey bees. J. Insect Physiol. 86, 11–16 (2016).CAS 
PubMed 

Google Scholar 
DesJardins, N. S. et al. A common fungicide, Pristine®, impairs olfactory associative learning performance in honey bees (Apis mellifera). Environ. Pollut. 288, 117720 (2021).CAS 
PubMed 

Google Scholar 
Fisher, A. et al. Colony field test reveals dramatically higher toxicity of a widely-used mito-toxic fungicide on honey bees (Apis mellifera). Environ. Pollut. 269, 115964 (2021).CAS 
PubMed 

Google Scholar 
Mahefarisoa, K. L., Simon Delso, N., Zaninotto, V., Colin, M. E. & Bonmatin, J. M. The threat of veterinary medicinal products and biocides on pollinators: A one health perspective. One Heal. 12, 100237 (2021).CAS 

Google Scholar 
Christen, V., Schirrmann, M., Frey, J. E. & Fent, K. Global transcriptomic effects of environmentally relevant concentrations of the neonicotinoids clothianidin, imidacloprid, and thiamethoxam in the brain of honey bees (Apis mellifera). Environ. Sci. Technol. 52, 7534–7544 (2018).ADS 
CAS 
PubMed 

Google Scholar 
Tison, L., Duer, A., Púčiková, V., Greggers, U. & Menzel, R. Detrimental effects of clothianidin on foraging and dance communication in honey bees. PLoS ONE 15, e0241134 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Tison, L., Rößner, A., Gerschewski, S. & Menzel, R. The neonicotinoid clothianidin impairs memory processing in honey bees. Ecotoxicol. Environ. Saf. 180, 139–145 (2019).CAS 
PubMed 

Google Scholar 
Morfin, N., Goodwin, P. H., Correa-Benitez, A. & Guzman-Novoa, E. Sublethal exposure to clothianidin during the larval stage causes long-term impairment of hygienic and foraging behaviours of honey bees. Apidologie 50, 595–605 (2019).CAS 

Google Scholar 
Yao, J., Zhu, Y. C. & Adamczyk, J. Responses of honey bees to lethal and sublethal doses of formulated clothianidin alone and mixtures. J. Econ. Entomol. 111, 1517–1525 (2018).CAS 
PubMed 

Google Scholar 
Bortolotti, L. et al. Effects of sub-lethal imidacloprid doses on the homing rate and foraging activity of honey bees. Bull. Insectology 56, 63–67 (2003).
Google Scholar 
Yang, E. C., Chuang, Y. C., Chen, Y. L. & Chang, L. H. Abnormal foraging behavior induced by sublethal dosage of imidacloprid in the honey bee (Hymenoptera: Apidae). J. Econ. Entomol. 101, 1743–1748 (2008).CAS 
PubMed 

Google Scholar 
Karahan, A., Cakmak, I., Hranitz, J. M., Karaca, I. & Wells, H. Sublethal imidacloprid effects on honey bee flower choices when foraging. Ecotoxicology 24, 2017–2025 (2015).CAS 
PubMed 

Google Scholar 
Dively, G. P., Embrey, M. S., Kamel, A., Hawthorne, D. J. & Pettis, J. S. Assessment of chronic sublethal effects of imidacloprid on honey bee colony health. PLoS ONE 10, 1–25 (2015).
Google Scholar 
Meikle, W. G. et al. Sublethal effects of imidacloprid on honey bee colony growth and activity at three sites in the U.S.. PLoS ONE https://doi.org/10.1371/journal.pone.0168603 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Wu, Y.-Y. et al. Sublethal effects of imidacloprid on targeting muscle and ribosomal protein related genes in the honey bee Apis mellifera L.. Sci. Rep. https://doi.org/10.1038/s41598-017-16245-0 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Wu, M.-C., Chang, Y.-W., Lu, K.-H. & Yang, E.-C. Gene expression changes in honey bees induced by sublethal imidacloprid exposure during the larval stage. Insect Biochem. Mol. Biol. 88, 12–20 (2017).CAS 
PubMed 

Google Scholar 
Peng, Y.-C. & Yang, E.-C. Sublethal dosage of imidacloprid reduces the microglomerular density of honey bee mushroom bodies. Sci. Rep. https://doi.org/10.1038/srep19298 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
de Almeida-Rossi, C., Roat, T. C., Tavares, D. A., Cintra-Socolowski, P. & Malaspina, O. Brain morphophysiology of africanized bee Apis mellifera exposed to sublethal doses of imidacloprid. Arch. Environ. Contam. Toxicol. 65, 234–243 (2013).PubMed 

Google Scholar 
Tosi, S., Burgio, G. & Nieh, J. C. A common neonicotinoid pesticide, thiamethoxam, impairs honey bee flight ability. Sci. Reports 7, 1–8 (2017).
Google Scholar 
Coulon, M. et al. Interactions between thiamethoxam and deformed wing virus can drastically impair flight behavior of honey bees. Front. Microbiol. 0, 766 (2020).
Google Scholar 
Shi, T.-F., Wang, Y.-F., Liu, F., Qi, L. & Yu, L.-S. Sublethal effects of the neonicotinoid insecticide thiamethoxam on the transcriptome of the honey bees (Hymenoptera: Apidae). J. Econ. Entomol. 110, 2283–2289 (2017).CAS 
PubMed 

Google Scholar 
Tesovnik, T. et al. Exposure of honey bee larvae to thiamethoxam and its interaction with Nosema ceranae infection in adult honey bees. Environ. Pollut. 256, 113443 (2020).CAS 
PubMed 

Google Scholar 
Friol, P. S., Catae, A. F., Tavares, D. A., Malaspina, O. & Roat, T. C. Can the exposure of Apis mellifera (Hymenoptera, Apiadae) larvae to a field concentration of thiamethoxam affect newly emerged bees?. Chemosphere 185, 56–66 (2017).ADS 
CAS 
PubMed 

Google Scholar 
Coulon, M. et al. Influence of chronic exposure to thiamethoxam and chronic bee paralysis virus on winter honey bees. PLoS ONE 14, e0220703 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Stoner, K. A. & Eitzer, B. D. Movement of soil-applied imidacloprid and thiamethoxam into nectar and pollen of squash (Cucurbita pepo). PLoS ONE 7, e39114 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Rolke, D., Persigehl, M., Peters, B., Sterk, G. & Blenau, W. Large-scale monitoring of effects of clothianidin-dressed oilseed rape seeds on pollinating insects in northern Germany: residues of clothianidin in pollen, nectar and honey. Ecotoxicology 25, 1691–1701 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wood, T. J., Kaplan, I., Zhang, Y. & Szendrei, Z. Honeybee dietary neonicotinoid exposure is associated with pollen collection from agricultural weeds. Proc. R. Soc. B 286, 1905 (2019).
Google Scholar 
Urlacher, E. et al. Measurements of chlorpyrifos levels in forager bees and comparison with levels that disrupt honey bee odor-mediated learning under laboratory conditions. J. Chem. Ecol. 42, 127–138 (2016).CAS 
PubMed 

Google Scholar 
Villalba, A., Maggi, M., Ondarza, P. M., Szawarski, N. & Miglioranza, K. S. B. Influence of land use on chlorpyrifos and persistent organic pollutant levels in honey bees, bee bread and honey: Beehive exposure assessment. Sci. Total Environ. 713, 136554 (2020).ADS 
CAS 
PubMed 

Google Scholar 
Traynor, K. S. et al. Pesticides in honey bee colonies: Establishing a baseline for real world exposure over seven years in the USA. Environ. Pollut. https://doi.org/10.1016/j.envpol.2021.116566 (2021).Article 
PubMed 

Google Scholar 
EPA Press Office. EPA Takes Action to Address Risk from Chlorpyrifos and Protect Children’s Health. (2021).Arena, M. & Sgolastra, F. A meta-analysis comparing the sensitivity of bees to pesticides. Ecotoxicology 233(23), 324–334 (2014).
Google Scholar 
Wright, G. A., Nicolson, S. W. & Shafir, S. Nutritional physiology and ecology of honey bees. Annu. Rev. Entomol 63, 327–344. https://doi.org/10.1146/annurev-ento-020117-043423 (2018).CAS 
Article 
PubMed 

Google Scholar 
Michener, C. D. The Social Behavior of the Bees: A Comparative Study (Belknap Press of Harvard University Press, 1974).
Google Scholar 
Stabler, D., Paoli, P. P., Nicolson, S. W. & Wright, G. A. Nutrient balancing of the adult worker bumblebee (Bombus terrestris) depends on the dietary source of essential amino acids. J. Exp. Biol. https://doi.org/10.1242/jeb.114249 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Böhme, F., Bischoff, G., Zebitz, C. P. W., Rosenkranz, P. & Wallner, K. Pesticide residue survey of pollen loads collected by honeybees (Apis mellifera) in daily intervals at three agricultural sites in South Germany. PLoS ONE 13, e0199995 (2018).PubMed 
PubMed Central 

Google Scholar 
Schilder, A. M. C., Hanson, E. J. & Hancock, J. F. An integrated approach to disease control in blueberries in Michigan. Acta Hortic. 715, 481–488 (2006).CAS 

Google Scholar 
Wise, J. C., Jenkins, P. E., Poppen, R. V. & Isaacs, R. Activity of broad-spectrum and reduced-risk insecticides on various life stages of cranberry fruitworm (Lepidoptera: Pyralidae) in Highbush Blueberry. J. Econ. Entomol. 103, 1720–1728 (2010).PubMed 

Google Scholar 
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2009).MATH 

Google Scholar  More

Fossil dataWe downloaded all fossil occurrences recorded for the order Scleractinia at the species level from the Paleobiology Database (PBDB – paleobiodb.org; accessed on 3 August 2021). This is the most comprehensive repository for palaeontological data in reef corals to date. Due to the nature of the data, no ethics approval was required. To minimize identification issues, we excluded taxa with uncertain generic and species assignments (i.e., classified as aff. and cf.) and only selected species that had accepted names. We also selected the variables classification and palaeoenvironment from the output options to facilitate taxonomic and environmental filters applied in downstream analyses. The full dataset consisted of 24,011 occurrences across 4235 species, spanning over 250 Myr of coral evolution from the Triassic to the present. Although our focus here lies on the Cenozoic, we used the complete fossil dataset (i.e., including all of the occurrences) to have estimates of the diversification dynamics in scleractinian corals throughout the whole timespan of their evolution.Evolutionary ratesWith the full palaeontological dataset, we estimated evolutionary rates through time in scleractinian corals using the Bayesian framework of the program PyRate (v3.0)12,36,37. This program uses fossil occurrence data to calculate the temporal variation in rates of preservation, speciation and extinction, while incorporating multiple sources of uncertainty12. At its core implementation, PyRate jointly estimates the times of origination (Ts) and extinction (Te) for each fossil lineage; the fossilization and sampling parameters that determine preservation rates (q); and the overall rates of speciation (λ) and extinction (μ) through time36. Recently, the program has been upgraded to include a reversible jump Markov Chain Monte Carlo (rjMCMC) algorithm to estimate diversification rate heterogeneity, which provides more accurate and precise estimates than other commonly used methods12. Therefore, despite the inherent bias of the fossil record (i.e., estimates are conditioned on sampled lineages), PyRate is a robust method to quantify speciation and extinction rates, and their respective temporal shifts, from fossil occurrence data.Extant taxa can also be included in the PyRate framework as long as they are also represented in the fossil record. This is done to extend the fossil geologic ranges to the recent times. Hence, the first step in our analysis was to identify which species in our dataset is still alive at the present. To do this, we matched the accepted species names in the PBDB dataset with those from the extant species dataset of Huang et al.38. Subsequently, we split our dataset into eleven independent subsets, with the goal of keeping each subset with an equal number of species. Each data subset included a random selection of species with their respective occurrences, which was enough to calculate Ts and Te (see below). This was done to avoid convergence issues, given the large size of our dataset and the consequent complexity of the model37. For each of our subsets, we generated fifty replicates by resampling the fossil ages from their temporal ranges to account for the uncertainty associated with the age of occurrences. We then used the maximum-likelihood test in PyRate to compare between three models of fossil preservation12: the homogeneous Poisson process (HPP; q is constant through time); the nonhomogeneous Poisson process (NHPP; q varies throughout the lifespan of a species); and the time-variable Poisson process (TPP; q varies across geological epochs). The latter model (TPP) was selected across all of our data subsets (Supplementary Table 2).After selecting the preservation model, we first focused on assessing the estimates of times of origination and extinction in each data subset, rather than using the full dataset to jointly estimate all parameters at once as in the original implementation of PyRate. This further reduced the complexity of the model and allowed for more precise parameter estimates. For each replicate in all of our data subsets, we approximated the posterior distribution of Ts and Te through a 50 million generation run of the rjMCMC algorithm under the TPP, sampling parameters every 40 thousand iterations. At the end of each run, we discarded 20% of the samples as burn-in and assessed chain convergence through the effective sample sizes of posterior parameter estimates, using the software Tracer39 (v1.7.1).From the results of this first set of models, we extracted the median estimates of Ts and Te across replicates, and we merged the estimates from the eleven independent data subsets. This merged data frame contained estimated times of origination and extinction for all coral lineages within our fossil dataset. We then used this merged Ts and Te data frame as input for another rjMCMC chain to finally estimate overall λ and μ through time, by applying the option -d in PyRate. In this option, Ts and Te for all fossil lineages are given as fixed values and, therefore, are not estimated by the model. The chain for this model was run for 100 million generations, sampling parameters at every 40 thousand iterations. Once again, we excluded 20% of the initial samples as burn-in and checked model convergence using Tracer. Finally, we calculated net diversification rates through time by subtracting the post burn-in samples of μ from λ.To explore the taxonomic idiosyncrasies in the evolutionary rates of reef corals, we selected the most abundant families on present-day coral reefs in terms of the number of colonies per area18 (Acroporidae, Agariciidae, Merulinidae, Mussidae, Pocilloporidae, and Poritidae). Altogether, species within these families account for ~40% of the total extant diversity in Scleractinia. These families also account for most of the occurrences in the PBDB fossil dataset (excluding extinct families, which are generally older and had little temporal overlap with extant ones): Acroporidae (1457 occ. in 165 spp.); Agariciidae (722 occ. in 89 spp.); Merulinidae (2464 occ. in 229 spp.); Mussidae (1146 occ. in 100 spp.); Pocilloporidae (615 occ. in 64 spp.); and Poritidae (1149 occ. in 91 spp.). Therefore, from our full dataset, we selected six independent ones encompassing all species in each of the selected families. We also selected only species that are classified as reef-associated within these families, since we were specifically interested in these environments. This selection had a negligible effect on the size of the individual datasets, given that the vast majority of fossil species within these families are reef-associated. In each family, we followed the same modelling steps described above to estimate μ and λ, and diversity trajectories. However, this time it was not necessary to split the datasets into subsets, given that each family has far less occurrences than the full dataset. We started by comparing models of preservation, which showed the TPP as the best supported for all families (Supplementary Table 3). Then we created fifty replicates by resampling fossil ages to accommodate the uncertainty associated with the time of occurrences. For each replicate, we ran the rjMCMC algorithm for 50 million generations under the TPP model, with a sampling frequency of 40 thousand iterations. We discarded initial 20% of the samples as burn-in, and assessed convergence through Tracer. We then combined all replicates, resampling 100 random samples from each replicate to assess the estimates of μ and λ through time for each family. Finally, we extracted diversity trajectories in each family for all of the replicates by applying the -ltt option in PyRate, which generates a table with estimated range-through diversity at every 0.1 Myr. From these trajectories, we calculated the mean difference in diversity (slope in species per 0.1 Myr) between subsequent time samples backwards from the present (i.e., diversity in time t was subtracted from diversity in time t-1) using the diff function in R (v4.0.3).As an alternative to PyRate, we also calculated the diversity dynamics of reef coral fossils using the R package divDyn40, which combines a range of published methods for quantifying fossil diversification rates. Differently from PyRate, the metrics applied in divDyn require that the fossil occurrences are split into discrete time bins. Therefore, these metrics treat the origination and extinction rates as independent parameters in each bin, while PyRate is designed to detect rate heterogeneity through a continuous time setting12. Our goal here, however, was not to compare models but to assess the robustness of our rate patterns and diversity trajectories using alternative methods. We divided our dataset into one-million-year time bins to have enough temporal resolution for rate calculations. To account for the uncertainty in the assignment of fossil ages, we created 50 binned replicates by sampling the age of each occurrence from a random uniform distribution, with bounds defined by the age ranges provided in the PBDB dataset. We then used the divDyn function to calculate the per capita rates of origination and extinction through time (based on the rate equations by Foote41) for all scleractinians (Supplementary Fig. 5a) and for reef-associated acroporids alone (Supplementary Fig. 5b). We also used the same procedure to generate range-through diversity curves for each of the six families selected previously, to compare with the curves generated by PyRate (Fig. 2a). Although the rate results differed between the PyRate (Fig. 1) and the divDyn (Supplementary Fig. 5) approaches, the general patterns remained unchanged. Rates are more volatile through time in divDyn estimates, with larger confidence intervals, which is expected from the metrics applied in the package12,42. Yet, we found the same peaks in extinction for Scleractinia: at the Cretaceous-Paleogene and Eocene-Oligocene boundaries, and at the Pliocene-Pleistocene (Supplementary Fig. 5a). The recent peak in speciation in Acroporidae was also detected, although less strong (Supplementary Fig. 5b). Despite these slight differences in rate estimates, the diversity curves reconstructed through divDyn (Supplementary Fig. 6) mirrored almost exactly the ones found with PyRate (Fig. 2), demonstrating that the overall macroevolutionary trends described herein (Figs. 1 and 2) are robust to methodological choices.Diversity-dependent modelsTo assess the effects of diversity dependency on the evolution of reef coral lineages, we implemented the Multivariate Birth-Death model (MBD)11 within the PyRate framework. This method was first described as the Multiple Clade Diversity Dependence model (MCDD)19, in which rates of speciation and extinction are modelled as having linear correlations with the diversity trajectories of other clades. At its original implementation, the MCDD was developed to assess the effects of negative interactions, where increasing species diversity in one group can suppress speciation rates and/or promote extinction in itself or in other ecologically similar clades19. However, the model also incorporates the possibility of positive interactions, where increasing diversity in one clade can correlate with enhanced rates of speciation or buffered extinction. Through further model developments43, the MCDD was updated to also include a horseshoe prior44 on the diversity-dependence parameters, which helped controlling for overparameterization and enhanced the power of the model to recover true effects43. More recently, this model took its current form as the MBD11, with the additional possibilities of including environmental correlates and setting exponential, rather than just linear, correlations.We first applied the MBD to estimate the diversity-dependent effects of individual extant coral families (i.e., the ones selected in the previous analysis; see Evolutionary rates) in their combined diversity trajectories. From the rjMCMC model results for individual families, we extracted estimates of Ts and Te in each of the fifty replicates and merged them across families. This merged dataset with fifty replicates of Ts and Te was then used as input for the MBD model, where we set the relative diversity trajectories of each individual family as predictors. We also included three key environmental predictors—paleotemperature, sea level and rate of sea-level change—to assess their influence in overall evolutionary rates. The paleotemperature data was obtained from Westerhold et al.45, and consists of global mean temperature estimates for the last 66 million years, averaged across 0.1 Myr time bins. Eustatic sea-level data was downloaded from Miller et al.46, and contains estimates of sea level for the last 100 million years in comparison to present-day levels, also split in ~0.1 Myr time bins. With this dataset, we calculated the average rate of sea-level change per million years, as measured from the absolute difference between subsequent sea-level values backwards in time (i.e., sea level in time t was subtracted from sea level in time t-1). These environmental factors were rescaled between 0 and 1 to maintain all predictors on the same relative scale.Under our MBD model, the speciation and extinction rates of all families combined could change through time and through correlations with the relative diversity of individual families or environmental factors. The strength and directionality (positive or negative) of the correlations are also jointly estimated for each predictor within the model11. We ran both linear and exponential correlation models (see formulas in Lehtonen et al.11) in each of our fifty replicates for 25 million generations, sampling parameters at every 25 thousand iterations. We then compared the linear and exponential models through the posterior harmonic means of their log likelihoods, which supported the exponential one as having a better fit. From the posterior estimates, we summarized the speciation (Fig. 3c) and extinction (Fig. 3d) correlation parameters (i.e., the strength of the effect) by calculating their median and 95% Highest Posterior Density (HPD) interval across replicates. Finally, we also summarized the effect of families on lineage turnover (Fig. 3e), which we conceptualize as the sum of the effects on speciation and extinction.The MBD model also provides posterior samples of the weight of the correlation parameters, which is estimated through the horseshoe prior11. In essence, this prior is able to reliably distinguish correlation parameters that should be considered noise from those that represent a true signal in the data11. The parameterization of the horseshoe prior contains local and global Bayesian shrinkage parameters44 from which shrinkage weights (w) can be calculated (see formulas in Lehtonen et al.11). These shrinkage weights associated with each correlation parameter in the MBD model vary between 0 and 1, with values closer to 0 representing noise and values closer to 1 representing a true signal. Through simulations, it has been shown that values of w  > 0.5 indicate that the correlation parameter in question significantly differs from the background noise, being the correlation positive or negative43. However, as a conservative way to infer the weight of correlation parameters, here we use a value of w  > 0.7 to detect significance. This value was calculated for each diversity-dependence parameter (speciation, extinction and turnover) from the median values drawn from the model posteriors.The spatial distribution of reef-associated taxa varied considerably throughout the Cenozoic, with biodiversity hotspots moving halfway across the globe47. Therefore, the best way to capture this dynamic biogeographic history in reef corals is by analysing global diversity patterns like we did in our main MBD model. However, to assess the robustness of our diversity-dependent results against the influence of geographic scale and site co-occurrences, we repeated all the modelling steps described above with two data subsets. First, we selected only fossil species that have occurrences in the Indo-Pacific Ocean (i.e., 30°W–180°W) within the six families. Second, we excluded sites in which the Acroporidae did not co-occur with the other families. In each of these data subsets, we calculated diversity trajectories and used them as predictors in a separate MBD model. These models had a merged dataset of Ts and Te of all species included in each case (Indo-Pacific and co-occurrences) as a response variable.Finally, we followed the same modelling procedures described above to investigate the diversity-dependent effects in family pairwise analyses. We applied the MBD model to assess the effects of all other families in each individual family at a time, while also estimating correlations with the key environmental predictors. From the rjMCMC model results for individual families, we extracted the fifty replicates of estimated Ts and Te. Each replicate was then used as input for an MBD run using the relative diversity trajectories of each other individual family as predictors, along with the environmental variables. Once again, we ran 25 million generations of the MBD, with a sampling frequency of 25 thousand, using both linear and exponential correlation models in each age replicate. For all families, we found that the exponential model had a better fit. We then summarized the correlation parameters and the shrinkage weights (Supplementary Fig. 7) derived from the exponential models per family by calculating the median and 95% confidence intervals across replicates.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Brown, R., Derksen, C, & Wang, L. A multi‐data set analysis of variability and change in Arctic spring snow cover extent, 1967–2008. J. Geophys. Res. Atmos. 115(D16), 1–16 (2010).Tan, A., Adam J. C., & Lettenmaier, D. P. Change in spring snowmelt timing in Eurasian Arctic rivers. J. Geophys. Rese. Atmos. 116(D3), 1–12 (2011).St. Jacques, J. M., & Sauchyn, D. J. Increasing winter baseflow and mean annual streamflow from possible permafrost thawing in the Northwest Territories, Canada. Geophys. Res. Lett. 36(1), 1–6 (2009).Liljedahl, A. K. et al. Pan-Arctic ice-wedge degradation in warming permafrost and its influence on tundra hydrology. Nat. Geosci. 9(4), 312–318 (2016).ADS 
CAS 
Article 

Google Scholar 
Farquharson, L. M. et al. Climate change drives widespread and rapid thermokarst development in very cold permafrost in the Canadian High Arctic. Geophys. Res. Lett. 46(12), 6681–6689 (2019).ADS 
Article 

Google Scholar 
Nitze, I. et al. Remote sensing quantifies widespread abundance of permafrost region disturbances across the Arctic and Subarctic. Nat. Commun. 9(1), 1–11 (2018).Article 

Google Scholar 
Lewkowicz, A. G. & Way, R. G. Extremes of summer climate trigger thousands of thermokarst landslides in a High Arctic environment. Nat. Commun. 10(1), 1–11 (2019).CAS 
Article 

Google Scholar 
Jones, M. K. W., Pollard, W. H. & Jones, B. M. Rapid initialization of retrogressive thaw slumps in the Canadian high Arctic and their response to climate and terrain factors. Environ. Res. Lett. 14(5), 055006 (2019).ADS 
Article 

Google Scholar 
Schuur, E. A. G. et al. Climate change and the permafrost carbon feedback. Nature 520(7546), 171–179 (2015).ADS 
CAS 
Article 

Google Scholar 
Sturm, M., Racine, C. & Tape, K. D. Increasing shrub abundance in the Arctic. Nature 411(6837), 546–547 (2001).ADS 
CAS 
Article 

Google Scholar 
Berner, L. T. et al. Summer warming explains widespread but not uniform greening in the Arctic tundra biome. Nat. Commun. 11(1), 1–12 (2020).ADS 
MathSciNet 
Article 

Google Scholar 
Tape, K. D. et al. Range expansion of moose in Arctic Alaska linked to warming and increased shrub habitat. PLoS ONE 11(4), e0152636 (2016).Article 

Google Scholar 
Ward, D. H. et al. Multi-decadal trends in spring arrival of avian migrants to the central Arctic coast of Alaska: Effects of environmental and ecological factors. J. Avian Biol. 47(2), 197–207 (2016).Article 

Google Scholar 
Tape, K. D. et al. Tundra be dammed: beaver colonization of the Arctic. Glob. Change Biol. 24(10), 4478–4488 (2018).ADS 
Article 

Google Scholar 
Whitfield, C. J. et al. Beaver-mediated methane emission: the effects of population growth in Eurasia and the Americas. Ambio 44(1), 7–15 (2015).CAS 
Article 

Google Scholar 
Westbrook, C. J., Cooper, D. J., & Baker, B. W. Beaver dams and overbank floods influence groundwater–surface water interactions of a Rocky Mountain riparian area. Water Resour. Res. 42(6), 1–12 (2006).
Bunn, S. E. & Arthington, A. H. Basic principles and ecological consequences of altered flow regimes for aquatic biodiversity. Environ. Manage 30(4), 492–507 (2002).Article 

Google Scholar 
Naiman, R. J. & Rogers, K. H. Large animals and system-level characteristics in river corridors. Bioscience 47(8), 521–529 (1997).Article 

Google Scholar 
Kemp, P. S. et al. Qualitative and quantitative effects of reintroduced beavers on stream fish. Fish Fish. 13(2), 158–181 (2012).Article 

Google Scholar 
Furbearer Reports, various authors. Alaska Department of Fish & Game, Division of Wildlife Conservation. Juneau, Alaska (1965–2017)Young, S.B., et al. The Environment of the Noatak River basin, Alaska. Center For Northern Studies. Wolcott, VT (1974)Melchior, H. R. ed. Terrestrial Mammals of the Chukchi-Imuruk Area. In Biological Survey of the Bering Land Bridge National Monument: Revised Final Report. Biology and Resource Management Program, Alaska Cooperative Park Studies Unit, University of Alaska (1979)Georgette, S. & Shiedt, A. Whitefish: traditional ecological knowledge and subsistence fishing in the Kotzebue Sound Region, Alaska. No. 290. Alaska Department of Fish and Game, Division of Subsistence (2005).Brubaker, M. et al. Climate change and health effects in Northwest Alaska. Glob. Health Action 4(1), 8445 (2011).MathSciNet 
Article 

Google Scholar 
Rabung, S. & Norton sound bering straight regional planning team. Norton Sound Bering Straight Regional Comprehensive Salmon Plan: Phase II. Alaska Department of Fish & Game, 1–217 (2015)Jones, B. M. et al. Increase in beaver dams controls surface water and thermokarst dynamics in an Arctic tundra region, Baldwin Peninsula, northwestern Alaska. Environ. Res. Lett. 15(7), 075005 (2020).ADS 
Article 

Google Scholar 
Brewer, M. C. The thermal regime of an arctic lake. EOS Trans. Am. Geophys. Union 39(2), 278–284 (1958).Article 

Google Scholar 
Jorgenson, M. T. et al. Resilience and vulnerability of permafrost to climate change. Can. J. For. Res. 40(7), 1219–1236 (2010).Article 

Google Scholar 
Lachenbruch, A. H. et al. Temperatures in permafrost. Temp. Meas. Control Sci. Ind. 1, 791 (1962).
Google Scholar 
Smith, M. W. Microclimatic influences on ground temperatures and permafrost distribution, Mackenzie Delta, Northwest Territories. Can. J. Earth Sci. 12(8), 1421–1438 (1975).ADS 
Article 

Google Scholar 
Langer, M. et al. Rapid degradation of permafrost underneath waterbodies in tundra landscapes—toward a representation of thermokarst in land surface models. J. Geophys. Res. Earth Surf. 121(12), 2446–2470 (2016).ADS 
Article 

Google Scholar 
Jones, B. M. et al. Identification of unrecognized tundra fire events on the north slope of Alaska. J. Geophys. Res. Biogeosci. 118(3), 1334–1344 (2013).Article 

Google Scholar 
Kantner, S. Swallowed by the Great Land: And Other Dispatches from Alaska’s Frontier. Mountaineers Books (2015)Pastick, N. J. et al. Distribution of near-surface permafrost in Alaska: Estimates of present and future conditions. Remote Sens. Environ. 168, 301–315 (2015).ADS 
Article 

Google Scholar 
Obu, J. et al. Northern Hemisphere permafrost map based on TTOP modelling for 2000–2016 at 1 km2 scale. Earth Sci. Rev. 193, 299–316 (2019).ADS 
Article 

Google Scholar 
Willby, N. J. et al. Rewilding wetlands: beaver as agents of within-habitat heterogeneity and the responses of contrasting biota. Philos. Trans. R. Soc. B Biol. Sci. 373(1761), 20170444 (2018).Article 

Google Scholar 
Kivinen, S., Nummi, P. & Kumpula, T. Beaver-induced spatiotemporal patch dynamics affect landscape-level environmental heterogeneity. Environ. Res. Lett. 15(9), 094065 (2020).ADS 
Article 

Google Scholar 
Pollock, M. M. et al. The importance of beaver ponds to coho salmon production in the Stillaguamish River basin, Washington, USA. North Am. J. Fish. Manag. 24(3), 749–760 (2004).Article 

Google Scholar 
Weber, N. et al. Alteration of stream temperature by natural and artificial beaver dams. PLoS ONE 12(5), e0176313 (2017).Article 

Google Scholar 
Nicieza, A. G. & Metcalfe, N. B. Growth compensation in juvenile Atlantic salmon: responses to depressed temperature and food availability. Ecology 78(8), 2385–2400 (1997).Article 

Google Scholar 
Deegan, L. A. et al. Influence of environmental variability on the growth of age-0 and adult Arctic grayling. Trans. Am. Fish. Soc. 128(6), 1163–1175 (1999).Article 

Google Scholar 
Jones, B. M. et al. Multi-dimensional remote sensing analysis documents beaver-induced permafrost degradation, Seward Peninsula, Alaska. Remote Sens. 13(23), 4863 (2021).ADS 
Article 

Google Scholar 
Turetsky, M. R. et al. Permafrost collapse is accelerating carbon release. Nature 569, 32–34 (2019).ADS 
CAS 
Article 

Google Scholar 
Jung, T. S. et al. Colonization of the Beaufort coastal plain by Beaver (Castor canadensis): a response to shrubification of the Tundra?. Can. Field-Naturalist 130(4), 332–335 (2016).Article 

Google Scholar 
Halley, D. J., Saveljev, A. P. & Rosell, F. Population and distribution of beavers Castor fiber and Castor canadensis in Eurasia. Mammal Rev. 51(1), 1–24 (2021).Article 

Google Scholar 
Ecke, F. et al. Meta-analysis of environmental effects of beaver in relation to artificial dams. Environ. Res. Lett. 12(11), 113002 (2017).ADS 
Article 

Google Scholar 
Raynolds, M. K. et al. A raster version of the Circumpolar Arctic Vegetation Map (CAVM). Remote Sens. Environ. 232, 111297 (2019).ADS 
Article 

Google Scholar  More

While rapid urbanization has been an engine of economic growth, it has also resulted in the encroachment on agricultural land and ecologically sensitive areas. In many places, it has also resulted in unequal access to opportunities and the inadequate provision of public services. Furthermore, climate change, biodiversity loss, economic crises, and other global processes are increasing the frequency and severity of disturbances, such as floods, droughts, storms, heatwaves, and sea level rise1, as well as the spread of epidemics and economic volatility. These shocks, in conjunction with unsustainable urban developments, pose challenges for cities, exposing their vulnerabilities and inadequate management paradigms2,3. For example, extreme water shortages in São Paulo have resulted from a confluence of increasing water demand, encroachment on urban watersheds, and deforestation in the Amazon4,5. These shortages have different impacts on urban dwellers in central and peripheral neighborhoods, disproportionately affecting already marginalized communities6. They are also a harbinger of future climate impacts on the city as a whole7.The functioning of cities depends on complex and interdependent social, ecological and technological systems. Citizens access water, electricity, food, health, and other services through complex networks of technological infrastructure and institutional arrangements that obscure the relationship between urban resource provision and the ecosystems that sustain them. The interdependencies in these systems mean that shocks can reverberate broadly, creating unanticipated and cascading crises. Shocks can propagate along supply chains, mobility networks, and across sectors, such as finance and health, and from local to global and back to local levels through non-linear feedback mechanisms8. For example, hurricanes hitting cities on the southern US coastline reveal the interdependencies between extreme weather events exacerbated by climate change and critical infrastructure systems and the sensitive ecological context in which they are situated. New Orleans has repeatedly seen the simultaneous breakdown of water supply, drainage services, and transportation networks, due to electricity system interruptions in response to coastal and riverine flooding or other storm impacts9. Furthermore, the poisoning of soil, water, and air from flooded and damaged oil refineries and other heavy industries disproportionally impacts marginalized communities because of pre-existing housing and insurance regulations10, as well as the region’s ecosystems. Similarly, the Covid-19 pandemic left cities scrambling to respond to an unprecedented health crisis, but also knock-on economic effects that led to shortages of jobs, affordable housing, and food11. These multi-hazard events and cascading crises, and their intersection with systemic racism and other environmental justice issues, highlight the need for coordination across sectors and scales. Furthermore, cities are embedded into international trade networks that deliver goods and services from around the globe12. Decision-making at the local scale therefore has the potential to influence resilience and sustainability both locally and globally13,14,15.Coordination across complex and interdependent urban systems requires adequate forms of governance. Since almost all future population growth is projected to be in urban areas, identifying which governance arrangements allow for inclusive decision-making, management, and planning, while allowing for system-wide transformations, is crucial not only for the well-being of urban dwellers but also for achieving climate goals and maintaining the biosphere16. Despite the urgent need for governance aimed at increasing urban sustainability1, there is limited research on the complexities of governing urban systems across sectors and scales17. Adequate frameworks and guidelines for how to govern such systems could help close the gap between the normative goal of urban sustainability across scales and the need to ensure resilience to shocks with on-the-ground practices18,19.While the concepts of resilience and sustainability are often used interchangeably, it is important to make the distinction between the two3,15,20. Resilience describes the perseverance of a system and its functions, especially its ability to recover from and reorganize in response to disturbances3,21. Resilience can therefore also describe systems in an undesirable regime (e.g., resource-intensive or socially unjust), where overcoming unsustainable practices is impeded due to the system’s resilience. Resilience approaches focus on the process of responding to disturbances, and building adaptive capacity to cope with unknown shocks and stressors22,23. Sustainability, on the other hand, is grounded in the maintenance of the biosphere, which is critical for human well-being and societal development3, and the equitable distribution of resources and services for current and future generations. Sustainability requires balancing societal and ecological goals across scales and over time20,24. It also requires transformative capacity—the ability to create fundamentally different human-environment interactions when a current system becomes untenable25. Sustainability approaches integrate normative values and anticipatory thinking to identify desirable future options and strategies to attain those options24,26,27.Governance refers to the ways in which individuals and organizations, public and private, manage their common affairs28. The process of governance (or the act of ‘governing’) encompasses “the activities of social, political, and administrative actors that can be seen as purposeful efforts to guide, steer, control, or manage the pursuance of public goods”29. The governance process is not limited to public goods and services, but also private goods, common-pool resources, and club goods. Governance of sustainability transformations must combine resilience and sustainability approaches24. This requires identifying and setting sustainability goals, facilitating the process of achieving these goals through proactive planning, learning, unlocking, and overcoming undesirable system states and path-dependencies30, while also anticipating and responding to unknown shocks and stressors22,23. Under certain conditions, socio-political shocks and natural disasters may provide a window of opportunity for transformative change31,32. Proactive governance may be able to harness these opportunities by responding to changing pressures and steering urban transformations towards sustainability31.In this review, we synthesize recent literature on the governance of urban sustainability transformations from a social-ecological-technological systems (SETS) perspective33,34. The conceptualization of cities as urban SETS highlights their complex and interdependent nature, the particular governance challenges this poses, and the importance of identifying governance structures that are able to coordinate across distinct sectors, spatial areas, and time. We bring together literature on sustainable resource management and service provision, social processes and transformative capacity, and governance and power relations among actors. The academic fields that this literature review integrates typically operate in isolation. However, their integration is crucial for enabling a deeper understanding of systemic urban sustainability transformations, as they provide insights into different aspects of the transformation of urban SETS.The remainder of this article is organized according to the elements of Fig. 1. In section 2, we introduce the concept of interdependent urban SETS across scales (Fig. 1a), and the role of governance in determining the evolution of urban SETS. Sections 3-5 represent the three elements in the outer circle of Fig. 1b: In section 3, we review innovations in sustainable resource management and service provision, which are fundamental to urban livelihoods, and require effective coordination across sectors and scales. In section 4, we turn to the role of ‘bottom-up’ experimentation for building sustainable and resilient communities, and the interplay between bottom-up and top-down efforts in initiating and enabling transformation processes. In section 5, we discuss the role of polycentric and multilevel governance arrangements, power-relations among governance actors, and the importance of governance ‘embeddedness’ in sustainability transformations. We summarize the elements of Fig. 1 and how they relate to the sustainable functioning and transformation of urban SETS in section 6. We conclude in section 7.Fig. 1: Governance of Urban Social-Ecological-Technological Systems.a Dependencies across time and space. Shocks and stressors have the potential to impact one or all elements of the SETS across spatial scales, and determine the evolution of the urban system (arrow ‘Temporal Dynamics’). Figure adapted from175. b Embedded Governance: Embedded actors and SETS elements interact to shape the governance of urban SETS, negotiate the interfaces of the SETS elements, and coordinate transformation processes towards sustainability across spatial and temporal scales. We focus on three interacting elements of particular importance: (1) Equitable and reliable access to resources in urban areas (‘Provision of goods and services’). These are an important basis for urban livelihoods and for inclusive social, economic, and political processes to take place. (2) The interactions of bottom-up processes and economic activities with top-down efforts by government actors (‘Bottom-up and top-down processes’). A balance between ‘top-down’ and ‘bottom-up’ dynamics can enable local innovation and experimentation to result in systemic transformations. (3) ‘Power-relations among actors’ determine the choice of policies, the distribution of resources and social influence, and thus, the trajectories of urban SETS. Arrows in (a) and (b) represent the embeddedness of governance.Full size imageUrban social-ecological-technological systems across scalesOur definition of ‘urban’ includes formal and informal settlements within a metropolitan area, as well as the footprints that an urban area produces through its supply chains and emissions. Thus, the spatial and institutional boundaries of ‘urban’ vary depending on the particular governance issue being considered. For example, the planning of bicycle lanes in a city has different boundaries than the issue of urban food security, which includes regional and global trade networks.Considering cities as SETS, and governance as embedded within these SETS, foregrounds the interdependencies, tradeoffs, and feedbacks across system elements at different scales and the implications of these dynamics for sustainability transformations16,34,35, as depicted in Fig. 1. Ecosystems include forests, oceans, rivers, lakes, groundwater, soils, and the atmosphere, which provide resources, including water, food, fuel, fiber, and other ecosystem services. The ecosystem dimension of SETS also includes the space or land that support urban economies and on which cities are built. The social system represents relevant actors and their interactions, including state authorities and agencies, individuals and civil society organizations, academics, business actors, and others, who gather information and make decisions over how to access and distribute resources, manage ecosystems, organize SETS in space, and consume or utilize the goods and services provided by ecosystems. The interactions among actors are mediated not only by norms, rules, culture, and power dynamics, but also by their ecological and technological environment. The technological system, or physical infrastructure, includes buildings, facilities, and infrastructure networks, which transport flows of resources, people, and information, and which mediate the interaction between ecosystems and social actors.One example of social-ecological-technological interactions is flood protection infrastructure (technological), which is planned and designed by decision-makers and built to protect urban populations (social) from storm surges, riverine flooding, and rising sea levels in coastal and riparian areas (ecological). Another example of SETS interactions that highlights the broad spatial dependencies are food supply chains. These include soils, water, plants, animals, and energy from the sun (ecological system), infrastructure deployed for planting and harvesting, for processing and storage, and for transportation and sale (technological), as well as the people involved along the supply chains and the citizens who depend on the delivery of food into urban areas (social system). Sustainable governance of these systems must ensure the inclusivity of the design and planning process, as well as the equitable and ecologically viable production and distribution of these goods and services.Governance, as defined here, determines the evolution of the SETS as a whole: how people live together and how they live within their social-ecological-technological environment, or, to take the specific examples from above: what types of flood protection measures are chosen, who is involved in the decision-making process, what type of food is grown where, for and by whom, using what technologies. Governance tools extend beyond formal laws and regulations to include informal and collective action mechanisms, such as access to reliable information, effective communication, monitoring and sanctioning of rule-breakers, and processes for resolving conflicts36,37. They also include nudges or messaging to promote more sustainable behaviors, the creation of ‘enabling environments’ for bottom-up processes and social movements, the building of adaptive capacities, and participatory mechanisms of co-production and management24,38,39,40,41,42,43, among others.Governance is also highly context- and path-dependent44,45, which is reflected in the diversity of sociopolitical and governance regimes around the world, the way in which cities are planned and built, the inclusiveness of the political process, and the distribution of resources, rights and services to citizens. That said, quantitative research has also identified regularities in how cities evolve across scales. This body of work has identified a set of scaling relations that can be used to predict the social, spatial and infrastructural properties of cities of different size and at different levels of development, and to inform urban planning strategies46,47,48.A sustainable governance system is able to anticipate and buffer temporal dynamics that arise from internal change processes or external pressures, and can adapt the urban system to ensure inclusive processes in determining how goods and services, such as transportation and health49, are provided. When critical links within the governance system and essential mechanisms of governance are missing, or when shocks are unanticipated, they can lead to catastrophic and cascading effects. For example, adequate anticipation, proactive and inclusive planning, and learning could make coastal cities, such as New Orleans, more prepared to respond to repeated hurricane threats. Tighter cross-scale and cross-level feedback loops, for example, between different cities, countries, and sectors, as well as anticipatory planning could have avoided some of the detrimental impacts of the Covid-19 pandemic and measures introduced to mitigate its spread50.Furthermore, efforts that do not take a holistic, cross-scale, and cross-level view but endeavor to transform one part of the system, may meet resistance, tradeoffs or unexpected responses from other system elements and scales35. For example, decisions by governments and private actors to reduce greenhouse gas (GHG) emissions has led to remarkable technological developments along the entire supply chain of production. As a case in point, in 2015, Apple partnered with its suppliers in China to install two gigawatts of renewable energy, reducing its supply-chain GHG emissions51. This shift in energy supply may entail other externalities. For example, the mineral resources required for a renewable energy transition are unevenly spread across the globe and require intensive, disruptive, and often exploitative mining practices52. The company, while improving its carbon footprint, has also been criticized for massive human rights abuses, including labor forces drawn from so-called ‘re-education camps’ for Uyghur people53. So, while a more conscious (urban) customer base is pushing powerful actors to address climate pressures through technological innovations, the abuse of labor rights and the existence of low labor standards remain unaddressed due to the disconnect between urban consumption and the production processes occurring in distant locations.The imbalance across different domains and dimensions of sustainability goals can be exacerbated by policies or incentives that are narrowly focused on one of the SETS elements rather than simultaneously considering sustainability across social, ecological, and technological systems. Therefore, sustainability transformation goals and processes should consider all of the SETS elements, including how they are embedded across sectors and scales.Sustainable urban resource management and service provisionHome to the majority of the global population, cities are the largest consumers of natural resources, and the largest producers of CO2 emissions and other environmental pollutants54. Thus, the way in which cities manage their resource demands and flows is decisive for the sustainability of urban areas, including the global ecosystems they rely on for the delivery of their goods and services. Conventional urban systems operate linearly: resources are extracted, used, and released as waste into the environment55. For example, cities access energy through extensive networks56; extract industrial and construction materials, such as metals and sand57,58; import food from across the globe59; reach hundreds of miles to extract water from distant rivers, lakes, and groundwater bodies60; and release pollution, such as domestic and industrial wastewater (with limited treatment) back into the environment61,62. These linear urban supply systems lack the necessary feedback between ecosystems and consumers to signal when ecosystems are being degraded or overexploited63,64.In contrast, cross-sectoral approaches focus on the synergies between the inputs and outputs of different sectors65. Systems of reuse and recycling66,67 have been deployed to address the problem of scarcity, such as implemented for Windhoek’s water supply system (Namibia)68, and in Singapore’s NEWater recycling system, both of which redistribute treated wastewater to customers69. Networked systems of water, energy, and food (WEF) can produce mutual benefits from their coordinated management17,70. WEF nexus approaches link the wastewater sector with agriculture through treated wastewater used in irrigation and recovered nutrients for the fertilization of crops71. Wastewater-energy co-benefits occur, for example, through the methane released from sewage sludge, which can be used for electricity production, such as implemented in Berlin (Germany), Melbourne (Australia)72,73 and urban China74. Food-energy co-benefits include the biodiesel produced from used cooking oil75 by several companies in North Carolina (USA)76. These cross-sector systems require intense coordination among social actors, ecosystems, and technologies16,35. For example, water suppliers and wastewater managers need to coordinate to ensure that water brought to the city and turned into wastewater can be adequately treated and discharged, and that the nutrients and biosolids extracted from sewage sludge adhere to quality standards to avoid soil and crop contamination (e.g., heavy metals contained in wastewater can contaminate agricultural soils, rivers, and groundwater).The circular economy (CE) and related concepts, such as urban metabolism, industrial ecology, cradle-to-cradle, material flows and life-cycle analysis, are recurring themes in the urban sustainability literature that stand in contrast to conventional, linear systems17,77,78. These CE approaches are gaining attention in both Global North and Global South cities79,80,81. Narrowly defined, the CE describes systems that reduce, recover, reuse, and recycle resources and materials across sectors19, similar to the networked WEF nexus. Other interpretations of the CE include well-being and inclusiveness as goals, and the need for a radical transformation of society, including shifts in mindsets and norms, policy-making, and business models, in order to decouple economic growth from resource consumption18. Safe and just operating space models15,82,83 combine ecological limits (i.e., the Planetary Boundaries84) with societal well-being. These approaches recognize that climate change and degrading ecosystems have disproportionate impacts on the urban poor, and, like other shocks, have the potential to increase existing inequalities85. Thus, building transformative capacity for climate adaptation planning and implementation must address existing inequities, unmet needs and the adaptive capacity of the urban poor86.Different social and political contexts take different approaches to addressing sustainability challenges. For example, China and Europe are leading the way regarding the adoption of circular city concepts, albeit with different drivers (i.e., governmental or ‘top-down’ in China vs. societal or ‘bottom-up’ in Europe) and implementation strategies (i.e., technological innovation vs. socio-technological and socio-ecological approaches)67. In China, the introduction of laws and financial incentives related to the CE and the resulting implementation of streamlined technological processes increased resource efficiency by 35% and reduced waste intensity by 47% between 2005-201387. European strategies have developed from increasing material efficiency, similar to those implemented in China and elsewhere87, to broader and more diverse efforts initiated at the sub-city and other levels of governance. These include urban living lab initiatives, social and industrial start-ups promoting and implementing CE approaches in various domains, including urban agri-food, textiles, plastics, recycling of building materials, water, mobility, and others, supported by civil society organizations that provide advice about the CE67. European CE efforts and related discourse thus extend beyond material flows to include the development of social, inclusive and collaborative values, green identities, sharing economies, and community gardens and spaces for creative exchange88. Such initiatives involve actors across levels, including businesses, non-governmental organizations (NGOs), state-level actors, and civil society, and, depending on the specific goals of resource flow efficiencies or community development88, are characterized by different levels of integration of SETS elements.Different approaches to the CE also reflect distinct origins, with some approaches emerging from the natural/engineering sciences and others from the social sciences. There is an ongoing convergence of CE concepts across different disciplines. Thus, much of the recent academic literature adopts holistic and transformational concepts of circularity, which include material flows, equity, justice, inclusion, and a balanced distribution of power in their definitions18. In contrast, governments and practitioners often define CE in technocratic terms. In both interpretations, the ability of the CE to achieve sustainability goals faces several fundamental challenges, such as the decoupling of economic development from consumption; rebound effects resulting from increased resource efficiency;89 and resource tradeoffs in the bioeconomy, such as the use of biofuels, which in some cases can be more environmentally damaging than fossil fuels18,90.Many urban sustainability efforts focus on local improvements in ecosystem protection, which are often achieved through the outsourcing of production and polluting industrial processes to other regions or countries, such as the case for feed crop production and heavy industry59,91, or focus on only one aspect of sustainability (e.g. GHG emissions) at the expense of others, as in the example of Apple discussed in Section 2. However, goods also flow into urban areas as indirect resources. Water is embedded in imported energy, and water and energy are embedded in imported food14,92, resulting in large ecological and water footprints15 and climate impacts. CE and WEF nexus approaches can lead to certain efficiency improvements, and recent, high-level policy developments, such as the European Green Deal and the US Green New Deal, contain elements of the CE, with implications for urban resource management. However, many practical challenges remain, including social and technological barriers and a dearth of legal and political frameworks to guide implementation of the CE19,93. Furthermore, reducing the footprints of consumption will ultimately require behavioral changes, such as shifts from meat- to plant-based diets94,95,96. These will in turn require carefully and holistically designed policy and governance interventions, such as the realignment of subsidies, price floors, trade regulations, and pricing of externalities to favor plant-based diets, as well as efforts to emphasize the link between food production and consumption97. Governance mechanisms that prioritize inclusive policy and decision-making arrangements can bridge bottom-up behavioral changes with top-down incentive structures and regulations. Broader definitions of sustainability goals that account for cross-scale and cross-sector interactions are needed to avoid externalities resulting from outsourcing and problem-shifting.Bottom-up and top-down processes in urban sustainability transformationsDifferent governance contexts will favor technocratic or social innovations, and will determine whether these are implemented through top-down control, emerge from bottom-up initiatives98, or result from a mix of the two. Urban SETS planning influences not only the efficiency of land and resource use, CO2 emissions, public health outcomes99, and the resilience of cities to climate impacts100; it also has major implications for socio-political processes, such as the emergence of local interest groups and social movements, and the equity and inclusiveness of urban development101. Both spatial planning and large-scale technology implementation in the urban context, often top-down processes, can enable or suppress transformative forces in urban settings. Throughout history, this has been demonstrated by the important function of physical urban spaces as places for congregations and the voicing of public demands, but also for the suppression and segregation of social groups101. For example, the restructuring of Paris after the French Revolution aimed to make it easier to control the urban population and quench potential uprisings101. Similarly, after the demonstrations and political upheavals during the ‘Arab Spring’, Cairo’s Tahrir square went from being the main square for public congregation to a camera- and police-patrolled, traffic-regulated space in which gatherings are prohibited. Online media platforms (new technological systems) have also had a critical role in the emergence of socio-political movements, but also the control of political opinions102,103. This has been demonstrated by the role of the Facebook and Twitter platforms in the mobilization of civil society in the Arab Spring movements and global Black-Lives-Matter protests, as well as right-wing supporters of former US president Trump. The role of space (‘ecological systems’) and technology in mediating social relations is critical for understanding the emergence of bottom-up processes, as well as the interaction between bottom-up and top-down processes, in overcoming the status quo and driving sustainability transformations.Lock-in is often related to the durability of urban infrastructure systems and the legacy effects of the built environment, including buildings, facilities, roads, and water and electricity supply networks; but it can also characterize entrenched political systems, formal and informal institutions, and social discourse104. This observation has prompted renewed attention to the question of which social dynamics can lead to transformative change, and the relationship between social tipping and ecological Earth system processes40,105. Social tipping points describe abrupt shifts in social norms or patterns of behavior40. The tendency for social conformity106 means that, under certain circumstances, localized policy interventions, messaging by formal institutions38, or bottom-up experimentation and innovations can trigger widespread changes in values or behaviors40,105,107. However, it must be recognized that social conformity can also work to undermine collective benefits, as was demonstrated in the fascist era during the first half of the twentieth century.One focus of the urban sustainability transformations literature is on identifying enabling environments for local experimentation as a potential way to overcome path-dependencies and lock-ins that rigidify the status-quo108. This ‘urban tinkering’ approach emphasizes the role of bottom-up, distributed change and the need to engage a range of social actors in the continuously unfolding process of urban SETS transformations. This perspective of the urban transformations and sustainability science literature24 has parallels with adaptive ecological systems, which are characterized by diversity, open exchange, modular organization and redundancy109—features that can increase adaptive capacity. Bottom-up initiatives generate these characteristics by building on the creativity and diverse knowledge of unique sets of local actors working together to develop context-specific approaches and solutions. These local initiatives and institutions create modular structures, which can operate separately, but work best when they are interlinked to a moderate degree to facilitate coordination and exchange110. Diversity, redundancy and modularity contribute to general resilience by enhancing buffering capacity; that is, by dampening the spread of disturbances across modules111,112.Different capacities are needed during distinct phases of an urban transformation, including, among others, the ability to dismantle drivers of unsustainable practices (‘unlocking’), and the ability to create innovations that contribute to sustainability and resilience, and to embed these innovations in governance structures, practices and discourse (‘transforming’)113. Furthermore, social learning and reflexive stakeholder action are important to ensure that transformative processes remain adaptive and aligned with evolving sustainability goals114. Social or collective learning is particularly important for managing the multidimensional contingencies of urban SETS, which require the integration of knowledge across scales and sectors, and are difficult for any single actor or entity to grasp115.Enabling conditions and capacities for local governance at the urban scale vary vastly, even under seemingly similar democratic regimes. For example, cities in India, Brazil, and South Africa, the three largest democracies in the Global South, had strikingly similar spatial and social forms in the 1980’s, marked by high levels of inequality and state-sanctioned differentiated citizenship with unequal access to formal housing and basic services116. In the 1990’s, all three countries implemented reforms aimed at decentralization, which shifted power from the central state to the municipal level, in order to promote more participatory forms of governance116. The degree of state-civil society embeddedness and different dynamics between top-down and bottom-up processes have led to significantly diverging developmental outcomes. India remains largely elite-dominated and controlled by the central state with very weak local capacity for governance and autonomy, and unequal and often poor outcomes for basic service provision, including access to housing. South Africa is characterized by strong capacities of local government, but remains technocratic, centrally controlled and fragmented by design of the apartheid regime, including only certain segments of the population in governance processes. In contrast, Brazil is characterized by strong, autonomous local capacity and broad-based, participatory governance processes. While not without challenges, and with significant heterogeneity across cities, more inclusive and empowering policies, such as new legal instruments that regularize ‘squatted’ land in informal settlements, have enabled the extension of infrastructure services and the improvement of living conditions in historically neglected neighborhoods116.Certain environments or junctures in time may create enabling conditions for the emergence of small-scale experimentation or interventions. These may, incrementally or abruptly, alter established paths and act as building blocks towards transformative change40,117,118. The existence of public spaces as meeting points, including community gardens, public parks and squares, such as the Tahrir Square prior to its closure to the public, and central Paris until its transformation after the French Revolution, may stimulate sociopolitical debate creating new narratives and ideas for change, and, through a sense of shared concern or interest, promote the emergence of community organizations and a local public sphere119,120. Resulting efforts, locally adapted, socially inclusive, decentralized initiatives can play a key role in urban sustainability transformations121.The COVID-19 pandemic has shown that rapid and large scale reorganization of daily life is possible, at least in the short-term. Many governments worldwide enforced social distancing measures, and closed their borders and ‘non-essential’ parts of their economies to mitigate the spread of the virus. These were enforced as top-down measures. While economic activities and associated carbon emissions are quickly returning to normal, such shocks can temporarily destabilize entrenched practices and narratives, potentially creating a critical juncture or ‘window of opportunity’ for institutional reform and systems change32,122. For example, COVID-19 drew heightened attention to many pre-existing social inequities and the inadequacy of existing social welfare, as well as the lack of preparedness by many urban areas to address a crisis on this scale. This, as well as compounding pressures from extreme weather events, have led to calls for improved healthcare, social safety nets, and broader sustainability transformations123,124. At the same time, the pandemic itself has also created growing inequality within countries and globally125. Unfortunately, systemic shocks can result in increased inequality, unless specific governance structures and mechanisms for reducing inequality are in place before and during such ‘windows of opportunity’85. Whether disruptions such as the pandemic become critical junctures for positive systemic transformations or instead accelerate existing unsustainable trends depends critically on proactive governance responses, and the structures that allow or prevent those responses from effecting wholistic change8,31,126 and avoiding maladaptation31,127. In contrast, transformational change can also be continuous, rather than spontaneous and disruptive, but requires a transformation of the governance system itself128. This relates back to the discussion of capacities and enabling conditions for sustainability transformations discussed above. Rather than transformed through an outside force, such as an epidemic, or through technical interventions, transformational change from within the governance system is achieved through deliberative acts that focus on changing the governance process that enables emergence, adaptation, and implementation of new paradigms, i.e., how things are done, rather than what is done128.Theories about the enabling conditions for urban transformation often depart from the on-the-ground reality of urban governance. Sectors are often managed by siloed entities, disconnected across sectors and spatial scales113; local sustainability initiatives (bottom-up) are often fragmented, lack funding, and rely on engaged individuals, which creates obstacles for lasting impacts and system-wide transformations119; and government initiatives (top-down) are typically short-term and responsive rather than long-term and mitigative, anticipatory or preparatory113. In the following section, we discuss the importance of governance structures that are able to create bridges across sectors, spatial scales, and institutional levels.Power relations and embeddedness of governance in urban SETSPower relations and embeddedness can be critical factors in determining whether bottom-up and top-down initiatives are able to transform urban SETS116,119,129. In this section, we focus on the relationship between the embeddedness of governance, power relations among actors and the sustainability of urban systems. The scholarly literature describes embeddedness as the interdependence and relations among societal actors, and links the concept of embeddedness with those of capacity and agency in governance processes116,130. Embeddedness forms the basis for coordination, conflict resolution and addressing collective action problems. The degree of embeddedness and the connectivity of actor relationships across sectors and levels of governance can be crucial for ensuring an inclusive political process in urban settings116.The relational quality of embeddedness is tightly connected to the relational concept of power130. Relational power is understood to be incorporated in social practices, embedded in social-ecological-technological relations, and productive, meaning that it “shapes, creates and transforms social relationships, practices and institutional arrangements”130,131. Recent scholarly work integrating complex adaptive systems thinking with the governance of sustainability transitions recognizes the need for a concept of power-relations that is based not only in the interactions among social agents, but also in SETS interactions130. We use the term embeddedness to refer to the relations among governance actors across levels and sectors and their relationship with ecological and technological elements in the urban system.Governance can be structured in diverse ways. Polycentric (PG) and multi-level governance (MLG) structures have received substantial attention in climate governance discourses29,37,110. Both PG and MLG are used to describe the interdependent, multi-level and multi-scalar nature of governance and the overlapping decision-making authority of relevant actors. The focus of these approaches on the interdependence between stakeholders at different levels and scales is particularly important for the conception of cities as dynamic SETS. Locally embedded, self-organizing structures are emphasized in PG. Bottom-up, modular PG systems with cross-level and cross-unit linkages can benefit sustainable and equitable resource management by increasing institutional fit, and enabling mechanisms of mutual monitoring, learning, accountability and conflict resolution, and increasing adaptive capacity37,110. Thus, PG is particularly important in the context of resilience and sustainability of urban SETS due to its focus on interactions across units of governance and on local contexts, which enables anticipation, learning and adaptation. MLG accounts for the increasingly international context of policy and administration, and emphasizes the redistribution of power and control away from the state. This redistribution occurs in three directions: (1) upwards to international actors and organizations, (2) downwards to cities, regions, and communities, and (3) outwards to civil society and non-state actors29. The process of globalization contributes to the redistribution of power, a change in the focal scale of governance from local to global, and can lead to a disconnect between policy design and implementation and the local context. Thus, the relationship between MLG and PG is especially relevant for urban governance, which must take into consideration both local and global sustainability issues.Governance of urban sustainability transformations is often associated with democratic, inclusive, multi-level, and multi-scale characteristics132. The degree of embeddedness varies across contexts. In Europe local initiatives and authorities collaborate with state, business, non-governmental, and academic actors and organizations at the local, sub-national, national, and European scale88, ensuring a certain extent of governance embeddedness at horizontal and vertical scales. In the federal system of the U.S., states are important actors in the design of environmental policies, some of which can be implemented and enforced at the local level. States with large enough markets can also shape industries and drive national standards, as in the case of California leading automobile emissions standards, as well as pioneering air pollution and plastic waste policies133,134. In India, cities have limited capacity for local self-governance in part due to a lack of cross-level embeddedness between the central state, municipalities and civil society and a highly skewed distribution of power concentrated in the central state116,135.Embedding governance vertically, from within the city up to the regional and global levels, and horizontally, across sectors, cities and regions, can help ensure that local sustainability improvements are not achieved at the cost of global sustainability15,136. Furthermore, voluntary coalitions—such as international city networks like ICLEI, C40, or the World Mayors Council on Climate Change—enable the sharing of information, the spread of innovations, and strengthen governance capacity and the potential impact of bottom-up or decentralized efforts137. City networks have existed since the early 1900s, but have proliferated in the past 30 years into a broad ecosystem of over 200 global partnerships between local authorities138,139. These range from publicly to privately initiated and from public membership networks to ones that include multiple types of stakeholders. In recent years, they have begun to play a pivotal role in defining and making progress towards global sustainability agendas, in part due to the increasing acknowledgment that decisions at the urban scale drive many global change processes140,141. Theory based on model analyses also predicts that these voluntary coalitions will be more successful in addressing climate change and other global sustainability challenges than efforts to achieve global consensus among national players142,143. These networked urban governance systems support sustainability initiatives through mutual efforts towards agreed upon goals, such as emissions reductions, but also by facilitating knowledge co-creation and dissemination through the systematic organization of learning across locales, and by offering institutional support144,145. In general, decentralized initiatives that are embedded in broader coalitions can generate a diversity of responses to sustainability challenges, allowing those that are most successful to spread to other locations41,119,146. However, disconnects between civil society initiatives and the local state or other municipalities can limit the potential for the diffusion of innovation across locations119.Shifts in power and decision-making authority from public to non-public actors147, typical of multilevel governance settings148, can both strengthen and undermine the agency of citizens and cities to influence sustainability149,150. For instance, critical infrastructure systems around the world have been transferred from public authorities to private actors. This has often compromised social inclusion, for example, by limiting access to healthcare for lower income groups in Europe and South America151,152,153. Similarly, the management of local infrastructure services, such as water supply, by transnational corporations has increased social inequities and negative ecological impacts154,155,156,157. Thus, while diverse public and private governance actors have the ability to usher in sustainable change51,158, shifting decision-making authority can also result in the concentration of power in private, political and economic elites at different scales. This can lead to policy capture and vested interests, which entrench high-carbon trajectories, instead of enabling transformations through participation and inclusion24,129. Strengthening governance embeddedness between the state, civil society and private entities can help maintain the balance of power among actors that is needed for sustainable urban development and transformations116.While governance interventions may aim to address climate risk, the sustainable management of natural resources, and social inclusion, their effectiveness depends on the vertical embeddedness of local governments with civil society and with higher level authorities at sub-national, national and international levels. For example, local actors in the Aburra Valley (Colombia) created an inter-city risk response network to address the increasing vulnerability of growing numbers of urban dwellers in areas at high risk of flooding and landslides. This horizontally and vertically embedded PG structure includes civil society, local and national government, international city networks, and international organizations and requires high levels of trust among these actors159.Governance inevitably requires attention to multiple issues and constituencies, and the prioritization of certain goals and values over others160. Inclusive, deliberative decision-making processes, such as multi-criteria decision approaches, can be used to explore synergies and tradeoffs between technological, social and ecological objectives161. A critique of such decision-making processes is that addressing climate change requires urgency and pragmatism. In addition, accounting for multiple interest groups increases the complexity of decision-making, increasing the time-, resource- and data-intensity as well as delays due to contested (political) interests that must be addressed by more transparent processes161. However, attention to power dynamics and measures that increase accountability and protect vulnerable populations are crucial. Participatory governance approaches aim to redistribute power by equitably involving diverse actors in the planning, decision-making and implementation of policies162. For such inclusive systems of governance to be successful, the incorporation of a range of judicial, popular, and political safeguards are needed to ensure lasting cooperation among actors and avoid opportunism and power imbalances163. For example, New York City’s water supply system relies on a set of institutions that ensure the balance of power among diverse actors, including local civil society, municipal, state, and federal government actors. Within this institutional framework the involved actors must continuously negotiate their various interests163.The sustainability of urban SETS requires the social actors involved, including state and non-state actors, to assume critical roles of governance164. Central to the role of the state and other actors in a regulatory or management position is the ability to manage uncertainty, maintain the rule of law and order, coordinate policies, and redistribute wealth, resources and opportunities165. To ensure that this ‘state power’ remains within certain limits, non-state actors and civil society need to monitor and check the state and elites through political and public action, which creates feedbacks between civil society and the state116,164. When this state-society embeddedness is lost, power concentrated at the national or state level can lead to hollow local states, even in democratic contexts116. This can result in (and be reinforced by) inequitable public services, including water, transport, and public health, and the degradation of natural resources and the built environment, despite significant economic growth116,166. Indian cities, such as Chennai, are typical examples, in which only half of the urban population has access to water on their premises, soil and water pollution and resource decline are accelerating, leaving taps dry and rivers contaminated167. In US cities, too, such as East Chicago (IL) and Flint (MI), out-of-balance power relations and the lack of response to feedback between SETS elements have exacerbated social inequalities168. In Flint (MI), sustained lead poisoning resulted from corroding lead water pipes after the change of the water source. One the other hand, efforts that promote sustainability locally, such as investments in public transit, bike- and pedestrian-friendly neighborhoods, and urban green space, can drive up property values and lead to gentrification through private property investment. This can, again, have major implications for low-income residents169,170,171. In such cases, state regulation (top-down) can help maintain the social cohesion of neighborhoods undergoing sustainability transformations. The urban governance cycle shown in Fig. 1b illustrates the two-way links between the equitable provision of urban services, top-down and bottom-up processes, and the power relations of the governance system.The provisioning of critical infrastructure services produces and maintains human well-being. It also conditions the ability of the state to provide legitimate order, facilitate effective economic development, and promote social inclusion165. In contrast, mistrust in the state and a rejection of state power and government can impede sustainability efforts51. In extreme cases, fragmentation of state-civil society contact can lead to the emergence of competing sovereignties, which can weaken the legitimacy of a central state172. For example, in many of the informal settlements housing the urban poor in Mexico City, a flourishing informal sector has replaced the state in providing sources of income, as well as basic infrastructure services, such as water and electricity. The state, focused on rapid macroeconomic growth rather than socio-political inclusion, tolerated the rise of the informal sector, whose leaders gained support and power by providing the services that the state failed to provide. Accelerated by increased violence and coercion emanating from the informal sector, the loss of legitimacy and erosion of the formal state has created a parallel rule by drug cartels172. Once eroded, regaining an inclusive governance system is a steep, uphill struggle. In contexts where trust, consent and reciprocity have been eroded, and ‘relational tipping points’ have been reached, governance must balance the urgency of ecological sustainability transformations with the slow process of rebuilding social relations173. These are crucial for sustainability efforts that are inclusive of marginalized populations, such as indigenous people173.In all of these examples, it is the governance context, its embeddedness with SETS interactions, and power dynamics, which provide the enabling (or disabling) environment for sustainability transformations to take place from local to global scales.SETS transformations and the urban governance cycleThe growing gap between the normative goals and responsibilities of urban governance systems described in sustainability narratives and on-the-ground practice highlights the need for a more systematic understanding of how urban governance can facilitate coordination and transformation of urban SETS. In this last section, we discuss how the literature included in this review contributes to understanding the interplay of sustainable resource management and service provision, bottom-up processes and top-down interventions, and power relations among actors in the transformation pathways of urban SETS.Sustainable provision of urban services. The way in which social actors employ technologies to interact with ecosystems and to extract and allocate resources are the basis for the provision of basic services, including water, energy, food, health, transportation, etc. As discussed in section 3, much of the empirical literature on the sustainable provision of services and resource management in linear, circular and nexus approaches is focused on achieving either environmental or social goals. Combining these goals requires addressing the pressures and stressors on natural ecosystems, including climate change and the degradation of global ecosystems, as well as the need for adequate access and just distribution of resources and services. Sustainable service provision is therefore not merely a question of circular material flows and the protection of ecological systems, but must also address social equity issues and cross-scale impacts on social, ecological, and technological system dynamics174. The multiple, entangled resource inflows and outflows of the urban metabolism, as well as various stakeholder interests, make defining and achieving sustainability goals a political challenge for SETS governance.Bottom-up and top-down processes. Identifying promising paths for sustainability transformations requires processes of experimentation in the social, ecological and technological spheres. Bottom-up processes are particularly relevant in this context as they generate diverse and locally adapted responses. Embedded governance allows successful experimentation to be transmitted through social learning dynamics. As such, these bottom-up initiatives and voluntary coalitions can be crucial for overcoming lock-in and breaking path-dependencies. Local initiatives can engender widespread behavioral shifts and norm dynamics, which can in turn spread along horizontally interlinked city networks. Furthermore, as modular building blocks of polycentric governance structures, local, bottom-up initiatives that are interlinked across scales also provide system-wide redundancies and buffering capacities that are crucial during periods of change, including the recovery from shocks and transformation processes.Governmental and other forms of top-down support are crucial for ensuring that adequate spaces for social exchange are maintained, and for promoting and upscaling promising local innovations that may help the transition to sustainability. Governments must also develop and uphold regulations to constrain power imbalances among actors and ensure the separation of legislative, judicial and executive powers. The built environment and (digital) technologies interact with these social processes, both through bottom-up movements and through top-down control. Widespread sustainability transformations can be triggered from the bottom-up if they are embedded in a SETS governance context and supported by structures that are built through bottom-up, top-down and multilevel interactions.Power Relations among actors. Urban SETS are complex, cross-boundary, multi-scale systems, whose governance is best described as multilevel. The diverse ecology of actors in multilevel governance creates power dynamics that can easily become unbalanced if they are not safeguarded and protected in an anticipatory manner. The authority of cities in determining how to best adapt or transform in response to the diverse pressures resulting from internal change and external shocks depends on the distribution of power and the degree of embeddedness into several levels of vertical governance, including international, national, and sub-national levels. Within cities, vertical embeddedness and power relations describe the inclusivity of decision-making processes, which affects the degree of trust, consent, and the wellbeing of residents. How SETS interactions, service provision, and social processes evolve thus depends on the embeddedness of governance across the SETS elements and from local to global levels. More

Study sites, sample collection and preparationAdult female mosquitoes used in this study had previously been collected from three areas: Kerio Valley (Baringo county), Rabai (Kilifi county) and Nguruman (Kajiado county) (Fig. 1), as part of vector-borne disease surveillance project and stored at – 80 °C at the International Centre of Insect Physiology and Ecology (icipe). The mosquitoes were surveyed between August 2019 and May 2020. Nguruman is an agropastoral area located in Kajiado county at the southern end of the Kenyan Rift Valley bordering Tanzania. The area has a semi-arid climate characterized by erratic rains, extreme temperatures, and cyclic and prolonged droughts30. The vegetation is dominated by bushland, grassland and open woodlands along seasonal river valleys. Specific indicator data for malaria is not available for Nguruman except for estimates pertaining to the larger Kajiado county which as of 2019 indicates a malaria incidence rate of 5 per 1000 population31. Collections in Kerio Valley (Baringo county within the Rift Valley) were conducted in Kapluk and Barwesa, both agro-pastoral areas with arid and semi-arid ecology. Malaria is a major vector-borne disease in the areas with report of perennially occurrence in neighboring riverine areas32. Rabai is one of the seven administrative sub-counties of Kilifi county in the coastal region of Kenya where malaria is endemic. The main economic activities in the area include subsistence agriculture, casual labor, crafts and petty trading. The weather patterns at the sites during the sampling period were as follows: Kerio Valley (mean daily temperature: 21.2 °C, mean daily rainfall: 4.1 mm, mean relative humidity: 73.4%); Rabai (mean daily temperature: 26.4 °C, mean daily rainfall: 2.1 mm; mean relative humidity: 78.1%) and Nguruman (mean daily temperature: 22.5 °C, mean daily rainfall: 0.9 mm, mean relative humidity: 61.2%).Mosquito survey and processingHost seeking mosquitoes were trapped using CDC light traps baited with dry ice (carbon dioxide) attractive to several mosquitoes. Traps were set outdoors about 10–15 m away from randomly selected homesteads from 18:00 h to 06:00 h. After collection, the mosquitoes were anesthetized with trimethylamine and temporarily stored in liquid nitrogen before transportation to the Emerging Infectious Disease (EID) laboratory at icipe and later stored at − 80 °C. Anopheline mosquitoes were morphologically identified to species level using published taxonomic keys15,33.DNA extraction and Anopheles species discriminationDNA was extracted from the head/thorax of individual mosquitoes using ISOLATE II Genomic DNA Extraction kit (Bioline, UK) following the manufacturer’s instructions and used for species discrimination and screening for P. falciparum infection and Gste2 mutations (described below).Cryptic sibling species of the Anopheles funestus and Anopheles gambiae complexes were identified using conventional PCR34,35 and/or sequencing. PCR for An. funestus complex in a 15 µl reaction volume comprised 0.5 µM of each primer targeting: Anopheles funestus s.s, Anopheles vaneedeni, Anopheles rivulorum, Anopheles parensis, Anopheles leesoni, Anopheles longipalpis A and Anopheles longipalpis C, 3 µl of 5X HOT FIREPol Blend Master Mix Ready to Load (Solis BioDyne, Estonia) and 2 µl of DNA template. The cycling conditions were initial denaturation at 95 °C for 15 min, and then 30 cycles of denaturation at 95 °C for 30 s, annealing at 46 °C for 30 s and extension at 72 °C for 40 s and final extension at 72 °C for 10 min. Size fragments of each species were scored after separation in 1.5% agarose gel electrophoresis stained with ethidium bromide against a 1 Kb DNA ladder (HyperLadder, Bioline, London, UK).For An. gambiae s.l., PCR in a 10 µl volume consisted of 2 µl of 5X Evagreen HRM Master Mix (Solis BioDyne, Estonia), 1 µl of DNA template and 10 µM concentration of each primer targeting An. gambiae s.s and An. arabiensis. The thermal cycling conditions included initial denaturation for 15 min at 95 °C followed by 40 cycles of denaturation at 95 °C for 20 s, annealing at 61 °C for 15 s and extension at 72 °C for 20 s followed by final extension at 72 °C for 7 min.A subset of An. funestus s.l. samples that failed to amplify using the established protocol, was further amplified and sequenced targeting the internal transcribed spacer 2 (ITS2) region of the ribosomal DNA (rDNA)36. This target has shown utility in discriminating closely related mosquito species including anophelines12 and sequences from diverse species for this marker are well represented in reference databases (e.g. GenBank). PCR volumes for rDNA ITS2 were 15 µl containing 0.5 µM of the forward and reverse primers, 3 µl of 5X HOT FIREPol Blend Master Mix Ready to Load (Solis BioDyne, Estonia) and 2 µl of DNA template. The cycling conditions were initial denaturation at 95 °C for 15 min, followed by 40 cycles of denaturation at 95 °C for 30 s, annealing at 60 °C for 30 s and extension at 72 °C for 45 s and final extension at 72 °C for 7 min. ExoSAP IT rapid cleanup kit (Affymetrix Inc., Santa Clara, CA, USA) was used to clean the PCR product as per the manufacturer’s guideline, and then outsourced for bidirectional Sanger sequencing to Macrogen, South Korea.Detection of malaria parasitesPlasmodium falciparum sporozoites in individual mosquitoes (head/thorax) were detected by analyzing high resolution melting (HRM) profiles generated from real time PCR products of non-coding mitochondrial sequence (ncMS)37. A P. falciparum DNA from National Institute for Biological Standards and Control (NIBSC; London, UK) was used as a reference positive control. PCR was carried out in a 10 µl volume consisting of 2 µl of 5X Evagreen HRM Master Mix (Solis BioDyne, Estonia), 1 µl of DNA template and 10 µM of each primer. PCR cycling conditions were initial denaturation for 15 min at 95 °C followed by 40 cycles of denaturation at 95 °C for 20 s, annealing at 61 °C for 15 s and extension at 72 °C for 20 s followed by final extension at 72 °C for 7 min. A fraction of RT-PCR-HRM positive samples were further analyzed using conventional PCR in a 10 µl volume consisting of 2 µl of 5X HOT FIREPol Blend Master Mix Ready to Load (Solis BioDyne, Estonia), 1 µl of DNA template and 10 µM of each primer. The cycling conditions comprised initial denaturation for 15 min at 95 °C followed by 40 cycles of denaturation at 95 °C for 20 s, annealing at 61 °C for 15 s and extension at 72 °C for 20 s followed by final extension at 72 °C for 7 min. PCR product of samples positive by RT-PCR were purified using ExoSAP- IT (USB Corporation, Cleveland, OH, USA) and outsourced for sequencing to Macrogen, South Korea. All sporozoite-positive mosquitoes were molecularly identified to species by PCR of the ITS2 region as described above.Genotyping for L119F-GSTe2 mutation and sequencingTwo outer and two inner primers in a PCR assay were used to genotype the L119F-GSTe2 mutations that confer resistance of An. funestus mosquitoes to pyrethroids/DDT19 as described previously28. Thus, only An. funestus s.l. was screened using this assay. Briefly, PCR in a 15 µl reaction volume consisted of 10 µM of each primer, 3 µl of 5X HOT FIREPol Blend Master Mix Ready to Load (Solis BioDyne, Estonia), and 2 µl of DNA template. The cycling conditions were initial denaturation at 95 °C for 15 min, followed by 40 cycles of denaturation at 95 °C for 30 s, annealing at 59 °C for 30 s and extension at 72 °C for 40 s and final extension at 72 °C for 7 min. Amplicons were resolved in a 1.5% agarose gel stained with ethidium bromide (Sigma-Aldrich, GmbH, Germany) against a 1 Kb DNA ladder (HyperLadder, Bioline, London, UK). The amplicons were scored as either homozygous susceptible (SS) at 312 bp, homozygous resistant (RR) at 523 bp or heterozygous (RS) when both bands were visualized.Representative GSTe2 allele positive samples were sequenced for the GSTe2 gene using the Gste2F and Gste2R primers as described previously38. PCR comprised a reaction volume of 15 µl in MyTaq DNA Polymerase Kit (Bioline, London, UK) containing 10 µM of each primer, 5X My Taq reaction buffer, 2 µl of My taq DNA polymerase and 1 µl of DNA template. PCR conditions were: initial denaturation of 5 min at 95 °C, followed by 30 cycles of 94 °C for 30 s, 58 °C for 30 s and 72 °C for 1 min, with a final extension at 72 °C for 10 min. Cleaning and sequencing of amplicons were performed as described above.Sequence and polymorphism analysisSequences (mosquito, P. falciparum, GSTe2) were viewed and cleaned in Geneious Prime39 and queried in GenBank using Basic Local Alignment Search Tool (BLastn). Parasite sequences were assigned as P. falciparum after  > 98% percentage identity. MAFFT in Geneious Prime39 was used to perform multiple sequence alignments with default parameters. Maximum likelihood (ML) trees were inferred for mosquito ITS2 sequences using the best fit model of sequence evolution with nodal support for different groupings evaluated through 1000 bootstrap replications. GSTe2 gene polymorphism analysis was performed in Geneious Prime and ML tree reconstructed from MAFFT alignment using PhyML v. 2.2.4. Haplotype distribution network was constructed using Templeton-Crandall Sing (TCS) program v. 1.2140.Statistical analysisRelative abundance was used to estimate the outdoor composition of the anopheline mosquitoes. Daily counts of female mosquito/trap/night for An. funestus s.l. and An. gambiae s.l. were compared for each area using generalized linear models (GLM) with negative binomial error structure based on best-fit model residuals. The mean catches/trap/night was computed for each of the species complexes. The P. falciparum sporozoite infection rates (Pfsp) were expressed as the number of positive specimens of the total number of specimens examined. The distribution of L119F-GSTe2 mutations was assessed by determining allelic frequencies in different species. Infection status among the resistant mosquitoes was compared using the Fisher’s Exact Test. Data were analyzed using R v 4.1.0 software at 95% confidence limit.Ethical considerationsEthical review and approval of the study was granted by the Scientific and Ethical Review Unit (SERU) of the Kenya Medical Research Institute (KEMRI) (Protocol No. SSC 2787). Prior to data collection, the purpose of the study, procedures and associated benefits/risks were provided to the local leadership at county and community levels. Additionally, informed verbal consent to trap mosquitoes around homesteads was obtained from household heads. More

Cooke, S. J. et al. Knowledge co-production: A pathway to effective fisheries management, conservation, and governance. Fisheries 46, 89–97 (2021).Article 

Google Scholar 
Kobluk, H. M. et al. Indigenous knowledge of key ecological processes confers resilience to a small-scale kelp fishery. People Nat. 3, 723–739 (2021).Article 

Google Scholar 
Lee, L. C. et al. Drawing on indigenous governance and stewardship to build resilient coastal fisheries: People and abalone along Canada’s northwest coast. Mar. Policy 109, 103701 (2019).Article 

Google Scholar 
Reid, A. J. et al. “Two-Eyed Seeing”: An Indigenous framework to transform fisheries research and management. Fish. Fish. 22, 243–261 (2021).Article 

Google Scholar 
Toniello, G., Lepofsky, D., Lertzman-Lepofsky, G., Salomon, A. K. & Rowell, K. 11,500 y of human–clam relationships provide long-term context for intertidal management in the Salish Sea, British Columbia. Proc. Natl Acad. Sci. 116, 22106–22114 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ahn, J. E. & Ronan, A. D. Development of a model to assess coastal ecosystem health using oysters as the indicator species. Estuar., Coast. Shelf Sci. 233, 106528 (2020).CAS 
Article 

Google Scholar 
Skilbeck, C. G., Heap, A. D. & Woodroffe, C. D. Geology and sedimentary history of modern estuaries. in Applications of Paleoenvironmental Techniques in Estuarine Studies (eds. Weckström, K., Saunders, K. M., Gell, P. A. & Skilbeck, C. G.) 45–74 (Springer Netherlands, 2017). https://doi.org/10.1007/978-94-024-0990-1_3.Durham, S. R., Gillikin, D. P., Goodwin, D. H. & Dietl, G. P. Rapid determination of oyster lifespans and growth rates using LA-ICP-MS line scans of shell Mg/Ca ratios. Palaeogeogr., Palaeoclimatol., Palaeoecol. 485, 201–209 (2017).Article 

Google Scholar 
Lockwood, R. & Mann, R. A conservation palaeobiological perspective on Chesapeake Bay oysters. Philos. Trans. R. Soc. B 374, 20190209 (2019).CAS 
Article 

Google Scholar 
Rick, T. C. et al. Millennial-scale sustainability of the Chesapeake Bay native American oyster fishery. Proc. Natl Acad. Sci. 113, 6568–6573 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Thompson, V. D. et al. Ecosystem stability and Native American oyster harvesting along the Atlantic Coast of the United States. Sci. Adv. 6, eaba9652 (2020).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zimmt, J. B., Lockwood, R., Andrus, C. F. T. & Herbert, G. S. Sclerochronological basis for growth band counting: A reliable technique for life-span determination of Crassostrea virginica from the mid-Atlantic United States. Palaeogeogr. Palaeoclimatol. Palaeoecol. 516, 54–63 (2019).Article 

Google Scholar 
Alleway, H. K. & Connell, S. D. Loss of an ecological baseline through the eradication of oyster reefs from coastal ecosystems and human memory. Conserv Biol. 29, 795–804 (2015).PubMed 
Article 

Google Scholar 
Beck, M. W. et al. Oyster reefs at risk and recommendations for conservation, restoration, and management. Bioscience 61, 107–116 (2011).Article 

Google Scholar 
Kirby, M. X. Fishing down the coast: Historical expansion and collapse of oyster fisheries along continental margins. Proc. Natl Acad. Sci. 101, 13096 (2004).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lotze, H. K. et al. Depletion, degradation, and recovery potential of estuaries and coastal seas. Science 312, 1806 (2006).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Zu Ermgassen, P. S. et al. Historical ecology with real numbers: Past and present extent and biomass of an imperilled estuarine habitat. Proc. R. Soc. B: Biol. Sci. 279, 3393–3400 (2012).Article 

Google Scholar 
Carranza, A., Defeo, O. & Beck, M. Diversity, conservation status and threats to native oysters (Ostreidae) around the Atlantic and Caribbean coasts of South America. Aquat. Conserv.: Mar. Freshw. Ecosyst. 19, 344–353 (2009).Article 

Google Scholar 
Pluckhahn, T. J. & Thompson, V. D. Woodland-period mound building as historical tradition: Dating the mounds and monuments at Crystal River (8CI1). J. Archaeological Sci.: Rep. 15, 73–94 (2017).Article 

Google Scholar 
Waselkov, G. A. Shellfish gathering and shell midden archaeology. Adv. Archaeol. Method Theory 10, 93–210 (1987).Article 

Google Scholar 
McNiven, I. J. Ritualized middening practices. J. Archaeol. Method Theory 20, 552–587 (2013).Article 

Google Scholar 
Hawkes, A. D. et al. Relative sea-level change in northeastern Florida (USA) during the last ~8.0 ka. Quat. Sci. Rev. 142, 90–101 (2016).ADS 
Article 

Google Scholar 
Kelley, J. T., Belknap, D. F. & Claesson, S. Drowned coastal deposits with associated archaeological remains from a sea-level “slowstand”: Northwestern Gulf of Maine, USA. Geology 38, 695–698 (2010).ADS 
Article 

Google Scholar 
Khan, N. S. et al. Drivers of Holocene sea-level change in the Caribbean. Quat. Sci. Rev. 155, 13–36 (2017).ADS 
Article 

Google Scholar 
Love, R. et al. The contribution of glacial isostatic adjustment to projections of sea-level change along the Atlantic and Gulf coasts of North America. Earth’s Future 4, 440–464 (2016).ADS 
Article 

Google Scholar 
Shugar, D. H. et al. Post-glacial sea-level change along the Pacific coast of North America. Quat. Sci. Rev. 97, 170–192 (2014).ADS 
Article 

Google Scholar 
Dougherty, A. J. et al. Redating the earliest evidence of the mid-Holocene relative sea-level highstand in Australia and implications for global sea-level rise. PLoS ONE. 14, e0218430 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bailey, G. N. The role of molluscs in coastal economies: The results of midden analysis in Australia. J. Archaeol. Sci. 2, 45–62 (1975).Article 

Google Scholar 
Habu, J., Matsui, A., Yamamoto, N. & Kanno, T. Shell midden archaeology in Japan: Aquatic food acquisition and long-term change in the Jomon culture. Quat. Int. 239, 19–27 (2011).Article 

Google Scholar 
Hale, J. C. et al. Submerged landscapes, marine transgression and underwater shell middens: Comparative analysis of site formation and taphonomy in Europe and North America. Quat. Sci. Rev. 258, 106867 (2021).Article 

Google Scholar 
Erlandson, J. M. et al. Shellfish, geophytes, and sedentism on Early Holocene Santa Rosa Island, Alta California, USA. J. Isl. Coast. Archaeol. 15, 504–524 (2020).Article 

Google Scholar 
Rick, T. C. Early to Middle Holocene estuarine shellfish collecting on the islands and mainland coast of the Santa Barbara Channel, California, USA. Open Quaternary 6, 9 (2020).Sanger, D. & Sanger, M. J. Boom and bust on the river: The story of the Damariscotta oyster shell heaps. Archaeol. East. North Am. 14, 65–78 (1986).
Google Scholar 
Moss, M. L. Shellfish gender, and status on the Northwest Coast: Reconciling archaeological, ethnographic, and ethnohistoric records of the Tlingit. Am. Anthropologist 95, 631–652 (1993).Article 

Google Scholar 
Cannon, A., Burchell, M. & Bathurst, R. Trends and strategies in shellfish gathering on the Pacific Northwest Coast of North America. in Early Human Impact on Megamolluscs (eds. Antczak, A. & Cipriani, R.) 7–22 (Archaeopress, 2008).Grier, C., Angelbeck, B. & McLay, E. Terraforming and monumentality as long-term social practice in the Salish Sea region of the Northwest Coast of North America. Hunt. Gatherer Res. 3, 107–132 (2017).Article 

Google Scholar 
Pluckhahn, T. J. & Thompson, V. D. New Histories of Village Life at Crystal River. (University Press of Florida, 2018).Thompson, V. D. et al. Ancient engineering of fish capture and storage in southwest Florida. Proc. Natl Acad. Sci. 117, 8374–8381 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sassaman, K. E. Complex hunter–gatherers in evolution and history: A North American perspective. J. Archaeol. Res. 12, 227–280 (2004).Article 

Google Scholar 
Luby, E. M. & Gruber, M. F. The dead must be fed: Symbolic meanings of the shellmounds of the San Francisco Bay area. Camb. Archaeol. J. 9, 95–108 (1999).Article 

Google Scholar 
Lightfoot, K. G. & Luby, E. M. Mound building by California hunter-gatherers. in The Oxford Handbook of North American Archaeology (ed. Pauketat, T.) 212–223 (Oxford University Press, 2012).Smith, A. D. T. Archaeological expressions of Holocene cultural and environmental change in coastal Southeast Queensland. (The University of Queensland, 2016).Reeder-Myers, L., Rick, T., Lowery, D., Wah, J. & Henkes, G. Human ecology and coastal foraging at Fishing Bay, Maryland, USA. J. Ethnobiol. 36, 595–616 (2016).Article 

Google Scholar 
Petrie, C. C. Tom Petrie’s reminiscences of Early Queensland (dating from 1837). (Watson, Ferguson & Company, 1904).Eipper, C. Statement of the Origin, Condition and Prospects, of the German Mission to the Aborigines at Moreton Bay, etc. (James Reading, 1841).Watkins, G. Notes on the Aboriginals of Stradbroke and Moreton Islands. Proc. R. Soc. Qld. 8, 40–50 (1891).
Google Scholar 
Ross, A. & with members of the Quandamooka Aboriginal Land Council. Aboriginal approaches to cultural heritage management: A Quandamooka case study. in Australian Archaeology ’95: Proceedings of the 1995 Australian Archaeological Association Annual Conference (eds. Ulm, S., Lilley, I. & Ross, A.) vol. Tempus 6 107–112 (Anthropology Museum, University of Queensland, 1996).Jenkins, J. A. & Gallivan, M. D. Shell on earth: Oyster harvesting, consumption, and deposition practices in the Powhatan Chesapeake. J. Isl. Coast. Archaeol. 15, 384–406 (2020).Article 

Google Scholar 
Hatch, M. B. A. & Wyllie-Echeverria, S. Historic distribution of Ostrea lurida (Olympia oyster) in the San Juan Archipelago. Wash. State Tribal Coll. Univ. Res. J. 1, 38–45 (2016).
Google Scholar 
Swanton, J. R. Social Organization and Social Usages of the Indians of the Creek Confederacy. (Bureau of American Ethnology, 1928).Hening, W. W. The Statutes at Large of Virginia. (1809).Wharton, J. The Bounty of the Chesapeake: Fishing in Colonial Virginia. (Virginia 350th Anniversary Celebration Corporation, 1957).Denys, N. Description géographique et historique des Costes de l’Amérique Septentrionale. Avec l’Histoire naturelle du Pais. (Chez Claude Barbin, 1672).Nicolar, J. The Life and Traditions of the Red Man. (Duke University Press, 2007 Print, 1893).Speck, F. G. Penobscot Man: The Life History of a Forest Tribe in Maine. (University of Pennsylvania Press, 1940).Washburn, K. Passamaquoddy tribe conducts oyster project. Bangor Daily News (1979).Kennedy, V. S. Shifting Baselines in the Chesapeake Bay: An Environmental History. (Johns Hopkins University Press, 2018).de Charlevoix, P. F. X. Journal of a Voyage to North America, Vollume II. Translated by Louise Phelps Kellogg. (The Caxton Club, 1923).Ingersoll, E. The Oyster Industry. (United States Bureau of Fisheries, United States Census Office, Government Printing Office, 1881).Brice, J. J. Report on the fish and fisheries of the coastal waters of Florida. in Report of the Commissioner for the Year Ending June 30, 1896 263–242 (U.S. Commission of Fish and Fisheries, U.S. Government Printing Office, 1896).Blake, B. & Zu Ermgassen, P. S. E. The history and decline of Ostrea lurida in Willapa Bay, Washington. J. Shellfish Res. 34, 273–280 (2015).Article 

Google Scholar 
Thurstan, R. H. et al. Charting two centuries of transformation in a coastal social-ecological system: A mixed methods approach. Global Environmental Change 61, 102058 (2020).Schulte, D. M. History of the Virginia oyster fishery, Chesapeake Bay, USA. Front. Mar. Sci. 4, 127 (2017).Fletcher, M.-S., Hamilton, R., Dressler, W. & Palmer, L. Indigenous knowledge and the shackles of wilderness. Proc. Natl Acad. Sci. 118, e2022218118 (2021).Ross, A., Coghill, S. & Coghill, B. Discarding the evidence: The place of natural resources stewardship in the creation of the Peel Island Lazaret Midden, Moreton Bay, southeast Queensland. Quat. Int. 385, 177–190 (2015).Article 

Google Scholar 
Reeder-Myers, L. A. & Rick, T. C. Kayak surveys in estuarine environments: addressing sea level rise and climate change. Antiquity 93, 1040–1051 (2019).Article 

Google Scholar 
Savarese, M., Walker, K. J., Stingu, S., Marquardt, W. H. & Thompson, V. The effects of shellfish harvesting by aboriginal inhabitants of Southwest Florida (USA) on productivity of the eastern oyster: Implications for estuarine management and restoration. Anthropocene 16, 28–41 (2016).Article 

Google Scholar 
Lulewicz, I. H., Thompson, V. D., Cramb, J. & Tucker, B. Oyster paleoecology and native American subsistence practices on Ossabaw Island, Georgia, USA. J. Archaeol. Sci.: Rep. 15, 282–289 (2017).
Google Scholar 
Hesterberg, S. G. et al. Prehistoric baseline reveals substantial decline of oyster reef condition in a Gulf of Mexico conservation priority area. Biol. Lett. 16, 20190865 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Cannarozzi, N. R. & Kowalewski, M. Seasonal oyster harvesting recorded in a Late Archaic period shell ring. PloS ONE. 14, e0224666 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Cook-Patton, S. C., Weller, D., Rick, T. C. & Parker, J. D. Ancient experiments: Forest biodiversity and soil nutrients enhanced by Native American middens. Landsc. Ecol. 29, 979–987 (2014).Article 

Google Scholar 
Stalter, R. & Kincaid, D. The vascular flora of five Florida shell middens. J. Torre. Botanical Soc. 131, 93–103 (2004).Article 

Google Scholar 
Kirby, M. X. & Miller, H. M. Response of a benthic suspension feeder (Crassostrea virginica Gmelin) to three centuries of anthropogenic eutrophication in Chesapeake Bay. Estuar. Coast. Shelf Sci. 62, 679–689 (2005).ADS 
Article 

Google Scholar 
Suttles, W. Variation in habitat and culture on the Northwest Coast. in Coastal Salish Essays 26–44 (University of Washington Press, 1987).Bliege Bird, R. & Nimmo, D. Restore the lost ecological functions of people. Nat. Ecol. Evolution 2, 1050–1052 (2018).Article 

Google Scholar 
Berkes, F. Indigenous ways of knowing and the study of environmental change. J. R. Soc. N.Z. 39, 151–156 (2009).Article 

Google Scholar 
Tengö, M., Malmer, P., Elmqvist, T. & Brondizio, E. S. A Framework for Connecting Indigenous, Local and Scientific Knowledge Systems. (2012).Ellis, E. C. et al. People have shaped most of terrestrial nature for at least 12,000 years. Proc. Natl Acad. Sci.118, e2023483118 (2021).Roberts, P. et al. Reimagining the relationship between Gondwanan forests and Aboriginal land management in Australia’s “Wet Tropics”. Iscience 24, 102190 (2021).Ogburn, D. M., White, I. & McPhee, D. P. The disappearance of oyster reefs from eastern Australian estuaries—impact of colonial settlement or mudworm invasion? Coast. Manag. 35, 271–287 (2007).Article 

Google Scholar 
Diggles, B. K. Historical epidemiology indicates water quality decline drives loss of oyster (Saccostrea glomerata) reefs in Moreton Bay, Australia. N.Z. J. Mar. Freshw. Res. 47, 561–581 (2013).CAS 
Article 

Google Scholar 
Pritchard, C., Shanks, A., Rimler, R., Oates, M. & Rumrill, S. The Olympia oyster Ostrea lurida: Recent advances in natural history, ecology, and restoration. J. Shellfish Res. 34, 259–271 (2015).Article 

Google Scholar 
Trimble, A. C., Ruesink, J. L. & Dumbauld, B. R. Factors preventing the recovery of a historically overexploited shellfish species, Ostrea lurida Carpenter 1864. J. Shellfish Res. 28, 97–106 (2009).Article 

Google Scholar 
White, J., Ruesink, J. L. & Trimble, A. C. The nearly forgotten oyster: Ostrea lurida Carpenter 1864 (Olympia oyster) history and management in Washington State. J. Shellfish Res. 28, 43–49 (2009).Article 

Google Scholar 
Harding, J. M., Spero, H. J., Mann, R., Herbert, G. S. & Sliko, J. L. Reconstructing early 17th century estuarine drought conditions from Jamestown oysters. Proc. Natl Acad. Sci. 107, 10549–10554 (2010).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mann, R., Harding, J. M. & Southworth, M. J. Reconstructing pre-colonial oyster demographics in the Chesapeake Bay, USA. Estuar., Coast. Shelf Sci. 85, 217–222 (2009).ADS 
Article 

Google Scholar 
Bayne, B. L. Biology of Oysters. (Elsevier Science & Technology, 2017).Galtsoff, P. S. The American Oyster Crassostrea virginica Gmelin. (United States Government Printing Office, 1964).Kennedy, V. S., Newell, R. I. E. & Eble, A. F. The Eastern Oyster: Crassostrea virginica. (University of Maryland Sea Grant Publications, 1996).Grabowski, J. H., Powers, S. P., Peterson, C. H., Gaskill, D. & Summerson, H. C. Growth and survivorship of non-native (Crassostrea gigas and Crassostrea ariakensis) versus native eastern (Crassostrea virginica) oysters. J. Shellfish Res. 23, 781–793 (2004).
Google Scholar 
Shumway, S. Natural environmental factors. in The eastern oyster Crassostrea virginica (eds. Kennedy, V., Newell, R. & Eble, A.) 467–513 (Maryland Sea Grant, 1996).Lyman, R. L. Paleoenvironmental reconstruction from faunal remains: Ecological basics and analytical assumptions. J. Archaeol. Res. 25, 315–371 (2017).MathSciNet 
Article 

Google Scholar 
Claasen, C. Shells. (Cambridge University Press, 1990).Giovas, C. M. The shell game: Analytic problems in archaeological mollusc quantification. J. Archaeol. Sci. 36, 1557–1564 (2009).Article 

Google Scholar 
Peltier, W. R., Argus, D. F. & Drummond, R. Space geodesy constrains ice age terminal deglaciation: The global ICE-6G_C (VM5a) model: Global Glacial Isostatic Adjustment. J. Geophys. Res.: Solid Earth 120, 450–487 (2015).ADS 
Article 

Google Scholar  More

Argonema gen. nov. Skoupý et Dvořák.Type species: Argonema galeatum.Morphology: Filamentous cyanobacterium, colonies macroscopic, growing in round bulbs and tufts. The filaments are dark green to blue-green, grey-green or brown-green in color. Cells are wider than they are long. Filaments sheathed, sheaths are colorless to light brown, distinct, and variable in length. The filament can protrude from the sheath or the sheath can exceed filament. Trichomes are cylindrical, not attenuated to slightly attenuated towards the end, slightly or not constricted at cell walls. The apical cell can be concave, dark brown, purple-brown to almost black. Cell content often granulated. Necridic cells present, reproduction by hormogonia. The morphological description was based on both culture and fresh material.Etymology: The genus epithet (Argonema) is derived from greek Argo – slow, latent (αργός) and nema – thread (νήμα).A. galeatum sp. nov. Skoupý et Dvořák.Morphology: The cells of A. galeatum are 6.5–9.1 µm (mean 7.81 µm) wide and 1.1–2.5 µm (mean 1.83 µm) long (Figs. 1–5). Filaments are straight, blue-green to gray-green in color. The sheaths are colorless to light brown, distinct, and variable in length. The filament can protrude from the sheath or the sheath can exceed filament. No true branching was observed. Trichomes are cylindrical, not attenuated or slightly attenuated towards the end, slightly or not constricted at cell walls. Some filaments have a concave apical cell that is dark brown, purple-brown to almost black (Fig. 11b). Cell content often granulated. Reproduction by necridic cells and subsequent breaking of the filaments into hormogonia (Fig. 11a,c). The morphological description was based on both culture and fresh material.Figures 1-8Microphotographs of Argonema galeatum (Figs 1–5) and Argonema antarcticum (Figs. 6–8) Trichomes of A. galeatum appear more straight (Fig 2), while trichomes of A. antarcticum form waves (Fig 6) and loops (Fig 7). Scale = 10 µm, wide arrow = necridic cells, arrowhead = granules, asterisk = colored apical cell, circle = empty sheath.Full size imageFigures 9 and 10Histograms of cell dimensions constructed using PAST software. Fig. 9 – Histogram of cell width frequencies in A. galeatum (blue) and A. antarcticum (red). Fig. 10 – Histogram of cell length frequencies in A. galeatum (blue) and A. antarcticum (red).Full size imageHolotype: 38,057, Herbarium of the Department of Botany (OL), Palacký University Olomouc, Czech Republic.Reference strain: Argonema galeatum A003/A1.Type locality: James Ross Island, Western Antarctica, 63.80589S, 57.92147 W.Habitat: Well-developed soil crust.Etymology: Species epithet A. galeatum was derived from latin galea – helmet.A. antarcticum sp. nov. Skoupý et Dvořák.Morphology: The cells are 7.6–9.2 µm (mean 8.52 µm) wide and 1.2–2.8 µm (mean 1.72 µm) long (Figs. 5–8). Filaments are wavy, gray-green to brown-green in color. The sheaths are colorless to light brown, distinct, and variable in length. The filament can protrude from the sheath or the sheath can exceed filament. No true branching was observed. Trichomes are cylindrical, not attenuated or slightly attenuated towards the end with a concave apical cell, slightly or not constricted at cell walls (Fig. 11d). Necridic cells present (Fig. 11e), reproduction by hormogonia. The morphological description was based on both culture and fresh material.Holotype: 38,058, Herbarium of the Department of Botany (OL), Palacký University, Olomouc, Czech Republic.Reference strain: Argonema antarcticum A004/B2.Type locality: James Ross Island, Western Antarctica, 63.89762S, 57.79743 W.Habitat: Well-developed soil crust.Etymology: Species epithet A. antarcticum was derived from the original sampling site.Morphological variabilityWe used light microscopy to assess the morphology of Argonema from soil crust samples and cultured strains. Argonema is morphologically similar to other Oscillatoriales, such as Lyngbya, Phormidium, and Oscillatoria. In culture, the morphology of A. galeatum and A. antarcticum differed slightly. Filaments of A. antarcticum are wider than cells of A. galeatum, averaging at 8.52 µm (A. galeatum – 7.81 µm). The average cell width/length ratio is 4.54 for A.galeatum and 4.89 for A. antarcticum. The cell width was significantly different between the two species (Nested ANOVA, p  More

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

Google Scholar 
Jones, K. E. et al. Global trends in emerging infectious diseases. Nature 451, 990–993 (2008).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Keesing, F. et al. Impacts of biodiversity on the emergence and transmission of infectious diseases. Nature 468, 647–652 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Carlson, C. J. et al. Climate change will drive novel cross-species viral transmission. Preprint at bioRxiv https://doi.org/10.1101/2020.01.24.918755 (2020).Gibb, R. et al. Zoonotic host diversity increases in human-dominated ecosystems. Nature https://doi.org/10.1038/s41586-020-2562-8 (2020).Loh, E. H. et al. Targeting transmission pathways for emerging zoonotic disease surveillance and control. Vector Borne Zoonotic Dis. 15, 432–437 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Hassell, J. M., Begon, M., Ward, M. J. & Fèvre, E. M. Urbanization and disease emergence: dynamics at the wildlife–livestock–human interface. Trends Ecol. Evol. 32, 55–67 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Cohen, J. M., Sauer, E. L., Santiago, O., Spencer, S. & Rohr, J. R. Divergent impacts of warming weather on wildlife disease risk across climates. Science 370, eabb1702 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Murray, M. H. et al. City sicker? A meta-analysis of wildlife health and urbanization. Front. Ecol. Environ. 17, 575–583 (2019).Article 

Google Scholar 
Becker, D. J., Hall, R. J., Forbes, K. M., Plowright, R. K. & Altizer, S. Anthropogenic resource subsidies and host–parasite dynamics in wildlife. Phil. Trans. R. Soc. B 373, 20170086 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Werner, C. S. & Nunn, C. L. Effect of urban habitat use on parasitism in mammals: a meta-analysis. Proc. Biol. Sci. 287, 20200397 (2020).PubMed 
PubMed Central 

Google Scholar 
Becker, D. J., Streicker, D. G. & Altizer, S. Linking anthropogenic resources to wildlife–pathogen dynamics: a review and meta-analysis. Ecol. Lett. 18, 483–495 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Becker, D. J. et al. Macroimmunology: the drivers and consequences of spatial patterns in wildlife immune defense. J. Anim. Ecol. 89, 972–995 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Albery, G. F. & Becker, D. J. Fast-lived hosts and zoonotic risk. Trends Parasitol. https://doi.org/10.1016/j.pt.2020.10.012 (2021).Seto, K. C., Güneralp, B. & Hutyra, L. R. Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools. Proc. Natl Acad. Sci. USA 109, 16083–16088 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Chen, G. et al. Global projections of future urban land expansion under shared socioeconomic pathways. Nat. Commun. 11, 537 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gao, J. & O’Neill, B. C. Mapping global urban land for the twenty-first century with data-driven simulations and shared socioeconomic pathways. Nat. Commun. 11, 2302 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Santini, L. et al. One strategy does not fit all: determinants of urban adaptation in mammals. Ecol. Lett. 22, 365–376 (2019).PubMed 
Article 

Google Scholar 
Ostfeld, R. S. et al. Life history and demographic drivers of reservoir competence for three tick-borne zoonotic pathogens. PLoS ONE 9, e107387 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Olival, K. J. et al. Host and viral traits predict zoonotic spillover from mammals. Nature 546, 646–650 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mollentze, N. & Streicker, D. G. Viral zoonotic risk is homogenous among taxonomic orders of mammalian and avian reservoir hosts. Proc. Natl Acad. Sci. USA 117, 9423–9430 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gutiérrez, J. S., Piersma, T. & Thieltges, D. W. Micro- and macroparasite species richness in birds: the role of host life history and ecology. J. Anim. Ecol. 88, 1226–1239 (2019).PubMed 
Article 

Google Scholar 
Teitelbaum, C. S. et al. A comparison of diversity estimators applied to a database of host–parasite associations. Ecography 43, 1316–1328 (2019).Article 

Google Scholar 
Jorge, F. & Poulin, R. Poor geographical match between the distributions of host diversity and parasite discovery effort. Proc. R. Soc. B 285, 20180072 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Allen, T. et al. Global hotspots and correlates of emerging zoonotic diseases. Nat. Commun. 8, 1124 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Gibb, R. et al. Mammal virus diversity estimates are unstable due to accelerating discovery effort. Biol. Lett. https://doi.org/10.1098/rsbl.2021.0427 (2022).Hughes, A. et al. Sampling biases shape our view of the natural world. Ecography 44, 1259–1269 (2021).Article 

Google Scholar 
Estes, L. et al. The spatial and temporal domains of modern ecology. Nat. Ecol. Evol. 2, 819–826 (2018).PubMed 
Article 

Google Scholar 
Titley, M. A., Snaddon, J. L. & Turner, E. C. Scientific research on animal biodiversity is systematically biased towards vertebrates and temperate regions. PLoS ONE 12, e0189577 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Lloyd-Smith, J. O. et al. Should we expect population thresholds for wildlife disease? Trends Ecol. Evol. 20, 511–519 (2005).PubMed 
Article 

Google Scholar 
Cummings, C. R. et al. Foraging in urban environments increases bactericidal capacity in plasma and decreases corticosterone concentrations in white ibises. Front. Ecol. Evol. 8, 575980 (2020).Article 

Google Scholar 
Hwang, J. et al. Anthropogenic food provisioning and immune phenotype: association among supplemental food, body condition, and immunological parameters in urban environments. Ecol. Evol. 8, 3037–3046 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Strandin, T., Babayan, S. A. & Forbes, K. M. Reviewing the effects of food provisioning on wildlife immunity. Phil. Trans. R. Soc. B 373, 20170088 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Downs, C. J., Dochtermann, N. A., Ball, R., Klasing, K. C. & Martin, L. B. The effects of body mass on immune cell concentrations of mammals. Am. Nat. 195, 107–114 (2020).PubMed 
Article 

Google Scholar 
Downs, C. J. et al. Extreme hyperallometry of mammalian antibacterial defenses. Preprint at bioRxiv https://doi.org/10.1101/2020.09.04.242107 (2020).Becker, D. J., Seifert, S. N. & Carlson, C. J. Beyond infection: integrating competence into reservoir host prediction. Trends Ecol. Evol. 35, 1062–1065 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Hanson, D. A., Britten, H. B., Restani, M. & Washburn, L. R. High prevalence of Yersinia pestis in black-tailed prairie dog colonies during an apparent enzootic phase of sylvatic plague. Conserv. Genet. 8, 789–795 (2007).CAS 
Article 

Google Scholar 
Gecchele, L. V., Pedersen, A. B. & Bell, M. Fine-scale variation within urban landscapes affects marking patterns and gastrointestinal parasite diversity in red foxes. Ecol. Evol. 10, 13796–13809 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Albery, G. F., Sweeny, A. R., Becker, D. J. & Bansal, S. Fine-scale spatial patterns of wildlife disease are common and understudied. Funct. Ecol. https://doi.org/10.1111/1365-2435.13942 (2021).Jones, K. E. et al. PanTHERIA: a species-level database of life history, ecology, and geography of extant and recently extinct mammals. Ecology 90, 2648–2648 (2009).Article 

Google Scholar 
Fritz, S. A., Bininda-Emonds, O. R. P. & Purvis, A. Geographical variation in predictors of mammalian extinction risk: big is bad, but only in the tropics. Ecol. Lett. 12, 538–549 (2009).PubMed 
Article 

Google Scholar 
Albery, G. F., Eskew, E. A., Ross, N. & Olival, K. J. Predicting the global mammalian viral sharing network using phylogeography. Nat. Commun. https://doi.org/10.1038/s41467-020-16153-4 (2020).IUCN Red List of Threatened Species Version 2019-2 (IUCN, 2019); https://www.iucnredlist.orgBecker, D. J. et al. Optimising predictive models to prioritise viral discovery in zoonotic reservoirs. Lancet Microbe https://doi.org/10.1016/S2666-5247(21)00245-7 (2022).Mason, P. Parasites of deer in New Zealand. N. Zeal. J. Zool. 21, 39–47 (1994).Article 

Google Scholar 
Wilman, H. et al. EltonTraits 1.0: species-level foraging attributes of the world’s birds and mammals. Ecology 95, 2027 (2014).Article 

Google Scholar 
Plourde, B. T. et al. Are disease reservoirs special? Taxonomic and life history characteristics. PLoS ONE 12, e0180716 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Gibb, R. et al. Data proliferation, reconciliation, and synthesis in viral ecology. Bioscience https://doi.org/10.1101/2021.01.14.426572 (2021).Stephens, P. R. et al. Global mammal parasite database version 2.0. Ecology 98, 1476 (2017).PubMed 
Article 

Google Scholar 
Wardeh, M., Risley, C., Mcintyre, M. K., Setzkorn, C. & Baylis, M. Database of host–pathogen and related species interactions, and their global distribution. Sci. Data 2, 150049 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Shaw, L. P. et al. The phylogenetic range of bacterial and viral pathogens of vertebrates. Mol. Ecol. 29, 3361–3379 (2020).PubMed 
Article 

Google Scholar 
Chamberlain, S. A. & Szöcs, E. taxize: taxonomic search and retrieval in R. F1000Res https://doi.org/10.12688/f1000research.2-191.v2 (2013).Carlson, C. J. et al. The Global Virome in One Network (VIRION): an atlas of vertebrate–virus associations. mBio 13, e0298521 (2022).Article 

Google Scholar 
Lindgren, F. & Rue, H. Bayesian spatial modelling with R-INLA. J. Stat. Softw. 63, 1–25 (2015).Article 

Google Scholar 
Lindgren, F., Rue, H. & Lindstrom, J. An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach. J. R. Stat. Soc. B 73, 423–498 (2011).Article 

Google Scholar 
Hadfield, J. D. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J. Stat. Softw. 33, 1–22 (2010).Article 

Google Scholar 
Winter, D. J. rentrez: an R package for the NCBI eUtils API. R J. 9, 520–526 (2017).Article 

Google Scholar 
Shipley, B. Confirmatory path analysis in a generalized multilevel context. Ecology 90, 363–368 (2009).PubMed 
Article 

Google Scholar 
Carlson, C. J., Dallas, T. A., Alexander, L. W., Phelan, A. L. & Phillips, A. J. What would it take to describe the global diversity of parasites? Proc. R. Soc. B 287, 20201841 (2020).PubMed 
PubMed Central 
Article 

Google Scholar  More