Philippa Kaur

Naef-Daenzer, B. Patch time allocation and patch sampling by foraging great and blue tits. Anim. Behav. 59, 989–999 (2000).CAS 
PubMed 
Article 

Google Scholar 
Kotler, B. P., Brown, J. S. & Bouskila, A. Apprehension and time allocation in gerbils: The effects of predatory risk and energetic state. Ecology 85, 917–922 (2004).Article 

Google Scholar 
Wajnberg, E., Bernhard, P., Hamelin, F. & Boivin, G. Optimal patch time allocation for time-limited foragers. Behav. Ecol. Sociobiol. 60, 1–10 (2006).Article 

Google Scholar 
Embar, K., Kotler, B. P. & Mukherjee, S. Risk management in optimal foragers: The effect of sightlines and predator type on patch use, time allocation, and vigilance in gerbils. Oikos 120, 1657–1666 (2011).Article 

Google Scholar 
Lima, S. L. & Bednekoff, P. A. Temporal variation in danger drives antipredator behavior: The predation risk allocation hypothesis. Am. Nat. 153, 649–659 (1999).PubMed 
Article 

Google Scholar 
Beauchamp, G. & Ruxton, G. D. A reassessment of the predation risk allocation hypothesis: A comment on Lima and Bednekoff. Am. Nat. 177, 143–146 (2011).PubMed 
Article 

Google Scholar 
Ferrari, M. C. O., Sih, A. & Chivers, D. P. The paradox of risk allocation: A review and prospectus. Anim. Behav. 78, 579–585 (2009).Article 

Google Scholar 
Wolf, L. L. & Hainsworth, F. R. Foraging efficiencies and time budgets in nectar-feeding birds. Ecology 56, 117–128 (1975).Article 

Google Scholar 
Litzow, M. A. & Piatt, J. F. Variance in prey abundance influences time budgets of breeding seabirds: Evidence from pigeon guillemots Cepphus columba. J. Avian Biol. 34, 54–64 (2003).Article 

Google Scholar 
Rishworth, G. M., Tremblay, Y. & Green, D. B. Drivers of time-activity budget variability during breeding in a pelagic seabird. PLoS One 9, e116544 (2014).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Stephens, D. W., Brown, J. S. & Ydenberg, R. C. Foraging: Behavior and Ecology. (The University of Chicago Press, 2007).Orians, G. & Pearson, N. On the theory of central place foraging. In Analysis of Ecological Systems (eds. Horn, D., Mitchell, R. & Stairs, G.) 154–177 (The Ohio State University Press, 1979).Chaurand, T. & Weimerskirch, H. The regular alternation of short and long foraging trips in the blue petrel Halobaena caerulea: A previously undescribed strategy of food provisioning in a pelagic seabird. J. Anim. Ecol. 63, 275–282 (1994).Article 

Google Scholar 
Weimerskirch, H. et al. Alternate long and short foraging trips in pelagic seabird parents. Anim. Behav. 47, 472–476 (1994).Article 

Google Scholar 
Welcker, J., Beiersdorf, A., Varpe, Ø. & Steen, H. Mass fluctuations suggest different functions of bimodal foraging trips in a central-place forager. Behav. Ecol. 23, 1372–1378 (2012).Article 

Google Scholar 
Welcker, J. et al. Flexibility in the bimodal foraging strategy of a high Arctic alcid, the little auk Alle alle. J. Avian Biol. 40, 388–399 (2009).Article 

Google Scholar 
Jakubas, D., Wojczulanis-Jakubas, K., Iliszko, L. M. & Kidawa, D. Flexibility of little auks foraging in various oceanographic features in a changing Arctic. Sci. Rep. https://doi.org/10.1038/s41598-020-65210-x (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Shoji, A. et al. Dual foraging and pair coordination during chick provisioning by Manx shearwaters: Empirical evidence supported by a simple model. J. Exp. Biol. 218, 2116–2123 (2015).PubMed 
PubMed Central 

Google Scholar 
Phillips, R. A., Wakefield, E. D., Croxall, J. P., Fukuda, A. & Higuchi, H. Albatross foraging behaviour: No evidence for dual foraging, and limited support for anticipatory regulation of provisioning at South Georgia. Mar. Ecol. Prog. Ser. 391, 279–292 (2009).ADS 
Article 

Google Scholar 
Brown, Z. W., Welcker, J., Harding, A. M. A., Walkusz, W. & Karnovsky, N. J. Divergent diving behavior during short and long trips of a bimodal forager, the little auk Alle alle. J. Avian Biol. 43, 215–226 (2012).Article 

Google Scholar 
Baduini, C. L. & Hyrenbach, K. D. Biogeography of procellariiform foraging strategies: Does ocean productivity influence provisioning?. Mar. Ornithol. 31, 101–112 (2003).
Google Scholar 
Navarro, J. & González-Solís, J. Environmental determinants of foraging strategies in Cory’s shearwaters Calonectris diomedea. Mar. Ecol. Prog. Ser. 378, 259–267 (2009).ADS 
CAS 
Article 

Google Scholar 
Ochi, D., Oka, N. & Watanuki, Y. Foraging trip decisions by the streaked shearwater Calonectris leucomelas depend on both parental and chick state. J. Ethol. 28, 313–321 (2010).Article 

Google Scholar 
Congdon, B. C., Krockenberger, A. K. & Smithers, B. V. Dual-foraging and co-ordinated provisioning in a tropical Procellariiform, the wedge-tailed shearwater. Mar. Ecol. Prog. Ser. 301, 293–301 (2005).ADS 
Article 

Google Scholar 
Peck, D. R. & Congdon, B. C. Colony-specific foraging behaviour and co-ordinated divergence of chick development in the wedge-tailed shearwater Puffinus pacificus. Mar. Ecol. Prog. Ser. 299, 289–296 (2005).ADS 
Article 

Google Scholar 
Weimerskirch, H. How can a pelagic seabird provision its chick when relying on a distant food resource? Cyclic attendance at the colony, foraging decision and body condition in sooty shearwaters. J. Anim. Ecol. 67, 99–109 (1998).Article 

Google Scholar 
Stempniewicz, L. BWP update. Little Auk (Alle alle). J. Birds West. Palearct. 3, 175–201 (2001).
Google Scholar 
Wojczulanis-Jakubas, K. & Jakubas, D. When and why does my mother leave me? The question of brood desertion in the Dovekie (Alle Alle). Auk 129, 632–637 (2012).Article 

Google Scholar 
Harding, A. M. A., Van Pelt, T. I., Lifjeld, J. T. & Mehlum, F. Sex differences in little auk Alle alle parental care: Transition from biparental to paternal-only care. Ibis (Lond. 1859). 146, 642–651 (2004).Article 

Google Scholar 
Wojczulanis-Jakubas, K. et al. Duration of female parental care and their survival in the little auk Alle alle—Are these two traits linked ?. Behav. Ecol. Sociobiol. 74, 1–11 (2020).Article 

Google Scholar 
Wojczulanis, K., Dariusz, J. & Lech, S. The Little Auk Alle alle: An ecological indicator of a changing Arctic and a model organism. Polar Biol. https://doi.org/10.1007/s00300-021-02981-7 (2021).Article 

Google Scholar 
Steen, H., Vogedes, D., Broms, F., Falk-Petersen, S. & Berge, J. Little auks (Alle alle) breeding in a High Arctic fjord system: Bimodal foraging strategies as a response to poor food quality?. Polar Res. 26, 118–125 (2007).Article 

Google Scholar 
Wojczulanis-Jakubas, K., Jakubas, D., Karnovsky, N. J. & Walkusz, W. Foraging strategy of little auks under divergent conditions on feeding grounds. Polar Res. 29, 22–29 (2010).Article 

Google Scholar 
Jakubas, D., Wojczulanis-Jakubas, K., Iliszko, L., Darecki, M. & Stempniewicz, L. Foraging strategy of the little auk Alle alle throughout breeding season—switch from unimodal to bimodal pattern. J. Avian Biol. 45, 551–560 (2014).Article 

Google Scholar 
Jakubas, D., Iliszko, L., Wojczulanis-Jakubas, K. & Stempniewicz, L. Foraging by little auks in the distant marginal sea ice zone during the chick-rearing period. Polar Biol. 35, 73–81 (2012).Article 

Google Scholar 
Jakubas, D. et al. Intra-seasonal variation in zooplankton availability, chick diet and breeding performance of a high Arctic planktivorous seabird. Polar Biol. 391, 1547–1561 (2016).Article 

Google Scholar 
Jakubas, D. et al. Foraging closer to the colony leads to faster growth in little auks. Mar. Ecol. Prog. Ser. 489, 263–278 (2013).ADS 
Article 

Google Scholar 
Stempniewicz, L. Predator-prey interactions between Glaucous Gull Larus hyperboreus and Little Auk Alle alle in Spitsbergen. Acta Ornithol. 29, 155–170 (1995).
Google Scholar 
Wojczulanis-Jakubas, K., Jakubas, D. & Stempniewicz, L. Changes in the glaucous gull predatory pressure on little auks in Southwest Spitsbergen. Waterbirds 28, 430–435 (2005).Article 

Google Scholar 
Kharitonov, S. Methods and Theoretical Aspects of Seabird Studies. (Proc 5 All-Russian Mar Biol School, Marine Biological Institute, 2007).Wojczulanis-Jakubas, K., Jakubas, D. & Stempniewicz, L. Avifauna of Hornsund area, SW Spitsbergen: Present state and recent changes. Polish Polar Res. 29, 187–197 (2008).
Google Scholar 
Keslinka, K. L., Wojczulanis-Jakubas, K., Jakubas, D. & Neubauer, G. Determinants of the little auk (Alle alle) breeding colony location and size in W and NW coast of Spitsbergen. PLoS One 14, 1–20 (2019).
Google Scholar 
Kidawa, D., Barcikowski, M. & Palme, R. Parent-offspring interactions in a long-lived seabird, the Little Auk (Alle alle): Begging and provisioning under simulated stress. J. Ornithol. 158, 145–157 (2017).Article 

Google Scholar 
Welcker, J., Beiersdorf, A., Varpe, Ø. & Steen, H. Mass fluctuations suggest different functions of bimodal foraging trips in a central-place forager. Behav. Ecol. https://doi.org/10.1093/beheco/ars131 (2012).Article 

Google Scholar 
Jakubas, D. & Wojczulanis, K. Predicting the sex of Dovekies by discriminant analysis. Waterbirds 30, 92–96 (2007).Article 

Google Scholar 
Grissot, A. et al. Parental coordination of chick provisioning in a planktivorous arctic seabird under divergent conditions on foraging grounds. Front. Ecol. Evol. 7, 349 (2019).Article 

Google Scholar 
Stoffel, M. A., Nakagawa, S. & Schielzeth, H. rptR: Repeatability estimation and variance decomposition by generalized linear mixed-effects models. Methods Ecol. Evol. 8, 1639–1644 (2017).Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing. R. (2019).Wojczulanis-Jakubas, K., Jakubas, D. & Stempniewicz, L. Sex-specific parental care by incubating Little Auks (Alle alle). Ornis Fenn. 86, 140–148 (2009).
Google Scholar 
Welcker, J., Steen, H., Harding, A. M. A. & Gabrielsen, G. W. Sex-specific provisioning behaviour in a monomorphic seabird with a bimodal foraging strategy. Ibis (Lond. 1859). 151, 502–513 (2009).Article 

Google Scholar 
Kidawa, D. et al. Parental efforts of an Arctic seabird, the little auk Alle alle under variable foraging conditions. Mar. Biol. Res. 11, 349–360 (2015).Article 

Google Scholar 
Wickham, H. Hadley Wickham. Media 35, 211 (2009).
Google Scholar 
Karnovsky, N. J. et al. Inter-colony comparison of diving behavior of an Arctic top predator: Implications for warming in the Greenland Sea. Mar. Ecol. Prog. Ser. 440, 229–240 (2011).ADS 
Article 

Google Scholar 
Karnovsky, N. et al. Foraging distributions of little auks Alle alle across the Greenland Sea: Implications of present and future Arctic climate change. Mar. Ecol. Prog. Ser. 415, 283–293 (2010).ADS 
Article 

Google Scholar 
Gremillet, D. et al. Little auks buffer the impact of current Arctic climate change. Mar. Ecol. Prog. Ser. 454, 197–206 (2012).ADS 
Article 

Google Scholar 
Harding, A. M. A. et al. Flexibility in the parental effort of an Arctic-breeding seabird. Funct. Ecol. 23, 348–358 (2009).Article 

Google Scholar 
Jakubas, D. et al. Foraging effort does not influence body condition and stress level in little auks. Mar. Ecol. Prog. Ser. 432, 277–290 (2011).ADS 
Article 

Google Scholar 
Jakubas, D., Wojczulanis-Jakubas, K., Iliszko, L. M., Strøm, H. & Stempniewicz, L. Habitat foraging niche of a High Arctic zooplanktivorous seabird in a changing environment. Sci. Rep. 7, 1–14 (2017).CAS 
Article 

Google Scholar  More

MacLachlan, N. J. & Guthrie, A. J. Re-emergence of bluetongue, African horse sickness, and other Orbivirus diseases. Vet. Res. https://doi.org/10.1051/vetres/2010007 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Koenraadt, C. J. M. et al. Bluetongue, Schmallenberg—What is next? Culicoides-borne viral diseases in the 21st Century. BMC Res. Notes 10, 77 (2014).
Google Scholar 
Dennis, S. J., Meyers, A. E., Hitzeroth, I. I. & Rybicki, E. P. African horse sickness: A review of current understanding and vaccine development in the. Viruses 11, 844 (2019).CAS 
PubMed Central 
Article 

Google Scholar 
Collins, Á. B., Doherty, M. L., Barrett, D. J. & Mee, J. F. Schmallenberg virus: A systematic international literature review (2011–2019) from an Irish perspective. Ir. Vet. J. 72, 1–22 (2019).Article 

Google Scholar 
Tkuwet, G. & Firesbhat, A. A review on African horse sickness. Eur. J. Appl. Sci. 7, 213–219 (2015).CAS 

Google Scholar 
Mellor, P. S. & Hamblin, C. African horse sickness. Vet. Res. 35, 445–466 (2004).PubMed 
Article 

Google Scholar 
Coetzee, P., Stokstad, M., Venter, E. H., Myrmel, M. & Van Vuuren, M. Bluetongue: A historical and epidemiological perspective with the emphasis on South Africa. Virol. J. 9, 198 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Cagienard, A., Griot, C., Mellor, P. S., Denison, E. & Stärk, K. D. Bluetongue vector species of Culicoides in Switzerland. Med. Vet. Entomol. 20, 239–247 (2006).CAS 
PubMed 
Article 

