Ecology

Contents of Cd and Ca in Panax notogensing rootsThe Ca content of P. notoginseng roots increased significantly with the increase of lime application rates under the same concentration of oxalic acid sprayed on leaves (Table 2). Compared with no lime application, the Ca content was the highest increased by 212% under 3750 kg hm−2 lime without spraying oxalic acid. The content of Ca slightly increased with the increase of oxalic acid spraying concentrations under the same rate of lime application.Table 2 Effects of foliar spraying of oxalic acid on contents of Cd and Ca in roots of Panax notoginseng under Cd stress.Full size tableThe contents of Cd in roots ranged from 0.22 to 0.70 mg kg−1. The content of 2250 kg hm−2 Cd decreased greatly with the increase of lime application rates under the same spraying concentration of oxalic acid. Compared with the control, the root Cd contents decreased by 68.57% under the application of 2250 kg hm−2 lime and 0.1 mol L−1 oxalic acid spraying. The Cd contents of P. notoginseng roots decreased significantly with the increase of oxalic acid spraying concentrations under application of non-lime and 750 kg hm−2 lime. The root Cd contents decreased at first and then increased with the increase of oxalic acid concentrations under the application of 2250 kg hm−2 lime and 3750 kg hm−2 lime. In addition, the Bivariate analysis showed that the Ca content of P. notoginseng roots was significantly affected by lime (F = 82.84**), and the Cd content of P. notoginseng roots was significantly affected by lime (F = 74.99**) and oxalic acid (F = 7.72*).MDA contents and relative antioxidase activitiesThe content of MDA decreased greatly with the increase of the rates of lime application and oxalic acid spraying concentrations. There was no significant difference in the content of MDA in the roots of P. notoginseng with non-lime and 3750 kg hm−2 lime application. Under 750 kg hm−2, 2250 kg hm−2 lime application, the MDA content with 0.2 mol L−1 oxalic acid spraying concentration treatment decreased by 58.38% and 40.21% comparing with non-oxalic acid spraying application, respectively. The content of MDA (7.57 nmol g−1) was the lowest under 750 kg hm−2 lime application and 0.2 mol L−1 oxalic acid spraying treatment (Fig. 1).Figure 1Effects of foliar spraying of oxalic acid on contents of malondialdehyde in roots of Panax notoginseng under Cd stress. Notes The figure legend showed the spray concentration of oxalic acid (mol L−1), different lowercase letters indicate significant differences between treatments at the same lime application rate (P  Rb1  > R1. The contents of the three saponins had no significant difference with increase of the concentrations of oxalic acid spraying and no application of lime (Table 4).Table 4 Effects of foliar oxalate application on the percentages of three saponins in roots of Panax notoginseng under Cd stress.Full size tableThe contents of R1 with 0.2 mol L−1 oxalic acid spraying was significantly lower than that without oxalic acid spraying and rates of 750 or 3750 kg hm−2 lime application. Under the concentration of 0 or 0.1 mol L−1 oxalic acid spraying, there was no significant difference in contents of R1 with increase of rates of lime application. Under the concentration of 0.2 mol L−1 oxalic acid spraying, the contents of R1 with 3750 kg hm−2 lime was significantly lower 43.84% than that without lime application (Table 4).The contents of Rg1 increased at first and then decreased with the increase of oxalic acid spraying concentrations and 750 kg hm−2 lime application. Under the application rates of 2250 or 3750 kg hm−2 lime, the contents of Rg1 decreased with the increase of oxalic acid spraying concentration. With the same concentration of oxalic acid spraying, the Rg1 content increased at first and then decreased with the increase of lime application rates. Compared with the control, except that the Rg1 content with three concentrations of oxalic acid spraying and 750 kg hm−2 lime was higher than that of the control, the contents of Rg1 in the roots of P. notoginseng under other treatments was lower than that of the control. The Rg1 content was the highest with 750 kg hm−2 lime and 0.1 mol L−1 oxalic acid spraying treatment, which was higher 11.54% than that of the control (Table 4).The contents of Rb1 increased first and then decreased with the increase of oxalic acid spraying concentration and 2250 kg hm−2 lime application. The content of Rb1 with 0.1 mol L−1 oxalic acid spraying reached the maximum value of 3.46%, which was higher 74.75% than that without oxalic acid spraying treatment. Under other lime application treatments, there was no significant difference among different oxalic acid spraying concentrations. With 0.1 and 0.2 mol L−1 oxalic acid spraying treatments, the contents of Rb1 decreased at first and then decreased with the increase of lime application rates (Table 4).Contents of flavonoidsWith the same concentration of oxalic acid spraying, the content of flavonoids increased at first and then decreased with the increase of the amounts of lime application. There was no significant difference in the content of flavonoids under different concentrations of oxalic acid spraying without the application of lime or 3750 kg hm−2 lime. Under 750 and 2250 kg hm−2 lime application, the content of flavonoids increased at first and then decreased with the increase of the concentration of oxalic acid spraying. Under the treatment of 750 kg hm−2 application and 0.1 mol L−1 oxalic acid spraying, the content of flavonoids was the highest, which was 4.38 mg g−1, which was higher 18.38% than that of the same rate of lime application and without spraying oxalic acid. The content of flavonoids with 0.1 mol L−1 oxalic acid spraying treatment increased by 21.74% compared with that without oxalic acid spraying treatment and 2250 kg hm−2 lime application (Fig. 5).Figure 5Effects of foliar spraying of oxalate on the contents of flavonoids in roots of Panax notoginseng under Cd stress.Full size imageBivariate analysis showed that the content of soluble sugar in P. notoginseng root was significantly relationship with the amount of lime application and the concentration of oxalic acid spraying. The content of soluble protein in root was significantly relationship with lime application rates, both of lime and oxalic acid. The contents of free amino acid and proline in roots were significantly relationship with lime application rates, oxalic acid spraying concentrations, both of lime and oxalic acid (Table 5).Table 5 Variance analysis of the effects of oxalic acid, calcium and cadmium on the contents of multiple medicinal ingredients in the roots of Panax notoginseng (F value).Full size tableThe content of R1 in the root of P. notoginseng was significantly relationship with oxalic acid spraying concentrations, lime application rates, both of lime and oxalic acid. The content of flavonoids was significantly relationship with oxalic acid spraying concentrations, lime application rates. More

Pachauri, R. K. et al. Climate change 2014: synthesis report. Contribution of Working Groups I, II and III to the fifth assessment report of the Intergovernmental Panel on Climate Change. (IPCC, 2014).Thackeray, S. J. et al. Phenological sensitivity to climate across taxa and trophic levels. Nature 535, 241–245 (2016).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Walther, G.-R. et al. Ecological responses to recent climate change. Nature 416, 389–395 (2002).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Root, T. L., MacMynowski, D. P., Mastrandrea, M. D. & Schneider, S. H. Human-modified temperatures induce species changes: Joint attribution. Proc. Natl. Acad. Sci. USA 102, 7465–7469 (2005).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Menzel, A. et al. European phenological response to climate change matches the warming pattern. Glob. Change Biol. 12, 1969–1976 (2006).ADS 
Article 

Google Scholar 
Rosenzweig, C. et al. Attributing physical and biological impacts to anthropogenic climate change. Nature 453, 353–357 (2008).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Myneni, R. B., Keeling, C. D., Tucker, C. J., Asrar, G. & Nemani, R. R. Increased plant growth in the northern high latitudes from 1981 to 1991. Nature 386, 698–702 (1997).ADS 
CAS 
Article 

Google Scholar 
Walkovszky, A. Changes in phenology of the locust tree (Robinia pseudoacacia L.) in Hungary. Int. J. Biometeorol. 41, 155–160 (1998).ADS 
Article 

Google Scholar 
Crick, H. Q. P. & Sparks, T. H. Climate change related to egg-laying trends [8]. Nature 399, 423–424 (1999).ADS 
CAS 
Article 

Google Scholar 
Both, C., Van Asch, M., Bijlsma, R. G., Van Den Burg, A. B. & Visser, M. E. Climate change and unequal phenological changes across four trophic levels: Constraints or adaptations?. J. Anim. Ecol. 78, 73–83 (2009).PubMed 
Article 

Google Scholar 
Brown, J. L., Li, S. H. & Bhagabati, N. Long-term trend toward earlier breeding in an American bird: A response to global warming?. Proc. Natl. Acad. Sci. USA. 96, 5565–5569 (1999).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Miller-Rushing, A. J., Lloyd-Evans, T. L., Primack, R. B. & Satzinger, P. Bird migration times, climate change, and changing population sizes. Glob. Chang. Biol. 14, 1959–1972 (2008).ADS 
Article 

Google Scholar 
Hughes, L. Biological consequences of global warming: Is the signal already apparent?. Trends Ecol. Evol. 15, 56–61 (2000).CAS 
PubMed 
Article 

Google Scholar 
Kudo, G., Nishikawa, Y., Kasagi, T. & Kosuge, S. Does seed production of spring ephemerals decrease when spring comes early?. Ecol. Res. 19, 255–259 (2004).Article 

Google Scholar 
Doi, H., Gordo, O. & Katano, I. Heterogeneous intra-annual climatic changes drive different phenological responses at two trophic levels. Clim. Res. 36, 181–190 (2008).Article 

Google Scholar 
Both, C., Bouwhuis, S., Lessells, C. M. & Visser, M. E. Climate change and population declines in a long-distance migratory bird. Nature 441, 81–83 (2006).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Post, E. & Forchhammer, M. C. Climate change reduces reproductive success of an Arctic herbivore through trophic mismatch. Philos. Trans. R. Soc. B Biol. Sci. 363, 2369–2375 (2008).Article 

Google Scholar 
Thackeray, S. J. et al. Trophic level asynchrony in rates of phenological change for marine, freshwater and terrestrial environments. Glob. Chang. Biol. 16, 3304–3313 (2010).ADS 
Article 

Google Scholar 
Visser, M. E. & Both, C. Shifts in phenology due to global climate change: The need for a yardstick. Proc. R. Soc. B Biol. Sci. 272, 2561–2569 (2005).Article 

Google Scholar 
Hegland, S. J., Nielsen, A., Lázaro, A., Bjerknes, A. L. & Totland, Ø. How does climate warming affect plant-pollinator interactions?. Ecol. Lett. 12, 184–195 (2009).PubMed 
Article 

Google Scholar 
Brown, C. J. et al. Quantitative approaches in climate change ecology. Glob. Chang. Biol. 17, 3697–3713 (2011).ADS 
PubMed Central 
Article 

Google Scholar 
Parmesan, C., Duarte, C., Poloczanska, E., Richardson, A. J. & Singer, M. C. Overstretching attribution. Nat. Clim. Chang. 1, 2–4 (2011).ADS 
Article 

Google Scholar 
Damos, P. & Savopoulou-Soultani, M. Temperature-driven models for insect development and vital thermal requirements. Psyche (London) 2012, (2012).Osawa, T. et al. Climate-mediated population dynamics enhance distribution range expansion in a rice pest insect. Basic Appl. Ecol. 30, 41–51 (2018).Article 

Google Scholar 
Kiritani, K. Predicting impacts of global warming on population dynamics and distribution of arthropods in Japan. Popul. Ecol. 48, 5–12 (2006).Article 

Google Scholar 
Ozawa, A., Uchiyama, T. & Kasai, A. Estimating the day on which the numbers of adult tea spiny whiteflies (Aleurocanthus camelliae Kanmiya and Kasai) peaked and the number of generations produced in a year based on the effective cumulative temperature in tea fields. Annu. Rep. Kansai Plant Prot. Soc. 58, 57–64 (2016) (in Japanese).Article 

Google Scholar 
Ozawa, A., Saito, T. & Ikeda, F. Effect of host plant and temperature on reproduction of American Serpentine Leafminer, Liriomyza trifolii (Burgess). Japan. J. Appl. Entomol. Zool. 43, 41–48 (1999) (in Japanese).Article 

Google Scholar 
Baier, P., Pennerstorfer, J. & Schopf, A. PHENIPS—A comprehensive phenology model of Ips typographus (L.) (Col., Scolytinae) as a tool for hazard rating of bark beetle infestation. For. Ecol. Manage. 249, 171–186 (2007).Article 

Google Scholar 
Karuppaiah, V. & Sujayanad, G. K. Impact of climate change on population dynamics of insect pests. (2012).Lipper, L. et al. Climate-smart agriculture for food security. Nat. Clim. Chang. 4, 1068–1072 (2014).ADS 
Article 

Google Scholar 
Campbell, B. M. et al. Reducing risks to food security from climate change. Glob. Food Sec. 11, 34–43 (2016).Article 

Google Scholar 
Kiritani, K. The low development threshold temperature and the thermal constant in insects and mites in Japan (2nd edition). Bull. Nationai Instirute Agro-Environmental Sci. 31, 1–74 (2012) (in Japanese).
Google Scholar 
Nagasawa, A., Takahashi, A. & Higuchi, H. Host plant use for oviposition by Trigonotylus caelestialium (Hemiptera: Miridae) and Stenotus rubrovittatus (Hemiptera: Miridae). Appl. Entomol. Zool. 47, 331–339 (2012).Article 

Google Scholar 
Higuchi, H. Ecology and management of rice bugs causing pecky rice. Japan. J. Appl. Entomol. Zool. 54, 171–188 (2010) (in Japanese).Article 

Google Scholar 
Ohtomo, R. Occurrence and control of Stenotus rubrovittatus (Hemiptera: Miridae) in Touhoku area in Japan. Japan. J. Appl. Entomol. Zool. 57, 137–149 (2013) (in Japanese).Article 

Google Scholar 
Kiritani, K. The impact of global warming and land-use change on the pest status of rice and fruit bugs (Heteroptera) in Japan. Glob. Chang. Biol. 13, 1586–1595 (2007).ADS 
Article 

Google Scholar 
Nagasawa, A. & Higuchi, H. Suitability of poaceous plants for nymphal growth of the pecky rice bugs Trigonotylus caelestialium and Stenotus rubrovittatus (Hemiptera: Miridae) in Niigata, Japan. Appl. Entomol. Zool. 47, 421–427 (2012).Article 

Google Scholar 
Gordo, O. & Sanz, J. J. Temporal trends in phenology of the honey bee Apis mellifera (L.) and the small white Pieris rapae (L.) in the Iberian Peninsula (1952–2004). Ecol. Entomol. 31, 261–268 (2006).Article 

Google Scholar 
Sparks, T. H. & Yates, T. J. The effect of spring temperature on the appearance dates of British butterflies 1883–1993. Ecography (Cop.) 20, 368–374 (1997).Article 

Google Scholar 
Stefanescu, C., Penuelas, J. & Filella, I. Effects of climatic change on the phenology of butterflies in the northwest Mediterranean Basin. Glob. Chang. Biol. 9, 1494–1506 (2003).ADS 
Article 

Google Scholar 
Dell, D., Sparks, T. H. & Dennis, R. L. H. Climate change and the effect of increasing spring temperatures on emergence dates of the butterfly Apatura iris (Lepidoptera: Nymphalidae). Eur. J. Entomol. 102, 161–167 (2005).Article 

Google Scholar 
Tabuchi, K. et al. Rice bugs in the Tohoku Region: Their occurrence and damage from 2003 to 2013. Bull. Natl. Agric. Res. Cent. Tohoku Reg. 117, 63–115 (2015) (in Japanese).
Google Scholar 
Jolly, W. M. & Running, S. W. Effects of precipitation and soil water potential on drought deciduous phenology in the Kalahari. Glob. Chang. Biol. 10, 303–308 (2004).ADS 
Article 

Google Scholar 
Kriticos, D. J. et al. CliMond: Global high-resolution historical and future scenario climate surfaces for bioclimatic modelling. Methods Ecol. Evol. 3, 53–64 (2012).Article 

Google Scholar 
Takeda, A. & Shimizu, K. Characteristics of the damaged grain caused by the sorghum plant bug, Stenotus rubrovittatus (Hemiptera: Miridae) at different infection periods. Annu. Rep. Kanto-Tosan Plant Prot. Soc. 2009, 85–87 (2009) (in Japanese).
Google Scholar 
Ishimoto, M. Seasonal Prevalence of Occurrence of the Rice Leaf Bug, Trigonotylus caelestialium (Kirkaldy) (Heteroptera: Miridae) on Paddy Rice Plants. Jpn. J. Appl. Entomol. Zool. 48, 79–85 (2004) (in Japanese).
Katase, M., Shimizu, K., Siina, S., Hagiwara, K. & Iwai, H. Seasonal occurrence of rice bugs in the northern part of Chiba Prefecture. Annu. Rep. Kanto-Tosan Plant Prot. Soc. 2007, 99–104 (2007) (in Japanese).
Google Scholar 
Shintani, Y. Effect of seasonal variation in host–plant quality on the rice leaf bug, Trigonotylus caelestialium. Entomol. Exp. Appl. 133, 128–135 (2009).Article 

Google Scholar 
Takada, M. B., Yoshioka, A., Takagi, S., Iwabuchi, S. & Washitani, I. Multiple spatial scale factors affecting mirid bug abundance and damage level in organic rice paddies. Biol. Control 60, 169–174 (2012).Article 

Google Scholar 
Yoshioka, A., Takada, M. B. & Washitani, I. Landscape effects of a non-native grass facilitate source populations of a native generalist bug, Stenotus rubrovittatus, in a heterogeneous agricultural landscape. J. Insect Sci. 14, 1–14 (2014).Article 

Google Scholar 
Watanabe, T. & Higuchi, H. Recent occurrence and problem of rice bugs. Plant Prot. 60, 201–203 (2006) (in Japanese).
Google Scholar 
Seino, H. An estimation of distribution of meteorological elements using GIS and AMeDAS data. J. Agric. Meteorol. 48, 379–383 (1993) (in Japanese).Article 

Google Scholar 
Ishigooka, Y., Tsuneo, K., Nishimori, M., Hasegawa, T. & Ohno, H. Spatial characterization of recent hot summers in Japan with agro-climatic indices related to rice production. J. Agric. Meteorol. 67, 209–224 (2011).Article 

Google Scholar 
Yamasaki, K. et al. Intraspecific variations in life history traits of two pecky rice bug species from Japan: Mapping emergence dates and number of annual generations. Ecol. Evol. (2021).Sakagami, Y. & Korenaga, R. Triangle method—A simple method for the estimation of total effective temperature. Japan. J. Appl. Entomol. Zool. 25, 52–54 (1981) (in Japanese).Article 

Google Scholar 
Shigehisa, S. Seasonal changes in egg diapause induction and effects of photoperiod and temperature on egg diapause in the sorghum plant bug, Stenotus rubrovittatus (Matsumura) (Heteroptera: Miridae). Japan. J. Appl. Entomol. Zool. 52, 229–232 (2008) (in Japanese).Article 

Google Scholar 
Ishimoto, M. Oviposition of sorghum plant bug, Stenotus rubrovittatus (Matsumura) (Heteroptera: Miridae) on rice plants. Japan. J. Appl. Entomol. Zool. 55, 193–197 (2011) (in Japanese).Article 

Google Scholar 
Ogata, M. et al. Fecundity and longevity in adults of the sorghum plant bug, Stenotus rubrovittatus (Matsumura) (Heteroptera: Miridae) under laboratory conditions. Japan. J. Entomol. 13, 129–132 (2010) (in Japanese).
Google Scholar 
Niiyama, T. & Iitomi, A. Optimal control timing of the rice leaf bug, Trigonotylus caelestialium (Heteroptera: Miridae). Annu. Rep. Soc. Plant Prot. North Japan. (in Japanese) (2003).Takeda, A., Shimizu, K., Shiina, S., Hagiwara, K. & Katase, M. Seasonal prevalence of Stenotus rubrovittatus (Hemiptera: Miridae) in a gramineous weed field and a rice field. Annu. Rep. Kanto-Tosan Plant Prot. Soc. 2008, 97–102 (2008) (in Japanese).
Google Scholar 
Niiyama, T. Studies on the ecology of Trigonotylus caelestialium and establishment of pesticide control techniques. Annu. Rep. Akita Prefect. Agric. Inst. 49, 147–180 (2009) (in Japanese).
Google Scholar 
Sato, M. Latest infestation period of Miridae causing pecky rice damage in Aomori Prefecture. Annu. Rep. Soc. Plant Prot. North Japan 2014, 129–134 (2014) (in Japanese).
Google Scholar 
Johnson, J. B. & Omland, K. S. Model selection in ecology and evolution. Trends Ecol. Evol. 19, 101–108 (2004).PubMed 
Article 

Google Scholar 
Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest package: tests in linear mixed effects models. J. Stat. Softw. 82, 1–26 (2017).Article 

Google Scholar  More

There are two phases involved in supporting the technical quality of the dataset: (i) the process used by Taronga when collecting data for each shark-bite incident, and (ii) the consistency modifications that we made during manuscript development.Phase oneFor each shark-bite incident, Taronga attempts to contact the most-relevant person involved in the event as possible — e.g., victim, victim’s family, witnesses. Those contacted are asked to complete a questionnaire with information relating to the shark bite. If applicable, questionnaires can also be completed by a fisheries officer in the relevant State and sent to Taronga. Taronga works closely with experts in each State’s fisheries department to validate information sourced in media reports. Each shark-bite case is unique, so the validation process varies depending on details specific to each incident. Forensic shark scientists within each State department are contacted after each shark bite to confirm details related to the incident. For example, validating species responsible for the bite often requires forensic analysis through expert examination of bite marks or artefacts21,22,23,24,25. When available, video footage is analysed for validation of information, such as confirmation of shark species and length.The database is cross-checked annually for the previous year with the International Shark-Attack File in Florida, USA, as well as with fisheries officers to ensure consistency. The database is cross-checked with the acknowledgement that there are discrepancies between versions due to differences in inclusion criteria. For example, the International Shark-Attack File includes bites where the victim was bitten aboard a boat, whereas the database we present here does not include bites aboard a boat.We acknowledge the limitations associated with this database, such as differences in reporting over time. For example, incidents might be reported more in recent decades due to technological advances making reporting more accessible or media publicising these events more widely18. There might also be reporting biases, for example, victims could be more likely to report a bite by a large, potentially dangerous shark (e.g., white, tiger, or bull shark) rather than a smaller, less-dangerous shark species (e.g., wobbegong shark). We also completed a quality assessment of the original database fields and redesigned the data acquisition and entry process (see Phase two) to allow exploration of shark-bite trends and patterns in Australia.Phase twoWe identified errors and inconsistencies in database fields. To avoid additional errors and inconsistencies, and to obtain a quality-controlled database, we redesigned the process for gaining and entering information into the database. This included creating a data descriptor (Supplementary File 2) used as a protocol to inform which questions to ask in the questionnaire. The data descriptor also directs the format of database entries by specifying information required in each field and by indicating the format of each entry (i.e., numeric, descriptive, or categorical). We manually inspected all previously entered data and adapted them to match the data descriptor. We checked each entry using the filter function in Microsoft Excel to identify any spelling and grammatical errors in the fields and ensure that all categories were grammatically identical. We standardised all metrics during this process (e.g., the data descriptor now stipulates that all length measurements should be recorded in metres).We validated and standardised the geographical locations of shark bites by converting all coordinates into decimal degrees using Microsoft Excel. We subsequently plotted all coordinates using the ggmap library27 in R (Version 4.0.2) (R Core Team 2020) (Fig. 2). We corrected any unusual coordinates (e.g., outside of Australian waters) and crosschecked them with site descriptions and states to ensure validity.Fig. 2Geographical locations of 1,196 shark bites in Australia. Each shark-bite incident is indicated by a red dot. (a) all shark-bite incidents; (b) bites most likely inflicted by bull sharks, (c) tiger sharks, and (d) white sharks. Two bite incidents that occurred at the Australian external territory, Cocos (Keeling) Islands, are not included on these maps. Background layers show elevation and major, perennial watercourses.Full size imageDuring the development of the data-descriptor protocol, we converted some previous descriptive columns into categorical columns. These columns included (but are not limited to) victim activity, attractant, injury location, injury severity, and weather condition. Converting these columns into categories facilitates analysis to investigate shark-bite patterns. For example, we converted victim activity into a categorical field to restrict answers to the following: snorkelling, motorised boating, unmotorised boating, boarding, swimming, standing, diving, fishing, or other, rather than allowing answers in any format. We used this information to create a time series to show the activity of shark-bite victims in Australia over time (Figs. 3 and 4). Shark bites have increased for boarders (including surfboarding, bodyboarding, kiteboarding, sailboarding, wakeboarding, and stand-up paddle boarding) over time, particularly since 1960 (Figs. 3b, 4). This is likely due to the increase in popularity of board sports, particularly surfing, since the 1960s28. This trend is likely not reflected in Fig. 3a because shark bites are unlikely to be classed as ‘provoked’ during board riding.Fig. 3Number of shark bites (black, dashed line) and proportion of activity done by the victim at the time of shark-bite incidents in Australia from 1900 to 2022. Panels represent shark bites in Australia that are; (a) provoked or (b) unprovoked. Boarding includes surfboarding, bodyboarding, kiteboarding, sailboarding, wakeboarding, and stand-up paddle boarding. Swimming includes snorkelling, spearfishing, freediving, body surfing, clinging to an object, falling into water, floating, or wading. Diving includes scuba-diving, hookah diving, or hard-hat diving. Fishing includes cleaning fish. No data for years 1908 and 1970.Full size imageFig. 4Number of shark bites (black, dashed line) and proportion of activity done by the victim at the time of shark-bite incidents in Australia from 1900 to 2022. Panels represent shark bites in Australia that are; (a) fatal or (b) non-fatal. Boarding includes surfboarding, bodyboarding, kiteboarding, sailboarding, wakeboarding, and stand-up paddle boarding. Swimming includes snorkelling, spearfishing, freediving, body surfing, clinging to an object, falling into water, floating, or wading. Diving includes scuba-diving, hookah diving, or hard-hat diving. Fishing includes cleaning fish. No data for years 1908 and 1970.Full size imageThe column representing a shark-bite victim’s recovery status also requires a categorical response, restricting answers to: fatal, injured, or uninjured (Fig. 5). Since 1900, the proportion of shark-bite-related fatalities has decreased (Fig. 5). This trend is also true for the three species most attributed to shark-bite-related fatalities, white (Carcharodon carcharias), tiger (Galeocerdo cuvier), and bull (Carcharhinus leucas) sharks. The decrease in shark-bite-related deaths is likely due to advancements in medical responses to shark-bite victims over time14 and better understanding among surfers about using tourniquets to stem bleeding following increased certification as first responders in workplace occupational health and safety requirements. Bites resulting in an uninjured victim includes interactions where the shark might have bitten the victim’s equipment (i.e., surfboard, bodyboard, kayak) rather than biting the person.Fig. 5Proportion of victim-recovery status (fatal = grey; red = injured; blue = uninjured) resulting from all unprovoked shark bites in Australia between 1900–2022. Blank years represent years without any reported occurrences. (a) all shark-bite incidents, (b) bites most likely inflicted by bull sharks, (c) tiger sharks, or (d) white sharks. No data for years 1908 and 1970.Full size imagePreviously, entries in the injury location column were descriptive. We converted the column to a categorical field restricting answers to the following: arm, hand, lower arm, upper arm, shoulder, neck, head, torso, leg, foot, calf, thigh, pelvic region, or other. During the analyses process, we further categorised these injury locations into four body areas (head, arm, torso, and leg) to assess how injury location affects recovery status (fatal or injured) (Fig. 6). Fatality most often occurred following shark bites to the torso (Fig. 6). This is likely due to the injuries to organs and major arteries resulting in blood loss, which is a leading cause of shark-bite fatalities29. This is the first time that the location of a shark bite on the body has been assessed relative to recovery status.Fig. 6Proportion of Australian shark bites resulting in either fatality or injury categorised by injury location on the victim’s body (left panel; 250 bites resulting in fatal injury, 723 bites resulting in non-fatal injury) and by species (right panel; bites by 201 tiger, 170 bull, 258 white, and 303 bites other sharks).Full size imageUnderstanding how the location of a shark-bite wound relates to victim recovery has value in informing the development of shark-bite mitigations. For example, the development of shark-bite-resistant wetsuits could potentially result in higher survival rates of the user if the fabric is concentrated around the torso region30,31. Redesigning data acquisition and entry process to allow for categorical columns permits these types of analyses.Some detail of a shark-bite incident might be lost by converting previously descriptive columns into categorical columns. We addressed this by complementing categorical columns with accompanying fields to retain details of the incident. For example, injury location and injury severity columns are both categorical and allow for user-friendly data analysis, whereas the injury description column is descriptive and provides added detail about the victim’s injuries if applicable. Individually, all three columns address certain aspects of the victim’s injuries, and together, all three columns comprehensively summarise injuries to the shark-bite victim.Our analysis of the Australian Shark-Incident Database suggests that tiger sharks are proportionally responsible for the most fatalities of all shark species in Australia (38% of all tiger shark bites result in fatality), followed by bull sharks (32% of all bull shark bites result in fatality), and white sharks (25% of all white shark bites result in fatality) (Fig. 6). We emphasise that these figures represent the overall percentage of bites resulting in fatality since 1791 and do not account for possible changes over time. These figures are also proportional to the number of bites by each respective species. In Australia, white sharks are responsible for the largest number of bites on humans (361 total) compared to tiger (229 total) and bull sharks (197 total). At the time of publication, white and tiger sharks were each responsible for 91 and 86 total fatalities on humans in Australia, respectively.There were 540 incidents in which time of day was recorded. We used these data to assess whether particular shark-bite incidents are more likely to occur at specific times of day (Fig. 7). To standardise reported 24-hour times, we took the location (latitude, longitude) information and reported time and date of each incident using the getSunlightTimes function in the suncalc library in R32 to calculate whether the incident occurred in one of four light-availability categories: dawn, day, dusk, or night. Shark bites occur mostly during the day, which likely reflects time of day when there are more ocean users present. However, there is a slightly higher proportion of bites at dusk for bull sharks compared to tiger and white sharks (Fig. 7). Identifying these trends can assist authorities in developing data-driven educational messaging as a shark-bite mitigation measure. This is important considering that enhanced education is the preferred mitigation measure scored by ocean users in New South Wales33.Fig. 7Period of day (corrected for local time) distribution of provoked and unprovoked shark bites in Australia by shark species from 1791 to 2022 (n = 540 bites from all species, 70 from bull, 59 from tiger, and 217 from white sharks).Full size imageThese examples demonstrate that the Australian Shark-Incident Database will be useful for scientists to analyse environmental, social, and biological related shark-bite patterns in Australia. Use of the newly developed data descriptor to standardise future applications and account for quality assurance and control will aid in keeping the database consistent for ease of analysis and interpretation. Ultimately, the publishing of this database will improve our understanding of shark-bite incidents in Australia and will equip us with the knowledge to aim to avoid or predict these events in the future. More

For this initial compilation, we focused on gathering traits from field guides and species accounts rather than the primary research literature because each represents the culmination of a comprehensive effort to describe a regional flora/fauna by local experts25. Authors of these guides have already done the hard work of scouring the literature, corresponding with fellow naturalists, and compiling occurrence records to support range, phenology, and habitat associations26. We began by performing a comprehensive review of all the holdings in the Florida Museum of Natural History’s McGuire Center for Lepidoptera and Biodiversity library, at the University of Florida. This, and subsequent searches in online databases, allowed us to compile a list of references that currently has more than 800 relevant resources.We initially identified the categories of trait information available in each resource and its format to target volumes for trait extraction and processing. Given the unequal availability of resources among regions, we had the explicit goal of identifying a corpus that would maximize the number of extractable trait data from as many butterfly species as evenly across the globe as possible. This led to our choice of 117 volumes within several global regions (Fig. 2, Supplementary Material S1) and a focus on measurements (wingspan/forewing length), phenology (months of adult flight and total duration of flight in months) and voltinism (the number of adult flight periods per year), habitat affinities, and host plants as traits (Table 1, Supplementary Material S2).Table 1 The total number of species represented by each trait in LepTraits 1.0.Full size tableTo process these resources, we developed a protocol to scan each volume, extract verbatim natural language descriptions, provide quality control for extraction, and then resolve given taxonomic names to a standardized list27. This provided a database of trait information in which each “cell” included all text from a single resource relevant to one trait category of a single taxon. In order to “atomize” the raw text into standardized metrics or a controlled list of descriptive terms, we developed a methodology appropriate to each trait. This resulted in a more fine-grained dataset in which each “cell” included a single, standardized trait value. Since the values of these taxon-specific traits frequently differed among resources, we then calculated “consensus” traits for each species, for example, the average forewing length (Table 1). A graphical representation of this process with an example trait is illustrated in Fig. 1.Fig. 1A graphical illustration of the processing workflow used to compile, scan, digitize, extract, atomize, and compile species trait records from literature resources. (1) Literature resources were examined for potential trait data and compiled into a single library; (2) each literature resource was scanned into.pdf format so that text could be readily copy and pasted from species accounts; (3) each.pdf file was uploaded to an online database with associated metadata for each literature resource; (4) trait extractors utilized an online interface to extract verbatim, raw text from designated resources; (5) verbatim, raw text extracts were either automatically (via regular-expressions and keyword searches) or manually atomized to a controlled vocabulary; (6) species consensus traits were calculated by aggregating resource-level records by name-normalized taxonomy. Rulesets were used for consensus trait building and are detailed in the supplementary material. Both resource-level and species consensus traits are presented in the dataset.Full size imageResource compilation and ingestionText sources from the master list were digitized by multiple participating institutions. They scanned each page of the book and converted the images to editable text with Abbyy FineReader optical character recognition (OCR) software (abbyy.com). These PDFs with copy-and-pastable text were then uploaded to a secure, online database that included citation information about each resource. The geographic breadth covered by each resource was designated using the World Geographic Scheme (WGS)28; this information was used to assess geographic evenness of our trait compilation efforts. Resource metadata, including the WGS scheme, were kept with each resource in an online database where individuals could access scanned copies of the resource for trait extraction.Verbatim data extractionIndividual workers were assigned to a resource and instructed to copy verbatim trait information from the original source. They then pasted that text into the relevant data field in a standardized, electronic form on an online portal designed to facilitate extraction and processing. Most field guides and other book-length resources are organized within a taxonomic hierarchy to describe traits of a family with a contiguous block of text, for example, family, then genus, species, and finally subspecies within species. We call these text blocks describing a single taxon “accounts” (e.g., family account, species account), and we recorded data at the taxonomic resolution provided in the original source. These taxonomic ranks included family, subfamily, tribe, genus, species, and subspecies. When information for a taxon was encountered outside its own account, the “extractor” (project personnel trained to manually extract verbatim text) assigned to glean data from the book entered this text into a separate entry for the taxon. Trait information from figure captions and tables were also extracted from the resource. Graphical representations of phenology and voltinism were common, and these visual data were converted to text descriptions. Each resource was extracted in stages, and each stage was subjected to a quality assurance and control process (see Technical Validation). This process corrected mistakes and attempted to find unextracted data overlooked by the extractor. These problems were corrected before the extractor could proceed with further trait extraction from the resource and were also used for training purposes.AtomizationVerbatim text extracts were subjected to an “atomization” process in which raw text was standardized into disaggregated, readily computable data. This conversion into the final trait data format (numerical, categorical, etc.) was two-pronged and involved both manual editing and semi-automated atomization of verbatim text. Regular expressions were used for most semi-automated atomization, including extraction of wing measurements, which were converted into centimeters. Keyword searches were also performed in the semi-automated pipeline for phenology, voltinism, and oviposition traits. For example, “univoltine” or “uni*” was searched for across the voltinism raw text, along with other search terms. All semi-automated atomization outputs were subject to quality assurance and control detailed further in Technical Validation. Manual atomization tasks were performed by multiple team members for traits which presented higher complexity. For example, habitat affinities and host plant associations were atomized manually along with a quality control protocol based on predefined rule sets that are described further in the Supplementary Material S3.Normalization and consensus traitsTo provide consensus traits at the species (and sometimes genus) level, we standardized nomenclature through a process we called “name-normalization,” which harmonizes taxonomy across all of our resources29. This name-normalization procedure relied on a comprehensive catalog of valid names and synonyms27. Following taxonomic harmonization, we compiled consensus traits based on rule sets specified in the metadata of each trait. For example, species-level consensus of primary and secondary host plant families required that at least one-third of the records for a given taxon list a particular family of plants (when multiple records were available).Categorical traits such as voltinism list all known voltinism patterns for a species regardless of geographic context. To this end, it is important that users of these data are aware that not all traits may be applicable to their study region. For example, some species may be univoltine at higher latitudes or elevations, but bivoltine elsewhere. We therefore present both the resource-level records as well as the species consensus traits for use in analysis.For this initial synopsis of butterfly species traits, we extracted records from 117 literature/web-based resources, resulting in 75,103 individual trait extraction records across 12,448 unique species, out of the ca. 19,200 species described to date27. Figure 2 indicates the geographic regions covered by our 117 resources, mapped at the resolution level-two regions in the World Geographic Scheme28. A full list of resources can be found in the Supplemental Material S1 as a bibliography. Similarly, the geographic distribution of trait records is indicated in Fig. 3. Resource and consensus species trait records varied in number and in the scope of taxonomic coverage. Table 1 indicates the number of unique records and species level records for each trait. Table 2 indicates the number of species-level records by family. Measurement traits, including wingspan and forewing length, were the most comprehensive traits extracted from our resource set. This represents one of the largest trait datasets and the most comprehensive dataset for butterflies to date.Fig. 2Geographic breadth of our butterfly trait resources. Using a global map of level-two regions (World Geographic Scheme, Brummitt 2001), we have indicated the total number of resources available within each geographic area). Grey areas indicate that no resources were extracted from that region.Full size imageFig. 3Geographic breadth of our butterfly trait records. Using a global map of level-two regions(World Geographic Scheme, Brummitt 2001), we have indicated the total number of trait records from each geographic region). Grey areas indicate that trait records were not extracted from that region.Full size imageTable 2 The number of species represented within each family in LepTraits 1.0.Full size table More

