Philippa Kaur

Hutchinson, G. E. Concluding remarks. Cold Spring Harb. Symp. Quant. Biol. 22, 415–427. https://doi.org/10.1101/sqb.1957.022.01.039 (1957).Article 

Google Scholar 
Manenti, R., Melotto, A., Guillaume, O., Ficetola, G. F. & Lunghi, E. Switching from mesopredator to apex predator: How do responses vary in amphibians adapted to cave living?. Behav. Ecol. Sociobiol. 74, 126. https://doi.org/10.1007/s00265-020-02909-x (2020).Article 

Google Scholar 
Pekár, S., García, L. F. & Viera, C. Behaviour and Ecology of Spiders (Springer, 2017).
Google Scholar 
Nawrocki, B., Colborne, S. F., Yurkowski, D. J. & Fisk, A. T. Foraging ecology of Bowfin (Amia calva), in the Lake Huron-Erie Corridor of the Laurentian Great Lakes: Individual specialists in generalist populations. J. Great Lakes Res. 42, 1452–1460. https://doi.org/10.1016/j.jglr.2016.08.002 (2016).Article 

Google Scholar 
Nifong, J. C. Living on the edge: Trophic ecology of Alligator mississippiensis (American alligator) with access to a shallow estuarine impoundment. Bull. Fla. Mus. Nat. Hist. 54, 13–49 (2016).
Google Scholar 
Stuart, S. N. et al. Status and trends of amphibian declines and extinctions worldwide. Science 306, 1783–1786. https://doi.org/10.1126/science.1103538 (2004).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Jaeger, A. et al. Age, sex, and breeding status shape a complex foraging pattern in an extremely long-lived seabird. Ecology 95, 2324–2333 (2014).Article 
PubMed 

Google Scholar 
Salwiczek, L. H. et al. Adult cleaner wrasse outperform capuchin monkeys, chimpanzees and orang-utans in a complex foraging task derived from cleaner: Client reef fish cooperation. PLoS ONE 7, e49068. https://doi.org/10.1371/journal.pone.0049068 (2012).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Juáres, M. A., Santos, M., Mennucci, J. A., Coria, N. R. & Mariano-Jelicich, R. Diet composition and foraging habitats of Adélie and gentoo penguins in three different stages of their annual cycle. Mar. Biol. 163, 105. https://doi.org/10.1007/s00227-016-2886-y (2016).Article 
CAS 

Google Scholar 
Rosenblatt, A. E. et al. Factors affecting individual foraging specialization and temporal diet stability across the range of a large “generalist” apex predator. Oecologia 178, 5–16. https://doi.org/10.1007/s00442-014-3201-6 (2015).Article 
ADS 
PubMed 

Google Scholar 
Casper, R. M. et al. The influence of diet on foraging habitat models: A case study using nursing Antarctic fur seals. Ecography 33, 748–759. https://doi.org/10.1111/j.1600-0587.2009.06155.x (2010).Article 

Google Scholar 
Pagani-Núñez, E., Barnett, C. A., Gu, H. & Goodale, E. The need for new categorizations of dietary specialism incorporating spatio-temporal variability of individual diet specialization. J. Zool. 300, 1–7. https://doi.org/10.1111/jzo.12364 (2016).Article 

Google Scholar 
Quevedo, M., Svanbäck, R. & Eklöv, P. Intrapopulation niche partitioning in a generalist predator limits food web connectivity. Ecology 90, 2263–2274. https://doi.org/10.1890/07-1580.1 (2009).Article 
PubMed 

Google Scholar 
Ćirović, D., Penezić, A., Milenković, M. & Paunović, M. Winter diet composition of the Golden jackal (Canis aureus L. 1758) in Serbia. Mamm. Biol. 79, 132–137. https://doi.org/10.1016/j.mambio.2013.11.003 (2014).Article 

Google Scholar 
Moser, C. F., de Avila, F. R., de Oliveira, M. & Tozetti, A. M. Diet composition and trophic niche overlap between two sympatric species of Physalaemus (Anura, Leptodactylidae, Leiuperinae) in a subtemperate forest of southern Brazil. Herpeto. Notes 10, 9–15 (2017).
Google Scholar 
Lunghi, E. et al. What shapes the trophic niche of European plethodontid salamanders?. PLoS ONE 13, e0205672. https://doi.org/10.1371/journal.pone.0205672 (2018).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Evangelista, C., Boiche, A., Lecerf, A. & Cucherousset, J. Ecological opportunities and intraspecific competition alter trophic niche specialization in an opportunistic stream predator. J. Anim. Ecol. 83, 1025–1034. https://doi.org/10.1111/1365-2656.12208 (2014).Article 
PubMed 

Google Scholar 
Cloyed, C. S. & Eason, P. K. Niche partitioning and the role of intraspecific niche variation in structuring a guild of generalist anurans. R. Soc. Open Sci. 4, 170060. https://doi.org/10.1098/rsos.170060 (2017).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Dehnhard, N. et al. Is individual consistency in body mass and reproductive decisions linked to individual specialization in foraging behavior in a long-lived seabird?. Ecol. Evol. 6, 4488–4501. https://doi.org/10.1002/ece3.2213 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Jirka, K. J. & Kraft, C. E. Diet niche width and individual specialization of Brook trout in Adirondack lakes. Trans. Am. Fish Soc. 146, 716–731. https://doi.org/10.1080/00028487.2017.1290680 (2017).Article 

Google Scholar 
Reading, C. & Jofré, G. Diet composition changes correlated with body size in the Smooth snake, Coronella austriaca, inhabiting lowland heath in southern England. Amphib. Reptil. 34, 463–470. https://doi.org/10.1163/15685381-00002899 (2013).Article 

Google Scholar 
Novak, T. et al. Niche partitioning in orbweaving spiders Meta menardi and Metellina merianae (Tetragnathidae). Acta Oecol. 36, 522–529. https://doi.org/10.1016/j.actao.2010.07.005 (2010).Article 
ADS 

Google Scholar 
Stamp, N. E. & Bowers, M. D. Presence of predatory wasps and stinkbugs alters foraging behavior of cryptic and non-cryptic caterpillars on plantain (Plantago lanceolata). Oncologic 95, 376–384 (1993).ADS 

Google Scholar 
Magnusson, W. E. & Lima, A. P. The ecology of a cryptic predator, Paleosuchus tigonatus, in a tropical rainforest. J. Herpetol. 25, 41–48 (1991).Article 

Google Scholar 
Riesch, R., Tobler, M. & Plath, M. Extremophile Fishes Ecology, Evolution, and Physiology of Teleosts in Extreme Environments (Springer, 2015).
Google Scholar 
Horikoshi, K. Barophiles: Deep-sea microorganisms adapted to an extreme environment. Curr. Opin. Microbiol. 1, 291–295 (1998).Article 
CAS 
PubMed 

Google Scholar 
Mammola, S. et al. Collecting eco-evolutionary data in the dark: Impediments to subterranean research and how to overcome them. Ecol. Evol. 11, 5911–5926. https://doi.org/10.1002/ece3.7556 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Crovetto, F., Romano, A. & Salvidio, S. Comparison of two non-lethal methods for dietary studies in terrestrial salamanders. Wildl. Res. 39, 266–270. https://doi.org/10.1071/WR11103 (2012).Article 

Google Scholar 
Wake, D. B. The enigmatic history of the European. Asian and American plethodontid salamanders. Amphib-reptile 34, 323–336 (2013).Article 

Google Scholar 
Lanza, B., Pastorelli, C., Laghi, P. & Cimmaruta, R. A review of systematics, taxonomy, genetics, biogeography and natural history of the genus Speleomantes Dubois, 1984 (Amphibia Caudata Plethodontidae). Atti Mus. Civ. Stor. Nat. Trieste 52, 5–135 (2006).
Google Scholar 
Rondinini, C., Battistoni, A., Peronace, V. & Teofili, C. Lista Rossa IUCN dei Vertebrati Italiani (Comitato Italiano IUCN e Ministero dell’Ambiente e della Tutela del Territorio e del Mare, 2013).European Community. Council directive 92/43/EEC of 21 May 1992 on the conservation of natural habitats and of wild fauna and flora. J. Eur. Union 206(7), 1–44 (1992).
Google Scholar 
Salvidio, S., Palumbi, G., Romano, A. & Costa, A. Safe caves and dangerous forests? Predation risk may contribute to salamander colonization of subterranean habitats. Sci. Nat. 104, 20. https://doi.org/10.1007/s00114-017-1443-y (2017).Article 
CAS 

Google Scholar 
Lunghi, E., Manenti, R. & Ficetola, G. F. Do cave features affect underground habitat exploitation by non-troglobite species?. Acta Oecol. 55, 29–35. https://doi.org/10.1016/j.actao.2013.11.003 (2014).Article 
ADS 

Google Scholar 
Salvidio, S., Oneto, F., Ottonello, D., Costa, A. & Romano, A. Trophic specialization at the individual level in a terrestrial generalist salamander. Can. J. Zool. 93, 79–83. https://doi.org/10.1139/cjz-2014-0204 (2015).Article 

Google Scholar 
Lunghi, E. et al. Environmental suitability models predict population density, performance and body condition for microendemic salamanders. Sci. Rep. 8, 7527. https://doi.org/10.1038/s41598-018-25704-1 (2018).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ficetola, G. F. et al. Differences between microhabitat and broad-scale patterns of niche evolution in terrestrial salamanders. Sci. Rep. 8, 10575. https://doi.org/10.1038/s41598-018-28796-x (2018).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Oneto, F., Ottonello, D., Pastorino, M. V. & Salvidio, S. in Scripta Herpetologica. Studies on Amphibians and Reptiles in honour of Benedetto Lanza (eds M. Capula & C. Corti) (Edizioni Belvedere, 2014).Lunghi, E., Mascia, C., Mulargia, M. & Corti, C. Is the Sardinian grass snake (Natrix natrix cetti) an active hunter in underground environments?. Spixiana 41, 160 (2018).
Google Scholar 
Ficetola, G. F., Canedoli, C. & Stock, F. The Racovitzan impediment and the hidden biodiversity of unexplored environments. Conserv. Biol. 33, 214–216. https://doi.org/10.1111/cobi.13179 (2019).Article 
PubMed 

Google Scholar 
Lunghi, E. et al. Field-recorded data on the diet of six species of European Hydromantes cave salamanders. Sci. Data 5, 180083. https://doi.org/10.1038/sdata.2018.83 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Lunghi, E. et al. Updating salamander datasets with phenotypic and stomach content information for two mainland Speleomantes. Sci. Data 8, 150. https://doi.org/10.1038/s41597-021-00931-w (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Deban, S. M. & Dicke, U. Activation patterns of the tongue-projector muscle during feeding in the imperial cave salamander Hydromantes imperialis. J. Exp. Biol. 207, 2071–2081. https://doi.org/10.1242/jeb.00978 (2004).Article 
PubMed 

Google Scholar 
Deban, S. M., O’Reilly, J. C., Dicke, U. & van Leeuwen, J. L. Extremely high-power tongue projection in plethodontid salamanders. J. Exp. Biol. 210, 655–667. https://doi.org/10.1242/jeb.02664 (2007).Article 
PubMed 

Google Scholar 
Vignoli, L., Caldera, F. & Bologna, M. A. Trophic niche of cave populations of Speleomantes italicus. J. Nat. Hist. 40, 1841–1850 (2006).Article 

Google Scholar 
Salvidio, S. et al. Consistency in trophic strategies between populations of the Sardinian endemic salamander Speleomantes imperialis. Anim. Biol. 67, 1–16. https://doi.org/10.1163/15707563-00002517 (2017).Article 

Google Scholar 
Lunghi, E., Manenti, R. & Ficetola, G. F. Seasonal variation in microhabitat of salamanders: Environmental variation or shift of habitat selection?. PeerJ 3, e1122. https://doi.org/10.7717/peerj.1122 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Lunghi, E. et al. Thermal equilibrium and temperature differences among body regions in European plethodontid salamanders. J. Therm. Biol. 60, 79–85. https://doi.org/10.1016/j.jtherbio.2016.06.010 (2016).Article 
PubMed 

Google Scholar 
Spotila, J. R. Role of temperature and water in the ecology of lungless salamanders. Ecol. Monogr. 42, 95–125 (1972).Article 

Google Scholar 
Manenti, R., Lunghi, E. & Ficetola, G. F. Distribution of spiders in cave twilight zone depends on microclimatic features and trophic supply. Invertebr. Biol. 134, 242–251. https://doi.org/10.1111/ivb.12092 (2015).Article 

Google Scholar 
Yurkowski, D. J. et al. Latitudinal variation in ecological opportunity and intraspecific competition indicates differences in niche variability and diet specialization of Arctic marine predators. Ecol. Evol. 6, 1666–1678. https://doi.org/10.1002/ece3.1980 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Bolnick, D. I. et al. The ecology of individuals: Incidence and implications of individual specialization. Am. Nat. 161, 1–28. https://doi.org/10.1086/343878 (2003).Article 
MathSciNet 
PubMed 

Google Scholar 
Araújo, M. S., Bolnick, D. L. & Layman, C. A. The ecological causes of individual specialisation. Ecol. Lett. 14, 948–958. https://doi.org/10.1111/j.1461-0248.2011.01662.x (2011).Article 
PubMed 

Google Scholar 
Lunghi, E. et al. Same diet, different strategies: Variability of individual feeding habits across three populations of Ambrosi’s cave salamander (Hydromantes ambrosii). Diversity 12, 180. https://doi.org/10.3390/d12050180 (2020).Article 

Google Scholar 
Costa, A., Crovetto, F. & Salvidio, S. European plethodontid salamanders on the forest floor: Local abundance is related to fine-scale environmental factors. Herpetol. Conserv. Biol. 11, 344–349 (2016).
Google Scholar 
Salvidio, S., Romano, A., Oneto, F., Ottonello, D. & Michelon, R. Different season, different strategies: Feeding ecology of two syntopic forest-dwelling salamanders. Acta Oecol. 43, 42–50 (2012).Article 
ADS 

Google Scholar 
Culver, D. C. & Pipan, T. The Biology of Caves and Other Subterranean Habitats 2nd edn. (Oxford University Press, 2019).Book 

Google Scholar 
Lunghi, E., Manenti, R. & Ficetola, G. F. Cave features, seasonality and subterranean distribution of non-obligate cave dwellers. PeerJ 5, e3169. https://doi.org/10.7717/peerj.3169 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Lunghi, E., Ficetola, G. F., Zhao, Y. & Manenti, R. Are the neglected Tipuloidea crane flies (Diptera) an important component for subterranean environments?. Diversity 12, 333. https://doi.org/10.3390/d12090333 (2020).Article 

Google Scholar 
Manenti, R. et al. The stenoendemic cave-dwelling planarians (Platyhelminthes, Tricladida) of the Italian Alps and Apennines: Conservation issues. J. Nat. Conserv. 45, 90–97. https://doi.org/10.1016/j.jnc.2018.08.001 (2018).Article 

Google Scholar 
Lunghi, E. et al. Ecological observations on hybrid populations of European plethodontid salamanders, genus Speleomantes. Diversity 13, 285. https://doi.org/10.3390/d13070285 (2021).Article 

Google Scholar 
Lunghi, E., Guillaume, O., Blaimont, P. & Manenti, R. The first ecological study on the oldest allochthonous population of European cave salamanders (Hydromantes sp.). Amphib-Reptile 39, 113–119. https://doi.org/10.1163/15685381-00003137 (2018).Article 

Google Scholar 
Bolnick, D. I. et al. Ecological release from interspecific competition leads to decoupled changes in population and individual niche width. Proc. R. Soc. B 277, 1789–1797. https://doi.org/10.1098/rspb.2010.0018 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Lunghi, E. et al. Comparative reproductive biology of European cave salamanders (genus Hydromantes): Nesting selection and multiple annual breeding. Salamandra 54, 101–108 (2018).
Google Scholar 
Ficetola, G. F., Lunghi, E. & Manenti, R. Microhabitat analyses support relationships between niche breadth and range size when spatial autocorrelation is strong. Ecography 43, 724–734. https://doi.org/10.1111/ecog.04798 (2020).Article 

Google Scholar 
Lormée, H., Jouventin, P., Trouve, C. & Chastel, O. Sex-specific patterns in baseline corticosterone and body condition changes in breeding Red-footed Boobies Sula sula. Ibis 145, 212–219 (2003).Article 

Google Scholar 
Du Plessis, K. L., Martin, R. O., Hockey, P. A. R., Cunningham, S. J. & Ridley, A. R. The costs of keeping cool in a warming world: Implications of high temperatures for foraging, thermoregulation and body condition of an arid-zone bird. Glob Change Biol 18, 3063–3070. https://doi.org/10.1111/j.1365-2486.2012.02778.x (2012).Article 
ADS 

Google Scholar 
Lunghi, E. & Corti, C. Predation of European cave salamanders (Speleomantes) by the spider Meta bourneti. Spixiana 44, 54 (2021).
Google Scholar 
Lunghi, E. Doubling the lifespan of European plethodontid salamanders. Ecology 103, e03581. https://doi.org/10.1002/ecy.3581 (2022).Article 
PubMed 

