Ecology

Parker, J. & Johnston, L. A. The proximate determinants of insect size. J. Biol. 5, 15 (2006).PubMed 
PubMed Central 
Article 

Google Scholar 
Honěk, A. Intraspecific variation in body size and fecundity in insects: A general relationship. Oikos 66, 483 (1993).Article 

Google Scholar 
Kingsolver, J. G. & Huey, R. B. Size, temperature, and fitness: Three rules. Evol. Ecol. Res. 10, 251–268 (2008).
Google Scholar 
Beukeboom, L. W. Size matters in insects—An introduction. Entomol. Exp. Appl. 166, 2–3 (2018).Article 

Google Scholar 
West, S. A., Flanagan, K. E. & Godfray, H. C. J. The relationship between parasitoid size and fitness in the field, a study of Achrysocharoides zwoelferi (Hymenoptera: Eulophidae). J. Anim. Ecol. 65, 631–639 (1996).Article 

Google Scholar 
Sagarra, L. A., Vincent, C. & Stewart, R. K. Body size as an indicator of parasitoid quality in male and female Anagyrus kamali (Hymenoptera: Encyrtidae). Bull. Entomol. Res. 91, 363–367 (2001).CAS 
PubMed 
Article 

Google Scholar 
Ellers, J., Alphen, J. J. M. V. & Sevenster, J. G. A field study of size–fitness relationships in the parasitoid Asobara tabida. J. Anim. Ecol. 67, 318–324 (1998).Article 

Google Scholar 
Armbruster, P. & Hutchinson, R. A. Pupal mass and wing length as indicators of fecundity in Aedes albopictus and Aedes geniculatus (Diptera: Culicidae). J. Med. Entomol. 39, 699–704 (2002).PubMed 
Article 

Google Scholar 
Tantawy, A. O. & Vetukhiv, M. O. Effects of size on fecundity, longevity and viability in populations of Drosophila pseudoobscura. Am. Nat. 94, 395–403 (1960).Article 

Google Scholar 
Lefranc, A. & Bundgaard, J. The influence of male and female body size on copulation duration and fecundity in Drosophila melanogaster. Hereditas 132, 243–247 (2004).Article 

Google Scholar 
Atkinson, D. Temperature and organism size: A biological law for ectotherms? Adv. Ecol. Res. 25, 1–58 (1994).Article 

Google Scholar 
Poças, G. M., Crosbie, A. E. & Mirth, C. K. When does diet matter? The roles of larval and adult nutrition in regulating adult size traits in Drosophila melanogaster. J. Insect Physiol. 139, 104051. https://doi.org/10.1016/j.jinsphys.2020.104051 (2020).CAS 
Article 
PubMed 

Google Scholar 
Tammaru, T. Determination of adult size in a folivorous moth: constraints at instar level? Ecol. Entomol. 23, 80–89 (1998).Article 

Google Scholar 
Miller, R. S. & Thomas, J. L. The effects of larval crowding and body size on the longevity of adult Drosophila melanogaster. Ecology 39, 118–125 (1958).Article 

Google Scholar 
Nijhout, H. F. The control of body size in insects. Dev. Biol. 261, 1–9 (2003).CAS 
PubMed 
Article 

Google Scholar 
Shingleton, A. W., Mirth, C. K. & Bates, P. W. Developmental model of static allometry in holometabolous insects. Proc. R. Soc. B 275, 1875–1885 (2008).PubMed 
PubMed Central 
Article 

Google Scholar 
Koenraadt, C. J. M. Pupal dimensions as predictors of adult size in fitness studies of Aedes aegypti (Diptera: Culicidae). J. Med. Entomol. 45, 331–336 (2008).CAS 
PubMed 
Article 

Google Scholar 
Stillwell, R. C., Dworkin, I., Shingleton, A. W. & Frankino, W. A. Experimental manipulation of body size to estimate morphological scaling relationships in Drosophila. JoVE 56, 3162. https://doi.org/10.3791/3162 (2011).Article 

Google Scholar 
Shin, S.-M., Akram, W. & Lee, J.-J. Effect of body size on energy reserves in Culex pipiens pallens females (Diptera: Culicidae). Entomol. Res. 42, 163–167 (2012).Article 

Google Scholar 
Mirth, C. K. & Riddiford, L. M. Size assessment and growth control: How adult size is determined in insects. BioEssays 29, 344–355 (2007).CAS 
PubMed 
Article 

Google Scholar 
Chown, S. L. & Gaston, K. J. Body size variation in insects: A macroecological perspective. Biol. Rev. 85, 139–169 (2010).PubMed 
Article 

Google Scholar 
Beadle, G. W., Tatum, E. L. & Clancy, C. W. Food level in relation to rate of development and eye pigmentation in Drosophila melanogaster. Biol. Bull. 75, 447–462 (1938).Article 

Google Scholar 
Gayon, J. History of the concept of allometry1. Am. Zool. 40, 748–758 (2000).
Google Scholar 
Takken, W. et al. Larval nutrition differentially affects adult fitness and Plasmodium development in the malaria vectors Anopheles gambiae and Anopheles stephensi. Parasit. Vectors 6, 345 (2013).PubMed 
PubMed Central 
Article 

Google Scholar 
Briegel, H. Metabolic relationship between female body size, reserves, and fecundity of Aedes aegypti. J. Insect Physiol. 36, 165–172 (1990).Article 

Google Scholar 
Ellers, J. Fat and eggs: An alternative method to measure the trade-off between survival and reproduction in insect parasitoids. Neth. J. Zool. 3, 227–235 (1996).
Google Scholar 
González-Tokman, D. et al. Energy storage, body size and immune response of herbivore beetles at two different elevations in Costa Rica. Rev. Biol. Trop. 67, 608–620 (2019).
Google Scholar 
Timmermann, S. E. & Briegel, H. Larval growth and biosynthesis of reserves in mosquitoes. J. Insect Physiol. 45, 461–470 (1999).CAS 
PubMed 
Article 

Google Scholar 
Strohm, E. Factors affecting body size and fat content in a digger wasp. Oecologia 123, 184–191 (2000).PubMed 
Article 
ADS 

Google Scholar 
Lease, H. M. & Wolf, B. O. Lipid content of terrestrial arthropods in relation to body size, phylogeny, ontogeny and sex. Physiol. Entomol. 36, 29–38 (2011).CAS 
Article 

Google Scholar 
Arrese, E. L. & Soulages, J. L. Insect fat body: Energy, metabolism, and regulation. Annu. Rev. Entomol. 55, 207–225 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kühnlein, R. P. Lipid droplet-based storage fat metabolism in Drosophila. J. Lipid Res. 53, 1430–1436 (2012).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Church, R. B. & Robertson, F. W. A biochemical study of the growth of Drosophila melanogaster. J. Exp. Zool. 162, 337–351 (1966).Article 

Google Scholar 
Merkey, A. B., Wong, C. K., Hoshizaki, D. K. & Gibbs, A. G. Energetics of metamorphosis in Drosophila melanogaster. J. Insect Physiol. 57, 1437–1445 (2011).CAS 
PubMed 
Article 

Google Scholar 
Nestel, D., Tolmasky, D., Rabossi, A. & Quesada-Allué, L. A. Lipid, carbohydrates and protein patterns during metamorphosis of the Mediterranean fruit fly, Ceratitis capitata (Diptera: Tephritidae). Ann. Entomol. Soc. Am. 96, 237–244 (2003).CAS 
Article 

Google Scholar 
Lee, K. P. & Jang, T. Exploring the nutritional basis of starvation resistance in Drosophila melanogaster. Funct. Ecol. 28, 1144–1155 (2014).Article 

Google Scholar 
Hahn, D. A. & Denlinger, D. L. Meeting the energetic demands of insect diapause: Nutrient storage and utilization. J. Insect Physiol. 53, 760–773 (2007).CAS 
PubMed 
Article 

Google Scholar 
Tejeda, M. T. et al. Effects of size, sex and teneral resources on the resistance to hydric stress in the tephritid fruit fly Anastrepha ludens. J. Insect Physiol. 70, 73–80 (2014).CAS 
PubMed 
Article 

Google Scholar 
Hoffmann, A. A., Hallas, R., Anderson, A. R. & Telonis-Scott, M. Evidence for a robust sex-specific trade-off between cold resistance and starvation resistance in Drosophila melanogaster. J. Evol. Biol. 18, 804–810 (2005).CAS 
PubMed 
Article 

Google Scholar 
Alaux, C., Ducloz, F., Crauser, D. & Le Conte, Y. Diet effects on honeybee immunocompetence. Biol. Lett. 6, 562–565 (2010).PubMed 
PubMed Central 
Article 

Google Scholar 
Bryk, B., Hahn, K., Cohen, S. M. & Teleman, A. A. MAP4K3 regulates body size and metabolism in Drosophila. Dev. Biol. 344, 150–157 (2010).CAS 
PubMed 
Article 

Google Scholar 
Gasser, M., Kaiser, M., Berrigan, D. & Stearns, S. C. Life-history correlates of evolution under high and low adult mortality. Evolution 54, 1260–1272 (2000).CAS 
PubMed 
Article 

Google Scholar 
Chippindale, A. K., Chu, T. J. F. & Rose, M. R. Complex trade-offs and the evolution of starvation resistance in Drosophila melanogaster. Evolution 50, 753 (1996).PubMed 
Article 

Google Scholar 
Kristensen, T. N., Overgaard, J., Loeschcke, V. & Mayntz, D. Dietary protein content affects evolution for body size, body fat and viability in Drosophila melanogaster. Biol. Lett. 7, 269–272 (2011).PubMed 
Article 

Google Scholar 
Juarez-Carreño, S. et al. Body-fat sensor triggers ribosome maturation in the steroidogenic gland to initiate sexual maturation in Drosophila. Cell Rep. 37, 109830 (2021).PubMed 
Article 
CAS 

Google Scholar 
Markow, T. A. The secret lives of Drosophila flies. Elife 4, e06793 (2015).PubMed Central 
Article 

Google Scholar 
Choma, M. A., Suter, M. J., Vakoc, B. J., Bouma, B. E. & Tearney, G. J. Physiological homology between Drosophila melanogaster and vertebrate cardiovascular systems. Dis. Model. Mech. 4, 411–420 (2011).CAS 
PubMed 
Article 

Google Scholar 
Morgan, T. H., Sturtevant, A. H., Muller, H. J. & Bridges, C. B. The Mechanism of Mendelian Heredity (H. Holt, 1923).
Google Scholar 
Dobzhansky, T. The influence of the quantity and quality of chromosomal material on the size of the cells in Drosophila melanogaster. Wilhelm Roux Arch. Entwickl Mech. Org. 115, 363–379 (1929).PubMed 
Article 

Google Scholar 
Musselman, L. P. & Kühnlein, R. P. Drosophila as a model to study obesity and metabolic disease. J. Exp. Biol. 221, 163881 (2018).Article 

Google Scholar 
DiAngelo, J. R. & Birnbaum, M. J. Regulation of fat cell mass by insulin in Drosophila melanogaster. Mol. Cell. Biol. 29, 6341–6352 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rovenko, B. M. et al. High sucrose consumption promotes obesity whereas its low consumption induces oxidative stress in Drosophila melanogaster. J. Insect Physiol. 79, 42–54 (2015).CAS 
PubMed 
Article 

Google Scholar 
Hardy, C. M. et al. Obesity-associated cardiac dysfunction in starvation-selected Drosophila melanogaster. Am. J. Physiol.-Regul. Integr. Compar. Physiol. 309, R658–R667 (2015).CAS 
Article 

Google Scholar 
Hardy, C. M. et al. Genome-wide analysis of starvation-selected Drosophila melanogaster—A genetic model of obesity. Mol. Biol. Evol. 35, 50–65 (2018).CAS 
PubMed 
Article 

Google Scholar 
Musselman, L. P. et al. A high-sugar diet produces obesity and insulin resistance in wild-type Drosophila. Dis. Model. Mech. 4, 842–849 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Henry, Y., Renault, D. & Colinet, H. Hormesis-like effect of mild larval crowding on thermotolerance in Drosophila flies. J. Exp. Biol. 221, 169342 (2018).Article 

Google Scholar 
Bulletin, E. P. P. O. Drosophila suzukii. EPPO Bull. 43, 417–424 (2013).Article 

Google Scholar 
Bächli, G., Vilela, C. R., Escher, S. A. & Saura, A. The Drosophilidae (Diptera) of Fennoscandia and Denmark (Brill Academic Publishers, 2004).Book 

Google Scholar 
Markow, T. A. & O’Grady, P. M. Drosophila: A Guide to Species Identification and Use (Elsevier, 2006).
Google Scholar 
Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).CAS 
PubMed 
Article 

Google Scholar 
Visser, B. et al. Variation in lipid synthesis, but genetic homogeneity, among Leptopilina parasitic wasp populations. Ecol. Evol. 8, 7355–7364 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Williams, C. M., Thomas, R. H., MacMillan, H. A., Marshall, K. E. & Sinclair, B. J. Triacylglyceride measurement in small quantities of homogenised insect tissue: Comparisons and caveats. J. Insect Physiol. 57, 1602–1613 (2011).CAS 
PubMed 
Article 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, 2020).Fox, J. & Weisberg, S. An R Companion to Applied Regression 2nd edn. (Sage, 2011).
Google Scholar 
Lenth, R., Singmann, H., Love, J., Buerkner, P. & Herve, M. Emmeans: Estimated marginal means, aka least-squares means. R Package Version 1, 3 (2018).
Google Scholar 
Burnham, K. P. & Anderson, D. R. A practical information-theoretic approach. In Model Selection and Multimodel Inference (ed. Burnham, K. P.) (Springer, 2002).MATH 