Model descriptionThe system dynamics model divided an urban dog population into the following subpopulations: (i) free-roaming dogs (both owned and unowned free-roaming, i.e. unrestricted dogs found on streets), (ii) shelter dogs (unowned restricted dogs living in shelters), and (iii) owned dogs (owned home-dwelling restricted dogs) (Fig. 1). The subpopulations change in size by individuals flowing between the different subpopulations or from flows extrinsically modelled (i.e. flows from subpopulations not included in the systems model; the acquisition of dogs from breeders and friends to the owned dog population, and the immigration/emigration of dogs from other neighbourhoods).Ordinary differential equations were used to describe the dog population dynamics. The models were written in R version 3.6.128, and numerically solved using the Runge–Kutta fourth order integration scheme with a 0.01 step sizes using the package “deSolve”29,30. For the baseline model, Eqs. (1–3) were used to describe the rates of change of dog subpopulations in the absence of management.Baseline free-roaming dog population (S):$$frac{dS}{dt}={r}_{s}times Stimes left(1-frac{S}{{K}_{s}}right)+alpha times O-delta times S$$
(1)
In the baseline model, the free-roaming dog population (Eq. 1) increases through the free-roaming dog intrinsic growth rate (rs), and the abandonment and roaming of dogs from the owned dog population (α) and decreases through adoption to the owned dog population (δ). The intrinsic growth rate is the sum of the effects of births, deaths, immigration, and emigration, which are not modelled separately. In this model, the growth rate of the free-roaming dog population is reduced depending on the population size in relation to the carrying capacity, through the logistic equation (rreal = rmax(1 − S/Ks))31. In the baseline simulation, the free-roaming dog population rises over time, until it stabilises at an equilibrium size.Baseline shelter dog population (H):$$frac{dH}{dt}=gamma times O-beta times H- {mu }_{h}times H$$
(2)
The shelter dog population (Eq. 2) increases through relinquishment of owned dogs (γ) and decreases through the adoption of shelter dogs to the owned dog population (β), and through the shelter dog death rate (µh). There is no carrying capacity for the shelter dog population as we assumed that more housing would be created as the population increases. This allowed calculation of the resources required to house shelter dogs.Baseline owned dog population (O),$$frac{dO}{dt}={r}_{o} times Otimes (1-frac{O}{{K}_{o}})+beta times H+delta times S-alpha times O-gamma times O$$
(3)
The owned dog population (Eq. 3) increases through the owned dog growth rate (ro), adoption of shelter dogs (β), and adoption of free-roaming dogs (δ); and decreases through abandonment/roaming (α) and relinquishment (γ) of owned dogs to the shelter dog population. The growth rate of the owned dog population (ro) combines the birth, death, and acquisition rates from sources other than the street or shelters (e.g. breeders, friends) and was modelled as density dependent by the limit to growth logistic formula (1 − O/Ko).Parameter estimatesDetailed descriptions of parameter estimates are provided in the supplementary information. The simulated environment was based on the city of Lviv, Ukraine. This city has an area of 182 km2 and a human population size of 717,803. Parameters were estimated from literature, where possible, and converted to monthly rates (Table 1). Initial sizes of the dog populations were estimated for the baseline simulation, based on our previous research in Lviv32. The carrying capacity depends on the availability of resources (i.e. food, shelter, water, and human attitudes and behaviour33) and is challenging to estimate. We assumed the initial free-roaming and owned dog populations were at carrying capacity. Initial population sizes for simulations including interventions were determined by the equilibrium population sizes from the baseline simulation (i.e. the stable population size, the points at which the populations were no longer increasing/decreasing).Table 1 Parameter description, parameter value, and minimum and maximum values used in the sensitivity analysis for the systems model.Full size tableEstimating the rate at which owned dogs are abandoned is difficult, as abandonment rates are often reported per dog-owning lifetime32,34 and owners are likely to under-report abandonment of dogs. Similarly, it is challenging to estimate the rate that owned dogs move from restricted to unrestricted (i.e. free-roaming). For simplicity, we modelled a combined abandonment/roaming rate (α) of 0.003 per month, estimated based on our previous research in Lviv and from literature34,35,36. We derive the owned dog relinquishment rate (γ) from New et al.37. We estimated shelter (β) and free-roaming adoption rates (δ) from shelter data in Lviv. We set the maximum intrinsic growth rate for the free-roaming dogs (rs) at 0.03 per month, similar to that reported in literature17,19,38. We assumed that demand for dogs was met quickly through a supply of dogs from births, breeders and friends and set a higher growth rate for the owned dog population (ro) at 0.07 per month.We assumed shelters operated with a “no-kill” policy (i.e. dogs were not killed in shelters as part of population management) and included a shelter dog death rate (µh) of 0.008 per month to incorporate deaths due to euthanasia for behavioural problems or health problems, or natural mortality. We modelled neutered free-roaming dog death rate (µn) explicitly for the CNR intervention at a minimum death rate of 0.02 per month38,39,40,41.InterventionsSix intervention scenarios were modelled (Table 2): sheltering; culling; CNR; responsible ownership; combined CNR and responsible ownership; and combined CNR and sheltering, representing interventions feasibly applied and often reported27. Table 2 outlines the equations describing each intervention. To simulate a sheltering intervention, a proportion of the free-roaming dog population was removed and added to the shelter dog population at sheltering rate (σ). In culling interventions, a proportion of the free-roaming dog population was removed through culling (χ).Table 2 Description of intervention parameters and coverages for simulations applied at continuous and annual periodicities.Full size tableFree-roaming dog population with sheltering intervention:$$frac{dS}{dt}={r}_{s}times Stimes left(1-frac{S}{{K}_{s}}right)+alpha times O-delta times S-sigma times S$$
(4)
Shelter dog population with sheltering intervention:$$frac{dH}{dt}=gamma times O-beta times H- {mu }_{h}times H+sigma times S$$
(5)
Free-roaming dog population with a culling intervention:$$frac{dS}{dt}={r}_{s}times Stimes left(1-frac{S}{{K}_{s}}right)+alpha times O-delta times S-chi times S$$
(6)
To simulate a CNR intervention, an additional subpopulation was added to the system (Eq. 7): (iv) the neutered free-roaming dog population (N; neutered, free-roaming). In this simulation, a proportion of the intact (I) free-roaming dog population was removed and added to the neutered free-roaming dog population. A neutering rate (φ) was added to the differential equations describing the intact free-roaming and the neutered free-roaming dog populations. Neutering was assumed to be lifelong (e.g. gonadectomy); a neutered free-roaming dog could not re-enter the intact free-roaming dog subpopulation. Neutered free-roaming dogs were removed from the population through the density dependent neutered dog death rate (µn); death rate increased when the population was closer to the carrying capacity. The death rate was a non-linear function of population size and carrying capacity modelled using a table lookup function (Fig. S1). Neutered free-roaming dogs were also removed through adoption to the owned dog population, and we assumed that adoption rates did not vary between neutered and intact free-roaming dogs.Neutered free-roaming dog population:$$frac{dN}{dt}=varphi times I-{mu }_{n}times N-delta times N$$
(7)
Intact free-roaming dog population with neutering intervention.$$frac{dI}{dt}={r}_{s}times Itimes left(1-frac{(I+N)}{{K}_{s}}right)+alpha times O-delta times I-varphi times I$$
(8)
To simulate a responsible ownership intervention, the baseline model was applied with decreased rate of abandonment/roaming (α) and increased rate of shelter adoption (β). To simulate combined CNR and responsible ownership, a proportion of the intact free-roaming dog population was removed through the neutering rate (φ), abandonments/roaming decreased (α) and shelter adoptions increased (β). In combined CNR and sheltering interventions, a proportion of the intact free-roaming dog population (I) was removed through neutering (φ) and added to the neutered free-roaming dog population (N), and a proportion was removed through sheltering (σ) and added to the shelter dog population (H).Intact free-roaming dog population with combined CNR and sheltering interventions:$$frac{dI}{dt}={r}_{s}times Stimes left(1-frac{(I+N)}{{K}_{s}}right)+alpha times O-delta times I-varphi times I- sigma times I$$
(9)
Intervention length, periodicity, and coverageAll simulations were run for 70 years to allow populations to reach equilibrium. It is important to note that this is a theoretical model; running the simulations for 70 years allows us to compare the interventions, but does not accurately predict the size of the dog subpopulations over this long time period. Interventions were applied for two lengths of time: (i) the full 70-year duration of the simulation; and (ii) a five-year period followed by no further intervention, to simulate a single period of investment in population management. In each of these simulations, we modelled the interventions as (i) continuous (i.e. a constant rate of e.g. neutering) and (ii) annual (i.e. intervention applied once per year). Interventions were run at low, medium, and high coverages (Table 2). As the processes are not equivalent, we apply different percentages for the intervention coverage (culling/neutering/sheltering) and the percent increase/decrease in parameter rates for the responsible ownership intervention. Intervention coverage refers to the proportion of dogs that are culled/neutered/sheltered per year (i.e. 20%, 40% and 70% annually) and, for responsible ownership interventions, the decrease in abandonment/roaming rate and increase in the adoption rate of shelter dogs (30%, 60% and 90% increase/decrease from baseline values). To model a low (20%), medium (40%) and high (70%) proportion of free-roaming dogs caught, but where half of the dogs were sheltered and half were neutered-and-returned, combined CNR and sheltering interventions were simulated at half-coverage (e.g. intervention rate of 0.7 was simulated by 0.35 neutered and 0.35 sheltered). For continuous interventions, sheltering (σ), culling (χ), and CNR (φ) were applied continuously during the length of the intervention. For annual interventions, σ, χ, and φ were applied to the ordinary differential equations using a forcing function applied at 12-month intervals. In simulations that included responsible ownership interventions, the decrease in owned dog abandonment/roaming (α) and the increase in shelter adoption (β) was assumed instantaneous and continuous (i.e. rates did not change throughout the intervention).Model outputsThe primary outcome of interest was the impact of interventions on free-roaming dog population size. For interventions applied for the duration of the simulation, we calculated: (i) equilibrium population size for each population; (ii) percent decrease in free-roaming dog population; (iii) costs of intervention in terms of staff-time; and (iv) an overall welfare score. For interventions applied for a five-year period, we also calculated: (v) minimum free-roaming dog population size and percent reduction from initial population size; and (vi) the length of time between the end of the intervention and time-point at which the free-roaming dog population reached above 20,000 dogs (the assumed initial free-roaming dog population size of Lviv, based on our previous research32, see Supplementary Information for detail).The costs of population management interventions vary by country (e.g. staff salaries vary between countries) and by the method of application (e.g. method of culling, or resources provided in a shelter). To enable a comparison of the resources required for each intervention, the staff time (staff working-months) required to achieve the intervention coverage was calculated. While this does not incorporate the full costs of an intervention, as equipment (e.g. surgical equipment), advertising campaigns, travel costs for the animal care team, and facilities (e.g. clinic or shelter costs) are not included, it can be used as a proxy for intervention cost. Using data provided from VIER PFOTEN International, we estimated the average number of staff required to catch and neuter the free-roaming dog population and to house the shelter dog population in each intervention, using this data as a proxy for catching and sheltering/culling. The number of dogs that can be cared for per shelter staff varies by shelter. To account for this, we estimated two staff-to-dog ratios (low and high). Table 3 describes the staff requirements for the different interventions.Table 3 Staff required for interventions and the number of dogs processed per staff per day.Full size tableUsing the projected population sizes, the staff time required for each staff type (e.g. number of veterinarian-months of work required) was calculated for each intervention. Relative salaries for the different staff types were estimated (Table 3). The relative salaries were used to calculate the cost of the interventions by:[Staff time required × relative salary ] × €20,000.Where €20,000 was the estimated annual salary of a European veterinarian, allowing relative staff-time costs to be compared between the different interventions. Average annual costs were reported.To provide overall welfare scores for each of the interventions, we apply the estimated welfare scores on a one to five scale, for each of the dog subpopulations, as determined by Hogasen et al. (2013)22. This scale is based on the Five Freedoms (freedom from hunger and thirst; freedom from discomfort; freedom from pain, injury, or disease; freedom to express normal behaviour; freedom from fear and distress42,43) and was calculated using expert opinions from 60 veterinarians in Italy22. The scores were weighted by the participants’ self-reported knowledge of different dog subpopulations, which resulted in the following scores: 2.8 for shelter dogs (WH); 3.5 for owned dogs (WO); 3.1 for neutered free-roaming dogs (WN); and 2.3 for intact free-roaming dogs (WI)22.Using these estimated welfare scores, we calculated an average welfare score for the total dog population based on the model’s projected population sizes for each subpopulation (Eq. 10). For interventions running for the duration of the simulation, the welfare score was calculated at the time point (t) when the population reached an equilibrium size. For interventions running for five years, the welfare score was calculated at the end of the five-year intervention. The percentage change in welfare scores from the baseline simulation were reported.$$Welfare score= frac{{H}_{t}times {W}_{H}+{O}_{t}times {W}_{O}+{N}_{t}times {W}_{N}+{I}_{t}times {W}_{I}}{{H}_{t}+{O}_{t}+{N}_{t}+{I}_{t}}$$
(10)
Model validation and sensitivity analysisA global sensitivity analysis was conducted on all parameters described in the baseline simulation and all interventions applied continuously, at high coverage, for the full duration of the simulation. A Latin square design algorithm was used in package “FME”44 to sample the parameters within their range of values (Table 1). For the global sensitivity analysis on interventions, all parameter values were varied, apart from the parameters involved in the intervention (e.g. culling, neutering, abandonment/roaming rates). The effects of altering individual parameters (local sensitivity analysis) on the population equilibrium was also examined for the baseline simulation using the Latin square design algorithm to sample each parameter, individually, within their range of values. Sensitivity analyses were run for 100 simulations over 50 years solved with 0.01 step sizes. More