Google Scholar 
Ficetola, G. F., Pennati, R. & Manenti, R. Spatial segregation among age classes in cave salamanders: Habitat selection or social interactions?. Popul Ecol 55, 217–226 (2013).Article 

Google Scholar 
Lunghi, E. et al. Interspecific and inter-population variation in individual diet specialization: Do environmental factors have a role?. Ecology 101, e03088. https://doi.org/10.1002/ecy.3088 (2020).Article 
PubMed 

Google Scholar 
Blamires, S. J. Plasticity in extended phenotypes: Orb web architectural responses to variations in prey parameters. J. Exp. Biol. 213, 3207–3212. https://doi.org/10.1242/jeb.045583 (2010).Article 
PubMed 

Google Scholar 
Costa, A. et al. Generalisation within specialization: Inter-individual diet variation in the only specialized salamander in the world. Sci. Rep. 5, 13260. https://doi.org/10.1038/srep13260 (2015).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lunghi, E. et al. Capture-mark-recapture data on the strictly protected Speleomantes italicus. Ecology 103, e3641. https://doi.org/10.1002/ecy.3641 (2022).Article 
PubMed 

Google Scholar 
Lunghi, E. & Bruni, G. Long-term reliability of visual implant elastomers in the Italian cave salamander (Hydromantes italicus). Salamandra 54, 283–286 (2018).
Google Scholar 
Lunghi, E., Bacci, F. & Zhao, Y. How can we record reliable information on animal colouration in the wild?. Diversity 13, 356. https://doi.org/10.3390/d13080356 (2021).Article 

Google Scholar 
Lunghi, E. et al. On the stability of the dorsal pattern of European cave salamanders (genus Hydromantes). Herpetozoa 32, 249–253. https://doi.org/10.3897/herpetozoa.32.e39030 (2019).Article 

Google Scholar 
Oksanen, J. et al. Vegan: Community Ecology Package. R package version 2.5-7. https://cran.r-project.org, https://github.com/vegandevs/vegan (2020).R Development Core Team. R: A language and environment for statistical computing. (R Foundation for Statistical Computing, 2021) http://www.R-project.org/.Anderson, M. J. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46. https://doi.org/10.1111/j.1442-9993.2001.01070.pp.x (2001).Article 

Google Scholar 
Băncilă, R. I., Hartel, T. R. P., Smets, J. & Cogălniceanu, D. Comparing three body condition indices in amphibians: A case study of yellow-bellied toad Bombina variegata. Amphib-Reptile 31, 558–562. https://doi.org/10.1163/017353710X518405 (2010).Article 

Google Scholar 
Labocha, M. K., Schutz, H. & Hayes, J. P. Which body condition index is best?. Oikos 123, 111–119. https://doi.org/10.1111/j.1600-0706.2013.00755.x (2014).Article 

Google Scholar 
Lunghi, E. et al. Photographic database of the European cave salamanders, genus Hydromantes. Sci. Data 7, 171. https://doi.org/10.1038/s41597-020-0513-8 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Lunghi, E., Corti, C., Manenti, R. & Ficetola, G. F. Consider species specialism when publishing datasets. Nat. Ecol. Evol. 3, 319. https://doi.org/10.1038/s41559-019-0803-8 (2019).Article 
PubMed 

Google Scholar  More

Systematic reviewWe searched relevant publications through Web of Science (all databases), Google Scholar, and the China National Knowledge Infrastructure Database between 1945 and March 2021 with the following combinations of keywords: (drain* OR lower* water table OR standing water depth OR ground water table level drawdown OR decline OR drought OR dry*) with (peatland* OR mire* OR fen OR bog OR swamp OR marsh*) with (soil respiration OR heterotrophic respiration OR microbial respiration OR soil CO2 OR soil carbon decompos* OR soil carbon minerali* or peat subsidence). Using these search terms, we initially identified 2120 different publications. To reliably evaluate WT decline impacts on SR and peat subsidence-associated soil CO2 emissions, the following further criteria were applied:1) Only paired studies with pristine peatland (i.e., undrained, near-natural peatland without direct drainage history) as a control and pristine peatland with direct WT decline (due to drainage and land use or climate-induced drying) as a treatment were included by carefully checking the descriptions of field conditions from the publications. For the pristine peatlands, we included the peatland only if the peat soil had at least 30% dry organic matter, a peat depth of >40 cm1, and did not have any direct drainage history2. We acknowledge that few, if any, untouched and completely pristine peatlands currently exist, particularly in Europe.2) WT decline in peatlands referred to only the WT depth lowered by drainage or climate-induced drying and/or additional management practices related to C or N input (e.g., manure/N fertilizers); treatments in which WT decline was combined with manipulated warming, elevated CO2, N deposition, etc., were excluded, while individual treatments (i.e., peatlands affected by WT decline without additional warming, elevated CO2, N deposition treatments, etc.) were included, as the primary objective of this study was to evaluate the responses of peatland C decomposition to WT decline.3) Each individual study included SR or at least one of its components (HR and AR), and the measurement intervals were at least monthly. The in situ measurements of SR or its components (HR and AR) covered at least the growing or nongrowing season in temperate/boreal climate zones and the whole wet or dry season in (sub)tropical climate zones.4) Both in situ and soil core/microcosm/mesocosm measurements of SR or its components (HR and AR) were included. SR and its components were exclusively measured using the chamber method. The results of the latter group were used to test the results of the former.Finally, 386 paired in situ and 21 paired soil core incubation measurements of SR or its components (HR and AR) were extracted from 63 in situ studies and 9 soil core studies, respectively (see Supplementary Data A). Furthermore, to estimate HR emissions from global drained peatlands, the in situ measured paired peat subsidence rate (Rps, cm yr–1) and drainage duration (i.e., years since first drainage) and the proportion of peat subsidence rate attributed to oxidation (Po, %) and drainage duration, as well as the soil (0–30 cm) organic C and bulk density in pristine peatlands, were extracted from peer-reviewed publications. In drained boreal and temperate peatlands, most studies measured the total subsidence (in meter) during a certain drainage period, therefore the average Rps was calculated as the ratio of total subsidence and drainage years. It was assumed that the Rps was faster at the beginning and lower at the end of drainage duration, so the average subsidence rate is the rate for the middle year of the drainage duration41. The remaining studies directly showed the in situ measured Rps at the ith year of drainage. A similar procedure was applied for the Po in the ith year of drainage. In sum, 230 paired Rps–drainage duration observations and 49 paired Po–drainage duration observations, as well as 76 SOC and 63 BD in pristine peatlands, were taken from 80, 25, 58, and 44 studies, respectively (see Supplementary Data B).Data compilationTo systematically evaluate the impacts of WT decline on SR in pristine peatlands and clarify the underlying mechanisms, we obtained data related to SR and its components (HR and AR) together with environmental variables such as the mean annual temperature [MAT], mean annual precipitation [MAP], peat depth [PD], WT depth [WTD], soil water content [SWC], soil temperature [Ts], soil redox potential [Eh], soil air oxygen level [O2], soil bulk density [BD], soil pH [pH], soil organic carbon [SOC], soil total nitrogen [TN], soil total phosphorus [TP], soil ammonium [({{{{rm{NH}}}}}_{4}^{+})], soil nitrate [({{{{rm{NO}}}}}_{3}^{-})], soil dissolved organic carbon [DOC], microbial biomass carbon [MBC], microbial biomass nitrogen [MBN], dissolved total phosphorus [DTP], belowground biomass [BGB], iron [Fe3+, Fe2+] and sulfate [({{{{rm{SO}}}}}_{4}^{2-})] when possible. If available, other important information, such as geographic location (latitude, longitude), climate and WT decline driver and duration, intensity, peatland type, Rps, Po, nutrient type, inundated condition, microtopography, and plant functional types, was recorded. For WT decline intensity, net WT declines greater and less than 30 cm were defined as deep and shallow declines, respectively, according to the IPCC wetland report42. The abovementioned information about pristine peatlands and peatlands affected by WT decline is compiled in Supplementary Data A and B.We subsequently extracted the mean ((bar{X})), standard deviation (SD) and replicates (n) from different publications. If studies reported standard error (SE) rather than SD, then SD was calculated by SE (sqrt{n}). If studies reported only the median, maximum, minimum, and 25th and 75th percentiles, then the mean and SD were estimated following the mathematical equations recommended by ref. 60. If neither SD nor SE was reported, then the missing SD was estimated by multiplying the reported mean by the average coefficient of variation (CV) obtained from the remaining observations, resulting in both the mean and SD being reported61. The data were either obtained directly from tables and texts or extracted by digitizing graphs using Getdata Graph Digitizer software (version 2.26, Russia).The final database consisted of 250 paired SR, 101 paired HR and 35 paired AR in situ observations. Only 35 paired observations simultaneously reported SR, HR, and AR. Twenty-one paired SR soil core incubation measurements were also collected to test the results of the in situ measurements. The dataset mainly originated from Europe, North America, and Southeast Asia, and most studies ( >70%) were conducted in temperate and boreal peatlands in the Northern Hemisphere (Fig. 1a). Moreover, 230 paired Rps–drainage duration observations and 49 paired Po–drainage duration observations (Fig. 5a, b) and an additional 485 drainage year (Supplementary Fig. 9) observations classified by climate zone (i.e., boreal, temperate and tropical) and land use (i.e., agriculture, forestry, and grassland) were collected. A total of 76 SOC and 63 BD measurements from pristine peatlands categorized by climate zone (i.e., boreal, temperate, and tropical) were extracted to estimate Rps by oxidation and associated soil HR from global pristine peatlands due to drainage activities (Supplementary Fig. 10 and Supplementary Data B). In this study, we were unable to estimate climate drying-induced net CO2 emissions through soil HR, as the areas of pristine peatlands affected by climate drying currently remain unknown.Meta-analysisTo assess the relative changes in SR and its components (HR and AR), as well as environmental variables (e.g., SOC, BD, Ts, etc.) due to WT decline, the log-transformed response ratio (RR) was used:62$${{{mathrm{ln}}}}({{{rm{RR}}}})=,{{{mathrm{ln}}}}({X}_{{{{rm{t}}}}}/{X}_{{{{rm{c}}}}})$$
(1)
The results are presented as the percent change ((elnRR  – 1) × 100). The variance (v) of RR was estimated using the following equation:$$v=frac{{{{{rm{SD}}}}}_{{{{rm{t}}}}}^{2}}{{n}_{{{{rm{t}}}}},{X}_{{{{rm{t}}}}}^{2}}+frac{{{{{rm{SD}}}}}_{{{{rm{c}}}}}^{2}}{{n}_{{{{rm{c}}}}},{X}_{{{{rm{c}}}}}^{2}}$$
(2)
where Xt and Xc indicate the means of the treatment and control, SDt and SDc indicate the SDs of the treatment and control and nt and nc indicate the numbers of replicates in the treatment and control, respectively.However, in our study, approximately 60% of the WTD and Eh observations for the peatlands in pristine condition (control) and affected by WT decline (treatment) showed opposite signs; e.g., the pristine peatlands generally exhibited positive WTDs (higher than the peat surface) and negative Eh values, while those affected by WT decline exhibited negative WTDs (lower than the peat surface) and positive Eh values. Since it is impossible to calculate the logarithm of negative values, we introduced a new study index (net changes) for these two variables in our meta-analysis according to ref. 63:$$D={X}_{{{{rm{t}}}}}-{X}_{{{{rm{c}}}}}$$
(3)
where Xt and Xc indicate the paired annual mean WTD and Eh for the treatment and control, respectively, and D indicates the difference between the treatment and control.The SD and variance (v) of D were estimated using the following equation:$${{{rm{SD}}}}=sqrt{frac{({n}_{{{{rm{c}}}}}-1);{{{{rm{SD}}}}}_{{{{rm{c}}}}}^{2}+({n}_{{{{rm{t}}}}}-1);{{{{rm{SD}}}}}_{{{{rm{t}}}}}^{2}}{{n}_{{{{rm{c}}}}}+{n}_{{{{rm{t}}}}}-2}}$$
(4)
$$v=frac{{{{{rm{SD}}}}}_{{{{rm{t}}}}}^{2}}{{n}_{{{{rm{t}}}}}}+frac{{{{{rm{SD}}}}}_{{{{rm{c}}}}}^{2}}{{n}_{{{{rm{c}}}}}}$$
(5)
where SDt and SDc indicate the SD of the treatment and control and nt and nc indicate the number of replicates for the treatment and control, respectively.The weighted mean RR or D was calculated by individual RR or D with bias-corrected 95% confidence intervals (CIs) using the rma.mv function in the metafor package in R software (R core team, 2019)64, in which the variable “study” was regarded as a random effect to account for the dependence of observations derived from the same study. The impact of WT decline on a response variable was considered significant if the 95% CI did not overlap 065. Differences between subgroups (e.g., WT decline driver, climate zone, drainage duration) were considered significant if the 95% CIs did not overlap each other65. The frequency distribution of RR was calculated to test variability among individual studies using the Gaussian function (i.e., normal distribution)66.Estimation of peat subsidence rate by oxidation and associated HR rateDrainage has induced widespread peat subsidence and associated large CO2 release through soil HR and consequently reduced the sustainable utilization of drained peatlands and contributed to global warming11,12. In this study, we estimated the spatial patterns of Rps by oxidation and associated soil HR from global drained peatlands. Using the 230 paired Rps and drainage duration observations synthesized in this study, we first constructed empirical models between Rps and drainage duration for drained peatlands categorized by climate zone (boreal, temperate and tropical climate) and land use (i.e., agriculture, forestry and grassland) (Fig. 5a, b). The values of Rps for certain groups classified by climate zone and land use could be estimated by using the corresponding empirical models established in this study and reported drainage durations that were extracted from the literature. The empirical models categorized by climate zone and land use are listed below (Fig. 5a, b):$${R}_{{{{rm{ps}}}}}{mbox{-}}{{{rm{Bor}}}}{mbox{-}}{{{rm{Tem}}}}{mbox{-}}{{{rm{Agr}}}}=13.95,{{{{rm{Dur}}}}}^{-0.58},,n=48,,{R}_{{{{rm{adj}}}}.}^{2}=0.85,,p; < ; 0.0001$$ (6) $${R}_{{{{rm{ps}}}}}{mbox{-}}{{{rm{Bor}}}}{mbox{-}}{{{rm{Tem}}}}{mbox{-}}{{{rm{For}}}}=5.36,{{{{rm{Dur}}}}}^{-0.83},,n=21,,{R}_{{{{rm{adj}}}}.}^{2}=0.92,,p; < ; 0.0001$$ (7) $${R}_{{{{rm{ps}}}}}{mbox{-}}{{{rm{Bor}}}}{mbox{-}}{{{rm{Tem}}}}{mbox{-}}{{{rm{Gra}}}}=5.55,{{{{rm{Dur}}}}}^{-0.36},,n=40,,{R}_{{{{rm{adj}}}}.}^{2}=0.61,,p; < ; 0.0001$$ (8) $${R}_{{{{rm{ps}}}}}{mbox{-}}{{{rm{Tro}}}}{mbox{-}}{{{rm{Agr}}}}{mbox{-}}{{{rm{For}}}}{mbox{-}}{{{rm{Gra}}}}=6.63,{{{{rm{Dur}}}}}^{-0.37},,n=121,,{R}_{{{{rm{adj}}}}.}^{2}=0.55,,p; < ; 0.0001$$ (9) where Rps indicates the peat subsidence rate (cm yr–1), Dur is the drainage duration, and the numbers indicate coefficients for the established empirical models. Bor, Tem, and Tro indicate boreal, temperate, and tropical climate zones, respectively. Agr, For, and Gra represent agriculture, forestry, and grassland land uses, respectively. We note that it was not possible to further distinguish these models between boreal and temperate climate zones and among agriculture, forestry, or grassland land use in tropical climates, as there is currently a lack of sufficient measurements, which warrants more research.However, the Rps is triggered by a combination of processes such as physical compaction by heavy equipment or livestock trampling and shrinkage through contraction of organic fibers when drying, consolidation by loss of water from pores in the peat and oxidation owing to the breakdown of peat organic matter10,11,12. Therefore, to reliably estimate the soil HR rate from Rps due to oxidation, the proportion of Rps attributed to oxidation (Po, in %) should be considered12. Using the 49 paired Po and drainage duration observations synthesized in this study, we then constructed empirical models between Po and drainage duration for drained peatlands that were also categorized by climate zone (boreal, temperate, and tropical climate) and land use (agriculture, forestry, and grassland) (Fig. 5c, d). Similarly, the Po values of certain groups classified by climate zone and land use could be estimated by using the corresponding empirical models established in this study and reported drainage durations that were extracted from the literature. The empirical models categorized by climate zone and land use are shown below (Fig. 5c, d):$$ {P}_{{{{rm{o}}}}}{mbox{-}}{{{rm{Tem}}}}{mbox{-}}{{{rm{Bor}}}}{mbox{-}}{{{rm{Agr}}}}{mbox{-}}{{{rm{For}}}}{mbox{-}}{{{rm{Gra}}}}=12.05,{{{mathrm{Ln}}}}({{{rm{Dur}}}})+2.15,,n=30,\ {R}_{{{{rm{adj}}}}.}^{2}=0.89,,p; < ;0.0001$$ (10) $$ {P}_{{{{rm{o}}}}}{mbox{-}}{{{rm{Tro}}}}{mbox{-}}{{{rm{Agr}}}}{mbox{-}}{{{rm{For}}}}{mbox{-}}{{{rm{Gra}}}}=14.36,{{{mathrm{Ln}}}}({{{rm{Dur}}}})+37.05,,n=19,\ {R}_{{{{rm{adj}}}}.}^{2}=0.81,,p; < ;0.0001$$ (11) where Po indicates the proportion of Rps attributable to oxidation, Dur is the drainage duration, and the numbers indicate coefficients for the established empirical models. The abbreviations Bor, Tem, Tro, Agr, For, and Gra have been described previously. We note that the different land uses shared the same models across temperate and boreal climates and tropical climate due to a lack of sufficient global observations. This will also induce some uncertainties in our analysis.Furthermore, the soil HR (FHR, Mt C yr−1) due to peat oxidation induced by drainage was estimated using the following equation according to ref. 11:$${F}_{{{{rm{HR}}}}}=sum {R}_{{{{rm{ps}}}},i,j}times {P}_{{{{rm{o}}}},i,j}times {{{{rm{SOC}}}}}_{i}times {{{{rm{BD}}}}}_{i}times {A}_{i,j}$$ (12) where SOC (g kg–1) and BD (g cm–3) indicate the soil (0–30 cm) organic C concentration and bulk density of pristine peatlands, respectively; A (×103 km2) indicates the drained peatland area; i indicates the climate zone (boreal, temperate or tropical); j indicates the land use (agriculture, forestry or grassland); and Rps (cm yr–1) and Po (%) are described in Eqs. (6–11). Datasets of the SOC concentration and BD and Rps due to oxidation were systematically reviewed and bootstrapped and categorized by climate zones and land uses (see Supplementary Fig. 10 and Supplementary Data B). Regarding the large uncertainties for areas of drained peatlands, we combined two previously published datasets (72, 61, 22, 37, 43, 26, 94, 109, and 39 × 103 km2 by ref. 18, and 37, 55, 4, 109, 63, 58, 96, 72, and 1 × 103 km2 by ref. 20. for agriculture-, forestry- and grassland-drained peatlands in boreal, temperate and tropical climate zones, respectively) and obtained their mean values with 95% CIs (for details, see bootstrapping procedure in Data analysis). Uncertainties (i.e., 95% CI) in total HR (δFHR) were propagated according to the Gaussian random error propagation principle as follows:$${{{rm{delta }}}}{F}_{{{{rm{HR}}}}}=sqrt{sum sqrt{begin{array}{c}{(delta {R}_{{{{rm{ps}}}},i,j})}^{2}times {({P}_{{{{rm{o}}}},i,j}times {{{{rm{SOC}}}}}_{i}times {{{{rm{BD}}}}}_{i}times {A}_{i,j})}^{2}+\ {(delta {P}_{{{{rm{o}}}},i,j})}^{2}times {({R}_{{{{rm{ps}}}},i,j}times {{{{rm{SOC}}}}}_{i}times {{{{rm{BD}}}}}_{i}times {A}_{i,j})}^{2}+\ {(delta {{{{rm{SOC}}}}}_{i})}^{2}times {({R}_{{{{rm{ps}}}},i,j}times {P}_{{{{rm{o}}}},i,j}times {{{{rm{BD}}}}}_{i}times {A}_{i,j})}^{2}+\ {(delta {{{{rm{BD}}}}}_{i})}^{2}times {({R}_{{{{rm{ps}}}},i,j}times {P}_{{{{rm{o}}}},i,j}times {{{{rm{SOC}}}}}_{i}times {A}_{i,j})}^{2}+\ {(delta {A}_{i,j})}^{2}times {({R}_{{{{rm{ps}}}},i,j}times {P}_{{{{rm{o}}}},i,j}times {{{{rm{SOC}}}}}_{i}times {{{{rm{BD}}}}}_{i})}^{2}end{array}}}$$ (13) where δFHR, δRps, δPo, δSOC, δBD, and δA indicate the 95% CIs of total soil HR, Rps, Po, SOC, and BD and drained peatland area, respectively, and i and j indicate the climate zone (boreal, temperate, tropical) and land use (agriculture, forestry, or grassland), respectively.To further estimate the total SR (FSR, Mt C yr−1) and its uncertainty (δFSR) from global drained peatlands, the following equations were used:$${F}_{{{{rm{SR}}}}}=sum frac{{F}_{{{{rm{HR}}}},i,j}}{{C}_{{{{rm{HR}}}},i,j}}$$ (14) $$delta {F}_{{{{rm{SR}}}}}=sqrt{sum sqrt{{(frac{1}{{C}_{{{{rm{HR}}}},i,j}})}^{2}times delta {F}_{{{{rm{HR}}}},i,j}^{2}+{(-frac{{F}_{{{{rm{HR}}}},i,j}}{{C}_{{{{rm{HR}}}},i,j}^{2}})}^{2}times delta {C}_{{{{rm{HR}}}},i,j}^{2}}}$$ (15) where CHR (%) indicates the mean relative contribution of HR to SR from simultaneously measured SR, HR, and AR from our meta-analysis (see Supplementary Fig. 11 and Supplementary Data A) and i and j indicate the climate zone (boreal, temperate, tropical) and land use (agriculture, forestry, or grassland), respectively. FHR and δFHR are given in Eqs. (12, 13). We note that the CHR could be classified only by climate zone, as there is a lack of sufficient measurements of land use; that is, the different land uses under the same climate shared the same CHR value, which may induce uncertainties in estimating the total SR from global drained peatlands.Regarding the abovementioned lack of sufficient measurements for distinguishing between boreal and temperate drained peatlands, we also used another method to estimate the annual total HR and SR from global drained peatlands. Specifically, we obtained the mean values of Rps by oxidation across boreal and temperate drained peatlands for each land use (i.e., climate zones were classified as boreal+temperate or tropical) (Supplementary Fig. 14 and Supplementary Table 1). The estimation process was the same as that previously described. The different estimation methods were likely to provide results with greater convergence.Data analysisSignificant differences in observed variables were tested by performing nonparametric analysis. Specifically, tests with two independent samples (i.e., Mann–Whitney U test) were used for only two variables (i.e., to compare the contribution of HR to SR between pristine and drained peatlands), and tests with two or more independent samples (i.e., Kruskal–Wallis test and pairwise comparisons) were used if there were three or more variables (i.e., SOC, BD and Rps due to oxidation in the boreal, temperate and tropical climate zones or different land uses). Linear or nonlinear regression analysis was performed to examine the relationships between the responses of SR and its components with environmental variables or the peat subsidence rate with drainage duration.To reliably estimate the uncertainties in Rps by oxidation, SOC, BD, drained peatland area, and relative contribution of HR to SR, bootstrap resampling with 10000 iterations was conducted using the boot package, and 95% CIs were calculated using the “basic” type. The ggplot 2 package in R software (R core team, 2019) was used for statistical analysis. Data were expressed as the means with their 95% CIs, and significance of the regression analyses was indicated at the level of p  More