Google Scholar 
Oluwayelu, D., Adebiyi, A. & Tomori, O. Endemic and emerging arboviral diseases of livestock in Nigeria: A review. Parasit. Vectors 11, 1–12 (2018).Article 

Google Scholar 
Sibhat, B., Ayelet, G., Gebremedhin, E. Z., Skjerve, E. & Asmare, K. Seroprevalence of Schmallenberg virus in dairy cattle in Ethiopia. Acta Trop. 178, 61–67 (2018).PubMed 
Article 

Google Scholar 
Aklilu, N. et al. African horse sickness outbreaks caused by multiple virus types in Ethiopia. Transbound. Emerg. Dis. 61, 185–192 (2014).CAS 
PubMed 
Article 

Google Scholar 
Rojas, J. M., Rodríguez-Martín, D., Martín, V. & Sevilla, N. Diagnosing bluetongue virus in domestic ruminants: Current perspectives. Vet. Med. Res. Rep. 10, 17 (2019).
Google Scholar 
Gizaw, D., Sibhat, D., Ayalew, B. & Sehal, M. Sero-prevalence study of bluetongue infection in sheep and goats in selected areas of Ethiopia. Ethiop. Vet. J. 20, 105 (2016).Article 

Google Scholar 
Abera, T. et al. Bluetongue disease in small ruminants in south western Ethiopia: Cross-sectional sero-epidemiological study. BMC Res. Notes 11, 112 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Mellor, P. S., Boorman, J. & Baylis, M. Culicoides biting midges: Their role as arbovirus vectors. Annu. Rev. Entomol. 45, 307–340 (2000).CAS 
PubMed 
Article 

Google Scholar 
Carpenter, S., Groschup, M. H., Garros, C., Felippe-Bauer, M. L. & Purse, B. V. Culicoides biting midges, arboviruses and public health in Europe. Antivir. Res. 100, 102–113 (2013).CAS 
PubMed 
Article 

Google Scholar 
Sick, F., Beer, M., Kampen, H. & Wernike, K. Culicoides biting midges—Underestimated vectors for arboviruses of public health and veterinary importance. Viruses 11, 376 (2019).CAS 
PubMed Central 
Article 

Google Scholar 
Blanda, V. et al. Geo-statistical analysis of Culicoides spp. distribution and abundance in Sicily, Italy. Parasit. Vectors 11, 78 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Vasić, A. et al. Species diversity, host preference and arbovirus detection of Culicoides (Diptera: Ceratopogonidae) in south-eastern Serbia. Parasit. Vectors 12, 1–9 (2019).Article 

Google Scholar 
Martin, E. et al. Culicoides species community composition and infection status with parasites in an urban environment of east central Texas, USA. Parasit. Vectors 12, 1–10 (2019).Article 

Google Scholar 
Gusmão, G. M. C., Brito, G. A., Moraes, L. S., Bandeira, M. D. C. A. & Rebêlo, J. M. M. Temporal variation in species abundance and richness of Culicoides (Diptera: Ceratopogonidae) in a tropical equatorial area. J. Med. Entomol. https://doi.org/10.1093/jme/tjz015 (2019).Article 
PubMed 

Google Scholar 
Sghaier, S. et al. New species of the genus Culicoides (Diptera Ceratopogonidae) for Tunisia, with detection of Bluetongue viruses in vectors. Vet. Ital. 53, 357–366 (2017).PubMed 

Google Scholar 
Gordon, S. J. G. et al. The occurrence of Culicoides species, the vectors of arboviruses, at selected trap sites in Zimbabwe. Onderstepoort J. Vet. Res. 82, e1–e8 (2015).PubMed 
Article 
CAS 

Google Scholar 
Villard, P. et al. Modeling Culicoides abundance in mainland France: Implications for surveillance. Parasit. Vectors 12, 1–10 (2019).Article 

Google Scholar 
Diarra, M. et al. Spatial distribution modelling of Culicoides (Diptera: Ceratopogonidae) biting midges, potential vectors of African horse sickness and bluetongue viruses in Senegal. Parasit. Vectors 11, 341 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Calvete, C. et al. Spatial distribution of Culicoides imicola, the main vector of bluetongue virus, Spain. Vet. Rec. 158, 130–131 (2006).CAS 
PubMed 
Article 

Google Scholar 
Purse, B. V. et al. Modelling the distributions of Culicoides bluetongue virus vectors in Sicily in relation to satellite-derived climate variables. Med. Vet. Entomol. 18, 90–101 (2004).CAS 
PubMed 
Article 

Google Scholar 
Purse, B. V. et al. Spatial and temporal distribution of bluetongue and its Culicoides vectors in Bulgaria. Med. Vet. Entomol. 20, 335–344 (2006).CAS 
PubMed 
Article 

Google Scholar 
Leta, S. et al. Modeling the global distribution of Culicoides imicola: An ensemble approach. Sci. Rep. 9, 1–9 (2019).ADS 
CAS 
Article 

Google Scholar 
Mulatu, T. & Hailu, A. The occurrence and identification of Culicoides species in the Western Ethiopia. Acad. J. Entomol. 12, 40–43 (2019).
Google Scholar 
Khamala, C. P. M. & Kettle, D. S. The Culicoides Latreille (Diptera: Ceratopogonidae) of East Africa. Trans. R. Entomol. Soc. Lond. 123, 1–95 (1971).Article 

Google Scholar 
Venter, G. J. Specie di Culicoides (Diptera: Ceratopogonidae) vettori del virus della Bluetongue in Sud Africa. Vet. Ital. 51, 325–333 (2015).PubMed 

Google Scholar 
Mathieu, B. et al. Development and validation of IIKC: An interactive identification key for Culicoides (Diptera: Ceratopogonidae) females from the Western Palaearctic region. Parasit. Vectors 5, 137 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Thuiller, W., Lafourcade, B., Engler, R. & Araújo, M. B. BIOMOD—A platform for ensemble forecasting of species distributions. Ecography (Cop.) 32, 369–373 (2009).Article 

Google Scholar 
Baylis, M., Bouayoune, H., Touti, J. & El Hasnaoui, H. Use of climatic data and satellite imagery to model the abundance of Culicoides imicola, the vector of African horse sickness virus, in Morocco. Med. Vet. Entomol. 12, 255–266 (1998).CAS 
PubMed 
Article 

Google Scholar 
Diarra, M. et al. Modelling the abundances of two major culicoides (Diptera: Ceratopogonidae) species in the Niayes area of Senegal. PLoS One 10, e0131021 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ramilo, D. W., Nunes, T., Madeira, S., Boinas, F. & da Fonseca, I. P. Geographical distribution of Culicoides (DIPTERA: CERATOPOGONIDAE) in mainland Portugal: Presence/absence modelling of vector and potential vector species. PLoS One 12, e0180606 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ben Rais Lasram, F. et al. The Mediterranean Sea as a ‘cul-de-sac’ for endemic fishes facing climate change. Glob. Chang. Biol. 16, 3233–3245 (2010).ADS 
Article 

Google Scholar 
Tiffin, P. & Ross-Ibarra, J. Goal-oriented evaluation of species distribution models accuracy and precision: True Skill Statistic profile and uncertainty maps. PeerJ PrePints https://doi.org/10.7287/peerj.preprints.488v1 (2014).Article 

Google Scholar 
Graham, M. H. Confronting multicollinearity in ecological multiple regression. Ecology 84, 2809–2815 (2003).Article 

Google Scholar 
Demissie, G. H. Seroepidemiological study of African horse sickness in southern Ethiopia. Open Sci. Repos. Vet. Med. 10, e70081919 (2013).
Google Scholar 
Zeleke, A., Sori, T., Powel, K., Gebre-Ab, F. & Endebu, B. Isolation and identification of circulating serotypes of African horse sickness virus in Ethiopia. J. Appl. Res. Vet. Med. 3, 40–43 (2005).
Google Scholar 
Ayelet, G. et al. Outbreak investigation and molecular characterization of African horse sickness virus circulating in selected areas of Ethiopia. Acta Trop. 127, 91–96 (2013).PubMed 
Article 

Google Scholar 
Gulima, D. Seroepidemiological study of bluetongue in indigenous sheep in selected districts of Amhara National Regional State, north western Ethiopia. Ethiop. Vet. J. 13, 1–15 (2009).
Google Scholar 
Borkent, A. & Dominiak, P. Catalog of the biting midges of the world (Diptera: Ceratopogonidae). Zootaxa 4787, 1–377 (2020).Article 

Google Scholar 
Borkent, A. & Wirth, W. W. World species of biting midges (Diptera: Ceratopogonidae). Bull. Am. Museum Nat. Hist. 233, 5–195 (1997).
Google Scholar 
Guichard, S. et al. Worldwide niche and future potential distribution of Culicoides imicola, a major vector of bluetongue and African horse sickness viruses. PLoS One 9, e112491 (2014).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Becker, E. E. E., Venter, G. J., Labuschagne, K., Greyling, T. & van Hamburg, H. Occurrence of Culicoides species Diptera: Ceratopogonidae) in the Khomas region of Namibia during the winter months. Vet. Ital. 48, 45–54 (2012).PubMed 

Google Scholar 
Capela, R. et al. Spatial distribution of Culicoides species in Portugal in relation to the transmission of African horse sickness and bluetongue viruses. Med. Vet. Entomol. 17, 165–177 (2003).CAS 
PubMed 
Article 

Google Scholar 
Calvete, C. et al. Modelling the distributions and spatial coincidence of bluetongue vectors Culicoides imicola and the Culicoides obsoletus group throughout the Iberian peninsula. Med. Vet. Entomol. 22, 124–134 (2008).CAS 
PubMed 
Article 

Google Scholar 
Riddin, M. A., Venter, G. J., Labuschagne, K. & Villet, M. H. Culicoides species as potential vectors of African horse sickness virus in the southern regions of South Africa. Med. Vet. Entomol. 33, 498–511 (2019).CAS 
PubMed 
Article 

Google Scholar 
Foxi, C. et al. Role of different Culicoides vectors (Diptera: Ceratopogonidae) in bluetongue virus transmission and overwintering in Sardinia (Italy). Parasit. Vectors 9, 440 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Musuka, G. N., Mellor, P. S., Meiswinkel, R., Baylis, M. & Kelly, P. J. Prevalence of Culicoides imicola and other species (Diptera: Ceratopogonidae) ateight sites in Zimbabwe: To the editor. J. S. Afr. Vet. Assoc. 72, 62–63 (2001).CAS 
PubMed 
Article 

Google Scholar 
Meiswinkel, R. The 1996 outbreak of African horse sickness in South Africa—the entomological perspective. Arch. Virol. Suppl. 14, 69–83 (1998).CAS 
PubMed 

Google Scholar 
Jean Pierre, T. et al. Characteristics, classification and genesis of vertisols under seasonally contrasted climate in the Lake Chad Basin, Central Africa. J. Afr. Earth Sci. 150, 176–193 (2019).ADS 
CAS 
Article 

Google Scholar 
Elias, E. Characteristics of Nitisol profiles as affected by land use type and slope class in some Ethiopian highlands. Environ. Syst. Res. 6, 1–15 (2017).Article 

Google Scholar 
Nachtergaele, F. The classification of leptosols in the world reference base for soil resources.Veronesi, E., Venter, G. J., Labuschagne, K., Mellor, P. S. & Carpenter, S. Life-history parameters of Culicoides (Avaritia) imicola Kieffer in the laboratory at different rearing temperatures. Vet. Parasitol. 163, 370–373 (2009).CAS 
PubMed 
Article 

Google Scholar 
Verhoef, F. A. A., Venter, G. J. & Weldon, C. W. Thermal limits of two biting midges, Culicoides imicola Kieffer and C. bolitinos Meiswinkel (Diptera: Ceratopogonidae). Parasites Vectors 7, 384 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Conte, A., Goffredo, M., Ippoliti, C. & Meiswinkel, R. Influence of biotic and abiotic factors on the distribution and abundance of Culicoides imicola and the Obsoletus Complex in Italy. Vet. Parasitol. 150, 333–344 (2007).CAS 
PubMed 
Article 

Google Scholar 
Martinez-de la Puente, J., Navarro, J., Ferraguti, M., Soriguer, R. & Figuerola, J. First molecular identification of the vertebrate hosts of Culicoides imicola in Europe and a review of its blood-feeding patterns worldwide: Implications for the transmission of bluetongue disease and African horse sickness. Med. Vet. Entomol. 31, 333–339 (2017).CAS 
PubMed 
Article 

Google Scholar 
Purse, B. V. et al. Impacts of climate, host and landscape factors on Culicoides species in Scotland. Med. Vet. Entomol. 26, 168–177 (2012).CAS 
PubMed 
Article 

Google Scholar 
Leta, S. et al. Updating the global occurrence of Culicoides imicola, a vector for emerging viral diseases. Sci. Data 6, 1–8 (2019).CAS 
Article 

Google Scholar  More

Vázquez, D. P., Gianoli, E., Morris, W. F. & Bozinovic, F. Ecological and evolutionary impacts of changing climatic variability. Biol. Rev. 92, 22–42. https://doi.org/10.1111/brv.12216 (2017).Article 
PubMed 

Google Scholar 
Marshall, K. E. & Sinclair, B. J. The impacts of repeated cold exposure on insects. J. Exp. Biol. 215, 1607–1613. https://doi.org/10.1242/jeb.059956 (2012).Article 
PubMed 

Google Scholar 
Bale, J. & Hayward, S. Insect overwintering in a changing climate. J. Exp. Biol. 213, 980–994 (2010).CAS 
PubMed 
Article 

Google Scholar 
Kingsolver, J. G. Feeding, growth, and the thermal environment of cabbage white caterpillars, Pieris rapae L. Physiol. Biochem. Zool. 73, 621–628 (2000).CAS 
PubMed 
Article 

Google Scholar 
Stange, E. E. & Ayres, M. P. Climate change impacts: Insects (JohnWiley & Sons, 2010).
Google Scholar 
Chapman, A. D. Numbers of Living Species in Australia and the World: Report for the Department of the Environment and Heritage Canberra, Australia (Department of the Environment and Heritage, 2006).
Google Scholar 
Colinet, H., Sinclair, B. J., Vernon, P. & Renault, D. Insects in fluctuating thermal environments. Annu. Rev. Entomol. 60, 123–140. https://doi.org/10.1146/annurev-ento-010814-021017 (2015).CAS 
Article 
PubMed 

Google Scholar 
Hickling, R., Roy, D. B., Hill, J. K. & Thomas, C. D. A northward shift of range margins in British Odonata. Glob. Change Biol. 11, 502–506. https://doi.org/10.1111/j.1365-2486.2005.00904.x (2005).ADS 
Article 

Google Scholar 
Rumpf, S. B., Hülber, K., Zimmermann, N. E. & Dullinger, S. Elevational rear edges shifted at least as much as leading edges over the last century. Glob. Ecol. Biogeogr. 28, 533–543. https://doi.org/10.1111/geb.12865 (2019).Article 

