Ecology

Bonan, G. B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 320, 1444–1449 (2008).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Frey, S. J. K. et al. Spatial models reveal the microclimatic buffering capacity of old-growth forests. Sci. Adv. 2, e1501392. https://doi.org/10.1126/sciadv.1501392 (2016).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Davis, K. T., Dobrowski, S. Z., Holden, Z. A., Higuera, P. E. & Abatzoglou, J. T. Microclimatic buffering in forests of the future: The role of local water balance. Ecography 42, 1–11 (2019).Article 

Google Scholar 
de Frenne, P. et al. Global buffering of temperatures under forest canopies. Nat. Ecol. Evol. 3, 744–749 (2019).Article 
PubMed 

Google Scholar 
Zellweger, F. et al. Forest microclimate dynamics drive plant responses to warming. Science 368, 772–775 (2020).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Rahman, M. A., Moser, A., Rötzer, T. & Pauleit, S. Microclimatic differences and their influence on transpirational cooling of Tilia cordata in two contrasting street canyons in Munich, Germany. Agric. For. Meteorol. 232, 443–456 (2017).Article 
ADS 

Google Scholar 
Rahman, M. A., Moser, A., Rötzer, T. & Pauleit, S. Within canopy temperature differences and cooling ability of Tilia cordata trees grown in urban conditions. Build. Environ. 114, 118–128 (2017).Article 

Google Scholar 
Ehbrecht, M., Schall, P., Ammer, C., Fischer, M. & Seidel, D. Effects of structural heterogeneity on the diurnal temperature range in temperate forest ecosystems. For. Ecol. Manag. 432, 860–867 (2019).Article 

Google Scholar 
Richter, R., Hutengs, C., Wirth, C., Bannehr, L. & Vohland, M. Detecting tree species effects on forest canopy temperatures with thermal remote sensing: The role of spatial resolution. Remote Sens. 13, 135. https://doi.org/10.3390/rs13010135 (2021).Article 
ADS 

Google Scholar 
IPCC. Climate change 2021: The physical science basis. In Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds Masson-Delmotte, V. et al.) (Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, In Press, 2021).Rahman, M. A. et al. Traits of trees for cooling urban heat islands: A meta-analysis. Build. Environ. 170, 106606. https://doi.org/10.1016/j.buildenv.2019.106606 (2020).Article 

Google Scholar 
Rahman, M. A., Moser, A., Rötzer, T. & Pauleit, S. Comparing the transpirational and shading effects of two contrasting urban tree species. Urban Ecosyst. 22, 683–697 (2019).Article 

Google Scholar 
Joly, F.-X. et al. Tree species diversity affects decomposition through modified micro-environmental conditions across European forests. New Phytol. 214, 1281–1293 (2017).Article 
CAS 
PubMed 

Google Scholar 
Lindo, Z. & Winchester, N. Out on a limb: microarthropod and microclimate variation in coastal temperate rainforest canopies. Insect Conserv. Divers. 6, 513–521 (2013).Article 

Google Scholar 
Pincebourde, S., Murdock, C. C., Vickers, M. & Sears, M. W. Fine-scale microclimatic variation can shape the responses of organisms to global change in both natural and urban environments. Integr. Comp. Biol. 56, 45–61 (2016).Article 
PubMed 

Google Scholar 
Janssen, P., Fuhr, M. & Bouget, C. Beyond forest habitat qualities: Climate and tree characteristics as the major drivers of epiphytic macrolichen assemblages in temperate mountains. J. Veg. Sci. 30, 42–54 (2019).Article 

Google Scholar 
Welti, E. A. R. et al. Temperature drives variation in flying insect biomass across a German malaise trap network. Insect Conserv. Divers. https://doi.org/10.1111/icad.12555 (2021).Article 

Google Scholar 
Lin, Y.-S., Medlyn, B. E. & Ellsworth, D. S. Temperature responses of leaf net photosynthesis: The role of component processes. Tree Physiol. 32, 219–231 (2012).Article 
CAS 
PubMed 

Google Scholar 
Simon, H. et al. Modeling transpiration and leaf temperature of urban trees: A case study evaluating the microclimate model ENVI-met against measurement data. Landsc. Urban Plan. 174, 33–40 (2018).Article 

Google Scholar 
Eamus, D., Boulain, N., Cleverly, J. & Breshears, D. D. Global change-type drought-induced tree mortality: Vapor pressure deficit is more important than temperature per se in causing decline in tree health. Ecol. Evol. 3, 2711–2729 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Eichenberg, D. et al. The effect of microclimate on wood decay is indirectly altered by tree species diversity in a litterbag study. J. Plant Ecol. 10, 170–178 (2017).Article 

Google Scholar 
Brockerhoff, E. G. et al. Forest biodiversity, ecosystem functioning and the provision of ecosystem services. Biodivers. Conserv. 26, 3005–3035 (2017).Article 

Google Scholar 
Martínez Pastur, G., Perera, A. H., Peterson, U. & Iverson, L. R. In Ecosystem Services from Forest Landscapes (eds Perera, A. H. et al.) 1–10 (Springer International Publishing, 2018).
Google Scholar 
Smithers, R. J. et al. Comparing the relative abilities of tree species to cool the urban environment. Urban Ecosyst. 21, 851–862 (2018).Article 

Google Scholar 
Shashua-Bar, L., Tsiros, I. X. & Hoffman, M. Passive cooling design options to ameliorate thermal comfort in urban streets of a Mediterranean climate (Athens) under hot summer conditions. Build. Environ. 57, 110–119 (2012).Article 

Google Scholar 
Song, J. & Wang, Z.-H. Impacts of mesic and xeric urban vegetation on outdoor thermal comfort and microclimate in Phoenix, AZ. Build. Environ. 94, 558–568 (2015).Article 

Google Scholar 
Meehl, G. A. & Tebaldi, C. More intense, more frequent, and longer lasting heat waves in the 21st century. Science 305, 994–997 (2004).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Christidis, N., Jones, G. S. & Stott, P. A. Dramatically increasing chance of extremely hot summers since the 2003 European heatwave. Nat. Clim. Change 5, 46–50 (2015).Article 
ADS 

Google Scholar 
Pfleiderer, P., Schleussner, C.-F., Kornhuber, K. & Coumou, D. Summer weather becomes more persistent in a 2 °C world. Nat. Clim. Change 9, 666–671 (2019).Article 
ADS 

Google Scholar 
Selten, F. M., Bintanja, R., Vautard, R. & van den Hurk, B. J. J. M. Future continental summer warming constrained by the present-day seasonal cycle of surface hydrology. Sci. Rep. 10, 4721. https://doi.org/10.1038/s41598-020-61721-9 (2020).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gartner, K., Nadezhdina, N., Englisch, M., Čermak, J. & Leitgeb, E. Sap flow of birch and Norway spruce during the European heat and drought in summer 2003. For. Ecol. Manag. 258, 590–599 (2009).Article 

Google Scholar 
Speak, A., Montagnani, L., Wellstein, C. & Zerbe, S. The influence of tree traits on urban ground surface shade cooling. Landsc. Urban Plan. 197, 103748. https://doi.org/10.1016/j.landurbplan.2020.103748 (2020).Article 

Google Scholar 
Rahman, M. A., Armson, D. & Ennos, A. R. A comparison of the growth and cooling effectiveness of five commonly planted urban tree species. Urban Ecosyst. 18, 371–389 (2015).Article 

Google Scholar 
Bowden, J. D. & Bauerle, W. L. Measuring and modeling the variation in species-specific transpiration in temperate deciduous hardwoods. Tree Physiol. 28, 1675–1683 (2008).Article 
PubMed 

Google Scholar 
Panferov, O. et al. The role of canopy structure in the spectral variation of transmission and absorption of solar radiation in vegetation canopies. IEEE Trans. Geosci. Remote Sens. 39, 241–253 (2001).Article 
ADS 

Google Scholar 
Lin, H., Chen, Y., Zhang, H., Fu, P. & Fan, Z. Stronger cooling effects of transpiration and leaf physical traits of plants from a hot dry habitat than from a hot wet habitat. Funct. Ecol. 31, 2202–2211 (2017).Article 

Google Scholar 
Fauset, S. et al. Differences in leaf thermoregulation and water use strategies between three co-occurring Atlantic forest tree species. Plant Cell Environ. 41, 1618–1631 (2018).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Chen, L., Zhang, Z. & Ewers, B. E. Urban tree species show the same hydraulic response to vapor pressure deficit across varying tree size and environmental conditions. PloS ONE 7, e47882. https://doi.org/10.1371/journal.pone.0047882 (2012).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gallego, H. A., Rico, M., Moreno, G. & Santa Regina, I. Leaf water potential and stomatal conductance in Quercus pyrenaica Willd. forests: Vertical gradients and response to environmental factors. Tree Physiol. 14, 1039–1047 (1994).Article 
PubMed 

Google Scholar 
Hölscher, D., Koch, O., Korn, S. & Leuschner, C. Sap flux of five co-occurring tree species in a temperate broad-leaved forest during seasonal soil drought. Trees 19, 628–637 (2005).Article 

Google Scholar 
Li, S. et al. Leaf gas exchange performance and the lethal water potential of five European species during drought. Tree Physiol. 36, 179–192 (2016).CAS 
PubMed 

Google Scholar 
Schnabel, F. et al. Cumulative growth and stress responses to the 2018–2019 drought in a European floodplain forest. Glob. Change Biol. 28, 1870–1883 (2022).Article 
CAS 

Google Scholar 
Sastry, A., Guha, A. & Barua, D. Leaf thermotolerance in dry tropical forest tree species: Relationships with leaf traits and effects of drought. AoB Plants 10, plx070. https://doi.org/10.1093/aobpla/plx070 (2018).Article 
CAS 
PubMed 

Google Scholar 
Banerjee, T. & Linn, R. Effect of vertical canopy architecture on transpiration, thermoregulation and carbon assimilation. Forests 9, 198. https://doi.org/10.3390/f9040198 (2018).Article 

Google Scholar 
Leuzinger, S. & Körner, C. Tree species diversity affects canopy leaf temperatures in a mature temperate forest. Agric. For. Meteorol. 146, 29–37 (2007).Article 
ADS 

Google Scholar 
Yi, K. et al. High heterogeneity in canopy temperature among co-occurring tree species in a temperate forest. J. Geophys. Res. Biogeosci. 125, e05892. https://doi.org/10.1029/2020JG005892 (2020).Article 

Google Scholar 
Hagemeier, M. & Leuschner, C. Functional crown architecture of five temperate broadleaf tree species: Vertical gradients in leaf morphology, leaf angle, and leaf area density. Forests 10, 265. https://doi.org/10.3390/f10030265 (2019).Article 

Google Scholar 
Raabe, K., Pisek, J., Sonnentag, O. & Annuk, K. Variations of leaf inclination angle distribution with height over the growing season and light exposure for eight broadleaf tree species. Agric. For. Meteor. 214–215, 2–11 (2015).Article 

Google Scholar 
Kafuti, C. et al. Foliar and wood traits covary along a vertical gradient within the crown of long-lived light-demanding species of the Congo Basin semi-deciduous forest. Forests 11, 35. https://doi.org/10.3390/f11010035 (2020).Article 

Google Scholar 
Peiffer, M., Bréda, N., Badeau, V. & Granier, A. Disturbances in European beech water relation during an extreme drought. Ann. For. Sci. 71, 821–829 (2014).Article 

Google Scholar 
Stratópoulos, L. M. F. et al. Tree species from two contrasting habitats for use in harsh urban environments respond differently to extreme drought. Int. J. Biometeorol. 63, 197–208 (2019).Article 
ADS 
PubMed 

Google Scholar 
McGloin, R. et al. Available energy partitioning during drought at two Norway spruce forests and a European beech forest in Central Europe. J. Geophys. Res. Atmos. 124, 3726–3742 (2019).Article 
ADS 

Google Scholar 
Schwaab, J. et al. Increasing the broad-leaved tree fraction in European forests mitigates hot temperature extremes. Sci. Rep. 10, 14153. https://doi.org/10.1038/s41598-020-71055-1 (2020).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Hari, V., Rakovec, O., Markonis, Y., Hanel, M. & Kumar, R. Increased future occurrences of the exceptional 2018–2019 Central European drought under global warming. Sci. Rep. 10, 12207. https://doi.org/10.1038/s41598-020-68872-9 (2020).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lyon, T. L., Weil, R. R. & Brady, N. C. The Nature and Properties of Soils 15th edn. (Pearson, 2017).
Google Scholar 
Zweifel, R., Böhm, J. P. & Häsler, R. Midday stomatal closure in Norway spruce—reactions in the upper and lower crown. Tree Physiol. 22, 1125–1136 (2002).Article 
CAS 
PubMed 

Google Scholar 
Rahman, M. A., Moser, A., Gold, A., Rötzer, T. & Pauleit, S. Vertical air temperature gradients under the shade of two contrasting urban tree species during different types of summer days. Sci. Total Environ. 633, 100–111 (2018).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Schuldt, B. et al. A first assessment of the impact of the extreme 2018 summer drought on Central European forests. Basic Appl. Ecol. 45, 86–103 (2020).Article 