Brown, M. E. et al. In Climate Change, global food security, and the U.S. food system (2015).Masson-Delmotte, V. et al. AR6 Climate Change 2021: The Physical Science Basis—IPCC. In Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (2021).Douris, J. et al. WMO Atlas of Mortality and Economic Losses from Weather, Climate and Water Extremes (1970–2019) (WMO-No. 1267). In WMO Statement on the state of the Global Climate vol. 1267 (WMO, 2021).Smith, A. B. U.S. Billion-dollar Weather and Climate Disasters, 1980–present (NCEI Accession 0209268). In National Centers for Environmental Information (2020).Lobell, D. B. et al. The critical role of extreme heat for maize production in the United States. Nat. Clim. Change 3, 497–501 (2013).Article 
ADS 

Google Scholar 
Mann, M. E. et al. Projected changes in persistent extreme summer weather events: The role of quasi-resonant amplification. Sci. Adv. 4, 5 (2018).Article 

Google Scholar 
Li, Y., Guan, K., Schnitkey, G. D., DeLucia, E. & Peng, B. Excessive rainfall leads to maize yield loss of a comparable magnitude to extreme drought in the United States. Glob. Chang. Biol. 25, 2325–2337 (2019).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Daloz, A. S. et al. Direct and indirect impacts of climate change on wheat yield in the Indo-Gangetic plain in India. J. Agric. Food Res. 4, 100–132 (2021).
Google Scholar 
Leng, G. Maize yield loss risk under droughts in observations and crop models in the Unites States. Environ. Res. Lett. 16, 24016 (2021).Article 

Google Scholar 
Backlund, P., Janetos, A., & Schimel, D. In The Effects of Climate Change on Agriculture, Land Resources, Water Resources, and Biodiversity in the United States Synthesis and Assessment Product, vol. 4.3 (2008).Scheelbeek, P. F. D., Tuomisto, H. L., Bird, F. A., Haines, A. & Dangour, A. D. Effect of environmental change on yield and quality of fruits and vegetables: Two systematic reviews and projections of possible health effects. Lancet Glob. Health 5, S21 (2017).Article 

Google Scholar 
Drewnowski, A., Dwyer, J., King, J. C. & Weaver, C. M. A proposed nutrient density score that includes food groups and nutrients to better align with dietary guidance. Nutr. Rev. 77, 404–416 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Weinberger, K. & Lumpkin, T. A. Diversification into horticulture and poverty reduction: A research agenda. World Dev. 35, 1464–1480 (2007).Article 

Google Scholar 
Barnabás, B., Jäger, K. & Fehér, A. The effect of drought and heat stress on reproductive processes in cereals. Plant Cell Environ. 31, 11–38 (2008).PubMed 

Google Scholar 
Kazan, K. & Lyons, R. The link between flowering time and stress tolerance. J. Exp. Bot. 67, 47–60 (2016).Article 
CAS 
PubMed 

Google Scholar 
Messina, C. D. et al. On the dynamic determinants of reproductive failure under drought in maize. In Silico Plants 1, 1–14 (2019).Article 

Google Scholar 
Yang, X., Wang, B., Chen, L., Li, P. & Cao, C. The different influences of drought stress at the flowering stage on rice physiological traits, grain yield, and quality. Sci. Rep. 9, 3742 (2019).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Deryng, D., Conway, D., Ramankutty, N., Price, J. & Warren, R. Global crop yield response to extreme heat stress under multiple climate change futures. Environ. Res. Lett. 9, 034011 (2014).Article 
ADS 

Google Scholar 
Owen, P. C. Responses of a semi-dwarf wheat to temperatures representing a tropical dry season. II. Extreme temperatures. Exp. Agric. 7, 43–47 (1971).Article 

Google Scholar 
Liu, F., Jensen, C. R. & Andersen, M. N. A review of drought adaptation in crop plants: Changes in vegetative and reproductive physiology induced by ABA-based chemical signals. Aust. J. Agric. Res. 56, 1245–1252 (2005).Article 
CAS 

Google Scholar 
Turc, O., Bouteillé, M., Fuad-Hassan, A., Welcker, C. & Tardieu, F. The growth of vegetative and reproductive structures (leaves and silks) respond similarly to hydraulic cues in maize. New Phytol. 212, 377–388 (2016).Article 
PubMed 

Google Scholar 
Fuad-Hassan, A., Tardieu, F. & Turc, O. Drought-induced changes in anthesis-silking interval are related to silk expansion: A spatio-temporal growth analysis in maize plants subjected to soil water deficit. Plant Cell Environ. 31, 1349–1360 (2008).Article 
PubMed 

Google Scholar 
USDA–National Agricultural Statistics Service (2021). https://data.nal.usda.gov/dataset/nass-quick-stats, accessed 29 December 2021.Challinor, A. J., Parkes, B. & Ramirez-Villegas, J. Crop yield response to climate change varies with cropping intensity. Glob. Chang. Biol. 21, 1679–1688 (2015).Article 
ADS 
PubMed 

Google Scholar 
Kukal, M. S. & Irmak, S. Climate-driven crop yield and yield variability and climate change impacts on the U.S. great plains agricultural production. Sci. Rep. 8, 1–18 (2018).Article 
ADS 

Google Scholar 
Iizumi, T. & Sakai, T. The global dataset of historical yields for major crops 1981–2016. Sci. Data 7, 1–7 (2020).Article 

Google Scholar 
Thornton, M. M. et al. In Daymet: Daily surface weather data on a 1-km grid for North America, Version 4. ORNL DAAC (2020).Ritchie, S. W., Hanway, J. J., Benson, G. O., & Herman, J. C. How a corn plant develops: Special report no. 48. In Ames: Iowa State University of Science and Technology Cooperative Extension Service (1986).Nicholls, N. Increased Australian wheat yield due to recent climate trends. Nature 387, 484–485 (1997).Article 
ADS 
CAS 

Google Scholar 
Schlenker, W. & Roberts, M. J. Nonlinear temperature effects indicate severe damages to US crop yields under climate change. Proc. Natl. Acad. Sci. USA 106, 15594–15598 (2009).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Dhaliwal, D. S. & Williams, M. M. I. I. Understanding variability in optimum plant density and recommendation domains for crowding stress tolerant processing sweet corn. PLoS ONE 15, e0228809 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Siirila-Woodburn, E. R. et al. A low-to-no snow future and its impacts on water resources in the western United States. Nat. Rev. Earth Environ. 2, 800–819 (2021).Article 
ADS 

Google Scholar 
Gilmore, E. C. & Rogers, J. S. heat units as a method of measuring maturity in corn. Agron. J. 50, 611–615 (1958).Article 

Google Scholar 
Wang, J. Y. A critique of the heat unit approach to plant response studies. Ecology 41, 785–790 (1960).Article 

Google Scholar 
Cross, H. Z. & Zuber, M. S. Prediction of flowering dates in maize based on different methods of estimating thermal units. Agron. J. 64, 351–355 (1972).Article 

Google Scholar 
Lobell, D. B., Bänziger, M., Magorokosho, C. & Vivek, B. Nonlinear heat effects on African maize as evidenced by historical yield trials. Nat. Clim. Change 1, 42–45 (2011).Article 
ADS 

Google Scholar 
Díaz, E. L. et al. In Chapter 20: US Caribbean. Impacts, Risks, and Adaptation in the United States: The Fourth National Climate Assessment, Volume II,(2018).Wang, Y. et al. Flowering dynamics, pollen, and pistil contribution to grain yield in response to high temperature during maize flowering. Environ. Exp. Bot. 158, 80–88 (2019).Article 

Google Scholar 
Lohani, N., Singh, M. B. & Bhalla, P. L. High temperature susceptibility of sexual reproduction in crop plants. J. Exp. Bot. 71, 555–568 (2020).Article 
CAS 
PubMed 

Google Scholar 
Jagadish, S. V. K., Craufurd, P. Q. & Wheeler, T. R. High temperature stress and spikelet fertility in rice (Oryza sativa L.). J. Exp. Bot. 58, 1627–1635 (2007).Article 
CAS 
PubMed 

Google Scholar 
Gourdji, S. M., Sibley, A. M. & Lobell, D. B. Global crop exposure to critical high temperatures in the reproductive period: Historical trends and future projections. Environ. Res. Lett. 8, 024041 (2013).Article 
ADS 

Google Scholar 
Hedhly, A., Hormaza, J. I. & Herrero, M. Global warming and sexual plant reproduction. Trends Plant Sci. 14, 30–36 (2008).Article 
PubMed 

Google Scholar 
Zhao, F. et al. Effects of heat stress during grain filling on sugar accumulation and enzyme activity associated with sucrose metabolism in sweet corn. Acta Agron. Sin. 39, 1644–1651 (2013).Article 
CAS 

Google Scholar 
Lobell, D. B., Bonfils, C. J., Kueppers, L. M. & Snyder, M. A. Irrigation cooling effect on temperature and heat index extremes. Geophys. Res. Lett. 35, 9705 (2008).Article 
ADS 

Google Scholar 
Thiery, W. et al. Present-day irrigation mitigates heat extremes. J. Geophys. Res. Atmos. 122, 1403–1422 (2017).Article 
ADS 

Google Scholar 
Li, Y. et al. Quantifying irrigation cooling benefits to maize yield in the US Midwest. Glob. Chang. Biol. 26, 3065–3078 (2020).Article 
ADS 
PubMed 

Google Scholar  More

Adams, D. C., Rohlf, F. J. & Slice, D. E. A field comes of age: Geometric morphometrics in the 21st century. Hystrix 24, 7–14. https://doi.org/10.4404/hystrix-24.1-6283 (2013).Article 

Google Scholar 
Terray, L. et al. Skull morphological evolution in Malagasy endemic Nesomyinae rodents. PLoS ONE 17, e0263045. https://doi.org/10.1371/journal.pone.0263045 (2022).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Viacava, P., Baker, A. M., Blomberg, S. P., Phillips, M. J. & Weisbecker, V. Using 3D geometric morphometrics to aid taxonomic and ecological understanding of a recent speciation event within a small Australian marsupial (Antechinus: Dasyuridae). Zool. J. Linn. Soc. 1–16. https://doi.org/10.1093/zoolinnean/zlab048 (2021).Brassard, C. et al. Morphological and functional divergence of the lower jaw between native and invasive red foxes. J. Mamm. Evol. 29, 335–352. https://doi.org/10.1007/s10914-021-09593-2 (2022).Article 