Sala, E. et al. Protecting the global ocean for biodiversity, food and climate. Nature 592, 397–402 (2021).ADS 
CAS 
Article 

Google Scholar 
Gillispie, C. C., Gratton-Guinness, I. & Fox, R. Pierre-Simon Laplace, 1749-1827: A Life in Exact Science (Princeton Univ. Press, 1999).Dinmore, T. A., Duplisea, D. E., Rackham, B. D., Maxwell, D. L. & Jennings, S. Impact of a large-scale area closure on patterns of fishing disturbance and the consequences for benthic communities. ICES J. Mar. Sci. 60, 371–380 (2003).Article 

Google Scholar 
Hiddink, J. G., Hutton, T., Jennings, S. & Kaiser, M. J. Predicting the effects of area closures and fishing effort restrictions on the production, biomass, and species richness of benthic invertebrate communities. ICES J. Mar. Sci. 63, 822–830 (2006).Article 

Google Scholar 
Greenstreet, S. P. R., Fraser, H. M. & Piet, G. J. Using MPAs to address regional-scale ecological objectives in the North Sea: modelling the effects of fishing effort displacement. ICES J. Mar. Sci. 66, 90–100 (2009).Article 

Google Scholar 
Suuronen, P. et al. A path to a sustainable trawl fishery in Southeast Asia. Rev. Fish. Sci. Aquac. 28, 499–517 (2020).Article 