Google Scholar 
Crawley, M. J. The R Book (Wiley, 2007).MATH 
Book 

Google Scholar 
Borash, D. J. & Ho, G. T. Patterns of selection: Stress resistance and energy storage in density-dependent populations of Drosophila melanogaster. J. Insect Physiol. 47, 1349–1356 (2001).CAS 
PubMed 
Article 

Google Scholar 
Klepsatel, P., Procházka, E. & Gáliková, M. Crowding of Drosophila larvae affects lifespan and other life-history traits via reduced availability of dietary yeast. Exp. Gerontol. 110, 298–308 (2018).PubMed 
Article 

Google Scholar 
Henry, Y., Overgaard, J. & Colinet, H. Dietary nutrient balance shapes phenotypic traits of Drosophila melanogaster in interaction with gut microbiota. Comp. Biochem. Physiol. A: Mol. Integr. Physiol. 241, 110626 (2020).CAS 
Article 

Google Scholar 
Ireland, S. & Turner, B. The effects of larval crowding and food type on the size and development of the blowfly, Calliphora vomitoria. Forensic Sci. Int. 159, 175–181 (2006).PubMed 
Article 

Google Scholar 
Saunders, D. S. & Bee, A. Effects of larval crowding on size and fecundity of the blow fly, Calliphora vicina (Diptera: Calliphoridae). EJE 92, 615–622 (2013).
Google Scholar 
Ziegler, R. Changes in lipid and carbohydrate metabolism during starvation in adult Manduca sexta. J. Comp. Physiol. B 161, 125–131 (1991).CAS 
PubMed 
Article 

Google Scholar 
Ojeda-Avila, T., Arthur Woods, H. & Raguso, R. A. Effects of dietary variation on growth, composition, and maturation of Manduca sexta (Sphingidae: Lepidoptera). J. Insect Physiol. 49, 293–306 (2003).CAS 
PubMed 
Article 

Google Scholar 
Borash, D. J., Gibbs, A. G., Joshi, A. & Mueller, L. D. A genetic polymorphism maintained by natural selection in a temporally varying environment. Am. Nat. 151, 148. https://doi.org/10.1086/286108 (1998).CAS 
Article 
PubMed 

Google Scholar 
Klepsatel, P., Knoblochová, D., Girish, T. N., Dircksen, H. & Gáliková, M. The influence of developmental diet on reproduction and metabolism in Drosophila. BMC Evol. Biol. 20, 93 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Matzkin, L. M., Johnson, S., Paight, C., Bozinovic, G. & Markow, T. A. Dietary protein and sugar differentially affect development and metabolic pools in ecologically diverse Drosophila. J. Nutr. 141, 1127–1133 (2011).CAS 
PubMed 
Article 

Google Scholar 
Musselman, L. P. et al. Role of fat body lipogenesis in protection against the effects of caloric overload in Drosophila. J. Biol. Chem. 288, 8028–8042 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Reeve, M. W., Fowler, K. & Partridge, L. Increased body size confers greater fitness at lower experimental temperature in male Drosophila melanogaster. J. Evol. Biol. 13, 836–844 (2000).Article 

Google Scholar 
Lounibos, L. P. et al. Does temperature affect the outcome of larval competition between Aedes aegypti and Aedes albopictus?. J. Vector Ecol. 27, 86–95 (2002).CAS 
PubMed 

Google Scholar 
Bergland, A. O., Genissel, A., Nuzhdin, S. V. & Tatar, M. Quantitative trait loci affecting phenotypic plasticity and the allometric relationship of ovariole number and thorax length in Drosophila melanogaster. Genetics 180, 567–582 (2008).PubMed 
PubMed Central 
Article 

Google Scholar 
Holm, S. et al. A comparative perspective on longevity: The effect of body size dominates over ecology in moths. J. Evol. Biol. 29, 2422–2435 (2016).CAS 
PubMed 
Article 

Google Scholar 
Nunney, L. The response to selection for fast larval development in Drosophila melanogaster and its effect on adult weight: An example of a fitness trade-off. Evolution 50, 1193–1204 (1996).PubMed 
Article 

Google Scholar 
Partridge, L. & Farquhar, M. Lifetime mating success of male fruitflies (Drosophila melanogaster) is related to their size. Anim. Behav. 31, 871–877 (1983).Article 

Google Scholar 
Markow, T. A. & Ricker, J. P. Male size, developmental stability, and mating success in natural populations of three Drosophila species. Heredity 69, 122–127 (1992).PubMed 
Article 

Google Scholar 
Wikelski, M. & Romero, L. M. Body size, performance and fitness in galapagos marine iguanas. Integr. Comp. Biol. 43, 376–386 (2003).PubMed 
Article 

Google Scholar 
van Buskirk, J. & Crowder, L. B. Life-history variation in marine turtles. Copeia 1994, 66–81 (1994).Article 

Google Scholar 
Broderick, A. C., Glen, F., Godley, B. J. & Hays, G. C. Variation in reproductive output of marine turtles. J. Exp. Mar. Biol. Ecol. 288, 95–109 (2003).Article 

Google Scholar 
Wauters, L. A. et al. Effects of spatio-temporal variation in food supply on red squirrel Sciurus vulgaris body size and body mass and its consequences for some fitness components. Ecography 30, 51–65 (2007).Article 

Google Scholar 
Lindström, J. Early development and fitness in birds and mammals. Trends Ecol. Evol. 14, 343–348 (1999).PubMed 
Article 

Google Scholar 
Reim, C., Teuschl, Y. & Blanckenhorn, W. U. Size-dependent effects of temperature and food stress on energy reserves and starvation resistance in yellow dung flies. Evol. Ecol. Res. 8, 1215–1234 (2006).
Google Scholar 
Kölliker-Ott, U. M., Blows, M. W. & Hoffmann, A. A. Are wing size, wing shape and asymmetry related to field fitness of Trichogramma egg parasitoids? Oikos 100, 563–573 (2003).Article 

Google Scholar 
Knapp, M. Relative importance of sex, pre-starvation body mass and structural body size in the determination of exceptional starvation resistance of Anchomenus dorsalis (Coleoptera: Carabidae). PLoS ONE 11, e0151459 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Lue, C.-H. et al. DROP: Molecular voucher database for identification of Drosophila parasitoids. Mol. Ecol. Resour. 21, 2437–2454 (2021).CAS 
PubMed 
Article 

Google Scholar 
Visser, B. et al. Loss of lipid synthesis as an evolutionary consequence of a parasitic lifestyle. Proc. Natl. Acad. Sci. 107, 8677–8682 (2010).CAS 
PubMed 
PubMed Central 
Article 
ADS 

Google Scholar 
Visser B et al. Why do
many parasitoids lack adult triglyceride accumulation, despite functioning fatty acid biosynthesis machinery? EcoEvoRxiv:
https://doi.org/10.32942/osf.io/zpf4jArakawa, R., Miura, M. & Fujita, M. Effects of host species on the body size, fecundity, and longevity of Trissolcus mitsukurii (Hymenoptera: Scelionidae), a solitary egg parasitoid of stink bugs. Appl. Entomol. Zool. 39, 177–181 (2004).Article 

Google Scholar 
Visser, B., Alborn, H.T., Rondeaux, S. et al. Phenotypic plasticity explains apparent reverse evolution of fat synthesis in parasitic
wasps. Sci Rep 11, 7751 (2021). https://doi.org/10.1038/s41598-021-86736-8.Krüger, A. P. et al. Effects of irradiation dose on sterility induction and quality parameters of Drosophila suzukii (Diptera: Drosophilidae). J. Econ. Entomol. 111, 741–746 (2018).PubMed 
Article 

Google Scholar 
Nikolouli, K. et al. Sterile insect technique and Wolbachia symbiosis as potential tools for the control of the invasive species Drosophila suzukii. J. Pest Sci. 91, 1–15 (2017).
Google Scholar 
Nikolouli, K., Sassù, F., Mouton, L., Stauffer, C. & Bourtzis, K. Combining sterile and incompatible insect techniques for the population suppression of Drosophila suzukii. J. Pest Sci. 93, 647–661 (2020).CAS 
Article 

Google Scholar 
Calkins, C. O. & Parker, A. G. Sterile insect quality. In Sterile Insect Technique (eds Dyck, V. A. et al.) 269–296 (Springer, 2005).Chapter 

Google Scholar  More

IPCC. Summary for policymakers. In: Climate change 2021: the physical science basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Masson-Delmotte, V, Zhai P, Pirani A, Connors S L, Péan C, Berger S, et al. Editors. Cambridge, UK and New York, NY, USA: Cambridge University Press. 2021. https://doi.org/10.1017/9781009157896.Zhang J, Zhang S, Cheng M, Jiang H, Zhang X, Peng C, et al. Effect of drought on agronomic traits of rice and wheat: a meta-analysis. Int J Environ Res Public Health. 2018;15:839.Article 

Google Scholar 
FAO. Global agriculture towards 2050, high-level expert forum, how to feed the world 2050. Rome: Food and Agriculture Organization of United Nations FAO. 2009.Agoussar A, Yergeau E. Engineering the plant microbiota in the context of the theory of ecological communities. Curr Opin Biotechnol. 2021;70:220–5.CAS 
Article 

Google Scholar 
Quiza L, St-Arnaud M, Yergeau E. Harnessing phytomicrobiome signaling for rhizosphere microbiome engineering. Front Plant Sci. 2015;6:507.Article 

Google Scholar 
Mitter B, Pfaffenbichler N, Flavell R, Compant S, Antonielli L, Petric A, et al. A new approach to modify plant microbiomes and traits by introducing beneficial bacteria at flowering into progeny seeds. Front Microbiol. 2017;8. https://doi.org/10.3389/fmicb.2017.00011.Li X, Jousset A, de Boer W, Carrión VJ, Zhang T, Wang X, et al. Legacy of land use history determines reprogramming of plant physiology by soil microbiome. ISME J. 2019;13:738–51.CAS 
Article 

Google Scholar 
Nelson EB, Simoneau P, Barret M, Mitter B, Compant S. Editorial special issue: the soil, the seed, the microbes and the plant. Plant Soil. 2018;422:1–5.CAS 
Article 

Google Scholar 
Trivedi P, Leach JE, Tringe SG, Sa T, Singh BK. Plant–microbiome interactions: from community assembly to plant health. Nat Rev Microbiol. 2020;18:607–21.CAS 
Article 

Google Scholar 
Moroenyane I, Tremblay J, Yergeau E. Soybean microbiome recovery after disruption is modulated by the seed and not the soil microbiome. Phytobiomes J. 2021;5:418–31.Article 

Google Scholar 
Xiong C, Zhu Y-G, Wang J-T, Singh B, Han L-L, Shen J-P, et al. Host selection shapes crop microbiome assembly and network complexity. New Phytol. 2021;229:1091–104.CAS 
Article 

Google Scholar 
Schimel J, Balser TC, Wallenstein M. Microbial stress-response physiology and its implications for ecosystem function. Ecology. 2007;88:1386–94.Article 

Google Scholar 
Allison SD, Martiny JBH. Resistance, resilience, and redundancy in microbial communities. Proc Natl Acad Sci USA. 2008;105:11512–9.CAS 
Article 

Google Scholar 
Bouskill NJ, Lim HC, Borglin S, Salve R, Wood TE, Silver WL, et al. Pre-exposure to drought increases the resistance of tropical forest soil bacterial communities to extended drought. ISME J. 2013;7:384–94.CAS 
Article 

Google Scholar 
Evans SE, Wallenstein MD. Soil microbial community response to drying and rewetting stress: does historical precipitation regime matter? Biogeochemistry. 2012;109:101–16.Article 

Google Scholar 
Meisner A, Snoek BL, Nesme J, Dent E, Jacquiod S, Classen AT, et al. Soil microbial legacies differ following drying-rewetting and freezing-thawing cycles. ISME J. 2021;15:1207–21.CAS 
Article 

Google Scholar 
Azarbad H, Constant P, Giard-Laliberté C, Bainard LD, Yergeau E. Water stress history and wheat genotype modulate rhizosphere microbial response to drought. Soil Biol Biochem. 2018;126:228–36.CAS 
Article 

Google Scholar 
Jones DL, Nguyen C, Finlay RD. Carbon flow in the rhizosphere: carbon trading at the soil–root interface. Plant Soil. 2009;321:5–33.CAS 
Article 

Google Scholar 
Sasse J, Martinoia E, Northen T. Feed your friends: do plant exudates shape the root microbiome? Trends Plant Sci. 2018;23:25–41.CAS 
Article 

Google Scholar 
Holz M, Zarebanadkouki M, Kuzyakov Y, Pausch J, Carminati A. Root hairs increase rhizosphere extension and carbon input to soil. Ann Bot. 2018;121:61–9.CAS 
Article 

Google Scholar 
Marasco R, Rolli E, Ettoumi B, Vigani G, Mapelli F, Borin S, et al. A drought resistance-promoting microbiome is selected by root system under desert farming. PLoS ONE. 2012;7:e48479.CAS 
Article 