Google Scholar 
Boessneck, J. & von den Driesch, A. The significance of measuring animal bones from archaeological sites. In Approaches to Faunal Analysis in the Middle East (eds Meadows, R. H. & Zeder, M. A.) 5–39 (Peabody Museum Bulletin 2, 1978).
Google Scholar 
Serjeantson, D. ‘Science is measurement’; ABMAP, a database of domestic animal bone measurements. Environ. Archaeol. 10, 97–103. https://doi.org/10.1179/env.2005.10.1.97 (2005).Article 

Google Scholar 
Haruda, A. F. Separating sheep (Ovis aries L.) and goats (Capra hircus L.) using geometric morphometric methods: An investigation of astragalus morphology from late and final Bronze Age Central Asian contexts. Int. J. Osteoarchaeol. 27, 551–562 (2017).Article 

Google Scholar 
Davis, S. J. M. Towards a metrical distinction between sheep and goat astragali. In Economic Zooarchaeology: Studies in Hunting, Herding and Early Agriculture (eds Rowley-Conwy, P. et al.) 93–138 (Oxbow Books Limited, 2019).
Google Scholar 
Jeanjean, M. et al. Sorting the flock: Quantitative identification of sheep and goat from isolated third lower molars and mandibles through geometric morphometrics. J. Archaeol. Sci. 141, 105580. https://doi.org/10.1016/j.jas.2022.105580 (2022).Article 

Google Scholar 
Evin, A. et al. Phenotype and animal domestication: A study of dental variation between domestic, wild, captive, hybrid and insular Sus scrofa. BMC Evol. Biol. 15, 1–16. https://doi.org/10.1186/s12862-014-0269-x (2015).Article 

Google Scholar 
Harbers, H. et al. The mark of captivity: Plastic responses in the ankle bone of a wild ungulate (Sus scrofa). R. Soc. Open Sci. 7, 192039. https://doi.org/10.1098/rsos.192039 (2020).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Drake, A. G., Coquerelle, M. & Colombeau, G. 3D morphometric analysis of fossil canid skulls contradicts the suggested domestication of dogs during the late Paleolithic. Sci. Rep. 5, 8299. https://doi.org/10.1038/srep08299 (2015).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ventresca Miller, A. R., Haruda, A., Varfolomeev, V., Goryachev, A. & Makarewicz, C. A. Close management of sheep in ancient Central Asia: Evidence for foddering, transhumance, and extended lambing seasons during the Bronze and Iron Ages. Sci. Technol. Archaeol. Res. 6, 41–60. https://doi.org/10.1080/20548923.2020.1759316 (2020).Duval, C., Lepetz, S., Horard-Herbin, M.-P. & Cucchi, T. Did Romanization impact Gallic pig morphology? New insights from molar geometric morphometrics. J. Archaeol. Sci. 57, 345–354. https://doi.org/10.1016/j.jas.2015.03.004 (2015).Article 

Google Scholar 
Davis, S. J. M. Zooarchaeological evidence for Moslem and Christian improvements of sheep and cattle in Portugal. J. Archaeol. Sci. 35, 991–1010. https://doi.org/10.1016/j.jas.2007.07.001 (2008).Article 

Google Scholar 
Samper Carro, S. C., Louys, J. & Oonnor, S. Shape does matter: A geometric morphometric approach to shape variation in Indo-Pacific fish vertebrae for habitat identification. J. Archaeol. Sci. 99, 124–134. https://doi.org/10.1016/j.jas.2018.09.010 (2018).Stimpson, C. M. A 48,000 year record of swiftlets (Aves: Apodidae) in North-western Borneo: Morphometric identifications and palaeoenvironmental implications. Palaeogeogr. Palaeoclimatol. Palaeoecol. 374, 132–143. https://doi.org/10.1016/j.palaeo.2013.01.011 (2013).Article 

Google Scholar 
Medina, M. E., Picasso, M. B. J., Campos, M. R. & Avila, N. C. Tarsometatarsus, eggshells, and the species level identification of large-sized flightless birds from Boyo Paso 2 (Sierras of Córdoba, Argentina). Int. J. Osteoarchaeol. 29, 584–594. https://doi.org/10.1002/oa.2754 (2019).Article 

Google Scholar 
Weaver, L. N. & Grossnickle, D. M. Functional diversity of small-mammal postcrania is linked to both substrate preference and body size. Curr. Zool. 66, 539–553. https://doi.org/10.1093/cz/zoaa057 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Chen, X., Milne, N. & O’Higgins, P. Morphological variation of the thoracolumbar vertebrae in Macropodidae and its functional relevance. J. Morphol. 266, 167–181. https://doi.org/10.1002/jmor.10370 (2005).Article 
PubMed 

Google Scholar 
Etienne, C., Filippo, A., Cornette, R. & Houssaye, A. Effect of mass and habitat on the shape of limb long bones: A morpho-functional investigation on Bovidae ( Mammalia: Cetartiodactyla ). J. Anat. 238, 886–904. https://doi.org/10.1111/joa.13359 (2020).Article 
PubMed 

Google Scholar 
Bassarova, M., Janis, C. M. & Archer, M. The calcaneum-on the heels of marsupial locomotion. J. Mamm. Evol. 16, 1–23. https://doi.org/10.1007/s10914-008-9093-7 (2009).Article 

Google Scholar 
Janis, C. M., Buttrill, K. & Figueirido, B. Locomotion in extinct giant kangaroos: Were Sthenurines hop-less monsters?. PLoS ONE 9, e109888. https://doi.org/10.1371/journal.pone.0109888 (2014).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Argot, C. Functional-adaptive analysis of the hindlimb anatomy of extant marsupials and the paleobiology of the paleocene marsupials Mayulestes ferox and Pucadelphys andinus. J. Morphol. 253, 76–108. https://doi.org/10.1002/jmor.1114 (2002).Article 
PubMed 

Google Scholar 
Mein, E. & Manne, T. Identifying marsupials from Australian archaeological sites: Current methodological challenges and opportunities in zooarchaeological practice. Archaeol. Ocean. 56, 133–141. https://doi.org/10.1002/arco.5234 (2021).Article 

Google Scholar 
Woinarski, J. C. Z. et al. Reading the black book: The number, timing, distribution and causes of listed extinctions in Australia. Biol. Conserv. 239, 108261. https://doi.org/10.1016/j.biocon.2019.108261 (2019).Article 

Google Scholar 
Garvey, J. Preliminary zooarchaeological interpretations from Kutikina Cave, south-west Tasmania. Aust. Aborig. Stud. 1, 57–62 (2006).
Google Scholar 
Veth, P. et al. Montebello Islands Archaeology: Late Quaternary Foragers on an Arid Coastline. (BAR Publishing, 2007).Morse, K. Who can see the sea? Prehistoric Aboriginal occupation of the Cape Range peninsula. Rec. West. Aust. Mus. Suppl. 45, 227–242 (1993).
Google Scholar 
Warburton, N. M. & Prideaux, G. Functional pedal morphology of the extinct tree-kangaroo Bohra (Diprotodontia: Macropodidae). In Macropods: The Biology of Kangaroos, Wallabies, and Rat-Kangaroos (eds Coulson, G. & Eldridge, M.) 137–151 (CSIRO Publishing, 2010).
Google Scholar 
Bishop, N. Functional anatomy of the macropodid pes. Proc. Linn. Soc. New South Wales 117, 17–50 (1997).ADS 

Google Scholar 
Szalay, F. S. Evolutionary History of the Marsupials and an Analysis of Osteological Characters. (Cambridge University Press, 1994).Veth, P. et al. Early human occupation of a maritime desert, Barrow Island, north-west Australia. Quat. Sci. Rev. 168, 19–29. https://doi.org/10.1016/j.quascirev.2017.05.002 (2017).Article 
ADS 

Google Scholar 
Moro, D. & Lagdon, R. History and environment of Barrow Island. Rec. West. Aust. Mus. Suppl. 83, 1–8. https://doi.org/10.18195/issn.0313-122x.83.2013.001-008 (2013).Veth, P., Ditchfield, K. & Hook, F. Maritime deserts of the Australian northwest. Aust. Archaeol. 79, 156–166. https://doi.org/10.1080/03122417.2014.11682032 (2014).Article 

Google Scholar 
Morse, K. Coastwatch: Pleistocene resource use on the Cape Range peninsula. In Australian Coastal Archaeology (eds Hall, J. & McNiven, I. J.) 73–78 (ANH Publications, 1999).
Google Scholar 
Baynes, A. & McDowell, M. C. The original mammal fauna of the Pilbara biogeographic region of north-western Australia. Rec. West. Aust. Mus. Suppl. 78, 285–298. https://doi.org/10.18195/issn.0313-122x.78(1).2010.285-298 (2010).Article 

Google Scholar 
Shortridge, G. C. An account of the geographical distribution of the marsupials and monotremes of south-west Australia, having special reference to the specimens collected during the Balston expedition of 1904–1907. Proc. Zool. Soc. Lond. 74, 803–848. https://doi.org/10.1111/j.1469-7998.1910.tb06974.x (1909).Article 

Google Scholar 
Ballard, C. K. Use of Epiphyseal and Total Fusion Scores as Methods of Age Estimation and Evaluation of Morphological Indices in the Macropodidae. (Northern Illinois University, 2007).Rose, R. W. Age estimation of the Tasmanian bettong (Bettongia gaimardi) (Marsupialia: Potoroidae). Wildl. Res. 16, 251–261. https://doi.org/10.1071/WR9890251 (1989).Article 

Google Scholar 
Johnson, P. M. & Delean, S. Reproduction in the northern bettong, Bettongia tropica Wakefield (Marsupialia: Potoroidae), in captivity, with age estimation and development of pouch young. Wildl. Res. 28, 79–85. https://doi.org/10.1071/WR00007 (2001).Article 

Google Scholar 
Thompson, C. K., Wayne, A. F., Godfrey, S. S. & Andrew Thompson, R. C. Survival, age estimation and sexual maturity of pouch young of the brush-tailed bettong (Bettongia penicillata) in captivity. Aust. Mammal. 37, 29–38. https://doi.org/10.1071/AM14025 (2015).Article 

Google Scholar 
Janis, C. M. Correlation of cranial and dental variables with dietary preferences in mammals: A comparison of macropodoids and ungulates. Mem. – Queensl. Museum 28, 349–366 (1990).
Google Scholar 
Sharman, G. B., Frith, H. J. & Calaby, J. H. Growth of the pouch young, tooth eruption and age determination in the red kangaroo, Megaleia rufa. CSIRO Wildl. Res. 9, 20–49. https://doi.org/10.1071/cwr9640020 (1964).Article 

Google Scholar 
Newsome, A. E., Merchant, J. C., Bolton, B. L. & Dudziński, M. L. Sexual dimorphism in molar progression and eruption in the agile wallaby. Wildl. Res. 4, 1–5. https://doi.org/10.1071/WR9770001 (1977).Article 

Google Scholar 
Poole, W. E., Merchant, J. C., Carpenter, S. M. & Calaby, J. H. Reproduction, growth and age determination in the yellow-footed rock-wallaby Petrogale xanthopus Gray, in captivity. Wildl. Res. 12, 127–136. https://doi.org/10.1071/WR9850127 (1985).Article 

Google Scholar 
Delaney, R. & Marsh, H. Estimating the age of wild rock-wallabies by dental radiography: A basis for quantifying the age structure of a discrete colony of Petrogale assimilis. Wildl. Res. 22, 547–559. https://doi.org/10.1071/WR9950547 (1995).Article 

Google Scholar 
Kido, N., Tanaka, S., Wada, Y., Sato, S. & Omiya, T. Molar eruption and identification of the eastern grey kangaroo (Macropus giganteus) at different ages. J. Vet. Med. Sci. 80, 648–652. https://doi.org/10.1292/jvms.17-0069 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Koo, T. K. & Li, M. Y. A guideline of selecting and reporting intraclass correlation coefficients for reliability research. J. Chiropr. Med. 15, 155–163. https://doi.org/10.1016/j.jcm.2016.02.012 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Claude, J. Log-shape ratios, Procrustes superimposition, elliptic Fourier analysis: Three worked examples in R. Hystrix 24, 94–102. https://doi.org/10.4404/hystrix-24.1-6316 (2013).Article 

Google Scholar 
Mosimann, J. E. Size allometry: Size and shape variables with characterizations of the Lognormal and generalized gamma distributions. J. Am. Stat. Assoc. 65, 930–945. https://doi.org/10.2307/2284599 (1970).Article 
MATH 

Google Scholar 
Kovarovic, K., Aiello, L. C., Cardini, A. & Lockwood, C. A. Discriminant function analyses in archaeology: Are classification rates too good to be true ?. J. Archaeol. Sci. 38, 3006–3018. https://doi.org/10.1016/j.jas.2011.06.028 (2011).Article 

Google Scholar 
Ramayah, T. et al. Discriminant analysis: An illustrated example. Afr. J. Bus. Manag. 4, 1654–1667 (2010).
Google Scholar 
Sanchez, P. M. The unequal group size problem in discriminant analysis. J. Acad. Mark. Sci. 2, 629–633. https://doi.org/10.1007/BF02729456 (1974).Article 

Google Scholar 
Hair, J. F., Black, W. C., Babin, B. J. & Anderson, R. E. Multivariate Data Analysis. (Cengage, 2018).Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S. (Springer, 2002).Fox, J. & Weisberg, S. An R Companion to Applied Regression. (Sage 2019).Harrell, F. E. & Dupont, C. Harrell Miscellaneous. (2021).Oksanen, J. et al. Community Ecology Package (2020).Kassambara, A. Pipe-Friendly Framework for Basic Statistical Tests. (2021).Korkmaz, S., Goksuluk, D. & Zararsiz, G. MVN: An R Package for Assessing Multivariate Normality. R J. 6, 151–162 (2014).Article 

Google Scholar 
Revelle, W. Procedures for Psychological, Psychometric, and Personality Research. (2022).Weisbecker, V. et al. Individual variation of the masticatory system dominates 3D skull shape in the herbivory-adapted marsupial wombats. Front. Zool. 16, 41. https://doi.org/10.1186/s12983-019-0338-5 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Richards, J. D. et al. The biology of banded (Lagostrophus fasciatus) and rufous (Lagorchestes hirsutus) hare-wallabies (Diprotodontia: Macropodidae) on Dorre and Bernier Islands, Western Australia. Wildl. Res. 28, 311–322. https://doi.org/10.1071/WR99109 (2001).Article 
ADS 

Google Scholar 
Ingleby, S. & Westoby, M. Habitat requirements of the spectacled hare-wallaby (Lagorchestes conspicillatus) in the Northern Territory and Western Australia. Wildl. Res. 19, 721–741. https://doi.org/10.1071/WR9920721 (1992).Article 

Google Scholar 
Helgen, K. M. & Flannery, T. F. Taxonomy and historical distribution of the wallaby genus Lagostrophus. Aust. J. Zool. 51, 199–212. https://doi.org/10.1071/ZO02078 (2003).Article 

Google Scholar 
McDowell, M. C. et al. Morphological and molecular evidence supports specific recognition of the recently extinct Bettongia anhydra (Marsupialia: Macropodidae). J. Mammal. 96, 287–296. https://doi.org/10.1093/jmammal/gyv006 (2015).Article 

Google Scholar 
Ingleby, S. Distribution and status of the northern nailtail wallaby, Onychogalea unguífera (Gould, 1841). Wildl. Res. 18, 655–676. https://doi.org/10.1071/WR9910655 (1991).Article 

Google Scholar 
Peters, C. et al. Species identification of Australian marsupials using collagen fingerprinting. R. Soc. Open Sci. 8, 211229. https://doi.org/10.1098/rsos.211229 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Prince, R. I. T. Banded hare-wallaby. In Mammals of Australia (eds Strahan, R. & van Dyck, S.) 406–408 (Reed New Holland, 2008).
Google Scholar 
De Deckker, P., Barrows, T. T. & Rogers, J. Land-sea correlations in the Australian region: Post-glacial onset of the monsoon in northwestern Western Australia. Quat. Sci. Rev. 105, 181–194. https://doi.org/10.1016/j.quascirev.2014.09.030 (2014).Article 
ADS 

Google Scholar 
Ward, I. et al. 50,000 years of archaeological site stratigraphy and micromorphology in Boodie Cave, Barrow Island, Western Australia. J. Archaeol. Sci. Rep. 15, 344–369. https://doi.org/10.1016/j.jasrep.2017.08.012 (2017).Article 

Google Scholar 
Skippington, J., Manne, T. & Veth, P. Isotopic indications of late Pleistocene and Holocene paleoenvironmental changes at Boodie Cave archaeological site, Barrow Island, Western Australia. Molecules 26, 2585. https://doi.org/10.3390/molecules26092582 (2021).Article 
CAS 