Google Scholar 
Amoroso, R. O. et al. Bottom trawl fishing footprints on the world’s continental shelves. Proc. Natl Acad. Sci. USA 115, E10275–E10282 (2018).CAS 
Article 

Google Scholar 
Atwood, T. B., Witt, A., Mayorga, J., Hammill, E. & Sala, E. Global patterns in marine sediment carbon stocks. Front. Mar. Sci. 7, 165 (2020).Article 

Google Scholar 
Smeaton, C., Hunt, C. A., Turrell, W. R. & Austin, W. E. N. Marine sedimentary carbon stocks of the United Kingdom’s exclusive economic zone. Front. Earth Sci. 9, 593324 (2021).Article 

Google Scholar 
Legge, O. et al. Carbon on the northwest European shelf: contemporary budget and future influences. Front. Mar. Sci. 7, 143 (2020).Article 

Google Scholar 
Melnychuk, M. C. et al. Identifying management actions that promote sustainable fisheries. Nat. Sustain. 4, 440–449 (2021).Article 

Google Scholar  More

Bowler, D. E. et al. Mapping human pressures on biodiversity across the planet uncovers anthropogenic threat complexes. People Nat. 2, 380–394 (2020).Article 

Google Scholar 
Persson, L. et al. Outside the safe operating space of the planetary boundary for novel entities. Environ. Sci. Technol. https://doi.org/10.1021/acs.est.1c04158 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
United Nations Environment Programme (UNEP). 2019. Global Mercury Assessment 2018. UN Environment Programme, Chemicals and Health Branch Geneva, Switzerland. https://www.unep.org/resources/publication/global-mercury-assessment-2018Rimmer, C. C., Miller, E. K., McFarland, K. P., Taylor, R. J. & Faccio, S. D. Mercury bioaccumulation and trophic transfer in the terrestrial food web of a montane forest. Ecotoxicology 19, 697–709 (2010).CAS 
PubMed 
Article 