Google Scholar 
Halsch, C. A. et al. Insects and recent climate change. Proc. Natl. Acad. Sci. 118, e2002543117. https://doi.org/10.1073/pnas.2002543117 (2021).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
McCain, C. M. & Garfinkel, C. F. Climate change and elevational range shifts in insects. Curr. Opin. Insect Sci. 47, 111–118. https://doi.org/10.1016/j.cois.2021.06.003 (2021).Article 
PubMed 

Google Scholar 
Angilletta, M. J. Jr. & Angilletta, M. J. Thermal Adaptation: A Theoretical and Empirical Synthesis (Oxford University Press, 2009).Book 

Google Scholar 
Angilletta, M. J. & Dunham, A. E. The temperature-size rule in ectotherms: Simple evolutionary explanations may not be general. Am. Nat. 162, 332–342. https://doi.org/10.1086/377187 (2003).Article 
PubMed 

Google Scholar 
Jensen, J. L. W. V. Sur les fonctions convexes et les inégalités entre les valeurs moyennes. Acta Math. 30, 175–193 (1906).MathSciNet 
MATH 
Article 

Google Scholar 
Ruel, J. J. & Ayres, M. P. Jensen’s inequality predicts effects of environmental variation. Trends Ecol. Evol. 14, 361–366 (1999).CAS 
PubMed 
Article 

Google Scholar 
Kingsolver, J. G. & Woods, H. A. Thermal sensitivity of growth and feeding in Manduca sexta Caterpillars. Physiol. Zool. 70, 631–638. https://doi.org/10.1086/515872 (1997).CAS 
Article 
PubMed 

Google Scholar 
Gillooly, J. F., Brown, J. H., West, G. B., Savage, V. M. & Charnov, E. L. Effects of size and temperature on metabolic rate. Science 293, 2248–2251 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Bale, J. S. et al. Herbivory in global climate change research: Direct effects of rising temperature on insect herbivores. Glob. Change Biol. 8, 1–16 (2002).ADS 
Article 

Google Scholar 
Robinet, C. & Roques, A. Direct impacts of recent climate warming on insect populations. Integr. Zool. 5, 132–142 (2010).PubMed 
Article 

Google Scholar 
García-Robledo, C., Kuprewicz, E. K., Staines, C. L., Erwin, T. L. & Kress, W. J. Limited tolerance by insects to high temperatures across tropical elevational gradients and the implications of global warming for extinction. Proc. Natl. Acad. Sci. 113, 680–685. https://doi.org/10.1073/pnas.1507681113 (2016).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Vasseur, D. A. et al. Increased temperature variation poses a greater risk to species than climate warming. Proc. R. Soc. B Biol. Sci. 281, 20132612 (2014).Article 

Google Scholar 
Sandehson, D. E. The relation of temperature to the growth of insects. J. Econ. Entomol. 3, 113–140 (1910).Article 

Google Scholar 
Cook, W. C. Some Effects of Alternating Temperatures on the Growth and Metabolism of Cutworm Larvae (Oxford University Press, 1927).
Google Scholar 
Kingsolver, J. G., Ragland, G. J. & Diamond, S. E. Evolution in a constant environment: Thermal fluctuations and thermal sensitivity of laboratory and field populations of Manduca sexta. Evolution 63, 537–541. https://doi.org/10.1111/j.1558-5646.2008.00568.x (2009).Article 
PubMed 

Google Scholar 
Eldridge, W. H., Sweeney, B. W. & Law, J. M. Fish growth, physiological stress, and tissue condition in response to rate of temperature change during cool or warm diel thermal cycles. Can. J. Fish. Aquat. Sci. 72, 1527–1537 (2015).CAS 
Article 

Google Scholar 
Bernhardt, J. R., Sunday, J. M., Thompson, P. L. & O’Connor, M. I. Nonlinear averaging of thermal experience predicts population growth rates in a thermally variable environment. Proc. R. Soc. B Biol. Sci. 285, 20181076. https://doi.org/10.1098/rspb.2018.1076 (2018).Article 

Google Scholar 
Morissette, J., Swart, S., Maccormack, T. J., Currie, S. & Morash, A. J. Thermal variation near the thermal optimum does not affect the growth, metabolism or swimming performance in wild Atlantic salmon Salmo salar. J. Fish Biol. 98, 1585–1589. https://doi.org/10.1111/jfb.14348 (2021).CAS 
Article 
PubMed 

Google Scholar 
Bozinovic, F. et al. The mean and variance of environmental temperature interact to determine physiological tolerance and fitness. Physiol. Biochem. Zool. 84, 543–552 (2011).PubMed 
Article 

Google Scholar 
Boggs, C. L. The fingerprints of global climate change on insect populations. Curr. Opin. Insect Sci. 17, 69–73. https://doi.org/10.1016/j.cois.2016.07.004 (2016).Article 
PubMed 

Google Scholar 
Lemoine, N. P., Drews, W. A., Burkepile, D. E. & Parker, J. D. Increased temperature alters feeding behavior of a generalist herbivore. Oikos 122, 1669–1678. https://doi.org/10.1111/j.1600-0706.2013.00457.x (2013).Article 

Google Scholar 
Vangansbeke, D. et al. Prey consumption by phytoseiid spider mite predators as affected by diurnal temperature variations. Biocontrol 60, 595–603 (2015).Article 

Google Scholar 
Davies, C., Coetzee, M. & Lyons, C. L. Effect of stable and fluctuating temperatures on the life history traits of Anopheles arabiensis and An. quadriannulatus under conditions of inter-and intra-specific competition. Parasit. Vectors 9, 1–9 (2016).Article 

Google Scholar 
Delava, E., Fleury, F. & Gibert, P. Effects of daily fluctuating temperatures on the Drosophila-Leptopilina boulardi parasitoid association. J. Therm. Biol. 60, 95–102 (2016).PubMed 
Article 

Google Scholar 
Amarasekare, P. & Coutinho, R. M. Effects of temperature on intraspecific competition in ectotherms. Am. Nat. 184, E50–E65. https://doi.org/10.1086/677386 (2014).Article 
PubMed 

Google Scholar 
Jiang, L. & Morin, P. J. Temperature-dependent interactions explain unexpected responses to environmental warming in communities of competitors. J. Anim. Ecol. 73, 569–576 (2004).Article 

Google Scholar 
Novich, R. A., Erickson, E. K., Kalinoski, R. M. & DeLong, J. P. The temperature independence of interaction strength in a sit-and-wait predator. Ecosphere 5, 1–9 (2014).Article 

Google Scholar 
Fox, L. R. Cannibalism in natural populations. Annu. Rev. Ecol. Syst. 6, 87–106 (1975).Article 

Google Scholar 
Polis, G. A. The evolution and dynamics of intraspecific predation. Annu. Rev. Ecol. Syst. 12, 225–251 (1981).Article 

Google Scholar 
Nishimura, K. & Isoda, Y. Evolution of cannibalism: Referring to costs of cannibalism. J. Theor. Biol. 226, 293–302. https://doi.org/10.1016/j.jtbi.2003.09.007 (2004).ADS 
MathSciNet 
Article 
PubMed 
MATH 

Google Scholar 
Crumrine, P. W. Body size, temperature, and seasonal differences in size structure influence the occurrence of cannibalism in larvae of the migratory dragonfly, Anax junius. Aquat. Ecol. 44, 761–770 (2010).Article 

Google Scholar 
Reglero, P., Urtizberea, A., Torres, A. P., Alemany, F. & Fiksen, Ø. Cannibalism among size classes of larvae may be a substantial mortality component in tuna. Mar. Ecol. Prog. Ser. 433, 205–219 (2011).ADS 
Article 

Google Scholar 
Nilsson-Örtman, V., Stoks, R. & Johansson, F. Competitive interactions modify the temperature dependence of damselfly growth rates. Ecology 95, 1394–1406. https://doi.org/10.1890/13-0875.1 (2014).Article 
PubMed 

Google Scholar 
Pritchard, G. & Leggott, M. Temperature, incubation rates and the origins of dragonflies. Adv. Odonatol. 3, 121–126 (1987).
Google Scholar 
Hassall, C. & Thompson, D. J. The effects of environmental warming on Odonata: A review. Int. J. Odonatol. 11, 131–153 (2008).Article 

Google Scholar 
Johansson, F. & Crowley, P. H. Larval cannibalism and population dynamics of dragonflies. In Aquatic Insects: Challenges to Populations 36–54 (CABI, 2008).Rudolf, V. H. W. & Rasmussen, N. L. Ontogenetic functional diversity: Size structure of a keystone predator drives functioning of a complex ecosystem. Ecology 94, 1046–1056 (2013).PubMed 
Article 

Google Scholar 
Hyeun-Ji, L. & Johansson, F. Compensating for a bad start: Compensatory growth across life stages in an organism with a complex life cycle. Can. J. Zool. 94, 41–47 (2016).Article 

Google Scholar 
Sokolovska, N., Rowe, L. & Johansson, F. Fitness and body size in mature odonates. Ecol. Entomol. 25, 239–248. https://doi.org/10.1046/j.1365-2311.2000.00251.x (2000).Article 

Google Scholar 
Karl, T. R. Modern global climate change. Science 302, 1719–1723. https://doi.org/10.1126/science.1090228 (2003).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Meehl, G. A. et al. Climate Change 2007: the physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge University Press, 2007).
Google Scholar 
Meehl, G. A. More intense, more frequent, and longer lasting heat waves in the 21st century. Science 305, 994–997. https://doi.org/10.1126/science.1098704 (2004).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Khelifa, R. Spatiotemporal Pattern of Phenology Across Geographic Gradients in Insects, in Chapter 1 (Geographic Gradients in Climate Change Response Explained by Non-linear Thermal-Performance Curves) (University of Zurich, 2017).
Google Scholar 
Boudot, J. P. & Kalkman, V. Atlas of the European Dragonflies and Damselflies (KNNV Publishing, 2015).
Google Scholar 
Norling, U. Growth, winter preparations and timing of emergence in temperate zone Odonata: Control by a succession of larval response patterns. Int. J. Odonatol. 24, 1–36 (2021).Article 

Google Scholar 
Sniegula, S. & Johansson, F. Photoperiod affects compensating developmental rate across latitudes in the damselfly Lestes sponsa. Ecol. Entomol. 35, 149–157. https://doi.org/10.1111/j.1365-2311.2009.01164.x (2010).Article 

Google Scholar 
Sniegula, S., Golab, M. J. & Johansson, F. Size-mediated priority and temperature effects on intra-cohort competition and cannibalism in a damselfly. J. Anim. Ecol. 88, 637–648. https://doi.org/10.1111/1365-2656.12947 (2019).Article 
PubMed 

Google Scholar 
Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675. https://doi.org/10.1038/nmeth.2089 (2012).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Benke, A. C. A method for comparing individual growth rates of aquatic insects with special reference to the Odonata. Ecology 51, 328–331 (1970).Article 

Google Scholar 
Nilsson-Örtman, V., Stoks, R., De Block, M. & Johansson, F. Generalists and specialists along a latitudinal transect: Patterns of thermal adaptation in six species of damselflies. Ecology 93, 1340–1352. https://doi.org/10.1890/11-1910.1 (2012).Article 
PubMed 

Google Scholar 
Eklund, A. et al. Sveriges Framtida Klimat: Underlag Till Dricksvattenutredningen (SMHI, 2015).
Google Scholar 
McPeek, M. A. Determination of species composition in the Enallagma damselfly assemblages of permanent lakes. Ecology 71, 83–98. https://doi.org/10.2307/1940249 (1990).Article 

Google Scholar 
Kirillin, G. et al. FLake-global: Online lake model with worldwide coverage. Environ. Model. Softw. 26, 683–684. https://doi.org/10.1016/j.envsoft.2010.12.004 (2011).Article 

Google Scholar 
SMHI. Advanced Climate Change Scenario Service. https://www.smhi.se.R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019).Suhling, F., Suhling, I. & Richter, O. Temperature response of growth of larval dragonflies–An overview. Int. J. Odonatol. 18, 15–30 (2015).Article 

Google Scholar 
Padfield, D., O’Sullivan, H. & Pawar, S. rTPC and nls.multstart: A new pipeline to fit thermal performance curves in R. Methods Ecol. Evol. 1, 1. https://doi.org/10.1111/2041-210X.13585 (2021).Article 

Google Scholar 
Kleiber, M. Body size and metabolism. Hilgardia 6, 315–353 (1932).CAS 
Article 

Google Scholar 
Hemmingsen, A. Reports of the Steno Memorial Hospital and Nordisk Insulin Laboratorium. Energy Metab. Relat. Body Size Respir. Surf. Evol. 9, 6–110 (1960).
Google Scholar 
Dell, A. I., Pawar, S. & Savage, V. M. Systematic variation in the temperature dependence of physiological and ecological traits. Proc. Natl. Acad. Sci. 108, 10591–10596 (2011).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Logan, J. D., Wolesensky, W. & Joern, A. Temperature-dependent phenology and predation in arthropod systems. Ecol. Model. 196, 471–482. https://doi.org/10.1016/j.ecolmodel.2006.02.034 (2006).Article 

Google Scholar 
Pink, M. & Abrahams, M. V. Temperature and its impact on predation risk within aquatic ecosystems. Can. J. Fish. Aquat. Sci. 73, 869–876. https://doi.org/10.1139/cjfas-2015-0302 (2016).Article 

Google Scholar 
DeAngelis, D., Cox, D. & Coutant, C. Cannibalism and size dispersal in young-of-the-year largemouth bass: Experiment and model. Ecol. Model. 8, 133–148 (1980).Article 

Google Scholar 
Fagan, W. F. & Odell, G. M. Size-dependent cannibalism in praying mantids: Using biomass flux to model size-structured populations. Am. Nat. 147, 230–268 (1996).Article 

Google Scholar 
Dong, Q. & Deangelis, D. L. Consequences of cannibalism and competition for food in a smallmouth bass population: An individual-based modeling study. Trans. Am. Fish. Soc. 127, 174–191 (1998).Article 

Google Scholar 
Verheyen, J. & Stoks, R. Temperature variation makes an ectotherm more sensitive to global warming unless thermal evolution occurs. J. Anim. Ecol. 88, 624–636. https://doi.org/10.1111/1365-2656.12946 (2019).Article 
PubMed 

Google Scholar 
Starr, S. M. & McIntyre, N. E. Effects of water temperature under projected climate change on the development and survival of Enallagma civile (Odonata: Coenagrionidae). Environ. Entomol. 49, 230–237. https://doi.org/10.1093/ee/nvz138 (2020).CAS 
Article 
PubMed 

Google Scholar 
Culler, L. E., McPeek, M. A. & Ayres, M. P. Predation risk shapes thermal physiology of a predaceous damselfly. Oecologia 176, 653–660. https://doi.org/10.1007/s00442-014-3058-8 (2014).ADS 
Article 
PubMed 

