Ecology

Miyata, T., Miyazawa, S. & Yasunaga, T. Two types of amino acid substitutions in protein evolution. J. Mol. Evol. 12, 219–236 (1979).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Li, W.-H., Wu, C.-I. & Luo, C.-C. A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes. Mol. Biol. Evol. 2, 150–174 (1985).PubMed 

Google Scholar 
Bielawski, J. P. & Yang, Z. Positive and negative selection in the DAZ gene family. Mol. Biol. Evol. 18, 523–529 (2001).CAS 
PubMed 
Article 

Google Scholar 
Ohta, T. Slightly deleterious mutant substitutions in evolution. Nature 246, 96–98 (1973).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Ohta, T. The nearly neutral theory of molecular evolution. Annu. Rev. Ecol. Evol. Syst. 23, 263–286 (1992).Article 

Google Scholar 
Johnson, K. P. & Seger, J. Elevated rates of nonsynonymous substitution in island birds. Mol. Biol. Evol. 18, 874–881 (2001).CAS 
PubMed 
Article 

Google Scholar 
Woolfit, M. & Bromham, L. Population size and molecular evolution on islands. Proc. Biol. Sci. 272, 2277–2282 (2005).CAS 
PubMed 
PubMed Central 

Google Scholar 
Ross, L., Hardy, N. B., Okusu, A. & Normark, B. B. Large population size predicts the distribution of asexuality in scale insects. Evolution 67, 196–206 (2013).PubMed 
Article 

Google Scholar 
Weber, C. C., Nabholz, B., Romiguier, J. & Ellegren, H. Kr/Kc but not dN/dS correlates positively with body mass in birds, raising implications for inferring lineage-specific selection. Genome Biol. 15, 542 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Brandt, A. et al. Effective purifying selection in ancient asexual oribatid mites. Nat. Commun. 8, 873 (2017).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Figuet, E. et al. Life history traits, protein evolution, and the nearly neutral theory in amniotes. Mol. Biol. Evol. 33(6), 1517–1527 (2016).CAS 
PubMed 
Article 

Google Scholar 
Saclier, N. et al. Life history traits impact the nuclear rate of substitution but not the mitochondrial rate in isopods. Mol. Biol. Evol. 35, 2900–2912 (2018).CAS 
PubMed 
Article 

Google Scholar 
Hebert, P. D. The Daphnia of North America: An Illustrated Fauna (on CD-ROM) (CyberNatural Software, Guelph, 1995).
Google Scholar 
Colbourne, J. K. et al. Phylogenetics and evolution of a circumarctic species complex (Cladocera: Daphnia pulex). Biol. J. Linn. Soc. 65, 347–365 (1998).
Google Scholar 
Crease, T. J., Omilian, A. R., Costanzo, K. S. & Taylor, D. J. Transcontinental phylogeography of the Daphnia pulex species complex. PLoS ONE 7, e46620 (2012).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mergeay, J., Verschuren, D. & De Meester, L. Cryptic invasion and dispersal of an American Daphnia in East Africa. Limnol. Oceanogr. 50, 1278–1283 (2005).ADS 
CAS 
Article 

Google Scholar 
Ma, X. et al. Lineage diversity and reproductive modes of the Daphnia pulex group in Chinese lakes and reservoirs. Mol. Phylogenet. Evol. 130, 424–433 (2019).PubMed 
Article 

Google Scholar 
So, M. et al. Invasion and molecular evolution of Daphnia pulex in Japan. Limnol. Oceanogr. 60, 1129–1138 (2015).ADS 
Article 

Google Scholar 
Duggan, I. C. et al. Identifying invertebrate invasions using morphological and molecular analyses: North American Daphnia ‘pulex’ in New Zealand fresh waters. Aquat. Invasions 7, 585–590 (2012).Article 

Google Scholar 
Ye, Z. et al. The rapid, mass invasion of New Zealand by North American Daphnia “pulex”. Limnol. Oceanogr. 66, 2673–2683 (2021).ADS 
Article 

Google Scholar 
Paland, S., Colbourne, J. K. & Lynch, M. Evolutionary history of contagious asexuality in Daphnia pulex. Evolution 59, 800–813 (2005).CAS 
PubMed 
Article 

Google Scholar 
Muller, H. J. The relation of recombination to mutational advance. Mutat. Res. 106, 2–9 (1964).CAS 
PubMed 
Article 

Google Scholar 
Felsenstein, J. The evolutionary advantage of recombination. Genetics 78, 737–756 (1974).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Paland, S. & Lynch, M. Transitions to asexuality result in excess amino acid substitutions. Science 311, 990–992 (2006).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Johnson, S. G. & Howard, R. S. Contrasting patterns of synonymous and nonsynonymous sequence evolution in asexual and sexual freshwater snail lineages. Evolution 61, 2728–2735 (2007).PubMed 
Article 

Google Scholar 
Neiman, M. et al. Accelerated mutation accumulation in asexual lineages of a freshwater snail. Mol. Biol. Evol. 27, 954–963 (2010).CAS 
PubMed 
Article 

Google Scholar 
Henry, L., Schwander, T. & Crespi, B. J. Deleterious mutation accumulation in asexual Timema stick insects. Mol. Biol. Evol. 29, 401–408 (2012).CAS 
PubMed 
Article 

Google Scholar 
Tucker, A. E. et al. Population-genomic insights into the evolutionary origin and fate of obligately asexual Daphnia pulex. Proc. Natl. Acad. Sci. 110, 15740–15745 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Colbourne, J. K. et al. The ecoresponsive genome of Daphnia pulex. Science 331, 555–561 (2011).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ye, Z. et al. A new reference genome assembly for the microcrustacean Daphnia pulex. G3 (Bethesda) 7, 1405–1416 (2017).CAS 
Article 