Google Scholar 
Baynes, A. & Jones, B. The mammals of Cape Range Peninsula, north-western Australia. Rec. West. Aust. Mus. Suppl. 45, 207–255 (1993).
Google Scholar 
Piper, C. & Veth, P. Palaeoecology and sea level changes: Decline of mammal species richness during late Quaternary island formation in the Montebello Islands, north-western Australia. Palaeontol. Electron. 24, a20. https://doi.org/10.26879/1050 (2021).Article 

Google Scholar 
Lyman, R. L. The history of ‘laundry lists’ in North American zooarchaeology. J. Anthropol. Archaeol. 39, 42–50. https://doi.org/10.1016/j.jaa.2015.02.003 (2015).Article 

Google Scholar 
Guillaud, E., Cornette, R. & Béarez, P. Is vertebral form a valid species-specific indicator for salmonids? The discrimination rate of trout and Atlantic salmon from archaeological to modern times. J. Archaeol. Sci. 65, 84–92. https://doi.org/10.1016/j.jas.2015.11.010 (2016).Article 

Google Scholar 
Monchot, H. & Gendron, D. Disentangling long bones of foxes (Vulpes vulpes and Alopex lagopus) from artic archaeological sites. J. Archaeol. Sci. 37, 799–806. https://doi.org/10.1016/j.jas.2009.11.009 (2010).Article 

Google Scholar  More

Estes, J. et al. Megafaunal impacts on structure and function of ocean ecosystems. Annu. Rev. Environ. Resour. 41, 83–116 (2016).
Google Scholar 
Ferretti, F., Worm, B., Britten, G., Heithaus, M. & Lotze, H. Patterns and ecosystem consequences of shark declines in the ocean. Ecol. Lett. 13, 1055–1071 (2010).PubMed 

Google Scholar 
Heithaus, M. R. et al. Seagrasses in the age of sea turtle conservation and shark overfishing. Front. Mar. Sci. 1, 1–6 (2014).
Google Scholar 
McCauley, D. et al. Marine defaunation: Animal loss in the global ocean. Science 347, 1255641 (2015).PubMed 

Google Scholar 
Pereira, H. et al. Scenarios for global biodiversity in the 21st century. Science 330, 1496–1501 (2010).ADS 
CAS 
PubMed 

Google Scholar 
Oliver, S., Braccini, M., Newman, S. & Harvey, E. S. Global patterns in the bycatch of sharks and rays. Mar. Policy 54, 86–97 (2015).
Google Scholar 
Pacoureau, N. et al. Half a century of global decline in oceanic sharks and rays. Nature 589, 567–574 (2021).ADS 
CAS 
PubMed 

Google Scholar 
Gilman, E. et al. Shark interactions in pelagic longline fisheries. Mar. Policy 32, 1–18 (2008).
Google Scholar 
Worm, B. et al. Global catches, exploitation rates, and rebuilding options for sharks. Mar. Policy 40, 194–204 (2013).
Google Scholar 
Bowlby, H. & Gibson, A. Implications of life history uncertainty when evaluating status in the Northwest Atlantic population of white shark (Carcharodon carcharias). Ecol. Evol. 10, 4990–5000 (2020).PubMed 
PubMed Central 

Google Scholar 
Dulvy, N. et al. Overfishing drives over one third of all sharks and rays toward a global extinction crisis. Curr. Biol. 31, 4773-4787.e8 (2021).CAS 
PubMed 

Google Scholar 
Heino, M., Pauli, B. & Dieckmann, U. Fisheries-induced evolution. Annu. Rev. Ecol. Evol. Syst. 46, 461–480 (2015).
Google Scholar 
Mitchell, J., McLean, D., Collins, S. & Langlois, T. Shark depredation in commercial and recreational fisheries. Rev. Fish Biol. Fish 28, 715–748 (2018).
Google Scholar 
Jaiteh, V. F., Loneragan, N. & Warren, C. The end of shark finning? Impacts of declining catches and fin demand on coastal community livelihoods. Mar. Policy 82, 224–233 (2017).
Google Scholar 
Seidu, I. et al. Fishing for survival: Importance of shark fisheries for the livelihoods of coastal communities in Western Ghana. Fish. Res. 246, 106157 (2022).
Google Scholar 
Gilman, E., Weijerman, M. & Suuronen, P. Ecological data from observer programs underpin ecosystem-based fisheries management. ICES J. Mar. Sci. 74, 1481–1495 (2017).
Google Scholar 
Melnychuk, M. et al. Identifying management actions that promote sustainable fisheries. Nat. Sustain. https://doi.org/10.1038/s41893-020-00668-1 (2021).Article 

Google Scholar 
Musyl, M. & Gilman, E. Meta-analysis of post-release fishing mortality in apex predatory pelagic sharks and white marlin. Fish Fish. 20, 466–500 (2019).
Google Scholar 
Clarke, S. A status snapshot of key shark species in the western and central pacific and potential management options. in WCPFC-SC7-2011/EB-WP-04. Western and Central Pacific Fisheries Commission, Kolonia, Federated States of Micronesia (2011).Dapp, D., Walker, T., Huveneers, C. & Reina, R. Respiratory mode and gear type are important determinants of elasmobranch immediate and post-release mortality. Fish Fish. 17, 507–524 (2016).
Google Scholar 
ICES. Report of the working group on elasmobranch fishes. in ICES CM 2018/ACOM:16. International Council for the Exploration of the Sea, Copenhagen (2018).Dicks, L. et al. A transparent process for “evidence-informed” policy making. Conserv. Lett. 7, 119–125 (2014).
Google Scholar 
Nichols, J., Kendall, W. & Boomer, G. Accumulating evidence in ecology: Once is not enough. Ecol. Evol. 9, 13991–14004 (2019).PubMed 
PubMed Central 

Google Scholar 
Nakagawa, S. et al. Meta-analysis of variation: Ecological and evolutionary applications and beyond. Methods Ecol. Evol. 6, 143–152 (2015).
Google Scholar 
Pfaller, J., Chaloupka, M., Bolten, A. & Bjorndal, K. Phylogeny, biogeography and methodology: A meta-analytic perspective on heterogeneity in adult marine turtle survival rates. Sci. Rep. 8, 5852. https://doi.org/10.1038/s41598-018-24262-w (2018).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Godin, A., Carlson, J. & Burgener, V. The effect of circle hooks on shark catchability and at-vessel mortality rates in longlines fisheries. Bull. Mar. Sci. 88, 469–483 (2012).
Google Scholar 
Reinhardt, J. et al. Catch rate and at-vessel mortality of circle hooks versus J-hooks in pelagic longline fisheries: A global meta-analysis. Fish Fish. 19, 413–430 (2018).
Google Scholar 
Rosa, D., Santos, C. & Coelho, R. Assessing the effects of hook, bait and leader type as potential mitigation measures to reduce bycatch and mortality rates of shortfin mako: A meta-analysis with comparisons for target, bycatch and vulnerable fauna interactions. in ICCAT Collective Volume of Scientifics Papers 76, 247–278 (2020).Santos, C., Rosa, D. & Coelho, R. Hook, bait and leader type effects on surface pelagic longline retention and mortality rates: A meta-analysis with comparisons for target, bycatch and vulnerable fauna interactions. in IOTC-2019-WPEB15-39. Indian Ocean Tuna Commission, Mahe, Seychelles (2019).Santos, C., Rosa, D. & Coelho, R. Progress on a meta-analysis for comparing hook, bait and leader effects on target, bycatch and vulnerable fauna interactions. in Collective Volume of Scientifics Papers ICCAT 77, 182–217 (2020).Condamine, F., Romieu, J. & Guinot, G. Climate cooling and clade competition likely drove the decline of lamniform sharks. PNAS 116, 20584–20590 (2019).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Vehtari, A., Gelman, A., Simpson, D., Carpenter, B. & Bürkner, P. Rank-normalization, folding, and localization: An improved Rhat for assessing convergence of MCMC (with Discussion). Bayesian Anal. 16, 667–718 (2021).MathSciNet 

Google Scholar 
Cinar, O., Nakagawa, S. & Viechtbauer, W. Phylogenetic multilevel meta-analysis: A simulation study on the importance of modelling the phylogeny. Methods Ecol. Evol. 13, 383–395 (2022).
Google Scholar 
Lajeunesse, M. Meta-analysis and the comparative phylogenetic method. Am. Nat. 174, 369–381 (2009).PubMed 

Google Scholar 
Chamberlain, S. et al. Does phylogeny matter? Assessing the impact of phylogenetic information in ecological meta-analysis. Ecol. Lett. 15, 627–636 (2012).PubMed 

Google Scholar 
Burns, J. & Strauss, S. More closely related species are more ecologically similar in an experimental test. Proc. Natl. Acad. Sci. USA 108, 5302–5307 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Cachera, M. & Le Loc’h, F. Assessing the relationships between phylogenetic and functional singularities in sharks (Chondrichthyes). Ecol. Evol. 7, 6292–6303 (2017).PubMed 
PubMed Central 

Google Scholar 
Münkemüller, T., Boucher, F. C., Thuiller, W. & Lavergne, S. Phylogenetic niche conservatism—Common pitfalls and ways forward. Funct. Ecol. 29, 627–639 (2015).PubMed 
PubMed Central 

Google Scholar 
Bazzi, M., Campione, N., Kear, B., Pimiento, C. & Ahlberg, P. Feeding ecology has shaped the evolution of modern sharks. Curr. Biol. 31, 5138–5148 (2021).CAS 
PubMed 

Google Scholar 
Sepulveda, C., Wegner, N., Bernal, D. & Graham, J. The red muscle morphology of the thresher sharks (family Alopiidae). J. Exp. Biol. 208, 4255–4261 (2005).CAS 
PubMed 

Google Scholar 
Wosnick, N. et al. Multispecies thermal dynamics of air-exposed ectothermic sharks and its implications for fisheries conservation. J. Exp. Mar. Biol. Ecol. 513, 1–9 (2019).
Google Scholar 
French, R. et al. High survivorship after catch-and-release fishing suggests physiological resilience in the endothermic shortfin mako shark (Isurus oxyrinchus). Conserv. Physiol. https://doi.org/10.1093/conphys/cov044 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Davis, M. Key principles for understanding fish bycatch discard mortality. Can. J. Fish. Aquat. Sci. 59, 1834–1843 (2002).
Google Scholar 
Massey, Y., Sabarros, P., Rabearisoa, N. & Bach, P. Drivers of at-haulback mortality of sharks caught during pelagic longline fishing experiments. in IOTC-2019-WPEB15-14_Rev1. Indian Ocean Tuna Commission, Mahe, Seychelles (2019).Musyl, M., Moyes, C., Brill, R. & Fragoso, N. Factors influencing mortality estimates in post-release survival studies: Comment on Campana et al. (2009). Mar. Ecol. Prog. Ser. 396, 157–159 (2009).Pimiento, C., Cantalapiedra, J., Shimada, K., Field, D. & Smaers, J. Evolutionary pathways towards gigantism in sharks and rays. Evolution 73, 588–599 (2019).PubMed 

Google Scholar 
Musyl, M. & Gilman, E. Post-release fishing mortality of blue (Prionace glauca) and silky shark (Carcharhinus falciformes) from a Palauan-based commercial longline fishery. Rev. Fish Biol. Fish. 28, 567–658 (2018).
Google Scholar 
Childs, D., Sheldon, B. & Rees, M. The evolution of labile traits in sex- and age-structured populations. J. Anim. Ecol. 85, 329–342 (2016).PubMed 
PubMed Central 

Google Scholar 
Comte, L., Murienne, J. & Grenouillet, G. Species traits and phylogenetic conservatism of climate-induced range shifts in stream fishes. Nat. Commun. 5, 5053 (2014).ADS 
PubMed Central 

Google Scholar 
IUCN. The IUCN Red List of Threatened Species. Version 2021-3. www.iucnredlist.org. ISSN 2307-8235 (International Union for the Conservation of Nature, Gland, Switzerland, 2022).García, V., Lucifora, L. & Ransom, M. The importance of habitat and life history to extinction risk in sharks, skates, rays and chimaeras. Proc. R. Soc. B 275, 83–89 (2008).PubMed 

Google Scholar 
Cortes, E. Perspectives on the intrinsic rate of population growth. Methods Ecol. Evol. 7, 1136–1145 (2016).
Google Scholar 
Ellis, J. et al. A review of capture and post-release mortality of elasmobranchs. J. Fish Biol. 90, 653–722 (2017).CAS 
PubMed 

Google Scholar 
Gallagher, A., Orbesen, E., Hammerschlag, N. & Serafy, J. Vulnerability of oceanic sharks as pelagic longline bycatch. Glob. Ecol. Conserv. 1, 50–59 (2014).
Google Scholar 
Afonso, A., Santiago, R., Hazin, H. & Hazin, F. Shark bycatch and mortality and hook bite-offs in pelagic longlines: Interactions between hook types and leader materials. Fish. Res. 131–133, 9–14 (2012).
Google Scholar 
Gilman, E., Chaloupka, M. & Musyl, M. Effects of pelagic longline hook size on species- and size-selectivity and survival. Rev. Fish Biol. Fish. 28, 417–433 (2018).
Google Scholar 
Epperly, S., Watson, J., Foster, D. & Shah, A. Anatomical hooking location and condition of animals captured with pelagic longlines: The grand banks experiments 2002–2003. Bull. Mar. Sci. 88, 513–527 (2012).
Google Scholar 
Amorim, S., Santos, M., Coelho, R. & Fernandez-Carvalho, J. Effects of 17/0 circle hooks and bait on fish catches in a southern Atlantic swordfish longline fishery. Aquat. Conserv. 25, 518–533 (2014).
Google Scholar 
Coelho, R., Fernandez-Carvalho, J., Lino, P. & Santos, M. An overview of the hooking mortality of elasmobranchs caught in a swordfish pelagic longline fishery in the Atlantic Ocean. Aquat. Living Resour. 25, 311–319 (2012).
Google Scholar 
Gilman, E. et al. A decision support tool for integrated fisheries bycatch management. Rev. Fish Biol. Fish. https://doi.org/10.1007/s11160-021-09693-5 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
Pascoe, S. et al. Use of incentive-based management systems to limit bycatch and discarding. Int. Rev. Environ. Resour. Econ. 4, 123–161 (2010).
Google Scholar 
Somers, K., Pfeiffer, L., Miller, S. & Morrison, W. Using incentives to reduce bycatch and discarding: Results under the west coast catch share program. Coast. Manag. 46, 1–17 (2019).
Google Scholar 
Abbott, J. & Wilen, J. Regulation of fisheries bycatch with common-pool output quotas. J. Environ. Econ. Manag. 57, 195–204 (2009).MATH 

Google Scholar 
Gilman, E. et al. Increasing the functionalities and accuracy of fisheries electronic monitoring systems. Aquat. Conserv. 29, 901–926 (2019).
Google Scholar 
Watling, J. Fishing observers ‘intimidated and bribed by EU crews’. Quota checks allegedly being compromised aboard Northwest Atlantic fishery boats, as observers report surveillance and theft. The Guardian (2012, accessed 21 July 2022). https://www.theguardian.com/environment/2012/may/18/fishing-inspectors-intimidated-bribed-crews.Clarke, S., Harley, S., Hoyle, S. & Rice, J. Population trends in Pacific oceanic sharks and the utility of regulations on shark finning. Conserv. Biol. 27, 197–209 (2013).PubMed 

Google Scholar 
Tolotti, M. T. et al. Banning is not enough: The complexities of oceanic shark management by tuna regional fisheries management organizations. Glob. Ecol. Conserv. 4, 1–7 (2015).
Google Scholar 
Gilman, E., Chaloupka, M., Merrifield, M., Malsol, N. & Cook, C. Standardized catch and survival rates, and effect of a ban on shark retention, Palau pelagic longline fishery. Aquat. Conserv. 26, 1031–1062 (2016).
Google Scholar 
CITES. Appendices I, II and III. Valid from 22 June 2021. Convention on International Trade in Endangered Species of Wild Fauna and Flora, United Nations Environment Program, Geneva (2021).Ward-Paige, C. A global overview of shark sanctuary regulations and their impact on shark fisheries. Mar. Policy 82, 87–97 (2017).
Google Scholar 
E.U. Regulation (E.U.) No 1380/2013 of the European Parliament and of the Council of 11 December 2013 on the Common Fisheries Policy, amending Council Regulations (EC) No 1954/2003 and (EC) No 1224/2009 and repealing Council Regulations (EC) No 2371/2002 and (EC) No 639/2004 and Council Decision 2004/585/EC. Official Journal of the European Union L354, 22–61 (2013).FAO. International Guidelines on Bycatch Management and Reduction of Discards (Food and Agriculture Organization of the United Nations, Rome, 2011).CCSBT. Resolution to Align CCSBT’s Ecologically Related Species Measures with those of other Tuna RFMOs (Commission for the Conservation of Southern Bluefin Tuna, Deakin West, Australia, 2021).IATTC. Active Resolutions and Recommendations (Inter-American Tropical Tuna Commission, La Jolla, 2022).ICCAT. Compendium. Management Recommendations and Resolutions Adopted by ICCAT for the Conservation of Atlantic Tunas and Tuna-like Species (International Commission for the Conservation of Atlantic Tunas, Madrid, 2021).IOTC. Compendium of Active Conservation and Management Measures for the Indian Ocean Tuna Commission (Indian Ocean Tuna Commission, Mahe, 2021).WCPFC. Conservation and Management Measures and Resolutions of the Western and Central Pacific Fisheries Commission. Compiled 31 August 2021 (Western and Central Pacific Fisheries Commission, Kolonia, Federated States of Micronesia, 2021).Faith, D. Threatened species and the potential loss of phylogenetic diversity: Conservation scenarios based on estimated extinction probabilities and phylogenetic risk analysis. Conserv. Biol. 22, 1461–1470 (2008).PubMed 