Google Scholar 
Hochberg, U., Rockwell, F. E., Holbrook, N. M. & Cochard, H. Iso/anisohydry: A plant-environment interaction rather than a simple hydraulic trait. Trends Plant Sci. 23, 112–120 (2018).Article 
CAS 
PubMed 

Google Scholar 
Leuschner, C., Wedde, P. & Lübbe, T. The relation between pressure–volume curve traits and stomatal regulation of water potential in five temperate broadleaf tree species. Ann. For. Sci. 76, 60. https://doi.org/10.1007/s13595-019-0838-7 (2019).Article 

Google Scholar 
Bartlett, M. K., Scoffoni, C. & Sack, L. The determinants of leaf turgor loss point and prediction of drought tolerance of species and biomes: A global meta-analysis. Ecol. Lett. 15, 393–405 (2012).Article 
PubMed 

Google Scholar 
Hartmann, H., Link, R. M. & Schuldt, B. A whole-plant perspective of isohydry: Stem-level support for leaf-level plant water regulation. Tree Physiol. 41, 901–905 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Alonso-Forn, D. et al. Revisiting the functional basis of sclerophylly within the leaf economics spectrum of oaks: Different roads to Rome. Curr. For. Rep. 6, 260–281 (2020).
Google Scholar 
Hirons, A. D. & Thomas, P. A. Applied Tree Biology (John Wiley & Sons Ltd, 2017).Book 

Google Scholar 
Richter, R., Reu, B., Wirth, C., Doktor, D. & Vohland, M. The use of airborne hyperspectral data for tree species classification in a species-rich Central European forest area. Int. J. Appl. Earth Obs. Geoinf. 52, 464–474 (2016).ADS 

Google Scholar 
Qiu, G. et al. Effects of evapotranspiration on mitigation of urban temperature by vegetation and urban agriculture. J. Integr. Agric. 12, 1307–1315 (2013).Article 

Google Scholar 
Meier, F. & Scherer, D. Spatial and temporal variability of urban tree canopy temperature during summer 2010 in Berlin, Germany. Theor. Appl. Climatol. 110, 373–384 (2012).Article 
ADS 

Google Scholar 
Landsberg, J. J. & James, G. B. Wind profiles in plant canopies: Studies on an analytical model. J. Appl. Ecol. 8, 729–741 (1971).Article 

Google Scholar 
Gromke, C. & Ruck, B. Aerodynamic modelling of trees for small-scale wind tunnel studies. Forestry 81, 243–258 (2008).Article 

Google Scholar 
Baldocchi, D. D. Turbulent transfer in a deciduous forest. Tree Physiol. 5, 357–377 (1989).Article 
CAS 
PubMed 

Google Scholar 
Derby, R. W. & Gates, D. M. The temperature of tree trunks—Calculated and observed. Am. J. Bot. 53, 580–587 (1966).
Google Scholar 
Jayalakshmy, M. S. & Philip, J. Thermophysical properties of plant leaves and their influence on the environment temperature. Int. J. Thermophys. 31, 2295–2304 (2010).Article 
ADS 
CAS 

Google Scholar 
Pieruschka, R., Huber, G. & Berry, J. A. Control of transpiration by radiation. Proc. Natl. Acad. Sci. U.S.A. 107, 13372–13377 (2010).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Meili, N. et al. Tree effects on urban microclimate: Diurnal, seasonal, and climatic temperature differences explained by separating radiation, evapotranspiration, and roughness effects. Urban For. Urban Green. 58, 126970. https://doi.org/10.1016/j.ufug.2020.126970 (2021).Article 

Google Scholar 
Oogathoo, S., Houle, D., Duchesne, L. & Kneeshaw, D. Vapour pressure deficit and solar radiation are the major drivers of transpiration of balsam fir and black spruce tree species in humid boreal regions, even during a short-term drought. Agric. For. Meteorol. 291, 108063. https://doi.org/10.1016/j.agrformet.2020.108063 (2020).Article 
ADS 

Google Scholar 
Betts, M. G., Phalan, B., Frey, S. J. K., Rousseau, J. S. & Yang, Z. Old-growth forests buffer climate-sensitive bird populations from warming. Divers. Distrib. 24, 439–447 (2018).Article 

Google Scholar 
Pureswaran, D. S., Roques, A. & Battisti, A. Forest insects and climate change. Curr. For. Rep. 4, 35–50 (2018).
Google Scholar 
de Frenne, P. et al. Forest microclimates and climate change: Importance, drivers and future research agenda. Glob. Change Biol. 27, 2279–2297 (2021).Article 
ADS 

Google Scholar 
Woods, C. L., Cardelús, C. L. & DeWalt, S. J. Microhabitat associations of vascular epiphytes in a wet tropical forest canopy. J. Ecol. 103, 421–430 (2015).Article 

Google Scholar 
Nakamura, A. et al. Forests and their canopies: Achievements and horizons in canopy science. Trends Ecol. Evol. 32, 438–451 (2017).Article 
PubMed 

Google Scholar 
European State of the Climate 2020, Copernicus Climate Change Service, Full report: climate.copernicus.eu/ESOTC/2020Munzi, S. et al. Lichens as ecological indicators in urban areas: beyond the effects of pollutants. J. Appl. Ecol. 51, 1750–1757 (2014).Article 

Google Scholar 
Kaspari, M., Clay, N. A., Lucas, J., Yanoviak, S. P. & Kay, A. Thermal adaptation generates a diversity of thermal limits in a rainforest ant community. Glob. Change Biol. 21, 1092–1102 (2015).Article 
ADS 

Google Scholar 
Baudier, K. M., Mudd, A. E., Erickson, S. C. & O’Donnell, S. Microhabitat and body size effects on heat tolerance: implications for responses to climate change (army ants: Formicidae, Ecitoninae). J. Anim. Ecol. 84, 1322–1330 (2015).Article 
PubMed 

Google Scholar 
Merinero, S., Dahlberg, C. J., Ehrlén, J. & Hylander, K. Intraspecific variation influences performance of moss transplants along microclimate gradients. Ecology 101, e02999. https://doi.org/10.1002/ecy.2999 (2020).Article 
PubMed 

Google Scholar 
Ben-Yakir, D. & Fereres, A. The effects of UV radiation on arthropods: A review of recent publications (2010–2015). Acta Hortic. 1134, 335–342 (2016).Vanhaelewyn, L., van der Straeten, D., de Coninck, B. & Vandenbussche, F. Ultraviolet radiation from a plant perspective: The plant-microorganism context. Front. Plant Sci. 11, 597642. https://doi.org/10.3389/fpls.2020.597642 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Jansen, E. Das Naturschutzgebiet Burgaue; Staatliches Umweltfachamt: Leipzig, Germany (1999).Sächsisches Landesamt für Umwelt, Landwirtschaft und Geologie (LFULG) & DWD Deutscher Wetterdienst (2019) [ed.]: 2018 Wetter trifft auf Klima. Dresden, Leipzig. https://www.klima.sachsen.de/download/Jahresrueckblick2018_A5_OeA.pdf.Haase, D. & Gläser, J. Determinants of floodplain forest development illustrated by the example of the floodplain forest in the District of Leipzig. For. Ecol. Manag. 258, 887–894 (2009).Article 

Google Scholar 
Patzak, R., Richter, R., Engelmann, R. A. & Wirth, C. Tree crowns as meeting points of diversity generating mechanisms: A test with epiphytic lichens in a temperate forest. Preprint at: https://www.biorxiv.org/content/https://doi.org/10.1101/2020.01.03.894303v1.full (2020).Meinen, C., Leuschner, C., Ryan, N. T. & Hertel, D. No evidence of spatial root system segregation and elevated fine root biomass in multi-species temperate broad-leaved forests. Trees 23, 941–950 (2009).Article 

Google Scholar 
van der Zande, D., Stuckens, J., Verstraeten, W. W., Muys, B. & Coppin, P. Assessment of light environment variability in broadleaved forest canopies using terrestrial laser scanning. Remote Sens. 2, 1564–1574. https://doi.org/10.3390/rs2061564 (2010).Article 
ADS 

Google Scholar 
Köstner, B., Granier, A. & Cermák, J. Sapflow measurements in forest stands: Methods and uncertainties. Ann. For. Sci. 55, 13–27 (1998).Article 

Google Scholar 
Granier, A. Evaluation of transpiration in a Douglas-fir stand by means of sap flow measurements. Tree Physiol. 3, 309–320 (1987).Article 
CAS 
PubMed 

Google Scholar 
Metzger, J. M. & Oren, R. The effect of crown dimension on transparency and the assessment of tree health. Ecol. Appl. 11, 1634–1640 (2001).Article 

Google Scholar 
Lefcheck, J. S. piecewiseSEM: Piecewise structural equation modelling in r for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579 (2016).Article 

Google Scholar 
Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & R Core Team (2020). nlme: Linear and nonlinear mixed effects models. R package version 3.1-151, https://CRAN.R-project.org/package=nlme.Dornelas, M. et al. Quantifying temporal change in biodiversity: Challenges and opportunities. Proc. Biol. Sci. 280, 20121931. https://doi.org/10.1098/rspb.2012.1931 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Shipley, B. The AIC model selection method applied to path analytic models compared using a d-separation test. Ecology 94, 560–564 (2013).Article 
PubMed 