Google Scholar 
Dokulil, M. T. et al. Increasing maximum lake surface temperature under climate change. Clim. Change https://doi.org/10.1007/s10584-021-03085-1 (2021).Article 

Google Scholar 
Merritt, R. W. & Cummins, K. W. An Introduction to the Aquatic Insects of North America 2nd edn. (Kendall/Hunt Publishing Company, 1984).
Google Scholar 
Verdonschot, R. & Peeters, E. T. Preference of larvae of Enallagma cyathigerum (Odonata: Coenagrionidae) for habitats of varying structural complexity. Eur. J. Entomol. 109, 229–234 (2012).Article 

Google Scholar 
McCarty, J. P., Wolfenbarger, L. L. & Wilson, J. A. eLS 1–13 (Wiley, 2017).Book 

Google Scholar 
Holzmann, K. L. Challenges in a Changing Climate: The Effect of Temperature Variation on Growth and Competition in Damselflies Independent thesis Advanced level (degree of Master (Two Years) thesis, Uppsala University. http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-467582 (2022). More

Hughes, W. O. H., Oldroyd, B. P., Beekman, M. & Ratnieks, F. L. W. Ancestral monogamy shows kin selection is key to the evolution of eusociality. Science 320, 1213–1216 (2008).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Ross, K. G. & Matthews, R. W. The Social Biology of Wasps (Cornell University Press, 1991).Book 

Google Scholar 
Itô, Y. & Higashi, S. Spring behaviour of Ropalidia plebeiana (Hymenoptera: Vespidae) within a huge aggregation of nests. Appl. Entomol. Zool. 22, 519–527 (1987).Article 

Google Scholar 
Saito, F. & Kojima, J.-I. Colony cycle in the south-eastern coastal populations of Ropalidia plebeiana, the only Ropalidia wasp occurring in temperate Australia. Entomol. Sci. 8, 263–275 (2005).Article 

Google Scholar 
Richards, O. W. The Australian social wasps (Hymenoptera: Vespidae). Aust. J. Zool. Suppl. Ser. 26, 1–132 (1978).Article 

Google Scholar 
Makino, S., Yamane, S., Itô, Y. & Spradbery, J. P. Process of comb division of reused nests in the Australian paper wasp Ropalidia plebeinana (Hymenoptera, Vespidae). Ins. Soc. 41, 411–422 (1994).Article 

Google Scholar 
Goodnight, K. F. & Queller, D. C. Computer software for performing likelihood tests of pedigree relationship using genetic markers. Mol. Ecol. 8, 1231–1234 (1999).Article 

Google Scholar 
Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 14, 2611–2620 (2005).CAS 
PubMed 
Article 

Google Scholar 
Crochet, P.-A. Genetic structure of avian populations—Allozymes revisited. Mol. Ecol. 9, 1463–1469 (2000).CAS 
PubMed 
Article 

Google Scholar 
Meirmans, P. G. & Hedrick, P. W. Assessing population structure: FST and related measures. Mol. Ecol. Res. 11, 5–18 (2011).Article 

Google Scholar 
Hedrick, P. W. A standardized genetic differentiation measure. Evolution 59, 1633–1638 (2005).CAS 
PubMed 
Article 

Google Scholar 
Itô, Y., Yamane, S. & Spradbery, J. P. Population consequences of huge nesting aggregations of Ropalidia plebeiana (Hymenoptera: Vespidae). Res. Popul. Ecol. 30, 279–295 (1988).Article 

Google Scholar 
Boomsma, J. J. & d’Ettorre, P. Nice to kin and nasty to non-kin: revisiting Hamilton’s early insights on eusociality. Biol. Lett. 9, 20130444 (2013).PubMed 
PubMed Central 
Article 

Google Scholar 
Hannonen, M. & Sundström, L. Worker nepotism among polygynous ants. Nature 421, 910 (2003).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Parsons, P. J., Grinsted, L. & Field, J. Partner choice correlates with fine scale kin structuring in the paper wasp Polistes dominula. PLoS ONE 14, e0221701 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Strassmann, J. E., Queller, D. C., Solis, C. R. & Hughes, C. R. Relatedness and queen number in the Neotropical wasp, Parachartergus colobopterus. Anim. Behav. 42, 461–470 (1991).Article 

Google Scholar 
Leadbeater, E., Carruthers, J. M., Green, J. P., Rosser, N. S. & Field, J. Nest inheritance is the missing source of direct fitness in a primitively eusocial insect. Science 333, 874–876 (2011).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Field, J. & Leadbeater, E. Cooperation between non-relatives in a primitively eusocial paper wasp, Polistes dominula. Philos. Trans. R. Soc. B 371, 20150093 (2016).Article 

Google Scholar 
Bhadra, A. & Gadagkar, R. We know that the wasps “know”: cryptic successors to the queen in Ropalidia marginata. Biol. Lett. 4, 634–637 (2008).PubMed 
PubMed Central 
Article 

Google Scholar 
Bang, A. & Gadagkar, R. Reproductive queue without overt conflict in the primitively eusocial wasp Ropalidia marginata. Proc. Nat. Acad. Sci. USA 109, 14494–14499 (2012).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Frank, S. A. Hierarchical selection theory and sex ratios. II. On applying the theory, and a test with fig wasps. Evolution 39, 949–964 (1985).PubMed 
Article 

Google Scholar 
Silk, J. B. & Brown, G. R. Local resource competition and local resource enhancement shape primate birth sex ratios. Proc. R. Soc. Lond. B. 275, 1761–1765 (2008).
Google Scholar 
Schwarz, M. P. Local resource enhancement and sex ratios in a primitively social bee. Nature 331, 346–348 (1988).ADS 
Article 

Google Scholar 
Cronin, A. L. & Schwarz, M. P. Sex ratios, local fitness enhancement and eusociality in the allodapine bee Exoneura richardsoni. Evol. Ecol. 11, 567–577 (1997).Article 

Google Scholar 
Schwarz, M. P., Bull, N. J. & Hogendoorn, K. Evolution of sociality in the allodapine bees: A review of sex allocation, ecology and evolution. Ins. Soc. 45, 349–368 (1998).Article 

Google Scholar 
Gamboa, G. J., Wacker, T. L., Duffy, K. G., Dobson, S. W. & Fishwild, T. G. Defence against intraspecific usurpation by paper wasp cofoundresses (Polistes fuscatus, Hymenoptera: Vespidae). Can J. Zool. 70, 2369–2372 (1992).Article 

Google Scholar 
Katada, S. & Iwahashi, O. Characteristics of usurped colonies in the subtropical paper wasp, Ropalidia fasciata (Hymenoptera: Vespidae). Ins. Soc. 43, 247–253 (1996).Article 

Google Scholar 
Yamane, S. Ecological factors influencing the colony cycle of Polistes wasps. in Natural History and Evolution of Paper-Wasps (Turillazzi, S. & West-Eberhard, M. J. eds.). 75–97. (Oxford University Press, 1996).Clouse, R. Some effects of group size on the output of beginning nests of Mischocyttarus mexicanus (Hymenoptera: Vespidae). Flor. Entomol. 84, 418–425 (2001).Article 

Google Scholar 
Strassmann, J. E. Female-biased sex ratios in social insects lacking morphological castes. Evolution 38, 256–266 (1984).PubMed 

Google Scholar 
Suzuki, T. Production schedule of males and reproductive females, investment sex ratios, and worker-queen conflict in paper wasps. Am. Nat. 128, 366–378 (1986).Article 

Google Scholar 
Tsuchida, K. & Suzuki, T. Conflict over sex ratio and male production in paper wasps. Ann. Zool. Fenn. 43, 468–480 (2006).
Google Scholar 
Ohtsuki, H. & Tsuji, K. Adaptive reproduction schedule as a cause of worker policing in social Hymenoptera: A dynamic game analysis. Am. Nat. 173, 747–758 (2009).PubMed 
Article 

Google Scholar 
Walsh, P. S., Metzger, D. A. & Higuchi, R. Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques 10, 506–513 (1991).CAS 
PubMed 

Google Scholar 
Bassam, B. J., Caetano-Anolles, G. & Gresshoff, P. M. Fast and sensitive silver staining of DNA in polyacrylamide gels. Anal. Biochem. 196, 80–83 (1991).CAS 
PubMed 
Article 

Google Scholar 
van Oosterhout, C., Hutchinson, W. F., Wills, D. P. M. & Shipley, P. MICRO-CHECHER: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes. 4, 535–538 (2004).Article 
CAS 

Google Scholar 
Peakall, R. & Smouse, P. E. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes. 6, 288–295 (2006).Article 

Google Scholar 
Raymond, M. & Rousset, F. GENEPOP (version 1.2): Population genetics software for exact tests and ecumenicism. J. Hered. 86, 248–249 (1995).Article 

Google Scholar 
Meirmans, P. G. GenoDive version 3.0: Easy-to-use software for the analysis of genetic data of diploids and polyploids. Mol. Ecol. Res. 20, 1126–1131 (2020).CAS 
Article 

Google Scholar 
Michalakis, Y. & Excoffier, L. A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics 142, 1061–1064 (1996).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945–959 (2000).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Earl, D. & vonHoldt, B. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Res. 4, 359–361 (2012).Article 

Google Scholar 
Kopelman, N. M., Mayzel, J., Jakobsson, M., Rosenberg, N. A. & Mayrose, I. Clumpak: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Res. 15, 1179–1191 (2015).CAS 
Article 

Google Scholar 
Goodnight, K. F. Relatedness 4.2c Release (Rice University, 1996).
Google Scholar 
Queller, D. C. A method for detecting kin discrimination within natural colonies of social insects. Anim. Behav. 47, 569–576 (1994).Article 

Google Scholar  More

By pooling the results of the entire known range analysis of 14 GPS-collared elephants living in the Kinabatangan, our study suggests that this populations range covers at least 628 km2 (Table 3). Nine different locations were identified as hot spots, representing 266.9 km2 or 43% of this range, suggesting that just under half is highly used and/or frequented (Fig. 1). We found that the size of individual’s hot spots was positively related to the size of the entire range, meaning the larger the entire range the larger the summed area of an elephants hot spots. On average, hot spots represented a relatively small percent of an animal’s entire range (ranging from 4 to 20%, averaging 12%, Table 3). However, time spent within these hot spots ranged from 10 to 60% (averaging 34% across elephants, Table 5), with time spent in hot spots being related to the overall size of the hot spots (the larger the hot spot the more time elephants spent in them).Identifying the location of these hot spots is essential in designing appropriate management practices in collaboration with land users and identifying the best location for elephant corridors. In the last 25 years, forest cover in the Lower Kinabatangan has been drastically reduced and fragmented46, eroding the biodiversity value of this landscape. Today, this region has little remaining forests, and what is left is insufficient for sustaining the local elephant population10. Moreover, forests are highly fragmented along the Kinabatangan River, with a number of bottlenecks constraining elephant movements9. The situation in this landscape is getting worse because of further land clearances for agriculture, namely oil palm; as well as for the highly controversial Sukau Bridge and new road/highway that is planned for the region.Our analyses revealed a highly significant difference between the average proportions of protected area, unprotected forest, and oil palm estate extents within the elephant’s entire range; and a substantive, but not significant, difference across these land use/land cover types within hot spots (Table SI 4). At the individual level, there was a highly significant negative relationship between the proportion of protected areas and oil palm estates both within the elephant’s entire range and within the hot spots.At the pooled level, we found that around 45% of the entire known range and hot spots were within forested environments (280.44 km2 and 120.29 km2 respectively). Our results showed strong fidelity of certain elephants to these forested habitats. Our k-means cluster analysis found that within elephant entire ranges and hot spots, two out of the three cluster groups had high or very high usage of forests. Both cluster 1, for the entire range, and cluster 1 for hot spots extents, had five females that on average used forest environments 90% of their time, with protected areas being used 64% and 59%, and unprotected forested being used on average 26% and 31%, respectively (Table 7).Individuals in cluster 2, for the entire range analysis, on average, spent 73% of their time in forests (57% of this in protected areas and 16% in unprotected forests; Table 7). For the hot spot analysis, the individuals in cluster 2 spent on average 65% of their time in forests (52% of this in the unprotected forests and 13% in protected forests; Table 7). Elephants within these clusters were all females. Our results suggest that forest may be of particular importance for females as they had forest as their dominant land cover type within their entire range, hot spot extents and time spent analyses (Fig. 3, Table 5). Several studies have shown that adult females influence and guide the movement patterns and habitat utilization for their family group and that females in family units tend to inhabit less risky areas, such as within natural forest habitat60,61,62.However, the unprotected forest is at risk. We identified about 8% (or 49 km2) of forest identified within the pooled entire known range were not protected, with half potentially being on state land, and the remaining half on land titles of various types (Table SI 4). For the pooled hot spot areas, unprotected forest was proportionally higher, comprising of 11% (or 29 km2) of the total extent, with 54% being potentially on State land and 46% on land titles (Table SI 4). Protecting these forests would be an essential and efficient way to secure key elephant habitat since all collared individuals were using these forest fragments in their entire range (averaging 11%, and ranging from 8 to 18%), and hot spot extents (averaging 20%, and ranging from 0 to 53%) (Table SI 4, Fig. 3). On average, 24% of time was spent in unprotected forests within hot spots, though this varied widely from 0% (for the male elephant known as Gading) to 61% (for the female matriarch named Jasmine) (Table 5). In fact, five females had large proportions of their hot spot extents (24–53%) in unprotected forests, spending substantial periods of their time (33–61%) within these threatened areas.Our findings show that unprotected forests around the villages of Bilit and Sukau, were of particular significance (Figs. 1, 2). These unprotected forests largely consist of lowland dry forest, seasonally flooded swamp forest, and swamp forest, which are considered important habitats for elephants for feeding, resting and moving47,63. Within these forests, and along the forest margins and river banks there are also natural open grasslands that consist of Phragmites karka and Dinochloa scabrida that provide essential forage, mainly in the riparian areas for elephants9,21,23. Forested environments are also considered to be important in providing natural refugee from human activities and disturbance. For example, elephants have been documented to form significantly larger group sizes, as well as engaging in significantly more social interactions, in natural forest habitat compared to, for example, oil palm landscapes63. Adult females, generally, avoid areas considered unsafe for their respective social units, are more selective in the resources they use, and require regular access to water because of the presence of young64,65,66. This may be why our results, strongly suggest that forest habitats seem to be most important for adult females.Another significant issue faced by these elephants is the threat from the controversial planned Sukau bridge and road/highway that is set out in the Sabah Structure Plan, an overarching policy document for the State58. Currently, a new road/highway is under construction on the northern bank of the village of Sukau, and this has already cleared areas of unprotected forest. This public road could link to a potential new bridge that would cross over the Kinabatangan River, cutting through unprotected forest and a protected area (Lower Kinabatangan Wildlife Sanctuary), before going through oil palm estates then through another protected area to the south and through the Tabin elephant population range. For the Kinabatangan, creating a public highway will cut the elephant population range into two parts (Figs. 2, 3). All collared elephants use this area, as it is a key bottleneck and the only alternative option to pass around Sukau village9. We found that nine elephants have hot spots that intersect or meet up with the current road (which will be up-graded and get considerably busier) and/or the planned road/highway alignment (Figs. SI 1 and 2). For these elephants, we calculated that they spent from 2 to 44% (average 14%) of their time within these hot spots (Table 4). Our statistical analyses suggest that if the road/highway goes ahead it will have a significant impact on the elephants’ behaviour with respect to time spent in the hot spots. Indeed, this infrastructure project could have dire consequences for these elephants and their family groups, by disrupting their ranging patterns and segmenting the entire elephant range into two (Figs. 2, 4). The existing road in Batu Putih has already proven to be an impassable barrier for this elephant population, as no elephants have been observed crossing this road since the early 2000s14. For elephants that do try and cross, vehicle collisions may become a significant threat to elephants and drivers alike67, and potentially increasing human–elephant conflict in the nearby villages, as well as in plantations14,68,69, exacerbating an already difficult situation for this small and fragmented population.Results from the pooled analysis show that about 53% of the entire known population range is within oil palm estates; and 51% for the pooled hot spots (Fig. 3, Table SI 4). Our k-means clustering analysis grouped 6 elephants into cluster 3 that on average spent 57% of time in oil palm estates; and 7 elephants into cluster 2 within the hot spot analysis that on average spent 73% of their time in oil palm estates (Table 6). All the males, were clustered within these groups (Table 5). In fact, the three collared males were amongst the highest users of oil palm estates (Fig. 3, Table SI4, and 5). This could be related to a ‘‘high risk, high gain’’ strategy, often adopted by males to increase body size and enhance reproductive success32,33,60. However, it is interesting to see that three females (Ita, Ratu and Koyah) and their respective social units, also seemed to have high levels of oil palm use, while other individuals had zero or very little use of oil palm (e.g. Aqeela, Jasmin, Sandi, Kasih; Table SI 4, Fig. 3). Differential choices may result from differences in individual knowledge and experience with people during past encounters, for example70,71.We identified that collared elephants were ranging in 11 known oil palm estates, with the five most regularly used being Melangking Oil Palm Plantation (with 12 elephants entire range overlapping with this estate and six hot spots), IOI Corporation (with 11 overlapping entire ranges, and eight hot spots), Genting Plantations (14 and seven, respectively), Sime Darby Plantation (five and two, respectively), and Karangan Agriculture (8 and 2, respectively) (Table 6; Fig. 4). Presence of bottlenecks and barriers (e.g. electric fences) may explain hot spot occurrences in these estates, as well as feeding opportunities, management strategies of specific estates, and historical and seasonal ranges.Linear features like major highways, electric fences and drainage ditches hamper elephant movements within the Lower Kinabatangan9. A previous study identified 20 bottlenecks in the Lower Kinabatangan with the two main ones (of 9 km and 6.5 km in length) found around the village of Sukau9. In addition, the unplanned and chaotic erection of electric fences by large estates and smallholdings has disrupted significantly elephant movement patterns and resulted in artificial hot spots for certain individuals (e.g. Liun, Ita, Gading and Sejati)35,72. Electric fences have widely been used to mitigate human–elephant conflicts. The establishment of fences rarely consider the traditional elephant routes nor the location of existing fences in neighbouring estates. If elephants manage to enter such areas, they often become trapped and experience difficulties in returning to nearby forests, exacerbating conflicts with people35.Certain estates such as Melangking Oil Palm Plantation have allowed elephants to roam freely in their estate (Muhammad Al-Shafieq, personal communication). Since 2017, this plantation has shown a drastic reduction in damages to their oil palms following the removal of their permanent electric fences surrounding their entire estate. Instead, this plantation is using a temporary electric fencing regime around newly planted palm areas. Concurrently, they now do not push elephants out of their estate, which can explain why Melangking Oil Palm Plantation is a significant hotspot in the region.Another reason why elephant ranges incorporate oil palm estates is to move between forest patches that are becoming completely isolated following forest conversion, as is the case close to Sukau (Fig. SI1 and SI2; Fig. 1). Unlike other elephant species that increase their speed of movement rates in highly disturbed areas27,30,66, the Bornean elephant has been observed doing the opposite, which may explain some of the hot spots within oil palm estates. This movement strategy may allow for better vigilance as seen on a few occasions when elephants spent 2–5 days in the Bukit Melapi-Yu Kwang Corridor, near the village of Sukau, before leaving the area (Othman, personal observation).Hot spots in the oil palm landscape can also be explained by feeding opportunities, since elephants feed on palm shoots, leaves and hearts73. Elephants are known to eat the shoots of newly planted oil palms, often killing the palms and causing significant economic damages35. Since 2010, many estates located in the Lower Kinabatangan have started a new palm rotation. Palms are replanted every 25 years. A new rotation includes land clearing, bole and root mass removal, and the shredding or chipping of felled palms. Elephants are attracted to the shredded palm hearts since it gives them easy access to one of their favourite food72. This particular behaviour does not cause economic damage, and some estate managers allow the elephants to stay and forage in the chipping areas. This was documented for several collared elephants, whose hot spots and time spent were particularly high within oil palm (e.g. Gading and Sandy, two males; and Ratu and Ita, two females). Once the shredded palms have dried, however, elephants will leave these areas and move elsewhere. Within oil palm estates, some elephants have been found to travel more directly and rapidly suggesting ‘exploratory’ behaviour, which could be associated with searching for young palms or areas of palm felling and chipping of palm hearts15.Lastly, elephants may still be using their historical range that used to be covered with forest before conversion to oil palm. Other factors potentially explaining the relatively high use of oil palm estates include seasonal variations of ranging patterns. Indeed events of drought or floods limit the access to various parts of the floodplain and will tend to confine the animals in some areas9,63.In Sabah the state authorities have recorded at least 200 elephant deaths from the year 2010 to 2021 and most of these have occurred on, or near, oil palm estates14,74,75,76. Deaths from non-natural causes are largely due to poisoning (both accidental and intentional), gunshot wounds, poaching for tusks and other body parts, and snares35. Stopping killing and enabling a safe coexistence between people and elephants within multiple-use landscapes that are dominated by oil palm is one of the key strategies developed in the Bornean Elephant Action Plan for Sabah (2020–2029), which was endorsed by the State14. Based on our results in Lower Kinabatangan, a series of recommendations are proposed.This study underscores the importance of remaining forested areas for the Lower Kinabatangan elephant population. Full protection of all forest fragments left in the Lower Kinabatangan is urgently needed. Several official mechanisms are available to fulfil this request that has been proposed for the past 20 years by various organizations46.The current network of forests available in the Lower Kinabatangan is too small and fragmented to sustain a viable elephant population. Forest corridors must be created across the landscape through reforestation exercises, whilst concurrently undertaking enrichment planting of native understory forage within forested areas as this may minimize the need for elephants to search for easily accessible food in high-risk oil palm landscapes21,22,23.Current governmental plans to build a road bridge and public road/highway linking the southern bank of the Kinabatangan River to Tabin Wildlife Reserve to the south will irreversibly impact the Lower Kinabatangan elephant population by cutting the current range into two isolated parts. This will impact the elephants ranging patterns, potentially even fragmenting the already small population into two groups, and potentially leading to elephant deaths by vehicle collisions (which is becoming increasingly common in Peninsular Malaysia), and increase the risk of poaching activities, all resulting in a decrease in the genetic diversity of the, already small and isolated, population14,67.Eventually, the future of the Kinabatangan elephant population resides in improving land use and management practices within oil palm estates currently used by elephants. We recommend that priority should be given at improving elephant movements in oil palm estates by removing unnecessary man-made barriers and only cautiously installing temporary electric fences to protect sensitive areas. For example, the use of electric fences around mature oil palm and areas whereby palms are being removed and chipped could be prohibited, and electric fences permitted solely for protecting oil palm nurseries, new plantings and young oil palms (e.g. up to 7–8 years old), and staff and office quarters. This would greatly allow for landscape permeability for elephants, and other species that need to cross the landscape for their ecological and biological needs14.A handful of guidelines exist to assist oil palm managers and staff in managing elephant populations in their respective estates72,77. However, there is a need for a more comprehensive set of guidelines, which delineate better practices with the aim to increase the protection of people and elephants outside protected areas. Guidelines should specify “do’s” and “don’ts” (based on best available data and knowledge) of actions needed before, during and after elephants visit oil palm estates and smallholdings.Sabah now is in an interesting transition within their palm oil sector. On the 21st October 2015, the Sabah State Cabinet committed to produce 100% certified sustainable palm oil, by 2025, under the Roundtable for Sustainable palm Oil (RSPO) Jurisdictional Certification approach. Under this approach, areas of High Conservation Value and areas identified within the High Carbon Stock Approach need specific management and monitoring, in order to comply with RSPO principles and criteria78,79,80. Sabah government can use this platform to build an integrated landscape level approach to better manage landscapes within known elephant ranges (which is considered a High Conservation Value species) to allow for a safe and permeable movement through the landscape.Eventually, long-term survival of the Bornean elephant will mainly depend on how people and elephants can co-exist. It is our hope that this study illustrates the importance of protecting all forested habitat and effectively managing areas outside of protected areas to allow for long-term elephant coexistence with humans in this landscape. More

Paerl, H. W. & Huisman, J. Climate change: a catalyst for global expansion of harmful cyanobacterial blooms. Environ. Microbiol. Rep. 1, 27–37 (2009).CAS 
PubMed 
Article 

Google Scholar 
Paerl, H. W. & Paul, V. J. Climate change: links to global expansion of harmful cyanobacteria. Water Res. 46, 1349–1363 (2012).CAS 
PubMed 
Article 

Google Scholar 
Huisman, J. et al. Cyanobacterial blooms. Nat. Rev. Microbiol. 16, 471 (2018).CAS 
PubMed 
Article 

Google Scholar 
Reissig, M., Trochine, C., Queimaliños, C., Balseiro, E. & Modenutti, B. Impact of fish introduction on planktonic food webs in lakes of the Patagonian Plateau. Biol. Conserv. 132, 437–447 (2006).Article 

Google Scholar 
Britton, J. R., Davies, G. D. & Harrod, C. Trophic interactions and consequent impacts of the invasive fish Pseudorasbora parva in a native aquatic foodweb: a field investigation in the UK. Biol. Invasions 12, 1533–1542 (2010).Article 

Google Scholar 
Beaulieu, M., Pick, F. & Gregory-Eaves, I. Nutrients and water temperature are significant predictors of cyanobacterial biomass in a 1147 lakes data set. Limnol. Oceanogr. 58, 1736–1746 (2013).ADS 
CAS 
Article 

Google Scholar 
O’Neil, J. M., Davis, T. W., Burford, M. A. & Gobler, C. J. The rise of harmful cyanobacteria blooms: the potential roles of eutrophication and climate change. Harmful Algae 14, 313–334 (2012).Article 
CAS 

Google Scholar 
Paerl, H. W. & Huisman, J. Blooms like it hot. Science 320, 57–58 (2008).CAS 
PubMed 
Article 

Google Scholar 
Toxic Cyanobacteria in Water: A Guide to Their Public Health Consequences, Monitoring, and Management. (E & FN Spon, 1999).Sukenik, A., Quesada, A. & Salmaso, N. Global expansion of toxic and non-toxic cyanobacteria: effect on ecosystem functioning. Biodivers. Conserv. 24, 889–908 (2015).Article 

Google Scholar 
Ibelings, B. W., Bormans, M., Fastner, J. & Visser, P. M. CYANOCOST special issue on cyanobacterial blooms: synopsis—a critical review of the management options for their prevention, control and mitigation. Aquat. Ecol. 50, 595–605 (2016).CAS 
Article 

Google Scholar 
Paerl, H. W. Mitigating harmful cyanobacterial blooms in a human- and climatically-impacted world. Life 4, 988–1012 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Rastogi, R. P., Madamwar, D. & Incharoensakdi, A. Bloom dynamics of cyanobacteria and their toxins: environmental health impacts and mitigation strategies. Front. Microbiol. https://doi.org/10.3389/fmicb.2015.01254 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Ewing, H. A. et al. “New” cyanobacterial blooms are not new: two centuries of lake production are related to ice cover and land use. Ecosphere 11, e03170 (2020).Article 

Google Scholar 
McGlone, M. S. & Wilmshurst, J. M. Dating initial Maori environmental impact in New Zealand. Quat. Int. 59, 5–16 (1999).Article 

Google Scholar 
Brooking, A. P. D. of H. T. & Brooking, T. The History of New Zealand. (Greenwood Publishing Group, 2004).Wilmshurst, J. M., Anderson, A. J., Higham, T. F. G. & Worthy, T. H. Dating the late prehistoric dispersal of Polynesians to New Zealand using the commensal Pacific rat. PNAS 105, 7676–7680 (2008).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
McGlone, M. S. The Polynesian settlement of New Zealand in relation to environmental and biotic changes. N. Z. J. Ecol. 12, 115–129 (1989).
Google Scholar 
McGlone, M. S. Polynesian deforestation of New Zealand: a preliminary synthesis. Archaeol. Ocean. 18, 11–25 (1983).Article 

Google Scholar 
McWethy, D. B. et al. Rapid landscape transformation in South Island, New Zealand, following initial Polynesian settlement. PNAS 107, 21343–21348 (2010).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
McWethy, D. B., Wilmshurst, J. M., Whitlock, C., Wood, J. R. & McGlone, M. S. A high-resolution chronology of rapid forest transitions following Polynesian arrival in New Zealand. PLoS ONE 9, e111328 (2014).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Star, P. New Zealand environmental history: a question of attitudes. Environ. Hist. Camb. 9, 463–475 (2003).Article 

Google Scholar 
Clark, A. H. The Invasion of New Zealand by Plants, People, and Animals (Rutgers University Press, 1949).
Google Scholar 
Wilmshurst, J. M. Human effects on the environment: European impact. Te Ara: The Encyclopedia of New Zealand https://teara.govt.nz/en/human-effects-on-the-environment/page-3 (2007).Smol, J. P. The ratio of diatom frustules to chrysophycean statospores: a useful paleolimnological index. Hydrobiologia 123, 199–208 (1985).Article 

Google Scholar 
Rees, A. B. H., Cwynar, L. C. & Cranston, P. S. Midges (Chironomidae, Ceratopogonidae, Chaoboridae) as a temperature proxy: a training set from Tasmania, Australia. J. Paleolimnol. 40, 1159–1178 (2008).ADS 
Article 

Google Scholar 
Epp, L. S., Stoof, K. R., Trauth, M. H. & Tiedemann, R. Historical genetics on a sediment core from a Kenyan lake: intraspecific genotype turnover in a tropical rotifer is related to past environmental changes. J. Paleolimnol. 43, 939–954 (2010).ADS 
Article 

Google Scholar 
Buchaca, T. et al. Rapid ecological shift following piscivorous fish introduction to increasingly eutrophic and warmer Lake Furnas (Azores Archipelago, Portugal): a paleoecological approach. Ecosystems 14, 458–477 (2011).CAS 
Article 

Google Scholar 
Cristescu, M. E. & Hebert, P. D. N. Uses and misuses of environmental DNA in biodiversity science and conservation. Annu. Rev. Ecol. Evol. Syst. 49, 209–230 (2018).Article 

Google Scholar 
Giguet-Covex, C. et al. Long livestock farming history and human landscape shaping revealed by lake sediment DNA. Nat. Commun. 5, 3211 (2014).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Alsos, I. G. et al. Plant DNA metabarcoding of lake sediments: how does it represent the contemporary vegetation. PLoS ONE 13, e0195403 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Nelson-Chorney, H. T. et al. Environmental DNA in lake sediment reveals biogeography of native genetic diversity. Front. Ecol. Environ. 17, 313–318 (2019).
Google Scholar 
Capo, E. et al. Lake sedimentary DNA research on past terrestrial and aquatic biodiversity: overview and recommendations. Quaternary 4, 6 (2021).Article 

Google Scholar 
Shokralla, S., Spall, J. L., Gibson, J. F. & Hajibabaei, M. Next-generation sequencing technologies for environmental DNA research. Mol. Ecol. 21, 1794–1805 (2012).CAS 
PubMed 
Article 

Google Scholar 
Taberlet, P., Coissac, E., Pompanon, F., Brochmann, C. & Willerslev, E. Towards next-generation biodiversity assessment using DNA metabarcoding. Mol. Ecol. 21, 2045–2050 (2012).CAS 
PubMed 
Article 

Google Scholar 
Thomsen, P. F. & Willerslev, E. Environmental DNA: an emerging tool in conservation for monitoring past and present biodiversity. Biol. Conserv. 183, 4–18 (2015).Article 

Google Scholar 
Keeley, N., Wood, S. A. & Pochon, X. Development and preliminary validation of a multi-trophic metabarcoding biotic index for monitoring benthic organic enrichment. Ecol. Ind. 85, 1044–1057 (2018).CAS 
Article 

Google Scholar 
Monchamp, M.-E., Walser, J.-C., Pomati, F. & Spaak, P. Sedimentary DNA reveals cyanobacterial community diversity over 200 years in two perialpine lakes. Appl. Environ. Microbiol. 82, 6472–6482 (2016).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pal, S., Gregory-Eaves, I. & Pick, F. R. Temporal trends in cyanobacteria revealed through DNA and pigment analyses of temperate lake sediment cores. J. Paleolimnol. 54, 87–101 (2015).ADS 
Article 

Google Scholar 
Dodsworth, W. Temporal Trends in Cyanobacteria Through Paleo-Genetic Analyses. (Université d’Ottawa/University of Ottawa, 2020). https://doi.org/10.20381/ruor-24401.Rinta-Kanto, J. M. et al. The diversity and distribution of toxigenic Microcystis spp. in present day and archived pelagic and sediment samples from Lake Erie. Harmful Algae 8, 385–394 (2009).CAS 
Article 

Google Scholar 
Zastepa, A. et al. Reconstructing a long-term record of microcystins from the analysis of lake sediments. Sci. Total Environ. 579, 893–901 (2017).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Schallenberg, M. et al. Ecosystem services of lakes. In Ecosystem Services in New Zealand (ed. Dymond, J.) 23 (Manaaki Whenua Press, 2013).
Google Scholar 
Ministry for the Environment & Ausseil, A.-G. Our freshwater 2020. www.mfe.govt.nz (2020).Takiwa—Map Page. https://lernz.takiwa.co/map.Leathwick, J. et al. Freshwater ecosystems of New Zealand (FENZ) geodatabase. Users guide. (2010).Cochrane, L. Reconstructing Ecological Change, Catchment Disturbance, and Anthropogenic Impact over the last 3000 years at Lake Pounui, Wairarapa, New Zealand. (2017).Burns, C. W. & Mitchell, S. F. Seasonal succession and vertical distribution of phytoplankton in Lake Hayes and Lake Johnson, South Island, New Zealand. N. Z. J. Mar. Freshw. Res. 8, 167–209 (1974).Article 

Google Scholar 
Lawa. Land, Air, Water Aotearoa (LAWA) https://www.lawa.org.nz/ (2018).Bunny, T., Perrie, A., Milne, J. & Keenan, L. Lake water quality in the Ruamāhanga Whaitua. 17 (2014).McKinnon, M. Volcanic Plateau region: The lure of trout. Te Ara—The Encyclopedia of New Zealand https://teara.govt.nz/en/volcanic-plateau-region/page-8 (2015).Burns, C. W. & Mitchell, S. F. Seasonal succession and vertical distribution of zooplankton in Lake Hayes and Lake Johnson. N. Z. J. Mar. Freshw. Res. 14, 189–204 (1980).Article 

Google Scholar 
Schallenberg, M. & Schallenberg, L. Lake Hayes restoration and monitoring plan. 55 https://a234f952-dbf2-444e-983e-ef311d984ee7.filesusr.com/ugd/c1b10b_d2993ed023cd4bdbac7eef71a89c2de7.pdf (2017).NIWA. NIWA https://niwa.co.nz/.Mackereth, F. J. H. A portable core sampler for lake deposits. Limnol. Oceanogr. 3, 181–191 (1958).ADS 
Article 

Google Scholar 
Howarth, J. D., Fitzsimons, S. J., Norris, R. J. & Jacobsen, G. E. Lake sediments record cycles of sediment flux driven by large earthquakes on the Alpine fault, New Zealand. Geology 40, 1091–1094 (2012).ADS 
CAS 
Article 

Google Scholar 
Trodahl, M. I., Rees, A. B. H., Newnham, R. M. & Vandergoes, M. J. Late Holocene geomorphic history of Lake Wairarapa, North Island, New Zealand. N. Z. J. Geol. Geophys. 59, 330–340 (2016).CAS 
Article 

Google Scholar 
Khan, S., Puddick, J., Burns, C. W., Closs, G. & Schallenberg, M. Palaeolimnological evaluation of historical nutrient and food web contributions to the eutrophication of two monomictic lakes. Submitted for Journal Publication (2022).Rinta-Kanto, J. M. et al. Quantification of toxic Microcystis spp. during the 2003 and 2004 blooms in Western Lake Erie using quantitative real-time PCR. Environ. Sci. Technol. 39, 4198–4205 (2005).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Nübel, U., Garcia-Pichel, F. & Muyzer, G. PCR primers to amplify 16S rRNA genes from cyanobacteria. Appl. Environ. Microbiol. 63, 3327–3332 (1997).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, 2020).RStudio Team. RStudio: Integrated Development for R. RStudio, PBC, Boston, MA. (2020).Wickham, H. et al. Welcome to the Tidyverse. J. Open Sour. Softw. 4, 1686 (2019).ADS 
Article 