Google Scholar 
Dolce, J. & Wilga, C. Evolutionary and ecological relationships of gill slit morphology in extant sharks. Bull. Mus. Comp. 161, 79–109 (2013).
Google Scholar 
MacLeod, N. & Forey, P. Morphology, Shape and Phylogeny (CRC Press, 2002).
Google Scholar 
Haddaway, N., Macura, B., Whaley, P. & Pullin, A. ROSES RepOrting standards for Systematic Evidence Syntheses: pro forma, flow-diagram and descriptive summary of the plan and conduct of environmental systematic reviews and systematic maps. Environ. Evid. https://doi.org/10.1186/s13750-018-0121-7 (2018).Article 

Google Scholar 
Pullin, A., Frampton, G., Livoreil, B. & Petrokofsky, G., Eds. Section 5. Conducting a Search. Key CEE Standards for Conduct and Reporting. In Pullin, A., Frampton, G., Livoreil, B., Petrokofsky, G., Eds. Guidelines and Standards for Evidence Synthesis in Environmental Management. Version 5.0. Collaboration for Environmental Evidence (2020).Pullin, A., Frampton, G., Livoreil, B. & Petrokofsky, G. (eds) Section 3. Planning a CEE Evidence Synthesis. In Pullin, A., Frampton, G., Livoreil, B., Petrokofsky, G. (eds) Guidelines and Standards for Evidence Synthesis in Environmental Management. Version 5.0. Collaboration for Environmental Evidence (2021).Page, M. et al. The PRISMA statement: An updated guideline for reporting systematic reviews. BMJ https://doi.org/10.1136/bmj.n.71 (2020).Article 
PubMed 

Google Scholar 
Tuyl, F., Gerlach, R. & Mengersen, K. Comparison of Bayes-Laplace, Jeffreys, and other priors: The case of zero events. Am. Stat. 62, 40–44 (2008).MathSciNet 

Google Scholar 
Dorai-Raj, S. binom: Binomial confidence intervals for several parameterizations. R package version 1.1-1. https://CRAN.R-project.org/package=binom (2014).van Lissa, C. Small sample meta-analyses: Exploring heterogeneity using MetaForest. Chapter 13. In Small Sample Size Solutions: A Guide for Applied Researchers and Practitioners (eds Van De Schoot, R. & Miočević, M.) 186–202 (Routledge, Oxford, 2020).
Google Scholar 
Wright, M. & Ziegler, A. ranger: A fast implementation of random forests for high dimensional data in C++ and R. J. Stat. Softw. 77, 1–17 (2017).
Google Scholar 
Janitza, S., Celik, E. & Boulesteix, A. A computationally fast variable importance test for random forests for high-dimensional data. Adv. Data Anal. Classif. 12, 885–915 (2018).MathSciNet 
MATH 

Google Scholar 
Mayer, M. missRanger: Fast imputation of missing values. R package version 2.1.3. https://CRAN.R-project.org/package=missRanger (2021).Konstantopoulos, S. Fixed effects and variance components estimation in three-level meta-analysis. Research Synthesis. Methods 2, 61–76 (2011).
Google Scholar 
Amaral, C., Pereira, F., Silva, D., Amorim, A. & de Carvalho, E. The mitogenomic phylogeny of the Elasmobranchii (Chondrichthyes). Mitochondrial DNA A 29, 867–878 (2017).
Google Scholar 
Hara, Y. et al. Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates. Nat. Ecol. Evol. 2, 1761–1771 (2018).PubMed 

Google Scholar 
Naylor, G. et al. A DNA sequence-based approach to the identification of shark and ray species and its implications for global elasmobranch diversity and parasitology. Bull. Am. Mus. Nat. 367, 1–262 (2012).
Google Scholar 
Stein, R. et al. Global priorities for conserving the evolutionary history of sharks, rays and chimaeras. Nat. Ecol. Evol. 2, 288–298 (2018).PubMed 

Google Scholar 
Maddison, D., Swofford, D., Maddison, W. & Cannatella, D. Nexus: An extensible file format for systematic information. Syst. Biol. 46, 590–621 (1997).CAS 
PubMed 

Google Scholar 
Upham, N., Esselstyn, J. & Jetz, W. Inferring the mammal tree: Species-level sets of phylogenies for questions in ecology, evolution, and conservation. PLoS Biol. 17, e3000494 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Paradis, E. & Schliep, K. ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528 (2019).CAS 
PubMed 

Google Scholar 
Yu, G. Using ggtree to visualize data on tree-like structures. Curr. Protoc. Bioinform. 69, e96. https://doi.org/10.1002/cpbi.96 (2020).Article 

Google Scholar 
Hadfield, J. & Nakagawa, S. General quantitative genetic methods for comparative biology: Phylogenies, taxonomies and multi-trait models for continuous and categorical characters. J. Evol. Biol 23, 494–508 (2010).CAS 
PubMed 

Google Scholar 
Lin, L. & Chu, H. Meta-analysis of proportions using generalized linear mixed models. Epidemiology 31, 713–717 (2020).PubMed 
PubMed Central 

Google Scholar 
Carpenter, B. et al. Stan: A probabilistic programming language. J. Stat. Softw. 76, 1–32 (2017).
Google Scholar 
Bürkner, P. brms: An R Package for Bayesian multilevel models using Stan. J. Stat. Softw. 81, 1–28 (2017).
Google Scholar 
Günhan, B., Röver, C. & Friede, T. Random-effects meta-analysis of few studies involving rare events. Res. Synth. Methods 11, 74–90 (2020).PubMed 

Google Scholar 
Pappalardo, P. et al. Comparing traditional and Bayesian approaches to ecological meta-analysis. Methods Ecol. Evol. 11, 1286–1295 (2020).
Google Scholar 
Ott, M., Plummer, M. & Roos, M. How vague is vague? How informative is informative? Reference analysis for Bayesian meta-analysis. Stat. Med. 40, 4505–4521 (2021).MathSciNet 
PubMed 
PubMed Central 

Google Scholar 
Wood, S. Generalized Additive Models: An Introduction with R 2nd edn. (Chapman and Hall, 2017).MATH 

Google Scholar 
Kruschke, J. & Liddell, T. The Bayesian new statistics: Hypothesis testing, estimation, meta-analysis, and power analysis from a Bayesian perspective. Psychon. Bull. Rev. 25, 178–206 (2018).PubMed 

Google Scholar 
Kay, M. tidybayes: Tidy data and geoms for Bayesian models. R package version 2.1.1. https://doi.org/10.5281/zenodo.1308151 (2020).Makowski, D., Ben-Shachar, M. & Lüdecke, D. bayestestR: Describing effects and their uncertainty, existence and significance within the Bayesian framework. J. Open Source Softw. 4, 1541 (2019).ADS 

Google Scholar 
Searle, S., Speed, F. & Milliken, G. Population marginal means in the linear model: An alternative to least squares means. Am. Stat. 34, 216–221 (1980).MathSciNet 
MATH 

Google Scholar 
Lenth, R. emmeans: Estimated marginal means, aka least-squares means. R package version 1.5.2-1. https://CRAN.R-project.org/package=emmeans (2020).Pagel, M. Inferring the historical patterns of biological evolution. Nature 401, 877–884 (1999).ADS 
CAS 
PubMed 

Google Scholar 
Münkemüller, T. et al. How to measure and test phylogenetic signal. Methods Ecol. Evol. 3, 743–756 (2012).
Google Scholar 
Vehtari, A., Gelman, A. & Gabry, J. Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Stat. Comput. 27, 1413–1432 (2017).MathSciNet 
MATH 

Google Scholar 
Yao, Y., Vehtari, A., Simpson, D. & Gelman, A. Using stacking to average Bayesian predictive distributions (with Discussion). Bayesian Anal. 13, 917–1003 (2018).MathSciNet 
MATH 

Google Scholar 
Gabry, J., Simpson, D., Vehtari, A., Betancourt, M. & Gelman, A. Visualization in Bayesian workflow. J. R. Soc. Ser. A 182, 1–14 (2019).MathSciNet 

Google Scholar 
Lazic, S., Mellor, J., Ashby, M. & Munafo, M. A Bayesian predictive approach for dealing with pseudoreplication. Sci. Rep. 10, 2020. https://doi.org/10.1038/s41598-020-59384-7 (2020).Article 
CAS 

Google Scholar 
Page, M., Sterne, J., Higgins, J. & Egger, M. Investigating and dealing with publication bias and other reporting biases in meta-analyses of health research: A review. Res. Synth. Methods 12, 248–259 (2021).PubMed 

Google Scholar 
Peters, J., Sutton, A., Jones, D., Abrams, K. & Rushton, L. Contour-enhanced meta-analysis funnel plots help distinguish publication bias from other causes of asymmetry. J. Clin. Epidemiol. 61, 991–996 (2008).PubMed 

Google Scholar 
Viechtbauer, W. Conducting meta-analyses in R with the metafor package. J. Stat. Softw. 36, 1–48 (2010).
Google Scholar 
Gasparrini, A., Armstrong, B. & Kenward, M. Multivariate meta-analysis for non-linear and other multi-parameter associations. Stat. Med. 31, 3821–3839 (2012).MathSciNet 
CAS 
PubMed 
PubMed Central 

Google Scholar  More

FAO. The State of the World’s Forests 2018—Forest Pathways to Sustainable Development (FAO, 2018).
Google Scholar 
Pan, Y. et al. A large and persistent carbon sink in the world’s forests. Science 333(6045), 988–993. https://doi.org/10.1126/science.1201609 (2011).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Luyssaert, S. et al. Tradeoffs in using European forests to meet climate objectives. Nature 562(7726), 259–262. https://doi.org/10.1038/s41586-018-0577-1 (2018).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Büntgen, U. et al. Limited capacity of tree growth to mitigate the global greenhouse effect under predicted warming. Nat. Commun. https://doi.org/10.1038/s41467-019-10174-4 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Anderegg, W. R. L. et al. Climate-driven risks to the climate mitigation potential of forests. Science 368, 1327. https://doi.org/10.1126/science.aaz7005 (2020).Article 
CAS 

Google Scholar 
Jiang, M. et al. The fate of carbon in a mature forest under carbon dioxide enrichment. Nature 580(7802), 227–231. https://doi.org/10.1038/s41586-020-2128-9 (2020).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Naudts, K. et al. Europe’s forest management did not mitigate climate warming. Science 351(6273), 597–599. https://doi.org/10.1126/science.aad7270 (2016).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Tong, X. et al. Forest management in southern China generates short term extensive carbon sequestration. Nat. Commun. https://doi.org/10.1038/s41467-019-13798-8 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Yu, K. et al. Effects of stand age on soil respiration in Pinus massoniana plantations in the hilly red soil region of Southern China. CATENA 178, 313–321. https://doi.org/10.1016/j.catena.2019.03.038 (2019).Article 
CAS 

Google Scholar 
Mei, G., Sun, Y. & Sajjad, S. Models for predicting the biomass of Cunninghamia lanceolata trees and stands in southeastern China. PLoS ONE 12, e0169747. https://doi.org/10.1371/journal.pone.0169747 (2017).Article 
CAS 

Google Scholar 
Wu, H. et al. Soil phosphorus bioavailability and recycling increasedwith stand age in Chinese fir plantations. Ecosystems 23, 973–988. https://doi.org/10.1007/s10021-019-00450-1 (2019).Article 

Google Scholar 
State Forestry Administration. General situation of forest resources in China. The 8th National Forest Inventory (State Forestry Administration, 2014).Wang, X. et al. Vegetation carbon storage and density of forest ecosystems in China. Chin. J. Appl. Ecol. 12(1), 13–16 (2001) (in Chinese with English Abstract).ADS 
CAS 

Google Scholar 
Kang, H. et al. Simulating the impact of climate change on the growth of Chinese fir plantations in Fujian province, China. NZ J. For. Sci. 47(1), 20. https://doi.org/10.1186/s40490-017-0102-6 (2017).Article 

Google Scholar 
Lu, Y. et al. A process-based approach to estimate Chinese fir (Cunninghamia lanceolata) distribution and productivity in southern China under climate change. Forests 6, 360–379. https://doi.org/10.3390/f6020360 (2015).Article 

Google Scholar 
Zhang, X. et al. Relative contributions of competition, stand structure, age, and climate factors to tree mortality of Chinese fir plantations: Long-term spacing trials in southern China. For. Ecol. Manag. 465, 118103. https://doi.org/10.1016/j.foreco.2020.118103 (2020).Article 

Google Scholar 
You, R. et al. Variation in wood physical properties and effects of climate for different geographic sources of Chinese fir in subtropical area of China. Sci. Rep. https://doi.org/10.1038/s41598-021-83500-w (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Djomo, A. N., Ibrahima, A., Saborowski, J. & Gravenhorst, G. Allometric equations for biomass estimations in Cameroon and pan moist tropical equations including biomass data from Africa. For. Ecol. Manag. 260(10), 1873–1885. https://doi.org/10.1016/j.foreco.2010.08.034 (2010).Article 

Google Scholar 
Peng, D. et al. Estimating the aboveground biomass for planted forests based on stand age and environmental variables. Remote Sens. 11(19), 2270. https://doi.org/10.3390/rs11192270 (2019).Article 
ADS 

Google Scholar 
Zhou, X. et al. Dynamic allometric scaling of tree biomass and size. Nat. Plants. 7(1), 42–49. https://doi.org/10.1038/s41477-020-00815-8 (2021).Article 
PubMed 

Google Scholar 
Li, L. Study on the tree volume table compilation of Chinese fir in Kaihua Forest Farm (Beijing Forestry University, 2011) http://cdmd.cnki.com.cn/Article/CDMD-10022-1011134655.htm (in Chinese).Wang, J. P. et al. Study on the effect of Chinese fir volume formula on estimating the volume of fir standing trees in different sites. Guizhou For. Technol. 19(1), 26–29 (1991) (in Chinese).
Google Scholar 
Zeng, W. S. et al. Compatible tree volume and aboveground biomass equations for Chinese fir plantation in Guizhou. J. Beijing For. Univ. 33(4), 1–6 (2011) (in Chinese).
Google Scholar 
Xia, Z. S. et al. Construction of tree volume equations for Chinese fir plantation in Guizhou Province, southwestern China. J. Beijing For. Univ. 34(1), 1–5 (2012) (in Chinese).
Google Scholar 
Lin, H. Study on biomass and carbon storage of main coniferous forest in Jiangle state-owned forestry farm. J. Fujian For. Sci. Technol. 45(1), 30–34. https://doi.org/10.13428/j.cnki.fjlk.2018.01.007 (2018) (in Chinese with English Abstract).Article 
ADS 

Google Scholar 
Cai, Z. A study on biomass models of Cunninghamia lanceolata plantation in Fujian. (Beijing Forestry University, 2014), http://cdmd.cnki.com.cn/Article/CDMD-10022-1014327550.htm (in Chinese).Chen, G. et al. Carbon storage in a chronosequence of Chinese fir plantations in southern China. For. Ecol. Manag. 300, 68–76. https://doi.org/10.1016/j.foreco.2012.07.046 (2013).Article 

Google Scholar 
Zhang, G. et al. Biomass Characteristics of dominant tree species (group) at Lingnan forest farm in Anhui province. Scientia Silvae Sinicae. 48(5), 136–140. https://doi.org/10.1007/s11783-011-0280-z (2012) (in Chinese with English abstract).Article 
ADS 
CAS 

Google Scholar 
Shi, W. et al. Biomass model and carbon storage of Chinese fir plantation in Dabieshan Mountains in Anhui. Resour. Environ. Yangtze Basin. 24(5), 758–764. https://doi.org/10.11870/cjlyzyyhj201505007 (2015) (in Chinese with English abstract).Article 

Google Scholar 
Li, H. & Zhao, P. Improving the accuracy of tree-level aboveground biomass equations with height classification at a large regional scale. For. Ecol. Manag. 289, 153–163. https://doi.org/10.1016/j.foreco.2012.10.002 (2013).Article 

Google Scholar 
Zeng, W. & Tang, S. A new general allometric biomass model. Nat. Precedings. https://doi.org/10.1038/npre.2011.6704.1 (2011).Article 

Google Scholar 
Schumacher, F. X. & Hall, F. D. S. Logarithmic expression of timber-tree volume. J. Agric. Res. 47(9), 719–734 (1933).
Google Scholar 
Honer, T. G. A new total cubic foot volume function. For. Chron. 41(4), 476–493. https://doi.org/10.5558/tfc41476-4 (1965).Article 