Google Scholar 
Keith, N. et al. High mutational rates of large-scale duplication and deletion in Daphnia pulex. Genome Res. 26, 60–69 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hall, D. J. An experimental approach to the dynamics of a natural population of Daphnia galeata mendotae. Ecology 45, 94–112 (1964).Article 

Google Scholar 
McCauley, E., Murdoch, W. W. & Nisbet, R. M. Growth, reproduction, and mortality of Daphnia pulex Leydig: Life at low food. Ecology 4, 505–514 (1990).
Google Scholar 
Xu, S. et al. High mutation rates in the mitochondrial genomes of Daphnia pulex. Mol. Biol. Evol. 29, 763–769 (2012).CAS 
PubMed 
Article 

Google Scholar 
Zheng, Y., Peng, R., Kuro-o, M. & Zeng, X. Exploring patterns and extent of bias in estimating divergence time from mitochondrial DNA sequence data in a particular lineage: A case study of salamanders (Order Caudata). Mol. Biol. Evol. 28, 2521–2535 (2011).CAS 
PubMed 
Article 

Google Scholar 
Zaret, T. M. Predation and Freshwater Communities (Yale University Press, New Haven, 1980).
Google Scholar 
Lynch, M. Predation, competition, and zooplankton community structure: An experimental study. Limnol. Oceanogr. 24, 253–272 (1979).ADS 
Article 

Google Scholar 
Mills, E. L. & Forney, J. L. Impact on Daphnia pulex of predation by young yellow perch in Oneida Lake, New York. Trans. Am. Fish. Soc. 112(2A), 154–161 (1983).Article 

Google Scholar 
Craddock, D. R. Effects of increased water temperature on Daphnia pulex. Fish. Bull. 74, 403–408 (1976).
Google Scholar 
Maruoka, N. & Urabe, J. Inter and intraspecific competitive abilities and the distribution ranges of two Daphnia species in Eurasian continental islands. Popul. Ecol. 62, 353–363 (2020).Article 

Google Scholar 
Dodson, S. I. & Hanazato, T. Commentary on effects of anthropogenic and natural organic chemicals on development, swimming behavior, and reproduction of Daphnia, a key member of aquatic ecosystems. Environ. Health Perspect. 103(Suppl 4), 7–11 (1995).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Claska, M. E. & Gilbert, J. J. The effect of temperature on the response of Daphnia to toxic cyanobacteria. Freshw. Biol. 39, 221–232 (1998).Article 

Google Scholar 
Bast, J. et al. Consequences of asexuality in natural populations: Insights from stick insects. Mol. Biol. Evol. 35, 1668–1677 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hartfield, M. Evolutionary genetic consequences of facultative sex and outcrossing. J Evol Biol 29, 5–22 (2016).CAS 
PubMed 
Article 

Google Scholar 
Hörandl, E. et al. Genome evolution of asexual organisms and the paradox of sex in eukaryotes. In Evolutionary Biology—A Transdisciplinary Approach (ed. Pontarotti, P.) (Springer, Cham, 2020). https://doi.org/10.1007/978-3-030-57246-4_7.Chapter 

Google Scholar 
Lynch, M., Bürger, R., Butcher, D. & Gabriel, W. The mutational meltdown in asexual populations. J. Hered. 84, 339–344 (1993).CAS 
PubMed 
Article 

Google Scholar 
Gordo, I. & Charlesworth, B. The degeneration of asexual haploid populations and the speed of Muller’s ratchet. Genetics 154, 1379–1387 (2000).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Downing, J. A. et al. The global abundance and size distribution of lakes, ponds, and impoundments. Limnol. Oceanogr. 51, 2388–2397 (2006).ADS 
Article 

Google Scholar 
McDonald, C. P., Rover, J. A., Stets, E. G. & Striegl, R. G. The regional abundance and size distribution of lakes and reservoirs in the United States and implications for estimates of global lake extent. Limnol. Oceanogr. 57, 597–606 (2012).ADS 
Article 

Google Scholar 
De Meester, L., Góme, A., Okamura, B. & Schwenk, K. The monopolization hypothesis and the dispersal-gene flow paradox in aquatic organisms. Acta Oecol. 23, 121–135 (2002).ADS 
Article 

Google Scholar 
Fukami, T., Bezemer, T. M., Mortimer, S. R. & Van Der Putten, W. H. Species divergence and trait convergence in experimental plant community assembly. Ecol. Lett. 8, 1283–1290 (2005).Article 

Google Scholar 
Makino, T. & Kawata, M. Invasive invertebrates associated with highly duplicated gene content. Mol. Ecol. 28, 1652–1663 (2019).CAS 
PubMed 
Article 

Google Scholar 
Kondrashov, F. A. Gene duplication as a mechanism of genomic adaptation to a changing environment. Proc. R. Soc. Lond. B Biol. Sci. 279, 5048–5057 (2012).
Google Scholar 
Barrick, J. E. & Lenski, R. E. Genome dynamics during experimental evolution. Nat. Rev. Genet. 14, 827–839 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rocha, E. P. C. Neutral theory, microbial practice: Challenges in bacterial population genetics. Mol. Biol. Evol. 35, 1338–1347 (2018).CAS 
PubMed 
Article 

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

Google Scholar 
Tanabe, A. S. Kakusan4 and Aminosan: Two programs for comparing nonpartitioned, proportional and separate models for combined molecular phylogenetic analyses of multilocus sequence data. Mol. Ecol. Resour. 11, 914–921 (2011).PubMed 
Article 

Google Scholar 
Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tian, X., Ohtsuki, H. & Urabe, J. Evolution of asexual Daphnia pulex in Japan: Variations and covariations of the digestive, morphological and life history traits. BMC Evol. Biol. 19, 122 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Chen, Y. et al. SOAPnuke: A MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data. Gigascience 7, 1–6 (2018).ADS 
PubMed 
PubMed Central 

Google Scholar 
Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2019).Article 
CAS 