Google Scholar 
Cristol, D. A. et al. The movement of aquatic mercury through terrestrial food webs. Science 320, 335 (2008).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Evers, D. The effects of methylmercury on wildlife: A comprehensive review and approach for interpretation. Encycl. Anthropocene 5, 181–194 (2018).Article 

Google Scholar 
Whitney, M. C. & Cristol, D. A. Impacts of sublethal mercury exposure on birds: a detailed review. Rev. Environ. Contam. Toxicol. 244, 113–163 (2017).
Google Scholar 
Seewagen, C. L. Threats of environmental mercury to birds: Knowledge gaps and priorities for future research. Bird Conserv. Int. 20, 112–123 (2010).Article 

Google Scholar 
Seewagen, C. L. The threat of global mercury pollution to bird migration: Potential mechanisms and current evidence. Ecotoxicology 29, 1254–1267 (2020).CAS 
PubMed 
Article 

Google Scholar 
Ma, Y., Branfireun, B. A., Hobson, K. A. & Guglielmo, C. G. Evidence of negative seasonal carry-over effects of breeding ground mercury exposure on survival of migratory songbirds. J. Avian Biol. 49, jav-01656 (2018).Article 

Google Scholar 
Newton, I. Can conditions experienced during migration limit the population levels of birds?. J. Ornithol. 147, 146–166 (2006).Article 

Google Scholar 
Klaassen, M., Hoye, B. J., Nolet, B. A. & Buttemer, W. A. Ecophysiology of avian migration in the face of current global hazards. Philos. Trans. R. Soc. B 367, 1719–1732 (2020).Article 

Google Scholar 
Zurell, D., Graham, C. H., Gallien, L., Thuiller, W. & Zimmermann, N. E. Long-distance migratory birds threatened by multiple independent risks from global change. Nat. Clim. Chang. 8, 992–996 (2018).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Seewagen, C. L., Ma, Y., Morbey, Y. E. & Guglielmo, C. G. Stopover departure behavior and flight orientation of spring-migrant Yellow-rumped Warblers (Setophaga coronata) experimentally exposed to methylmercury. J. Ornithol. 160, 617–624 (2019).Article 

Google Scholar 
Seewagen, C. L. Blood mercury levels and the stopover refueling performance of a long-distance migratory songbird. Can. J. Zool. 91, 41–45 (2013).CAS 
Article 

Google Scholar 
Adams, E. M., Williams, K. A., Olsen, B. J. & Evers, D. C. Mercury exposure in migrating songbirds: Correlations with physical condition. Ecotoxicology 29, 1240–1253 (2020).CAS 
PubMed 
Article 

Google Scholar 
Ma, Y., Perez, C. R., Branfireun, B. A. & Guglielmo, C. G. Dietary exposure to methylmercury affects flight endurance in a migratory songbird. Environ. Pollut. 234, 894–901 (2018).CAS 
PubMed 
Article 

Google Scholar 
Gerson, A. R., Cristol, D. A. & Seewagen, C. L. Environmentally relevant methylmercury exposure reduces the metabolic scope of a model songbird. Environ. Pollut. 246, 790–796 (2019).CAS 
PubMed 
Article 

Google Scholar 
Jenni, L. & Jenni-Eiermann, S. Fuel supply and metabolic constraints in migrating birds. J. Avian Biol. 29, 521–552 (1998).Article 

Google Scholar 
McWilliams, S. R., Guglielmo, C., Pierce, B. & Klaassen, M. Flying, fasting, and feeding in birds during migration: A nutritional and physiological ecology perspective. J. Avian Biol. 35, 377–393 (2004).Article 

Google Scholar 
Guglielmo, C. G. Move that fatty acid: Fuel selection and transport in migratory birds and bats. Integr. Comp. Biol. 50, 336–345 (2010).PubMed 
Article 

Google Scholar 
Guglielmo, C. G. Obese super athletes: Fat-fueled migration in birds and bats. J. Exp. Biol. 221(Suppl_1), 165753 (2018).Article 

Google Scholar 
Kawakami, T. et al. Differential effects of cobalt and mercury on lipid metabolism in the white adipose tissue of high-fat diet-induced obesity mice. Toxicol. Appl. Pharmacol. 258, 32–42 (2012).CAS 
PubMed 
Article 

Google Scholar 
Yadetie, F. et al. Global transcriptome analysis of Atlantic cod (Gadus morhua) liver after in vivo methylmercury exposure suggests effects on energy metabolism pathways. Aquat. Toxicol. 126, 314–325 (2013).CAS 
PubMed 
Article 

Google Scholar 
Park, K. & Seo, E. Association between toenail mercury and metabolic syndrome is modified by selenium. Nutrients 8, 424 (2016).PubMed Central 
Article 
CAS 

Google Scholar 
Caito, S. W., Newell-Caito, J., Martell, M., Crawford, N. & Aschner, M. Methylmercury induces metabolic alterations in Caenorhabditis elegans: Role for C/EBP transcription factor. Toxicol. Sci. 174, 112–123 (2020).CAS 
PubMed 
Article 

Google Scholar 
Edmonds, S. T., O’Driscoll, N. J., Hillier, N. K., Atwood, J. L. & Evers, D. C. Factors regulating the bioavailability of methylmercury to breeding rusty blackbirds in northeastern wetlands. Environ. Pollut. 171, 148–154 (2012).CAS 
PubMed 
Article 

Google Scholar 
Rowse, L. M., Rodewald, A. D., Mažeika, S. & Sullivan, P. Pathways and consequences of contaminant flux to Acadian flycatchers (Empidonax virescens) in urbanizing landscapes of Ohio, USA. Sci. Total Environ. 485, 461–467 (2014).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Marsh, R. L. Catabolic enzyme activities in relation to premigratory fattening and muscle hypertrophy in the gray catbird (Dumetella carolinensis). J. Comp. Physiol. 141, 417–423 (1981).CAS 
Article 

Google Scholar 
Guglielmo, C. G., Haunerland, N. H., Hochachka, P. W. & Williams, T. D. Seasonal dynamics of flight muscle fatty acid binding protein and catabolic enzymes in a migratory shorebird. Am. J. Physiol.-Regul. Integr. Comp. Physiol. 282(5), R1405–R1413 (2002).CAS 
PubMed 
Article 

Google Scholar 
Maillet, D. & Weber, J. M. Relationship between n-3 PUFA content and energy metabolism in the flight muscles of a migrating shorebird: Evidence for natural doping. J. Exp. Biol. 210, 413–420 (2007).CAS 
PubMed 
Article 

Google Scholar 
Weber, J. M. Metabolic fuels: Regulating fluxes to select mix. J. Exp. Biol. 214, 286–294 (2011).CAS 
PubMed 
Article 

