Philippa Kaur

Canadell, J. G. et al. Multi-decadal increase of forest burned area in Australia is linked to climate change. Nat. Commun. 12, 6921 (2021).Nolan, R. H. et al. What do the Australian Black Summer fires signify for the global fire crisis? Fire 4, 97 (2021).Article 

Google Scholar 
Levin, N., Yebra, M. & Phinn, S. Unveiling the factors responsible for Australia’s Black Summer fires of 2019/2020. Fire 4, 58 (2021).Article 

Google Scholar 
Abram, N. J. et al. Connections of climate change and variability to large and extreme forest fires in southeast Australia. Commun. Earth Environ. 2, 8 (2021).Keenan, R. et al. No evidence that timber harvesting increased the scale or severity of the 2019/20 bushfires in south-eastern Australia. Aust. For. 84, 133–138 (2021).Article 

Google Scholar 
Fire Severity in Harvested Areas (New South Wales Department of Primary Industry, 2020); https://www.dpi.nsw.gov.au/__data/assets/pdf_file/0020/1222391/fire-severity-in-harvested-areas.pdfLindenmayer, D. B., Kooyman, R. M., Taylor, C., Ward, M. & Watson, J. E. Recent Australian wildfires made worse by logging and associated forest management. Nat. Ecol. Evol. 4, 898–900 (2020).Bowman, D. M., Williamson, G. J., Gibson, R. K., Bradstock, R. A. & Keenan, R. J. The severity and extent of the Australia 2019–20 Eucalyptus forest fires are not the legacy of forest management. Nat. Ecol. Evol. 5, 1003–1010 (2021).Poulos, H. M., Barton, A. M., Slingsby, J. A. & Bowman, D. M. Do mixed fire regimes shape plant flammability and post-fire recovery strategies? Fire 1, 39 (2018).Article 

Google Scholar 
Lindenmayer, D. B. et al. Logging elevated the probability of high-severity fire in the 2019–20 Australian forest fires. Nat. Ecol. Evol. https://doi.org/10.1038/s41559-022-01716-z (2022).Peterson, D. A. et al. Australia’s Black Summer pyrocumulonimbus super outbreak reveals potential for increasingly extreme stratospheric smoke events. NPJ Clim. Atmos. Sci. 4, 38 (2021).Gibson, R., Danaher, T., Hehir, W. & Collins, L. A remote sensing approach to mapping fire severity in south-eastern Australia using sentinel 2 and random forest. Remote Sens. Environ. 240, 111702 (2020).Article 

Google Scholar 
Taylor, C., McCarthy, M. A. & Lindenmayer, D. B. Nonlinear effects of stand age on fire severity. Conserv. Lett. 7, 355–370 (2014).Article 

Google Scholar 
Price, O. F. & Bradstock, R. A. The efficacy of fuel treatment in mitigating property loss during wildfires: insights from analysis of the severity of the catastrophic fires in 2009 in Victoria, Australia. J. Environ. Manag. 113, 146–157 (2012).Article 

Google Scholar 
Lindenmayer, D., Taylor, C. & Blanchard, W. Empirical analyses of the factors influencing fire severity in southeastern Australia. Ecosphere 12, e03721 (2021).Article 

Google Scholar 
Bowman, D. M., Williamson, G. J., Prior, L. D. & Murphy, B. P. The relative importance of intrinsic and extrinsic factors in the decline of obligate seeder forests. Global Ecol. Biogeogr. 25, 1166–1172 (2016).Article 

Google Scholar 
Taylor, C., Blanchard, W. & Lindenmayer, D. B. Does forest thinning reduce fire severity in Australian eucalypt forests? Conserv. Lett. 14, e12766 (2021).Article 

Google Scholar 
Lindenmayer, D. B. & Taylor, C. New spatial analyses of Australian wildfires highlight the need for new fire, resource, and conservation policies. Proc. Natl Acad. Sci. USA 117, 12481–12485 (2020).CAS 
Article 

Google Scholar 
Cruz, M., Alexander, M. & Plucinski, M. The effect of silvicultural treatments on fire behaviour potential in radiata pine plantations of South Australia. For. Ecol. Manag. 397, 27–38 (2017).Article 

Google Scholar 
Lindenmayer, D. B., Hobbs, R. J., Likens, G. E., Krebs, C. J. & Banks, S. C. Newly discovered landscape traps produce regime shifts in wet forests. Proc. Natl Acad. Sci. USA 108, 15887–15891 (2011).CAS 
Article 

Google Scholar  More