Google Scholar 
Moroenyane I, Mendes L, Tremblay J, Tripathi B, Yergeau É. Plant compartments and developmental stages modulate the balance between niche-based and neutral processes in soybean Microbiome. Microb Ecol. 2021;82:416–28. https://doi.org/10.1007/s00248-021-01688-w.Chaparro JM, Badri DV, Vivanco JM. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 2014;8:790–803.CAS 
Article 

Google Scholar 
Azarbad H, Tremblay J, Giard-Laliberté C, Bainard LD, Yergeau E. Four decades of soil water stress history together with host genotype constrain the response of the wheat microbiome to soil moisture. FEMS Microbiol Ecol. 2020;96. https://doi.org/10.1093/femsec/fiaa098.Chen S, Cade-Menun BJ, Bainard LD, St. Luce M, Hu Y, Chen Q. The influence of long-term N and P fertilization on soil P forms and cycling in a wheat/fallow cropping system. Geoderma. 2021;404:115274.CAS 
Article 

Google Scholar 
Smith EG, Zentner RP, Campbell CA, Lemke R, Brandt K. Long-term crop rotation effects on production, grain quality, profitability, and risk in the northern great plains. Agron J. 2017;109:957–67.Article 

Google Scholar 
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.Article 

Google Scholar 
Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest package: tests in linear mixed effects models. J Stat Softw. 2017;82:1–26.Article 

Google Scholar 
Yang S. otuSummary: summarizing OTU table regarding the composition, abundance and beta diversity of abundant and rare biospheres. 2018.Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package. R package version 2.5-6. 2019. Available online at: https://CRAN.R-project.org/package=vegan.Wagner MR, Lundberg DS, del Rio TG, Tringe SG, Dangl JL, Mitchell-Olds T. Host genotype and age shape the leaf and root microbiomes of a wild perennial plant. Nat Commun. 2016;7:12151.CAS 
Article 

Google Scholar 
Hardoim PR, Hardoim CCP, Overbeek LS, van, Elsas JD. Dynamics of seed-borne rice endophytes on early plant growth stages. PLoS ONE. 2012;7:e30438.CAS 
Article 

Google Scholar 
Grady KL, Sorensen JW, Stopnisek N, Guittar J, Shade A. Assembly and seasonality of core phyllosphere microbiota on perennial biofuel crops. Nat Commun. 2019;10:4135.Article 

Google Scholar 
Walsh CM, Becker-Uncapher I, Carlson M, Fierer N. Variable influences of soil and seed-associated bacterial communities on the assembly of seedling microbiomes. ISME J. 2021;15:2748–62.Morales Moreira ZP, Helgason BL, Germida JJ Environment has a Stronger Effect than Host Plant Genotype in Shaping Spring Brassica napus Seed Microbiomes. Phytobiomes J. 2021:PBIOMES-08-20-0059-R.Abdullaeva Y, Ambika Manirajan B, Honermeier B, Schnell S, Cardinale M. Domestication affects the composition, diversity, and co-occurrence of the cereal seed microbiota. J Adv Res. 2021;31:75–86.CAS 
Article 

Google Scholar 
Chandrasekaran M, Chanratana M, Kim K, Seshadri S, Sa T. Impact of arbuscular mycorrhizal fungi on photosynthesis, water status, and gas exchange of plants under salt stress–a meta-analysis. Front Plant Sci. 2019;10. https://doi.org/10.3389/fpls.2019.00457.Begum N, Qin C, Ahanger MA, Raza S, Khan MI, Ashraf M, et al. Role of arbuscular mycorrhizal fungi in plant growth regulation: implications in abiotic stress tolerance. Front Plant Sci. 2019;10. https://doi.org/10.3389/fpls.2019.01068.Rashid MI, Mujawar LH, Shahzad T, Almeelbi T, Ismail IMI, Oves M. Bacteria and fungi can contribute to nutrients bioavailability and aggregate formation in degraded soils. Microbiol Res. 2016;183:26–41.CAS 
Article 

Google Scholar 
Crowther TW, Thomas SM, Maynard DS, Baldrian P, Covey K, Frey SD, et al. Biotic interactions mediate soil microbial feedbacks to climate change. Proc Natl Acad Sci USA. 2015;112:7033–8.CAS 
Article 

Google Scholar 
Evans SE, Wallenstein MD. Climate change alters ecological strategies of soil bacteria. Ecol Lett. 2014;17:155–64.Article 

Google Scholar 
Santos-Medellín C, Edwards J, Liechty Z, Nguyen B, Sundaresan V. Drought stress results in a compartment-specific restructuring of the rice root-associated microbiomes. mBio. 2017;8:e00764–17.Article 

Google Scholar 
Naylor D, Coleman-Derr D. Drought stress and root-associated bacterial communities. Front Plant Sci. 2018;8:2223.Article 

Google Scholar 
Zhou C, Ma Z, Zhu L, Xiao X, Xie Y, Zhu J, et al. Rhizobacterial Strain Bacillus megaterium BOFC15 induces cellular polyamine changes that improve plant growth and drought resistance. Int J Mol Sci. 2016;17:E976.Article 

Google Scholar 
Bokhari A, Essack M, Lafi FF, Andres-Barrao C, Jalal R, Alamoudi S, et al. Bioprospecting desert plant Bacillus endophytic strains for their potential to enhance plant stress tolerance. Sci Rep. 2019;9:18154.CAS 
Article 

Google Scholar 
Xu L, Naylor D, Dong Z, Simmons T, Pierroz G, Hixson KK, et al. Drought delays development of the sorghum root microbiome and enriches for monoderm bacteria. Proc Natl Acad Sci USA. 2018;115:E4284–93.CAS 
PubMed 
PubMed Central 

Google Scholar 
Na X, Cao X, Ma C, Ma S, Xu P, Liu S, et al. Plant stage, not drought stress, determines the effect of cultivars on bacterial community diversity in the rhizosphere of broomcorn millet (Panicum miliaceum L.). Front Microbiol. 2019;10. https://doi.org/10.3389/fmicb.2019.00828. More

Experiment 1To study the acceptance and tag retention rates of honey bees under different introduction conditions, we set up three two-frame observation hives with (sim 1500) adult bees and a queen. Observation hives were set up in a shed, with an entrance tube that connected to the outside 4 ft above the ground so bees could freely forage in the surrounding fields (Fig. 1A). We put emerging brood from healthy source colonies in an incubator ((33.5^{circ }) C, (ge 55)% RH) and tagged individuals that emerged overnight. Each hive had two vent holes (1” Dia.) through which we could introduce bees (Fig. 1A).Figure 1(A) Observation hive with introduction holes (red), through which bees were introduced via funnel or introduction cage. (B) Plastic tags, silicon tags, and sucrose spray. (C) Photograph of a tagged bee foraging (Photo by Greg Yauney).Full size imageIn our initial experiment, there were seven treatment groups with 20 bees each per colony (n = 420 bees total). The seven treatments were: Control (C), No Sucrose (NS), Plastic (P), Wood Glue (WG), Not Incubated (NI), No Cage (NC), and Day (D). The Control treatment was designed to be a positive control, where we applied all the techniques we thought might increase acceptance of bees into a colony and tag retention rates. All other treatments had a single difference in the tag, tagging process, or introduction process to distinguish it from the control group, as detailed in Table 1. We glued a tag to the thorax of each bee (Fig. 1B) and marked the abdomen with a paint pen (Posca) to distinguish among treatment groups. In order to glue tags on, we picked each bee up, placed a small amount of glue on the thorax, and placed a tag on top of the glue with a pair of forceps (see Video 1 which details the tagging process). All bees except those in the Plastic group were tagged with 1.7 mm(^2) silicon tags (3.4 mm area). Silicon was chosen because it is a material representative of ASICs, which you would expect in a custom chip designed to track bee foraging flights. Plastic tags were 3 mm Dia. plastic discs (7.07 mm area), which are the commercially available bee tags commonly used in honey bee tracking and behavior experiments (Betterbee). All tags were glued on with shellac glue, the glue that comes with commerical honey bee marking kits (Betterbee), except for in the wood glue group, where they were glued on with wood glue (Titebond III). Next, bees were placed in a container with a bit of honey and stored until they were ready to be introduced. All bees except those in the Not Incubated group were placed in the incubator ((33.5^{circ }) C, (ge 55)% RH). The Not Incubated group was stored in a room environment, with variation between 21–27(^{circ }) C and 35–42% RH until introduction. Bees in the Day treatment spent 5 h in the incubator and then were sprayed with a light sucrose syrup (1 sucrose: 1 water (v/v)) and introduced at 4pm while the hives were still actively foraging. The rest of the bees spent between 5 and 8 h in the incubator or room environment before being introduced at 10:30 pm, after foraging had concluded. All except the No Sucrose group were sprayed with a light sucrose syrup before being introduced. The No Cage bees were rapidly introduced through one of the vent holes on the top of the hive using a funnel. The rest of the bees were placed in a cage together, which we connected to the introduction holes at the top of the colony, allowing them to move freely between the cage and the hive.Table 1 Experimental design used for preparation and introduction of treatment groups.Full size tableBeginning on day 2 (07/09/2020), we observed each hive in the morning on days 2-4 and 6-9 to see how many bees per group were present, hereinafter referred to as presence, and how many bees per group were present with tags, hereinafter referred to as success. We selected a random order in which to observe the three hives and a random order in which to observe the treatment groups for each hive. Each side of each hive had a grid drawn on it that divided it into nine squares. We scanned each side of each colony by eye for each treatment, starting with the lower left square of the grid on the first side, moving across the row, and then moving up to the next row, counting presence and success, using a tally counter when needed. We then moved rapidly to the other side and started at the top left of the grid, scanning row by row until we had observed each square in the grid. After an initial scan for each treatment, we placed the covers on the hives and shook for 10s to encourage bees to move around in the hive, and then waited for at least 15 minutes before a second observation. The maximum presence and success from the two daily observations were used for each treatment group and hive for analysis. Since we collected data by scanning each colony, we sometimes found more bees from a group in an observation hive than we had found in the same hive on previous day(s), even though more time had passed. Over the course of the experiment, our hives grew in size, and we believed we were seeing less tagged bees in part because they made up a smaller proportion of the hive population, and so decided to do a destructive sampling before the tagged bees reached foraging age. After dark on day 14 (7/21/2020), we made sure no tagged bees were dead on the bottom of the hives. We blocked the entrances, vacuumed all bees at the entrances into containers, and froze vacuumed bees and the three colonies, so that we could do a destructive sampling of all 3 colonies. This allowed us to get a final count of the presence and success for each of the seven treatment groups. We dissected each frozen colony, removing and inspecting each dead bee, and recorded the presence and success of each treatment group.Experiments 2 and 3We set up three two-frame observation hives in the same shed used for experiment 1 to conduct follow-up experiments in August 2020. The goal of experiment 2 was to compare Gorillaglue gel, an easily accessible Superglue (SG), to Titebond III, a readily accessible Wood Glue (WG2) used in experiment 1. We placed frames of capped brood in an incubator overnight to produce one day old nurse bees. We picked up each bee, placed a small dot of either superglue or wood glue on the thorax, and then placed 1.7 mm(^2) silicon tags on top of the glue. Bees were stored in the incubator with honey for 5–6 h until after dark. Then, we sprayed the bees with a light sucrose syrup and connected their cages to the vent holes at the top of the observation hives, allowing the bees to freely move between their cage and the hives. These details are summarized in Table 1.Some honey bee tagging projects may benefit from tagging foragers as opposed to nurse bees, because nurse bees are the youngest workers and if you tag them you must wait for them to reach foraging age, during which time they may lose their tags. Specifically, tagging foragers as opposed to nurses will be advantageous when the tag price is extremely high or the project is very time constrained, and knowing the exact age of tagged bees is not important for the project goals. Since foragers are older workers that have already acquired the colony scent and learned to navigate the area surrounding their hive, the optimal methods for introducing nurses and foragers may differ. It is not easy to use bees from a source colony, because if they are within foraging range of their maternal colony, they will attempt to fly back home. The goal of experiment 3 was to apply a treatment that had high success with nurse bees (Experiment 1: WG) to foragers, and compare with releasing foragers near their colony and allowing them to return freely. We call these treatments Hive Introduced (HI) and Natural Release (NR), respectively. All foragers for this experiment were collected from the observation hives and were introduced back to the same observation hive after tagging, either through the vent holes at the top of the hive or by releasing the bees near the entrance of the hive. We collected foragers from each colony entrance into a cage with an insect vacuum (Hand-Held DC Vac/Aspirator, Bioquip), specifically aspirating bees that were arriving from foraging trips or had nectar loads, and placed them in the fridge to anesthetize them. We then selected those with intact wings, placed a dot of wood glue on their thoraxes, and placed silicon tags on top of the glue. Both treatment groups were stored in the incubator ((33.5^{circ }) C, (ge 55)% RH) and given honey to feed on. After 2 h in the incubator, the containers with NR bees were sprayed with a light sucrose syrup and placed on the ground 5 ft in front of their respective hive entrances and opened, allowing the bees to fly back to their hives unaided. At 10PM, when it was dark and foraging had concluded, the HI bees were sprayed with a light sucrose syrup. Their cages were then connected to the vent holes at the top of the observation hives, allowing them to freely move between their cage and the hives.Experiments 2 and 3 were conducted in the same hives simultaneously, but were considered separate experiments because experiment 2 was conducted with nurses of known age and experiment 3 was conducted with foragers of unknown age. Nurses and foragers typically have an age difference and experience different levels of risk due to the behaviors they engage in, and so we analyzed these data separately in order to not confound our results. Beginning on day two (08/26/2020), we observed each hive on days 2–11 and 15–21 to determine introduction presence and success for experiment 2 and experiment 3. Forager observations (experiment 3) were always done early in the morning, before foraging activity commenced. As in experiment 1, we randomized observation order, scanned colonies for each treatment group before and after shaking, and used the maximum presence and success from the two observations for analysis. Since we collected data on multiple days by scanning each colony, we occasionally found more bees in a group than we had found on previous day(s), even though more time had passed.Statistical methodsStatistical analyses were performed in R 4.0.520. To determine which preparation and introduction techniques were associated with the highest presence and success, we built generalized linear mixed-effects models (glmms)21 for the proportion of present and success bees to introduced bees respectively, with treatment and sampling day as fixed effects, and colony as a random effect. For experiment 1, treatment was a categorical variable, where the Control bees were the reference group. We assessed the significance of the full models using Wald likelihood ratio chi-square tests on each glmm (‘Anova’ function in the ‘car’ package with test set to ‘Chisq’)22. In all statistical tests, (alpha) was set to 0.05. The destructive data from experiment 1 were analyzed separately from hive observation data. We ran a correlation test to determine the relationship between hive observation data from the final observation day, day 9, and the destructive sampling on day 14 using the ggpubr package23. More