Google Scholar 
Feige, J. N., Gelman, L., Michalik, L., Desvergne, B. & Wahli, W. From molecular action to physiological outputs: Peroxisome proliferator-activated receptors are nuclear receptors at the crossroads of key cellular functions. Prog. Lipid. Res. 45, 120–159 (2006).CAS 
PubMed 
Article 

Google Scholar 
Bensinger, S. J. & Tontonoz, P. Integration of metabolism and inflammation by lipid-activated nuclear receptors. Nature 454, 470–477. https://doi.org/10.1038/nature07202 (2008).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Ynalvez, R., Gutierrez, J. & Gonzalez-Cantu, H. Mini-review: Toxicity of mercury as a consequence of enzyme alteration. Biometals 29, 781–788 (2016).CAS 
PubMed 
Article 

Google Scholar 
Gerson, A. R. & Guglielmo, C. G. Energetics and metabolite profiles during early flight in American robins (Turdus Migratorius). J. Comp. Physiol. B. 183, 983–991 (2013).CAS 
PubMed 
Article 

Google Scholar 
Price, E. R., McFarlan, J. T. & Guglielmo, C. G. Preparing for migration? The effects of photoperiod and exercise on muscle oxidative enzymes, lipid transporters, and phospholipids in white-crowned sparrows. Physiol. Biochem. Zool. 83, 252–262 (2010).CAS 
PubMed 
Article 

Google Scholar 
Bradley, S. S., Dick, M. F., Guglielmo, C. G. & Timoshenko, A. V. Seasonal and flight-related variation of galectin expression in heart, liver and flight muscles of yellow-rumped warblers (Setophaga coronata). Glycoconj. J. 34, 603–611 (2017).CAS 
PubMed 
Article 

Google Scholar 
McFarlan, J. T., Bonen, A. & Guglielmo, C. G. Seasonal upregulation of fatty acid transporters in flight muscles of migratory white-throated sparrows (Zonotrichia albicollis). J. Exp. Biol. 212, 2934–2940 (2009).CAS 
PubMed 
Article 

Google Scholar 
Zhang, Y., King, M. O., Harmon, E., Eyster, K. & Swanson, D. L. Migration-induced variation of fatty acid transporters and cellular metabolic intensity in passerine birds. J. Comp. Physiol. B. 185, 797–810 (2015).CAS 
PubMed 
Article 

Google Scholar 
Dick, M. F. & Guglielmo, C. G. Dietary polyunsaturated fatty acids influence flight muscle oxidative capacity but not endurance flight performance in a migratory songbird. Am. J. Physiol.-Regul. Integr. Compar. Physiol. 316(4), R362–R375 (2019).CAS 
Article 

Google Scholar 
Schmittgen, T. D. & Livak, K. J. Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 3, 1101–1108 (2008).CAS 
PubMed 
Article 

Google Scholar 
Bittencourt, L. O. et al. Oxidative biochemistry disbalance and changes on proteomic profile in salivary glands of rats induced by chronic exposure to methylmercury. Oxid. Med. Cell. Longev. https://doi.org/10.1155/2017/5653291 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Shi, Q., Sun, N., Kou, H., Wang, H. & Zhao, H. Chronic effects of mercury on Bufo gargarizans larvae: Thyroid disruption, liver damage, oxidative stress and lipid metabolism disorder. Ecotoxicol. Environ. Saf. 164, 500–509 (2018).CAS 
PubMed 
Article 

Google Scholar 
Nøstbakken, O. J. et al. Dietary methylmercury alters the proteome in Atlantic salmon (Salmo salar) kidney. Aquat. Toxicol. 108, 70–77 (2012).PubMed 
Article 
CAS 

Google Scholar 
Zink, E. M. Comparison of the mercury induced proteomes of Escherichia coli MG1655 with and without the NR1 plasmid. MSc thesis, Washington State University, Pullman, WA (2009).Lundgren, B. O. & Kiessling, K. H. Seasonal variation in catabolic enzyme activities in breast muscle of some migratory birds. Oecologia 66, 468–471 (1985).ADS 
PubMed 
Article 

Google Scholar 
Banerjee, S. & Chaturvedi, C. M. Migratory preparation associated alterations in pectoralis muscle biochemistry and proteome in Palearctic-Indian emberizid migratory finch, red-headed bunting, Emberiza bruniceps. Comp. Biochem. Physiol. D Genom. Proteom. 17, 9–25 (2016).CAS 

Google Scholar 
Dick, M. F. The long haul: migratory flight preparation and performance in songbirds. Ph.D. dissertation, University of Western Ontario, London, Canada (2017).Driedzic, W. R., Crowe, H. L., Hicklin, P. W. & Sephton, D. H. Adaptations in pectoralis muscle, heart mass, and energy metabolism during premigratory fattening in semipalmated sandpipers (Calidris pusilla). Can. J. Zool. 71, 1602–1608 (1993).Article 

Google Scholar 
De Moranville, K. J. et al. PPAR expression, muscle size and metabolic rates across the gray catbird’s annual cycle are greatest in preparation for fall migration. J. Exper. Biol. 222, 198028 (2019).Article 

Google Scholar 
Zajac, D. M., Cerasale, D. J., Landman, S. & Guglielmo, C. G. Behavioral and physiological effects of photoperiod-induced migratory state and leptin on Zonotrichia albicollis: II. Effects on fatty acid metabolism. Gen. Comp. Endocrinol. 174, 269–275 (2011).CAS 
PubMed 
Article 

Google Scholar 
Tinant, G. et al. Methylmercury displays pro-adipogenic properties in rainbow trout preadipocytes. Chemosphere 263, 127917 (2021).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Cambier, S. et al. At environmental doses, dietary methylmercury inhibits mitochondrial energy metabolism in skeletal muscles of the zebra fish (Danio rerio). Int. J. Biochem. Cell Biol. 41, 791–799 (2009).CAS 
PubMed 
Article 

Google Scholar 
Ferain, A. et al. Transcriptional effects of phospholipid fatty acid profile on rainbow trout liver cells exposed to methylmercury. Aquat. Toxicol. 199, 174–187 (2018).CAS 
PubMed 
Article 

Google Scholar 
Börchers, T., Højrup, P., Nielsen, S. U., Roepstorff, P., Spener, F., Knudsen, J. Revision of the amino acid sequence of human heart fatty acid-binding protein. In Cellular Fatty Acid-binding Proteins 127–133 (Springer, Boston, 1990).Dörmann, P. et al. Amino acid exchange and covalent modification by cysteine and glutathione explain isoforms of fatty acid-binding protein occurring in bovine liver. J. Biol. Chem. 268, 16286–16292 (1993).PubMed 
Article 

Google Scholar 
Su, X. & Abumrad, N. A. Cellular fatty acid uptake: A pathway under construction. Trends Endocrinol. Metab. 20(2), 72–77 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
van Oort, M. M. et al. Each of the four intracellular cysteines of CD36 is essential for insulin-or AMP-activated protein kinase-induced CD36 translocation. Arch. Physiol. Biochem. 120, 40–49 (2014).PubMed 
Article 
CAS 

Google Scholar 
Wang, G., Bonkovsky, H. L., de Lemos, A. & Burczynski, F. J. Recent insights into the biological functions of liver fatty acid binding protein 1. J. Lipid Res. 56, 2238–2247 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Vallee, B. L. & Ulmer, D. D. Biochemical effects of mercury, cadmium, and lead. Annu. Rev. Biochem. 41, 91–128 (1972).CAS 
PubMed 
Article 

Google Scholar 
Aschner, M. & Syversen, T. Methylmercury: Recent advances in the understanding of its neurotoxicity. Ther. Drug Monit. 27, 278–283 (2005).PubMed 
PubMed Central 
Article 

Google Scholar 
Kenow, K. P., Meyer, M. W., Hines, R. K. & Karasov, W. H. Distribution and accumulation of mercury in tissues of captive-reared common loon (Gavia immer) chicks. Environ. Toxicol. Chem. 26, 1047–1055 (2007).CAS 
PubMed 
Article 

Google Scholar 
Varian-Ramos, C. W., Whitney, M., Rice, G. W. & Cristol, D. A. Form of dietary methylmercury does not affect total mercury accumulation in the tissues of zebra finch. Bull. Environ. Contam. Toxicol. 99, 1–8 (2017).CAS 
PubMed 
Article 

Google Scholar 
Rizzetti, D. A. et al. Chronic mercury at low doses impairs white adipose tissue plasticity. Toxicology 418, 41–50 (2019).CAS 
PubMed 
Article 

Google Scholar 
Richter, C. A. et al. Methylmercury-induced changes in gene transcription associated with neuroendocrine disruption in largemouth bass (Micropterus salmoides). Gen. Comp. Endocrinol. 203, 215–224 (2014).CAS 
PubMed 
Article 

Google Scholar 
Barnes, D. M., Hanlon, P. R. & Kircher, E. A. Effects of inorganic HgCl2 on adipogenesis. Toxicol. Sci. 75(2), 368–377 (2003).CAS 
PubMed 
Article 

Google Scholar 
Corder, K. R., DeMoranville, K. J., Russell, D. E., Huss, J. M. & Schaeffer, P. J. Annual life-stage regulation of lipid metabolism and storage and association with PPARs in a migrant species: the gray catbird (Dumetella carolinensis). J. Exp. Biol. 219, 3391–3398 (2016).PubMed 

Google Scholar 
DeMoranville, K. J., Carter, W. A., Pierce, B. J. & McWilliams, S. R. Flight training in a migratory bird drives metabolic gene expression in the flight muscle but not liver, and dietary fat quality influences select genes. Am. J. Physiol.-Regul. Integr. Compar. Physiol. 319(6), R637–R652 (2020).CAS 
Article 

Google Scholar 
Gavrilova, O. et al. Liver peroxisome proliferator-activated receptor γ contributes to hepatic steatosis, triglyceride clearance, and regulation of body fat mass. J. Biol. Chem. 278(36), 34268–34276 (2003).CAS 
PubMed 
Article 

Google Scholar 
Bedoucha, M., Atzpodien, E. & Boelsterli, U. A. Diabetic KKAy mice exhibit increased hepatic PPARγ1 gene expression and develop hepatic steatosis upon chronic treatment with antidiabetic thiazolidinediones. J. Hepatol. 35, 17–23 (2001).CAS 
PubMed 
Article 

Google Scholar 
Egeler, O., Williams, T. D. & Guglielmo, C. G. Modulation of lipogenic enzymes, fatty acid synthase and Δ 9-desaturase, in relation to migration in the western sandpiper (Calidris mauri). J. Comp. Physiol. B 170, 169–174 (2000).CAS 
PubMed 
Article 

Google Scholar 
Klaper, R. et al. Use of a 15k gene microarray to determine gene expression changes in response to acute and chronic methylmercury exposure in the fathead minnow (Pimephales promelas). J. Fish Biol. 72, 2207–2280 (2008).CAS 
Article 

Google Scholar 
Calow, P. Physiological costs of combating chemical toxicants: Ecological implications. Comp. Biochem. Physiol. C 100, 3–6 (1991).CAS 
PubMed 
Article 

Google Scholar 
Spalding, M. G. et al. Histologic, neurologic, and immunologic effects of methylmercury in captive great egrets. J. Wildl. Dis. 36, 423–435 (2000).CAS 
PubMed 
Article 

Google Scholar 
Carlson, J. R., Cristol, D. & Swaddle, J. P. Dietary mercury exposure causes decreased escape takeoff flight performance and increased molt rate in European starlings (Sturnus vulgaris). Ecotoxicology 23, 1464–1473 (2014).CAS 
PubMed 
Article 

Google Scholar 
Faaborg, J. et al. Conserving migratory land birds in the New World: Do we know enough?. Ecol. Appl. 20, 398–418 (2010).PubMed 
Article 

Google Scholar 
Duijns, S. et al. Body condition explains migratory performance of a long-distance migrant. Proc. R. Soc. B https://doi.org/10.1098/rspb.2017.1374 (2017).Article 
PubMed 
PubMed Central 

Google Scholar  More

Neufeld JD, Wagner M, Murrell JC. Who eats what, where and when? Isotope-labelling experiments are coming of age. ISME J. 2007;1:103–10.CAS 
PubMed 
Article 