Obtaining bacterial strains with varied resistance levelsWe experimentally evolved 24 replicate lineages of E. coli to tolerate increasing concentrations of chloramphenicol. By serially passaging bacterial cultures through 14 increasing chloramphenicol levels, we obtained 336 (24 lineages × 14 concentrations) populations of E. coli across a gradient of resistance levels (Fig. 2). The 24 replicate lineages enabled us to study the variability arising from the stochastic nature of mutation acquisition. We refer to these populations as “cultures” rather than “strains” due to the possible coexistence of multiple genotypes.Resistance incurs costs in both thermal tolerance and maximum growth rateWe measured growth rates of experimentally evolved E. coli cultures at three different temperatures: their historic temperature of 37 °C, and the novel temperatures of 32 °C and 42 °C. We hypothesized that growth rate costs of resistance would be larger in the novel temperatures, consistent with reduced thermal niche breadth.Overall, we found the growth rates decreased strongly with increasing antibiotic resistance (Fig. 3A). We then calculated relative growth rates for each lineage by dividing the growth rate at each timepoint by the growth rate of the culture at timepoint 1 (T1) at the appropriate temperature (e.g., all cultures at 32 °C were standardized by the ancestral growth rate at 32 °C). Analysis of these relative growth rates showed that there was both a fitness cost in maximum growth rate and a fitness cost in thermal niche breadth; the linear model showed a strong negative effect of increasing resistance on growth rate at 37 °C (F1, 974 = 988.2, p  More

Oppel, S. et al. Spatial scales of marine conservation management for breeding seabirds. Mar. Policy 98, 37–46 (2018).Article 

Google Scholar 
Lewison, R. et al. Research priorities for seabirds: improving conservation and management in the 21st century. Endanger. Species Res 17, 93–121 (2012).Article 

Google Scholar 
Hasegawa, H. & DeGange, A. R. The Short-tailed Albatross, Diomedea albatrus, its status, distribution and natural history. Am. Birds 36, 806–814 (1982).
Google Scholar 
Tickell, W. L. N. Albatrosses (Pica Press, 2000).BirdLife International. Phoebastria albatrus. The IUCN Red List of Threatened Species, e.T22698335A132642113 https://doi.org/10.2305/IUCN.UK.2018-2.RLTS.T22698335A132642113.en (2018).Japan Ministry of the Environment. Ministry of the Environment Red List (Government of Japan, 2020).COSEWIC. COSEWIC Assessment and Status Report on the Short-tailed Albatross Phoebastria albatrus in Canada (Committee on the Status of Endangered Wildlife in Canada, 2013).Environment Canada. Recovery Strategy for the Short-tailed Albatross (Phoebastria albatrus) and the Pink-footed Shearwater (Puffinus creatopus) in Canada (Environment Canada, 2008).United States of America Fish and Wildlife Service. Endangered and Threatened Wildlife and Plants; Final Rule to List the Short-tailed Albatross as Endangered in the United States. 65 FR 46643, 46643–4654, Document Number 00–19123 (2000).United States of America Fish and Wildlife Service. Short-tailed Albatross (Phoebastria albatrus) 5-Year Review: Summary and Evaluation (United States of America Fish and Wildlife Service, 2020).United States of America Fish and Wildlife Service. Short-tailed Albaross Recovery Plan (United States of America Fish and Wildlife Service, 2008).Orben, R. A. et al. Ontogenetic changes in at-sea distributions of immature short-tailed albatrosses Phoebastria albatrus. Endanger. Species Res 35, 23–37 (2018).Article 

Google Scholar 
Orben, R. A. et al. Across borders: external factors and prior behaviour influence North Pacific albatross associations with fishing vessels. J. Appl. Ecol. 58, 1272–1283 (2021).Article 

Google Scholar 
Fox, C. H., Robertson, C., O’Hara, P. D., Tadey, R. & Morgan, K. H. Spatial assessment of albatrosses, commercial fisheries, and bycatch incidents on Canada’s Pacific coast. Mar. Ecol. Prog. Ser. 672, 205–222 (2021).Article 

Google Scholar 
Piatt, J. F. et al. Predictable hotspots and foraging habitat of the endangered short-tailed albatross (Phoebastria albatrus) in the North Pacific: implications for conservation. Deep Sea Res. Part II 53, 387–398 (2006).Article 

Google Scholar 
Suryan, R. M. et al. Migratory routes of short-tailed albatrosses: use of exclusive economic zones of North Pacific Rim countries and spatial overlap with commercial fisheries in Alaska. Biol. Conserv. 137, 450–460 (2007).Article 

Google Scholar 
Suryan, R. M. & Fischer, K. N. Stable isotope analysis and satellite tracking reveal interspecific resource partitioning of nonbreeding albatrosses off Alaska. Can. J. Zool. 88, 299–305 (2010).CAS 
Article 

Google Scholar 
Zador, S. G., Punt, A. E. & Parrish, J. K. Population impacts of endangered short-tailed albatross bycatch in the Alaskan trawl fishery. Biol. Conserv. 141, 872–882 (2008).Article 

Google Scholar 
Geernaert, T. O., Gilroy, H. L., Kaimmer, S. M., Williams, G. H. & Trumble, R. J. A Feasibility Study that Investigates Options for Monitoring Bycatch of the Short-tailed Albatross in the Pacific Halibut Fishery off Alaska (International Pacific Halibut Commission, 2001).Guy, T. J. et al. Overlap of North Pacific albatrosses with the U.S. west coast groundfish and shrimp fisheries. Fish. Res. 147, 222–234 (2013).Article 

Google Scholar 
Bolnick, D. I. et al. The ecology of individuals: incidence and implications of individual specialization. Am. Natural 161, 1–28 (2003).Article 

Google Scholar 
Votier, S. C. et al. Individual responses of seabirds to commercial fisheries revealed using GPS tracking, stable isotopes and vessel monitoring systems. J. Appl. Ecol. 47, 487–497 (2010).Article 

Google Scholar 
Wakefield, E. D. et al. Long-term individual foraging site fidelity—why some gannets don’t change their spots. Ecology 96, 3058–3074 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
Votier, S. C. et al. Effects of age and reproductive status on individual foraging site fidelity in a long-lived marine predator. Proc. R. Soc. B 284, 20171068 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
Sztukowski, L. A. et al. Sex differences in individual foraging site fidelity of Campbell albatross. Mar. Ecol. Prog. Ser. 601, 227–238 (2018).Article 

Google Scholar 
Gutowsky, S. E. et al. Divergent post-breeding distribution and habitat associations of fledgling and adult Black-footed Albatrosses Phoebastria nigripes in the North Pacific. Ibis 156, 60–72 (2014).Article 

Google Scholar 
Weimerskirch, H., Åkesson, S. & Pinaud, D. Postnatal dispersal of wandering albatrosses Diomedea exulans: implications for the conservation of the species. J. Avian Biol. 37, 23–28 (2006).
Google Scholar 
Olson, S. L. & Hearty, P. J. Probable extirpation of a breeding colony of Short-tailed Albatross (Phoebastria albatrus) on Bermuda by Pleistocene sea-level rise. Proc. Natl Acad. Sci. 100, 12825–12829 (2003).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Dall, W. H. Notes on pre-historic remains in the Aleutian islands. Proc. Calif. Acad. Sci. 4, 283–287 (1872).
Google Scholar 
Eda, M. et al. Inferring the ancient population structure of the vulnerable albatross Phoebastria albatrus, combining ancient DNA, stable isotope, and morphometric analyses of archaeological samples. Conserv. Genet. 13, 143–151 (2012).Article 

Google Scholar 
Cousins, K. L., Dalzell, P. & Gilman, E. Managing pelagic longline-albatross interactions in the North Pacific Ocean. Mar. Ornithol 28, 159–174 (2000).
Google Scholar 
Hobson, K. A. & Montevecchi, W. A. Stable isotopic determinations of trophic relationships of great auks. Oecologia 87, 528–531 (1991).PubMed 
Article 
PubMed Central 

Google Scholar 
Fuller, B. T. et al. Pleistocene paleoecology and feeding behavior of terrestrial vertebrates recorded in a pre-LGM asphaltic deposit at Rancho La Brea, California. Palaeogeogr. Palaeoclimatol. Palaeoecol. 537, 109383 (2020).Article 

Google Scholar 
Hobson, K. A. & Clark, R. G. Assessing avian diets using stable isotopes I: turnover of 13C in tissues. Condor 94, 181–188 (1992).Article 

Google Scholar 
Hyland, C., Scott, M. B., Routledge, J. & Szpak, P. Stable carbon and nitrogen isotope variability of bone collagen to determine the number of isotopically distinct specimens. J. Archaeol. Method Theory https://doi.org/10.1007/s10816-021-09533-7 (2021).Article 

Google Scholar 
Hedges, R. E. M., Clement, J. G., Thomas, D. L. & O’Connell, T. C. Collagen turnover in the adult femoral mid‐shaft: modeled from anthropogenic radiocarbon tracer measurements. Am. J. Phys. Anthropol. 133, 808–816 (2007).PubMed 
Article 
PubMed Central 

Google Scholar 
Guiry, E. J., Orchard, T. J., Royle, T. C. A., Cheung, C. & Yang, D. Y. Dietary plasticity and the extinction of the passenger pigeon (Ectopistes migratorius). Quat. Sci. Rev. 233, 106225 (2020).Article 

Google Scholar 
Minagawa, M. & Wada, E. Stepwise enrichment of 15N along food chains: further evidence and the relation between δ15N and animal age. Geochim. Cosmochim. Acta 48, 1135–1140 (1984).CAS 
Article 

Google Scholar 
DeNiro, M. J. & Epstein, S. Influence of diet on the distribution of carbon isotopes in animals. Geochim. Cosmochim. Acta 42, 495–506 (1978).CAS 
Article 

Google Scholar 
Hobson, K. A., Ambrose, W. G. Jr & Renaud, P. E. Sources of primary production, benthic-pelagic coupling, and trophic relationships within the Northeast Water Polynya: insights from δ13C and δ15N analysis. Mar. Ecol. Prog. Ser. 128, 1–10 (1995).Article 

Google Scholar 
Sigman, D., Karsh, K. & Casciotti, K. Ocean process tracers: nitrogen isotopes in the ocean in Encyclopedia of Ocean Science (eds Steele, J. H. et al.) 4139–4152 (Academic Press, 2009).Guiry, E. Complexities of stable carbon and nitrogen isotope biogeochemistry in ancient freshwater ecosystems: implications for the study of past subsistence and environmental change. Front. Ecol. Evol 7, 313 (2019).Article 

Google Scholar 
Rau, G. H., Takahashi, T. & Des Marais, D. J. Latitudinal variations in plankton δ13C: implications for CO2 and productivity in past oceans. Nature 341, 516–518 (1989).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Popp, B. N. et al. Effect of phytoplankton cell geometry on carbon isotopic fractionation. Geochim. Cosmochim. Acta 62, 69–77 (1998).CAS 
Article 

Google Scholar 
Laws, E. A., Popp, B. N., Bidigare, R. R., Kennicutt, M. C. & Macko, S. A. Dependence of phytoplankton carbon isotopic composition on growth rate and (CO2) aq: theoretical considerations and experimental results. Geochim. Cosmochim. Acta 59, 1131–1138 (1995).CAS 
Article 

Google Scholar 
Vokhshoori, N. L. et al. Broader foraging range of ancient short-tailed albatross populations into California coastal waters based on bulk tissue and amino acid isotope analysis. Mar. Ecol. Prog. Ser. 610, 1–13 (2019).CAS 
Article 

Google Scholar 
Sherwood, O. A., Lehmann, M. F., Schubert, C. J., Scott, D. B. & McCarthy, M. D. Nutrient regime shift in the western North Atlantic indicated by compound-specific δ15N of deep-sea gorgonian corals. Proc. Natl Acad. Sci. 108, 1011–1015 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Szpak, P., Savelle, J. M., Conolly, J. & Richards, M. P. Variation in late holocene marine environments in the Canadian Arctic Archipelago: evidence from ringed seal bone collagen stable isotope compositions. Quat. Sci. Rev. 211, 136–155 (2019).Article 

Google Scholar 
Guiry, E. J. et al. Deforestation caused abrupt shift in Great Lakes nitrogen cycle. Limnol. Oceanogr. 65, 1921–1935 (2020).CAS 
Article 

Google Scholar 
Wiley, A. E. et al. Millennial-scale isotope records from a wide-ranging predator show evidence of recent human impact to oceanic food webs. Proc. Natl Acad. Sci. 110, 8972–8977 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Keeling, C. D. The Suess effect: 13Carbon-14Carbon interrelations. Environ. Int. 2, 229–300 (1979).CAS 
Article 

Google Scholar 
McMahon, K. W., Thorrold, S. R., Elsdon, T. S. & McCarthy, M. D. Trophic discrimination of nitrogen stable isotopes in amino acids varies with diet quality in a marine fish. Limnol. Oceanogr. 60, 1076–1087 (2015).CAS 
Article 

Google Scholar 
Chikaraishi, Y. et al. Determination of aquatic food‐web structure based on compound‐specific nitrogen isotopic composition of amino acids. Limnol. Oceanogr. Methods 7, 740–750 (2009).CAS 
Article 

Google Scholar 
Jackson, A. L., Inger, R., Parnell, A. C. & Bearhop, S. Comparing isotopic niche widths among and within communities: SIBER–stable isotope Bayesian ellipses in R. J. Anim. Ecol. 80, 595–602 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
Ambrose, S. H. Preparation and characterization of bone and tooth collagen for isotopic analysis. J. Archaeol. Sci. 17, 431–451 (1990).Article 

Google Scholar 
Guiry, E. J. & Szpak, P. Improved quality control criteria for stable carbon and nitrogen isotope measurements of ancient bone collagen. J. Archaeol. Sci. 132, 105416 (2021).CAS 
Article 

Google Scholar 
Thompson, D. R. & Furness, R. W. Stable-isotope ratios of carbon and nitrogen in feathers indicate seasonal dietary shifts in Northern Fulmars. Auk 112, 493–498 (1995).Article 

Google Scholar 
Carter, H. R. & Sealy, S. G. Historical occurrence of the short-tailed Albatross in British Columbia and Washington. 1841–1958. Wildl. Afield 11, 24–38 (2014).
Google Scholar 
Crockford, S. The Archaeological History of Short-tailed Albatross in British Columbia: A Review and Summary of STAL Skeletal Remains, as Compared to Other Avian Species, Identified from Historic and Prehistoric Midden Deposits. Report on file, Canadian Wildlife Service (2003).Borrmann, R. M., Phillips, R. A., Clay, T. A. & Garthe, S. High foraging site fidelity and spatial segregation among individual great black-backed gulls. J. Avian Biol. 50, e02156 (2019).Article 

Google Scholar 
Wilkinson, B. P., Haynes-Sutton, A. M., Meggs, L. & Jodice, P. G. High spatial fidelity among foraging trips of Masked Boobies from Pedro Cays, Jamaica. PLoS ONE 15, e0231654 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Araújo, M. S., Bolnick, D. I. & Layman, C. A. The ecological causes of individual specialisation. Ecol. Lett. 14, 948–958 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
Grémillet, D. et al. Offshore diplomacy, or how seabirds mitigate intra-specific competition: a case study based on GPS tracking of Cape gannets from neighbouring colonies. Mar. Ecol. Prog. Ser. 268, 265–279 (2004).Article 

Google Scholar 
Irons, D. B. Foraging area fidelity of individual seabirds in relation to tidal cycles and flock feeding. Ecology 79, 647–655 (1998).Article 

Google Scholar 
Piper, W. H. Making habitat selection more “familiar”: a review. Behav. Ecol. Sociobiol. 65, 1329–1351 (2011).Article 

Google Scholar 
Davoren, G. K., Montevecchi, W. A. & Anderson, J. T. Search strategies of a pursuit‐diving marine bird and the persistence of prey patches. Ecol. Monogr. 73, 463–481 (2003).Article 

Google Scholar 
Hazen, E. L. et al. Marine top predators as climate and ecosystem sentinels. Front. Ecol. Environ. 17, 565–574 (2019).Article 

Google Scholar 
Dall, S. R. X., Bell, A. M., Bolnick, D. I. & Ratnieks, F. L. An evolutionary ecology of individual differences. Ecol. Lett. 15, 1189–1198 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
McAllister, N. M. Avian fauna from the Yuquot excavation in The Yuquot Project, Volume 2 (eds. Folan, W. J. & Dewhirst, J.) 103–174 (National Historic Parks and Sites Branch, 1980).Drucker, P. I. The Northern and Central Nootkan tribes. Bureau of American Ethnology Bulletin 144, 1–480 (1951).
Google Scholar 
Lepofsky, D. & Caldwell, M. Indigenous marine resource management on the Northwest Coast of North America. Ecol. Process 2, 12 (2013).Article 

Google Scholar 
Dewhirst, J. The Indigenous Archaeology of Yuquout, a Nootkan Outside Village (National Historic Parks and Sites Branch, 1980).Longin, R. New method of collagen extraction for radiocarbon dating. Nature 230, 241–242 (1971).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Guiry, E. J. & Hunt, B. P. V. Integrating fish scale and bone isotopic compositions for ‘deep time’ retrospective studies. Mar. Environ. Res. 160, 104982 (2020).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Hobson, K. A., Atwell, L. & Wassenaar, L. I. Influence of drinking water and diet on the stable-hydrogen isotope ratios of animal tissues. Proc. Natl Acad. Sci. 96, 8003–8006 (1999).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Qi, H., Coplen, T. B., Geilmann, H., Brand, W. A. & Böhlke, J. K. Two new organic reference materials for δ13C and δ15N measurements and a new value for the δ13C of NBS 22 oil. Rapid Commun. Mass Spectrom 17, 2483–2487 (2003).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Qi, H. et al. A new organic reference material, l-glutamic acid, USGS41a, for δ13C and δ15N measurements − a replacement for USGS41. Rapid Commun. Mass Spectrom 30, 859–866 (2016).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
Szpak, P., Metcalfe, J. Z. & Macdonald, R. A. Best practices for calibrating and reporting stable isotope measurements in archaeology. J. Archaeol. Sci. Rep 13, 609–616 (2017).
Google Scholar 
Hammer, Ø., Haper, A. T. & Ryan, P. D. PAST: paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4, 4 (2001).
Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019).RStudio Team. RStudio: Integrated Development for R (RStudio, PBC, 2019). More