Google Scholar 
Lee, T. H. et al. SNPhylo: A pipeline to construct a phylogenetic tree from huge SNP data. BMC Genomics 15, 162 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
R Core Team, R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna, Austria, 2019). https://www.R-project.org/ More

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Google Scholar  More

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

Wäckers, F. L. & van Rijn, P. C. J. Food for Protection: An Introduction. In Plant-Provided Food for Carnivorous Insects: A Protective Mutualism and its Applications (eds Wäckers, F. L. et al.) 1–14 (Cambridge University Press, 2005). https://doi.org/10.1017/CBO9780511542220.002.Chapter 

Google Scholar 
Benelli, G. et al. The impact of adult diet on parasitoid reproductive performance. J. Pest Sci. 90, 807–823. https://doi.org/10.1007/s10340-017-0835-2 (2017).Article 

Google Scholar 
Wäckers, F. Assessing the suitability of flowering herbs as parasitoid food sources: Flower attractiveness and nectar accessibility. Biol. Control. 29, 307–314. https://doi.org/10.1016/j.biocontrol.2003.08.005 (2004).Article 

Google Scholar 
Heimpel, G. E. & Jervis, M. A. Does Floral Nectar Improve Biological Control by Parasitoids? In Plant-Provided Food for Carnivorous Insects: A Protective Mutualism and its Applications (eds Wäckers, F. L. et al.) 267–304 (Cambridge University Press, 2009). https://doi.org/10.1017/CBO9780511542220.010.Chapter 

Google Scholar 
Wäckers, F. L. Suitability of (extra-)Floral Nectar, Pollen, and Honeydew as Insect Food Sources. In Plant-Provided Food for Carnivorous Insects: A Protective Mutualism and its Applications (eds Wäckers, F. L. et al.) 17–74 (Cambridge University Press, 2005). https://doi.org/10.1017/CBO9780511542220.003.Chapter 

Google Scholar 
Wäckers, F. L., van Rijn, P. C. & Heimpel, G. E. Honeydew as a food source for natural enemies: Making the best of a bad meal?. Biol. Control. 45, 176–184. https://doi.org/10.1016/j.biocontrol.2008.01.007 (2008).Article 

Google Scholar 
Jervis, M. A., Ellers, J. & Harvey, J. A. Resource acquisition, allocation, and utilization in parasitoid reproductive strategies. Annu. Rev. Entomol. 53, 361–385. https://doi.org/10.1146/annurev.ento.53.103106.093433 (2008).CAS 
Article 

Google Scholar 
Rosenheim, J. A. An evolutionary argument for egg limitation. Evolution 50, 2089–2094 (1996).Article 

Google Scholar 
Rosenheim, J. A. The relative contributions of time and eggs to the cost of reproduction. Evolution 53, 376–385 (1999).Article 

Google Scholar 
Rosenheim, J. A., Jepsen, S. J., Matthews, C. E., Smith, D. S. & Rosenheim, M. R. Time limitation, egg limitation, the cost of oviposition, and lifetime reproduction by an insect in nature. Am. Nat. 172, 486–496 (2008).Article 

Google Scholar 
Rosenheim, J. A., Heimpel, G. E. & Mangel, M. Egg maturation, egg resorption and the costliness of transient egg limitation in insects. Proc. Royal Soc London. Ser. B Biol. Sci. 267, 1565–1573 (2000).CAS 
Article 

Google Scholar 
Takasu, K. & Hirose, Y. Host searching behavior in the parasitoid Ooencyrtus nezarae Ishii (Hymenoptera: Encyrtidae) as influenced by non-host food deprivation. Appl. Entomol. Zool. 26, 415–417. https://doi.org/10.1303/aez.26.415 (1991).Article 

Google Scholar 
Sisterson, M. S. & Averill, A. L. Costs and benefits of food foraging for a braconid parasitoid. J. Insect Behav. 15, 571–588. https://doi.org/10.1023/A:1016389402543 (2002).Article 

Google Scholar 
Jacob, H. S. & Evans, E. W. Influence of food deprivation on foraging decisions of the parasitoid Bathyplectes curculionis (Hymenoptera: Ichneumonidae). Ann. Entomol. Soc. Am. 94, 605–611. https://doi.org/10.1603/0013-8746(2001)094[0605:iofdof]2.0.co;2 (2001).Article 

Google Scholar 
Siekmann, G., Keller, M. A. & Tenhumberg, B. The sweet tooth of adult parasitoid cotesia rubecula: Ignoring hosts for nectar?. J. Insect Behav. 17, 459–476. https://doi.org/10.1023/b:joir.0000042535.76279.c7 (2004).Article 

Google Scholar 
Williams, L., Deschodt, P., Pointurier, O. & Wyckhuys, K. A. Sugar concentration and timing of feeding affect feeding characteristics and survival of a parasitic wasp. J. Insect Physiol. 79, 10–18. https://doi.org/10.1016/j.jinsphys.2015.05.004 (2015).CAS 
Article 

Google Scholar 
Talamas, E. J. et al. A maximalist approach to the systematics of a biological control agent: Gryon aetherium Talamas, sp. nov. (Hymenoptera, Scelionidae). J. Hymenopt. Res. 87, 323–480. https://doi.org/10.3897/jhr.87.72842 (2021).Article 

Google Scholar 
Straser, R. K., Daane, K. M., Talamas, E. & Wilson, H. Evaluation of egg parasitoid Hadronotus pennsylvanicus as a prospective biocontrol agent of the leaffooted bug Leptoglossus zonatus. Biocontrol https://doi.org/10.1007/s10526-022-10131-z (2022).Article 

Google Scholar 
Mitchell, P. L. & Mitchell, F. L. Parasitism and predation of leaffooted bug (Hemiptera: Heteroptera: Coreidae) eggs. Ann. Entomol. Soc. Am. 79, 854–860. https://doi.org/10.1093/aesa/79.6.854 (1986).Article 