Brierley AS, Kingsford MJ. Impacts of climate change on marine organisms and ecosystems. Curr Biol. 2009;19:R602–R614.CAS 
PubMed 
Article 

Google Scholar 
Gissi E, Manea E, Mazaris AD, Fraschetti S, Almpanidou V, Bevilacqua S, et al. A review of the combined effects of climate change and other local human stressors on the marine environment. Sci Total Environ. 2021;755:142564.CAS 
PubMed 
Article 

Google Scholar 
Carella F, Antuofermo E, Farina S, Salati F, Mandas D, Prado P, et al. In the wake of the ongoing mass mortality events: co-occurrence of Mycobacterium, Haplosporidium and other pathogens in Pinna nobilis collected in Italy and Spain (Mediterranean Sea). Front Mar Sci. 2020;7:48.Article 

Google Scholar 
Seuront L, Nicastro KR, Zardi GI, Goberville E. Decreased thermal tolerance under recurrent heat stress conditions explains summer mass mortality of the blue mussel Mytilus edulis. Sci Rep. 2019;9:17498.PubMed 
PubMed Central 
Article 

Google Scholar 
Fey SB, Siepielski AM, Nussle S, Cervantes-Yoshida K, Hwan JL, Huber ER, et al. Recent shifts in the occurrence, cause, and magnitude of animal mass mortality events. Proc Natl Acad Sci USA. 2015;112:1083–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Scarpa F, Sanna D, Azzena I, Mugetti D, Cerruti F, Hosseini S, et al. Multiple non-species-specific pathogens possibly triggered the mass mortality in Pinna nobilis. Life. 2020;10:238.CAS 
PubMed Central 
Article 

Google Scholar 
Bradley M, Kutz SJ, Jenkins E, O’Hara TM. The potential impact of climate change on infectious diseases of Arctic fauna. Int J Circumpolar Health. 2005;64:468–77.PubMed 
Article 

Google Scholar 
Beyer J, Green NW, Brooks S, Allan IJ, Ruus A, Gomes T, et al. Blue mussels (Mytilus edulis spp.) as sentinel organisms in coastal pollution monitoring: a review. Mar Environ Res. 2017;130:338–65.CAS 
PubMed 
Article 

Google Scholar 
Siravegna G, Marsoni S, Siena S, Bardelli A. Integrating liquid biopsies into the management of cancer. Nat Rev Clin Oncol. 2017;14:531–48.CAS 
PubMed 
Article 

Google Scholar 
Wan JCM, Massie C, Garcia-Corbacho J, Mouliere F, Brenton JD, Caldas C, et al. Liquid biopsies come of age: towards implementation of circulating tumour DNA. Nat Rev Cancer. 2017;17:223–38.CAS 
PubMed 
Article 

Google Scholar 
Mandel P, Metais P. Nuclear acids in human blood plasma. Comptes Rendus Séances Soc Biol Filiales. 1948;142:241–3.CAS 

Google Scholar 
Bronkhorst AJ, Ungerer V, Holdenrieder S. The emerging role of cell-free DNA as a molecular marker for cancer management. Biomol Detect Quantif. 2019;17:100087.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ignatiadis M, Sledge GW, Jeffrey SS. Liquid biopsy enters the clinic – implementation issues and future challenges. Nat Rev Clin Oncol. 2021;18:297–312.PubMed 
Article 

Google Scholar 
Lo YM, Corbetta N, Chamberlain PF, Rai V, Sargent IL, Redman CW, et al. Presence of fetal DNA in maternal plasma and serum. Lancet. 1997;350:485–7.CAS 
PubMed 
Article 

Google Scholar 
Moufarrej MN, Wong RJ, Shaw GM, Stevenson DK, Quake SR. Investigating pregnancy and its complications using circulating cell-free RNA in women’s blood during gestation. Front Pediatr. 2020;8:605219.PubMed 
PubMed Central 
Article 

Google Scholar 
Oellerich M, Sherwood K, Keown P, Schutz E, Beck J, Stegbauer J, et al. Liquid biopsies: donor-derived cell-free DNA for the detection of kidney allograft injury. Nat Rev Nephrol. 2021;17:591–603.CAS 
PubMed 
Article 

Google Scholar 
Wong FC, Lo YM. Prenatal diagnosis innovation: genome sequencing of maternal plasma. Annu Rev Med. 2016;67:419–32.CAS 
PubMed 
Article 

Google Scholar 
Gu W, Deng X, Lee M, Sucu YD, Arevalo S, Stryke D, et al. Rapid pathogen detection by metagenomic next-generation sequencing of infected body fluids. Nat Med. 2021;27:115–24.CAS 
PubMed 
Article 

Google Scholar 
Huang YF, Chen YJ, Fan TC, Chang NC, Chen YJ, Midha MK, et al. Analysis of microbial sequences in plasma cell-free DNA for early-onset breast cancer patients and healthy females. BMC Med Genom. 2018;11:16.Article 
CAS 

Google Scholar 
Goggs R, Jeffery U, LeVine DN, Li RHL. Neutrophil-extracellular traps, cell-free DNA, and immunothrombosis in companion animals: a review. Vet Pathol. 2020;57:6–23.CAS 
PubMed 
Article 

Google Scholar 
Kowarsky M, De Vlaminck I, Okamoto J, Neff NF, LeBreton M, Nwobegabay J, et al. Cell-free DNA reveals potential zoonotic reservoirs in non-human primates. BioRxiv. 2018;481093.Caza F, Bernet E, Veyrier FJ, Betoulle S, St-Pierre Y. Hemocytes released in seawater act as Trojan horses for spreading of bacterial infections in mussels. Sci Rep. 2020;10:19696.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Andrew S. FastQC: a quality control tool for high throughput sequence data. 2010. http://www.bioinformatics.babraham.ac.uk/projects/fastqc.Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Morgulis A, Gertz EM, Schäffer AA, Agarwala R. A fast and symmetric DUST implementation to mask low-complexity DNA sequences. Comput Biol. 2006;13:1028–40.CAS 
Article 

Google Scholar 
Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.CAS 
PubMed 
Article 

Google Scholar 
Wood DE, Salzberg SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014;15:R46.PubMed 
PubMed Central 
Article 

Google Scholar 
Cuccuru G, Orsini M, Pinna A, Sbardellati A, Soranzo N, Travaglione A, et al. Orione, a web-based framework for NGS analysis in microbiology. Bioinformatics. 2014;30:1928–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ondov BD, Bergman NH, Phillippy AM. Interactive metagenomic visualization in a Web browser. BMC Bioinform. 2011;12:385.Article 

Google Scholar 
Lüskow F, Riisgård H. In situ filtration rates of blue mussels (Mytilus edulis) measured by an open-top chamber method. OJMS. 2018;8:395–406.Article 

Google Scholar 
Szpechcinski A, Struniawska R, Zaleska J, Chabowski M, Orlowski T, Roszkowski K, et al. Evaluation of fluorescence-based methods for total vs. amplifiable DNA quantification in plasma of lung cancer patients. J Physiol Pharmacol. 2008;59:675–81.PubMed 

Google Scholar 
Tissot C, Toffart AC, Villar S, Souquet PJ, Merle P, Moro-Sibilot D, et al. Circulating free DNA concentration is an independent prognostic biomarker in lung cancer. Eur Respir J. 2015;46:1773–80.CAS 
PubMed 
Article 

Google Scholar 
Kustanovich A, Schwartz R, Peretz T, Grinshpun A. Life and death of circulating cell-free DNA. Cancer Biol Ther. 2019;20:1057–67.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Prouteau A, Denis JA, De Fornel P, Cadieu E, Derrien T, Kergal C, et al. Circulating tumor DNA is detectable in canine histiocytic sarcoma, oral malignant melanoma, and multicentric lymphoma. Sci Rep. 2021;11:877.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Vandewoestyne M, Van Hoofstat D, Franssen A, Van Nieuwerburgh F, Deforce D. Presence and potential of cell free DNA in different types of forensic samples. For Sci Int Genet. 2013;7:316–20.CAS 
Article 

Google Scholar 
Kowarsky M, Camunas-Soler J, Kertesz M, De Vlaminck I, Koh W, Pan W, et al. Numerous uncharacterized and highly divergent microbes which colonize humans are revealed by circulating cell-free DNA. Proc Natl Acad Sci USA. 2017;114:9623–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Meddeb R, Dache ZAA, Thezenas S, Otandault A, Tanos R, Pastor B, et al. Quantifying circulating cell-free DNA in humans. Sci Rep. 2019;9:5220.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Li YF, Yang N, Liang X, Yoshida A, Osatomi K, Power D, et al. Elevated seawater temperatures decrease microbial diversity in the gut of Mytilus coruscus. Front Physiol. 2018;9:839.PubMed 
PubMed Central 
Article 

Google Scholar 
Musella M, Wathsala R, Tavella T, Rampelli S, Barone M, Palladino G, et al. Tissue-scale microbiota of the Mediterranean mussel (Mytilus galloprovincialis) and its relationship with the environment. Sci Total Environ. 2020;717:137209.CAS 
PubMed 
Article 

Google Scholar 
Thompson JR, Randa MA, Marcelino LA, Tomita-Mitchell A, Lim E, Polz MF. Diversity and dynamics of a north atlantic coastal Vibrio community. Appl Environ Microbiol. 2004;70:4103–10.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pfister CA, Meyer F, Antonopoulos DA. Metagenomic profiling of a microbial assemblage associated with the California mussel: a node in networks of carbon and nitrogen cycling. PLoS One. 2010;5:e10518.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Galand PE, Casamayor EO, Kirchman DL, Potvin M, Lovejoy C. Unique archaeal assemblages in the Arctic Ocean unveiled by massively parallel tag sequencing. ISME J. 2009;3:860–9.CAS 
PubMed 
Article 

Google Scholar 
Korzhenkov AA, Toshchakov SV, Bargiela R, Gibbard H, Ferrer M, Teplyuk AV, et al. Archaea dominate the microbial community in an ecosystem with low-to-moderate temperature and extreme acidity. Microbiome. 2019;7:11.PubMed 
PubMed Central 
Article 

Google Scholar 
Spain EA, Johnson SC, Hutton B, Whittaker JM, Lucieer V, Watson SJ, et al. Shallow seafloor gas emissions near Heard and McDonald Islands on the Kerguelen Plateau, southern Indian Ocean. Earth Space Sci. 2020;7:e2019EA000695.Article 

Google Scholar 
Farías L, Florez-Leiva L, Besoain V, Sarthou G, Fernández C. Dissolved greenhouse gases (nitrous oxide and methane) associated with the naturally iron-fertilized Kerguelen region (KEOPS 2 cruise) in the Southern Ocean. Biogeosciences. 2015;12:1925–40.Article 

Google Scholar 
Legendre M, Bartoli J, Shmakova L, Jeudy S, Labadie K, Adrait A, et al. Thirty-thousand-year-old distant relative of giant icosahedral DNA viruses with a pandoravirus morphology. Proc Natl Acad Sci USA. 2014;111:4274–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Levasseur A, Andreani J, Delerce J, Bou Khalil J, Robert C, La Scola B, et al. Comparison of a modern and fossil pithovirus reveals its genetic conservation and evolution. Genome Biol Evol. 2016;8:2333–9.PubMed 
PubMed Central 
Article 

Google Scholar 
Kelley JL, Brown AP, Therkildsen NO, Foote AD. The life aquatic: advances in marine vertebrate genomics. Nat Rev Genet. 2016;17:523–34.CAS 
PubMed 
Article 

Google Scholar 
Colmer SF, Luethy D, Abraham M, Stefanovski D, Hurcombe SD. Utility of cell-free DNA concentrations and illness severity scores to predict survival in critically ill neonatal foals. PLoS One. 2021;16:e0242635.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rushton JG, Ertl R, Klein D, Tichy A, Nell B. Circulating cell-free DNA does not harbour a diagnostic benefit in cats with feline diffuse iris melanomas. J Feline Med Surg. 2019;21:124–32.PubMed 
Article 

Google Scholar 
Tagawa M, Shimbo G, Inokuma H, Miyahara K. Quantification of plasma cell-free DNA levels in dogs with various tumors. J Vet Diagn Investig. 2019;31:836–43.CAS 
Article 