Google Scholar 
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag, 2016).MATH 
Book 

Google Scholar 
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17, 10–12 (2011).Article 

Google Scholar 
Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).CAS 
PubMed 
Article 

Google Scholar 
Yilmaz, P. et al. The SILVA and “all-species living tree project (LTP)” taxonomic frameworks. Nucleic Acids Res. 42, D643–D648 (2013).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Glöckner, F. O. et al. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J. Biotechnol. 261, 169–176 (2017).PubMed 
Article 
CAS 

Google Scholar 
McMurdie, P. J. & Holmes, S. phyloseq: an R Package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 8, e61217 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Oksanen, J. et al. vegan: Community Ecology Package. R package version 2.5–6. 2019. (2019).Williams, P. A. & Cameron, E. K. Creating gardens: the diversity and progression of European plant introductions. In Biological Invasions in New Zealand Vol. 186 (eds Allen, R. B. & Lee, W. G.) 33–47 (Springer-Verlag, 2006).Chapter 

Google Scholar 
Simpson, G. L. Modelling palaeoecological time series using generalised additive models. Front. Ecol. Evol. https://doi.org/10.3389/fevo.2018.00149 (2018).Article 

Google Scholar 
Chen, H. & Boutros, P. C. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinform. 12, 35 (2011).Article 