Google Scholar 
Boschker HTS, Nold SC, Wellsbury P, Bos D, de Graaf W, Pel R, et al. Direct linking of microbial populations to specific biogeochemical processes by 13C-labelling of biomarkers. Nature. 1998;392:801–5Jehmlich N, Schmidt F, von Bergen M, Richnow H-H, Vogt C. Protein-based stable isotope probing (Protein-SIP) reveals active species within anoxic mixed cultures. ISME J. 2008;2:1122–33.CAS 
PubMed 
Article 

Google Scholar 
Radajewski S, Ineson P, Parekh NR, Colin Murrell J. Stable-isotope probing as a tool in microbial ecology. Nature. 2000;403:646–9.CAS 
PubMed 
Article 

Google Scholar 
Manefield M, Whiteley AS, Griffiths RI, Bailey MJ. RNA stable isotope probing, a novel means of linking microbial community function to phylogeny. Appl Environ Microbiol. 2002;68:5367–73.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Berry D, Mader E, Lee TK, Woebken D, Wang Y, Zhu D, et al. Tracking heavy water (D2O) incorporation for identifying and sorting active microbial cells. Proc Natl Acad Sci USA. 2015;112:E194–203.CAS 
PubMed 

Google Scholar 
Jehmlich N, Vogt C, Lünsmann V, Richnow HH, von Bergen M. Protein-SIP in environmental studies. Curr Opin Biotechnol. 2016;41:26–33.CAS 
PubMed 
Article 

Google Scholar 
Haichar, FEZ, Achouak W, Christen R, Heulin T, et al. Identification of cellulolytic bacteria in soil by stable isotope probing. Environ Microbiol. 2007;9:625–34Rangel-Castro JI, Ignacio Rangel-Castro J, Killham K, Ostle N, Nicol GW, Anderson IC, et al. Stable isotope probing analysis of the influence of liming on root exudate utilization by soil microorganisms. Environ Microbiol. 2005;7:828–38.CAS 
PubMed 
Article 

Google Scholar 
Wang Y, Song Y, Tao Y, Muhamadali H, Goodacre R, Zhou N-Y, et al. Reverse and multiple stable isotope probing to study bacterial metabolism and interactions at the single cell level. Anal Chem. 2016;88:9443–50.CAS 
PubMed 
Article 

Google Scholar 
Sharma K, Palatinszky M, Nikolov G, Berry D, Shank EA. Transparent soil microcosms for live-cell imaging and non-destructive stable isotope probing of soil microorganisms. Elife. 2020;9:e56275.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lee KS, Landry Z, Pereira FC, Wagner M, Berry D, Huang WE, et al. Raman microspectroscopy for microbiology. Nat. Rev. Methods Primers. 2021;1:80.CAS 
Article 

Google Scholar 
Hatzenpichler R, Krukenberg V, Spietz RL, Jay ZJ. Next-generation physiology approaches to study microbiome function at single cell level. Nat Rev Microbiol. 2020;18:241–56.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wagner M. Single-cell ecophysiology of microbes as revealed by raman microspectroscopy or secondary ion mass spectrometry imaging. Ann Rev Microbiol. 2009;63:411–29Lennon JT, Jones SE. Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat Rev Microbiol. 2011;9:119–30.CAS 
PubMed 
Article 

Google Scholar 
Lewis K. Persister cells, dormancy and infectious disease. Nat Rev Microbiol. 2006;5:48–56.PubMed 
Article 
CAS 

Google Scholar 
Nielsen KM, Johnsen PJ, Bensasson D, Daffonchio D. Release and persistence of extracellular DNA in the environment. Environ Biosafety Res. 2007;6:37–53.CAS 
PubMed 
Article 

Google Scholar 
Blazewicz SJ, Barnard RL, Daly RA, Firestone MK. Evaluating rRNA as an indicator of microbial activity in environmental communities: limitations and uses. ISME J. 2013;7:2061–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Nocker A, Sossa-Fernandez P, Burr MD, Camper AK. Use of propidium monoazide for live/dead distinction in microbial ecology. Appl Environ Microbiol. 2007;73:5111–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tawakoli PN, Al-Ahmad A, Hoth-Hannig W, Hannig M, Hannig C. Comparison of different live/dead stainings for detection and quantification of adherent microorganisms in the initial oral biofilm. Clin Oral Investig. 2013;17:841–50.CAS 
PubMed 
Article 

Google Scholar 
Netuschil L, Auschill TM, Sculean A, Arweiler NB. Confusion over live/dead stainings for the detection of vital microorganisms in oral biofilms-which stain is suitable? BMC Oral Health. 2014;14:2.PubMed 
PubMed Central 
Article 

Google Scholar 
Hatzenpichler R, Connon SA, Goudeau D, Malmstrom RR, Woyke T, Orphan VJ. Visualizing in situ translational activity for identifying and sorting slow-growing archaeal−bacterial consortia. Proc Natl Acad Sci USA. 2016;113:E4069–78.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kuru E, Hughes HV, Brown PJ, Hall E, Tekkam S, Cava F, et al. In Situ probing of newly synthesized peptidoglycan in live bacteria with fluorescent D-amino acids. Angew Chem Int Ed Engl. 2012;51:12519–23.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kopf SH, McGlynn SE, Green-Saxena A, Guan Y, Newman DK, Orphan VJ. Heavy water and15N labelling with NanoSIMS analysis reveals growth rate-dependent metabolic heterogeneity in chemostats. Environ Microbiol. 2015;17:2542–56Kopf SH, Sessions AL, Cowley ES, Reyes C, Van Sambeek L, Hu Y, et al. Trace incorporation of heavy water reveals slow and heterogeneous pathogen growth rates in cystic fibrosis sputum. Proc Natl Acad Sci USA. 2016;113:E110–6.CAS 
PubMed 
Article 

Google Scholar 
Neubauer C, Kasi AS, Grahl N, Sessions AL, Kopf SH, Kato R, et al. Refining the Application of Microbial Lipids as Tracers of Staphylococcus aureus Growth Rates in Cystic Fibrosis Sputum. J Bacteriol. 2018;200:e00365–18.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Haider S, Wagner M, Schmid MC, Sixt BS, Christian JG, Häcker G, et al. Raman microspectroscopy reveals long-term extracellular activity of Chlamydiae. Mol Microbiol. 2010;77:687–700.CAS 
PubMed 
Article 

Google Scholar 
Kloehn J, Boughton BA, Saunders EC, O’Callaghan S, Binger KJ, McConville MJ. Identification of Metabolically Quiescent Leishmania mexicana Parasites in Peripheral and Cured Dermal Granulomas Using Stable Isotope Tracing Imaging Mass Spectrometry. mBio. 2021;12:e00129–21.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kong L, Setlow P, Li Y-Q. Direct analysis of water content and movement in single dormant bacterial spores using confocal Raman microspectroscopy and Raman imaging. Anal Chem. 2013;85:7094–101.CAS 
PubMed 
Article 

Google Scholar 
Knudsen SM, Cermak N, Delgado FF, Setlow B, Setlow P, Manalis SR. Water and small-molecule permeation of dormant Bacillus subtilis spores. J Bacteriol. 2016;198:168–77.CAS 
PubMed 
Article 

Google Scholar 
Chen D, Huang S-S, Li Y-Q. Real-time detection of kinetic germination and heterogeneity of single Bacillus spores by laser tweezers Raman spectroscopy. Anal Chem. 2006;78:6936–41.CAS 
PubMed 
Article 

Google Scholar 
Devictor V, Clavel J, Julliard R, Lavergne S, Mouillot D, Thuiller W, et al. Defining and measuring ecological specialization. J Appl Ecol. 2010;47:15–25.Article 

Google Scholar 
Pereira FC, Berry D. Microbial nutrient niches in the gut. Environ Microbiol. 2017;19:1366–78.PubMed 
PubMed Central 
Article 

Google Scholar 
Shakya M, Lo C-C, Chain PSG. Advances and challenges in metatranscriptomic analysis. Front Genet. 2019;10:904.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Berry D, Loy A. Stable-Isotope probing of human and animal microbiome function. Trends Microbiol. 2018;26:999–1007.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Terrado R, Pasulka AL, Lie AA-Y, Orphan VJ, Heidelberg KB, Caron DA. Autotrophic and heterotrophic acquisition of carbon and nitrogen by a mixotrophic chrysophyte established through stable isotope analysis. ISME J. 2017;11:2022–34.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Dekas AE, Parada AE, Mayali X, Fuhrman JA, Wollard J, Weber PK, et al. Characterizing Chemoautotrophy and Heterotrophy in Marine Archaea and Bacteria With Single-Cell Multi-isotope NanoSIP. Front Microbiol. 2019;10:2682.PubMed 
PubMed Central 
Article 

Google Scholar 
Wegener G, Bausch M, Holler T, Thang NM, Mollar XP, Kellermann MY, et al. Assessing sub-seafloor microbial activity by combined stable isotope probing with deuterated water and 13C-bicarbonate. Environ Microbiol. 2019;14:1517–27Jing X, Gou H, Gong Y, Su X, Xu L, Ji Y, et al. Raman-activated cell sorting and metagenomic sequencing revealing carbon-fixing bacteria in the ocean. Environ Microbiol. 2018;20:2241–55.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Xu J, Zhu D, Ibrahim AD, Allen CCR, Gibson CM, Fowler PW, et al. Raman deuterium isotope probing reveals microbial metabolism at the single-cell level. Anal Chem. 2017;89:13305–12.CAS 
PubMed 
Article 

Google Scholar 
Zhang M, Hong W, Abutaleb NS, Li J, Dong P-T, Zong C, et al. Rapid determination of antimicrobial susceptibility by stimulated Raman scattering imaging of D2O metabolic incorporation in a single bacterium. Adv Chem Microsc Life Sci Transl Med. 2021.Lima C, Muhamadali H, Xu Y, Kansiz M, Goodacre R. Imaging Isotopically Labeled Bacteria at the Single-Cell Level Using High-Resolution Optical Infrared Photothermal Spectroscopy. Anal Chem. 2021;93:3082–8.CAS 
PubMed 
Article 

Google Scholar 
Ackermann M. A functional perspective on phenotypic heterogeneity in microorganisms. Nat Rev Microbiol. 2015;13:497–508.CAS 
PubMed 
Article 

Google Scholar 
Balaban NQ, Merrin J, Chait R, Kowalik L, Leibler S. Bacterial persistence as a phenotypic switch. Science. 2004;305:1622–5.CAS 
PubMed 
Article 

Google Scholar 
Maamar H, Raj A, Dubnau D. Noise in gene expression determines cell fate in Bacillus subtilis. Science. 2007;317:526–9.CAS 
PubMed 
Article 

Google Scholar 
Emonet T, Cluzel P. Relationship between cellular response and behavioral variability in bacterial chemotaxis. Proc Natl Acad Sci USA. 2008;105:3304–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ozbudak EM, Thattai M, Lim HN, Shraiman BI, Van Oudenaarden A. Multistability in the lactose utilization network of Escherichia coli. Nature. 2004;427:737–40.CAS 
PubMed 
Article 

Google Scholar 
Kiviet DJ, Nghe P, Walker N, Boulineau S, Sunderlikova V, Tans SJ. Stochasticity of metabolism and growth at the single-cell level. Nature. 2014;514:376–9.CAS 
PubMed 
Article 

Google Scholar 
Kotte O, Volkmer B, Radzikowski JL, Heinemann M. Phenotypic bistability in Escherichia coli’s central carbon metabolism. Mol Syst Biol. 2014;10:736.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
New AM, Cerulus B, Govers SK, Perez-Samper G, Zhu B, Boogmans S, et al. Different levels of catabolite repression optimize growth in stable and variable environments. PLoS Biol. 2014;12:e1001764.PubMed 
PubMed Central 
Article 

Google Scholar 
Solopova A, van Gestel J, Weissing FJ, Bachmann H, Teusink B, Kok J, et al. Bet-hedging during bacterial diauxic shift. Proc Natl Acad Sci USA. 2014;111:7427–32.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Schreiber F, Littmann S, Lavik G, Escrig S, Meibom A, Kuypers MMM, et al. Phenotypic heterogeneity driven by nutrient limitation promotes growth in fluctuating environments. Nat Microbiol. 2016;1:16055.CAS 
PubMed 
Article 

Google Scholar 
Nikolic N, Schreiber F, Dal Co A, Kiviet DJ, Bergmiller T, Littmann S, et al. Cell-to-cell variation and specialization in sugar metabolism in clonal bacterial populations. PLoS Genet. 2017;13:e1007122.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Takhaveev V, Heinemann M. Metabolic heterogeneity in clonal microbial populations. Curr Opin Microbiol. 2018;45:30–8.CAS 
PubMed 
Article 