Identification of two Zn/Co responsive proteins in diatomsZn and Co growth rate experiments in which Zn or Co (omitting the other) were added to the growth media were conducted and harvested for proteomic analysis. Growth rates of the marine diatom species Thalassiosira pseudonana CCMP1335, Phaeodactylum tricornutum CCMP632, Pseudo-nitzschia delicatissima UNC1205 and Chaetoceros sp. RS19 (Chaetoceros RS19 herein) were conducted in a consistent media composition to allow for intercomparison among species (see “Methods”). The onset of growth limitation by Zn and Co was evident by decreased growth rates under low [Zn2+] and [Co2+], and the ability to use Co to restore Zn-limited growth was species-specific and consistent with prior results for the diatoms T. pseudonana, P. tricornutum and P. delicatissima (Fig. 1a, b)9 and for other eukaryotic algae2,8,10. Growth rates of Chaetoceros RS19 were not stimulated by increasing [Co2+] up to 23.5 pM in the absence of added Zn. This inability to substitute Co for Zn in Chaetoceros RS19 was clearly distinct from that of other diatoms, but was consistent with previous observations in Chaetoceros calcitrans10, implying a genus-wide attribute.Fig. 1: Growth responses of diatoms to varying [Zn2+] and [Co2+] and initial detection of ZCRPs in T. pseudonana.Growth rates of four diatoms over a range of a [Zn2+] and b [Co2+]. Data are presented as mean values of biological duplicate cultures. Data is available in Supplementary Table 1. Global proteomic analyses comparing the proteomes of pooled biological duplicate cultures (n = 2) of T. pseudonana in c high vs. low added Zn and d high vs. low added Co. Each point is an identified protein with the mean of technical triplicate abundance scores in one treatment plotted against the mean of abundance scores in another treatment. The solid line denotes 1:1 abundance. Error bars in c are the standard deviation of technical triplicate measurements.Full size imageThe proteome as a function of Zn2+ and Co2+ was explored in the marine diatom T. pseudonana harvested during log phase growth. Global proteomic analysis comparing low (1.1 pM) versus high (10.2 pM) added [Zn2+] and low (2.3 pM) versus high (23.4 pM) added [Co2+] revealed two uncharacterized diatom proteins that greatly increased in abundance at low [Zn2+] or [Co2+] (Fig. 1c, d). These proteins were annotated as a CobW/HypB/UreG, nucleotide binding domain and a bacterial extra-cellular solute binding domain, respectively, within the manually curated JGI Thaps3 T. pseudonana genome17 and were identified in T. pseudonana cultures with high confidence (≥9 exclusive unique peptides, 100% protein probability; Supplementary Fig. 1). BLAST sequence alignments showed these proteins to be homologous with CobW-like proteins (with 31.69% identity relative to Pseudomonas denitrificans CobW) and with the bacterial nickel transport protein NikA (with 30.5% identity relative to E. coli NikA), respectively. Based on their clear response to Zn and Co in the proteomes of multiple diatom species (Fig. 2a–d), the lack of definitive annotations in diatoms, and their genetic distance from bacterial homologs, these proteins are referred to as ZCRP-A and ZCRP-B (Zn/Co Responsive Protein A and B) in this study. Abundance patterns of these proteins were also investigated in P. tricornutum, P. delicatissima and Chaetoceros RS19. ZCRP-A spectral abundance counts were significantly (Kendall correlation, p  10 times. j Topology predictions from five sub-methods (OCTOPUS, Philius, PolyPhobius, SCAMPI, and SPOCTOPUS), consensus prediction (TOPCONS), and predicted ΔG values for P. tricornutum ZCRP-B generated using the TOPCONS webserver (https://topcons.cbr.su.se/)27,28. k Extent of Co uptake after 24 h for wild-type (WT), ZCRPA-knockout (KO), and ZCRPA-overexpression (OE) lines of P. tricornutum normalized to fluorescence units (fsu). Data are presented as mean values ± the standard deviation of biological triplicate cultures (n = 3). Individual data points are overlaid as white circles. The extent of Co uptake was found to be significantly larger in the ZCRPA-OE line compared to the wild-type via one-way ANOVA (f(3) = 23.16, p = 0.000268) and post hoc Dunnett test (p = 0.00048).Full size imageTo date, connections between COG0523 proteins and utilization of Zn and Co have been explored primarily in prokaryotic organisms. For example, the COG0523 protein CobW has a role in vitamin B12 biosynthesis and thus Co use19,21. In contrast, a subgroup of other COG0523 proteins (YjiA, YeiR, ZigA, and ZagA) have been implicated in Zn2+ metabolism8,13,14,15,16, and a client protein to the metallochaperone ZagA in Bacillus subtilis has been identified22.Compared to bacteria, less is known about the function of COG0523 proteins in marine phytoplankton, though COG0523 protein family members are known to occur in all kingdoms8,23. A recent study described the presence of COG0523 domain proteins upregulated under low Zn in the coccolithophore Emiliania huxleyi, but without further functional characterization24, implying a potential Zn-related function of a COG0523 protein in a marine alga distinct from the marine diatoms included in this study.Although various proteins belonging to the COG0523 subgroup share similar conserved domains, they possess different metal binding abilities and thus likely have different functions among the diverse organisms in which they are found. For example, recent work has established that CobW preferentially binds Co2+ as the cognate metal and acts as a Co2+ chaperone ultimately supplying vitamin B12 in bacteria, whereas the closely related putative metal chaperones YeiR and YjiA (homologs of CobW) bind Zn2+19. We can infer from homology and the response to low Zn and low Co in the present study that Zn2+ and Co2+ are likely both cognate metals for diatom ZCRP-A. Further metal binding and affinity assays can confirm and characterize metal binding in this protein.Frustule morphologyPhenotypic plasticity in P. tricornutum is well documented. Two basic cell morphotypes, fusiform and triradiate, are found in natural liquid environments. It is thought that by adopting the triradiate form, a cell increases its surface area and thus the area of membrane available for enzymatic activity or molecular diffusion of dissolved inorganic carbon (DIC) into the cell. The triradiate form is known to be more common under DIC limiting conditions, which supports this hypothesis25. Distinct morphological differences resulted from the knockout (KO) of the ZCRP-A gene. In P. tricornutum, ZCRP-A knockout cells consistently adopted a triradiate shape while wild-type cells were fusiform (Fig. 4i). Normally, triradiate cells of P. tricornutum spontaneously revert to fusiform across generations26, thus it is notable that ZCRP-A knockout cells have consistently maintained their triradiate shape for 4+ years in culture irrespective of media [Zn2+]. As Zn2+ is the predominant metal cofactor used in diatom CAs, the adoption of the triradiate form in knockout P. tricornutum cells may be a response to a disruption of the carbon concentrating mechanism caused by a reduction in Zn acquisition capability due to ZCRP-A knockout. This is consistent with the observed relative increase in Mn2+-utilizing CA (ι-CA) in the knockout line compared to the wild-type (Supplementary Fig. 5).ZCRP-B sequence analysis and cellular localizationUnlike COG0523 proteins, the relationship of ZCRP-B abundance to environmental Zn and Co concentrations does not appear to have been previously described. Topology predictions of P. tricornutum ZCRP-B using TOPCONS27,28 revealed a single predicted transmembrane domain near the N-terminus, with the majority of the protein predicted to be oriented outside the membrane (Fig. 4j). Overexpression and fluorescent tagging of ZCRP-B confirmed localization to the cell membrane (Fig. 4e–h; Supplementary Fig. 3b). A single predicted transmembrane domain contrasts with the Zrt/Irt-like divalent metal transporters (ZIPs) in eukaryotic algae, which have 7+ transmembrane domains and are key Zn transporters in many organisms29,30. It is therefore most likely that ZCRP-B is not a transporter itself, but one part of a multi-protein membrane complex and potentially interacts with the ZIP system. A sequence database similarity search (BLASTp, NCBI) found the ZCRP-B protein to be homologous with NikA, a protein subunit of the bacterial ATP-binding cassette (ABC) type Ni transport system protein Nik (30.5% identity with E. coli NikA, E = 7e−49, Supplementary Fig. 6). This transporter is well characterized in bacteria and is comprised of five subunits NikA-E. NikB and NikC are two pore-forming integral inner membrane proteins, NikD and NikE are two inner membrane-associated proteins with ATPase activity, and NikA is the periplasmic component that functions as the initial metal receptor31. No proteins with homology to NikB nor NikC were detected in the P. tricornutum proteomes generated in this study. Two uncharacterized P. tricornutum proteins were homologous with NikD (28.8% identity, E = 1e−14) and NikE (34.9% identity, E = 1.33e−8), though neither had abundance trends similar to ZCRP-B, implying that their function and regulation are independent of ZCRP-B.The sequence of a functionally similar bacterial ABC transport complex, CntABCDF (cobalt nickel transporter, also known as Opp1) from Staphylococcus aureus was also compared to NikA and ZCRP-B (Supplementary Fig. 6). CntA shares 25.6% identity with ZCRP-B (E = 3e−28), and similar to NikA, is an extra-cytoplasmic solute-binding protein that transports Ni, Zn and Co. CntA functions as a Ni/Co acquisition system in Zn-limited S. aureus32. Although the Nik and Cnt systems serve Ni and Co transport in bacteria, ZCRP-B responds to Zn and Co in marine diatoms, which have a significant Zn demand. This may imply a recruitment and repurposing of this bacterial Ni transporter component as part of the Zn acquisition systems during the evolution of marine diatoms.ZCRP-B as a putative high-affinity ligandSequence similarity to the extracellular transport components NikA and CntA (Supplementary Fig. 6), localization to the plasma membrane (Fig. 4b; Supplementary Fig. 3b), and increased abundance under low Zn and Co conditions (Fig. 2b) of P. tricornutum ZCRP-B suggests a metal-binding role as part of a high-affinity transport complex. The induction of ZCRP-B expression at low [Zn2+] (Fig. 2a–c) fits the description of a high-affinity Zn uptake system observed in marine algae that is known to be induced at low free [Zn2+]33,34, suggesting that this protein is involved in an adaptive response to extremely scarce Zn availability. Furthermore, ZCRP-B could contribute to the pool of high-affinity organic ligands that complex dissolved Zn, either by dissociation from living cells or upon cell death by viral lysis and grazing, in the upper water column12,35.The identification of a membrane-associated Zn-Co responsive protein-containing putative metal-binding sites allows us to reconsider the mechanisms of cellular metal uptake in diatoms. Prior physiological experiments observed Zn uptake in marine diatoms to approach the limits of diffusion33, and predicted kinetic control with fast cell surface metal binding and uptake relative to dissociation and release back to the seawater environment36. To enable this transport capability, it was postulated that transporters might be so abundant that the membrane becomes crowded37. Here, the observation of a putative Zn-binding, membrane-associated protein with only 1 predicted transmembrane domain instead implies a separation of the Zn concentrating function at the cell surface relative to its transport into the cell. In this scenario when Zn is scarce, biosynthesis of ZCRP-B increases and is tethered to the cell surface to compete Zn away from natural dissolved Zn ligands35 and/or chelate Zn atoms that make it through the diffusive boundary layer to the membrane. In this manner, ZCRP-B would increase the surface Zn concentration in the vicinity of Zn transporters, and multiple ZCRP-B proteins could supply nearby surface ZIP transporters or be endocytosed, avoiding the predicted membrane crowding of transporters problem. Aristilde and colleagues have previously demonstrated that weak natural Zn-binding ligands containing cysteine do indeed enhance cellular Zn uptake within the diatom Thalassiosira weissflogii, with heightened effects in Zn-limited compared to Zn-replete cells38. They proposed the formation of a transient tertiary complex between the Zn-bound ligand and Zn transporters (ZIPs and heavy metal P-type ATPases) at the cell surface, which could be mediated by a surface-tethered Zn binding ligand such as ZCRP-B. Future studies could examine the mechanism of Zn exchange between ZCRP-B and Zn/Co transporters such as the ZIPs in eukaryotic algae, which were also detected at lower Zn and Co abundances in P. tricornutum but with relatively lower spectral counts (Supplementary Fig. 7a, b), consistent with this model. Furthermore, the proposed mechanism of ZCRP-B binding is similar to that of the high-affinity Fe3+ binding protein ISIP2a, previously characterized in marine algae as an iron starvation-induced protein39. ISIP2a has been characterized as a phytotransferrin involved in endocytosis-mediated high-affinity Fe uptake in P. tricornutum that acts to concentrate Fe at the cell surface and is an extracellular protein anchored to the membrane with one transmembrane domain39. As the protein sequences of P. tricornutum ZCRP-B and ISIP2a share no significant similarity, it is possible that the uptake mechanism of ZCRP-B is similar to that of ISIP2a, but specific to high-affinity Zn and Co uptake rather than Fe. This suggests a common strategy of using extracellular membrane-anchored metal acquisition proteins in marine algae faced with metal limitation.Co uptake in wild-type and mutant diatom strainsAs ZCRP-A and ZCRP-B abundance is related to media [Co2+] (Fig. 2a–d), we investigated differences in the extent of Co uptake after 24 h among Zn/Co-limited wild-type, ZCRP-A knockout, ZCRP-A overexpression, and ZCRP-B overexpression lines of P. tricornutum via addition of the radiotracer 57Co (see methods). The extent of Co uptake among genetically modified P. tricornutum lines was observed to be significantly different via one-way ANOVA (f(3) = 23.16, p = 0.000268). A Dunnet post hoc test revealed that uptake was significantly greater (2.6× larger) in the ZCRP-A overexpression line compared to wild-type (p = 0.00048, Fig. 4k). We interpret this result as the overexpression of ZCRP-A creating a larger intracellular binding capacity for Co, thus protecting it from intracellular sensor or regulatory systems and/or efflux pumps. In contrast, no significant difference in Co uptake rates was observed when comparing ZCRP-A knockout, ZCRP-B overexpression, and wild-type lines, suggesting that P. tricornutum ZCRP-A knockout cells are capable of compensating for knockout to maintain Co metabolism, perhaps through the use of low-affinity transporters33. This is consistent with these uptake experiments being conducted using seawater media with a relatively abundant concentration of Zn (background of 7.7 pM Co and 4.0 nM Zn in the absence of EDTA), thus the use of low-affinity transporters was likely sufficient to acquire Zn and Co for growth, and neither ZCRP-A knockout nor ZCRP-B overexpression would be expected to add any metabolic benefit (Fig. 4k). Moreover, if ZCRP-B is only one part of a multi-protein acquisition and transport complex as hypothesized, overexpression of the single protein may not result in enhanced functionality.Abundance patterns of CAs in two diatomsCarbonic anhydrase enzymes constitute a major reservoir of Zn and Co within marine diatoms7. Within the stroma, intracellular chloroplastic CAs are essential in supplying CO2 to RUBISCO as they convert HCO3−, the predominant species of inorganic carbon in the pyrenoid, into CO240,41. Seven subclasses of CAs have been identified in marine diatoms to date and are designated as alpha, beta, gamma, delta, zeta, theta, and iota (α, β, γ, δ, ζ, θ, and ι). While Zn2+ is the cofactor most commonly used in algal CAs, utilization of both cadmium (Cd2+) and cobalt (Co2+) in place of Zn2+ at the active site of ζ-CA (CDCA) and a δ-CA, respectively, has been previously documented2,5,42. Overall, Zn-utilizing CAs increased in abundance with increasing Zn, consistent with the need for rapid HCO3− conversion at faster growth rates (Fig. 5; Supplementary Fig. 7). Specifically, spectral abundance counts of two β-type CAs, PtCA1 and PtCA2, became abundant in high [Co2+] (23.4 pM) and [Zn2+] ( > 1.1 pM) and were inversely related to ZCRP-A abundance (Supplementary Fig. 7). Both PtCA1 and PtCA2 are known to localize to the chloroplast pyrenoid41,43. Moreover, the increasing abundance trends of the Zn-utilizing α-CAs (CA-II and CA-VI) and the θ-CA Pt43233, which localize to the periplastidial compartment, chloroplast endoplasmic reticulum, and thylakoid lumen, respectively, at higher and Zn/Co provide further evidence for this strategy of increasing CA use under Zn-replete and higher growth rate conditions (Fig. 5; Supplementary Fig. 7)43,44.Fig. 5: Comparison of α-CA, ι-CA, and ZCRP abundances.Spectral counting abundance scores of a alpha CA, iota CA, and b ZCRP-A and ZCRP-B detected in Zn and Co treatments of P. tricornutum measured by global proteomic analysis. Data are plotted as means ± the standard deviation of technical triplicate measurements of pooled biological duplicate cultures (n = 2). Protein names are shown with their corresponding JGI protein ID.Full size imageIn contrast, abundance trends of the recently discovered ι-CA were inversely related to Zn2+ (Fig. 5). Originally identified in T. pseudonana, ι-CA was found to localize to the inner chloroplast membrane surrounding the stroma and is unusual in that it prefers Mn2+ to Zn2+ as a cofactor45. In the present study, spectral counts of P. tricornutum ι-CA decreased as metal concentrations increased, similar to that observed for ZCRP-A and ZCRP-B (Fig. 5). This ι-CA response was consistent with a Zn sparing strategy under low [Zn2+] and [Co2+] used to prioritize the use of Zn2+ for other metalloenzyme functions.Due to the inverse relationship between the abundances of ZCRP-A and chloroplastic Zn2+-requiring CAs in P. tricornutum (that is, all CAs detected with the exception of ι-CA) and the various types of CAs in T. pseudonana (Supplementary Fig. 7), it seems unlikely ZCRP-A directly interacts with CAs. These results are instead consistent with the hypothesis that ZCRP-A functions as a Zn2+ allocation and prioritization mechanism during Zn limitation. The role of Zn2+ in key transcriptional and translational proteins such as RNA polymerase and ribosomal proteins is well known, and major reservoirs of Zn are associated with these transcription and translation systems in the fast-growing copiotrophic bacterium Pseudoalteromonas6. The availability of Zn in ribosomes and the ER is therefore likely also a cellular priority in diatoms, and could benefit from utilizing the putative chaperone and trafficking capability of ZCRP-A when Zn is scarce. We, therefore, posit that ZCRP-A may serve as a Zn2+ trafficking or storage protein that contributes to the prioritization and movement of Zn2+ to the ER or CER, while the Mn-utilizing Mn ι-CA compensates for the lowered Zn availability in the chloroplast. The increased biosynthesis of ZCRP-A may be an important function to shift Zn homeostasis, competing for intracellular Zn and trafficking it towards the ER or CER.Distribution of putative ZCRP homologs among oceanic taxaPutative ZCRP homologs among eukaryotic oceanic taxa were identified by BLAST searching the P. tricornutum ZCRP-A and ZCRP-B protein sequences against all available transcriptomes in the Marine Microbial Eukaryotic Transcriptome Sequencing Project (MMETSP) database, which includes over 650 assembled and annotated transcriptomes of oceanic microbial eukaryotes46. Phylogenetic analysis revealed the presence of putative ZCRP-A and ZCRP-B homologs in a wide variety of organisms belonging to the Chromista kingdom that could be further categorized into Bacillariophyceae, Dinophyceae, and Prymnesiophyceae classes (Supplementary Figs. 8 and  9). Notably, the Chaetoceros RS-19 ZCRP-A homolog did not phylogenetically cluster with the other diatoms (Bacillariophyceae), but instead appears to be more closely related to E. coli YjiA (Supplementary Fig. 8). Furthermore, the lack of the conserved G2/Switch I region in the Chaetoceros RS-19 homolog (Fig. 3) is anomalous in comparison to other putative homologs identified within the MMETSP database. Overall, ZCRPs are not exclusive to oceanic diatoms, but rather are widely distributed amongst oceanic taxa.Metaproteomic detection of ZCRP-A and ZCRP-BTo investigate the use of ZCRP-A and ZCRP-B in the natural environment, we searched metaproteomic data collected during the KM1128 METZYME (Metals and Enzymes in the Pacific) research expedition on the R/V Kilo Moana October 1–25, 2011 from Oʻahu, Hawaiʻi, to Apia, Samoa (Fig. 6a). dZn followed a nutrient-like distribution as described previously, with an average surface (40 m) dZn concentration of 1.21 nM and average deep water (3000 m) concentration of 10.37 nM47 (Fig. 6b). dCo was highly depleted in the upper photic zone as the result of biological uptake48,49 (Fig. 6c). Eukaryotic homologs of ZCRP-A and ZCRP-B were detected at multiple stations at surface ( More

Bowman, D., Williamson, G. J., Gibson, R. K., Bradstock, R. A. & Keenan, R. J. The severity and extent of the Australia 2019–20 Eucalyptus forest fires are not the legacy of forest management. Nat. Ecol. Evol. 5, 1003–1010 (2021).Article 

Google Scholar 
Lindenmayer, D. B., Kooyman, R., Taylor, C., Ward, M. & Watson, J. Recent Australian wildfires made worse by logging and associated forest management. Nat. Ecol. Evol. 4, 898–900 (2020).Article 

Google Scholar 
Gould, J. S., Knight, I. & Sullivan, A. L. Physical modelling of leaf scorch height from prescribed fires in young Eucalyptus sieberi regrowth forests in South-Eastern Australia. Int. J. Wildl. Fire 7, 7–20 (1997).Article 

Google Scholar 
Keith, D. Ocean Shores to Desert Dunes: The Native Vegetation of NSW and the ACT (Department of Environment and Conservation NSW, 2004).Burrows, N. Predicting canopy scorch height in jarrah forests. CALM Sci. 2, 267–274 (1997).
Google Scholar 
Penney, G., Habibi, D. & Cattani, M. Firefighter tenability and its influence on wildfire suppression. Fire Saf. J. 106, 38–51 (2019).Article 

Google Scholar 
Sharples, J. J. et al. Natural hazards in Australia: extreme bushfire. Clim. Change 139, 85–99 (2016).Article 

Google Scholar 
Attiwill, P. M. et al. Timber harvesting does not increase fire risk and severity in wet eucalypt forests of Southern Australia. Conserv. Lett. 7, 341–354 (2014).Article 

Google Scholar 
Lindenmayer, D., Taylor, C. & Blanchard, W. Empirical analyses of the factors influencing fire severity in southeastern Australia. Ecosphere 12, e03721 (2021).Article 

Google Scholar 
Taylor, C., Blanchard, W. & Lindenmayer, D. B. Does forest thinning reduce fire severity in Australian eucalypt forests? Conserv. Lett. 14, e12766 (2020).
Google Scholar 
Taylor, C., McCarthy, M. A. & Lindenmayer, D. B. Non-linear effects of stand age on fire severity. Conserv. Lett. 7, 355–370 (2014).Article 

Google Scholar 
Furlaud, J. M., Prior, L. D., Williamson, G. J. & Bowman, D. M. J. S. Fire risk and severity decline with stand development in Tasmanian giant Eucalyptus forest. For. Ecol. Manag. 502, 119724 (2021).Article 

Google Scholar 
Price, O. F. & Bradstock, R. A. The efficacy of fuel treatment in mitigating property loss during wildfires: insights from analysis of the severity of the catastrophic fires in 2009 in Victoria, Australia. J. Environ. Manag. 113, 146–157 (2012).Article 

Google Scholar 
Taylor, C. & Lindenmayer, D. B. The adequacy of Victoria’s protected areas for conserving its forest-dependent fauna. Austral Ecol. 44, 1076–1090 (2019).Article 

Google Scholar 
Taylor, C., Blanchard, W. & Lindenmayer, D. B. What are the relationships between thinning and fire severity? Austral Ecol. https://doi.org/10.1111/aec.13096 (2021).La Sala, A. Thinning Regrowth Eucalypts Native Forest Silviculture Technical Bulletin No. 13 (Forestry Tasmania, 2001).Cary, G. J., Blanchard, W., Foster, C. N. & Lindenmayer, D. B. Effects of altered fire intervals on critical timber production and conservation values. Int. J. Wildl. Fire 30, 322–328 (2021).Article 

Google Scholar 
Filkov, A. I. et al. The determinants of crown fire runs during extreme wildfires in broadleaf forests in Australia. Adv. For. Fire Res. https://doi.org/10.14195/978-989-26-16-506_190; http://hdl.handle.net/10316.2/44517 (2018).Lindenmayer, D. B., Hobbs, R. J., Likens, G. E., Krebs, C. J. & Banks, S. C. Newly discovered landscape traps produce regime shifts in wet forests. Proc. Natl Acad. Sci. USA 108, 15887–15891 (2011).CAS 
Article 

Google Scholar  More

Boyle, T. H. & Anderson, E. in Cacti: Biology and Uses (ed. Nobel, P. S.) 125–141 (Univ. California Press, 2002).Gibson, A. C. & Nobel, P. S. The Cactus Primer (Harvard Univ. Press, 1986).Bravo Hollis, H. & Sánchez Mejorada, H. Las Cactáceas de México (Univ. Nacional Autónoma de México, 1978).Goettsch, B. et al. High proportion of cactus species threatened with extinction. Nat. Plants 1, 15142 (2015).CAS 
PubMed 

Google Scholar 
Benavides, E., Breceda, A. & Anadón, J. D. Winners and losers in the predicted impact of climate change on cacti species in Baja California. Plant Ecol. 222, 29–44 (2021).
Google Scholar 
Nobel, P. S. Responses of some North American CAM plants to freezing temperatures and doubled CO2 concentrations: implications of global climate change for extending cultivation. J. Arid. Environ. 34, 187–196 (1996).
Google Scholar 
Reyes-García, C. & Andrade, J. L. Crassulacean acid metabolism under global climate change. N. Phytol. 181, 754–757 (2009).
Google Scholar 
Smith, S. D., Didden-Zopfy, B. & Nobel, P. S. High-temperature responses of North American cacti. Ecology 65, 643–651 (1984).
Google Scholar 
Larios, E., González, E. J., Rosen, P. C., Pate, A. & Holm, P. Population projections of an endangered cactus suggest little impact of climate change. Oecologia 192, 439–448 (2020).PubMed 

Google Scholar 
Esparza-Olguı́n, L., Valverde, T. & Vilchis-Anaya, E. Demographic analysis of a rare columnar cactus (Neobuxbaumia macrocephala) in the Tehuacan Valley, Mexico. Biol. Conserv. 103, 349–359 (2002).
Google Scholar 
Seal, C. E. et al. Thermal buffering capacity of the germination phenotype across the environmental envelope of the Cactaceae. Glob. Change Biol. 23, 5309–5317 (2017).
Google Scholar 
Huang, J., Yu, H., Guan, X., Wang, G. & Guo, R. Accelerated dryland expansion under climate change. Nat. Clim. Change 6, 166–171 (2016).
Google Scholar 
Gurvich, D. E. et al. Combined effect of water potential and temperature on seed germination and seedling development of cacti from a mesic Argentine ecosystem. Flora 227, 18–24 (2017).
Google Scholar 
Nuzhyna, N., Baglay, K., Golubenko, A. & Lushchak, O. Anatomically distinct representatives of Cactaceae Juss. family have different response to acute heat shock stress. Flora 242, 137–145 (2018).
Google Scholar 
Andrade, J. L. & Nobel, P. S. Microhabitats and water relations of epiphytic cacti and ferns in a lowland neotropical forest. Biotropica 29, 261–270 (1997).
Google Scholar 
Williams, D. G., Hultine, K. R. & Dettman, D. L. Functional trade-offs in succulent stems predict responses to climate change in columnar cacti. J. Exp. Bot. 65, 3405–3413 (2014).PubMed 

Google Scholar 
Aragón-Gastélum, J. L. et al. Induced climate change impairs photosynthetic performance in Echinocactus platyacanthus, an especially protected Mexican cactus species. Flora Morphol. Distrib. Funct. Ecol. Plants 209, 499–503 (2014).
Google Scholar 
Martorell, C., Montañana, D. M., Ureta, C. & Mandujano, M. C. Assessing the importance of multiple threats to an endangered globose cactus in Mexico: cattle grazing, looting and climate change. Biol. Conserv. 181, 73–81 (2015).
Google Scholar 
Dávila, P., Téllez, O. & Lira, R. Impact of climate change on the distribution of populations of an endemic Mexican columnar cactus in the Tehuacán-Cuicatlán Valley, Mexico. Plant Biosyst. 147, 376–386 (2013).
Google Scholar 
Conver, J. L., Foley, T., Winkler, D. E. & Swann, D. E. Demographic changes over >70 yr in a population of saguaro cacti (Carnegiea gigantea) in the northern Sonoran Desert. J. Arid. Environ. 139, 41–48 (2017).
Google Scholar 
Carrillo-Angeles, I. G., Suzán-Azpiri, H., Mandujano, M. C., Golubov, J. & Martínez-Ávalos, J. G. Niche breadth and the implications of climate change in the conservation of the genus Astrophytum (Cactaceae). J. Arid. Environ. 124, 310–317 (2016).
Google Scholar 
de Cavalcante, A. M. B. & de Duarte, A. S. Modeling the distribution of three cactus species of the Caatinga biome in future climate scenarios. Int. J. Ecol. Environ. Sci. 45, 191–203 (2019).
Google Scholar 
de Cavalcante, A. M. B., de Duarte, A. S. & Ometto, J. P. H. B. Modeling the potential distribution of Epiphyllum phyllanthus (L.) Haw. under future climate scenarios in the Caatinga biome. An. Acad. Bras. Cienc. 92, 351–358 (2020).
Google Scholar 
Tellez-Valdes, O. & DiVila-Aranda, P. Protected areas and climate change: a case study of the cacti in the Tehuacan-Cuicatlan biosphere reserve, Mexico. Conserv. Biol. 17, 846–853 (2003).
Google Scholar 
dos Santos Simões, S., Zappi, D., da Costa, G. M., de Oliveira, G. & Aona, L. Y. S. Spatial niche modelling of five endemic cacti from the Brazilian Caatinga: past, present and future. Austral Ecol. 45, 1–13 (2019).
Google Scholar 
Gorostiague, P., Sajama, J. & Ortega-Baes, P. Will climate change cause spatial mismatch between plants and their pollinators? A test using Andean cactus species. Biol. Conserv. 226, 247–255 (2018).
Google Scholar 
Butler, C. J., Wheeler, E. A. & Stabler, L. B. Distribution of the threatened lace hedgehog cactus (Echinocereus reichenbachii) under various climate change scenarios. J. Torre. Bot. Soc. 139, 46–55 (2012).
Google Scholar 
Johnson, C. N. Species extinction and the relationship between distribution and abundance. Nature 394, 272–274 (1998).CAS 

Google Scholar 
Thuiller, W., Lavorel, S. & Araújo, M. B. Niche properties and geographical extent as predictors of species sensitivity to climate change. Glob. Ecol. Biogeogr. 14, 347–357 (2005).
Google Scholar 
Enquist, B. J. Cyberinfrastructure for an integrated botanical information network to investigate the ecological impacts of global climate change on plant biodiversity. Preprint at PeerJ https://doi.org/10.7287/peerj.preprints.2615v2 (2016).Buisson, L., Thuiller, W., Casajus, N., Lek, S. & Grenouillet, G. Uncertainty in ensemble forecasting of species distribution. Glob. Change Biol. 16, 1145–1157 (2010).
Google Scholar 
Thuiller, W., Guéguen, M., Renaud, J., Karger, D. N. & Zimmermann, N. E. Uncertainty in ensembles of global biodiversity scenarios. Nat. Commun. 10, 1446 (2019).PubMed 
PubMed Central 

Google Scholar 
Goettsch, B., Durán, A. P. & Gaston, K. J. Global gap analysis of cactus species and priority sites for their conservation. Conserv. Biol. 33, 369–376 (2018).PubMed 

Google Scholar 
Maitner, B. S. et al. The bien R package: A tool to access the Botanical Information and Ecology Network (BIEN) database. Methods Ecol. Evol. 9, 373–379 (2018).
Google Scholar 
Karger, D. N. et al. Climatologies at high resolution for the Earth’s land surface areas. Sci. Data 4, 170122 (2017).PubMed 
PubMed Central 

Google Scholar 
Sanderson, B. M., Knutti, R. & Caldwell, P. A representative democracy to reduce interdependency in a multimodel ensemble. J. Clim. 28, 5171–5194 (2015).
Google Scholar 
Brodzik, M. J., Billingsley, B., Haran, T., Raup, B. & Savoie, M. H. EASE-Grid 2.0: Incremental but significant improvements for Earth-gridded data sets. ISPRS Int. J. Geo-Inf. 1, 32–45 (2012).
Google Scholar 
Venter, O. et al. Global terrestrial human footprint maps for 1993 and 2009. Sci. Data 3, 160067 (2016).PubMed 
PubMed Central 

Google Scholar 
Phillips, S. maxnet: Fitting ‘maxent’ species distribution models with ‘glmnet’. R package version 0.1.4. https://CRAN.R-project.org/package=maxnet (2017).Friedman, J., Hastie, T. & Tibshirani, R. Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33, 1–22 (2010).PubMed 
PubMed Central 

Google Scholar 
Dormann, C. F. et al. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography 36, 27–46 (2013).
Google Scholar 
Franklin, S. B., Gibson, D. J., Robertson, P. A., Pohlmann, J. T. & Fralish, J. S. Parallel analysis: a method for determining significant principal components. J. Veg. Sci. 6, 99–106 (1995).
Google Scholar 
Roberts, D. R. et al. Cross-validation strategies for data with temporal, spatial, hierarchical, or phylogenetic structure. Ecography 40, 913–929 (2017).
Google Scholar 
Merow, C., Smith, M. J. & Silander, J. A. A practical guide to MaxEnt for modeling species’ distributions: what it does, and why inputs and settings matter. Ecography 36, 1058–1069 (2013).
Google Scholar 
Allouche, O., Tsoar, A. & Kadmon, R. Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 43, 1223–1232 (2006).
Google Scholar 
Calabrese, J. M., Certain, G., Kraan, C. & Dormann, C. F. Stacking species distribution models and adjusting bias by linking them to macroecological models. Glob. Ecol. Biogeogr. 23, 99–112 (2014).
Google Scholar 
R Core Team R: A Language and Environment for Statistical Computing Version 3.6.0 (R Foundation for Statistical Computing, 2019). https://www.R-project.org/ More

Xu, Y. et al. Annual oil palm plantation maps in Malaysia and Indonesia from 2001 to 2016. Earth Syst. Sci. Data 12, 847–867 (2020).Article 

Google Scholar 
Meijaard, E. et al. The environmental impacts of palm oil in context. Nat. Plants 6, 1418–1426 (2020).Article 

Google Scholar 
Guillaume, T. et al. Carbon costs and benefits of Indonesian rainforest conversion to plantations. Nat. Commun. 9, 2388 (2018).Article 

Google Scholar 
Ordway, E. M. & Asner, G. P. Carbon declines along tropical forest edges correspond to heterogeneous effects on canopy structure and function. Proc. Natl Acad. Sci. USA 117, 7863–7870 (2020).CAS 
Article 

Google Scholar 
Hansen, M. C. et al. High-resolution global maps of 21st-century forest cover change. Science 342, 850 (2013).CAS 
Article 

Google Scholar 
Santoro, M. et al. The global forest above-ground biomass pool for 2010 estimated from high-resolution satellite observations. Earth Syst. Sci. Data 13, 3927–3950 (2021).Article 

Google Scholar 
The World Database on Protected Areas (WDPA) (UNEP-WCMC and IUCN, accessed 12 February 2020); www.protectedplanet.netMahmud, A., Rehrig, M. & Hills, G. Improving the Livelihoods of Palm Oil Smallholders: The Role of the Private Sector (FSG, 2010).Lasco, R. Forest carbon budgets in Southeast Asia following harvesting and land cover change. Sci. China 45, 55–64 (2002).Article 

Google Scholar 
Historical Greenhouse Gas Emissions (Climate Watch, accessed 6 October 2021); https://www.climatewatchdata.org/Euler, M., Schwarze, S., Siregar, H. & Qaim, M. Oil palm expansion among smallholder farmers in Sumatra, Indonesia. J. Agric. Econ. 67, 658–676 (2016).Article 

Google Scholar 
Donofrio, S., Rothrock, P. & Leonard, J. J. F. T. Supply Change: Tracking Corporate Commitments to Deforestation-free SupplyChains, 2017 (Forest Trends, 2017).Rist, L., Feintrenie, L. & Levang, P. The livelihood impacts of oil palm: smallholders in Indonesia. Biodivers. Conserv. 19, 1009–1024 (2010).Article 

Google Scholar 
Saadun, N. et al. Socio-ecological perspectives of engaging smallholders in environmental-friendly palm oil certification schemes. Land Use Policy 72, 333–340 (2018).Article 

Google Scholar 
Hansen, M. C., Stehman, S. V. & Potapov, P. V. Quantification of global gross forest cover loss. Proc. Natl Acad. Sci. USA 107, 8650 (2010).CAS 
Article 

Google Scholar 
Santoro, M. & Cartus, O. ESA Biomass Climate Change Initiative (Biomass_cci): Global datasets of forest above-ground biomass for the year 2017 v.1 (Centre for Environmental Data Analysis, 2019); https://doi.org/10.5285/bedc59f37c9545c981a839eb552e4084Busch, J. et al. Reductions in emissions from deforestation from Indonesia’s moratorium on new oil palm, timber, and logging concessions. Proc. Natl Acad. Sci. USA 112, 1328–1333 (2015).CAS 
Article 

Google Scholar 
McGarigal, K., Cushman, S. A. & Ene, E. FRAGSTATS v.4: spatial pattern analysis program for categorical and continuous maps (Univ. Massachusetts, 2012). More

Arneth, A. et al. Terrestrial biogeochemical feedbacks in the climate system. Nat. Geosci. 3, 525–532 (2010).CAS 
Article 

Google Scholar 
Green, J. K. et al. Regionally strong feedbacks between the atmosphere and terrestrial biosphere. Nat. Geosci. 10, 410–414 (2017).CAS 
Article 

Google Scholar 
Heimann, M. & Reichstein, M. Terrestrial ecosystem carbon dynamics and climate feedbacks. Nature 451, 289–292 (2008).CAS 
Article 

Google Scholar 
Beer, C. et al. Temporal and among-site variability of inherent water use efficiency at the ecosystem level. Glob. Biogeochem. Cycles 23, 1–13 (2009).Article 

Google Scholar 
Keenan, T. F. et al. Increase in forest water-use efficiency as atmospheric carbon dioxide concentrations rise. Nature 499, 324–327 (2013).CAS 
Article 

Google Scholar 
Frank, D. C. et al. Water-use efficiency & transpiration across European forests during the Anthropocene. Nat. Clim. Change 5, 579–583 (2015).CAS 
Article 

Google Scholar 
Mastrotheodoros, T. et al. Linking plant functional trait plasticity and the large increase in forest water use efficiency. J. Geophys. Res. Biogeosci. 122, 2393–2408 (2017).Article 

Google Scholar 
Lavergne, A. et al. Observed and modelled historical trends in the water-use efficiency of plants and ecosystems. Glob. Change Biol. 25, 2242–2257 (2019).Article 

Google Scholar 
Huxman, T. E. et al. Convergence across biomes to a common rain-use efficiency. Nature 429, 651–654 (2004).CAS 
Article 

Google Scholar 
Yang, Y. et al. Contrasting responses of water use efficiency to drought across global terrestrial ecosystems. Sci. Rep. 6, 23284 (2016).CAS 
Article 

Google Scholar 
Huang, L. et al. A global examination of the response of ecosystem water-use efficiency to drought based on MODIS data. Sci. Total Environ. 601–602, 1097–1107 (2017).Article 

Google Scholar 
Reichstein, M. et al. Severe drought effects on ecosystem CO2 and H2O fluxes at three Mediterranean evergreen sites: revision of current hypotheses? Glob. Change Biol. 8, 999–1017 (2002).Article 

Google Scholar 
Reichstein, M. et al. Inverse modeling of seasonal drought effects on canopy CO2/H2O exchange in three Mediterranean ecosystems. J. Geophys. Res. Atmos. 108, 4726 (2003).Article 

Google Scholar 
Cooley, S. S. et al. Assessing regional drought impacts on vegetation and evapotranspiration: a case study in Guanacaste, Costa Rica. Ecol. Appl. 29, e01834 (2019).Article 

Google Scholar 
Medrano, H., Flexas, J. & Galmés, J. Variability in water use efficiency at the leaf level among Mediterranean plants with different growth forms. Plant Soil 317, 17–29 (2008).Article 

Google Scholar 
Soh, W. K. et al. Rising CO2 drives divergence in water use efficiency of evergreen and deciduous plants. Sci. Adv. 5, eaax7906 (2019).CAS 
Article 

Google Scholar 
Wang, M., Chen, Y., Wu, X. & Bai, Y. Forest-type-dependent water use efficiency trends across the northern hemisphere. Geophys. Res. Lett. 45, 8283–8293 (2018).Article 

Google Scholar 
Enquist, B. et al. Scaling from traits to ecosystems: developing a general trait driver theory via integrating trait-based and metabolic scaling theories. Adv. Ecol. Res. 52, 249–318 (2015).Article 

Google Scholar 
Gross, N. et al. Functional trait diversity maximizes ecosystem multifunctionality. Nat. Ecol. Evol. 1, 0132 (2017).Article 

Google Scholar 
Bagousse‐Pinguet, Y. L. et al. Testing the environmental filtering concept in global drylands. J. Ecol. 105, 1058–1069 (2017).Article 

Google Scholar 
Ponce Campos, G. E. et al. Ecosystem resilience despite large-scale altered hydroclimatic conditions. Nature 494, 349–352 (2013).CAS 
Article 

Google Scholar 
Fisher, J. B. et al. The future of evapotranspiration: global requirements for ecosystem functioning, carbon and climate feedbacks, agricultural management, and water resources. Water Resour. Res. 53, 2618–2626 (2017).Article 

Google Scholar 
Xue, B.-L. et al. Global patterns, trends, and drivers of water use efficiency from 2000 to 2013. Ecosphere 6, art174 (2015).Article 

Google Scholar 
Fisher, J. B. et al. ECOSTRESS: NASA’s next generation mission to measure evapotranspiration from the International Space Station. Water Resour. Res. 56, e2019WR026058 (2020).Article 

Google Scholar 
Higgins, M. A. et al. Geological control of floristic composition in Amazonian forests. J. Biogeogr. 38, 2136–2149 (2011).Article 

Google Scholar 
De Kauwe, M. G., Keenan, T. F., Medlyn, B. E., Prentice, I. C. & Terrer, C. Satellite based estimates underestimate the effect of CO2 fertilization on net primary productivity. Nat. Clim. Change 6, 892–893 (2016).Article 

Google Scholar 
Huang, M. et al. Seasonal responses of terrestrial ecosystem water-use efficiency to climate change. Glob. Change Biol. 22, 2165–2177 (2016).Article 

Google Scholar 
Lin, Y.-S. et al. Optimal stomatal behaviour around the world. Nat. Clim. Change 5, 459–464 (2015).CAS 
Article 

Google Scholar 
Medlyn, B. E. et al. Reconciling the optimal and empirical approaches to modelling stomatal conductance. Glob. Change Biol. 17, 2134–2144 (2011).Article 

Google Scholar 
Peters, W. et al. Increased water-use efficiency and reduced CO2 uptake by plants during droughts at a continental scale. Nat. Geosci. 11, 744–748 (2018).CAS 
Article 

Google Scholar 
Cheng, L. et al. Recent increases in terrestrial carbon uptake at little cost to the water cycle. Nat. Commun. 8, 110 (2017).Article 

Google Scholar 
Fisher, J. B. ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS): Level-3 Evapotranspiration L3(ET_PT-JPL) Algorithm Theoretical Basis Document. Jet Propulsion Laboratory, California Institute of Technology (2018).Running, S. W. et al. A continuous satellite-derived measure of global terrestrial primary production. BioScience 54, 547–560 (2004).Article 

Google Scholar 
Heinsch, F. et al. Evaluation of remote sensing based terrestrial productivity from MODIS using regional tower eddy flux network observations. IEEE Trans. Geosci. Remote Sens. 44, 1908–1925 (2006).Article 

Google Scholar 
Zhao, M., Heinsch, F., Nemani, R. & Running, S. Improvements of the MODIS terrestrial gross and net primary production global data set. Remote Sens. Environ. 95, 164–176 (2005).Article 

Google Scholar 
Ryu, Y. et al. Integration of MODIS land and atmosphere products with a coupled-process model to estimate gross primary productivity and evapotranspiration from 1 km to global scales. Glob. Biogeochem. Cycles 25, GB4017 (2011).Article 

Google Scholar  More