Google Scholar 
Shi J, Zhang R, Li J, Zhang R. Size profile of cell-free DNA: a beacon guiding the practice and innovation of clinical testing. Theranostics. 2020;10:4737–48.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Fernando MR, Jiang C, Krzyzanowski GD, Ryan WL. Analysis of human blood plasma cell-free DNA fragment size distribution using EvaGreen chemistry based droplet digital PCR assays. Clin Chim Acta. 2018;483:39–47.CAS 
PubMed 
Article 

Google Scholar 
Findlay AJ. Microbial impact on polysulfide dynamics in the environment. FEMS Microbiol Lett. 2016;363:fnw103.PubMed 
Article 
CAS 

Google Scholar 
Jørgensen BB, Findlay AJ, Pellerin A. The biogeochemical sulfur cycle of marine sediments. Front Microbiol. 2019;10:849.PubMed 
PubMed Central 
Article 

Google Scholar 
Teske A, Brinkhoff T, Muyzer G, Moser DP, Rethmeier J, Jannasch HW. Diversity of thiosulfate-oxidizing bacteria from marine sediments and hydrothermal vents. Appl Environ Microbiol. 2000;66:3125–33.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zhang X, Du Z, Zheng R, Luan Z, Qi F, Cheng K, et al. Development of a new deep-sea hybrid Raman insertion probe and its application to the geochemistry of hydrothermal vent and cold seep fluids. Deep Sea Res Part I Oceanogr Res Pap. 2017;123:1–12.Article 
CAS 

Google Scholar 
Egger M, Riedinger N, Mogollón JM, Jørgensen BB. Global diffusive fluxes of methane in marine sediments. Nat Geosci. 2018;11:421–5.CAS 
Article 

Google Scholar 
Ansorge R, Romano S, Sayavedra L, Kupczok A, Tegetmeyer HE, Dubilier N, et al. Functional diversity enables multiple symbiont strains to coexist in deep-sea mussels. Nat Microbiol. 2019;4:2487–97.PubMed 
Article 
CAS 

Google Scholar 
Russell SL, Pepper-Tunick E, Svedberg J, Byrne A, Ruelas Castillo J, Vollmers C, et al. Horizontal transmission and recombination maintain forever young bacterial symbiont genomes. PLoS Genet. 2020;16:e1008935.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Angly FE, Felts B, Breitbart M, Salamon P, Edwards RA, Carlson C, et al. The marine viromes of four oceanic regions. PLoS Biol. 2006;4:e368.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Li Z, Pan D, Wei G, Pi W, Zhang C, Wang JH, et al. Deep sea sediments associated with cold seeps are a subsurface reservoir of viral diversity. ISME J. 2021;15:2366–78.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Thongsripong P, Chandler JA, Kittayapong P, Wilcox BA, Kapan DD, Bennett SN. Metagenomic shotgun sequencing reveals host species as an important driver of virome composition in mosquitoes. Sci Rep. 2021;11:8448.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Koonin EV, Krupovic M, Agol VI. The Baltimore classification of viruses 50 years later: how does it stand in the light of virus evolution? Microbiol Mol Biol Rev. 2021;85:e0005321.PubMed 
Article 

Google Scholar 
Koonin EV, Dolja VV, Krupovic M, Varsani A, Wolf YI, Yutin N, et al. Global organization and proposed megataxonomy of the virus world. Microbiol Mol Biol Rev. 2020;84:e00061–19.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Breitbach S, Tug S, Simon P. Circulating cell-free DNA: an up-coming molecular marker in exercise physiology. Sports Med. 2012;42:565–86.PubMed 
Article 

Google Scholar 
Preissner KT, Herwald H. Extracellular nucleic acids in immunity and cardiovascular responses: between alert and disease. Thromb Haemost. 2017;117:1272–82.PubMed 
Article 

Google Scholar 
Schwarzenbach H. Circulating nucleic acids as biomarkers in breast cancer. Breast Cancer Res. 2013;15:211.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Murphy DJ. Freezing resistance in intertidal invertebrates. Annu Rev Physiol. 1983;45:289–99.CAS 
PubMed 
Article 

Google Scholar 
Robledo JAF, Yadavalli R, Allam B, Espinosa EP, Gerdol M, Greco S, et al. From the raw bar to the bench: bivalves as models for human health. Dev Comp Immunol. 2019;92:260–82.Article 

Google Scholar 
Cowart DA, Murphy KR, Cheng CC. Metagenomic sequencing of environmental DNA reveals marine faunal assemblages from the West Antarctic Peninsula. Mar Genom. 2018;37:148–60.Article 

Google Scholar 
Parducci L, Bennett KD, Ficetola GF, Alsos IG, Suyama Y, Wood JR, et al. Ancient plant DNA in lake sediments. New Phytol. 2017;214:924–42.CAS 
PubMed 
Article 

Google Scholar 
Mariani S, Baillie C, Giuliano C, Riesgo A. Sponges as natural environmental DNA samplers. Curr Biol. 2019;29:R401–R402.CAS 
PubMed 
Article 

Google Scholar 
Weber S, Brink L, Wörner M, Künzel S, Veith M, Teubner D, et al. Molecular diet analysis in zebra and quagga mussels (Dreissena spp.) and an assessment of the utility of aquatic filter feeders as biological eDNA filters. BioRxiv. 2021; 432951.Caza F, Joly de Boissel PG, Villemur R, Betoulle S, St-Pierre Y. Liquid biopsies for omics-based analysis in sentinel mussels. Plos One. 2019;14:e0223525.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hunter ME, Ferrante JA, Meigs-Friend G, Ulmer A. Improving eDNA yield and inhibitor reduction through increased water volumes and multi-filter isolation techniques. Sci Rep. 2019;9:5259.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Burkhardt W III, Calci KR. Selective accumulation may account for shellfish-associated viral illness. Appl Environ Microbiol. 2000;66:1375–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Di Girolamo R, Liston J, Matches J. Ionic bonding, the mechanism of viral uptake by shellfish mucus. Appl Environ Microbiol. 1977;33:19–25.PubMed 
PubMed Central 
Article 

Google Scholar 
Metzger MJ, Reinisch C, Sherry J, Goff SP. Horizontal transmission of clonal cancer cells causes leukemia in soft-shell clams. Cell. 2015;161:255–63.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Metzger MJ, Villalba A, Carballal MJ, Iglesias D, Sherry J, Reinisch C, et al. Widespread transmission of independent cancer lineages within multiple bivalve species. Nature. 2016;534:705–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Canesi L, Gallo G, Gavioli M, Pruzzo C. Bacteria–hemocyte interactions and phagocytosis in marine bivalves. Microsc Res Tech. 2002;57:469–76.PubMed 
Article 

Google Scholar 
Andruszkiewicz EA, Koseff JR, Fringer OB, Ouellette NT, Lowe AB, Edwards CA, et al. Modeling environmental DNA transport in the coastal ocean using Lagrangian particle tracking. Front Mar Sci. 2019;6:477.Article 

Google Scholar 
Wood ZT, Lacoursière-Roussel A, LeBlanc F, Trudel M, Kinnison MT, Garry McBrine C, et al. Spatial heterogeneity of eDNA transport improves stream assessment of threatened salmon presence, abundance, and location. Front Ecol Evol. 2021;9:650717.Article 

Google Scholar 
Rand AC, Jain M, Eizenga JM, Musselman-Brown A, Olsen HE, Akeson M, et al. Mapping DNA methylation with high-throughput nanopore sequencing. Nat Methods. 2017;14:411–3.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Simpson JT, Workman RE, Zuzarte PC, David M, Dursi LJ, Timp W. Detecting DNA cytosine methylation using nanopore sequencing. Nat Methods. 2017;14:407–10.CAS 
PubMed 
Article 

Google Scholar 
Cavalli G, Heard E. Advances in epigenetics link genetics to the environment and disease. Nature. 2019;571:489–99.CAS 
PubMed 
Article 

Google Scholar 
Fan G, Song Y, Yang L, Huang X, Zhang S, Zhang M, et al. Initial data release and announcement of the 10,000 Fish Genomes Project (Fish10K). Gigascience. 2020;9:giaa080.PubMed 
PubMed Central 
Article 

Google Scholar  More