Google Scholar 
Yasuda, K. Function of the male pheromone of the leaf-footed plant bug, Leptoglossus australis (Fabricius) (Heteroptera: Coreidae) and its kairomonal effect. Jpn. Agric. Res. Q. 32, 161 (1998).CAS 

Google Scholar 
Bates, S. L. & Borden, J. H. Parasitoids of Leptoglossus occidentalis Heidemann (Heteroptera: Coreidae) in British Columbia. J. Entomol. Soc. Br. Columbia 101, 143–144 (2004).
Google Scholar 
Maltese, M., Caleca, V., Guerrieri, E. & Strong, W. B. Parasitoids of Leptoglossus occidentalis Heidemann (Heteroptera: Coreidae) recovered in western North America and first record of its egg parasitoid Gryon pennsylvanicum (Ashmead) (Hymenoptera: Platygastridae) in California. The Pan-Pacific Entomol. 88, 347–355. https://doi.org/10.3956/2012-23.1 (2012).Article 

Google Scholar 
Roversi, P. F. et al. Pre-release risk assessment of the egg-parasitoid Gryon pennsylvanicum for classical biological control of Leptoglossus occidentalis. J. Appl. Entomol. 138, 27–35. https://doi.org/10.1111/jen.12062 (2013).Article 

Google Scholar 
Nechols, J. R., Tracy, J. L. & Vogt, E. A. Comparative ecological studies of indigenous egg parasitoids (Hymenoptera: Scelionidae: Encyrtidae) of the squash bug, Anasa tristis (Hemiptera: Coreidae). J. Kansas Entomol. Soc. 62, 177–188 (1989).
Google Scholar 
Cornelius, M. L., Buffington, M. L., Talamas, E. J. & Gates, M. W. Impact of the egg parasitoid, Gryon pennsylvanicum (Hymenoptera: Scelionidae), on sentinel and wild egg masses of the squash bug (Hemiptera: Coreidae) in Maryland. Environ. Entomol. 45, 367–375. https://doi.org/10.1093/ee/nvv228 (2016).Article 

Google Scholar 
Cornelius, M. L., Hu, J. S. & Vinyard, B. T. Comparative study of egg parasitism by Gryon pennsylvanicum (Hymenoptera: Scelionidae) on two squash bug species Anasa tristis and Anasa armigera (Hemiptera: Coreidae). Environ. Entomol. https://doi.org/10.1093/ee/nvy145 (2018).Article 

Google Scholar 
Daane, K. M. et al. Stink bugs and leaffooted bugs. Pistachio Prod. Man. Publ. 3545, 225–238 (2016).
Google Scholar 
Joyce, A. L., Higbee, B. S., Haviland, D. R. & Brailovsky, H. Genetic variability of two leaffooted bugs, Leptoglossus clypealis and Leptoglossus zonatus (Hemiptera: Coreidae) in the Central Valley of California. J. Econ. Entomol. 110, 2576–2589. https://doi.org/10.1093/jee/tox222 (2017).CAS 
Article 

Google Scholar 
Zalom, F. G., Haviland, D. R., Symmes, E. T. & Tollerup, K. Almonds: Insects and Mites. University of California, Agriculture and Natural Resources, Oakland, CA, USA, University of California IPM Pest Management Guidelines, Publication 3431 ed. (2018).Michailides, T. J., Rice, R. E. & Ogawa, J. M. Succession and significance of several hemipterans attacking a pistachio orchard. J. Econ. Entomol. 80, 398–406. https://doi.org/10.1093/jee/80.2.398 (1987).Article 

Google Scholar 
Michailides, T. The ‘Achilles heel’of pistachio fruit. Calif. Agric. 43, 10–11 (1989).
Google Scholar 
Michailides, T. J. & Morgan, D. P. Association of botryosphaeria panicle and shoot blight of pistachio with injuries of fruit caused by hemiptera insects and birds. Plant Dis. 100, 1405–1413. https://doi.org/10.1094/pdis-09-15-1077-re (2016).Article 

Google Scholar 
Daane, K. et al. Large bugs damage pistachio nuts most severely during midseason. Calif. Agric. 59, 95–102 (2005).Article 

Google Scholar 
Haviland, D., Bentley, W., Beede, R. & Daane, K. Pistachios: Insects and mites. Univ. California IPM Pest Manag. Guidel. Publ. 3461 (2018).Joyce, A. L., Barman, A. K., Doll, D. & Higbee, B. S. Assessing feeding damage from two leaffooted bugs, Leptoglossus clypealis Heidemann and Leptoglossus zonatus (Dallas) (Hemiptera: Coreidae), on four almond varieties. Insects 10, 333. https://doi.org/10.3390/insects10100333 (2019).Article 

Google Scholar 
Stahl, J. M., Scaccini, D., Pozzebon, A. & Daane, K. M. Comparing the feeding damage of the invasive brown marmorated stink bug to a native stink bug and leaffooted bug on California pistachios. Insects 11, 688. https://doi.org/10.3390/insects11100688 (2020).Article 

Google Scholar 
Olson, D. L. & Nechols, J. R. Effects of squash leaf trichome exudates and honey on adult feeding, survival, and fecundity of the squash bug (Heteroptera: Coreidae) egg parasitoid Gryon pennsylvanicum (Hymenoptera: Scelionidae). Environ. Entomol. 24, 454–458. https://doi.org/10.1093/ee/24.2.454 (1995).Article 

Google Scholar 
Sabbatini Peverieri, G., Furlan, P., Simoni, S., Strong, W. & Roversi, P. Laboratory evaluation of Gryon pennsylvanicum (Ashmead) (Hymenoptera: Platygastridae) as a biological control agent of Leptoglossus occidentalis Heidemann (Heteroptera: Coreidae). Biol. Control. 61, 104–111. https://doi.org/10.1016/j.biocontrol.2012.01.005 (2012).Article 