Google Scholar 
Altschuler SJ, Wu LF. Cellular heterogeneity: do differences make a difference? Cell. 2010;141:559–63.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Beaumont HJE, Gallie J, Kost C, Ferguson GC, Rainey PB. Experimental evolution of bet hedging. Nature. 2009;462:90–3.CAS 
PubMed 
Article 

Google Scholar 
Calabrese F, Voloshynovska I, Musat F, Thullner M, Schlömann M, Richnow HH, et al. Quantitation and comparison of phenotypic heterogeneity among single cells of monoclonal microbial populations. Front Microbiol. 2019;10:2814.PubMed 
PubMed Central 
Article 

Google Scholar 
Zimmermann M, Escrig S, Hübschmann T, Kirf MK, Brand A, Inglis RF, et al. Phenotypic heterogeneity in metabolic traits among single cells of a rare bacterial species in its natural environment quantified with a combination of flow cell sorting and NanoSIMS. Front Microbiol. 2015;6:243.PubMed 
PubMed Central 
Article 

Google Scholar 
Zimmermann M, Escrig S, Lavik G, Kuypers MMM, Meibom A, Ackermann M, et al. Substrate and electron donor limitation induce phenotypic heterogeneity in different metabolic activities in a green sulphur bacterium. Environ Microbiol Rep. 2018;10:179–83.CAS 
PubMed 
Article 

Google Scholar 
Sheik AR, Muller EE, Audinot J-N, Lebrun LA, Grysan P, Guignard C, et al. In situ phenotypic heterogeneity among single cells of the filamentous bacterium Candidatus Microthrix parvicella. ISME J. 2016;10:1274–9.CAS 
PubMed 
Article 

Google Scholar 
Seyedsayamdost MR, Case RJ, Kolter R, Clardy J. The Jekyll-and-Hyde chemistry of Phaeobacter gallaeciensis. Nat Chem. 2011;3:331–5.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ferrier-Pagès C, Leal MC. Stable isotopes as tracers of trophic interactions in marine mutualistic symbioses. Ecol Evol. 2019;9:723–40.PubMed 
Article 

Google Scholar 
Pasulka AL, Thamatrakoln K, Kopf SH, Guan Y, Poulos B, Moradian A, et al. Interrogating marine virus-host interactions and elemental transfer with BONCAT and nanoSIMS-based methods. Environ Microbiol. 2018;20:671–92.CAS 
PubMed 
Article 

Google Scholar 
Kopp C, Domart-Coulon I, Escrig S, Humbel BM, Hignette M, Meibom A. Subcellular investigation of photosynthesis-driven carbon assimilation in the symbiotic reef coral Pocillopora damicornis. mBio. 2015;6:e02299–14.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rädecker N, Pogoreutz C, Gegner HM, Cárdenas A, Roth F, Bougoure J, et al. Heat stress destabilizes symbiotic nutrient cycling in corals. Proc Natl Acad Sci U S A. 2021;118:e2022653118.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Krueger T, Bodin J, Horwitz N, Loussert-Fonta C, Sakr A, Escrig S, et al. Temperature and feeding induce tissue level changes in autotrophic and heterotrophic nutrient allocation in the coral symbiosis – a NanoSIMS study. Sci Rep. 2018;8:12710.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Gibbin E, Gavish A, Krueger T, Kramarsky-Winter E, Shapiro O, Guiet R, et al. Vibrio coralliilyticus infection triggers a behavioural response and perturbs nutritional exchange and tissue integrity in a symbiotic coral. ISME J. 2019;13:989–1003.CAS 
PubMed 
Article 

Google Scholar 
Rix L, Ribes M, Coma R, Jahn MT, de Goeij JM, van Oevelen D, et al. Heterotrophy in the earliest gut: a single-cell view of heterotrophic carbon and nitrogen assimilation in sponge-microbe symbioses. ISME J. 2020;14:2554–67.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Thomas T, Moitinho-Silva L, Lurgi M, Björk JR, Easson C, Astudillo-García C, et al. Diversity, structure and convergent evolution of the global sponge microbiome. Nat Commun. 2016;7:11870.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mills MM, Turk-Kubo KA, van Dijken GL, Henke BA, Harding K, Wilson ST, et al. Unusual marine cyanobacteria/haptophyte symbiosis relies on N2 fixation even in N-rich environments. ISME J. 2020;14:2395–406.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Turk-Kubo KA, Mills MM, Arrigo KR, van Dijken G, Henke BA, Stewart B, et al. UCYN-A/haptophyte symbioses dominate N2 fixation in the Southern California Current System. ISME Commun. 2021;1:1–13.Article 

Google Scholar 
Moore CM, Mills MM, Arrigo KR, Berman-Frank I, Bopp L, Boyd PW, et al. Processes and patterns of oceanic nutrient limitation. Nat Geosci. 2013;6:701–10.CAS 
Article 

Google Scholar 
Scheller S, Yu H, Chadwick GL, McGlynn SE, Orphan VJ. Artificial electron acceptors decouple archaeal methane oxidation from sulfate reduction. Science. 2016;351:703–7.CAS 
PubMed 
Article 

Google Scholar 
Pereira FC, Wasmund K, Cobankovic I, Jehmlich N, Herbold CW, Lee KS, et al. Rational design of a microbial consortium of mucosal sugar utilizers reduces Clostridiodes difficile colonization. Nat Commun. 2020;11:5104.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Mooshammer M, Kitzinger K, Schintlmeister A, Ahmerkamp S, Nielsen JL, Nielsen PH, et al. Flow-through stable isotope probing (Flow-SIP) minimizes cross-feeding in complex microbial communities. ISME J. 2021;15:348–53.CAS 
PubMed 
Article 

Google Scholar 
Drescher K, Nadell CD, Stone HA, Wingreen NS, Bassler BL. Solutions to the public goods dilemma in bacterial biofilms. Curr Biol. 2014;24:50–5.CAS 
PubMed 
Article 

Google Scholar 
Słomka J, Alcolombri U, Secchi E, Stocker R, Fernandez VI. Encounter rates between bacteria and small sinking particles. New J Phys. 2020;22:043016.Article 

Google Scholar 
Alcolombri U, Peaudecerf FJ, Fernandez VI, Behrendt L, Lee KS, Stocker R. Sinking enhances the degradation of organic particles by marine bacteria. Nat Geosci. 2021;14:775–80.CAS 
Article 

Google Scholar 
University of Massachusetts Amherst Massachusetts Lynn Margulis, Margulis L, Fester R. Symbiosis as a source of evolutionary innovation: speciation and morphogenesis. MIT Press; 1991. 454 p.Legin AA, Schintlmeister A, Sommerfeld NS, Eckhard M, Theiner S, Reipert S, et al. Nano-scale imaging of dual stable isotope labeled oxaliplatin in human colon cancer cells reveals the nucleolus as a putative node for therapeutic effect. Nanoscale Adv. 2021;3:249–62.CAS 
Article 

Google Scholar 
Schaible GA, et al. Correlative SIP-FISH-Raman-SEM-NanoSIMS links identity, morphology, biochemistry, and physiology of environmental microbes. ISME COMMUN. 2022;2:52.Article 

Google Scholar 
Yu G-H, Chi Z-L, Kappler A, Sun F-S, Liu C-Q, Teng HH, et al. Fungal nanophase particles catalyze iron transformation for oxidative stress removal and iron acquisition. Curr Biol. 2020;30:2943–50.e4.CAS 
PubMed 
Article 

Google Scholar 
Subirana MA, Riemschneider S, Hause G, Dobritzsch D, Schaumlöffel D, Herzberg M. High spatial resolution imaging of subcellular macro and trace element distribution during phagocytosis. Metallomics. 2022;14:mfac011.PubMed 
Article 

Google Scholar 
Bonnin EA, Fornasiero EF, Lange F, Turck CW, Rizzoli SO. NanoSIMS observations of mouse retinal cells reveal strict metabolic controls on nitrogen turnover. BMC Mol Cell Biol. 2021;22:5.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Jo MC, Liu W, Gu L, Dang W, Qin L. High-throughput analysis of yeast replicative aging using a microfluidic system. Proc Natl Acad Sci U S A. 2015;112:9364–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Anggraini D, Ota N, Shen Y, Tang T, Tanaka Y, Hosokawa Y, et al. Recent advances in microfluidic devices for single-cell cultivation: methods and applications. Lab Chip. 2022;22:1438–68.CAS 
PubMed 
Article 

Google Scholar 
Eriksen R, Daria V, Gluckstad J. Fully dynamic multiple-beam optical tweezers. Opt Express. 2002;10:597–602.PubMed 
Article 

Google Scholar 
Dai X, Fu W, Chi H, Mesias VSD, Zhu H, Leung CW, et al. Optical tweezers-controlled hotspot for sensitive and reproducible surface-enhanced Raman spectroscopy characterization of native protein structures. Nat Commun. 2021;12:1292.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Collins DJ, Morahan B, Garcia-Bustos J, Doerig C, Plebanski M, Neild A. Two-dimensional single-cell patterning with one cell per well driven by surface acoustic waves. Nat Commun. 2015;6:8686.CAS 
PubMed 
Article 

Google Scholar 
Hu F, Shi L, Min W. Biological imaging of chemical bonds by stimulated Raman scattering microscopy. Nat Methods. 2019;16:830–42.CAS 
PubMed 
Article 

Google Scholar 
Ge X, Pereira FC, Mitteregger M, Berry D, Zhang M, Hausmann B, et al. SRS-FISH: A high-throughput platform linking microbiome metabolism to identity at the single-cell level. Proc Natl Acad Sci U S A. 2022;119:e2203519119.PubMed 
PubMed Central 
Article 

Google Scholar 
Vandergrift GW, Kew W, Lukowski JK, Bhattacharjee A, Liyu AV, Shank EA, et al. Imaging and direct sampling capabilities of nanospray desorption electrospray ionization with absorption-mode 21 Tesla Fourier transform ion cyclotron resonance mass spectrometry. Anal Chem. 2022;94:3629–36.CAS 
PubMed 
Article 

Google Scholar 
Harrison JP, Berry D. Vibrational spectroscopy for imaging single microbial cells in complex biological samples. Front Microbiol. 2017;8:675.PubMed 
PubMed Central 
Article 

Google Scholar 
Mayali X. NanoSIMS: microscale quantification of biogeochemical activity with large-scale impacts. Ann Rev Mar Sci. 2020;12:449–67.PubMed 
Article 

Google Scholar 
Alexandrov T. Spatial metabolomics and imaging mass spectrometry in the age of artificial intelligence. Annu Rev Biomed Data Sci. 2020;3:61–87.PubMed 
PubMed Central 
Article 

Google Scholar 
Boschker HTS, Middelburg JJ. Stable isotopes and biomarkers in microbial ecology. FEMS Microbiol Ecol. 2002;40:85–95.CAS 
PubMed 
Article 

Google Scholar 
Mayali X, Weber PK, Nuccio E, Lietard J, Somoza M, Blazewicz SJ, et al. Chip-SIP: Stable Isotope Probing analyzed with rRNA-targeted microarrays and nanoSIMS. Methods Mol Biol. 2019;2046:71–87.PubMed 
Article 

Google Scholar 
Chokkathukalam A, Kim D-H, Barrett MP, Breitling R, Creek DJ. Stable isotope-labeling studies in metabolomics: new insights into structure and dynamics of metabolic networks. Bioanalysis. 2014;6:511–24.CAS 
PubMed 
Article 

Google Scholar 
Hiller K, Metallo CM, Kelleher JK, Stephanopoulos G. Nontargeted elucidation of metabolic pathways using stable-isotope tracers and mass spectrometry. Anal Chem. 2010;82:6621–8.CAS 
PubMed 
Article 

Google Scholar 
Rusconi R, Garren M, Stocker R. Microfluidics expanding the frontiers of microbial ecology. Annu Rev Biophys. 2014;43:65–91.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lee KS, Pereira FC, Palatinszky M, Behrendt L, Alcolombri U, Berry D, et al. Optofluidic Raman-activated cell sorting for targeted genome retrieval or cultivation of microbial cells with specific functions. Nat Protoc. 2021;16:634–76.CAS 
PubMed 
Article 

Google Scholar 
Wagner M, Haider S. New trends in fluorescence in situ hybridization for identification and functional analyses of microbes. Curr Opin Biotechnol. 2012;23:96–102.CAS 
PubMed 
Article 

Google Scholar  More