Google Scholar 
Juggins, S. rioja: analysis of quaternary science data. (2020).de Vries, A. & Ripley, B. D. ggdendro: create dendrograms and tree diagrams using ‘ggplot2’. (2022).Anderson, M. J. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46 (2001).
Google Scholar 
Anderson, M. J. Distance-based tests for homogeneity of multivariate dispersions. Biometrics 62, 245–253 (2006).MathSciNet 
PubMed 
MATH 
Article 

Google Scholar 
Soo, R. M. et al. An expanded genomic representation of the phylum Cyanobacteria. Genome Biol. Evol. 6, 1031–1045 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
MacKeigan, P. W. et al. Comparing microscopy and DNA metabarcoding techniques for identifying cyanobacteria assemblages across hundreds of lakes. Harmful Algae 113, 102187 (2022).CAS 
PubMed 
Article 

Google Scholar 
Wood, S. A. et al. Trophic state and geographic gradients influence planktonic cyanobacterial diversity and distribution in New Zealand lakes. FEMS Microbiol. Ecol. https://doi.org/10.1093/femsec/fiw234 (2017).Article 
PubMed 

Google Scholar 
Becker, S., Richl, P. & Ernst, A. Seasonal and habitat-related distribution pattern of Synechococcus genotypes in Lake Constance. FEMS Microbiol. Ecol. 62, 64–77 (2007).CAS 
PubMed 
Article 

Google Scholar 
Sánchez-Baracaldo, P., Handley, B. A. & Hayes, P. K. Picocyanobacterial community structure of freshwater lakes and the Baltic Sea revealed by phylogenetic analyses and clade-specific quantitative PCR. Microbiology (Reading) 154, 3347–3357 (2008).Article 
CAS 

Google Scholar 
Pilon, S. et al. Contrasting histories of microcystin-producing cyanobacteria in two temperate lakes as inferred from quantitative sediment DNA analyses. Lake Reserv. Manag. 35, 102–117 (2019).CAS 
Article 

Google Scholar 
Queenstown’s Pioneering Beginnings. https://www.queenstownnz.co.nz/stories/post/queenstowns-pioneer-beginnings/ (2017).Fish, G. R. A limnological study of four lakes near Rotorua. N. Z. J. Mar. Freshw. Res. 4, 165–194 (1970).Article 

Google Scholar 
Lake Rotoehu—Lakes Water Quality Society. https://lakeswaterquality.co.nz/lake-rotoehu/.Bay of Plenty Regional Council, Rotorua District Council, & Te Arawa Lakes Trust. Lake Rotoehu Action Plan. 61 http://www.rotorualakes.co.nz/vdb/document/76 (2007).Hobbs, W. O. et al. Using a lake sediment record to infer the long-term history of cyanobacteria and the recent rise of an anatoxin producing Dolichospermum sp.. Harmful Algae 101, 101971 (2021).CAS 
PubMed 
Article 

Google Scholar 
de la Escalera, G. M., Antoniades, D., Bonilla, S. & Piccini, C. Application of ancient DNA to the reconstruction of past microbial assemblages and for the detection of toxic cyanobacteria in subtropical freshwater ecosystems. Mol. Ecol. 23, 5791–5802 (2014).Article 
CAS 

Google Scholar 
Retrolens—Historical Imagery Resource. https://retrolens.co.nz/.Strayer, D. L. Alien species in fresh waters: ecological effects, interactions with other stressors, and prospects for the future. Freshw. Biol. 55, 152–174 (2010).Article 

Google Scholar 
Hall, S. R. & Mills, E. L. Exotic species in large lakes of the world. Aquat. Ecosyst. Health Manag. 3, 105–135 (2000).Article 

Google Scholar 
Gehrke, P. C. & Harris, J. H. The role of fish in cyanobacterial blooms in Australia. Mar. Freshw. Res. 45, 905–915 (1994).Article 

Google Scholar 
Burns, C. W. & Schallenberg, M. Impacts of nutrients and zooplankton on the microbial food web of an ultra-oligotrophic lake. J. Plankton Res. 20, 1501–1525 (1998).Article 

Google Scholar 
Rowe, D. K. & Schallenberg, M. Food webs in lakes. In Freshwaters of New Zealand (ed. Harding, J. S.) 23 (Wellington, N.Z.: New Zealand Hydrological Society, 2004).Gliwicz, Z. M. & Pijanowska, J. The role of predation in zooplankton succession. In Plankton Ecology: Succession in Plankton Communities (ed. Sommer, U.) 253–296 (Springer, 1989).Chapter 

Google Scholar 
Vanni, M. J. & Findlay, D. L. Trophic cascades and phytoplankton community structure. Ecology 71, 921–937 (1990).Article 

Google Scholar 
Smith, K. F. & Lester, P. J. Trophic interactions promote dominance by cyanobacteria (Anabaena spp.) in the pelagic zone of lower Karori reservoir, Wellington, New Zealand. N. Z. J. Mar. Freshw. Res. 41, 143–155 (2007).Article 

Google Scholar 
Smith, K. F. & Lester, P. J. Cyanobacterial blooms appear to be driven by top-down rather than bottom-up effects in the Lower Karori Reservoir (Wellington, New Zealand). N. Z. J. Mar. Freshw. Res. 40, 53–63 (2006).CAS 
Article 

Google Scholar 
Caroppo, C. Ecology and biodiversity of picoplanktonic cyanobacteria in coastal and brackish environments. Biodivers. Conserv. 24, 949–971 (2015).Article 

Google Scholar 
Pulina, S. et al. Picophytoplankton seasonal dynamics and interactions with environmental variables in three Mediterranean coastal lagoons. Estuaries Coasts 40, 469–478 (2017).CAS 
Article 

Google Scholar 
Callieri, C. Picophytoplankton in freshwater ecosystems: the importance of small-sized phototrophs. Freshw. Rev. 1, 1–28 (2008).Article 

Google Scholar 
Keefer, D. K. Investigating landslides caused by earthquakes: a historical review. Surv. Geophys. 23, 473–510 (2002).ADS 
Article 

Google Scholar 
Fan, X. et al. Earthquake-induced chains of geologic hazards: patterns, mechanisms, and impacts. Rev. Geophys. 57, 421–503 (2019).ADS 
Article 