Google Scholar 
Burkhart, H. E. Cubic-foot volume of loblolly pine to any merchantable top limit. South. J. Appl. For. 2, 7–9. https://doi.org/10.1093/sjaf/1.2.7 (1977).Article 

Google Scholar 
Lee, D., Seo, Y. & Choi, J. Estimation and validation of stem volume equations for Pinus densiflora, Pinus koraiensis, and Larix kaempferi in South Korea. For. Sci. Technol. 13(2), 77–82. https://doi.org/10.1080/21580103.2017.1315963 (2017).Article 

Google Scholar 
Chen, B. H. & Chen, C. Y. A preliminary study on the biomass and productivity of Picea koraiensis forests in the dunes. Scientia Silvae Sinicae 4, 269–278 (1980) (in Chinese).
Google Scholar 
Niklas, K. J. Plant Allometry: The Scaling of Form and Process (University of Chicago Press, 1994).
Google Scholar 
Ketterings, Q. M. et al. Reducing uncertainty in the use of allometric biomass equations for predicting above-ground tree biomass in mixed secondary forests. For. Ecol. Manag. 146, 199–209. https://doi.org/10.1016/S0378-1127(00)00460-6 (2001).Article 

Google Scholar 
Chen, X. G. The biomass and allometric equation of a 20-years-old Cunninghamia lanceolata plantation. Prot. For. Sci. Technol. 4, 28–29, 40. https://doi.org/10.3969/j.issn.1005-5215.2007.04.010.(inChinese) (2007).Article 

Google Scholar 
Wang, X. P. et al. Climatic control of primary forest structure and DBH–height allometry in Northeast China. For. Ecol. Manag. 234, 264–274. https://doi.org/10.1016/j.foreco.2006.07.007 (2006).Article 

Google Scholar 
Peng, C. et al. Developing and evaluating tree height–diameter models at three geographic scales for black spruce in Ontario. N. J. Appl. For. 21(2), 83–92. https://doi.org/10.1093/njaf/21.2.83 (2004).Article 

Google Scholar 
López-Serrano, F. R. et al. Site and weather effects in allometries: A simple approach to climate change effect on pines. For. Ecol. Manag. 215(1–3), 251–270. https://doi.org/10.1016/j.foreco.2005.05.014 (2005).Article 

Google Scholar 
Zhang, C. et al. Developing aboveground biomass equations both compatible with tree volume equations and additive systems for single trees in Poplar plantations in Jiangsu Province, China. Forests 7, 32. https://doi.org/10.3390/f7020032 (2016).Article 

Google Scholar 
Liu, J. C. et al. Comparing non-destructive methods to estimate volume of three tree taxa in Beijing, China. Forests 10, 92. https://doi.org/10.3390/f10020092 (2019).Article 

Google Scholar 
Thangjam, U. et al. Developing tree volume equation for Parkia timoriana grown in home gardens and shifting cultivation areas of North-East India. For. Trees Livelihoods 28(12), 1–13. https://doi.org/10.1080/14728028.2019.1624200 (2019).Article 

Google Scholar 
Dutcă, I. et al. Does slope aspect affect the aboveground tree shape and volume allometry of European Beech (Fagus sylvatica L.) trees?. Forests 13, 1071. https://doi.org/10.3390/f13071071 (2022).Article 

Google Scholar 
Segura, M. & Kanninen, M. Allometric models for tree volume and total aboveground biomass in a tropical humid forest in Costa Rica. Biotropica 37(1), 2–8. https://doi.org/10.2307/30045500 (2005).Article 

Google Scholar 
Wang, X. W. et al. Additive tree biomass equations for Betula platyphylla Suk. plantations in Northeast China. Ann. For. Sci. 75, 60. https://doi.org/10.1007/s13595-018-0738-2 (2018).Article 

Google Scholar 
Niklas, K. J. & Enquist, B. J. Canonical rules for plant organ biomass partitioning and annual allocation. Am. J. Bot. 89(5), 812–819. https://doi.org/10.3732/ajb.89.5.812 (2002).Article 
PubMed 

Google Scholar 
Xiang, W. H. et al. General allometric equations and biomass allocation of Pinus massoniana trees on a regional scale in southern China. Ecol. Res. 26, 697–711. https://doi.org/10.1007/s11284-011-0829-0 (2011).Article 

Google Scholar 
Brown, S. Measuring carbon in forests: Current status and future challenges. Environ. Pollut. 116, 363–372. https://doi.org/10.1016/s0269-7491(01)00212-3 (2002).Article 
CAS 
PubMed 

Google Scholar 
Brassard, B. W. et al. Influence of environmental variability on root dynamics in northern forests. Crit. Rev. Plant Sci. 28, 179–197. https://doi.org/10.1080/07352680902776572 (2009).Article 

Google Scholar 
Montagu, K. D. et al. Developing general allometric relationship for regional estimates of carbon sequestration—An example using Eucalyptus pilularis from seven contrasting sites. For. Ecol. Manag. 204, 113–127. https://doi.org/10.1016/j.foreco.2004.09.003 (2005).Article 

Google Scholar 
Williams, R. J. et al. Allometry for estimating aboveground tree biomass in tropical and subtropical eucalypt woodlands: Towards general predictive equations. Aust. J. Bot. 53, 607–619. https://doi.org/10.1071/BT04149 (2005).Article 

Google Scholar 
Ouimet, R. et al. Estimation of coarse root biomass and nutrient content for sugar maple, jack pine, and black spruce using stem diameter at breast height. Can. J. For. Res. 38, 92–100. https://doi.org/10.1139/x07-134 (2008).Article 

Google Scholar 
Peichl, M. & Arain, M. A. Allometry and partitioning of above-and belowground tree biomass in an age-sequence of white pine forests. For. Ecol. Manag. 253, 68–80. https://doi.org/10.1016/j.foreco.2007.07.003 (2007).Article 

Google Scholar 
Bond-Lamberty, B. et al. Aboveground and below-ground biomass and sapwood area allometric equations for six boreal tree species of northern Manitoba. Can. J. For. Res. 32, 1441–1450. https://doi.org/10.1139/x02-063 (2002).Article 

Google Scholar 
King, J. S. et al. Biomass partitioning in red pine (Pinus resinosa) along a chronosequence in the Upper Peninsula of Michigan. Can. J. For. Res. 37(1), 93–102. https://doi.org/10.1139/x06-217 (2007).Article 

Google Scholar 
Ziania, D. & Mencuccini, M. Aboveground biomass relation-ships for beech (Fagus moesiaca Cz.) trees in Vermio Mountain, northern Greece, and generalised equations for Fagus sp. Ann. For. Sci. 60(5), 439–448. https://doi.org/10.1051/forest:2003036 (2003).Article 

Google Scholar 
Martin, J. G. et al. Aboveground biomass and nitrogen allocation of ten deciduous southern Appalachian tree species. Can. J. For. Res. 28(11), 1648–1659. https://doi.org/10.1139/x98-146 (1998).Article 

Google Scholar 
Wang, C. Biomass allometric equations for 10 co-occurring tree species in Chinese temperate forests. For. Ecol. Manag. 222, 9–16. https://doi.org/10.1016/j.foreco.2005.10.074 (2006).Article 

Google Scholar  More

Darwin C. The structure and distribution of coral reefs, 3rd edn. D. Appleton & Company: New York, NY, USA, 1889.Lajeunesse TC, Parkinson JE, Gabrielson PW, Jeong HJ, Reimer JD, Voolstra CR, et al. Systematic Revision of Symbiodiniaceae Highlights the Antiquity and Diversity of Coral Endosymbionts. Curr Biol. 2018;28:2570–80.e6.CAS 
PubMed 

Google Scholar 
Muscatine L, Porter JW. Reef corals: mutualistic symbioses adapted to nutrient-poor environments. Bioscience. 1977;27:454–60.
Google Scholar 
Rohwer F, Seguritan V, Azam F, Knowlton N. Diversity and distribution of coral-associated bacteria. Mar Ecol Prog Ser. 2002;243:1–10.
Google Scholar 
Rosenberg E, Koren O, Reshef L, Efrony R, Zilber-Rosenberg I. The role of microorganisms in coral health, disease and evolution. Nat Rev Microbiol. 2007;5:355–62.CAS 
PubMed 

Google Scholar 
Muscatine L. The role of symbiotic algae in carbon and energy flux in reef corals. Coral Reefs. 1990;25:75–87.
Google Scholar 
Falkowski PG, Dubinsky Z, Muscatine L, McCloskey L. Population control in symbiotic corals. Bioscience. 1993;43:606–11.
Google Scholar 
Baker DM, Freeman CJ, Wong JCY, Fogel ML, Knowlton N. Climate change promotes parasitism in a coral symbiosis. ISME J. 2018;12:921–30.CAS 
PubMed 
PubMed Central 

Google Scholar 
Rädecker N, Pogoreutz C, Gegner HM, Cárdenas A, Perna G, Geißler L, et al. Heat stress reduces the contribution of diazotrophs to coral holobiont nitrogen cycling. ISME J. 2022;16:1110–8.PubMed 

Google Scholar 
Rädecker N, Pogoreutz C, Voolstra CR, Wiedenmann J, Wild C. Nitrogen cycling in corals: the key to understanding holobiont functioning? Trends Microbiol. 2015;23:490–7.PubMed 

Google Scholar 
Bourne DG, Webster NS. Coral Reef Bacterial Communities. In: Rosenberg E, DeLong EF, editors. The Prokaryotes. Springer: Berlin Heidelberg; 2013. pp. 163–87.Ainsworth DT, Krause L, Bridge T, Torda G, Raina J-B, Zakrzewski M, et al. The coral core microbiome identifies rare bacterial taxa as ubiquitous endosymbionts. ISME J. 2015;9:2261–74.CAS 

Google Scholar 
Pernice M, Raina J-B, Rädecker N, Cárdenas A, Pogoreutz C, Voolstra CR. Down to the bone: the role of overlooked endolithic microbiomes in reef coral health. ISME J. 2020;14:325–34.PubMed 

Google Scholar 
Pogoreutz C, Oakley CA, Rädecker N, Cárdenas A, Perna G, Xiang N, et al. Coral holobiont cues prime Endozoicomonas for a symbiotic lifestyle. ISME J. 2022;16:1883–95.CAS 
PubMed 
PubMed Central 

Google Scholar 
Pogoreutz C, Voolstra CR, Rädecker N, Weis V. The coral holobiont highlights the dependence of cnidarian animal hosts on their associated microbes. In: Bosch TCG, Hadfield MG, editors. Cellular Dialogues in the Holobiont. Boca Raton: CRC Press; 2020. pp. 91–118.Babbin AR, Tamasi T, Dumit D, Weber L, Rodríguez MVI, Schwartz SL, et al. Discovery and quantification of anaerobic nitrogen metabolisms among oxygenated tropical Cuban stony corals. ISME J. 2021;15:1222–35.CAS 
PubMed 

Google Scholar 
Glaze TD, Erler DV, Siljanen HMP. Microbially facilitated nitrogen cycling in tropical corals. ISME J. 2022;16:68–77.CAS 
PubMed 

Google Scholar 
Lesser MP, Morrow KM, Pankey SM, Noonan SHC. Diazotroph diversity and nitrogen fixation in the coral Stylophora pistillata from the Great Barrier Reef. ISME J. 2018;12:813–24.CAS 
PubMed 

Google Scholar 
Cardini U, Bednarz VN, Naumann MS, van Hoytema N, Rix L, Foster RA, et al. Functional significance of dinitrogen fixation in sustaining coral productivity under oligotrophic conditions. Proc R Soc B. 2015;282:20152257.PubMed 
PubMed Central 

Google Scholar 
Pogoreutz C, Rädecker N, Cárdenas A, Gärdes A, Wild C, Voolstra CR. Nitrogen fixation aligns with nifH abundance and expression in two coral trophic functional groups. Front Microbiol. 2017;8:1187.PubMed 
PubMed Central 

Google Scholar 
Pogoreutz C, Rädecker N, Cárdenas A, Gärdes A, Voolstra CR, Wild C. Sugar enrichment provides evidence for a role of nitrogen fixation in coral bleaching. Glob Chang Biol. 2017;23:3838–48.PubMed 

Google Scholar 
Bednarz VN, van de Water JA, Rabouille S, Maguer JF, Grover R, Ferrier‐Pagès C. Diazotrophic community and associated dinitrogen fixation within the temperate coral Oculina patagonica. Environ Microbiol. 2019;21:480–95.CAS 
PubMed 

Google Scholar 
Lema KA, Willis BL, Bourne DG. Corals form characteristic associations with symbiotic nitrogen-fixing bacteria. Appl Environ Microbiol. 2012;78:3136–44.CAS 
PubMed 
PubMed Central 

Google Scholar 
Lema KA, Clode PL, Kilburn MR, Thornton R, Willis BL, Bourne DG. Imaging the uptake of nitrogen-fixing bacteria into larvae of the coral Acropora millepora. ISME J. 2016;10:1804–8.CAS 
PubMed 

Google Scholar 
Santos HF, Carmo FL, Duarte G, Dini-Andreote F, Castro CB, Rosado AS, et al. Climate change affects key nitrogen-fixing bacterial populations on coral reefs. ISME J. 2014;8:2272–9.PubMed 
PubMed Central 

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 USA. 2021;118:e2022653118.PubMed 
PubMed Central 

Google Scholar 
Braker G, Fesefeldt A, Witzel K-P. Development of PCR primer systems for amplification of nitrite reductase genes (nirK and nirS) to detect denitrifying bacteria in environmental samples. Appl Environ Microbiol. 1998;64:3769–75.CAS 
PubMed 
PubMed Central 

Google Scholar 
Tilstra A, El-Khaled YC, Roth F, Rädecker N, Pogoreutz C, Voolstra CR, et al. Denitrification aligns with N2 fixation in Red Sea corals. Sci Rep. 2019;9:1–9.Tilstra A, Roth F, El-Khaled YC, Pogoreutz C, Rädecker N, Voolstra CR, et al. Relative abundance of nitrogen cycling microbes in coral holobionts reflects environmental nitrate availability. R Soc Open Sci. 2021;8:201835.CAS 
PubMed 
PubMed Central 

Google Scholar 
Xiang N, Hassenrück C, Pogoreutz C, Rädecker N, Simancas-Giraldo SM, Voolstra CR, et al. Contrasting microbiome dynamics of putative denitrifying bacteria in two octocoral species exposed to dissolved organic carbon (DOC) and warming. Appl Environ Microbiol. 2022;88:e01886-21.El-Khaled YC, Roth F, Tilstra A, Rädecker N, Karcher DB, Kürten B, et al. In situ eutrophication stimulates dinitrogen fixation, denitrification, and productivity in Red Sea coral reefs. Mar Ecol Prog Ser. 2020;645:55–66.CAS 

Google Scholar 
Beauchamp EG, Trevors JT, Paul JW. Carbon sources for bacterial Denitrification. In: Stewart BA. Advances in Soil Science. Springer: New York, NY; 1989. pp. 113–42.Baker AC. Flexibility and Specificity in Coral-Algal Symbiosis: Diversity, Ecology, and Biogeography of Symbiodinium. Ann Rev Ecol Evol Syst. 2003;34:661–89.
Google Scholar 
Wang J-T, Chen Y-Y, Tew KS, Meng P-J, Chen CA. Physiological and Biochemical Performances of Menthol-Induced Aposymbiotic Corals. PLoS ONE. 2012;7:e46406.CAS 
PubMed 
PubMed Central 

Google Scholar 
Cui G, Liew YJ, Li Y, Kharbatia N, Zahran NI, Emwas A-H, et al. Host-dependent nitrogen recycling as a mechanism of symbiont control in Aiptasia. PLoS Genet. 2019;15:e1008189.CAS 
PubMed 
PubMed Central 

Google Scholar 
Rädecker N, Raina J-B, Pernice M, Perna G, Guagliardo P, Kilburn MR, et al. Using Aiptasia as a Model to Study Metabolic Interactions in Cnidarian-Symbiodinium Symbioses. Front Physiol. 2018;9:214.PubMed 
PubMed Central 

Google Scholar 
Voolstra CR. A journey into the wild of the cnidarian model systemAiptasiaand its symbionts. Mol Ecol. 2013;22:4366–8.PubMed 

Google Scholar 
Sunagawa S, Wilson EC, Thaler M, Smith ML, Caruso C, Pringle JR, et al. Generation and analysis of transcriptomic resources for a model system on the rise: the sea anemone Aiptasia pallida and its dinoflagellate endosymbiont. BMC Genom. 2009;10:258.
Google Scholar 
Xiang T, Hambleton EA, DeNofrio JC, Pringle JR, Grossman AR. Isolation of clonal axenic strains of the symbiotic dinoflagellate Symbiodinium and their growth and host specificity1. J Phycol. 2013;49:447–58.CAS 
PubMed 

Google Scholar 
Thornhill DJ, Lewis AM, Wham DC, Lajeunesse TC. Host‐specialist lineages dominate the adaptive radiation of reef coral endosymbionts. Evolution. 2014;68:352–67.CAS 
PubMed 

Google Scholar 
Bieri T, Onishi M, Xiang T, Grossman AR, Pringle JR. Relative Contributions of Various Cellular Mechanisms to Loss of Algae during Cnidarian Bleaching. PLoS ONE. 2016;11:e0152693.PubMed 
PubMed Central 