Google Scholar 
R Core Team (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/. More

Bacterial strains and mediaA collection of 185 genome-sequenced bacterial isolates, described previously14, was utilized to assemble the synthetic community used in this work. These isolates were obtained from surface-sterilized Brassicaceae roots, primarily Arabidopsis thaliana, grown in two soils from North Carolina, USA35. This isolate collection includes strains V. paradoxus CL14, Arthrobacter CL28, Acinetobacter CL69 and Acinetobacter CL71, which are also used in this work in individual strain contexts. V. paradoxus CL14 ΔHS33, which has a clean deletion of genes with gene ID 2643613677 through 2643613653 was constructed previously14 and used here. Additional strains were obtained from the American Type Culture Collection (ATCC): E. soli LF7 (ATCC BAA-2102), R. pomeroyi (ATCC 700808) and B. japonicum (ATCC 10324). P. phytofirmans PsJN (DSMZ 17436) was obtained from the DSMZ-German Collection of Microorganisms and Cell Cultures. P. putida strain 1290 was generously provided by Johan Leveau (University of California Davis). Pseudomonas strain Root 562 was generously provided by Paul Schulze-Lefert (Max-Planck-Gesellschaft). All bacteria, with exceptions noted below, were routinely grown on LB agar plates (10 g l−1 tryptone, 5 g l−1 yeast extract, 10 g l−1 NaCl, 1.5% (w/v) agar) and in 2xYT liquid medium (16 g l−1 tryptone, 10 g l−1 yeast extract, 5 g l−1 NaCl) at 28 °C. The 175-member (185-member minus 10 Variovorax strains) synthetic community (SC185-10V) was grown on KB medium as was done previously to culture this synthetic community14. B. japonicum (ATCC 10324) was routinely grown on liquid and solidified YM medium (1 g l−1 yeast extract, 10 g l−1 mannitol, 0.5 g l−1 dipotassium phosphate, 0.2 g l−1 magnesium sulfate, 0.1 g l−1 NaCl, 1 g l−1 CaCO4, pH 6.8, solidified with 1.5% agar as necessary) at 28 °C. R. pomeroyi (ATCC 700808) was routinely grown on liquid and solidified LB medium supplemented with 2% sea salt (Millipore Sigma S9883) and solidified with 1.5% (w/v) agar as necessary. M9 base medium was formulated using 1x M9 minimal salts medium (Sigma M6030) supplemented with 2 mM MgSO4, 0.1 mM CaCl2 and 10 µM FeSO4. A carbon source or sources were added to this M9 base medium to support bacterial growth. Unique strains constructed in this study are available upon request.Bacterial 16S rRNA sequencingBacterial colonization of Arabidopsis roots was assessed using a method similar to the previous study14. Roots from 8–10 plants were collected into sterilized 2 ml tubes containing three 4 mm glass beads and root fresh weight in each tube was obtained. Five such samples were collected for each bacterial treatment. The roots were washed three times with sterile distilled water and stored at −80 °C until further processing. The roots were then lyophilised for 48 h using a Labconco freeze-dry system and pulverized using an MPBio tissue homogenizer. DNA was extracted from the root samples and bacterial cell pellets saved from the bacteria for input into the experiment using the DNeasy PowerSoil HTP 96 kit (Qiagen) according to manufacturer instructions. The V3-V4 region of the bacterial 16S rRNA gene was amplified and sequenced as previously described14.16S amplicon sequence data processingThe 16S sequencing data from synthetic community experiments were processed as previously described14. Briefly, usable read output from MT-Toolbox36 (reads with 100% primer sequences that successfully merged with their pair) were filtered for quality with Sickle37 by not allowing any window with Q score under 20. The resulting sequences were globally aligned to a 16S rRNA gene sequence reference dataset from genome assemblies of the synthetic community members. For strains that do not have an intact 16S rRNA sequence in their assembly, Sanger sequencing was used to obtain the 16S rRNA gene sequence of the strains for inclusion in the reference dataset. The reference dataset also included sequences from Arabidopsis organellar sequences and known bacterial contaminants. Sequence alignment was performed with USEARCH v.7.109038 using the optional usearch_global at a 98% identity threshold. On average, 85% of read sequences matched an expected isolate. The 185 isolates of our 185-member synthetic community could not all be distinguished from one another on the basis of the V3-V4 sequence. They were thus classified into 97 unique sequences encompassing a set of identical (clustered at 100%) V3-V4 sequences coming from a single or multiple isolate strains. An isolate abundance table was created from the sequence mapping results.We estimated 16S rRNA absolute abundance using a plasmid spike-in method39. Synthetic DNA was spiked at known quantities into samples before DNA extraction and the ratio of added to recovered synthetic DNA served as a conversion factor by which the total number of 16S rRNA molecules in a given sample was estimated. We designed a plasmid which included 16S V3-V4 primer binding sequences flanking a randomly generated DNA sequence matching the most frequent length and Guanine + Cytosine (GC) content of amplicons generated using the same primer sequences from wild soil. These sequences were synthesized by Geneart (Invitrogen) and supplied cloned in plasmid pMA-T. The plasmid was transformed into E. coli and isolated using a miniprep spin kit (Qiagen). Specific volumes of this isolated plasmid were then added to individual samples before DNA extraction to spike-in approximately 20% of the predicted 16S copies occurring within the sample. We performed colony-forming units (c.f.u.) counting using similarly treated plant samples (that is, growth on SynCom-inoculated agar plates) to obtain an estimate of the 16S copy number per mg fresh weight of plant roots. We plated serial dilutions of plant root samples ground in MgCl2 on LB to perform c.f.u. counts. The c.f.u. count multiplied by a given sample’s fresh weight were used to calculate sample-specific predicted 16S copy numbers.Plant growth conditions and root growth inhibition assayA. thaliana ecotype Col-0 seeds were sterilized in 70% household bleach, 0.2% Tween-20 for 10 min with vigorous agitation and then rinsed 10 times with sterile distilled water. Seeds were stratified at 4 °C in sterile distilled water for 1–2 d. Plants were germinated for 7 d on 0.5x MS agar medium (2.22 g l−1 PhytoTech Labs M-404: Murashige & Skoog modified basal medium with Gamborg vitamins, 0.5 g l−1 MES hydrate, pH adjusted to 5.7, solidified with 1% (w/v) agar) supplemented with 0.5% (w/v) sucrose in vertical 12 ×12 cm square plates under long-day conditions (21 °C/18 °C, 16 h light/8 h dark, day/night cycle). Then 8 to 10 plants were aseptically transferred to 12 ×12 cm plates containing 0.5x MS agar medium without sucrose where the medium surface was spread with the bacterial inoculum. For assays with IAA addition, 100 nM IAA was added to the medium before pouring the plates. The plant root tip location was marked on plates after transfer to record the initial root tip position. The plates containing the plants and bacteria were incubated vertically under short-day conditions (22 °C/18 °C, 9 h light/15 h dark, day/night cycle) for an additional 11 d. Plates were imaged on a document scanner and primary root elongation was determined using imageJ to quantify the change in root tip position from the initial to the final position.Bacterial inoculation of plantsIndividual bacterial strains were grown on agar plates of the media types specified above at 28 °C. Before plant inoculation, a single colony was picked into the appropriate liquid medium and grown at 28 °C to late exponential or early stationary phase. To remove the medium from the bacteria before inoculation, strains were washed three times in sterile 10 mM MgCl2. The optical density at 600 nm (OD600) was measured for each washed strain and normalized to OD600 of 0.01 in 10 mM MgCl2. For plant experiments with mono-association of an individual strain, 100 µl of OD600 = 0.01 washed bacteria was spread on the 12 ×12 cm plate before plant transfer. For experiments in duo-association with Arthrobacter CL28, 100 µl of OD600 = 0.01 washed Arthrobacter CL28 was spread along with 100 µl of OD600 = 0.01 of the second strain.The 175-member synthetic community (SC185-10V) was prepared as described for the 185-member synthetic community used previously14 by leaving out the 10 isolates from the genus Variovorax. Briefly, 7 d before plant transfer, strains were inoculated individually into 600 µl KB medium in a 96-well plate and grown at 28 °C for 5 d. At 2 d before plant transfer, 20 ul from these 5-day-old cultures were transferred to 380 ul fresh KB medium in a new set of 96-well plates and both sets of plates were returned to the incubator for 2 d. This resulted in two cultures of each strain, one 7 d old and the other 2 d old, which were combined. The OD600 of the strains in each well was measured and the strains were combined while normalizing the OD600 of each strain in the pool. This pool was washed twice with 10 mM MgCl2 and diluted to OD600 = 0.2. For experiments with the SC185-10V SynCom, 100 µl of this OD600 = 0.2 washed pool was spread on 12 ×12 cm plates. For treatments where an additional strain was added to the SC185-10V SynCom, the individual strain was washed as described above, diluted to OD600 = 0.0034 in 10 mM MgCl2, and 100 µl of this dilution was spread on the plates with the SC185-10V SynCom. This addition of the individual strain corresponded to an OD600 three times that of a single strain in the SC185-10V SynCom (0.0034 = (0.2/175) × 3). For the addition of the 10 Variovorax strains to the SC185-10V SynCom experiment, the 10 Variovorax strains were grown individually in 2xYT medium from colonies grown on plates. The OD600 values of the 10 cultures were measured and the 10 strains were pooled while normalizing the OD600 of each strain to the same value. This mixture of the 10 Variovorax strains was then treated as the individual strains for washing and addition of 100 µl of OD600 = 0.0034 to the SC185-10V SynCom on plates.Construction of vectors with Variovorax CL14 iad gene insertsPortions of the V. paradoxus CL14 IAA degradation locus were subcloned into broad host range vector pBBR1MCS-232. Primers JMC579 through JMC604 (Supplementary Table 8) were used to amplify 3–5 kb segments of the locus by PCR using Q5 DNA polymerase (New England Biolabs). These primers were designed to amplify sections beginning and ending at gene start codons and with appropriate overlapping sequences for Gibson assembly either into the pBBR1MCS-2 backbone or to the adjacent section to make larger vector inserts, as appropriate. The pBBR1MCS-2 vector backbone was prepared for Gibson assembly by amplifying the vector by PCR using primers JMC577 and JMC578 (Supplementary Table 8) and subsequently treating with DpnI to remove circular vector template. PCR fragments were cleaned up as necessary using the QIAquick PCR purification kit (Qiagen). Appropriate fragments were mixed to construct the vectors by Gibson assembly using HiFi DNA Assembly Mastermix (New England Biolabs) according to manufacturer instructions. Gibson assembly products were transformed into NEB 10beta chemically competent E. coli (New England Biolabs) and selected on LB plates supplemented with 50 µg ml−1 kanamycin. Vectors were miniprepped using either the ZR plasmid miniprep classic kit or Zymo BAC DNA miniprep kit (Zymo Research) and confirmed via restriction mapping with PstI-HF (New England Biolabs) and Sanger sequencing (Genewiz).To construct vectors that are derivatives of pBBR1::70–66, the Q5 site-directed mutagenesis kit (New England Biolabs) was used for gene deletion. Briefly, vector pBBR1::70–66 was used as a PCR template and portions of this vector were amplified by PCR using primers JMC641 through JMC650 (Supplementary Table 8) and Q5 DNA polymerase (New England Biolabs). PCR products were cleaned up and circularized using KLD Mastermix (New England Biolabs). The product was transformed into NEB 10beta chemically competent E. coli (New England Biolabs) and selected on LB plates supplemented with 50 µg ml−1 kanamycin. Vectors were miniprepped and Sanger sequenced as described above to confirm the construction of the correct vectors.Conjugation of vectors to V. paradoxus CL14 ΔHS33Vectors were conjugated into V. paradoxus CL14 ΔHS33 using tri-parental mating. The helper E. coli strain carrying plasmid pRK201340 and donor NEB 10beta E. coli strains containing the pBBR1MCS-2-based vectors with Variovorax IAA degradation locus gene inserts were cultured in LB media containing 50 µg ml−1 kanamycin at 37 °C. V. paradoxus CL14 ΔHS33 was grown in 2xYT medium containing 100 µg ml−1 ampicillin at 28 °C. V. paradoxus CL14 wild type and derivative strains such as ΔHS33 are naturally resistant to ampicillin and this ampicillin selection allows for recovery of only Variovorax from the conjugation reaction. To prepare for conjugation, all bacteria were pelleted by centrifugation at 5,000 × g for 5 min and washed 3 times in 2xYT medium without antibiotics. For each conjugation reaction, equal volumes (100–300 µl) of each of the three washed bacteria: recipient V. paradoxus CL14 ΔHS33, donor NEB 10beta E. coli containing a pBBR1MCS-2-based vector, and helper E. coli pRK2013 were mixed. Control conjugation mixtures of each pair of strains and individual strains alone were performed in parallel to ensure successful selection of exconjugants only from mixtures of all three strains together. Conjugation mixtures were pelleted by centrifugation at 5,000 × g for 5 min, resuspended in 50 µl 2xYT media, transferred to LB media plates without antibiotics and allowed to dry in a laminar flow hood. These conjugation plates were incubated overnight at 28 °C. After 18–24 h, exconjugants were selected by streaking from the pooled conjugation mixtures on the LB plate without antibiotics to LB plates containing 50 µg ml−1 kanamycin and 100 µg ml−1 ampicillin. This selects for only V. paradoxus CL14 ΔHS33 (ampicillin resistant) containing the pBBR1MCS-2-based vector (kanamycin resistant). Individual colonies were picked into and subsequently cultured in 2xYT medium containing 50 µg ml−1 kanamycin and 100 µg ml−1 ampicillin at 28 °C.Construction of V. paradoxus CL14 gene deletionsUnmarked gene deletions in V. paradoxus CL14 were constructed as described previously14 using the suicide vector backbone pMo130 originally developed for gene knockouts in Burkholderia spp.41. Primers JMC203 and JMC204 (Supplementary Table 8) were used to amplify the pMO130 vector backbone by PCR. This product was subsequently treated with DpnI (New England Biolabs) to digest circular template DNA. Primers JMC605 through JMC612 and JMC671 through JMC677 (Supplementary Table 8) were used to amplify flanking regions for the gene deletion targets from V. paradoxus CL14 genomic DNA. All PCR was performed using Q5 DNA polymerase (New England Biolabs) and products were cleaned up, as appropriate, with the QIAquick PCR purification kit (Qiagen). These PCR products were assembled into suicide vectors using HiFi Gibson Assembly Mastermix (New England Biolabs), transformed into chemically competent NEB 5alpha E. coli (New England Biolabs), and selected on LB plates with 50 µg ml−1 kanamycin. Vectors were miniprepped using the ZR plasmid miniprep classic kit (Zymo Research) and confirmed by Sanger sequencing (Genewiz). Confirmed vectors were transformed into the chemically competent bi-parental mating strain E. coli WM3064. Transformants were selected at 37 °C on LB media supplemented with 50 µg ml−1 kanamycin and 0.3 mM diaminopimelic acid (DAP), and single colonies picked into LB medium also with 50 µg ml−1 kanamycin and 0.3 mM DAP.Bi-parental mating was performed by growing E. coli WM3064 containing the appropriate suicide vector as described above, and V. paradoxus CL14 was grown in 2xYT medium containing 100 µg ml−1 ampicillin at 28 °C. Both E. coli and Variovorax were washed separately three times using 2xYT medium, then mixed in a 1:1 ratio and pelleted. All centrifugation steps were performed at 5,000 × g for 5 min. The pelleted conjugation mixtures were resuspended in 1/10 the volume of 2xYT, plated on LB agar with 0.3 mM DAP and grown at 28 °C overnight. Exconjugants from these plates were streaked out and grown on LB agar with 100 µg ml−1 ampicillin, 50 µg ml−1 kanamycin, and no DAP at 28 °C. These strains were purified by streaking and growing on plates of the same medium once more. These strains with suicide vector integration were then grown once in liquid LB containing 100 µg ml−1 ampicillin and 1 mM isopropyl 1-thio-d-galactopyranoside (IPTG) at 28 °C and then streaked on plates containing media with 10 g l−1 tryptone, 5 g l−1 yeast extract, 100 g l−1 sucrose, 1.5% agar, 100 µg ml−1 ampicillin and 1 mM IPTG. Colonies from these plates were picked and grown in the same liquid media. These strains were then assessed for gene deletion by PCR using primers JMC657 through JMC660 and JMC697 through JMC699 (Supplementary Table 8). The Quick-DNA miniprep kit (Zymo Research) was used to isolate all genomic DNA for PCR screening. To purify the knockout strains, they were streaked and grown out three times on LB plates containing 100 µg ml−1 ampicillin before a final PCR confirmation. To check the purity of the final strains, PCR was performed with one primer outside the deletion region and one inside the deleted gene to ensure no product is produced for the knockout strain. The sequences for the primers used for this PCR reaction (JMC691, JMC717, JMC718, JMC693 and JMC694) can be found in Supplementary Table 8.Measurement of bacterial growth and IAA degradationIndividual strains were grown in 5 ml cultures in various media types supplemented with IAA at 28 °C and 250 r.p.m. To screen the V. paradoxus CL14 ΔHS33 pBBR1 vector complemented mutants, 2xYT medium supplemented with 0.1 mg ml−1 IAA was used. For comparison of other V. paradoxus CL14 mutants, M9 medium with 15 mM succinate and 0.1 mg ml−1 (0.57 mM) IAA was used. For comparison of IAA-degrading strains from diverse genera, M9 medium with 0.1% (w/v) casamino acids (Bacto) and 0.1 mg ml−1 IAA was used. For R. pomeroyi, 2% (w/v) sea salts were added to this M9 medium with casamino acids and IAA. The pBBR1 vector library in E. coli NEB 10beta was screened in LB medium supplemented with 0.04 mg ml−1 IAA and grown at 37 °C and 250 r.p.m. For all media types, IAA was solubilized in 100% ethanol at 20 mg ml−1 and diluted to 0.1 mg ml−1 in the media, resulting in 0.5% (v/v) ethanol in the media.To measure growth, a 200 µl sample was taken from the growing cultures and OD600 was determined on an Infinite M200 Pro plate reader (Tecan). Subsequently, cells were pelleted by centrifugation at 4,200 × g for 15 min and 50 µl of supernatant was transferred to a new 96-well plate and frozen at −80 °C until further analysis. IAA degradation was determined by thawing the plates containing 50 µl aliquots of culture supernatant and combining this with 100 µl of Salkowski reagent (10 mM ferric chloride and 35% perchloric acid)42. This was performed alongside mixing 50 µl of IAA standards with 100 µl of Salkowski reagent in the same 96-well plate format. Colour development was allowed to proceed for 1 h and absorbance was read at 530 nm on the Infinite M200 Pro plate reader (Tecan). The absorbances measured were converted to IAA concentration on the basis of the absorbances measured for the IAA standards.Liquid Chromatography Dual Mass Spectroscopy (LC–MS/MS) metabolomics on Variovorax IAA degradationV. paradoxus CL14 was grown in 5 ml cultures of M9 minimal medium supplemented with either 0.1 mg ml−1 IAA, 0.1 mg ml−1 13C6-IAA (with the 6 carbons of the benzene ring of the indole labelled, Cambridge Isotope Laboratories CLM-1896-PK), and/or 15 mM succinate. Cultures and parallel media controls were incubated at 28 °C with shaking at 250 r.p.m. Cultures and media controls were centrifuged (4,200 × g for 15 min at 4 °C) to pellet cells; supernatants were transferred to new tubes and both pellets (intracellular fraction) and supernatants (extracellular fraction) were stored frozen at −80 °C until extraction. All subsequent work was performed over dry ice or in chilled cold blocks. Frozen pellets from the intracellular fraction were thawed for 3 h at 4 °C, then 800 µl of cold LCMS-grade water was added to the pellets with repeated pipetting to break up the pellet until visually homogeneous. Samples were then quickly returned to −80 °C to freeze the suspension. Frozen pellet suspensions and extracellular solutions were lyophilised until dry. The cells from the dried pellet suspensions were lysed and homogenized with a bead mill (BioSpec Mini-Beadbeater-96) using one sterile 3.2 mm steel ball in each tube for 3 rounds of 5 s each with 10 s breaks in between to reduce heat production. Dried extracellular samples were concentrated by resuspension in 100 µl LCMS- grade methanol, vortexed 3 times for 10 s each, water bath sonicated for 20 min, incubated at 4 °C overnight, centrifuged (1,000 × g, 4 °C, 5 min), and the methanol supernatant was dried using a speed vacuum concentrator. On the day of LC–MS/MS analysis, homogenized dry material was suspended in LCMS-grade methanol with internal standard mix (100 µM U-13C/15N-labelled amino acids, SIGMA 767964). Intracellular samples were suspended at 11.1 µl mg−1 of original sample cell pellet wet weight; extracellular samples were suspended at 38.9 µl mg−1 of corresponding cell pellet wet weight from the culture. The solutions were vortexed 3 times for 10 s each, bath sonicated in ice water for 10 min, chilled at −20 °C for 10 min, then centrifuged (10,000 × g, 5 min, 10 °C) to pellet insoluble material. Supernatants containing the methanol extracts were filtered through 0.22 µm PVDF microcentrifuge filtration tubes (10,000 × g, 5 min, 10 °C); filtrates were transferred to glass vials and immediately capped. Filtrates were then analysed by LC–MS/MS using an Agilent 1290 UHPLC system connected to a Thermo Q Exactive Hybrid Quadrupole-Orbitrap mass spectrometer equipped with a heated electrospray ionization (HESI-II) source probe. Extracts were chromatographically separated on a ZORBAX RRHD Eclipse Plus C18, 95 Å, 2.1 × 50 mm, 1.8 µm column (Agilent) for non-polar metabolomics. Separation, ionization, fragmentation and data acquisition parameters are specified in Supplementary Table 7. Briefly, metabolites were separated by gradient elution followed by MS1 and data-dependent (top 2 most abundant MS1 ions not previously fragmented in last 7 s) MS2 collection; targeted data analysis was performed by comparison of sample peaks to a library of analytical standards analysed under the same conditions or by searching the raw data files for predicted m/z values based on structural information of compounds of interest. Three parameters were compared: matching m/z, retention time and fragmentation spectra using Metabolite Atlas (https://github.com/biorack/metatlas)43,44. Identification and standard reference comparison details are provided in Supplementary Table 6. Raw and processed data are available for download at the JGI Joint Genome Portal under ID 1340427. Statistical comparisons were performed using R version 3.6.2, using package agricolae 1.3–5 and stats 3.6.245; boxplots were generated with base R graphics using the boxplot function.Phylogenomic analysisTo guide the delineation of the IAA degradation operons across the bacterial tree of life, we constructed two Hidden Markov Model (HMM) profiles of the genes iacC and iacD by subsetting all homologous genes from the previously validated operons (Extended Data Fig. 4). In parallel, we downloaded the assembly files for all available complete genomes deposited in the NCBI RefSeq 202 repository46. For the 220,000 assembly files downloaded, we performed open reading frame (ORF) prediction using prodigal. We then used the two HMM profiles described above to query the predicted ORFs. Utilizing ad hoc scripts, we constructed a table of HMM hits along the genomes scanned and subset genomic loci where both iacC and iacD genes appeared adjacent to one another. The logic of using the iacC and iacD genes as anchor genes for our search is that the adjacent physical location of both iacC and iacD homologues is a conserved feature across all previously experimentally validated IAA-degrading operons (Extended Data Fig. 4). Next, for each region containing the adjacent iacC and iacD homologue genes, we extracted the gene neighbourhood adjacent to the anchor hit by extracting the amino acid sequence of ORFs +10 kb and −10kb with respect to the anchor hit. Using hmmscan from the Hmmer v3.1.b2 suite47, we performed HMM profiling in all ORFs extracted via our neighbourhood delineation against the COG database version 2003. Finally, we used the COG profiles across the neighbourhoods to create a matrix describing the prevalence of COGs across the regions (candidate regions) with the adjacent iacC and iacD homologue genes.For each genome containing at least one candidate region, we performed taxonomic classification using the GTDB database48. Due to the size of our estimated genomic matrix and to reduce potential biases due to over-representation of certain lineages within RefSeq, we performed principal coordinate analysis (PCoA) using a reduced matrix containing one representative candidate region per species. Species labelling was obtained from the GTDB taxonomic classification described above. PCoA was performed using the oh.pco function from the ohchibi package49, taking as input a binary version of the gene matrix described above. We classified candidate reads into the two types of IAA-degrading operon (iac-like and iad-like), utilizing a majority count-based approach using marker COGs conserved between the previously experimentally validated IAA-degrading operons (Extended Data Fig. 4). Specifically, for each potential operon, we determined the prevalence of COGs that a priori (Extended Data Fig. 4) showed differential prevalence across the two degrading operons (for example, iacA, iacB and iacI are markers of the iac operon, while iorB/iadB and iotA/iadA are exclusive markers of the iad-like operon). Hybrid gene clusters were defined as operons that exhibited the hallmark COGs of both operons.In parallel, we performed phylogenetic inference over all the genomes belonging to genera with at least one representative strain harbouring any of the two types of IAA-degrading operon. This phylogenetic tree was constructed using a super-matrix-based approach as previously described35. Finally, for each genus with at least one assembly harbouring a positive IAA-degrading operon, we estimated the prevalence of the trait across the genus by dividing the total number of isolates with detectable IAA degradation locus by the total number of isolates belonging to that genus in the dataset. In addition, to see the phylogenetic evenness of the distribution of the IAA degradation trait across each genus, we calculated the phylogenetic ratio by calculating the ratio between the average phylogenetic distance (computed via the cophenetic.phylo function from the ape R package50) of isolates with a detectable IAA degradation locus and the total average phylogenetic distance of all isolates within that genus. We constructed the MarR phylogeny using the MarR sequences from candidate regions with 100% markers of one of the two types of IAA-degrading operon. Amino acid sequences of the MarR homologues were aligned using MAFFT51 and phylogenetic inference was performed using FastTree 252.RNA-seq on Variovorax strainsV. paradoxus CL14 was grown in 5 ml cultures of M9 minimal medium supplemented with 15 mM succinate and 0.5% (v/v) ethanol alone or containing IAA. IAA was at a final concentration of 0.1 mg ml−1 in the medium to which it was added. Cultures were prepared at a starting OD600 of 0.02 and incubated at 28 °C, shaking at 250 r.p.m. Cells from all samples were collected for RNA-seq at 18 h to ensure IAA was still present in the cultures of strains that degraded IAA most rapidly. Cells were pelleted by centrifuging the culture at 4,200 × g for 15 min and removing the supernatant. Cell pellets were frozen at −80 °C before RNA extraction. To extract RNA, cells were lysed in TRIzol reagent (Invitrogen) according to manufacturer instructions for lysis and phase separation. After these steps, RNA was purified from the aqueous phase using the RNeasy mini kit (Qiagen) including the optional on-column DNase digestion with RNase-free DNase set (Qiagen). Total RNA was quantified using the Qubit 2.0 fluorometer (Invitrogen) and RNA-seq libraries were prepared using the Universal Prokaryotic RNA-Seq Prokaryotic AnyDeplete kit (Tecan) according to manufacturer instructions. The resulting libraries were pooled and sequenced on the Illumina HiSeq4000 to generate 50 bp single-end reads.RNA-seq data analysisThe V. paradoxus CL14 RNA-seq sequence data were analysed as described previously14. Briefly, the raw reads were mapped to the V. paradoxus CL14 genome (fasta file available at https://github.com/isaisg/variovoraxRGI/blob/master/rawdata/2643221508.fna) using bowtie253 with the ‘very sensitive’ flag. Hits to each individual coding sequence were counted and annotated using the function featureCounts from the R package Rsubread54, inputting the V. paradoxus CL14 gff file (available at https://github.com/isaisg/variovoraxRGI/blob/master/rawdata/2643221508.gff) and using the default parameters with the flag allowMultiOverlap = FALSE. Finally, DESeq255 was used to estimate Differentially Expressed Genes (DEGs) between treatments, with the corresponding fold-change estimates and False Discovery Rate (FDR) adjusted P values. For visualization purposes, we performed z-score standardization of each gene across samples and we visualized this standardized expression values utilizing a heat map constructed using ggplot256. These data can be found in Supplementary Table 9.MarR protein expression and purificationThe coding sequence for each gene can be found in Supplementary Table 10. MarR expression plasmids were synthesized as codon-optimized genes for E. coli expression by BioBasic in the pLIC-His N-term vector (pMCSG7) and transformed into E. coli BL21 (DE3) Gold cells for expression. Cells were grown in the presence of ampicillin in LB medium with shaking at 225 r.p.m. at 37 °C to an OD600 of 0.5, at which point the temperature was reduced to 18 °C. At an OD600 of 0.8, protein expression was induced by the addition of 0.1 mM IPTG and incubation continued overnight. Cells were collected by centrifugation at 4,500 × g for 20 min at 4 °C in a Sorvall (model RC-3B) swinging bucket centrifuge. Cell pellets were resuspended in buffer A (20 mM potassium phosphate, pH 7.4, 50 mM imidazole, 500 mM NaCl), DNase, lysozyme and a Roche Complete EDTA-free protease inhibitor tablet. Resuspended cells were sonicated and clarified via centrifugation at 17,000 × g for 60 min in a Sorvall (model RC-5B) swinging bucket centrifuge. The lysate was applied to a nickel-nitrilotriacetic acid HP column (GE Healthcare) on an Aktaxpress Fast Performance Liquid Chromatography (FPLC) system (Amersham Bioscience) and washed with buffer A. Protein was eluted with buffer B (20 mM potassium phosphate, pH 7.4, 500 mM imidazole, 500 mM NaCl). Fractions containing the protein of interest were combined and passed over a HiLoadTM 16/60 SuperdexTM 200 gel filtration column. Proteins were eluted in S200 buffer (20 mM HEPES, pH 7.4, 300 mM NaCl). Fractions were combined and concentrated for long-term storage at −80 °C.MarR mutant proteins were created by site-directed mutagenesis using primers from Integrated DNA Technologies. The mutant plasmids were sequenced to confirm the mutations. The mutants were produced and purified using E. coli BL21 (DE3) Gold as described above.Ligand binding studies by isothermal titration calorimetry (ITC)All ITC measurements were performed at 25 °C using an Auto-ITC200 microcalorimeter (MicroCal/GE Healthcare). The buffer employed was 20 mM HEPES, pH 7.4, 50 mM NaCl and 0.5% dimethly sulfoxide (DMSO) for protein/ligand binding and 20 mM HEPES, pH 7.4 and 300 mM NaCl for DNA/protein binding experiments. For ligand binding experiments, the calorimetry cell (volume 200 ml) was loaded with MarR wild-type, mutant or homologue protein at a concentration of 50 μM. The syringe was loaded with a ligand concentration of 0.5 or 2 mM. For DNA binding experiments, wild-type MarR_73 did not bind any of the DNA oligos examined; however, we hypothesized that this arose from the ability of this native receptor to remain bound to ligands retained from its recombinant expression in E. coli. Thus, we employed the MarR_73 S28A protein with reduced ligand binding capacity. Here, the calorimetry cell was loaded with duplex oligo at a concentration of 25 μM and the syringe was loaded with MarR S28A mutant protein, which was necessary to prevent ligand binding during expression and purification, at a concentration of 0.5 mM. A typical injection protocol included a single 0.2 μl first injection followed by 20 1.5 μl injections of the syringe sample into the calorimetry cell. The spacing between injections was kept at 180 s and the reference power at 8 μcal s−1. The data were analysed using Origin for ITC version 7.0 software supplied by the manufacturer and fit well to a one-site binding model. Two independent ITC measurements were performed for each condition. A non-integer N value (for example, 0.73 in Fig. 2a) indicates that some protein monomers may not be in an active conformation, and thus do not bind ligand. Additionally, small measurement errors in assessing the protein or ligand concentrations may also contribute to non-integer N values in ITC. To confirm that 300 mM NaCl did not negatively impact DNA binding, MarR_73 S28A was examined by ITC in 150 mM NaCl. In this condition, the KD for the 22 bp duplex was 0.428 ± 0.002 μM (N = 1.75 ± 0.014), while the KD for 24 bp duplex was 0.151 ± 0.025 μM (N = 2.51 ± 0.26).Protein crystallographyV. paradoxus MarR_73 was crystallized using the sitting drop vapour diffusion method at 20 °C in conditions outlined in Supplementary Table 4. Crystallization drops were set up using the Oryx4 protein crystallization robot (Douglas Instruments) and contained 0.15 μl protein and 0.15 μl well solution. For all V. paradoxus MarR_73 wild-type conditions, ligands were added at 10-fold molar excess before crystallization trials and crystals appeared within 2–5 d. V. paradoxus MarR_73 with the S28A and R46A mutations was crystallized in similar conditions as the wild-type protein. Similarly, P. putida MarR_iacR, B. japonicum MarR_Bj1, A. baumannii MarR_Ab and E. soli MarR_Es were crystallized using vapour diffusion methods in sitting drop trays at 20 °C and crystals appeared within 3–5 d. All crystallization conditions are outlined in Supplementary Table 4. Crystal specimens were cryoprotected with the well solution supplemented with glycerol to 20% (v/v) (Supplementary Table 4). X-ray diffraction data were collected at the Advanced Photon Source beamline 23-ID-D (Supplementary Table 3). Diffraction images were reduced using either XDS or Denzo and scaled with either Aimless or Scalepack57,58,59. The V. paradoxus MarR_73 structure in complex with IAA was determined by molecular replacement using the structure of 3CDH as a search model in Phaser60. All subsequent structures of V. paradoxus MarR_73 were determined using the V. paradoxus MarR_73 IAA complex structure (PDB: 7KFO) as a search model. The P. putida MarR_iacR and B. japonicum MarR_Bj1 structures were determined by molecular replacement using the structure of 3CJN as the search model. P. putida MarR_iacR (PDB: 7KUA) was subsequently used as the search model for molecular replacement to solve A. baumannii MarR_Ab and E. soli MarR_Es. A nickel ion was placed in the model of MarR_Ab. The following ions or molecules were examined and refined in this location in the MarR_Ab structure: water, Na, Mg, K, Ca, Mn, Fe, Co, Ni, Cu, Zn and Ba. Water, Na, Mg, K, Ca and Ba were deemed unacceptable in this site due to poor difference density. Of the remaining ions considered, there were no sources of Mn, Fe, Co, Cu or Zn in the protein expression media, protein purification buffers, protein storage buffer, crystallization condition or cryoprotectant solutions. Thus, we concluded that the ion present in this structure is Ni due to the use of a nickel-affinity column during the protein’s purification. It is unclear why this ion remained bound to MarR_Ab even after the subsequent size exclusion chromatography purification step, or why such an ion is only observed in this structure of the proteins examined. All structures were refined with either Phenix.refine or Refmac using iterative model building in Coot to the final parameters outlined in Supplementary Table 361,62. MarR_73 is a dimer with one protein monomer in the asymmetric unit and the dimer generated by crystallographic symmetry. PDB accession codes and associated crystallographic data are reported in Supplementary Table 3.Statistics and reproducibilityNo statistical method was used to predetermine sample size, but our sample sizes are similar to those reported in previous publications14,63,64. No data were excluded from the analyses. The experiments were randomized. The investigators were not blinded to allocation during experiments and outcome assessment. Where not stated, data distribution was assumed to be normal, but this was not formally tested.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Anderson, M. The evolution of eusociality. Annu. Rev. Ecol. Syst. 15, 165–189 (1984).Article 

Google Scholar 
Wilkinson, G. S. Reciprocal food sharing in the vampire bat. Nature 308, 181–184 (1984).Article 

Google Scholar 
Feistner, A. & Mcgrew, W. Food-sharing in primates: a critical review. Perspect. Primate Biol 3, (1989).Hoelzel, A. R. Killer whale predation on marine mammals at Punta Norte, Argentina; food sharing, provisioning and foraging strategy. Behav. Ecol. Sociobiol. 29, 197–204 (1991).Article 

Google Scholar 
Behmer, S. T. Animal behaviour: feeding the superorganism. Curr. Biol. 19, R366–R368 (2009).Article 
CAS 
PubMed 

Google Scholar 
Cassill, D. L. & Tschinkel, W. R. Information flow during social feeding in ant societies. in Information Processing in Social Insects (eds. Detrain, C., Deneubourg, J. L. & Pasteels, J. M.) 69–81 (Birkhäuser, 1999). https://doi.org/10.1007/978-3-0348-8739-7_4.Hunt, J. H. Trophallaxis and the evolution of eusocial Hymenoptera. in The Biology of Social Insects (CRC Press, 1982).Sorensen, A. A., Busch, T. M. & Vinson, S. B. Trophallaxis by temporal subcastes in the fire ant, Solenopsis invicta, in response to honey. Physiol. Entomol. 10, 105–111 (1985).Article 

Google Scholar 
Meurville, M.-P. & LeBoeuf, A. C. Trophallaxis: the functions and evolution of social fluid exchange in ant colonies (Hymenoptera: Formicidae). Myrmecol. News 31, 1–30 (2021).Bodner, L. et al. Nutrient utilization during male maturation and protein digestion in the Oriental hornet. Biology 11, 241 (2022).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Sorensen, A. A., Kamas, R. S. & Vinson, S. B. The influence of oral secretions from larvae on levels of proteinases in colony members of Solenopsis invicta Buren (Hymenoptera: Formicidae). J. Insect Physiol. 29, 163–168 (1983).Article 
CAS 

Google Scholar 
Erthal, M., Peres Silva, C. & Ian Samuels, R. Digestive enzymes in larvae of the leaf cutting ant, Acromyrmex subterraneus (Hymenoptera: Formicidae: Attini). J. Insect Physiol. 53, 1101–1111 (2007).Article 
CAS 
PubMed 

Google Scholar 
Went, F. W., Wheeler, J. & Wheeler, G. C. Feeding and digestion in some ants (Veromessor and Manica). BioScience 22, 82–88 (1972).Article 

Google Scholar 
Ishay, J. & Ikan, R. Food exchange between adults and larvae in Vespa orientalis F. Anim. Behav. 16, 298–303 (1968).Article 
CAS 
PubMed 

Google Scholar 
Hunt, J. H. The evolution of social wasps. (Oxford University Press, USA, 2007).Abe, T., Yoshiya, T., Hiromitsu, M. & Kawasaki, Y. Y. Comparative study of the composition of hornet larval saliva, its effect on behaviour and role of trophallaxis. Comp. Biochem. Physiol. Part C: Comp. Pharmacol. 99, 79–84 (1991).Article 

Google Scholar 
Ishay, J. & Ikan, R. Gluconeogenesis in the Oriental hornet Vespa orientalis F. Ecology 49, 169–171 (1968).Article 

Google Scholar 
Brock, R. E., Cini, A. & Sumner, S. Ecosystem services provided by aculeate wasps. Biol. Rev. 96, 1645–1675 (2021).Article 
PubMed 

Google Scholar 
Ueno, T. Flower-visiting by the invasive hornet Vespa velutina nigrithorax (Hymenoptera: Vespidae). Int. J. Chem., Environ. Biol. Sci. 3, 444–448 (2015).
Google Scholar 
Käfer, H., Kovac, H. & Stabentheiner, A. Respiration patterns of resting wasps (Vespula sp.). J. Insect Physiol. 59, 475–486 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Bodner, L., Bouchebti, S. & Levin, E. Allocation and metabolism of naturally occurring dietary amino acids in the Oriental hornet. Insect Biochem. Mol. Biol. 139, 103675 (2021).Article 
CAS 
PubMed 

Google Scholar 
Levin, E., Lopez-Martinez, G., Fane, B. & Davidowitz, G. Hawkmoths use nectar sugar to reduce oxidative damage from flight. Science 355, 733–735 (2017).Article 
CAS 
PubMed 

Google Scholar 
Hunt, J. H., Baker, I. & Baker, H. G. Similarity of amino acids in nectar and larval saliva: the nutritional basis for trophallaxis in social wasps. Evolution 36, 1318–1322 (1982).Article 
CAS 
PubMed 

Google Scholar 
Hunt, J. H., Jeanne, R. L., Baker, I. & Grogan, D. E. Nutrient dynamics of a swarm-founding social wasp species, Polybia occidentalis (Hymenoptera: Vespidae). Ethology 75, 291–305 (1987).Article 

Google Scholar 
Cassill, D. L. & Tschinkel, W. R. Allocation of liquid food to larvae via trophallaxis in colonies of the fire ant, Solenopsis invicta. Anim. Behav. 50, 801–813 (1995).Article 

Google Scholar 
Buffin, A., Denis, D., Simaeys, G. V., Goldman, S. & Deneubourg, J.-L. Feeding and stocking up: radio-labelled food reveals exchange patterns in ants. PLOS ONE 4, e5919 (2009).Article 
PubMed 
PubMed Central 

Google Scholar 
Quque, M. et al. Hierarchical networks of food exchange in the black garden ant Lasius niger. Insect Sci. 28, 825–838 (2021).Article 
PubMed 

Google Scholar 
Cassill, D. L. & Tschinkel, W. R. A duration constant for worker-to-larva trophallaxis in fire ants. Ins. Soc. 43, 149–166 (1996).Article 

Google Scholar 
Cassill, D. L. & Tschinkel, W. R. Regulation of diet in the fire ant, Solenopsis invicta. J. Insect Behav. 12, 307–328 (1999).Article 

Google Scholar 
Wilson, E. O. & Eisner, T. Quantitative studies of liquid food transmission in ants. Ins. Soc. 4, 157–166 (1957).Article 

Google Scholar 
Markin, G. P. Food distribution within laboratory colonies of the argentine ant, Tridomyrmex humilis (Mayr). Ins. Soc. 17, 127–157 (1970).Article 

Google Scholar 
Howard, D. F. & Tschinkel, W. R. The flow of food in colonies of the fire ant, Solenopsis invicta: a multifactorial study. Physiol. Entomol. 6, 297–306 (1981).Article 

Google Scholar 
Suryanarayanan, S. & Jeanne, R. L. Antennal drumming, trophallaxis, and colony development in the social wasp Polistes fuscatus (Hymenoptera: Vespidae). Ethology 114, 1201–1209 (2008).Article 

Google Scholar 
Greenwald, E., Segre, E. & Feinerman, O. Ant trophallactic networks: simultaneous measurement of interaction patterns and food dissemination. Sci. Rep. 5, 12496 (2015).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Baltiansky, L., Sarafian-Tamam, E., Greenwald, E. & Feinerman, O. Dual-fluorescence imaging and automated trophallaxis detection for studying multi-nutrient regulation in superorganisms. Methods Ecol. Evol. 12, 1441–1457 (2021).Article 

Google Scholar 
Feldhaar, H. et al. Stable isotopes: past and future in exposing secrets of ant nutrition (Hymenoptera: Formicidae). Myrmecol. N. 13, 3–13 (2010).
Google Scholar 
Bouchebti, S., Bodner, L., Bergman, M., Magory Cohen, T. & Levin, E. The effects of dietary proline, β-alanine, and γ-aminobutyric acid (GABA) on the nest construction behavior in the Oriental hornet (Vespa orientalis). Sci. Rep. 12, 7449 (2022).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Motro, M., Motro, U., Ishay, J. S. & Kugler, J. Some social and dietary prerequisites of oocyte development in Vespa orientalis L. workers. Ins. Soc. 26, 155–164 (1979).Article 

Google Scholar 
Levin, E., McCue, M. D. & Davidowitz, G. More than just sugar: allocation of nectar amino acids and fatty acids in a Lepidopteran. Proc. R. Soc. B: Biol. Sci. 284, 20162126 (2017).Article 

Google Scholar 
Wright, G. A., Nicolson, S. W. & Shafir, S. Nutritional physiology and ecology of honey bees. Annu. Rev. Entomol. 63, 327–344 (2018).Article 
CAS 
PubMed 

Google Scholar 
Helm, B. R. et al. The geometric framework for nutrition reveals interactions between protein and carbohydrate during larval growth in honey bees. Biol. Open 6, 872–880 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Paoli, P. P. et al. Nutritional balance of essential amino acids and carbohydrates of the adult worker honeybee depends on age. Amino Acids 46, 1449–1458 (2014).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Stabler, D., Paoli, P. P., Nicolson, S. W. & Wright, G. A. Nutrient balancing of the adult worker bumblebee (Bombus terrestris) depends on the dietary source of essential amino acids. J. Exp. Biol. 218, 793–802 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Arganda, S. et al. Parsing the life-shortening effects of dietary protein: effects of individual amino acids. Proc. R. Soc. B 284, 20162052 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Csata, E. & Dussutour, A. Nutrient regulation in ants (Hymenoptera: Formicidae): a review. Myrmecol. N. 29, 111–124 (2019).
Google Scholar 
Gottsberger, G., Schrauwen, J. & Linskens, H. F. Amino acids and sugars in nectar, and their putative evolutionary significance. Pl. Syst. Evol. 145, 55–77 (1984).Article 
CAS 

Google Scholar 
Ozimek, L. et al. Nutritive value of protein extracted from honey bees. J. Food Sci. 50, 1327–1329 (1985).Article 
CAS 

Google Scholar 
Nicolson, S. W. & Thornburg, R. W. Nectar chemistry. in Nectaries and Nectar (eds. Nicolson, S. W., Nepi, M. & Pacini, E.) 215–264 (Springer Netherlands, 2007). https://doi.org/10.1007/978-1-4020-5937-7_5.Contrera, F. A. L., Imperatriz-Fonseca, V. L. & Koedam, D. Trophallaxis and reproductive conflicts in social bees. Insect Soc. 57, 125–132 (2010).Article 

Google Scholar 
Carter, G. G. & Wilkinson, G. S. Food sharing in vampire bats: reciprocal help predicts donations more than relatedness or harassment. Proc. R. Soc. B: Biol. Sci. 280, 20122573 (2013).Article 

Google Scholar 
Nalepa, C. A. Origin of termite eusociality: trophallaxis integrates the social, nutritional, and microbial environments. Ecol. Entomol. 40, 323–335 (2015).Article 

Google Scholar 
Werenkraut, V., Arbetman, M. P. & Fergnani, P. N. The Oriental hornet (Vespa orientalis L.): a threat to the Americas? Neotrop. Entomol. 51, 330–338 (2022).Article 
PubMed 

Google Scholar 
Darchen, R. Biologie de Vespa orientalis. Les premiers stades de développement. Ins. Soc. 11, 141–157 (1964).Article 

Google Scholar 
Van Itterbeeck, J. et al. Rearing techniques for hornets with emphasis on Vespa velutina (Hymenoptera: Vespidae): A review. J. Asia-Pac. Entomol. 24, 103–117 (2021).Article 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/. (2020).Bouchebti, S., Bodner,L. & Levin, E. Continuous exchange of nectar nutrients in an Oriental hornet colony- Dataset [Data set]. Zenodo. https://doi.org/10.5281/zenodo.7135100 (2022). More

De Meyer, M., Robertson, M., Peterson, A. & Mansell, M. Ecological niches and potential geographical distributions of Mediterranean fruit fly (Ceratitis capitata) and Natal fruit fly (Ceratitis rosa). J. Biogeogr. 35, 270–281 (2008).
Google Scholar 
Nguyen, A. D. et al. Trade-offs in cold resistance at the northern range edge of the common woodland ant Aphaenogaster picea (Formicidae). Am. Nat. 194, E151–E163 (2019).Article 
PubMed 

Google Scholar 
Gilioli, G. et al. Non-linear physiological responses to climate change: the case of Ceratitis capitata distribution and abundance in Europe. Biol. Invasions 24, 261–279 (2022).Article 

Google Scholar 
Lancaster, L. T., Dudaniec, R. Y., Hansson, B. & Svensson, E. I. Latitudinal shift in thermal niche breadth results from thermal release during a climate-mediated range expansion. J. Biogeogr. 42, 1953–1963 (2015).Article 

Google Scholar 
Hallas, R., Schiffer, M. & Hoffmann, A. A. Clinal variation in Drosophila serrata for stress resistance and body size. Genet. Res. 79, 141–148 (2002).Article 
PubMed 

Google Scholar 
Hoffmann, A. A., Anderson, A. & Hallas, R. Opposing clines for high and low temperature resistance in Drosophila melanogaster. Ecol. Lett. 5, 614–618 (2002).Article 

Google Scholar 
Ragland, G. & Kingsolver, J. Influence of seasonal timing on thermal ecology and thermal reaction norm evolution in Wyeomyia smithii. J. Evol. Biol. 20, 2144–2153 (2007).Article 
CAS 
PubMed 

Google Scholar 
MacMillan, H. A. & Sinclair, B. J. Mechanisms underlying insect chill-coma. J. Insect Physiol. 57, 12–20 (2011).Article 
CAS 
PubMed 

Google Scholar 
Neilson, E. W. et al. There’sa storm a-coming: Ecological resilience and resistance to extreme weather events. Ecol. Evol. 10, 12147–12156 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Overgaard, J., Hoffmann, A. A. & Kristensen, T. N. Assessing population and environmental effects on thermal resistance in Drosophila melanogaster using ecologically relevant assays. J. Therm. Biol. 36, 409–416 (2011).Article 

Google Scholar 
Maysov, A. Chill coma temperatures appear similar along a latitudinal gradient, in contrast to divergent chill coma recovery times, in two widespread ant species. J. Exp. Biol. 217, 2650–2658 (2014).Article 
PubMed 

Google Scholar 
David, R. J. et al. Cold stress tolerance in Drosophila: analysis of chill coma recovery in D. melanogaster. J. therm. biol. 23, 291–299 (1998).Article 

Google Scholar 
Overgaard, J. & MacMillan, H. A. The integrative physiology of insect chill tolerance. Annu. Rev. Physiol. 79, 187–208 (2017).Article 
CAS 
PubMed 

Google Scholar 
Andersen, M. K. & Overgaard, J. The central nervous system and muscular system play different roles for chill coma onset and recovery in insects. Comp. Biochem. Physiol. A: Mol. Integr. Physiol. 233, 10–16 (2019).Article 
CAS 

Google Scholar 
Macdonald, S., Rako, L., Batterham, P. & Hoffmann, A. Dissecting chill coma recovery as a measure of cold resistance: evidence for a biphasic response in Drosophila melanogaster. J. Insect Physiol. 50, 695–700 (2004).Article 
CAS 
PubMed 

Google Scholar 
Gibert, P., Moreteau, B., Pétavy, G., Karan, D. & David, J. R. Chill-coma tolerance, a major climatic adaptation among Drosophila species. Evolution 55, 1063–1068 (2001).Article 
CAS 
PubMed 

Google Scholar 
Ayrinhac, A. et al. Cold adaptation in geographical populations of Drosophila melanogaster: phenotypic plasticity is more important than genetic variability. Funct. Ecol. 18, 700–706 (2004).Article 

Google Scholar 
Castañeda, L. E., Lardies, M. A. & Bozinovic, F. Interpopulational variation in recovery time from chill coma along a geographic gradient: a study in the common woodlouse, Porcellio laevis. J. Insect Physiol. 51, 1346–1351 (2005).Article 
PubMed 

Google Scholar 
Tonione, M. A., Cho, S. M., Richmond, G., Irian, C. & Tsutsui, N. D. Intraspecific variation in thermal acclimation and tolerance between populations of the winter ant Prenolepis imparis. Ecol. Evol. 10, 4749–4761 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Karl, I., Janowitz, S. A. & Fischer, K. Altitudinal life-history variation and thermal adaptation in the copper butterfly Lycaena tityrus. Oikos 117, 778–788 (2008).Article 

Google Scholar 
Ghalambor, C. K., Huey, R. B., Martin, P. R., Tewksbury, J. J. & Wang, G. Are mountain passes higher in the tropics? Janzen’s hypothesis revisited. Integr. Comp. Biol. 46, 5–17 (2006).Article 
PubMed 

Google Scholar 
Addo-Bediako, A., Chown, S. L. & Gaston, K. J. Thermal tolerance, climatic variability and latitude. In Proceedings of the Royal Society of London. Series B: Biological Sciences 267, 739–745 (2000).Poikela, N., Tyukmaeva, V., Hoikkala, A. & Kankare, M. Multiple paths to cold tolerance: the role of environmental cues, morphological traits and the circadian clock gene vrille. BMC ecol. Evol. 21, 1–20 (2021).
Google Scholar 
Andersen, J. L. et al. How to assess Drosophila cold tolerance: chill coma temperature and lower lethal temperature are the best predictors of cold distribution limits. Funct. Ecol. 29, 55–65 (2015).Article 

Google Scholar 
Papadopoulos, N., Katsoyannos, B., Carey, J. & Kouloussis, N. Seasonal and annual occurrence of the Mediterranean fruit fly (Diptera: Tephritidae) in northern Greece. Ann. Entomol. Soc. Am. 94, 41–50 (2001).Article 

Google Scholar 
Malacrida, A. et al. Globalization and fruitfly invasion and expansion: the medfly paradigm. Genetica 131, 1–9 (2007).Article 
CAS 
PubMed 

Google Scholar 
Egartner, A., Lethmayer, C., Gottsberger, R. A. & Blümel, S. In Joint Meeting of the IOBC-WPRS Working Groups “Pheromones and other semiochemicals in integrated production” & “Integrated Protection of Fruit Crops” at. 143–152.Nyamukondiwa, C., Kleynhans, E. & Terblanche, J. S. Phenotypic plasticity of thermal tolerance contributes to the invasion potential of mediterranean fruit flies (Ceratitis capitata). Ecol. Entomol. 35, 565–575 (2010).Article 

Google Scholar 
Weldon, C. W., Terblanche, J. S. & Chown, S. L. Time-course for attainment and reversal of acclimation to constant temperature in two Ceratitis species. J. Therm. Biol. 36, 479–485 (2011).Article 

Google Scholar 
Pujol-Lereis, L. M., Rabossi, A. & Quesada-Allué, L. A. Analysis of survival, gene expression and behavior following chill-coma in the medfly Ceratitis capitata: effects of population heterogeneity and age. J. Insect Physiol. 71, 156–163 (2014).Article 
CAS 
PubMed 

Google Scholar 
Pujol-Lereis, L. M., Fagali, N. S., Rabossi, A., Catalá, Á. & Quesada-Allué, L. A. Chill-coma recovery time, age and sex determine lipid profiles in Ceratitis capitata tissues. J. Insect Physiol. 87, 53–62 (2016).Article 
CAS 
PubMed 

Google Scholar 
Weldon, C. W., Nyamukondiwa, C., Karsten, M., Chown, S. L. & Terblanche, J. S. Geographic variation and plasticity in climate stress resistance among southern African populations of Ceratitis capitata (Wiedemann)(Diptera: Tephritidae). Sci. Rep. 8, 1–13 (2018).Article 
CAS 

Google Scholar 
Nyamukondiwa, C., Weldon, C. W., Chown, S. L., le Roux, P. C. & Terblanche, J. S. Thermal biology, population fluctuations and implications of temperature extremes for the management of two globally significant insect pests. J. Insect Physiol. 59, 1199–1211 (2013).Article 
CAS 
PubMed 

Google Scholar 
Mitchell, K. A., Boardman, L., Clusella-Trullas, S. & Terblanche, J. S. Effects of nutrient and water restriction on thermal tolerance: A test of mechanisms and hypotheses. Comp. Biochem. Physiol. A: Mol. Integr. Physiol. 212, 15–23 (2017).Article 
CAS 

Google Scholar 
Hoffmann, A. A. & Ross, P. A. Rates and patterns of laboratory adaptation in (mostly) insects. J. Econ. Entomol. 111, 501–509 (2018).Article 
PubMed 

Google Scholar 
Popa-Báez, Á. -D. et al. Climate stress resistance in male Queensland fruit fly varies among populations of diverse geographic origins and changes during domestication. BMC Genet. 21, 1–19 (2020).Article 

Google Scholar 
Beck, H. E. et al. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 5, 1–12 (2018).Article 

Google Scholar 
Kozak, K. H., Graham, C. H. & Wiens, J. J. Integrating GIS-based environmental data into evolutionary biology. Trends Ecol. Evol. 23, 141–148 (2008).Article 
PubMed 

Google Scholar 
Oyen, K. J. et al. Body mass and sex, not local climate, drive differences in chill coma recovery times in common garden reared bumble bees. J. Comp. Physiol. B. 191, 843–854 (2021).Article 
PubMed 

Google Scholar 
Angert, A. L., Bontrager, M. G. & Ågren, J. What do we really know about adaptation at range edges?. Annu. Rev. Ecol. Evol. Syst. 51, 341–361 (2020).Article 

Google Scholar 
Terblanche, J. S. & Hoffmann, A. A. Validating measurements of acclimation for climate change adaptation. Curr. Opin. insect sci. 41, 7–16 (2020).Article 
PubMed 

Google Scholar 
Kourti, A. Patterns of variation within and between Greek populations of Ceratitis capitata suggest extensive gene flow and latitudinal clines. J. Econ. Entomol. 97, 1186–1190 (2004).Article 
CAS 
PubMed 

Google Scholar 
Hangartner, S., Lasne, C., Sgrò, C. M., Connallon, T. & Monro, K. Genetic covariances promote climatic adaptation in Australian Drosophila. Evolution 74, 326–337 (2020).Article 
PubMed 

Google Scholar 
Bontrager, M. & Angert, A. L. Gene flow improves fitness at a range edge under climate change. Evol. Let. 3, 55–68 (2019).Article 

Google Scholar 
Liu, Q. et al. Extension of the growing season increases vegetation exposure to frost. Nat. Commun. 9, 1–8 (2018).
Google Scholar 
Schwartz, M. D., Ahas, R. & Aasa, A. Onset of spring starting earlier across the Northern Hemisphere. Glob. Change Biol. 12, 343–351 (2006).Article 

Google Scholar 
Ma, Q., Huang, J. G., Hänninen, H. & Berninger, F. Divergent trends in the risk of spring frost damage to trees in Europe with recent warming. Glob. Change Biol. 25, 351–360 (2019).Article 

Google Scholar 
Unterberger, C. et al. Spring frost risk for regional apple production under a warmer climate. PLoS ONE 13, e0200201 (2018).Article 
MathSciNet 
PubMed 
PubMed Central 

Google Scholar 
Manrakhan, A., Daneel, J.-H., Stephen, P. R. & Hattingh, V. Cold Tolerance of Immature Stages of Ceratitis capitata and Bactrocera dorsalis (Diptera: Tephritidae). J. Econ. Entomol. 115(2), 482–492 (2022).Article 
PubMed 

Google Scholar 
Papadopoulos, N. T., Carey, J. R., Katsoyannos, B. I. & Kouloussis, N. A. Overwintering of the mediterranean fruit fly (Diptera: Tephritidae) in Northern Greece. Ann. Entomol. Soc. Am. 89, 526–534 (1996).Article 

Google Scholar 
Papadopoulos, N. T., Katsoyannos, B. I. & Carey, J. R. Temporal changes in the composition of the overwintering larval population of the Mediterranean fruit fly (Diptera: Tephritidae) in Northern Greece. Ann. Entomol. Soc. Am. 91, 430–434 (1998).Article 

Google Scholar 
Katsoyannos, B. I., Kouloussis, N. A. & Carey, J. R. Seasonal and annual occurrence of Mediterranean fruit flies (Diptera: Tephritidae) on Chios Island, Greece: Differences between two neighboring citrus orchards. Ann. Entomol. Soc. Am. 91, 43–51 (1998).Article 

Google Scholar 
Mavrikakis, P. G., Economopoulos, A. P. & Carey, J. R. Continuous winter reproduction and growth of the mediterranean fruit fly (Diptera: Tephritidae) in Heraklion, crete Southern Greece. Environ. Entomol. 29, 1180–1187 (2000).Article 

Google Scholar 
Israely, N., Ziv, Y. & Oman, S. D. Spatiotemporal distribution patterns of Mediterranean fruit fly (Diptera: Tephritidae) in the central region of Israel. Ann. Entomol. Soc. Am. 98, 77–84 (2005).Article 

Google Scholar 
Bahrndorff, S., Lauritzen, J. M., Sørensen, M. H., Noer, N. K. & Kristensen, T. N. Responses of terrestrial polar arthropods to high and increasing temperatures. J. Exp. Biol. 224, jeb230797 (2021).Article 
PubMed 

Google Scholar 
Sinclair, B. J. & Roberts, S. P. Acclimation, shock and hardening in the cold. J. Therm. Biol. 30, 557–562 (2005).Article 

Google Scholar 
Bahrndorff, S., Gertsen, S., Pertoldi, C. & Kristensen, T. N. Investigating thermal acclimation effects before and after a cold shock in Drosophila melanogaster using behavioural assays. Biol. J. Lin. Soc. 117, 241–251 (2016).Article 

Google Scholar 
Sarmad, M., Ishfaq, A., Arif, H. & Zaka, S. M. Effect of short-term cold temperature stress on development, survival and reproduction of Dysdercus koenigii (Hemiptera: Pyrrhocoridae). Cryobiology 92, 47–52 (2020).Article 
PubMed 

Google Scholar 
Steyn, V. M., Mitchell, K. A., Nyamukondiwa, C. & Terblanche, J. S. Understanding costs and benefits of thermal plasticity for pest management: Insights from the integration of laboratory, semi-field and field assessments of Ceratitis capitata (Diptera: Tephritidae). Bull. Entomol. Res., 1–11 (2022).Davis, H. E., Cheslock, A. & MacMillan, H. A. Chill coma onset and recovery fail to reveal true variation in thermal performance among populations of Drosophila melanogaster. Sci. Rep. 11, 1–10 (2021).Article 

Google Scholar 
Noh, S., Everman, E. R., Berger, C. M. & Morgan, T. J. Seasonal variation in basal and plastic cold tolerance: Adaptation is influenced by both long-and short-term phenotypic plasticity. Ecol. Evol. 7, 5248–5257 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Bruins, H. J. Ancient desert agriculture in the Negev and climate-zone boundary changes during average, wet and drought years. J. Arid Environ. 86, 28–42 (2012).Article 

Google Scholar 
Hoffmann, A. A., Sørensen, J. G. & Loeschcke, V. Adaptation of Drosophila to temperature extremes: Bringing together quantitative and molecular approaches. J. Therm. Biol. 28, 175–216 (2003).Article 

Google Scholar 
Kawecki, T. J. & Ebert, D. Conceptual issues in local adaptation. Ecol. Lett. 7, 1225–1241 (2004).Article 

Google Scholar 
Nyamukondiwa, C. & Terblanche, J. S. Thermal tolerance in adult mediterranean and Natal fruit flies (Ceratitis capitata and Ceratitis rosa): Effects of age, gender and feeding status. J. Therm. Biol. 34, 406–414 (2009).Article 

Google Scholar 
Team, R. C. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria (2021).Fick, S. E. & Hijmans, R. J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).Article 

Google Scholar 
Mazerolle, M. J. Model selection and multimodel inference using the AICcmodavg package (2020).Therneau, T. A Package for Survival Analysis in R. R Package Version 3.2-13.(2021. (2021).Kassambara, A., Kosinski, M., Biecek, P. & Fabian, S. Survminer: Drawing Survival Curves using’ggplot2′. R package version 0.4. 9. 2021. (2021).Lenth, R. V. Emmeans: Estimated Marginal Means, aka Least-Squares Means. R package version 1.7.2. (2022). More

Sarwar, N. et al. Phytoremediation strategies for soils contaminated with heavy metals: Modifications and future perspectives. Chemosphere 171, 710–721 (2017).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Lin, H. M. et al. Cadmium-stress mitigation through gene expression of rice and silicon addition. Plant Growth Regul.: Int. J. Nat. Synthetic Regul. 81(1), 91–101 (2017).Article 
CAS 

Google Scholar 
Pan, F. S. et al. Enhanced Cd extraction of oilseed rape (Brassica napus) by plant growth-promoting bacteria isolated from Cd hyperaccumulator Sedum alfredii Hance. Int. J. Phytorem. 19(1/6), 281–289 (2017).Article 
CAS 

Google Scholar 
Puangprasert, S. & Prueksasit, T. Health risk assessment of airborne Cd, Cu, Ni and Pb for electronic waste dismantling workers in Buriram Province, Thailand. J. Environ. Manag. 252, 109601 (2019).Article 
CAS 

Google Scholar 
Tipu, M. I. et al. Growth and physiology of maize (Zea mays L.) in a nickel-contaminated soil and phytoremediation efficiency using EDTA. J. Plant Growth Regul. 40(2), 774–786 (2021).Article 
CAS 

Google Scholar 
Chaturvedi, N., Dhal, N. K. & Patra, H. K. EDTA and citric acid-mediated phytoextraction of heavy metals from iron ore tailings using Andrographis paniculata: A comparative study. Int. J. Min. Reclam. Environ. 29(1), 33–46 (2015).Article 
CAS 

Google Scholar 
Wang, G. Y. et al. Heavy metal removal by GLDA washing: Optimization, redistribution, recycling, and changes in soil fertility. Sci. Total Environ. 569–570, 557–568 (2016).Article 
ADS 
PubMed 

Google Scholar 
Kołodyńska, D. Cu(II), Zn(II), Co(II) and Pb(II) removal in the presence of the complexing agent of a new generation. Desalination 267(2–3), 175–183 (2011).Article 

Google Scholar 
Guo, X. F. et al. Mixed chelators of EDTA, GLDA, and citric acid as washing agent effectively remove Cd, Zn, Pb, and Cu from soils. J. Soils Sediments 18(2), 835–844 (2017).
Google Scholar 
Wang, X. et al. Subcellular distribution and chemical forms of cadmiun in Bechmeria nivea L. Gaud. Environ. Exp. Bot. 62(3), 389–395 (2008).Article 
CAS 

Google Scholar 
Gallego, S. M. et al. Unravelling cadmium toxicity and tolerance in plants: Insight into regulatory mechanisms. Environ. Exp. Bot. 83, 33–46 (2012).Article 
CAS 

Google Scholar 
Clemens, S., Aarts, M. G. M., Thomine, S. & Verbruggen, N. Plant science: The key to preventing slow cadmium poisoning. Trends Plant Sci. 18(2), 92–99 (2013).Article 
CAS 
PubMed 

Google Scholar 
Zhou, J. T. et al. Integration of cadmium accumulation, subcellular distribution, and physiological responses to understand cadmium tolerance in apple rootstocks. Front. Plant Sci. 8, 966 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Yang, L. P., Zhu, J., Wang, P., Lyu, D. G. & Li, H. F. Effect of Cd on growth, physiological response, Cd subcellular distribution and chemical forms of Koelreuteria paniculata. Ecotoxicol. Environ. Saf. 160, 10–18 (2018).Article 
CAS 
PubMed 

Google Scholar 
Wang, W. J., Zhang, M. Z. & Liu, J. N. Subcellular distribution and chemical forms of Cd in Bougainvillea spectabilis Willd. as an ornamental phytostabilizer: An integrated consideration. Int. J. Phytorem. 20(11), 1087–1095 (2017).Article 

Google Scholar 
Weigel, H. J. & Jäger, H. J. Subcellular distribution and chemical form of cadmium in bean plants. Plant Physiol. 65(3), 480–482 (1980).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Khanna, K., Kohli, S. K., Ohri, P., Bhardwaj, R. & Ahmad, P. Agroecotoxicological aspect of Cd in soil–plant system: Uptake, translocation and amelioration strategies. Environ. Sci. Pollut. Res. 29, 30908–30934 (2022).Article 
CAS 

Google Scholar 
Wei, Z. B., Chen, X. H., Wu, Q. T. & Tan, M. Biodegradable chelator GLDA induced remediation of heavy metal contaminated soil in Southeast Jingtian. Environ. Sci. 36(5), 1864–1869 (2015).CAS 

Google Scholar 
Wang, K., Liu, Y. H., Song, Z. G., Wang, D. & Qiu, W. W. Chelator complexes enhanced Amaranthus hypochondriacus L. phytoremediation efficiency in Cd-contaminated soils. Chemosphere 237, 124480 (2019).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Meng, N., Wang, M., Chen, L., Zheng, H. & Chen, S. B. Remediation effects of different herbaceous plants intercropping on Cd-contaminated soil. China Environ. Sci. 38(7), 2618–2624 (2018).CAS 

Google Scholar 
Jones, D. & Willett, V. Experimental evaluation of methods to quantify dissolved organic nitrogen (don) and dissolved organic carbon (doc) in soil. Soil Biol. Biochem. 38(5), 991–999 (2006).Article 
CAS 

Google Scholar 
Su, F. L. et al. The distribution and enrichment characteristics of copper in soil and Phragmites australis of Liao River estuary wetland. Environ. Monit. Assess.: Int. J. 190(6), 1–9 (2018).Article 
CAS 

Google Scholar 
Shahid, M., Dumat, C. & Khalid, S. Reviews of Environmental Contamination and Toxicology Vol. 241, 3–137 (Springer, 2016).
Google Scholar 
Yuliya, V. et al. Comparison of soil-to-root transfer and translocation coefficients of trace elements in vines of Chardonnay and Muscat white grown in the same vineyard. Sci. Hortic. 192, 89–96 (2015).Article 

Google Scholar 
Liu, Q. Q., Chen, Y. H., Shen, Z. G. & Zheng, L. Q. Roles of cell wall in plant heavy metal tolerance. Plant Physiol. J. 50(5), 605–611 (2014).
Google Scholar 
Zhen, S. et al. Foliar application of Zn reduces Cd accumulation in grains of late rice by regulating the antioxidant system, enhancing Cd chelation onto cell wall of leaves, and inhibiting Cd translocation in rice. Sci. Total Environ. 770, 145302 (2021).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Shi, Y. X. et al. Simulation of the absorption, migration and accumulation process of heavy metal elements in soil-crop system. Environ. Sci. 37(10), 3996–4003 (2016).
Google Scholar 
Yan, X. X. et al. Effect of foliar application of different manganese fertilizers on cadmium accumulation and subcellular distribution in pak choi. J. Agro Environ. Sci. 38(8), 1872–1881 (2019).
Google Scholar 
He, S., Wu, Q. & He, Z. Effect of DA-6 and EDTA alone or in combination on uptake, subcellular distribution and chemical form of Pb in Lolium perenne. Chemosphere 93(11), 2782–2788 (2013).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Li, C. C. et al. Integration of metal chemical forms and subcellular partitioning to understand metal toxicity in two lettuce (Lactuca sativa L.) cultivars. Plant Soil 384(1/2), 201–212 (2014).Article 
CAS 

Google Scholar 
Li, D., He, T., Saleem, M. & He, G. Metalloprotein-specific or critical amino acid residues: Perspectives on plant-precise detoxification and recognition mechanisms under cadmium stress. Int. J. Mol. Sci. 23(3), 1734 (2022).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Perriguey, J., Sterckeman, T. & Morel, J. L. Effect of rhizosphere and plantrelated factors on the cadmium uptake by maize(Zea mays L.). Environ. Exp. Bot. 63(1/3), 333–341 (2008).Article 
CAS 

Google Scholar 
Dai, S. et al. Effects of biochar amendments on speciation and bioavailability of heavy metals in coal-mine-contaminated soil. Hum. Ecol. Risk Assess. Int. J. 24(7), 1887–1900 (2018).Article 
CAS 

Google Scholar 
Hou, S., Zheng, N., Tang, L., Ji, X. F. & Li, Y. Y. Effect of soil pH and organic matter content on heavy metals availability in maize (Zea mays L.) rhizospheric soil of non-ferrous metals smelting area. Environ. Monit. Assess. 191(10), 634 (2019).Article 
PubMed 

Google Scholar 
Wu, H. J. et al. Effects of Astragalus smicuson cadmium effectiveness in paddy soil and cadmium accumulation in rice plant. Chin. Agric. Sci. Bull. 33(16), 105–111 (2017).ADS 

Google Scholar 
Jin, P. K., Liu, K. J. & Wang, X. B. Conversion and utilization of slowly biodegradable organic matter. Chin. J. Environ. Eng. 10(5), 2168–2174 (2016).CAS 

Google Scholar 
Kopáček, J. et al. Factors affecting the leaching of dissolved organic carbon after tree dieback in an unmanaged European mountain forest. Environ. Sci. Technol. 52(11), 6291–6299 (2018).Article 
ADS 
PubMed 

Google Scholar 
Anwar, S. et al. Impact of chelator-induced phytoextraction of cadmium on yield and ionic uptake of maize. Int. J. Phytorem. 19(6), 505–513 (2017).Article 
CAS 

Google Scholar 
Wu, J. M., Xi, M. & Kong, F. L. Review of researches on the factors influencing the dynamics of dissolved organic carbon in soils. Geol. Rev. 59(5), 953–961 (2013).CAS 

Google Scholar 
AkzoNobel. Dissolvine GL® Technichal Brochure 1–5 (AkzoNobel Amsterdam, 2010).
Google Scholar 
Beygi, M. & Jalali, M. Assessment of trace elements (Cd, Cu, Ni, Zn) fractionation and bioavailability in vineyard soils from the Hamedan, Iran. Geoderma 337, 1009–1020 (2019).Article 
ADS 
CAS 

Google Scholar 
Gul, I. et al. Comparative effectiveness of organic and inorganic amendments on cadmium bioavailability and uptake by Pelargonium hortorum. J. Soils Sediments 19(5), 2346–2356 (2019).Article 
CAS 

Google Scholar 
Wang, H., Sun, L. N., Li, H. B. & Sun, T. Y. Effect of different chelators application on Cd accumulation in metal polluted soils by Beta vulgaris var. cicla L. Ecol. Environ. 17(6), 2249–2252 (2008).
Google Scholar 
Zhang, G. X. et al. Effects of biochars on the availability of heavy metals to ryegrass in an alkaline contaminated soil. Environ. Pollut. 218, 513–522 (2016).Article 
CAS 
PubMed 

Google Scholar 
Gu, M. H. et al. Effects of manganese application on the formation of manganese oxides and cadmium fixation in soil. Ecol. Environ. Sci. 229(2), 360–368 (2020).
Google Scholar 
Bradl, H. B. Adsorption of heavy metal ions on soils and soils constituents. J. Colloid Interface Sci. 277(1), 1–18 (2004).Article 
ADS 
CAS 
PubMed 

Google Scholar  More

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.This is a summary of: Marin, F. R. et al. Protecting the Amazon forest and reducing global warming via agricultural intensification. Nat. Sustain. https://doi.org/10.1038/s41893-022-00968-8 (2022). More

Longline fishing boats such as these at Honolulu’s harbour in Hawaii must respect a large no-fishing zone off the western side of the archipelago.Credit: Sarah Medoff

Large no-fishing areas can drive the recovery of commercially valuable fish species, a study suggests. Ten years’ worth of fisheries data have shown that catch rates of two important types of tuna increased drastically in the vicinity of a marine protected area surrounding the northwestern Hawaiian islands.“It’s a win–win for fish and fishermen,” says Jennifer Raynor, an economist at the University of Wisconsin–Madison and a co-author of the study, which was published on 20 October in Science1.The results highlight the value of large-scale marine protected areas — a type of environmental management that has emerged in the past two decades, mostly in the Pacific Ocean, says Kekuewa Kikiloi, who studies Hawaiian culture at the University of Hawaii at Mānoa. Countries around the world have committed to protecting 30% of their land and oceans by 2030.Previous research showed that marine protected areas can help to restore populations of creatures that don’t move around much or at all, such as corals2 and lobsters3. Raynor and her colleagues wanted to test whether the areas could also drive the recovery of migratory species and provide spillover benefits for fisheries. The researchers looked at one of the largest such areas in the world, the 1.5-million-square-kilometre Papahānaumokuākea Marine National Monument, which was created in 2006 and expanded in 2016 to protect biological and cultural resources.The team focused on the Hawaiian ‘deep-set’ longline fishery, which mainly targets yellowfin tuna (Thunnus albacares) and bigeye tuna (Thunnus obesus).The researchers analysed catch data collected on fishing vessels between 2010 and late 2019. Then, they compared catch rates at various distances up to 600 nautical miles (1,111 kilometres) from the protected area, before and after its expansion in 2016. (The protected area itself currently extends for 200 nautical miles from the northwestern part of the Hawaiian archipelago.) They found that after the expansion, catch rates — defined as the number of fish caught for every 1,000 hooks deployed — went up, and that the increases were greater the closer the boats were to the no-fishing zone. At distances of up to 100 nautical miles, the catch rate for yellowfin tuna increased by 54%, and that for bigeye tuna by 12%. Some other types of catch rate also increased, but not by equally significant margins.The size of the Papahānaumokuākea Marine National Monument — more than three times the surface area of California — probably played a part in the positive effects, as did its shape. It spans about 2,000 kilometres from west to east, protecting large swathes of ocean waters at tropical latitudes. This means that tropical fish such as yellowfin and bigeye tuna — which tend to move along an east–west axis to stay in their preferred temperature range — can travel a long way and still stay in the no-fishing zone.What’s more, says Raynor, Papahānaumokuākea is a spawning ground for yellowfin tuna. Because the animals don’t travel far from their birthplace, the no-take zone provides refuge from fishing, helping tuna to aggregate and reproduce.“It is exciting to see that there are benefits to the fishing industry from this marine protected area,” says David Kroodsma, director of research and innovation at Global Fishing Watch in Oakland, California, a US non-governmental organization that monitors fishing activity worldwide. However, he adds, it’s unclear whether the results can be generalized to other areas of the world.Regardless, the findings could help others to design marine protected areas so that benefits trickle down to fisheries, says Steve Gaines, a marine ecologist at the University of California, Santa Barbara. The study, he says, “provides a platform to definitively evaluate what is working and what isn’t”.Co-managed by Indigenous populations, the state of Hawaii and the US government, Papahānaumokuākea is an example of a collaborative management strategy that bridges Indigenous knowledge and modern science, Kikiloi says. The approach, he adds, “can work successfully in other places too, if given a chance”. More

Field samplingLake Masoko fish were chased into fixed gill nets and SCUBA by a team of professional divers at different target depths determined by diver depth gauge (12× male benthic, 12× male littoral). Riverine fish (11× Mbaka River and 1× Itupi river) were collected by local fishermen. On collection, all fish were euthanized using clove oil. Collection of wild fish was done in accordance with local regulations and permits in 2015, 2016, 2018 and 2019. On collection, fish were immediately photographed with color and metric scales, and tissues were dissected and stored in RNAlater (Sigma-Aldrich); some samples were first stored in ethanol. Only male specimens (showing bright nuptial coloration) were used in this study for the practical reason of avoiding any misassignment of individuals to the wrong population (only male individuals show clear differences in phenotypes and could therefore be reliably assigned to a population). Furthermore, we assumed that any epigenetic divergence relevant to speciation should be contributing to between-population differences in traits possessed by both sexes (habitat occupancy, diet). To investigate the role of epigenetics in phenotypic diversification and adaptation to different diets, homogenized liver tissue – a largely homogenous and key organ involved in dietary metabolism, hormone production and hematopoiesis – was used for all RNA-seq and WGBS experiments.Common-garden experimentCommon-garden fish were bred from wild-caught fish specimens, collected and imported at the same time by a team of professional aquarium fish collectors according to approved veterinary regulations of the University of Bangor, UK. Wild-caught fish were acclimatized to laboratory tanks and reared to produce first-generation (G1) common-garden fish, which were reared under the same controlled laboratory conditions in separate tanks (light–dark cycles, diet: algae flakes daily, 2–3 times weekly frozen diet) for approximately 6 months (post hatching). G1 adult males showing bright nuptial colors were culled at the same biological stages (6 months post hatching) using MS222 in accordance with the veterinary regulations of the University of Bangor, UK. Immediately on culling, fish were photographed and tissues collected and snap-frozen in tubes.Stable isotopesTo assess dietary/nutritional profiles in the three ecomorph populations, carbon (δ13C) and nitrogen (δ15N) isotope analysis of muscle samples (for the same individuals as RRBS; 12, 12 and 9 samples for benthic, littoral and riverine populations, respectively) was undertaken by elemental analyzer isotope ratio mass spectrometry by Iso-Analytical Limited. It is important to note that stable isotope analysis does not depend on the use of the same tissue as the ones used for the RRBS/WGBS samples45. Normality tests (Shapiro–Wilk, using the R package rstatix v.0.7.0), robust for small sample sizes, were performed to assess sample deviation from a Gaussian distribution. Levene’s test for homogeneity of variance was then performed (R package carData v.3.0-5) to test for homogeneity of variance across groups. Finally, Welch’s ANOVA was performed followed by Games–Howell all-pairs comparison tests with adjusted P value using Tukey’s method (rstatix v.0.7.0). Mean differences in isotope measurements and 95% CI mean differences were calculated using Dabestr v.0.3.0 with 5,000 bootstrapped resampling.Throughout this manuscript, all box plots are defined as follows: centre line, median; box limits, upper and lower quartiles; whiskers, 1.5× interquartile range; points, outliers.RNA-seqNext-generation sequencing library preparationTotal RNA from liver tissues stored in RNAlater was extracted using a phenol/chloroform approach (TRIzol reagent; Sigma-Aldrich). Of note, when tissues for bisulphite sequencing samples were not available, additional wild-caught samples were used (Supplementary Table 3). The quality and quantity of RNA extraction were assessed using TapeStation (Agilent Technologies), Qubit and NanoDrop (Thermo Fisher Scientific). Next-generation sequencing (NGS) libraries were prepared using poly(A) tail-isolated RNA fraction and sequenced on a NovaSeq system (S4; paired-end 100/150 bp; Supplementary Table 3), yielding on average 32.9 ± 3.9 Mio reads.Read alignment and differential gene expression analysisAdaptor sequence in reads, low-quality bases (Phred score  More