Google Scholar 
Manighetti, I. et al. Repeated giant earthquakes on the Wairarapa fault, New Zealand, revealed by Lidar-based paleoseismology. Sci. Rep. 10, 2124 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
McSaveney, E. Historic earthquakes: the 1942 Wairarapa earthquakes. Te Ara Encyclopedia of New Zealand https://teara.govt.nz/en/historic-earthquakes/page-9 (2006).New Zealand’s environmental reporting series: our atmosphere and climate. (Ministry for the Environment & Stats NZ, 2020).Beng, K. C. & Corlett, R. T. Applications of environmental DNA (eDNA) in ecology and conservation: opportunities, challenges and prospects. Biodivers. Conserv. 29, 2089–2121 (2020).Article 

Google Scholar 
Freeland, J. R. The importance of molecular markers and primer design when characterizing biodiversity from environmental DNA. Genome https://doi.org/10.1139/gen-2016-0100 (2016).Article 
PubMed 

Google Scholar 
Barnes, M. A. et al. Environmental conditions influence eDNA persistence in aquatic systems. Environ. Sci. Technol. 48, 1819–1827 (2014).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Barnes, M. A. et al. Environmental conditions influence eDNA particle size distribution in aquatic systems. Environmental DNA https://doi.org/10.1002/edn3.160 (2020).Article 

Google Scholar 
Corinaldesi, C., Beolchini, F. & Dell’anno, A. Damage and degradation rates of extracellular DNA in marine sediments: implications for the preservation of gene sequences. Mol. Ecol. 17, 3939–3951 (2008).CAS 
PubMed 
Article 

Google Scholar 
Eichmiller, J. J., Best, S. E. & Sorensen, P. W. Effects of temperature and trophic state on degradation of environmental DNA in lake water. Environ. Sci. Technol. 50, 1859–1867 (2016).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Strickler, K. M., Fremier, A. K. & Goldberg, C. S. Quantifying effects of UV-B, temperature, and pH on eDNA degradation in aquatic microcosms. Biol. Conserv. 183, 85–92 (2015).Article 

Google Scholar 
Seymour, M. et al. Acidity promotes degradation of multi-species environmental DNA in lotic mesocosms. Commun. Biol. https://doi.org/10.1038/s42003-017-0005-3 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Dommain, R. et al. The challenges of reconstructing tropical biodiversity with sedimentary ancient DNA: a 2200-year-long metagenomic record from Bwindi Impenetrable Forest, Uganda. Front. Ecol. Evol. https://doi.org/10.3389/fevo.2020.00218 (2020).Article 

Google Scholar 
Jöhnk, K. D. et al. Summer heatwaves promote blooms of harmful cyanobacteria. Glob. Change Biol. 14, 495–512 (2008).ADS 
Article 

Google Scholar 
Sogin, M. L. et al. Microbial diversity in the deep sea and the underexplored “rare biosphere”. PNAS 103, 12115–12120 (2006).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