Google Scholar 
Baumgarten S, Simakov O, Esherick LY, Liew YJ, Lehnert EM, Michell CT, et al. The genome of Aiptasia, a sea anemone model for coral symbiosis. Proc Natl Acad Sci USA. 2015;112:11893–8.CAS 
PubMed 
PubMed Central 

Google Scholar 
Correa AMS, McDonald MD, Baker AC. Development of clade-specific Symbiodinium primers for quantitative PCR (qPCR) and their application to detecting clade D symbionts in Caribbean corals. Mar Biol. 2009;156:2403–11.CAS 

Google Scholar 
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25:402–8.CAS 
PubMed 

Google Scholar 
Lee JA, Francis CA. DeepnirSamplicon sequencing of San Francisco Bay sediments enables prediction of geography and environmental conditions from denitrifying community composition. Environ Microbiol. 2017;19:4897–912.CAS 
PubMed 

Google Scholar 
Huggett J, Dheda K, Bustin S, Zumla A. Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005;6:279–84.CAS 
PubMed 

Google Scholar 
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.CAS 
PubMed 
PubMed Central 

Google Scholar 
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBO J. 2011;17:10–2.
Google Scholar 
Boutet E, Lieberherr D, Tognolli M, Schneider M, Bairoch A. UniProtKB/Swiss-Prot: the manually annotated section of the UniProt KnowledgeBase. Methods Mol Biol. 2007;406:89–112.Abascal F, Zardoya R, Telford MJ. TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res. 2010;38:7–13.
Google Scholar 
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.CAS 
PubMed 
PubMed Central 

Google Scholar 
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.PubMed 
PubMed Central 

Google Scholar 
Fish JA, Chai B, Wang Q, Sun Y, Brown CT, Tiedje JM, et al. FunGene: the functional gene pipeline and repository. Front Microbiol. 2013;4:291.PubMed 
PubMed Central 

Google Scholar 
Wickham H. ggplot2. Wiley Interdiscip Rev Comput Stat. 2011;3:180–5.Oksanen J, Kindt R, Legendre P, O’Hara B, Stevens MHH, Oksanen MJ, et al. The vegan package. Commun Ecol Package. 2007;10:719.
Google Scholar 
McMurdie PJ, Holmes S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE. 2013;8:e61217.CAS 
PubMed 
PubMed Central 

Google Scholar 
Lin H, Peddada SD. Analysis of compositions of microbiomes with bias correction. Nat Commun. 2020;11:1–11.CAS 

Google Scholar 
Meunier V, Geissler L, Bonnet S, Rädecker N, Perna G, Grosso O, et al. Microbes support enhanced nitrogen requirements of coral holobionts in a high CO 2 environment. Mol Ecol. 2021;30:5888–99.CAS 
PubMed 

Google Scholar 
Geissler L, Meunier V, Rädecker N, Perna G, Rodolfo-Metalpa R, Houlbrèque F, et al. Highly Variable and Non-complex Diazotroph Communities in Corals From Ambient and High CO2 Environments. Front Mar Sci. 2021;8:754682.Thornhill DJ, Xiang Y, Pettay DT, Zhong M, Santos SR. Population genetic data of a model symbiotic cnidarian system reveal remarkable symbiotic specificity and vectored introductions across ocean basins. Mol Ecol. 2013;22:4499–515.CAS 
PubMed 

Google Scholar 
Röthig T, Costa RM, Simona F, Baumgarten S, Torres AF, Radhakrishnan A, et al. Distinct bacterial communities associated with the coral model Aiptasia in aposymbiotic and symbiotic states with Symbiodinium. Front Mar Sci. 2016;3:234.
Google Scholar 
Hartman LM, Blackall LL, van Oppen MJH. Antibiotics reduce bacterial load in Exaiptasia diaphana, but biofilms hinder its development as a gnotobiotic coral model. Access Microbiol. 2022;4:000314.CAS 
PubMed 
PubMed Central 

Google Scholar 
Lawson CA, Raina JB, Kahlke T, Seymour JR, Suggett DJ. Defining the core microbiome of the symbiotic dinoflagellate, Symbiodinium. Environ Microbiol Rep. 2018;10:7–11.CAS 
PubMed 

Google Scholar 
Matthews JL, Raina JB, Kahlke T, Seymour JR, van Oppen MJ, Suggett DJ. Symbiodiniaceae‐bacteria interactions: rethinking metabolite exchange in reef‐building corals as multi‐partner metabolic networks. Environ Microbiol. 2020;22:1675–87.PubMed 

Google Scholar 
Costa RM, Cárdenas A, Loussert-Fonta C, Toullec G, Meibom A, Voolstra CR. Surface Topography, Bacterial Carrying Capacity, and the Prospect of Microbiome Manipulation in the Sea Anemone Coral Model Aiptasia. Front Microbiol. 2021;12:637834.Pelve EA, Fontanez KM, DeLong EF. Bacterial succession on sinking particles in the ocean’s interior. Front Microbiol. 2017;8:2269.PubMed 
PubMed Central 

Google Scholar 
Welles L, Lopez-Vazquez CM, Hooijmans CM, Van Loosdrecht MCM, Brdjanovic D. Prevalence of ‘Candidatus Accumulibacter phosphatis’ type II under phosphate limiting conditions. AMB Express. 2016;6:1–12.Kaneko T. Complete Genomic Sequence of Nitrogen-fixing Symbiotic Bacterium Bradyrhizobium japonicum USDA110. DNA Res. 2002;9:189–97.PubMed 

Google Scholar 
Cziesielski MJ, Liew YJ, Cui G, Schmidt-Roach S, Campana S, Marondedze C, et al. Multi-omics analysis of thermal stress response in a zooxanthellate cnidarian reveals the importance of associating with thermotolerant symbionts. Proc R Soc B: Biol Sci. 2018;285:20172654.
Google Scholar 
Xiang T, Lehnert E, Jinkerson RE, Clowez S, Kim RG, Denofrio JC, et al. Symbiont population control by host-symbiont metabolic interaction in Symbiodiniaceae-cnidarian associations. Nat Commun. 2020;11:1–9.CAS 

Google Scholar  More

All animal care and experimental protocols complied with local and federal laws and guidelines and were approved by the appropriate governing body in Berlin, Germany, the Landesamt fur Gesundheit und Soziales (LaGeSo G-0224/20).Experimental breeding and designThe all-female Amazon molly (Poecilia formosa) is a naturally clonal, live-bearing fish species that gives birth to broods of genetically identical offspring. Like all unisexual vertebrates, Amazon mollies are the result of inter-specific hybridization44,45. As such, this ‘frozen hybrid’ has a heterozygous genome from its ancestral P. mexicana mother and P. latipinna father alleviating concerns about reduced genetic variation and the resulting inbreeding depression often associated with artificially selected isogenic animals. Additionally, despite their clonal nature, the Amazon’s genome shows no evidence of increased mutation accumulation, genomic decay or transposable element activity suggesting the genomes of these animals are evolving in similar ways as sexual species46. They reproduce through gynogenesis where the meiotic process is disrupted so that the eggs contain a full maternal genome. The egg must be fused with a sperm from one of their ancestral species to stimulate embryogenesis, but this paternal DNA is not incorporated into the egg. This provides the opportunity to control when reproduction occurs by controlling the females’ access to male sperm donors.We placed adult females, as potential mothers of experimental fish, in individual (5-gallon) breeding tanks with two Atlantic molly (P. mexicana) males for one week to act as sperm donors. Amazon mollies give birth to broods of generally ~8-30 individuals. A brood is born at once (i.e. all individuals are born within minutes of each other) and birth generally happens early in the day close to dawn. These parental fish were lab-bred and themselves sisters, so of the same age and lineage, and were kept at similar social densities and under standardized environmental conditions throughout their lives to further minimize potential variation in maternal experience. Each breeding tank contained an artificial plant as refuge and was checked frequently each day for the presence of offspring, especially during the morning hours when births are most likely. Newborn mollies were always found in the morning and then singly netted by trained animal caretakers, into individual experimental tanks where their behavior was automatically recorded for the next 70 days (see below). Moving the fish from the maternal tank to the experimental tanks was done in a standardized manner (i.e. individual fish were netted and placed into small dishes of water and then placed in the tracking tanks to limit exposure to the air) by the same caretakers to minimize variation in experience among individual fish. Altogether, eight mothers provided offspring that completed the entire 10-week experiment (Supplementary Table 1).Experimental tanks (27 x 27 cm), made of white Perspex, consisted of four equally sized compartments, and were evenly lit from below using 6500K-LEDs. Environmental conditions were highly standardized across tanks: all tanks were on the same 11:13 (L:D) light schedule, water depth was maintained at 10 cm depth, temperature was maintained at 25 ± 1 °C by a room air conditioning system, and fish received a standardized amount of powdered flake fish food (TetraMin™) twice daily. Opaque blinds surrounded the tanks to further limit outside disturbances. All experimental tanks were connected to the same filtration system where water could mix in the sump tank, allowing chemical cues to be shared across all experimental fish. Previous work has shown exposure to just chemical cues of conspecifics is sufficient in preventing the developmental of pathological behavior that could be associated with development in complete isolation14. We initially placed a total of 40 newborn individuals into the tracking tanks. At the end of the 10-week experiment, we were able to achieve complete tracking data on 26 individuals; camera malfunctions prevented data collection on four individuals, two individuals jumped into neighboring tanks causing the loss of data of all four individuals as we could not verify their identity; four newborn individuals escaped through holes in the water outlet of the tanks; and four individuals died as newborns. All results in the manuscript are on these 26 animals, though including data from all 40 (e.g. patterns of individual variation on the first day post birth) did not change the results or their interpretation (see Supplementary Table 2).Behavioral trackingWe developed a custom recording system using Raspberry Pi computers, which are an upcoming low-cost, highly adaptable solution for many applications in the biological sciences25. Specifically, we created a local network of Raspberry Pi 3B + ’s, each connected to a Raspberry Pi camera positioned exactly above an experimental tank, commanded by a lab computer, and connected to the server on the institute network (Supplementary Fig. 1). We programmed the Raspberry Pi’s using pirecorder26 to take timestamped photos every 3 s across the daily light period, each day, for 10 weeks, and store them automatically in dedicated, automatically named folders on the server. Image settings and resolution were thereby optimized to minimize file size while assuring image quality. After the experimental period, we created videos of all the recorded images of each fish of each day. These videos were subsequently tracked with the Biotracker software27, using background subtraction, providing the x, y coordinates of each fish in each frame. We then processed the data, including scaling and converting the coordinates to mm, and, for each frame, computed fish’s swimming speed (cm/s) and distance from the tank walls (cm). We then summarized these variables both on an hourly and daily basis to compute fish’s median swimming speed, inter-quartile range of swimming speeds, activity (proportion of time spent moving >0.5 cm/s), and median border distance. To quantify fish’s body size over time, we randomly selected five photos per week of each compartment, making sure the fish was away from the compartment walls and did not show strong body curvature, and then used ImageJ software to measure total body length (mm) from the tip of the snout to the end of the body. By averaging the measurements of the five images, we acquired one body size measurement per week.Error checkingWe collected up to 924,000 photos on each individual throughout the experimental period resulting in a total of over 24 million data points collected on our experimental animals (N = 26 individuals). To ensure that our tracking software accurately captured the behavior of our fish, we checked for potential tracking errors in two ways. First, we estimated overall error rates. To do this, we selected at random a starting frame from within a day; then we manually checked each of the subsequent 200 frames and identified whether an error was made (fish was not properly located by BioTracker) or not (fish was properly located) by visual inspection of the videos. We estimated the error rate as the number of errors divided by the total number of checked frames. The overall median error rate over the entire observation period was estimated to be 7%. Error rates increased earlier in the observation period when the fish were smaller (Supplementary Note I). As such, as a second step, we manually went through and corrected all frames for the very first day of tracking (i.e. day 1 post-birth) for all fish (~13,200 frames per individual) as this is a critical time period for one of our research questions. This ensured that the resulting behavioral data were completely accurate for this day. This manual correction allowed us the additional opportunity to compare how well our automatically tracked (i.e. not manually corrected) data performed compared to the manually corrected data. We found that the automatically tracked data re-created near identical estimates of among- and within-individual variance components and most importantly the among-individual correlation between the automatically tracked and manually corrected data was over 0.98 for our behavioral variables (Supplementary Note I). This strongly suggests that any errors introduced by our automated tracking software have minimal influence of our behavioral variables at best and do not affect our interpretation of the results.Statistical analysesWe used linear mixed, or hierarchical, models to partition the behavioral variation across different times periods into its among- and within-individual components. Throughout we focused our analysis on the 26 individuals for which we had complete data for the entire 10-week observation period to ensure comparable variation over time and across models.Our first question of interest was to test when individual differences in behavior first appeared over the course of the experiment. We started by investigating behavior on the first day post birth (Fig. 1A, Supplementary Table 2) and then planned to proceed in a day-by-day fashion until significant repeatability in behavior was apparent (Supplementary Table 3). We used hourly median swimming speed (11 observations for each of 26 individuals) as our response variable and included ‘hour’ and ‘total length (TL)’ as fixed effects and ‘individual’ was included as our random effect of interest. Including TL as a covariate allowed us to test whether behavior was related to an offspring’s body size on its first day of life. We set the first hour of the day as 0 and mean-centered TL as this would allow the among- (and within-) individual variance components to be estimated at these values (i.e. the earliest possible moment from when we could record behavior in the fish). We estimated the adjusted repeatability of median swimming speed as the variance attributable to individual identity over the total variance not explained by the fixed effects. We additionally estimated both marginal and conditional R-squared values which estimate the variance explained by the fixed effects only and the variance explained by the fixed and random effects combined, respectively. As our individual experimental fish came from different mothers, we first explored a number of different variance structures including random intercepts and slopes for both individual ID and maternal ID. This allowed us to test whether maternal identity explained variation in individual behavior. However, the most supported model included random intercepts and slopes for individual ID and not for mother ID, indicating that our methods to reduce variation among mothers were successful (Table 1). We used median swimming speed as our behavioral variable of interest throughout the main manuscript, as this behavior was tightly correlated with most of our other behavioral variables (Supplementary Fig. 2); though results using the other behavioral variables yielded the same interpretation (i.e. that significant individuality in (any) behavior was present on the very first day post-birth; Supplementary Table 2).Our second research question was to investigate how individual behavioral variance changed over the course of the entire observation period (70 days). Again, we first explored several different variance structures to test the importance of maternal identity and/or individual identity on behavioral variation. We found support for the inclusion of random slopes at the individual level, but not maternal level (Table 1). This indicates that levels of among- (and within-) individual variation may differ throughout the observation period. To investigate patterns of change in the variance components, we ran a series of models where we centered the observation covariate on different days. Individual intercepts are estimated when all covariates are set to zero, so this allowed us to ‘slice’ the data to estimate the among- and within-individual variance at different time points over the ten weeks. We ran 11 models as we chose to center the data every 7 days (first model was centered on observation 1; 11th model was centered on observation 70). The predicted individual intercepts (best linear unbiased predictors) and estimated variance components from each model are plotted in Fig. 3.We also closely investigated any potential influence of body size and/or growth rate differences on behavioral expression and individual behavioral variation in this entire 10-week data set. First, we estimated the repeatability of both weekly total length and weekly growth rates to determine if individuals consistently differed in these traits. Then, we ran a series of models with median weekly swimming speed as the response variable and included either weekly total length, weekly growth rate, and/or overall growth rate (estimated over the entire 10 weeks), as our fixed effects of interest. Each model also included the random effects of individual intercepts and slopes. Finally, because body size varies both among individuals (some individuals are on average larger than others) and within individuals (as they grow), we also performed within-individual centering of total length. In this fifth model, we included each individual’s average total length and their weekly deviation from their average length as the two fixed effects of interest. Individual identity and slopes were included as random effects. For all models, we estimated the variance explained by the fixed effects (marginal R2) and the fixed and random effects together (conditional R2). These results are reported in Table 2.For our third and final research question, we tested whether early-life behavior predicted later-life behavior. To test this, we estimated the among-individual correlation (including ‘individual ID’ as our random effect) in behavior using multivariate mixed models where the daily median swimming speeds in each week were the response variables (7 observations per week per individual; 10 weeks total; Fig. 4A). Then to investigate how the strength of these correlations may change over development, we used a linear model to test whether the correlation strength was predicted by the interaction between the first week included in the correlation and distance to the next week in the correlation (1, 2, 3, 4 or 5 weeks away in time; Fig. 4B).All models were performed using Markov Chain Monte Carlo estimation with the MCMCglmm package38 in R v3.6.139. We set our models to run 510,000 iterations with a 10,000 burn-in and thinning every 200 iterations. To ensure proper model mixing and convergence, we initially ran 5 independent chains and inspected posterior trace plots of parameter estimates (Supplementary Note II). In a preliminary analysis we tested three different prior settings (Supplementary Note II); results did not change with prior settings so we chose parameter-expanded priors for all models reported here as these are generally considered to be more robust. An R Markdown file with all the results presented here is included in Supplementary Note II.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More