Time dynamics of the mollusk communitiesIn this section, the presence-absence of the species recorded in the three periods (T1, T2, T3) are analyzed in relation to time, habitat, and human impact. The list of the 28 species of freshwater mollusks (17 gastropods and 11 bivalves) in T1–T3, their codes, and origins are given in Table 1.The number of mollusk species has increased in time as the river has shifted from lotic to a mixture of flowing and stagnant sectors due to the building of reservoirs. T1 was characterized mainly by rheophilic elements and prosobranchs. Some species became extinct during the hydro-technical works (before or during T2) and are unlikely to recover, such as the rheophilic Theodoxus transversalis and Lithoglyphus naticoides. Other rheophilic species disappeared between T1 and T2 but managed to survive in tributaries and repopulated some sectors during the last years. The most remarkable recovery is that of the thick-shelled river mussel Unio crassus, a species protected by EU legislation. T2 was characterized by some extinctions but also colonization by lentiphilic pulmonates and tolerant, resistant species such as some clams. A few lotic species also survived in the river sectors between the dams. In T3, we encountered a rich and diverse community, including some newly established populations of AIS and the discontinuous presence of both lentic and lotic communities. Overall, the present-day fauna is richer than in former periods, consisting of 15 species of gastropods and 8 bivalves. The AIS included the gastropods Physa acuta and Ferrissia californica, which arrived in the area most likely during the XXth century, Viviparus acerosus, which is native to the Danube, but unknown until after 2000 in the upper-middle Olt River basin, the bivalves Dreissena polymorpha, also native in the Danube but an invader in the middle Olt since 2008–2010, Sinanodonta woodiana, first found in 2015, and Corbicula fluminea, which was first discovered in the Olt (and also in Transylvania) during our survey in February 2020. The mean number of native species per river’s sector increases almost linearly (2.8 species per sector in T1, 3.3 in T2, and 4.6 in T3), while the corresponding values for AIS increase non-linearly (no AIS in T1, 0.6 species per sector in T2 and 3.2 in T3).In the CCA of freshwater mollusk community changes through time (Period as predictor), the adjusted explained variation was 23.6% (test on all axes, pseudo-F = 5.9, p = 0.001). The polygons delimiting the positions of the sites during the three periods of time show no overlap, and they were distinct and separated in the ordination space (Fig. 1a). T2 (the period with maximum human impact) is distinctly placed and separated from the period without impact (T1) along both ordination axes. Meantime, T3 is closer to T1, having an intermediate position between the other two periods, showing a trend of recovery, such as the return of some species. In the CCA of T1–T3 species presence-absence predicted by the selected environmental descriptors (Period, Habitat, and Impact) (Fig. 1b), the adjusted explained variation was 28.36% (test on all axes, pseudo-F = 4.2, p = 0.001). FD(Rao) computed on all FT was plotted as isolines by GAM on the ordination space (model AIC = -17.19, model test F = 5.1, p = 0.003; tests of non-linearity in predictor effects: F = 3.9, p = 0.03). The functional diversity decreased from T1 to T2, then increased sharply to T3; it also decreased from rivers (R) to lakes (L) and along the human impact gradient (Impact).Figure 1Canonical correspondence analysis (CCA) of mollusk communities: (a) classification diagram of sampling sites based on period (as predictors): T1—XIXth century, T2—1995–2000, T3—2020 (adjusted explained variation 23.6%; first axis accounts for 17.6% the second for 6.0%, both axes are significant, p = 0.001); (b) CCA diagram of species occurrence constrained by environmental predictors (period, habitat: L—lakes, lentic sector in reservoirs, R—river, lotic sectors, and Impact—human impact) with functional diversity expressed as Rao quadratic entropy index (FD (Rao)) isolines plotted by generalized additive models (GAM) on the ordination space (adjusted explained variation 28.36%; first axis accounts for 16.3%, the second for 6.0%, both axes are significant, p = 0.001) .Full size imageIn the dc-CA with the selected predictors on T1–T3 presence-absence data, the first two axes separate the communities by period, each positioned in a distinct quadrant (Fig. 2). After a decrease in diversity from T1 to T2, in T3, there were more species and higher functional diversity. In time, there was a reduction in body size, a switch from species with separate sexes to hermaphrodites, a transition of oviposition towards ovo-viviparity (in snails), and external fecundation (in bivalves), and a switch of the feeding type. The dc-CA adjusted explained variation was 16.47%; tests based on sectors and species showed significant relationships (combined test for all axes, pseudo-F = 2.6, p = 0.006), the dimensionality test based on case scores was significant for the first axis (pseudo-F = 4.2, p = 0.001) and marginally significant for the second one (pseudo-F = 1.1, p = 0.053). In contrast, the dimensionality test based on species scores was significant only for the first axis (pseudo-F = 1.6, p = 0.004). The adjusted variation explained by environmental predictors (Hab, Impact, and Period) was 28.36%, and by the selected functional traits (Sexes, FeedT, SizeM, and Ovipos) was 14.64%.Figure 2Double-constrained correspondence analysis (dc-CA) with selected predictors on presence-absence data in T1–T3. The selected functional traits (in blue) are Sexes (circles: H—hermaphrodite, S—separate sexes), Feeding type (squares: SCR—scraper, SS—scraper and sediment, SF—scraper and filter, F—filter, SEDF—suspension and deposit feeder), Oviposition (diamonds: OV—ovo-viviparity, CAP—capsule/eggmass, BE—parental care, juveniles in brood pouches of demibranchs, No—no oviposition, external fecundation), and mean body size (SizeM); the selected environmental predictors (in red) are time (Period, with levels T1—XIXth century, T2—1995–2000, T3—2020), habitat (R—river, lotic sector; L—lake, a lentic sector in reservoirs) and human impact (Impact). Species are coded by the first three letters of the genus and species names. The adjusted explained variation was 16.47%, the first axis accounts for 12.7% and the second for 2.2%. Native species have black labels, while aliens (AIS) are written in green.Full size imageWe have split the binary data describing communities into two parts: natives and AIS, using the latter as predictors. We partitioned the variation in native species composition explained by the three predictor groups (Period, Environment, and AIS) (Fig. 3), subjecting the explanatory variables to an interactive forward selection procedure. We used RDA with centered response variables (CCA can not be used because the empty rows in some tables hinder the use of a proper hierarchical permutation scheme). The adjusted explained variation was 39.6% (the simple effects: time accounted for 22.33%, habitat and impact 24.73%, and the selected AIS 20.82%). All simple and unique effects were significant (p  More

FAO (Food and Agriculture Organization of the United Nations, Fisheries and Aquaculture Department). The State of World Fisheries and Aquaculture 2020 (Food and Agriculture Organization of the United Nations, 2020).Costello, C. et al. The future of food from the sea. Nature 588(7836), 1–6 (2020).
Google Scholar 
Sarah, A. B. et al. Trends in the detection of aquatic non-indigenous species across global marine, estuarine and freshwater ecosystems: A 50-year perspective. Divers. Distrib. 26(12), 1780–1797 (2020).
Google Scholar 
FAMA (Fisheries Administration of the Ministry of Agriculture of the PRC). China Fishery Statistical Yearbook (China Academic Journal Electronic Publishing House, 1949–1975). https://www.cafs.ac.cn/kxyj/qgyytjnj.htm. (in Chinese).FAMA (Fisheries Administration of the Ministry of Agriculture of the PRC). China Fishery Statistical Yearbook (China Agriculture Press, 2021). (in Chinese).Shelton, W. L. & Rothbard, S. Exotic species in global aquaculture—a review. Isr. J. Aquac. 58(1), 3–28 (2006).
Google Scholar 
Ju, R. et al. Emerging risks of non-native species escapes from aquaculture: Call for policy improvements in China and other developing countries. J. Appl. Ecol. 57, 86–90 (2020).
Google Scholar 
Zhu, C. & Dong S. Aquaculture site selection and carrying capacity management in the People’s Republic of China. In Site Selection and Carrying Capacities for Inland and Coastal Aquaculture (eds Ross, L. G., Telfer, T. C., Falconer, L., Soto, D. & Aguilar Manjarrez, J.) 219–230 (FAO/Institute of Aquaculture, Expert Workshop, 6–8 December 2010, University of Stirling, UK, FAO, Rome, 2013).Falconer, L., Telfer, T. C. & Ross, L. G. Investigation of a novel approach for aquaculture site selection. J. Environ. Manag. 181, 791–804 (2016).
Google Scholar 
Liu, Y. et al. Spatiotemporal variations in suitable areas for Japanese scallop aquaculture in the Dalian coastal area from 2003 to 2012. Aquaculture 422, 172–183 (2014).
Google Scholar 
Reverter, M. et al. Aquaculture at the crossroads of global warming and antimicrobial resistance. Nat. Commun. 11, 1870 (2020).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Eve, G., Marc, B. & Montserrat, R. Immune response of the sea cucumber Parastichopus regalis to different temperatures: Implications for aquaculture purposes. Aquaculture 497, 357–363 (2018).
Google Scholar 
Gentry, R. et al. Mapping the global potential for marine aquaculture. Nat. Ecol. Evol. 1(9), 1317 (2017).PubMed 

Google Scholar 
Kim, B. et al. Impact of seawater temperature on Korean aquaculture under representative concentration pathways (RCP) scenarios. Aquaculture 542(3), 736893 (2021).
Google Scholar 
Wentz, F. J. et al. Satellite measurements of sea surface temperature through clouds. Science 288(5467), 847–850 (2000).ADS 
CAS 
PubMed 

Google Scholar 
Mediodia, H. Effects of sea surface temperature on tuna catch: Evidence from countries in the Eastern Pacific Ocean. Ocean. Coast. Manag. 209, 105657 (2021).
Google Scholar 
Liu, S., Zhang, Z., Wu, J. & Yu W. Spatial and temporal variations of potential habitats of jumbo flying squid Dosidicus gigas off Peru under increasing sea surface Temperature. Fish. Sci. (2020). (in Chinese with an English Abstract).Nian, R. et al. The identification and prediction in abundance variation of Atlantic cod via long short-term memory with periodicity, time–frequency co-movement, and lead-lag effect across sea surface temperature, sea surface salinity, catches, and prey biomass from 1919 to 2016. Front. Mar. Sci. 8, 665716 (2021).
Google Scholar 
Radiarta, I. & Saitoh, S. Biophysical models for Japanese scallop, Mizuhopecten yessoensis, aquaculture site selection in Funka Bay, Hokkaido, Japan, using remotely sensed data and geographic information system. Aquacult. Int. 17(5), 403 (2009).
Google Scholar 
Laama, C. & Bachar, N. Evaluation of site suitability for the expansion of mussel farming in the Bay of Souahlia (Algeria) using empirical models. J. Appl. Aquac. 31(4), 337–355 (2019).
Google Scholar 
Liu, Y. et al. Impact of iceanographic environmental shifts and atmospheric events on the sustainable development of coastal aquaculture: A case study of kelp and scallops in southern Hokkaido, Japan. Sustainability 7(2), 1263–1279 (2015).
Google Scholar 
Foo, S. & Gregory, P. Sea surface temperature in coral reef restoration outcomes. Environ. Res. Lett. 15(7), 074045 (2020).ADS 

Google Scholar 
Warren, D. & Seifert, S. Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria. Ecol. Appl. 21(2), 335–342 (2011).PubMed 

Google Scholar 
Warren, D., Glor, R. & Turelli, M. ENMTools: a toolbox for comparative studies of environmental niche models. Ecography 33(3), 607–611 (2010).
Google Scholar 
Sillero, N. What does ecological modelling model? A proposed classification of ecological niche models based on their underlying methods. Ecol. Model. 222(8), 1343–1346 (2011).
Google Scholar 
Bo, Z., Xin-Jun, C. & Gang, L. Relationship between the resource and fishing ground of mackerel and environmental factors based on GAM and GLM models in the East China Sea and Yellow Sea. Shuichan Xuebao 32(3), 379–386 (2008) (in Chinese with an English abstract).
Google Scholar 
Chen, P. & Chen, X. Analysis of habitat distribution of Argentine shortfin squid (Illex argentinus) in the southwest Atlantic Ocean using maximum entropy model. Shuichan Xuebao 40(6), 893–902 (2016) (in Chinese with an English abstract).
Google Scholar 
Zhang, S., Shi, Y., Li, F., Zhu, M. & Wei, Z. Prediction of potential fishing ground for Pacific saury (Cololabis saira) based on MAXENT model. J. Ocean. Univ. China 29(2), 280–286 (2020) (in Chinese with an English abstract).
Google Scholar 
Phillips, S. & Elith, J. On estimating probability of presence from use–availability or presence–background data. Ecology 94(6), 1409–1419 (2013).PubMed 

Google Scholar 
Yang, H. et al. Current advances and technological prospects of the sea cucumber seed industry in China. Mar. Sci. 7, 2–9 (2020) (in Chinese with an English abstract).CAS 

Google Scholar 
Chang, Y., Ding, J., Song, J. & Yang, W. Biology and Aquaculture of Sea Cucumbers and Sea Urchins (Ocean Press, Beijing, 2004) (in Chinese).
Google Scholar 
Li, C. & Hu, W. Status, trend and countermeasure in development of sea cucumber Apostichopus Japonicus Selenka industry in China. Mar. Econ. China. 1, 3–20 (2017) (in Chinese with an English abstract).
Google Scholar 
FAMA (Fisheries Administration of the Ministry of Agriculture of the PRC). China Fishery Statistical Yearbook. China Agriculture Press; 2003. (in Chinese).FAMA (Fisheries Administration of the Ministry of Agriculture of the PRC). China Fishery Statistical Yearbook (China Agriculture Press, 2012). (in Chinese).He, C. & Huang, G. On Apostichopus japonicus culture in China and major culture provinces. Fish. Inf. St. (2014). (in Chinese with an English abstract).Su, L., Zhou, C., Hu, L. & Xu, J. Development status and sustainable development of Apostichopus japonicus industry in south China. Fish. Sci. Technol. Inf. 2, 57–60 (2014) (in Chinese).
Google Scholar 
Guo, F. Research and analysis report on sea cucumber Apostichopus japonicus aquaculture industry in typical regions of North and South China: A case study of Wafangdian city and Xiapu county. Masteral dissertation, Dalian Ocean University. (2021). (in Chinese with an English abstract).Agatsuma, Y. Strongylocentrotus intermedius. In Sea Urchins: Biology and Ecology 4th edn (ed. Lawrence, J. M.) 609–621 (Elsevier, Amsterdam, 2020).
Google Scholar 
Wang, Z. & Chang, Y. Studies on hatching of Japanese sea urchin Strongylocentrotus intermedius. J. Fish. Sci. C 4(1), 60–67 (1997) (in Chinese with an English abstract).
Google Scholar 
Liao, Y. Fauna of China: Echinodermata: Holothuroidea (Science Press, 1997) (in Chinese).Merckx, B. et al. Null models reveal preferential sampling, spatial autocorrelation and overfitting in habitat suitability modelling. Ecol. Model. 222(3), 588–597 (2011).CAS 

Google Scholar 
Matthew, A. A method for implementing a statistically significant number of data classes in the Jenks algorithm. In Proceedings of the Sixth International Conference on Fuzzy Systems and Knowledge Discovery 35–38 (Tianjin, China. 2009).Zhao, G. Water environment analysis of two typical breeding patterns. Masteral dissertation, Hebei Agricultural University. (2019). (in Chinese with an English abstract).Liu, C. & Lan, Y. Situation and countermeasure of sea cucumber culturing industry in Fujian Province. J. Fujian Fish. (2013). (in Chinese with an English abstract).Fei, G. et al. Effect of water temperature on digestive enzyme activity and gut mass in sea cucumber Apostichopus japonicus (Selenka), with special reference to aestivation. J. Oceanol. Limnol. 27(4), 714–722 (2009).
Google Scholar 
Han, C., Lin, P., et al. A study on key technique of Stichopus japonicus Selenka farming in southern sea area. Mod. Fish. Inf. (2011). (in Chinese with an English abstract).Chang, Y., Wang, Z. & Wang, G. Effect of temperature and algae on feeding and growth in sea urchin, Strongylocentrotus intermedius. J. Fish. China. (1997). (in Chinese with an English abstract).Lawrence, J. et al. Temperature effect of feed consumption, absorption, and assimilation efficiencies and production of the sea urchin Strongylocentrotus intermedius. J. Shellfish Res. 28, 389–395 (2009).
Google Scholar 
Zhao, C. et al. Effects of temperature and feeding regime on food consumption, growth, gonad production and quality of the sea urchin Strongylocentrotus intermedius. J. Mar. Biolog. Assoc. 96(1), 185–195 (2015).
Google Scholar 
Lawrence, J., Zhao, C. & Chang, Y. Large-scale production of sea urchin (Strongylocentrotus intermedius) seed in a hatchery in China. Aquac. Int. 27(1), 1–7 (2019).CAS 

Google Scholar 
Chang, Y. et al. Aquaculture of Strongylocentrotus intermedius in Fujian coastal areas. South China Fish. Sci. 16(3), 1–9 (2020).
Google Scholar 
Yu, Z. Raft culture technique of sea urchin in south China. China Fish. 376(003), 57–57 (2007) (in Chinese).
Google Scholar  More

Site selectionWe selected 104 sites covering a range of ecological, land cover, and climate conditions across North America (Table 1). These sites were selected because they are part of either the National Ecological Observatory Network (NEON) or AmeriFlux network, all have PhenoCams, and each has at least one year of available flux data between 2017 and 2021. Among the included sites, 44 are part of the NEON.Table 1 List of AmeriFlux and NEON sites included in the dataset. .Full size tablePlanetScope image database compilationThe LSP metrics included in the dataset are derived from a database of daily 3 m PlanetScope imagery. To compile this database, a Python script was created to search, request, and download imagery using Planet’s RESTful API interface (https://developers.planet.com/docs/apis/data/). For each site, the area of interest (AOI) was defined using a GeoJSON file that prescribed a 10 by 10 km box centered over the flux tower at each site. Each GeoJSON was then used to submit search requests to the API. As part of the search process, the following filters were applied to ensure that good quality images with adequate clear sky views and high-accuracy geolocation were downloaded: (1) quality category identified as ‘standard’; (2) cloud cover less than or equal to 50%; and (3) ground control is ‘true’. Filtering was performed using all available PlanetScope ‘PSScene4Band’ imagery from 2016 to 2022. Once the API completed the search, the Python script read the search results, submitted orders, and the selected imagery was downloaded from Planet’s cloud-based system to local storage. During execution of the Python script, a log file was created to keep track of successful and failed orders. If an order failed, the script was run again targeting the specific order that failed. The resulting dataset included over 1.8 million unique files with, on average, 3,885 scene images for each site (i.e., the number of images, on average, that overlap part of each 10 by 10 km site), and had a total volume of 62.2 TB.Image processingTo ensure that high-quality image time series were used to generate LSP metrics, we used PlanetScope per-pixel quality assurance information to exclude pixels that had low quality in all 4 bands (i.e., blue, green, red, and near-infrared). Specifically, we excluded pixels where the Unusable Data Mask (layer ‘umd’) was not 0 (i.e., we retained pixels that were not cloud contaminated or located in non-image areas) and pixels where the Usable Data Mask (layer ‘umd2’) is 0 (i.e., we retained pixels that were not contaminated by snow, shadow, haze, or clouds). We then cropped all the images to exclude pixels outside of the 10 by 10 km window centered over each tower. We selected this window size based on published results showing that 80% of the average monthly footprint at eddy covariance towers ranges from 103 to 107 square meters22. Note that the swath for PlanetScope imagery often did not cover entire sites and some sites (e.g., the tall tower at US-Pfa) have larger footprints than other sites. Similarly, most sites had multiple PlanetScope image acquisitions on the same day. To create image time series, we mosaiced all available imagery at each site on each date, and, under the assumption that geolocation error was non-systematic and modest, we created a single image for each date using the mean surface reflectance for pixels with multiple values on the same day. The resulting database of daily surface reflectance images were sorted in chronological order, sub-divided into 200 sub-areas at each site (i.e., 0.5 km2 each), and then stored as image stacks to facilitate parallel processing to estimate LSP metrics, where each image stack included all images from July 1, 2016 through January 31, 2022.Creation of daily EVI2 time seriesTo estimate LSP metrics we adapted the algorithm described by Bolton et al.19, which was originally implemented to estimate LSP metrics from harmonized Landsat and Sentinel-2 (HLS) imagery, for use with PlanetScope imagery. Prior to LSP estimation, daily images of the two-band Enhanced Vegetation Index30 (EVI2) data were generated from PlanetScope imagery and then interpolated to create smooth time series of daily EVI2 values at each pixel in three main steps. First, sources of variation related to clouds, atmospheric aerosols, and snow were detected and removed from the EVI2 time series at each pixel based on data masks provided with PlanetScope imagery (described above) and outlier detection criteria (i.e., de-spiking and removing negative EVI2 values). Second, we identified the ‘background’ EVI2 value (the minimum EVI2 value outside of the growing season) based on the 10th percentile of snow-free EVI2 values at each pixel. Any dates with EVI2 values below the background value were replaced with the background EVI2. Third, penalized cubic smoothing splines were used to gap-fill and smooth the data to create daily EVI2 time series across all years of available data. Complete details on these steps are given in Bolton et al.19. This approach has been tested and shown to yield PlanetScope EVI2 time series that are consistent with both EVI2 time series from HLS imagery and time series of the Green Chromatic Coordinate (GCC) from PhenoCam imagery26. We used the EVI2 instead of other vegetation indices such as the Enhanced Vegetation Index (EVI) or the Normalized Difference Vegetation Index (NDVI) because EVI2 is less sensitive to noise from atmospheric effects relative to EVI and is less prone to saturation over dense canopies and noise from variation in soil background reflectance over sparse canopies relative to the NDVI30,32. Thus, phenological metrics from EVI2 time series tend to have better agreement with PhenoCam observations than corresponding metrics from NDVI33.Identifying phenological cyclesPrior to estimating LSP metrics, we first identity unique growth cycles by searching the period before and after each local peak in the daily PlanetScope EVI2 time series. To be considered a valid growth cycle, the difference in EVI2 between the local minimum and maximum was required to be at least 0.1 and greater than 35% of the total range in EVI2 over the 24-month period centered on the target year ± 6 months. The start of each growth cycle is restricted to occur within 185 days before the peak of the cycle and at least 30 days after the previous peak. The same procedure was applied in reverse at the end of the cycle to constrain the range of end dates for each growth cycle. This procedure is applied recursively over the time series until each local peak has been assessed and all growth cycles (with associated green-up period, peak greenness, and green-down period) are identified in the time series at each pixel. As part of this process, the algorithm provides the number of growth cycles identified for each year in the time series.Retrieving LSP metricsLSP metrics are estimated for each pixel in up to two growth cycles in each year. If no growth cycles are detected, the algorithm returns fill values for all timing metrics, but does report values for the four annual metrics: EVImax, EVIamp, EVIarea, and numObs (see below). If more than two growth cycles are detected, which is rare but does occur (e.g., alfalfa, which is harvested and regrows multiple times in a year), the algorithm records 7 LSP metrics for each of the two growth cycles with the largest EVI2 amplitudes. The resulting dataset includes seven ‘timing’ metrics that identify the timing of greenup onset, mid-greenup, maturity, peak EVI2, greendown onset, mid-greendown, and dormancy. These metrics record the day of year (DOY) when the EVI2 time series exceeds 15%, 50%, and 90% of EVI2 amplitude during the greenup phase, reaches its maximum, and goes below 90%, 50%, and 15% of EVI2 amplitude during the greendown phase. In addition, the algorithm records three complementary metrics that characterize the magnitude of seasonality and total ‘greenness’ at each pixel in each growth cycle: the EVI2 amplitude, the maximum EVI2, and the growing season integral of EVI2, which is calculated as the sum of daily EVI2 values between the growth cycle start- and end-dates (i.e., from greenup onset to dormancy).Quality assurance flagsQuality Assurance (QA) values are estimated at each pixel based on the density of observations and the quality of spline fits during each phenophase of the growing season. A QA value of 1 (high quality) is assigned if the correlation between observed versus fitted daily EVI2 values is greater than 0.75 and the maximum gap during each phase is less than 30 days. A QA value of 2 (moderate quality) is assigned if the correlation coefficient is less than 0.75 or the length of the maximum gap over the segment is greater than 30 days. A QA value of 3 (low quality) is assigned if the correlation coefficient is less than 0.75 and the length of the maximum gap over the segment is greater than 30 days. A QA value of 4 is assigned if no growth cycles were detected or insufficient data were available to run the algorithm. More

Turner, D. & Bateson, P. (eds) The Domestic Cat: The Biology of Its Behaviour (Cambridge Univ. Press, 2000).
Google Scholar 
Bradshaw, J. W. S., Goodwin, D., Legrand-Defrétin, V. & Nott, H. M. R. Food selection by the domestic cat, an obligate carnivore. Comp. Biochem. Physiol. A Physiol. 114, 205–209 (1996).CAS 
PubMed 
Article 

Google Scholar 
Trouwborst, A., McCormack, P. C. & Martínez Camacho, E. Domestic cats and their impacts on biodiversity: A blind spot in the application of nature conservation law. People Nat. 2, 235–250 (2020).Article 

Google Scholar 
Crowley, S. L., Cecchetti, M. & McDonald, R. A. Our wild companions: Domestic cats in the anthropocene. Trends Ecol. Evol. 35, 477–483 (2020).PubMed 
Article 

Google Scholar 
Driscoll, C. A. et al. The Near Eastern origin of cat domestication. Science 317, 519–523 (2007).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Van Neer, W., Linseele, V., Friedman, R. & De Cupere, B. More evidence for cat taming at the Predynastic elite cemetery of Hierakonpolis (Upper Egypt). J. Archaeol. Sci. 45, 103–111 (2014).Article 

Google Scholar 
Ottoni, C. et al. The palaeogenetics of cat dispersal in the ancient world. Nat. Ecol. Evol. 1, 0139 (2017).Article 

Google Scholar 
Baca, M. et al. Human-mediated dispersal of cats in the Neolithic Central Europe. Heredity 121, 557–563 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Vigne, J. The beginning of cat domestication in East and West Asia. Doc. Archaeobiol. 15, 343–354 (2019).
Google Scholar 
Krajcarz, M. et al. Ancestors of domestic cats in Neolithic Central Europe: Isotopic evidence of a synanthropic diet. Proc. Natl. Acad. Sci. USA 117, 17710–17719 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Piontek, A. M. et al. Analysis of cat diet across an urbanisation gradient. Urban Ecosyst. 24, 59–69 (2021).Article 

Google Scholar 
Medina, F. M. et al. A global review of the impacts of invasive cats on island endangered vertebrates. Glob. Chang. Biol. 17, 3503–3510 (2011).ADS 
Article 

Google Scholar 
Moseby, K. E., Peacock, D. E. & Read, J. L. Catastrophic cat predation: A call for predator profiling in wildlife protection programs. Biol. Conserv. 191, 331–340 (2015).Article 

Google Scholar 
Loss, S. R., Will, T. & Marra, P. P. The impact of free-ranging domestic cats on wildlife of the United States. Nat. Commun. 4, 1396 (2013).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Beaumont, M. et al. Genetic diversity and introgression in the Scottish wildcat. Mol. Ecol. 10, 319–336 (2001).CAS 
PubMed 
Article 

Google Scholar 
Beugin, M. P. et al. Hybridization between Felis silvestris silvestris and Felis silvestris catus in two contrasted environments in France. Ecol. Evol. 10, 263–276 (2020).PubMed 
Article 

Google Scholar 
Biró, Z., Lanszki, J., Szemethy, L., Heltai, M. & Randi, E. Feeding habits of feral domestic cats (Felis catus), wild cats (Felis silvestris) and their hybrids: Trophic niche overlap among cat groups in Hungary. J. Zool. 266, 187–196 (2005).Article 

Google Scholar 
Széles, G. L., Purger, J. J., Molnár, T. & Lanszki, J. Comparative analysis of the diet of feral and house cats and wildcat in Europe. Mammal. Res. 63, 43–53 (2018).Article 

Google Scholar 
Ottoni, C. & Van Neer, W. The dispersal of the domestic cat paleogenetic and zooarcheological evidence. Near East. Archaeol. 83, 38–45 (2020).Article 

Google Scholar 
Bitz-Thorsen, J. & Gotfredsen, A. B. Domestic cats (Felis catus) in Denmark have increased significantly in size since the Viking Age. Danish J. Archaeol. 7, 241–254 (2018).Article 

Google Scholar 
Faure, E. & Kitchener, A. C. An archaeological and historical review of the relationships between felids and people. Anthrozoos 22, 221–238 (2009).Article 

Google Scholar 
von den Driesch, A. Kulturgeschichte der Hauskatze. In Krankheiten der Katze, Bd. 1 (eds Schmidt, V. & Horzinek, M. C.) 17–40 (Fischer, 1992).
Google Scholar 
Głażewska, I. & Kijewski, T. A new view on the European feline population from mtDNA analysis in Polish domestic cats. Forensic Sci. Int. Genet. 27, 116–122 (2017).PubMed 
Article 
CAS 

Google Scholar 
Cucchi, T. et al. Tracking the Near Eastern origins and European dispersal of the western house mouse. Sci. Rep. 10, 8276 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Van Klinken, G. J., Richards, M. P. & Hedges, B. E. M. An overview of causes for stable isotopic variations in past European human populations: environmental, ecophysiological, and cultural effects. In Biogeochemical Approaches to Paleodietary Analysis (eds Ambrose, S. & Katzenberg, M.) 39–63 (Kluwer Academic Publishers, 2002). https://doi.org/10.1007/0-306-47194-9_3.Chapter 

Google Scholar 
Drucker, D. G., Bridault, A., Hobson, K. A., Szuma, E. & Bocherens, H. Can carbon-13 in large herbivores reflect the canopy effect in temperate and boreal ecosystems? Evidence from modern and ancient ungulates. Palaeogeogr. Palaeoclimatol. Palaeoecol. 266, 69–82 (2008).Article 

Google Scholar 
Koch, P. L. Isotopic study of the biology of modern and fossil vertebrates. In Stable Isotopes in Ecology and Environmental Science (eds Michener, R. & Lajtha, K.) 99–154 (Blackwell Publishing Ltd, 2007). https://doi.org/10.1002/9780470691854.ch5.Chapter 

Google Scholar 
Hofman-Kamińska, E. et al. Foraging habitats and niche partitioning of European large herbivores during the holocene—Insights from 3D dental microwear texture analysis. Palaeogeogr. Palaeoclimatol. Palaeoecol. 506, 183–195 (2018).Article 

Google Scholar 
Bocherens, H., Hofman-Kamińska, E., Drucker, D. G., Schmölcke, U. & Kowalczyk, R. European bison as a refugee species? Evidence from isotopic data on Early Holocene bison and other large herbivores in northern Europe. PLoS ONE 10, 1–19 (2015).Article 
CAS 

Google Scholar 
Hu, Y. et al. Earliest evidence for commensal processes of cat domestication. Proc. Natl. Acad. Sci. USA. 111, 116–120 (2014).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Haruda, A. F. et al. The earliest domestic cat on the Silk Road. Sci. Rep. 10, 11241 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Meckstroth, A. M., Miles, A. K. & Chandra, S. Diets of introduced predators using stable isotopes and stomach contents. J. Wildl. Manag. 71, 2387–2392 (2007).Article 

Google Scholar 
McDonald, B. W. et al. High variability within pet foods prevents the identification of native species in pet cats’ diets using isotopic evaluation. PeerJ 8, e8337 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Maeda, T., Nakashita, R., Shionosaki, K., Yamada, F. & Watari, Y. Predation on endangered species by human-subsidized domestic cats on Tokunoshima Island. Sci. Rep. 9, 16200 (2019).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Stewart, G. R., Aidar, M. P. M., Joly, C. A. & Schmidt, S. Impact of point source pollution on nitrogen isotope signatures (δ15N) of vegetation in SE Brazil. Oecologia 131, 468–472 (2002).ADS 
PubMed 
Article 

Google Scholar 
Graven, H., Keeling, R. F. & Rogelj, J. Changes to carbon isotopes in atmospheric CO2 over the industrial era and into the future. Glob. Biogeochem. Cycles 34, 1–21 (2020).Article 
CAS 

Google Scholar 
DeNiro, M. J. Postmortem preservation and alteration of in vivo bone collagen isotope ratios in relation to palaeodietary reconstruction. Nature 317, 806–809 (1985).ADS 
CAS 
Article 

Google Scholar 
Linderholm, A. & Kjellström, A. Stable isotope analysis of a medieval skeletal sample indicative of systemic disease from Sigtuna Sweden. J. Archaeol. Sci. 38, 925–933 (2011).Article 

Google Scholar 
Webb, E. C. et al. Compound-specific amino acid isotopic proxies for distinguishing between terrestrial and aquatic resource consumption. Archaeol. Anthropol. Sci. 10, 1–18 (2018).Article 

Google Scholar 
Müldner, G. & Richards, M. P. Stable isotope evidence for 1500 years of human diet at the city of York, UK. Am. J. Phys. Anthropol. 133, 682–697 (2007).PubMed 
Article 

Google Scholar 
Müldner, G. & Richards, M. P. Fast or feast: Reconstructing diet in later medieval England by stable isotope analysis. J. Archaeol. Sci. 32, 39–48 (2005).Article 

Google Scholar 
van der Sluis, L. G., Hollund, H. I., Kars, H., Sandvik, P. U. & Denham, S. D. A palaeodietary investigation of a multi-period churchyard in Stavanger, Norway, using stable isotope analysis (C, N, H, S) on bone collagen. J. Archaeol. Sci. Rep. 9, 120–133 (2016).
Google Scholar 
Polet, C. & Katzenberg, M. A. Reconstruction of the diet in a mediaeval monastic community from the coast of Belgium. J. Archaeol. Sci. 30, 525–533 (2003).Article 

Google Scholar 
Kosiba, S. B., Tykot, R. H. & Carlsson, D. Stable isotopes as indicators of change in the food procurement and food preference of Viking Age and Early Christian populations on Gotland (Sweden). J. Anthropol. Archaeol. 26, 394–411 (2007).Article 

Google Scholar 
Olsen, K. C. et al. Isotopic anthropology of rural German medieval diet: Intra- and inter-population variability. Archaeol. Anthropol. Sci. 10, 1053–1065 (2018).Article 

Google Scholar 
Benevolo, L. The European City (Blackwell Publishers, 1993).
Google Scholar 
Barrett, J. et al. Detecting the medieval cod trade: A new method and first results. J. Archaeol. Sci. 35, 850–861 (2008).Article 

Google Scholar 
Barrett, J. H. et al. Interpreting the expansion of sea fishing in medieval Europe using stable isotope analysis of archaeological cod bones. J. Archaeol. Sci. 38, 1516–1524 (2011).Article 

Google Scholar 
Bogaard, A., Heaton, T. H. E., Poulton, P. & Merbach, I. The impact of manuring on nitrogen isotope ratios in cereals: Archaeological implications for reconstruction of diet and crop management practices. J. Archaeol. Sci. 34, 335–343 (2007).Article 

Google Scholar 
Heaton, T. H. E. Spatial, species, and temporal variations in the 13C/12C ratios of C3 plants: Implications for palaeodiet studies. J. Archaeol. Sci. 26, 637–649 (1999).Article 

Google Scholar 
Bogaard, A. et al. Crop manuring and intensive land management by Europe’s first farmers. Proc. Natl. Acad. Sci. USA. 110, 12589–12594 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Styring, A. K. et al. Refining human palaeodietary reconstruction using amino acid δ15N values of plants, animals and humans. J. Archaeol. Sci. 53, 504–515 (2015).CAS 
Article 

Google Scholar 
Guiry, E. Complexities of stable carbon and nitrogen isotope biogeochemistry in ancient freshwater ecosystems: Implications for the study of past subsistence and environmental change. Front. Ecol. Evol. 7, 313 (2019).Article 

Google Scholar 
Fuller, B. T., Müldner, G., Van Neer, W., Ervynck, A. & Richards, M. P. Carbon and nitrogen stable isotope ratio analysis of freshwater, brackish and marine fish from Belgian archaeological sites (1st and 2nd millennium AD). J. Anal. At. Spectrom. 27, 807–820 (2012).CAS 
Article 

Google Scholar 
Robson, H. K. et al. Carbon and nitrogen stable isotope values in freshwater, brackish and marine fish bone collagen from Mesolithic and Neolithic sites in central and northern Europe. Environ. Archaeol. 21, 105–118 (2016).Article 

Google Scholar 
Hobson, K. A., Piatt, J. F. & Pitocchelli, J. Using stable isotopes to determine seabird trophic relationships. J. Anim. Ecol. 63, 786–798 (1994).Article 

Google Scholar 
Guiry, E. & Buckley, M. Urban rats have less variable, higher protein diets. Proc. R. Soc. B Biol. Sci. 285, 20181441 (2018).Article 
CAS 

Google Scholar 
Bicknell, A. W. J. et al. Stable isotopes reveal the importance of seabirds and marine foods in the diet of St Kilda field mice. Sci. Rep. 10, 1–12 (2020).Article 
CAS 

Google Scholar 
Hoffmann, R. C. Medieval fishing. In Working with Water in Medieval Europe. Technology and Resource-Use (ed. Squatriti, P.) 331–393 (Brill, 2000).
Google Scholar 
Gillies, C. & Clout, M. The prey of domestic cats (Felis catus) in two suburbs of Auckland City, New Zealand. J. Zool. 259, 309–315 (2003).Article 

Google Scholar 
Brickner-Braun, I., Geffen, E. & Yom-Tov, Y. The domestic cat as a predator of Israeli wildlife. Isr. J. Ecol. Evol. 53, 129–142 (2007).Article 

Google Scholar 
Flockhart, D. T. T., Norris, D. R. & Coe, J. B. Predicting free-roaming cat population densities in urban areas. Anim. Conserv. 19, 472–483 (2016).Article 

Google Scholar 
Castañeda, I., Zarzoso-Lacoste, D. & Bonnaud, E. Feeding behaviour of red fox and domestic cat populations in suburban areas in the south of Paris. Urban Ecosyst. 23, 731–743 (2020).Article 

Google Scholar 
Zhu, Y., Siegwolf, R. T. W., Durka, W. & Körner, C. Phylogenetically balanced evidence for structural and carbon isotope responses in plants along elevational gradients. Oecologia 162, 853–863 (2010).ADS 
PubMed 
Article 

Google Scholar 
Männel, T. T., Auerswald, K. & Schnyder, H. Altitudinal gradients of grassland carbon and nitrogen isotope composition are recorded in the hair of grazers. Glob. Ecol. Biogeogr. 16, 583–592 (2007).Article 

Google Scholar 
Pińska, K. & Badura, M. Warunki przyrodnicze i dieta roślinna mieszkańców Pucka w późnym średniowieczu. In Puck – kultura materialna małego miasta w późnym średniowieczu (ed. Starski, M.) 517 (Uniwersytet Warszawski, 2017).
Google Scholar 
Lefebvre, A. et al. Morphology of estuarine bedforms, Weser Estuary, Germany. Earth Surf. Process. Landforms 47, 242–256 (2022).ADS 
Article 

Google Scholar 
Bischop, D. & Von der Küchelmann, H. C. Küche in den Graben – Bremens Stadtgraben und die Essgewohnheiten seiner Anwohner an der Wende zur Frühen Neuzeit. In Lebensmittel im Mittelalter und in der frühen Neuzeit. Erzeugung, Verarbeitung, Versorgung. Beiträge des 16. Kolloquiums des Arbeitskreises zur archäologischen Erforschung des mittelalterlichen Handwerks, Soester Beiträge zur Archäologie 15 (ed. Melzer, W.) 137–151 (Mocker und Jahn, 2018).
Google Scholar 
Elmshäuser, K. & Pordzik, V. V. Lachsgarnen, Tomen und Kumpanen – Die älteste Bremer Fischeramtsrolle. Bremisches Jahrb. 98, 13–72 (2019).
Google Scholar 
Küchelmann, H. C. Viel Butter bei wenig Fisch. Zwei Fischknochenkomplexe des 12.–13. Jahrhunderts aus der Bremer Altstadt. In Grenzen überwinden. Archäologie zwischen Disziplin und Disziplinen. Festschrift für Uta Halle zum 65. Geburtstag, Internationale Archäologie Studia Honoraria 40 (eds Kahlow, S. et al.) 413–426 (Verlag Marie Leidorf GmbH, 2021).
Google Scholar 
Schwarcz, H. P. & Schoeninger, M. J. Stable isotope analyses in human nutritional ecology. Am. J. Phys. Anthropol. 34, 283–321 (1991).Article 

Google Scholar 
Wallace, M. et al. Stable carbon isotope analysis as a direct means of inferring crop water status and water management practices. World Archaeol. 45, 388–409 (2013).Article 

Google Scholar 
van der Merwe, N. J. & Medina, E. The canopy effect, carbon isotope ratios and foodwebs in amazonia. J. Archaeol. Sci. 18, 249–259 (1991).Article 

Google Scholar 
Ervynck, A. Orant, pugnant, laborant. The diet of the three orders in the feudal society of medieval north-western Europe. In Behaviour Behind Bones. The Zooarchaeology of Ritual, Religion, Status and Identity (eds O’Day, S. J. et al.) 215–223 (Oxbow Books, 2004).
Google Scholar 
von den Driesch, A. A guide to the measurement of animal bones from archaeological sites. Peabody Museum Bull. 1, 1–137 (1976).
Google Scholar 
O’Connor, T. P. Wild or domestic? Biometric variation in the cat Felis silvestris Schreber. Int. J. Osteoarchaeol. 17, 581–595 (2007).Article 

Google Scholar 
Kratochvíl, Z. Schadelkriterien der Wild- und Hauskatze (Felis silvestris silvestris Schreber 1777 und Felis s. f. catus L. 1758). Acta Sci. Nat. Brno 7, 1–50 (1973).
Google Scholar 
Kratochvíl, Z. Das Postkranialskelett der Wild- und Hauskatze (Felis silvestris und F. lybica f. catus). Acta Sci. Nat. Brno 10, 1–43 (1976).
Google Scholar 
Dyce, K. M., Sack, W. O. & Wensing, C. J. G. Textbook of Veterinary Anatomy (Saunders/Elsevier, 2010).
Google Scholar 
Krajcarz, M. et al. On the trail of the oldest domestic cat in Poland. An insight from morphometry, ancient DNA and radiocarbon dating. Int. J. Osteoarchaeol. 26, 912–919 (2016).Article 

Google Scholar 
Bronk Ramsey, C. Radiocarbon calibration and analysis of stratigraphy: The OxCal program. Radiocarbon 37, 425–430 (1995).CAS 
Article 

Google Scholar 
Bronk Ramsey, C., Dee, M., Lee, S., Nakagawa, T. & Staff, R. Developments in the calibration and modeling of radiocarbon dates. Radiocarbon 52, 953–961 (2010).Article 

Google Scholar 
Ferreira, J. P., Leitão, I., Santos-Reis, M. & Revilla, E. Human-related factors regulate the spatial ecology of domestic cats in sensitive areas for conservation. PLoS ONE 6, e25970 (2011).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pirie, T. J., Thomas, R. L. & Fellowes, M. D. E. Pet cats (Felis catus) from urban boundaries use different habitats, have larger home ranges and kill more prey than cats from the suburbs. Landsc. Urban Plan. 220, 104338 (2022).Article 

Google Scholar 
Bocherens, H. et al. Paleobiological implications of the isotopic signatures (13C, 15N) of fossil mammal collagen in Scladina cave (Sclayn, Belgium). Quat. Res. 48, 370–380 (1997).Article 

Google Scholar 
Longin, R. New method of collagen extraction for radiocarbon dating. Nature 230, 241–242 (1971).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Boudin, M., Boeckx, P., Vandenabeele, P. & Van Strydonck, M. Improved radiocarbon dating of contaminated protein-containing archaeological samples via cross-flow nanofiltrated amino acids. Rapid Commun. Mass Spectrom. 27, 2039–2050 (2013).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Wojcieszak, M., Van Den Brande, T., Ligovich, G. & Boudin, M. Pretreatment protocols performed at the Royal Institute for Cultural Heritage (RICH) prior to AMS 14C measurements. Radiocarbon 62, e14–e24 (2020).Article 

Google Scholar 
Hammer, Ø. PAST. PAleontological Statistics. Version 4.05 Reference manual (Natural History Museum University of Oslo, 2021).
Google Scholar 
Hammer, Ø., Harper, D. A. T. & Ryan, P. D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4, 1–9 (2001).
Google Scholar 
Rohland, N., Glocke, I., Aximu-Petri, A. & Meyer, M. Extraction of highly degraded DNA from ancient bones, teeth and sediments for high-throughput sequencing. Nat. Protoc. 13, 2447–2461 (2018).CAS 
PubMed 
Article 

Google Scholar 
Nguyen, L. T., Schmidt, H. A., Von Haeseler, A. & Minh, B. Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274 (2015).CAS 
PubMed 
Article 

Google Scholar  More