Our research involved work with animal subjects (unionid mussels and Round Goby fishes) and was conducted following relevant regulations and standard procedures. The field collections were carried out under Pennsylvania Fish and Boat Commission permits (# 2018-01-0136 and 2019-01-0026). The experimental protocols were approved by Penn State University’s Institutional Animal Care and Use Committee (IACUC# 201646941 and 201646962). All new DNA sequencing data are made publicly available in GenBank (with accession numbers provided in Table 1) and a BioProject (# PRJNA813547) of the National Center for Biotechnology Information40.Propensity of Round Goby to consume unionid mussels in a controlled lab settingStream table setupWe conducted lab experiments to observe the potential predation of juvenile freshwater mussels by the Round Goby, following standard research protocols for work with animal subjects (IACUC# 201646962, Penn State University). We constructed four artificial stream tables in an aquatic laboratory, each measuring 3 × 2 m and featuring two run and two pool sections (each 0.63 × 0.56 × 0.46 m). Water flow was produced using eight Homsay 920 GPH submersible water pumps, which pumped water from a central reservoir tub into each table at the start of each run section. The water flow direction was clockwise for stream tables 1 and 3, and counterclockwise for stream tables 2 and 4. Water pumped into the stream tables exited via two drains located medially of each run section, where it flowed back to the central reservoir tub. Each stream table was filled with a 6 mm layer of substrate consisting of a mixture of sand, gravel (4–6 mm), and crushed stone (size 2B, with an average size of ~ 19 mm). The day before each experiment, field technicians traveled to local streams and collected macroinvertebrates using one minute D-frame kick net samples for each of the four stream tables. The macroinvertebrates and associated substrate were transported back to the facility and were introduced into each stream table system.Preferential feeding experimentsBefore each experiment commenced, juvenile Plain Pocketbook mussels (Lampsillis cardium) were introduced into each stream table (with 165 mussel specimens for experiment 1 and 100 mussels for experiments 2 and 3). This  widespread and abundant species is not imperiled in Pennsylvania, and mussels were provided for this study by the White Sulphur Springs National Fish Hatchery located in southeast West Virgina. The mussels were allowed to acclimate in the stream tables for 2 h before commencing each experiment. Ten Round Gobies were introduced into each stream system (stream tables 1 and 2 for experiment 1, and all tables for experiments 2 and 3). The total length (from nose tip to caudal tip) of each fish was measured prior to introduction and after the termination of experiments 2 and 3. Experiment 1 was conducted for 3 weeks, while experiments 2 and 3 were conducted for 8 days. During these experiments, Round Gobies were allowed to exist in the systems and feed preferentially, on the mussels and macroinvertebrates, for the allotted time before each investigation concluded. We acknowledge that in these experiments, the mussel abundances are higher and macroinvertebrate densities lower and less rich than commonly occur in the natural stream environment. Further, the Round Goby fish densities used are much higher than currently in the French Creek watershed, though are comparable to what is currently seen in parts of the Great Lakes basin. Nonetheless, the experiment scenarios allowed us to observe if Round Gobies would consume the mussels when given the choice to feed on a variety of food items.Evaluation of unionids consumed by fishRound Gobies were removed from the stream tables upon completion of each experiment. They were euthanized using  > 250 mg/L buffered (pH ~ 7) tricaine-S (MS222) solution. The fish were submerged for 10 min beyond the cessation of opercular movement to ensure proper euthanasia, and tissues were collected after we confirmed complete euthanasia—compliant with AVMA guidelines and approved by the IACUC protocol. The Round Gobies were placed in a 10% solution of formalin for preservation, and after 2 weeks, they were rinsed with clean water and were placed in 70% ethanol for long-term storage. After fish were removed from the system, the water was drained, and the substrate was sifted to recover the remaining mussels. Mussels were counted, and live individuals were returned to holding tanks for use in subsequent experiments. To further assess whether Round Gobies had consumed mussels during the investigation, Round Gobies were x-rayed using a Bruker Skyscan 1176 micro-CT scanner. After that, the stomachs of each fish were excised, and the contents examined using a Leica CME dissection scope to confirm the identity of Plain Pocketbook mussels. Contents posterior to the stomach were not analyzed because they could not be reliably counted and identified.DNA metabarcoding to identify mussel species consumed by Round Goby in a stream settingFish and mussel sample acquisitionWe collected 39 Round Gobies directly from streams within the French Creek watershed—their newly invaded natural stream habitats—in June 2018. We aimed to quantify which species, if any, of unionid mussels they consumed. Fish collection locations included LeBoeuf Creek at Moore Road and 100 m below the confluence of French Creek and LeBoeuf Creek. The unionid mussel populations and the environmental field settings at these locations are detailed by Clark et al.19. A team of field technicians collected fish by kick seining (3 m × 1 m × 9.5 mm nylon mesh) while moving downstream. Seining was the sampling method of choice compared to electrofishing to avoid possible regurgitation of food items prior to excision of fishes’ stomachs. The stream reaches sampled at each location were between 100 and 200 m in length and included riffle, run, and pool habitats. In addition to fish samples, unionid mussel samples from French Creek were also collected for analysis (under Pennsylvania Fish & Boat Commission collectors permits # 2018-01-0136 and 2019-01-0026). Following standard research protocols (under IACUC# 201646941, Penn State University), the Round Gobies collected were euthanized using buffered Tricaine-S (MS222) solution; and stomachs were excised using sterilized utensils before being placed in sterilized tubes filled with 97% ethanol. After excision of stomachs, fishes were placed in a 10% formalin solution for preservation. After 2 weeks, fishes were rinsed with clean water and transferred to 70% ethanol for storage. The stomach samples were immediately placed in ethanol and on ice in the field. Samples were stored in a freezer before being shipped to the US Geological Survey’s Eastern Ecological Science Center for various molecular ecology analyses. Once the fish and mussel samples arrived at this lab, they were recorded and stored at four °C until analysis.Primer developmentSpecific primers targeting a moderately conserved region of the mitochondrial COI gene for 25 species of unionids inhabiting French Creek were designed. Previously a PCR-based amplification method utilizing restriction enzyme digests was used to identify genetic fingerprints of 25 unionid species inhabiting French Creek41. Here, we designed a new degenerate PCR primer set modified with sequencing overhangs to facilitate compatibility with a MiSeq amplicon sequencing method previously designed for 16S Amplicon sequencing. We targeted the locus of the mitochondrial COI gene of unionids known to inhabit the Atlantic Slope Drainage. Consensus sequences were derived using Multalin analysis and a tiling method to identify conserved primer binding regions flanking an ~ 300 bp region of the COI gene. This gene was targeted in part due to the availability of partial or complete sequences representing these target species in the NCBI reference database40. Cytochrome oxidase sequences were downloaded for the 25 unionid mussel species of interest. However, a COI sequence for the Rabbitsfoot (Theliderma cylindrica) mussel was absent from the NCBI database, which required us to sequence this region for an in-house reference (which is described later in the paper). We designed a degenerate primer cocktail specific to all mussel species of interest that amplified a ~ 289 bp product, with forward and reverse primers used for the amplification of unionid specific COI presented as supplemental information (see Table S-230. We evaluated the suitability of the primers using samples from field identified mussels. For primer optimization, PCR was run across a gradient of annealing temperatures to determine suitability. In addition, we used Round Goby DNA as a template to evaluate specificity. In addition to Round Goby stomach samples, mussel samples of several species collected from French Creek were included as positive controls.DNA extraction from tissue samplesFollowing the manufacturer’s protocols, tissue samples (including fish stomach and mussel tissue) were extracted with the Zymo Research ZymoBIOMICS 96 MagBead DNA Kit (San Diego, CA). Random samples of DNA extracts were analyzed on an Agilent 2100 Bioanalyzer using a high-sensitivity assay kit. Fragments in the target amplicon range were apparent (albeit not known to be of mussel origin). All samples were stored at − 20 °C until PCR was performed. DNA from both the T. cylindrica and L. complanata samples were analyzed for DNA quality.Rolling circle amplification of mitochondrial genomesTo acquire COI sequences for T. cylindrica and L. complanata, we subjected archived DNA samples to rolling circle amplification (RCA) followed by amplicon sequencing on the MiSeq. In short, 2 µl of DNA template was added to 2 µl Equiphi29 DNA polymerase reaction buffer containing 1 µl of Exonuclease-resistant random primers (ThermoFisher). Samples were denatured by heating to 95 °C for 3 min followed immediately by cooling on ice for more than 5 min. A volume of 5 µl was added to an RCA master mix containing 1.5 µl of 10 × Equiphi29 DNA polymerase reaction buffer, 0.2 µl of 100 mM dithiothreitol, 8 µl of 2.5 mM dNTPs, 1 µl of Eqiphi29 DNA polymerase (10U) and 4.3 µl of nuclease-free water. The samples were heated to 45 °C for 3 h and then 65 °C for 10 min. Samples were then placed in ice and then frozen at − 20 °C. All RCA products were normalized to 0.2 ng/µl in 10 mM Tris–HCl, pH 8.5. Normalized RCA product was utilized as a template for an Illumina Nextera XT library preparation. Sequencing libraries were prepared following the Nextera XT Library Preparation Reference Guide (CT# 15031942 v01) using the Nextera XT Library Preparation Kit (Illumina, San Diego, CA). Final libraries were analyzed for size and quality using the Agilent BioAnalyzer with the accompanying DNA 1000 Kit (Agilent, Santa Clara, CA). Libraries were quantified using the Qubit H.S. Assay Kit (Invitrogen, Carlsbad, CA) and normalized to 4 nM using 10 mM Tris, pH 8.5. Libraries were pooled and run on the Illumina MiSeq at a concentration of 10 pM with a 5% PhiX spike with run parameters of 1 × 150. Bioinformatic processing of this data is outlined below.Amplification of the cytochrome oxidase 1 geneExtracted genomic DNA was used as template for end-point PCR. Samples evaluated were from mussels and round gobies (see supporting Table S-330). The ~ 289 bp COI region was amplified with the mussel primers as follows. The amplification reaction contained 0.15 µM of each primer, 1 µL of the initial amplification product, and Promega Go Taq Green Master Mix following manufacturer recommendations for a 25 µL reaction. The thermocycler program consisted of an initial denaturing step of 95 °C for 3 min, followed by 30 cycles of 30 s at 95 °C, 30 s at 52 °C, and 1 min at 72 °C. Products were subjected to a final extension of 72 °C for 5 min then held until collection at 12 °C. An appropriately sized amplification product was confirmed for each reaction by electrophoresis of 5 µL of the reaction product through a 1.5% I.D. N.A. agarose gel (FMC Bioproducts) at 100 V for 45 min. PCR products were cleaned with the Qiagen Qiaquick PCR purification kit (Valencia, CA) and quantified using the Qubit dsDNA H.S. Assay Kit (Thermofisher Scientific, Grand Island, NY). Samples were diluted in 10 mM Tris buffer (pH 8.5) to a final concentration of 5 ng/µL.Generation of mock mussel samplesTo better understand and minimize sources of error or bias in the taxonomic assignment, we created a mock extraction by mixing sequences from known mussel taxa at defined concentrations. For each mussel, approximately 25-mg of tissue was extracted with the ZymoBIOMICS 96 MagBead DNA Kit (San Diego, CA) following the manufacturer’s protocol. The COI sequence was amplified from each species using the same primer-protocol combination described above. A total of 5 PCR products were mixed at equal concentration (mass/volume) to generate the mock sample (“Mock” hereafter). To confirm the identity of these inputs, each COI region was amplified and sequenced on the Illumina MiSeq during the same run as the Mock and samples.Sequencing library preparation and quality assessmentNext-generation sequencing was performed on the Illumina MiSeq platform to observe species-specific sequences and determine the diet of the Round Goby. Inclusion of the overhangs on the amplification primers allowed us to utilize the Illumina 16S Metagenomic Sequencing Library Preparation protocol42. Amplicon libraries were prepared following the same manufacturer’s protocol. All samples were indexed using the Illumina Nextera XT multiplex library indices. DNA read size spectra were determined with the Agilent 2100 Bioanalyzer using the Agilent DNA 1000 Kit (Santa Clara, Calif.). Libraries were quantified with the Qubit dsDNA H.S. Assay Kit (ThermoFisher Scientific, Grand Island, N.Y.) and normalized to 4 nM (nM) using 10 mM (mM) Tris (hydroxymethyl) aminomethane buffer pH 8.5. A final concentration of 10 picomolar library with a 6.5% PhiX control spike was created with the combined pool of all indexed libraries. All bioinformatic operations were completed on CLC Genomic Workbench v20 (Qiagen, Valencia, Calif.).Read filtering, trimming, and RNAseq metabarcoding assemblyFASTQ files from the sequencing runs were imported as paired-end reads into CLC Genomics Workbench v20.0.4 (Qiagen Bioinformatics, Redwood City, Calif.) for initial filtering of exogenous sequence adaptors and poor-quality base calls. The trimmed overlapping paired-end reads were mapped to the 25 target unionid sequences specific for the species of interest. Several mapping iterations were run using different levels of stringency. We utilized + 2/− 3 match-mismatch scoring and set the length fraction to 0.90. Analyses were iterated using different similarity fractions ranging from 0.90 to 0.99. Reads were annotated, and relative abundance was determined using a curated reference library (see supporting Datasets S-1 and S-230). More

Hughes, T. P. et al. Global warming and recurrent mass bleaching of corals. Nature 543, 373–377 (2017).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Glynn, P. W. Coral reef bleaching: facts, hypotheses and implications. Glob. Chang. Biol. 2, 495–509 (1996).ADS 
Article 

Google Scholar 
Riegl, B. & Piller, W. E. Possible refugia for reefs in times of environmental stress. Int. J. Earth Sci. 92, 520–531 (2003).Article 

Google Scholar 
Hinderstein, L. M. et al. Theme section on ‘Mesophotic Coral Ecosystems: Characterization, Ecology, and Management’. Coral Reefs 29, 247–251 (2010).ADS 
Article 

Google Scholar 
Bongaerts, P., Ridgway, T., Sampayo, E. M. & Hoegh-Guldberg, O. Assessing the ‘deep reef refugia’ hypothesis: Focus on Caribbean reefs. Coral Reefs 29, 309–327 (2010).Article 

Google Scholar 
Smith, T. B. et al. Caribbean mesophotic coral ecosystems are unlikely climate change refugia. Global Change Biol. 22, 2756–2765 (2016).ADS 
Article 

Google Scholar 
Frade, P. R. et al. Deep reefs of the Great Barrier Reef offer limited thermal refuge during mass coral bleaching. Nat. Commun. 9, 3447 (2018).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Holstein, D. M., Paris, C. B., Vaz, A. C. & Smith, T. B. Modeling vertical coral connectivity and mesophotic refugia. Coral Reefs 35, 23–37 (2016).ADS 
Article 

Google Scholar 
Prasetia, R., Sinniger, F., Hashizume, K. & Harii, S. Reproductive biology of the deep brooding coral Seriatopora hystrix: Implications for shallow reef recovery. PLoS ONE 12, e0177034 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Shlesinger, T., Grinblat, M., Rapuano, H., Amit, T. & Loya, Y. Can mesophotic reefs replenish shallow reefs? Reduced coral reproductive performance casts a doubt. Ecology 99, 421–437 (2018).PubMed 
Article 

Google Scholar 
Gleason, D. F. & Hofmann, D. K. Coral larvae: From gametes to recruits. J. Exp. Mar. Bio. Ecol. 408, 42–57 (2011).Article 

Google Scholar 
Hughes, T. P. & Tanner, J. E. Recruitment failure, life histories, and long-term decline of Caribbean corals. Ecology 81, 2250–2263 (2000).Article 

Google Scholar 
Bongaerts, P. et al. Deep reefs are not universal refuges: Reseeding potential varies among coral species. Sci. Adv. 3, e1602373 (2017).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
van Oppen, M. J. H., Bongaerts, P., Underwood, J. N., Peplow, L. M. & Cooper, T. F. The role of deep reefs in shallow reef recovery: An assessment of vertical connectivity in a brooding coral from west and east Australia. Mol. Ecol. 20, 1647–1660 (2011).PubMed 
Article 

Google Scholar 
Cohen, I. & Dubinsky, Z. Long term photoacclimation responses of the coral Stylophora pistillata to reciprocal deep to shallow transplantation: Photosynthesis and calcification. Front. Mar. Sci. 2, 45 (2015).Article 

Google Scholar 
Eyal, G. et al. Euphyllia paradivisa, a successful mesophotic coral in the northern Gulf of Eilat/Aqaba, Red Sea. Coral Reefs 35, 91–102 (2016).ADS 
Article 

Google Scholar 
Ben-Zvi, O. et al. Photophysiology of a mesophotic coral 3 years after transplantation to a shallow environment. Coral Reefs 39, 903–913 (2020).Article 

Google Scholar 
Murata, N., Takahashi, S., Nishiyama, Y. & Allakhverdiev, S. I. Photoinhibition of photosystem II under environmental stress. Biochim. Biophys. Acta Bioenerget. 1767, 414–421 (2007).CAS 
Article 

Google Scholar 
Takahashi, S. & Murata, N. How do environmental stresses accelerate photoinhibition?. Trends Plant Sci. 13, 178–182 (2008).CAS 
PubMed 
Article 

Google Scholar 
Cumbo, V. R., Baird, A. H. & van Oppen, M. J. H. The promiscuous larvae: Flexibility in the establishment of symbiosis in corals. Coral Reefs 32, 111–120 (2013).ADS 
Article 

Google Scholar 
Little, A. F., Van Oppen, M. J. H. & Willis, B. L. Flexibility in algal endosymbioses shapes growth in reef corals. Science 304, 1492–1494 (2004).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Sinniger, F., Morita, R. & Harii, S. ‘Locally extinct’ coral species Seriatopora hystrix found at upper mesophotic depths in Okinawa. Coral Reefs 32, 153 (2013).ADS 
Article 

Google Scholar 
Sinniger, F. et al. Overview of the mesophotic coral ecosystems around Sesoko Island, Okinawa, Japan. Galaxea J. Coral Reef Stud. 24, 69–76 (2022).Article 

Google Scholar 
Loya, Y. et al. Coral bleaching: the winners and the losers. Ecol. Lett. 4, 122–131 (2001).Article 

Google Scholar 
van Woesik, R., Sakai, K., Ganase, A. & Loya, Y. Revisiting the winners and the losers a decade after coral bleaching. Mar. Ecol. Prog. Ser. 434, 67–76 (2011).ADS 
Article 

Google Scholar 
Sinniger, F., Prasetia, R., Yorifuji, M., Bongaerts, P. & Harii, S. Seriatopora diversity preserved in upper mesophotic coral ecosystems in Southern Japan. Front. Mar. Sci. 4, 155 (2017).Article 

Google Scholar 
Atoda, K. The larva and postlarval development of some reef-building corals. V. Seriatopora hystrix. Sci. Rep. Tohoku Univ. 19, 33–39 (1951).
Google Scholar 
Hata, T. et al. Coral larvae are poor swimmers and require fine-scale reef structure to settle. Sci. Rep. 7, 2249 (2017).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Harii, S. & Kayanne, H. Larval dispersal, recruitment, and adult distribution of the brooding stony octocoral Heliopora coerulea on Ishigaki Island, southwest Japan. Coral Reefs 22, 188–196 (2003).Article 

Google Scholar 
Mulla, A. J., Lin, C. H., Takahashi, S. & Nozawa, Y. Photo-movement of coral larvae influences vertical positioning in the ocean. Coral Reefs 40, 1297–1306 (2021).Article 

Google Scholar 
Figueiredo, J., Baird, A. H., Harii, S. & Connolly, S. R. Increased local retention of reef coral larvae as a result of ocean warming. Nat. Clim. Chang. 4, 498–502 (2014).ADS 
Article 

Google Scholar 
Shanks, A. L., Largier, J., Brink, L., Brubaker, J. & Hooff, R. Demonstration of the onshore transport of larval invertebrates by the shoreward movement of an upwelling front. Limnol. Oceanogr. 45, 230–236 (2000).ADS 
Article 

Google Scholar 
Singh, T. et al. Long-term trends and seasonal variations in environmental conditions in Sesoko Island, Okinawa, Japan. Galaxea J. Coral Reef Stud. 24, 121–133 (2022).Article 

Google Scholar 
Roth, M. S., Fan, T.-Y. & Deheyn, D. D. Life history changes in coral fluorescence and the effects of light intensity on larval physiology and settlement in Seriatopora hystrix. PLoS ONE 8, e59476 (2013).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mundy, C. N. & Babcock, R. C. Role of light intensity and spectral quality in coral settlement: Implications for depth-dependent settlement?. J. Exp. Mar. Bio. Ecol. 223, 235–255 (1998).Article 

Google Scholar 
Nesa, B., Baird, A. H., Harii, S., Yakovleva, I. & Hidaka, M. Algal symbionts increase DNA damage in coral planulae exposed to sunlight. Zool. Stud. 51, 12–17 (2012).CAS 

Google Scholar 
Cunning, R. & Baker, A. C. Excess algal symbionts increase the susceptibility of reef corals to bleaching. Nat. Clim. Change 3, 259–262 (2013).ADS 
Article 

Google Scholar 
Nakamura, T. Mass coral bleaching event in Sekisei lagoon observed in the summer of 2016. J. Jpn. Coral Reef Soc. 19, 29–40 (2017).Article 

Google Scholar 
Sakai, K., Singh, T. & Iguchi, A. Bleaching and post-bleaching mortality of Acropora corals on a heat-susceptible reef in 2016. PeerJ 7, e8138 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Edmunds, P. J., Gates, R. D. & Gleason, D. F. The biology of larvae from the reef coral Porites astreoides, and their response to temperature disturbances. Mar. Biol. 139, 981–989 (2001).Article 

Google Scholar 
Baker, A. C. Reef corals bleach to survive change. Nature 411, 765–766 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Bongaerts, P. et al. Adaptive divergence in a scleractinian coral: Physiological adaptation of Seriatopora hystrix to shallow and deep reef habitats. BMC Evol. Biol. 11, 303 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Einbinder, S. et al. Novel adaptive photosynthetic characteristics of mesophotic symbiotic microalgae within the reef-building coral, Stylophora pistillata. Front. Mar. Sci. 3, 195 (2016).Article 

Google Scholar 
Rogers, C. S., Fitz, H. C., Gilnack, M., Beets, J. & Hardin, J. Scleractinian coral recruitment patterns at Salt River submarine canyon, St. Croix, U.S. Virgin Islands. Coral Reefs 3, 69–76 (1984).ADS 
Article 

Google Scholar 
Maida, M., Collb, J. C. & Sammarco, P. W. Shedding new light on scleractinian coral recruitment. J. Exp. Mar. Biol. Ecol. 180, 189–202 (1994).Article 

Google Scholar 
Sato, M. Mortality and growth of juvenile coral Pocillopora damicornis (Linnaeus). Coral Reefs 4, 27–33 (1985).ADS 
Article 

Google Scholar 
Nozawa, Y. Micro-crevice structure enhances coral spat survivorship. J. Exp. Mar. Biol. Ecol. 367, 127–130 (2008).Article 

Google Scholar 
Gleason, D. F. & Wellington, G. M. Ultraviolet radiation and coral bleaching. Nature 365, 836–838 (1993).ADS 
Article 

Google Scholar 
Shlesinger, T. & Loya, Y. Depth-dependent parental effects create invisible barriers to coral dispersal. Commun. Biol. 4, 1–10 (2021).Article 

Google Scholar 
Groves, S. H. et al. Growth rates of Porites astreoides and Orbicella franksi in mesophotic habitats surrounding St. Thomas, US Virgin Islands. Coral Reefs 37, 345–354 (2018).ADS 
Article 

Google Scholar 
Al-Horani, F. A., Al-Moghrabi, S. M. & De Beer, D. The mechanism of calcification and its relation to photosynthesis and respiration in the scleractinian coral Galaxea fascicularis. Mar. Biol. 142, 419–426 (2003).CAS 
Article 

Google Scholar 
Jiang, L. et al. Increased temperature mitigates the effects of ocean acidification on the calcification of juvenile Pocillopora damicornis, but at a cost. Coral Reefs 37, 71–79 (2018).ADS 
Article 

Google Scholar 
Jurriaans, S. & Hoogenboom, M. O. Thermal performance of scleractinian corals along a latitudinal gradient on the Great Barrier Reef. Philos. Trans. R. Soc. B Biol. Sci. 374, 20180546 (2019).CAS 
Article 

Google Scholar 
Brown, B. E. et al. Diurnal changes in photochemical efficiency and xanthophyll concentrations in shallow water reef corals: evidence for photoinhibition and photoprotection. Coral Reefs 18, 99–105 (1999).Article 

Google Scholar 
Salih, A., Larkum, A., Cox, G., Kühl, M. & Hoegh-Guldberg, O. Fluorescent pigments in corals are photoprotective. Nature 408, 850–853 (2000).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Matz, M. V., Marshall, N. J. & Vorobyev, M. Are corals colorful?. Photochem. Photobiol. 82, 345–350 (2006).CAS 
PubMed 
Article 

Google Scholar 
Haddock, S. H. D. & Dunn, C. W. Fluorescent proteins function as a prey attractant: Experimental evidence from the hydromedusa Olindias formosus and other marine organisms. Biol. Open 4, 1094–1104 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Eyal, G. et al. Spectral diversity and regulation of coral fluorescence in a mesophotic reef habitat in the Red Sea. PLoS ONE 10, 1–19 (2015).Article 
CAS 

Google Scholar 
Ben-Zvi, O., Eyal, G. & Loya, Y. Light-dependent fluorescence in the coral Galaxea fascicularis. Hydrobiologia 759, 15–26 (2015).Article 

Google Scholar 
Roth, M. et al. Fluorescent proteins in dominant mesophotic reef-building corals. Mar. Ecol. Prog. Ser. 521, 63–79 (2015).ADS 
CAS 
Article 

Google Scholar 
Ben-Zvi, O., Eyal, G. & Loya, Y. Response of fluorescence morphs of the mesophotic coral Euphyllia paradivisa to ultra-violet radiation. Sci. Rep. 9, 1–9 (2019).CAS 
Article 

Google Scholar 
Hughes, T. P. et al. Coral reefs in the Anthropocene. Nature 546, 82–90 (2017).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Oliver, E. C. J. et al. Longer and more frequent marine heatwaves over the past century. Nat. Commun. 9, 1324 (2018).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Nakamura, T., van Woesik, R. & Yamasaki, H. Photoinhibition of photosynthesis is reduced by water flow in the reef-building coral Acropora digitifera. Mar. Ecol. Prog. Ser. 301, 109–118 (2005).ADS 
Article 

Google Scholar  More