Google Scholar 
Cornelius, M. L., Vinyard, B. T., Mowery, J. D. & Hu, J. S. Ovipositional behavior of the egg parasitoid Gryon pennsylvanicum (Hymenoptera: Scelionidae) on two squash bug species Anasa tristis (Hemiptera: Coreidae) and Anasa armigera: Effects of parasitoid density, nutrition, and host egg chorion on parasitism rates. Environ. Entomol. 49, 1307–1315. https://doi.org/10.1093/ee/nvaa118 (2020).CAS 
Article 

Google Scholar 
Vogt, E. & Nechols, J. The influence of host deprivation and host source on the reproductive biology and longevity of the squash bug egg parasitoid Gryon pennsylvanicum (Ashmead) (Hymenoptera: Scelionidae). Biol. Control. 3, 148–154. https://doi.org/10.1006/bcon.1993.1022 (1993).Article 

Google Scholar 
Olson, D., Fadamiro, H., Lundgren, J. & Heimpel, G. E. Effects of sugar feeding on carbohydrate and lipid metabolism in a parasitoid wasp. Physiol. Entomol. 25, 17–26 (2000).CAS 
Article 

Google Scholar 
Jervis, M. A., Heimpel, G. E., Ferns, P. N., Harvey, J. A. & Kidd, N. A. C. Life-history strategies in parasitoid wasps: A comparative analysis of “ovigeny”. J. Animal Ecol. 70, 442–458. https://doi.org/10.1046/j.1365-2656.2001.00507.x (2001).Article 

Google Scholar 
Jervis, M. A. & Ferns, P. N. The timing of egg maturation in insects: Ovigeny index and initial egg load as measures of fitness and of resource allocation. Oikos 107, 449–461 (2004).Article 

Google Scholar 
Lee, J. C. & Heimpel, G. E. Effect of floral nectar, water, and feeding frequency on cotesia glomerata longevity. Biocontrol 53, 289–294 (2008).Article 

Google Scholar 
Wu, H., Meng, L. & Li, B. Effects of feeding frequency and sugar concentrations on lifetime reproductive success of Meteorus pulchricornis (Hymenoptera: Braconidae). Biol. Control. 45, 353–359. https://doi.org/10.1016/j.biocontrol.2008.01.017 (2008).CAS 
Article 

Google Scholar 
King, B. H. Offspring sex ratios in parasitoid wasps. Q. Rev. Biol. 62, 367–396. https://doi.org/10.1086/415618 (1987).Article 

Google Scholar 
Berndt, L. A. & Wratten, S. D. Effects of alyssum flowers on the longevity, fecundity, and sex ratio of the leafroller parasitoid Dolichogenidea tasmanica. Biol. Control. 32, 65–69. https://doi.org/10.1016/j.biocontrol.2004.07.014 (2005).Article 

Google Scholar 
Sabbatini Peverieri, G. et al. Host egg age of Leptoglossus occidentalis (Heteroptera: Coreidae) and parasitism by Gryon pennsylvanicum (Hymenoptera: Platygastridae). J. Econ. Entomol. 106, 633–640. https://doi.org/10.1603/ec12344 (2013).Article 

Google Scholar 
Abram, P. K., Brodeur, J., Urbaneja, A. & Tena, A. Nonreproductive effects of insect parasitoids on their hosts. Annu. Rev. Entomol. 64(1), 259–276 (2019).CAS 
Article 

Google Scholar 
Lewis, W. & Takasu, K. Use of learned odours by a parasitic wasp in accordance with host and food needs. Nature 348, 635–636 (1990).ADS 
Article 

Google Scholar 
Takasu, K. & Lewis, W. Importance of adult food sources to host searching of the larval parasitoid Microplitis croceipes. Biol. Control 5, 25–30 (1995).Article 

Google Scholar 
Wäckers, F. The effect of food deprivation on the innate visual and olfactory preferences in the parasitoid Cotesia rubecula. J. Insect Physiol. 40, 641–649 (1994).Article 

Google Scholar 
Lightle, D., Ambrosino, M. & Lee, J. C. Sugar in moderation: Sugar diets affect short-term parasitoid behaviour. Physiol. Entomol. 35, 179–185 (2010).CAS 
Article 

Google Scholar 
Varennes, Y.-D., Gonzalez Chang, M., Boyer, S. & Wratten, S. Nectar feeding increases exploratory behaviour in the aphid parasitoid Diaeretiella rapae (Mcintosh). J. Appl. Entomol. 140, 479–483 (2016).Article 

Google Scholar 
Takano, S. & Takasu, K. Food deprivation increases reproductive effort in a parasitoid wasp. Biol. Control. 133, 75–80. https://doi.org/10.1016/j.biocontrol.2019.03.010 (2019).Article 

Google Scholar 
Landis, D. A., Wratten, S. D. & Gurr, G. M. Habitat management to conserve natural enemies of arthropod pests in agriculture. Annu. Rev. Entomol. 45(1), 175–201 (2000).CAS 
Article 

Google Scholar 
Masner, L. A revision of gryon haliday in North America (Hymenoptera: Proctotrupoidea: Scelionidae). Can. Entomol. 115, 123–174. https://doi.org/10.4039/ent115123-2 (1983).Article 

Google Scholar 
Vogt, E. A. & Nechols, J. R. Diel activity patterns of the squash bug egg parasitoid Gryon pennsylvanicum (Hymenoptera: Scelionidae). Ann. Entomol. Soc. Am. 84, 303–308. https://doi.org/10.1093/aesa/84.3.303 (1991).Article 

Google Scholar 
Wiedemann, L. M., Canto-Silva, C. R., Romanowski, H. P. & Redaelli, L. R. Oviposition behavior of Gryon gallardoi (Hym.: Scelionidae) on eggs of Spartocera dentiventris (Hem.: Coreidae). Braz. J. Biol. 63, 133 (2003).CAS 
Article 

Google Scholar 
Friard, O. & Gamba, M. BORIS: a free, versatile open-source event-logging software for video/audio coding and live observations. Methods Ecol. Evol. 7, 1325–1330. https://doi.org/10.1111/2041-210x.12584 (2016).Article 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. https://www.r-project.org/ (2019).Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48. https://doi.org/10.18637/jss.v067.i01 (2015).Article 

Google Scholar 
Hothorn, T., Bretz, F. & Westfall, P. Simultaneous inference in general parametric models. Biom. J. 50, 346–363 (2008).MathSciNet 
Article 

Google Scholar  More

Study siteThe study area covers the campus of Suranaree University of Technology (SUT) and its surrounding landscape in Muang, Nakhon Ratchasima, Thailand (14.879° N, 102.018° E; Fig. 1). The university campus covers about 11.2 km2, and comprises a matrix of human modified lands interspersed with mixed deciduous forest fragments (at the onset of this study we identified there were 37 mixed deciduous forest fragments on campus, mean = 7.36 ± 1.48 ha, range = 0.45–45.6 ha [note, “±” is used for standard error throughout the text]). More than 15,000 students are enrolled at SUT, and there are numerous multi-story classrooms, laboratory and workshop buildings, residential housing, parking areas, eating and sports facilities, an elementary school, and a large hospital on the university campus. During the first term of the 2019 school year, 7622 students, as well as numerous SUT staff, lived in on-campus residential areas. The landscape surrounding the university is primarily dominated by agriculture, though there are also patches of less-disturbed areas as well as several densely populated villages and suburban housing divisions among the monoculture plots of upland crops (e.g., cassava, maize, and eucalyptus).Figure 1Study site map illustrating the land-use types spanning the area where the Malayan kraits (Bungarus candidus) were tracked in Muang Nakhon Ratchasima, Nakhon Ratchasima province, Thailand. Map created using QGIS v.3.8.2 (https://qgis.org/) in combination with Inkscape v.1.1.0 (https://inkscape.org/).Full size imageThe study site is located within the Korat Plateau region with an altitude range of 205–285 m above sea level. Northeast Thailand has a tropical climate, and the average daily temperature from 1 January 2018 to 31 December 2020 in Muang Nakhon Ratchasima was 28.29 °C, with daily averages ranging from 19.3 to 34.1 °C38. The region receives an average annual rainfall ranging from 1270 to 2000 mm39. There are three distinct seasons in northeast Thailand: cold, wet, and hot, each are classified by annual changes in temperature and rainfall. Cold season is typically between mid-October and mid-February, hot season is generally from mid-February to May, while the highly unpredictable rainfall of the wet season is predominantly concentrated between the months May to October39,40.Due to the representation of agriculture, semi-urban, and suburban areas with patches of more natural areas all within a relatively small area, we determined the university campus provided an ideal setting to examine how land-use features and human activity influence the movements of B. candidus. Additionally, past studies have indicated northeast Thailand hosts the most bites by B. candidus in Thailand29,33, making sites like ours ideal.Study animalsWe opportunistically sampled Malayan kraits captured as a result of notifications from locals and ad-hoc encounters during transit due to low detectability in visual encounter surveys, in addition to those discovered through unstandardized visual encounter surveys. Upon capture, we collected morphometric data, including snout-vent length (SVL), tail length (TL), mass, and sex (Table 1, Supp. Table 1). We measured body lengths with a tape measure, measured body mass with a digital scale, and determined sex via cloacal probing, all while the snakes were anesthetized via inhaling vaporized isoflurane. We then housed individuals with an SVL > 645 mm and mass > 50 g in plastic boxes (with refugia and water) prior to surgical transmitter implantation by a veterinarian from the Nakhon Ratchasima Zoo. We attempted to minimize the time snakes were in captivity awaiting implantation; however, delays arose due to the veterinarian’s availability, the snake being mid-ecdysis, or the snake having a bolus that needed to pass through the digestive tract before implantation (n = 21 implantations, mean = 5.02 ± 0.61 days, range = 0.60–13.02 days). The Nakhon Ratchasima Zoo veterinarian implanted radio transmitters (1.8 g BD-2 or 3.6 g SB-2 Holohil Inc, Carp, Canada) into the coelomic cavity using procedures described by Reinert and Cundall41, while the snake was anesthetized. We assigned each individual an ID according to sex and individual detection number (e.g., M02 = a male was the second B. candidus individual documented during the study). We released the implanted individuals as close as possible to their capture locations (mean = 65.31 m ± 13.7 m, range = 0–226.42 m), though on six occasions we moved individuals ≥ 100 m because the individual came from either residential areas or a busy road (all but one were moved  800 mm; thus, nine of the males were adults and four were juveniles (though two of the males had an SVL > 720 mm, and therefore likely sub-adults). The single telemetered female was an adult.Individual tracking durations varied (mean = 106.46 ± 15.36 days, range = 28.5–222.77 days; Supp. Fig. 1), as many individuals were lost due to unexpected premature transmitter failures (n = 5) or unsuccessful recapture efforts due to individuals sheltering under large buildings as the transmitter reached the end of its battery life (n = 4). We only recorded one confirmed mortality in the study, M01, who was killed by a motorized vehicle when crossing a road (n = 1). Another three individuals were lost due to unknown reasons, which may have been due to premature transmitter failure, mortality, or the animal moving beyond radio signal despite extensive search efforts. Thus, we only successfully recaptured and re-implanted five individuals (M01 once, M02 twice, M07 once, M27 once, and M33 twice). Transmitter batteries generally lasted approximately 90–110 days, so we aimed to replace transmitters after ≥ 90 days of use. At the end of the study, only one individual was successfully recaptured to remove the transmitter.Data collectionWe used very high frequency radio-telemetry to locate each telemetered individual on average every 24.20 h (SE ± 0.41, 0.17–410.0 h; see Supp. Fig. 2 for distribution of tracking time lags). We aimed to locate each individual’s shelter locations once each day during the daylight (06:00–18:00 h); however, we were occasionally (n = 34 days) unable to locate a snake for several consecutive days when we were unable to obtain radio signal due to an individual having moved far away or deep underneath a large structure. There were also a few occasions where we were unable to track snakes due to prolonged and heavy rainfall (n = 4 days), as the moisture damages equipment, or other reasons (n = 4 days). We additionally located snakes nocturnally (18:00–06:00 h) ad hoc and in an attempt to observe nocturnal behaviors and movement pathways when animals were active. We defined fixes as any time a telemetered individual was located, and relocations (i.e., moves) as the occasions where we located an individual > 5 m from its previous known location.Each day we manually honed in on signal via a radio receiver to locate individuals (as described by Amelon et al.42, and recorded locations in Universal Transverse Mercator (UTM; 47 N World Geodetic System 84) coordinate reference system with a handheld global positioning system (GPS) unit (Garmin 64S GPS, Garmin International, Inc., Olathe, Kansas) directly above the sheltered snake. We generally approached within one meter of sheltering snakes during daylight to precisely record shelter locations and identify shelter type. Since we could not visually confirm snake locations, we methodically eliminated all possible locations where the snake could possibly be while at close range with the minimum possible gain on the radio receiver.Telemetered kraits tended to be inactive and sheltering underground during the daylight, thus we were confident that our diurnal location checks would not affect their movements. However, in some cases we resorted to determining an individual’s location via triangulation, where multiple lines cast from different vantage points towards the snake intersect on the snake’s location on the GPS, allowing us to determine the animal’s coordinate location from approximately 10–30 m away. This helped ensure that we recorded locations with greater accuracy when snakes sheltered underneath large buildings, as it allowed us to move away from large structures that hindered the GPS accuracy. This technique was also implemented during some nocturnal location checks when a snake was believed to be active among dense vegetation, in an attempt to prevent disturbance of the animals’ natural behavior. While we did hope to gain visual observations of active individuals during the night, we exercised more caution during nocturnal location checks, typically maintaining a minimum distance of approximately 5 m in attempt to lessen the chances of disturbing an active individual’s behavior. If the animal was active we recorded the animal’s observed behavioral state (i.e., moving, feeding, or foraging). When the radio signal was stable and the individual was not visible, we recorded the animal’s behavior as “sheltering”. We strived for an accuracy of  5 m difference), and land-use type (e.g., mixed deciduous forest, human-settlement, semi-natural area, agriculture, plantation; see Supp. Figs. 3 and 4 for photos of land-use types), behavior (e.g., sheltering, moving, foraging, or feeding), and shelter type (e.g., anthropogenic, burrow, or unknown, note we also recorded if we suspected the shelter to be part of a termite tunnel complex due to a close proximity to a visible termite mound; Supp. Fig. 5).During each location check we recorded the straight-line distance between the current and previous locations (distance moved/step length) with the GPS device. We then used step-lengths to summarize their movements by estimating the mean daily displacement (MDD; the total distance moved divided by the number of days the snake was located) and mean movement distance (MMD; the mean relocation distance, excludes distances ≤ 5). In order to limit biases due to some snakes being located multiple times within a given day/night, we limited our sample for estimating MMD and MDD to only include a single location per day. This was accomplished by manually removing “extra” nocturnal location checks that occurred within the same day, making sure to have all shelter relocations present within the dataset. When calculating MDD, we used the total number of daily location checks rather than the number of days between the individual’s tracking start and stop date since there were some days where individuals were not tracked. We also used the same one location check per day dataset to calculate movement/relocation probabilities and to examine each individual’s MMD, MDD, and relocation probability for the overall tracking duration as well as for each season.When feasible, we positioned a Bushnell (Bushnell Corporation, Overland Park, Kansas) time lapse field camera (Trophy Cam HD Essential E3, Model:119837) with infrared night capability on a tripod spaced 2–5 m from occupied shelter sites. We positioned the cameras so that we may gather photos of the focal snake as it exited the shelter site and/or behaviors exhibited near the shelter. We programmed the cameras using a combined setting, including field scan, which continuously captured one photo every minute, along with a motion sensor setting, which took photos upon movement trigger outside of the regular 1-min intervals.Space use and site fidelityAll analyses and most visualizations were done in R v.4.0.5 using RStudio v.1.4.1106 43,44. We attempted to estimate home ranges for the telemetered B. candidus individuals using autocorrelated kernel density estimates (AKDEs) using R package ctmm v.0.6.045,46 in order to better understand the spatial requirements of B. candidus. However, examination of the variograms revealed that the majority of the variograms had not fully stabilized (i.e., limited evidence of range stability in our sample), and many individuals had extremely low effective sample sizes (21.82 ± 9.75, range = 1.49–135.75; Supp. Table 4). Therefore, we do not report home ranges in this text, as the AKDE estimates would violate the assumption of range residency and either underestimate or misrepresent B. candidus spatial requirements. We also examined the speed estimates resulting from fitted movement models. Resulting variograms and tentative home range estimates are included in a supplementary file for viewing only (Supp. Fig. 6, Supp. Table 4). The original code is from Montaño et al.47.Since our data was not sufficient to estimate home range size for the telemetered B. candidus, we instead used Dynamic Brownian Bridge Movement Models (dBBMMs) with the R package move v.4.0.648 to estimate within study occurrence distributions. We caution readers that these are not home range estimates but instead modeling the potential movement pathways animals could have traversed49. Use of dBBMMs not only allows us to estimate occurrence distributions for each individual, thus helping us better understand the animal’s movement pathways and resource use, but it also allows us to examine movement patterns through dBBMM derived motion variance50,51. We selected a window size of 19 and margin size of 5, to catch short resting periods with the margin, while the window size of 19 is long enough to get a valid estimate of motion variance when the animals exhibit activity/movement. Contours however are somewhat arbitrary; therefore, we used three different contours levels (90%, 95%, 99%) to estimate dBBMM occurrence distributions (using R packages adehabitatHR v.0.4.19, and rgeos v.0.5.5), and show the sensitivity to contour choice52,53.All movement data, either including initial capture locations or beginning with the first location check ~ 24 h post release, was used for production of both the AKDEs and dBBMMs for each individual. We also estimated dBBMM occurrence distributions for each telemetered individual with the exception of M29, which only made three small moves within a burrow complex during the short time he was radio-tracked before transmitter failure.We compared space use estimates to two previously published B. candidus tracking datasets34,36, and one unpublished dataset shared on the Zenodo data repository54, all originating from the Sakaerat Biosphere Reserve (approximately 41 km to the south of our study site): two adult males from within the forested area of the reserve [one tracked every 27.8 ± 0.99 h over a period of 103 days, the other tracked every 38.63 ± 11.2 h over a period of 30.58 days]34,54, and a juvenile male from agriculture on a forest boundary [tracked every 50.19 ± h for 66.91 days]36. The previous studies on B. candidus only tracked the movements of a single individual each, had coarser tracking regimes, and used traditional—fundamentally flawed methods55,56—to estimate space use34,36. Therefore, we ran dBBMMs with these previous datasets using the same window (19) and margin size (5).To quantify site reuse and time spent at sites (residency time) we used recursive analysis with the R package recurse v.1.1.257. We defined each site as a circular area with a radius of 5 m around each unique location (matching the targeted GPS accuracy). Then we calculated each individual’s overall number of relocations, each individual’s total number of relocations to each site, and each individual’s site revisit frequency and residency time at each unique site. Then we plotted revisited locations on a land-use map with space use estimates (95% and 99% dBBMM) in an attempt to help identify and highlight activity centers for telemetered individuals (see Supp. Figs. 7–13). All maps were created using Quantum Geographic Information System (QGIS v.3.8.2).Habitat selectionWe used Integrated Step Selection Function models (ISSF) to examine the influence of land-use features on the movements of B. candidus at both the individual and population levels. We included movement data from all male individuals that used more than one habitat feature in our ISSF analysis. Therefore, we excluded F16 and M29 who both only used settlement habitat. Excluding M29 was justified by the individual having been tracked for the shortest duration (19 days) and had the fewest number of moves (n = 3), thus there were not enough relocations for ISSF models to work effectively. Using modified code from Smith et al.51 that used ISSF with Burmese python radio-telemetry data, we used the package amt v.0.1.458 to run ISSF for each individual, with Euclidean distance to particular land-use features (natural areas, agriculture, settlement, buildings, and roads) to determine association or avoidance of features. Cameron Hodges created all land-use shape files in QGIS by digitizing features from satellite imagery and verified all questionable satellite land-use types via on-ground investigation.The semi-natural areas, plantations, mixed deciduous forest and water bodies (such as irrigation canals and ponds which have densely vegetated edges) were all combined into a single layer of less-disturbed habitats which we refer to as “natural areas”. All feature raster layers were then converted into layers with a gradient of continuous values of Euclidean distances to the land-use features, and were inverted in order to avoid zero-inflation of distance to feature values and to make the resulting model directional effects easier to more intuitive. We were able to generate 200 random steps per each observed step (following Smith et al.51), due to the coarse temporal resolution of manually collected radio-telemetry data (i.e., we were not computational limited when deciding the number of random locations). Higher numbers of random steps are preferable as they can aid in detecting smaller effects and rarer landscape features59.To investigate individual selection, we created nine different models testing for association to habitat features, with one being a null model which solely incorporated step-length and turning angle to predict movement60, five examining land-use features individually (agriculture, buildings, settlement, natural areas, roads), and the other three being multi-factor models. Each model considers distance to a land-use variable, step-length, and turn-angle as an aspect of the model. After running each of the nine models for each individual, we then examined the AIC for each model, point estimates (with lower and upper confidence intervals), and p-values in order to identify the best models for each individual and determine the strongest relationships and trends among the samples. We considered models with ∆ AIC  More

In this picture, taken in February at Copenhagen Zoo, I’m holding a vacuum device equipped with a tiny fan and filter. The devices — we call them air samplers — are designed to collect DNA samples from the air. We deployed three samplers at the zoo: one in a stable with two okapi (Okapia johnstoni) and two duikers, one in a rainforest house and one outside, near an exhibit of animals that live in the African savannah.At best, we had hoped to detect nearby animals in small enclosures — an okapi in the stable, for instance. But as we reported in Current Biology, the devices outperformed our expectations (C. Lynggaard Curr. Biol. 32, 701–707; 2022). They picked up identifiable DNA from 49 vertebrates, including guppies in the rainforest pool, ostriches and giraffes in the savannah area, and even cats and dogs in the park next door. Interestingly, we didn’t get any signal from turtles in the rainforest house. Maybe turtles mostly keep their DNA to themselves.Our analysis ultimately found that the sampler could detect animals from nearly 200 metres away. The giraffe in the picture is standing much closer than is necessary for collection of a sample.Airborne DNA is all around us. Birds release skin cells when they flap their wings. Saliva from all sorts of animals can become airborne. Animals release DNA when they defecate. In November 2021, I received a grant to start a research group whose goal is to collect airborne DNA in nature. This approach could transform conservation biology and species monitoring. We could detect rare animals and get a better understanding of diversity without disturbing an environment.We have so many lines of inquiry in this work. The location of the samplers, the rates of air flow, the time, the best methods for sorting DNA from the sample — we’re still trying to work all of these out. We hadn’t expected that the zoo experiment would ever work, so we’re scrambling to plan the next steps. It’s an exciting time. More

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Google Scholar  More