Ecology

Tilman, D., Balzer, C., Hill, J. & Befort, B. L. Global food demand and the sustainable intensification of agriculture. Proc. Natl Acad. Sci. USA 108, 20260–20264 (2011).CAS 
Article 

Google Scholar 
Mueller, N. D. et al. Closing yield gaps through nutrient and water management. Nature 490, 254–257 (2012).CAS 
Article 

Google Scholar 
Hong, C. et al. Global and regional drivers of land-use emissions in 1961–2017. Nature 589, 554–561 (2021).CAS 
Article 

Google Scholar 
Laurance, W. F., Sayer, J. & Cassman, K. G. Agricultural expansion and its impacts on tropical nature. Trends Ecol. Evol. 29, 107–116 (2014).Article 

Google Scholar 
Cassman, K. G. & Grassini, P. A global perspective on sustainable intensification research. Nat. Sustain. 3, 262–268 (2020).Article 

Google Scholar 
Hodge, I., Hauck, J. & Bonn, A. The alignment of agricultural and nature conservation policies in the European Union. Conserv. Biol. 29, 996–1005 (2015).Article 

Google Scholar 
Heilmayr, R., Rausch, L. L., Munger, J. & Gibbs, H. K. Brazil’s Amazon Soy Moratorium reduced deforestation. Nat. Food 1, 801–810 (2020).Article 

Google Scholar 
Diffenbaugh, N. S. et al. Quantifying the influence of global warming on unprecedented extreme climate events. Proc. Natl Acad. Sci. USA 114, 4881–4886 (2017).CAS 
Article 

Google Scholar 
Iizumi, T. & Ramankutty, N. How do weather and climate influence cropping area and intensity? Glob. Food Security 4, 46–50 (2015).Article 

Google Scholar 
Davis, K. F., Downs, S. & Gephart, J. A. Towards food supply chain resilience to environmental shocks. Nat. Food 2, 54–65 (2020).Article 

Google Scholar 
Wang, X. et al. Emergent constraint on crop yield response to warmer temperature from field experiments. Nat. Sustain. 3, 908–916 (2020).Article 

Google Scholar 
Lobell, D. B., Schlenker, W. & Costa-Roberts, J. Climate trends and global crop production since 1980. Science 333, 616–620 (2011).CAS 
Article 

Google Scholar 
Sloat, L. L. et al. Climate adaptation by crop migration. Nat. Commun. 11, 1243 (2020).CAS 
Article 

Google Scholar 
Afifi, T., Liwenga, E. & Kwezi, L. Rainfall-induced crop failure, food insecurity and out-migration in Same-Kilimanjaro, Tanzania. Clim. Dev. 6, 53–60 (2014).Article 

Google Scholar 
Stigter, K. in Applied Agrometeorology (ed. Stigter, K.) 531–534 (Springer, 2010).Seifert, C. A. & Lobell, D. B. Response of double cropping suitability to climate change in the United States. Environ. Res. Lett. 10, 024002 (2015).Article 

Google Scholar 
Kawasaki, K. Two harvests are better than one: double cropping as a strategy for climate change adaptation. Am. J. Agr. Econ. 101, 172–192 (2019).Article 

Google Scholar 
Ceglar, A., Zampieri, M., Toreti, A. & Dentener, F. Observed northward migration of agro‐climate zones in Europe will further accelerate under climate change. Earths Future 7, 1088–1101 (2019).Article 

Google Scholar 
Cohn, A. S., VanWey, L. K., Spera, S. A. & Mustard, J. F. Cropping frequency and area response to climate variability can exceed yield response. Nat. Clim. Change 6, 601–604 (2016).Article 

Google Scholar 
Challinor, A. J., Simelton, E. S., Fraser, E. D. G., Hemming, D. & Collins, M. Increased crop failure due to climate change: assessing adaptation options using models and socio-economic data for wheat in China. Environ. Res. Lett. 5, 034012 (2010).Article 

Google Scholar 
Ray, D. K. & Foley, J. A. Increasing global crop harvest frequency: recent trends and future directions. Environ. Res. Lett. 8, 044041 (2013).Article 

Google Scholar 
Wu, W. et al. Global cropping intensity gaps: increasing food production without cropland expansion. Land Use Policy 76, 515–525 (2018).Article 

Google Scholar 
Pugh, T. A. M. et al. Climate analogues suggest limited potential for intensification of production on current croplands under climate change. Nat. Commun. 7, 12608 (2016).CAS 
Article 

Google Scholar 
Scherer, L. A., Verburg, P. H. & Schulp, C. J. E. Opportunities for sustainable intensification in European agriculture. Glob. Environ. Change 48, 43–55 (2018).Article 

Google Scholar 
Qin, Y. et al. Agricultural risks from changing snowmelt. Nat. Clim. Change 10, 459–465 (2020).Article 

Google Scholar 
Waha, K. et al. Multiple cropping systems of the world and the potential for increasing cropping intensity. Glob. Environ. Change 64, 102131 (2020).Article 

Google Scholar 
Raderschall, C. A., Vico, G., Lundin, O., Taylor, A. R. & Bommarco, R. Water stress and insect herbivory interactively reduce crop yield while the insect pollination benefit is conserved. Glob. Chang. Biol. 27, 71–83 (2021).CAS 
Article 

Google Scholar 
Ding, M. et al. Variation in cropping intensity in Northern China from 1982 to 2012 based on GIMMS-NDVI data. Sustainability 8, 1123 (2016).Article 

Google Scholar 
Yu, Q., Xiang, M., Sun, Z. & Wu, W. The complexity of measuring cropland use intensity: an empirical study. Agr. Syst. 192, 103180 (2021).Article 

Google Scholar 
Moore, F. C. & Lobell, D. B. Adaptation potential of European agriculture in response to climate change. Nat. Clim. Change 4, 610–614 (2014).Article 

Google Scholar 
Agnolucci, P. et al. Impacts of rising temperatures and farm management practices on global yields of 18 crops. Nat. Food 1, 562–571 (2020).Article 

Google Scholar 
Zhu, P. & Burney, J. Temperature‐driven harvest decisions amplify US winter wheat loss under climate warming. Glob. Change Biol. 27, 550–562 (2021).CAS 
Article 

Google Scholar 
Ortiz-Bobea, A., Knippenberg, E. & Chambers, R. G. Growing climatic sensitivity of U.S. agriculture linked to technological change and regional specialization. Sci. Adv. 4, 4343 (2018).Article 

Google Scholar 
Duku, C., Zwart, S. J. & Hein, L. Impacts of climate change on cropping patterns in a tropical, sub-humid watershed. PLoS ONE 13, 0192642 (2018).Article 

Google Scholar 
Folberth, C. et al. The global cropland-sparing potential of high-yield farming. Nat. Sustain. 3, 281–289 (2020).Article 

Google Scholar 
Lobell, D. B. et al. The critical role of extreme heat for maize production in the United States. Nat. Clim. Change 3, 497–501 (2013).Article 

Google Scholar 
Yang, X. et al. Potential benefits of climate change for crop productivity in China. Agric. For. Meteorol. 208, 76–84 (2015).Article 

Google Scholar 
Burney, J., Woltering, L. & Burke, M. Solar-powered drip irrigation enhances food security in the Sudano–Sahel. Proc. Natl Acad. Sci. USA 107, 1848–1853 (2010).CAS 
Article 

Google Scholar 
You, L. et al. What is the irrigation potential for Africa? A combined biophysical and socioeconomic approach. Food Policy 36, 770–782 (2011).Article 

Google Scholar 
Zheng, B., Chenu, K., Fernanda Dreccer, M. & Chapman, S. C. Breeding for the future: what are the potential impacts of future frost and heat events on sowing and flowering time requirements for Australian bread wheat (Triticum aestivium) varieties? Glob. Change Biol. 18, 2899–2914 (2012).Article 

Google Scholar 
Flach, R., Fader, M., Folberth, C., Skalský, R. & Jantke, K. The effects of cropping intensity and cropland expansion of Brazilian soybean production on green water flows. Environ. Res. Commun. 2, 071001 (2020).Article 

Google Scholar 
Wood, S. A., Jina, A. S., Jain, M., Kristjanson, P. & DeFries, R. S. Smallholder farmer cropping decisions related to climate variability across multiple regions. Glob. Environ. Change 25, 163–172 (2014).Article 

Google Scholar 
Paola, A. D. et al. The expansion of wheat thermal suitability of Russia in response to climate change. Land Use Policy 78, 70–77 (2018).Article 

Google Scholar 
Brunelle, T. & Makowski, D. Assessing whether the best land is cultivated first: a quantile analysis. PLoS ONE 15, e0242222 (2020).CAS 
Article 

Google Scholar 
Lark, T. J., Spawn, S. A., Bougie, M. & Gibbs, H. K. Cropland expansion in the United States produces marginal yields at high costs to wildlife. Nat. Commun. 11, 4295 (2020).CAS 
Article 

Google Scholar 
Zabel, F., Putzenlechner, B. & Mauser, W. Global agricultural land resources—a high resolution suitability evaluation and its perspectives until 2100 under climate change conditions. PLoS ONE 9, e107522 (2014).Article 

Google Scholar 
Petkeviciene, B. The effects of climate factors on sugar beet early sowing timing. Agron. Res. 7, 436–443 (2009).
Google Scholar 
Ainsworth, E. A. & Long, S. P. 30 years of free-air carbon dioxide enrichment (FACE): what have we learned about future crop productivity and its potential for adaptation? Glob. Change Biol. 27, 27–49 (2021).CAS 
Article 

Google Scholar 
Collins, M. et al. in Climate Change 2013: The Physical Science Basis (eds Stocker, T. F. et al.) Ch. 11 (Cambridge Univ. Press, 2013).Pendergrass, A. G., Knutti, R., Lehner, F., Deser, C. & Sanderson, B. M. Precipitation variability increases in a warmer climate. Sci. Rep. 7, 17966 (2017).Article 

Google Scholar 
Asadieh, B. & Krakauer, N. Y. Global trends in extreme precipitation: climate models versus observations. Hydrol. Earth Syst. Sci. 19, 877–891 (2015).Article 

Google Scholar 
Zhang, Y., You, L., Lee, D. & Block, P. Integrating climate prediction and regionalization into an agro-economic model to guide agricultural planning. Clim. Change 158, 435–451 (2020).Article 

Google Scholar 
Turner, S. W. D., Hejazi, M., Yonkofski, C., Kim, S. H. & Kyle, P. Influence of groundwater extraction costs and resource depletion limits on simulated global nonrenewable water withdrawals over the twenty‐first century. Earths Future 7, 123–135 (2019).Article 

Google Scholar 
Zhu, W., Jia, S., Devineni, N., Lv, A. & Lall, U. Evaluating China’s water security for food production: the role of rainfall and irrigation. Geophys. Res. Lett. 46, 11155–11166 (2019).Article 

Google Scholar 
FAOSTAT (Food and Agriculture Organization of the United Nations, 1997).Egli, L., Schröter, M., Scherber, C., Tscharntke, T. & Seppelt, R. Crop asynchrony stabilizes food production. Nature 588, E7–E12 (2020).CAS 
Article 

Google Scholar 
Hersbach, H. et al. ERA5 Hourly Data on Single Levels from 1979 to Present (Copernicus Climate Change Service (C3S) Climate Data Store (CDS), accessed 1 August 2020); https://doi.org/10.24381/cds.adbb2d47 (2018).Feng, P. et al. Impacts of rainfall extremes on wheat yield in semi-arid cropping systems in eastern Australia. Clim. Change 147, 555–569 (2018).Article 

Google Scholar 
Teluguntla, P. et al. in Land Resources Monitoring, Modeling, and Mapping with Remote Sensing (ed. Thenkabail, P. S.) 849 (CRC Press, 2015).Hawkins, E. et al. Increasing influence of heat stress on French maize yields from the 1960s to the 2030s. Glob. Change Biol. 19, 937–947 (2013).Article 

Google Scholar 
Friedman, J., Hastie, T. & Tibshirani, R. Regularization paths for generalized linear models via coordinate descent. J. Stat. Softw. 33, 1–22 (2010).Article 

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

Google Scholar 
Deryng, D., Sacks, W. J., Barford, C. C. & Ramankutty, N. Simulating the effects of climate and agricultural management practices on global crop yield. Glob. Biogeochem. Cycles 25, GB2006 (2011).New, M., New, M., Lister, D., Hulme, M. & Makin, I. A high-resolution data set of surface climate over global land areas. Clim. Res. 21, 1–25 (2002).Article 

Google Scholar 
Willmott, C. J. Terrestrial Air Temperature and Precipitation: Monthly and Annual Time Series (1950–1996) (Center for Climatic Research, 2000); http://climate.geog.udel.edu/~climate/html_pages/README.ghcn_ts2.htmlVan Beveren, I. Total factor productivity estimation: a practical review. J. Econ. Surv. 26, 98–128 (2012).Article 

Google Scholar 
Xu, J. et al. Double cropping and cropland expansion boost grain production in Brazil. Nat. Food 2, 264–273 (2021).Article 

Google Scholar 
Friedl, M. & Gray, J. MCD12Q2 MODIS/Terra+ Aqua Land Cover Dynamics Yearly L3 Global 500 m SIN Grid V006 (NASA EOSDIS, 2019).Sulla-Menashe, D. & Friedl, M. A. User Guide to Collection 6 MODIS Land Cover (MCD12Q1 and MCD12C1) Product (USGS, 2018).Schwalm, C. R., Glendon, S. & Duffy, P. B. RCP8.5 tracks cumulative CO2 emissions. Proc. Natl Acad. Sci. USA 117, 19656–19657 (2020).CAS 
Article 

Google Scholar 
Lange, S. Trend-preserving bias adjustment and statistical downscaling with ISIMIP3BASD (v1.0). Geosci. Model Dev. 12, 3055–3070 (2019).Article 

Google Scholar 
Peng Zhu. Climate effects on caloric yield and cropping frequency. Zenodo https://doi.org/10.5281/zenodo.7038556 (2022). More

Knowlton, N. et al. in Life in the World’s Oceans 65–78 (Wiley-Blackwell, 2010).Fisher, R. et al. Species richness on coral reefs and the pursuit of convergent global estimates. Curr. Biol. 25, 500–505. https://doi.org/10.1016/j.cub.2014.12.022 (2015).CAS 
Article 
PubMed 

Google Scholar 
Brandl, S. J., Goatley, C. H. R., Bellwood, D. R. & Tornabene, L. The hidden half: Ecology and evolution of cryptobenthic fishes on coral reefs. Biol. Rev. 93, 1846–1873. https://doi.org/10.1111/brv.12423 (2018).Article 
PubMed 

Google Scholar 
Appeltans, W. et al. The magnitude of global marine species diversity. Curr. Biol. 22, 2189–2202. https://doi.org/10.1016/j.cub.2012.09.036 (2012).CAS 
Article 
PubMed 

Google Scholar 
Carvalho, S. et al. Beyond the visual: Using metabarcoding to characterize the hidden reef cryptobiome. Proc. R. Soc. B Biol. Sci. https://doi.org/10.1098/rspb.2018.2697 (2019).Article 

Google Scholar 
Kramer, M. J., Bellwood, O., Fulton, C. J. & Bellwood, D. R. Refining the invertivore: Diversity and specialisation in fish predation on coral reef crustaceans. Mar. Biol. 162, 1779–1786. https://doi.org/10.1007/s00227-015-2710-0 (2015).CAS 
Article 

Google Scholar 
Brandl, S. J. et al. Demographic dynamics of the smallest marine vertebrates fuel coral reef ecosystem functioning. Science 364, 1189–1192. https://doi.org/10.1126/science.aav3384 (2019).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Kramer, M. J., Bellwood, D. R. & Bellwood, O. Cryptofauna of the epilithic algal matrix on an inshore coral reef, Great Barrier Reef. Coral Reefs 31, 1007–1015. https://doi.org/10.1007/s00338-012-0924-x (2012).ADS 
Article 

Google Scholar 
Rocha, L. A. et al. Specimen collection: An essential tool. Science 344, 814–815. https://doi.org/10.1126/science.344.6186.814 (2014).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Berumen, M. L. et al. The status of coral reef ecology research in the Red Sea. Coral Reefs 32, 737–748. https://doi.org/10.1007/s00338-013-1055-8 (2013).ADS 
Article 

Google Scholar 
Paknia, O., Sh, H. R. & Koch, A. Lack of well-maintained natural history collections and taxonomists in megadiverse developing countries hampers global biodiversity exploration. Org. Divers. Evol. 15, 619–629. https://doi.org/10.1007/s13127-015-0202-1 (2015).Article 

Google Scholar 
Knowlton, N. & Leray, M. Censusing marine life in the twentyfirst Century. Genome 58, 238–238 (2015).
Google Scholar 
Yu, D. W. et al. Biodiversity soup: Metabarcoding of arthropods for rapid biodiversity assessment and biomonitoring. Methods Ecol. Evol. 3, 613–623. https://doi.org/10.1111/j.2041-210X.2012.00198.x (2012).Article 

Google Scholar 
Ransome, E. et al. The importance of standardization for biodiversity comparisons: A case study using autonomous reef monitoring structures (ARMS) and metabarcoding to measure cryptic diversity on Mo’orea coral reefs, French Polynesia. PLoS ONE https://doi.org/10.1371/journal.pone.0175066 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Coker, D. J., DiBattista, J. D., Sinclair-Taylor, T. H. & Berumen, M. L. Spatial patterns of cryptobenthic coral-reef fishes in the Red Sea. Coral Reefs 37, 193–199. https://doi.org/10.1007/s00338-017-1647-9 (2018).ADS 
Article 

Google Scholar 
Pearman, J. K. et al. Cross-shelf investigation of coral reef cryptic benthic organisms reveals diversity patterns of the hidden majority. Sci. Rep. 8, 8090. https://doi.org/10.1038/s41598-018-26332-5 (2018).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Pearman, J. K. et al. Disentangling the complex microbial community of coral reefs using standardized Autonomous Reef Monitoring Structures (ARMS). Mol. Ecol. 28, 3496–3507. https://doi.org/10.1111/mec.15167 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Selkoe, K. A. et al. The DNA of coral reef biodiversity: Predicting and protecting genetic diversity of reef assemblages. Proc. R. Soc. B-Biol. Sci. https://doi.org/10.1098/rspb.2016.0354 (2016).Article 

Google Scholar 
DiBattista, J. D. et al. Digging for DNA at depth: Rapid universal metabarcoding surveys (RUMS) as a tool to detect coral reef biodiversity across a depth gradient. PeerJ https://doi.org/10.7717/peerj.6379 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
DiBattista, J. D. et al. Assessing the utility of eDNA as a tool to survey reef-fish communities in the Red Sea. Coral Reefs 36, 1245–1252. https://doi.org/10.1007/s00338-017-1618-1 (2017).ADS 
Article 

Google Scholar 
Nester, G. M. et al. Development and evaluation of fish eDNA metabarcoding assays facilitate the detection of cryptic seahorse taxa (family: Syngnathidae). Environ. DNA 2, 614–626 (2020).Article 

Google Scholar 
West, K. M. et al. eDNA metabarcoding survey reveals fine-scale coral reef community variation across a remote, tropical island ecosystem. Mol. Ecol. 29, 1069–1086. https://doi.org/10.1111/mec.15382 (2020).CAS 
Article 
PubMed 

Google Scholar 
DiBattista, J. D. et al. Environmental DNA can act as a biodiversity barometer of anthropogenic pressures in coastal ecosystems. Sci. Rep. https://doi.org/10.1038/s41598-020-64858-9 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Worm, B. et al. Impacts of biodiversity loss on ocean ecosystem services. Science 314, 787–790. https://doi.org/10.1126/science.1132294 (2006).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Spalding, M. et al. Mapping the global value and distribution of coral reef tourism. Mar. Policy 82, 104–113. https://doi.org/10.1016/j.marpol.2017.05.014 (2017).Article 

Google Scholar 
Thomsen, P. F. & Willerslev, E. Environmental DNA – An emerging tool in conservation for monitoring past and present biodiversity. Biol. Cons. 183, 4–18. https://doi.org/10.1016/j.biocon.2014.11.019 (2015).Article 

Google Scholar 
Hughes, T. P. et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science 359, 80–83. https://doi.org/10.1126/science.aan8048 (2018).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Monroe, A. A. et al. In situ observations of coral bleaching in the central Saudi Arabian Red Sea during the 2015/2016 global coral bleaching event. PLoS ONE https://doi.org/10.1371/journal.pone.0195814 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Roth, F. et al. Coral reef degradation affects the potential for reef recovery after disturbance. Mar. Environ. Res. 142, 48–58. https://doi.org/10.1016/j.marenvres.2018.09.022 (2018).CAS 
Article 
PubMed 

Google Scholar 
Foster, T. & Gilmour, J. P. Seeing red: Coral larvae are attracted to healthy-looking reefs. Mar. Ecol. Prog. Ser. 559, 65–71. https://doi.org/10.3354/meps11902 (2016).ADS 
CAS 
Article 

Google Scholar 
Karcher, D. B. et al. Nitrogen eutrophication particularly promotes turf algae in coral reefs of the central Red Sea. PeerJ https://doi.org/10.7717/peerj.8737 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Pancrazi, I., Ahmed, H., Cerrano, C. & Montefalcone, M. Synergic effect of global thermal anomalies and local dredging activities on coral reefs of the Maldives. Marine Pollut. Bull. https://doi.org/10.1016/j.marpolbul.2020.111585 (2020).Article 

Google Scholar 
Vercelloni, J. et al. Forecasting intensifying disturbance effects on coral reefs. Glob. Change Biol. 26, 2785–2797. https://doi.org/10.1111/gcb.15059 (2020).ADS 
Article 

Google Scholar 
González-Barrios, F. J., Cabral-Tena, R. A. & Alvarez-Filip, L. Recovery disparity between coral cover and the physical functionality of reefs with impaired coral assemblages. Glob. Change Biol. 27, 640–651. https://doi.org/10.1111/gcb.15431 (2020).ADS 
Article 

Google Scholar 
Rice, M. M., Ezzat, L. & Burkepile, D. E. Corallivory in the anthropocene: Interactive effects of anthropogenic stressors and corallivory on coral reefs. Front. Marine Sci. https://doi.org/10.3389/fmars.2018.00525 (2019).Article 

Google Scholar 
Lin, Y.-J. et al. Long-term ecological changes in fishes and macro-invertebrates in the world’s warmest coral reefs. Sci. Total Environ. 750, 142254. https://doi.org/10.1016/j.scitotenv.2020.142254 (2021).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Loreau, M. & de Mazancourt, C. Biodiversity and ecosystem stability: A synthesis of underlying mechanisms. Ecol. Lett. 16, 106–115. https://doi.org/10.1111/ele.12073 (2013).Article 
PubMed 

Google Scholar 
Cardinale, B. J. et al. Biodiversity loss and its impact on humanity. Nature 486, 59–67. https://doi.org/10.1038/nature11148 (2012).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Handley, L. L. How will the “molecular revolution’ contribute to biological recording?. Biol. J. Lin. Soc. 115, 750–766. https://doi.org/10.1111/bij.12516 (2015).Article 

Google Scholar 
Ducklow, H. W., Doney, S. C. & Steinberg, D. K. Contributions of long-term research and time-series observations to marine ecology and biogeochemistry. Ann. Rev. Mar. Sci. 1, 279–302. https://doi.org/10.1146/annurev.marine.010908.163801 (2009).Article 
PubMed 

Google Scholar 
Hughes, B. B. et al. Long-term studies contribute disproportionately to ecology and policy. Bioscience 67, 271–281. https://doi.org/10.1093/biosci/biw185 (2017).Article 

Google Scholar 
Kraft, N. J. B. et al. Community assembly, coexistence and the environmental filtering metaphor. Funct. Ecol. 29, 592–599. https://doi.org/10.1111/1365-2435.12345 (2015).Article 

Google Scholar 
Leibold, M. A. et al. The metacommunity concept: A framework for multi-scale community ecology. Ecol. Lett. 7, 601–613. https://doi.org/10.1111/j.1461-0248.2004.00608.x (2004).Article 

Google Scholar 
Vellend, M. The Theory of Ecological Communities (MPB-57). (Princeton University Press, 2016).Condon, R. H. et al. Recurrent jellyfish blooms are a consequence of global oscillations. Proc. Natl. Acad. Sci. U.S.A. 110, 1000–1005. https://doi.org/10.1073/pnas.1210920110 (2013).ADS 
Article 
PubMed 

Google Scholar 
Boero, F., Kraberg, A. C., Krause, G. & Wiltshire, K. H. Time is an affliction: Why ecology cannot be as predictive as physics and why it needs time series. J. Sea Res. 101, 12–18. https://doi.org/10.1016/j.seares.2014.07.008 (2015).ADS 
Article 

Google Scholar 
Pearman, J. K., Anlauf, H., Irigoien, X. & Carvalho, S. Please mind the gap – Visual census and cryptic biodiversity assessment at central Red Sea coral reefs. Mar. Environ. Res. 118, 20–30. https://doi.org/10.1016/j.marenvres.2016.04.011 (2016).CAS 
Article 
PubMed 

Google Scholar 
David, R. et al. Lessons from photo analyses of autonomous reef monitoring structures as tools to detect (bio-)geographical, spatial, and environmental effects. Mar. Pollut. Bull. 141, 420–429. https://doi.org/10.1016/j.marpolbul.2019.02.066 (2019).CAS 
Article 
PubMed 

Google Scholar 
Pennesi, C. & Danovaro, R. Assessing marine environmental status through microphytobenthos assemblages colonizing the autonomous reef monitoring structures (ARMS) and their potential in coastal marine restoration. Mar. Pollut. Bull. 125, 56–65. https://doi.org/10.1016/j.marpolbul.2017.08.001 (2017).CAS 
Article 
PubMed 

Google Scholar 
Chang, J. J. M., Ip, Y. C. A., Bauman, A. G. & Huang, D. MinION-in-ARMS: Nanopore sequencing to expedite barcoding of specimen-rich macrofaunal samples from Autonomous Reef Monitoring Structures. Front. Marine Sci. https://doi.org/10.3389/fmars.2020.00448 (2020).Article 

Google Scholar 
Hazeri, G. et al. Latitudinal species diversity and density of cryptic crustacean (Brachyura and Anomura) in micro-habitat Autonomous Reef Monitoring Structures across Kepulauan Seribu, Indonesia. Biodivers. J. Biol. Divers. 20 (2019).Al-Rshaidat, M. M. D. et al. Deep COI sequencing of standardized benthic samples unveils overlooked diversity of Jordanian coral reefs in the northern Red Sea. Genome 59, 724–737. https://doi.org/10.1139/gen-2015-0208 (2016).CAS 
Article 
PubMed 

Google Scholar 
Pearman, J. K. et al. Pan-regional marine benthic cryptobiome biodiversity patterns revealed by metabarcoding Autonomous Reef Monitoring Structures. Mol. Ecol. https://doi.org/10.1111/mec.15692 (2020).Article 
PubMed 

Google Scholar 
Leray, M. & Knowlton, N. DNA barcoding and metabarcoding of standardized samples reveal patterns of marine benthic diversity. Proc. Natl. Acad. Sci. U.S.A. 112, 2076–2081. https://doi.org/10.1073/pnas.1424997112 (2015).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Obst, M. et al. A marine biodiversity observation network for genetic monitoring of hard-bottom communities (ARMS-MBON). Front. Marine Sci. https://doi.org/10.3389/fmars.2020.572680 (2020).Article 

Google Scholar 
Hughes, T. P. et al. Ecological memory modifies the cumulative impact of recurrent climate extremes. Nat. Clim. Chang. 9, 40–43. https://doi.org/10.1038/s41558-018-0351-2 (2019).ADS 
Article 

Google Scholar 
Hughes, T. P., Kerry, J. T. & Simpson, T. Large-scale bleaching of corals on the Great Barrier Reef. Ecology 99, 501–501. https://doi.org/10.1002/ecy.2092 (2018).CAS 
Article 
PubMed 

Google Scholar 
Furby, K. A., Bouwmeester, J. & Berumen, M. L. Susceptibility of central Red Sea corals during a major bleaching event. Coral Reefs 32, 505–513. https://doi.org/10.1007/s00338-012-0998-5 (2013).ADS 
Article 

Google Scholar 
Froehlich, C. Y. M., Klanten, O. S., Hing, M. L., Dowton, M. & Wong, M. Y. L. Uneven declines between corals and cryptobenthic fish symbionts from multiple disturbances. Sci. Rep. https://doi.org/10.1038/s41598-021-95778-x (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Bellwood, D. R. et al. Coral recovery may not herald the return of fishes on damaged coral reefs. Oecologia 170, 567–573. https://doi.org/10.1007/s00442-012-2306-z (2012).ADS 
Article 
PubMed 

Google Scholar 
Archana, A. & Baker, D. M. Multifunctionality of an urbanized coastal marine ecosystem. Front. Marine Sci. https://doi.org/10.3389/fmars.2020.557145 (2020).Article 

Google Scholar 
Servis, J. A., Reid, B. N., Timmers, M. A., Stergioula, V. & Naro-Maciel, E. Characterizing coral reef biodiversity: Genetic species delimitation in brachyuran crabs of Palmyra Atoll Central Pacific. Mitochondrial DNA Part A 31, 178–189. https://doi.org/10.1080/24701394.2020.1769087 (2020).CAS 
Article 

Google Scholar 
Chaves-Fonnegra, A. et al. Bleaching events regulate shifts from corals to excavating sponges in algae-dominated reefs. Glob. Change Biol. 24, 773–785. https://doi.org/10.1111/gcb.13962 (2018).ADS 
Article 

Google Scholar 
Perry, C. T. & Morgan, K. M. Post-bleaching coral community change on southern Maldivian reefs: Is there potential for rapid recovery?. Coral Reefs 36, 1189–1194. https://doi.org/10.1007/s00338-017-1610-9 (2017).ADS 
Article 

Google Scholar 
DeCarlo, T. M. The past century of coral bleaching in the Saudi Arabian central Red Sea. PeerJ https://doi.org/10.7717/peerj.10200 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Cortés, J. et al. in Coral Reefs of the Eastern Tropical Pacific: Persistence and Loss in a Dynamic Environment (eds Peter W. Glynn, Derek P. Manzello, & Ian C. Enochs) 203–250 (Springer Netherlands, 2017).Enochs, I. C. & Manzello, D. P. Species richness of motile cryptofauna across a gradient of reef framework erosion. Coral Reefs 31, 653–661. https://doi.org/10.1007/s00338-012-0886-z (2012).ADS 
Article 

Google Scholar 
Timmers, M. A. et al. Biodiversity of coral reef cryptobiota shuffles but does not decline under the combined stressors of ocean warming and acidification. Proc. Natl. Acad. Sci. 118, e2103275118. https://doi.org/10.1073/pnas.2103275118 (2021).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Khalil, M. T., Bouwmeester, J. & Berumen, M. L. Spatial variation in coral reef fish and benthic communities in the central Saudi Arabian Red Sea. PeerJ https://doi.org/10.7717/peerj.3410 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Roik, A. et al. Year-long monitoring of physico-chemical and biological variables provide a comparative baseline of coral reef functioning in the central Red Sea. PLoS ONE https://doi.org/10.1371/journal.pone.0163939 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Largier, J. L. Considerations in estimating larval dispersal distances from oceanographic data. Ecol. Appl. 13, S71–S89 (2003).Article 

Google Scholar 
Volkov, I., Banavar, J. R., Hubbell, S. P. & Maritan, A. Patterns of relative species abundance in rainforests and coral reefs. Nature 450, 45–49. https://doi.org/10.1038/nature06197 (2007).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Alsaffar, Z., Cúrdia, J., Borja, A., Irigoien, X. & Carvalho, S. Consistent variability in beta-diversity patterns contrasts with changes in alpha-diversity along an onshore to offshore environmental gradient: The case of Red Sea soft-bottom macrobenthos. Mar. Biodivers. 49, 247–262. https://doi.org/10.1007/s12526-017-0791-3 (2017).Article 

Google Scholar 
Alsaffar, Z. et al. The role of seagrass vegetation and local environmental conditions in shaping benthic bacterial and macroinvertebrate communities in a tropical coastal lagoon. Sci. Rep. https://doi.org/10.1038/s41598-020-70318-1 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Rocha, L. A. et al. Mesophotic coral ecosystems are threatened and ecologically distinct from shallow water reefs. Science 361, 281–284. https://doi.org/10.1126/science.aaq1614 (2018).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Soininen, J., Lennon, J. J. & Hillebrand, H. A multivariate analysis of beta diversity across organisms and environments. Ecology 88, 2830–2838. https://doi.org/10.1890/06-1730.1 (2007).Article 
PubMed 

Google Scholar 
Chust, G. et al. Dispersal similarly shapes both population genetics and community patterns in the marine realm. Sci. Rep. https://doi.org/10.1038/srep28730 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Gianuca, A. T., Declerck, S. A. J., Lemmens, P. & De Meester, L. Effects of dispersal and environmental heterogeneity on the replacement and nestedness components of beta-diversity. Ecology 98, 525–533. https://doi.org/10.1002/ecy.1666 (2017).Article 
PubMed 

Google Scholar 
Enochs, I. C., Toth, L. T., Brandtneris, V. W., Afflerbach, J. C. & Manzello, D. P. Environmental determinants of motile cryptofauna on an eastern Pacific coral reef. Mar. Ecol. Prog. Ser. 438, 105-U127. https://doi.org/10.3354/meps09259 (2011).ADS 
Article 

Google Scholar 
Hughes, T. P. et al. Coral reefs in the anthropocene. Nature 546, 82–90. https://doi.org/10.1038/nature22901 (2017).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Fabricius, K. E. Effects of terrestrial runoff on the ecology of corals and coral reefs: Review and synthesis. Mar. Pollut. Bull. 50, 125–146. https://doi.org/10.1016/j.marpolbul.2004.11.028 (2005).CAS 
Article 
PubMed 

Google Scholar 
Chaidez, V., Dreano, D., Agusti, S., Duarte, C. M. & Hoteit, I. Decadal trends in Red Sea maximum surface temperature. Sci. Rep. https://doi.org/10.1038/s41598-018-25731-y (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Hubbell, S. P. in Monographs in Population Biology. The unified neutral theory of biodiversity and biogeography Vol. 32 Monographs in Population Biology i-xiv, 1–375 (2001).Dornelas, M., Connolly, S. R. & Hughes, T. P. Coral reef diversity refutes the neutral theory of biodiversity. Nature 440, 80–82. https://doi.org/10.1038/nature04534 (2006).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Baselga, A. Partitioning the turnover and nestedness components of beta diversity. Glob. Ecol. Biogeogr. 19, 134–143. https://doi.org/10.1111/j.1466-8238.2009.00490.x (2010).Article 

Google Scholar 
Legendre, P. Interpreting the replacement and richness difference components of beta diversity. Glob. Ecol. Biogeogr. 23, 1324–1334. https://doi.org/10.1111/geb.12207 (2014).Article 

Google Scholar 
Hollander, M. & Wolfe, D. A. Nonparametric statistical methods. Ergonomics 18, 701–702 (1975).
Google Scholar 
Kohler, K. E. & Gill, S. M. Coral point count with excel extensions (CPCe): A visual basic program for the determination of coral and substrate coverage using random point count methodology. Comput. Geosci. 32, 1259–1269. https://doi.org/10.1016/j.cageo.2005.11.009 (2006).ADS 
Article 

Google Scholar 
Geller, J., Meyer, C., Parker, M. & Hawk, H. Redesign of PCR primers for mitochondrial cytochrome c oxidase subunit I for marine invertebrates and application in all-taxa biotic surveys. Mol. Ecol. Resour. 13, 851–861. https://doi.org/10.1111/1755-0998.12138 (2013).CAS 
Article 
PubMed 

Google Scholar 
Hao, X., Jiang, R. & Chen, T. Clustering 16S rRNA for OTU prediction: A method of unsupervised Bayesian clustering. Bioinformatics 27, 611–618. https://doi.org/10.1093/bioinformatics/btq725 (2011).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581. https://doi.org/10.1038/nmeth.3869 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2020).Ranwez, V., Harispe, S., Delsuc, F. & Douzery, E. J. P. MACSE: Multiple alignment of coding SEquences accounting for frameshifts and stop codons. PLoS ONE https://doi.org/10.1371/journal.pone.0022594 (2011).Article 
PubMed 
PubMed Central 

Google Scholar 
Machida, R. J., Leray, M., Ho, S. L. & Knowlton, N. Data Descriptor: Metazoan mitochondrial gene sequence reference datasets for taxonomic assignment of environmental samples. Sci. Data https://doi.org/10.1038/sdata.2017.27 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. https://doi.org/10.1128/aem.00062-07 (2007).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Generate High-Resolution Venn and Euler Plots v. 1.6.20 (2018).Ginestet, C. ggplot2: Elegant graphics for data analysis. J. R. Stat. Soc. Ser. Stat. Soc. 174, 245–245. https://doi.org/10.1111/j.1467-985X.2010.00676_9.x (2011).Article 

Google Scholar 
McMurdie, P. J. & Holmes, S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE https://doi.org/10.1371/journal.pone.0061217 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Anderson, M. J. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46. https://doi.org/10.1111/j.1442-9993.2001.01070.pp.x (2001).Article 

Google Scholar 
Hervé, M. Testing and plotting procedures for biostatistics v. 0.9-79. Retrieved from https://cran.r-project.org/web/packages/RVAideMemoire/index.html (2021).De Caceres, M. & Legendre, P. Associations between species and groups of sites: Indices and statistical inference. Ecology 90, 3566–3574. https://doi.org/10.1890/08-1823.1 (2009).Article 
PubMed 

Google Scholar 
Legendre, P. & Anderson, M. J. Distance-based redundancy analysis: Testing multispecies responses in multifactorial ecological experiments. Ecol. Monogr. 69, 1–24. https://doi.org/10.1890/0012-9615(1999)069[0001:dbratm]2.0.co;2 (1999).Article 

Google Scholar 
Roberts, D. Ordination and multivariate analysis for ecology v. 2.0-1. Retrieved from http://ecology.msu.montana.edu/labdsv/R (2019).Dray, S., Bauman, D., Blanchet, G., Borcard, D., Clappe, S., Guenard, G. & Wagner, H. Adespatial: Multivariate multiscale spatial analysis v. 0.3-13. Retrieved from https://cran.r-project.org/package=adespatial (2021). More

Overview of U.S. domestic migration flowsThe most recent ACS county-to-county flow dataset26 reports that about 45.6 million people migrated to the U.S. per year during the period 2015–2019, which corresponds to 14.2% of the U.S. population27. Approximately 43.5 million annual moves corresponded to domestic migration (moves within the U.S.28), while 2.1 million accounted for inflows of individuals from other countries (viz. international immigration).With respect to domestic migration, 25.7 million people per year migrated within the same county, thus showing that the highest share of domestic flows (59%) is intra-county. Annually, about 10.4 people moved between different counties within the same state, thus intra-state flows account for 24% of the domestic migration (Supplementary Fig. 1), mainly driven by the search for more affordable housing, better jobs, and for family reasons such as change in marital status29. Long distance moves, captured by inter-state flows, represent the remaining 17% of domestic flows, which comprises about 7.5 million moves per year. Here, we will refer to these domestic migration flows as inflows or outflows, and netflows (inflows-outflows).The United States Office of Management and Budget (OMB) classifies counties as metropolitan, micropolitan, or neither30. A metropolitan statistical area contains a core urban area of at least 50,000 population. A metro area represents a functional delineation of an urban area with a network of strong socioeconomic ties, and provision of infrastructure services31,32,33. A micropolitan statistical area contains an urban core of at lest 10,000 but less than 50,000 inhabitants. There are over 380 metropolitan statistical areas in the U.S., each composed of one or more counties, accounting for about 86% of the total U.S. population and comprising approximately 28% of the land area of the country. For this reason, our analysis focuses on the growth dynamics of MSA counties. Supplementary Fig. 2 shows the 3141 counties (administrative subdivisions of the states) in the U.S., comprising about 321 million inhabitants in the starting year of the ACS 5-Year survey period (2015–2019) of our analysis26.Population growth has two components, namely natural growth and migration. Natural growth accounts for births minus deaths, and migration comprises domestic and international migration. With recent trends showing that births and natural increase have declined in the U.S. and in recent years contribute less to overall city population growth34,35, migration patterns become more relevant to the study of city population growth. Because the ACS flow files contain international inflows only, the relative importance of migrations on population growth is here addressed by x = ∣Inflows−Outflows∣/∣Births−Deaths∣ (Supplementary Figs. 3, 4), which is the ratio between domestic netflows and natural growth. The statistical distribution of this quantity computed for all U.S. counties is well fitted by a lognormal distribution, and shows that x≥1 for 76.5% of counties. For most counties, domestic migration dominates population growth, and understanding the spatial structure of domestic netflows (and their distribution within a city) is crucial to the comprehension of the mechanisms behind the heterogeneity of city population growth.At this spatial granularity, we observe a strong heterogeneity among the U.S. counties (Supplementary Fig. 2) for the period 2015 − 2019, along with examples of specific MSAs. In particular, the relative dispersion of counties relative growth due to netflows is higher than one for about 85% of the metro areas, indicating a large heterogeneity within the same city and pointing towards the spatial structure of domestic migration. The observed difference in the netflows stresses the relevance of our approach: counties belonging to the same city may have specific growth rates due to population flow patterns, thus indicating preferential flow destinations and pinpointing the direction in which the city has expanded.Heterogeneity of inter- and intra-city flowsInter-city flows represent the major component of the total flows (~55%), while intra-city flows represent ~25%. Flows between metro and micro areas, and between metro and non-statistical areas are the smallest components, with ~13% and ~7%, respectively. Given that about 80% of the domestic migration are composed of intra- and inter-city flows, we will focus our attention on describing the structure of intra- and inter-city flows, but in the Supplementary Information we offer a brief analysis of flows between metro and micro areas, and between metro and non-statistical areas.Inter-city flows are not uniform across the U.S. cities. The most intense annual netflows ( >2000 people per year), accounting for approximately 17% of the entire inter-city U.S. netflows, are mainly from New York and Chicago to California and Florida (Fig. 2), and from Los Angeles to neighboring cities. Notably, netflows among the Midwestern cities are mostly negative and below the threshold we set. These flows are mainly responsible for increasing or decreasing the population of a given city. Intra-city flow patterns, illustrated with the 7 most populous U.S. cities with more than 5 counties, are also non-uniform.Fig. 2: Heterogeneity of inter- and intra-city netflows.The map (A) suggests that the domestic redistribution of people between different U.S. metro areas are non-uniform: the black arrows, indicating the direction of the most intense inter-city netflows (higher than 2000 people per year), reveal migration trends from northern and eastern cities to western and southern regions. Cities (composed of one or more counties) are colored according to the relative growth (viz. population growth adjusted by population) of the whole MSA during the 2015–2019 period, and the black intensity and the thickness of the arrows are proportional to the netflows. Alaska and Hawaii are not shown. Panels (B–H), which are close-up of New York (B), Chicago (C), Dallas (D), Houston (E), Washington D.C. (F), Philadelphia (G), Atlanta (H), suggest that the most intense intra-city netflows are oriented radially outwards: people are moving from core to external counties. Here, counties are colored according to their relative growth in the 2015–2019 period and the width of the arrows is proportional to the netflows between origin and destination counties.Full size imageOur analysis reveals that city centers (defined as the core county with the highest population density) are more likely to have negative netflows, indicating that people are leaving the central regions of cities. The arrows in Fig. 2 indicate the direction of the most intense netflows, supporting this finding and highlighting that there is a trend of people moving from internal to external regions, contributing to population growth and spatial expansion of U.S. cities. In fact, we found no correlation between relative population growth (viz. population growth by county size) and distance from the core county (Supplementary Fig. 5A) for the 46 cities with more than 5 counties, with relative growth about 0.03 ± 0.05. On the other hand, we found that relative natural growth (Supplementary Fig. 5B) is negatively correlated with the distance to core county, thus natural growth is less relevant as a component of growth in the outer regions of cities. Consequently, our results show that not only the contribution of each component of growth changes with distance to core county, but also that the internal redistribution of people is an important mechanism of growth, mainly in the external counties.We also examined variability in inter- and intra-city flows within the 50 states (Supplementary Fig. 6). Total flows within a state increase, as expected, with the state population. Two special cases are, however, of interest: (1) two states (Vermont and Rhode Island) with small populations have only one MSA, in which case within-state inter-city flows are zero; and (2) nearly 40%, or 149, of MSAs have only one county, in which case intra-city flows could not be estimated. For all other cases, we observe on average an equal split between inter- and intra-city flows, but with considerable variability among the states, with a mean about 0.5 and standard deviation about 0.2. A generalization of the intra- and inter-city migratory patterns for all 46 cities with more than 5 counties shows that the percentage of migrants from intra- and inter-city flows are of the same order of magnitude (Fig. 3).Fig. 3: Roles of intra- and inter-city flows in driving the heterogeneous population growth of cities.We define the core county as the one with the highest population density, and we plot the percentage of inflows due to intra- (A) and inter-city flows (B) of each county within a city as a function of its distance to the core county. The percentage of outflows due to intra- and inter-city flows are shown in (C) and (D), respectively. The positive correlation of the relative growth with distance due to intra-city flows in (E), along with the lack of correlation due to inter-city flows in (F), indicates that intra-city flows are mainly responsible for increasing the population in the external regions of cities. The sizes of red circles and blue squares are proportional to the city population. The range of distances is split into equally spaced bins. The number of counties n within each bin, from left to right, is 46, 1, 4, 7, 7, 17, 21, 31, 36, 38, 34, 31, 31, 30, 20, 20, 21, 14, 17, 9, 9, 6, 4, 2, 5, 5, 2, 1. The black dots and the error bars indicate the mean and the 90% interval, respectively, of the counties within the corresponding bin. We also show the Pearson correlation coefficient R and the p-value associated with the two-sided test of the null hypothesis of non-correlation.Full size imageApart from the core county, flows from the same city correspond to about 50% of the inflow of people in the counties, presenting a slightly positive correlation with their distance from the city center (Fig. 3A). The low percentage for the core county indicates that it is not the major destination of flows from the same city. The percentage of inflows from other cities is higher in the core county and decays as we move towards the suburbs (Fig. 3B). The moderate negative correlation of this percentage with the distance reveals that inflows from other cities are more likely to concentrate in the core regions of a city.The percentage of outflows directed from the core county to other counties within the same city has a slightly negative correlation with the distance of the origin county to the city center, so it is more likely to find intra-city flows with outflows from internal regions (Fig. 3C). The core county is an exception again, suggesting that it is less likely that someone leaving the core county will move to another county within the same city. The slightly negative correlation of the percentage of outflows directed to other cities suggests that there is a trend of people leaving the core county and the central regions to move to other cities (Fig. 3D). The high percentage of inflows (Fig. 3B) and outflows (Fig. 3D) in the central region due to inter-city flows implies that the central regions of cities are more dynamic and diverse and that people tend to move to counties with similar levels of urbanization. The same pattern is observed for flows between metro and micro areas, and for metro and non-statistical areas, allowing us to conclude that people moving from rural areas are more likely to move to the external regions of a city (Supplementary Fig. 7).The positive correlation of the relative growth with the distance due to intra-city flows (Fig. 3E) shows that the resulting intra-city redistribution of people, given by the difference between inflows and outflows, is such that there is a trend from core county to the external counties (viz. suburbs). When compared to the relative growth due to inter-city flows (Fig. 3F), which do not show any trend and that have negative values for the most distant counties, it becomes clear that intra-city flows play a major role in the population increase observed in outer regions of cities. Interestingly, large circle and square dots in Fig. 3E and F suggest that the loss of people due to inter-city netflows is more intense than the gain of people due to intra-city netflows in some external counties of the largest metro areas, thus explaining the population decline in some outer regions of New York and Chicago (as shown in Fig. 2B and C).The population growth due to intra-city flows is also depicted in Fig. 4. The concentration of flows below the diagonal captures the heterogeneity and the preferential destination of intra-city netflows. We observe that people are more likely to move to lower population density counties when moving from one place to another within the same city, as exemplified by 7 cities in panel A. Panel B summarizes this analysis for the 46 cities with more than 5 counties by showing the fraction ({{{{{{{mathcal{F}}}}}}}}) of intra-city netflows to lower density counties. We note that more than 93% of the cities have ({{{{{{{mathcal{F}}}}}}}} > 0.5) and that there is a positive correlation of ({{{{{{{mathcal{F}}}}}}}}) with the city population, and C shows the rank of cities according to the fraction of intra-city netflows to lower density counties.Fig. 4: People are moving to counties with lower population density.A The population density of the origin (ρo) and destination (ρd) counties of intra-city netflows for New York, Chicago, Dallas, Houston, Washington D.C., Philadelphia, Atlanta, reveal that the majority of the flows occur from high to low-density counties. The size of the symbols are proportional to the intensity of the netflow, and the black line corresponds to y = x. B The fraction of netflows to lower density counties ({{{{{{{mathcal{F}}}}}}}}) has a positive correlation with city population when we consider the 46 MSAs with more than 5 counties, suggesting that intra-city netflows to lower density counties are more frequent as the city size increases. We also show the Pearson correlation coefficient R and the p-value associated with the two-sided test of the null hypothesis of non-correlation. C The ranking of the cities according to ({{{{{{{mathcal{F}}}}}}}}).Full size imagePopulation density does not seem to play a major role in driving flows between counties of different cities. The fraction of inter-city netflows to lower density counties is about 57% when we consider all the 384 MSAs. The heterogeneity in the inter-city netflow pattern can be assessed by analyzing ({{{{{{{mathcal{F}}}}}}}}) versus the population of the destination city (Fig. 5A, B) and ({{{{{{{mathcal{F}}}}}}}}) versus the population of the origin city (Fig. 5C, D). The negative correlation of ({{{{{{{mathcal{F}}}}}}}}) with the population of the destination city in panel A indicates that inflows are more likely to come from lower density counties as the destination city size increases. The positive correlation of ({{{{{{{mathcal{F}}}}}}}}) with the population of the origin city in panel C reveals that outflows tend to be directed to lower density counties as the origin city size increases. The trends observed in panels A and C reveal that inter-city flows are more likely between counties with different population densities rather than between counties with similar population densities. Panels B and D show the rank order of cities according to a function of the destination city size and the origin city size, respectively.Fig. 5: Inter-city flow patterns depend on the population size of the origin and destination cities.Each point corresponds to a particular city. A Fraction ({{{{{{{mathcal{F}}}}}}}}) of netflows going to lower density counties versus the population of the destination city. Inflows to counties of large cities (with population greater than 106, dashed line) usually comes from counties with lower population densities. B Rank of cities according to the share of inflows from lower density counties. C Fraction ({{{{{{{mathcal{F}}}}}}}}) versus the population of the origin city. Outflows from counties of large cities usually go to cities with lower density counties. D The rank of cities according to the share of inter-city netflows to lower density counties is presented. The dots are colored according to the city population density (darker red means higher density). We also show the Pearson correlation coefficient R and the p-value associated with the two-sided test of the null hypothesis of non-correlation.Full size imageWe would expect that there might be preferential locations within a given city to which people move due to various factors such as lower costs of housing and employment opportunities. However, it seems that house prices have little to no effect on intra-city netflows (Supplementary Fig. 8). While the fraction of intra-city netflows to counties with less expensive houses is about 0.8 for cities like New York, Chicago and Washington, this fraction is about 0.2 for cities like Dallas, Houston and Philadelphia. The lack of a clear national pattern highlights the specificity of each city and the heterogeneity of the regional housing market in the U.S.36,37. On the other hand, the fraction of intra-city netflows to counties with lower unemployment rates is higher than 0.5 for about 2/3 of the cities (Supplementary Fig. 9), thus showing that people are more likely to move to counties with lower unemployment rates.Statistical structure of inter-city flowsIntra-city flows capture the internal redistribution of population, without altering the total city population. In this context, we focus on inter-city flows to investigate whether or not extreme flows play an important role in shaping the growth of counties as observed at the city level5. For cities, Verbavatz and Barthelemy5 introduce a stochastic equation to describe population growth, composed of two terms. The first term accounts for out-of-system growth, which includes natural growth and international migration, and the second term accounts for the growth due to domestic netflows. They find that total netflows adjusted by population size can be well approximated by a Lévy distribution, and this heavy-tailed distribution indicates that rare and extreme inter-city flows (viz. migratory shocks) dominate city population growth.Here, we find that, for counties, the distribution of total netflows adjusted by population size, which is represented by ζi and captures the intensity of inter-city migratory flows (see the section “Methods” for details), can be approximated by a Gaussian distribution (Fig. 6). The lack of a heavy tail in the empirical distribution of ζi suggests the absence of extreme flows at the county level, thus indicating that the growth of counties can be described by smoother migratory process than cities. Given that cities do experience migratory shocks5, our findings indicate that cities redistribute inflows among its different counties, leading to a spill-over effect that dampens flow shocks at the county level.Fig. 6: Extreme shocks are dissipated at the county level.The distribution of ζi, which is the sum of the netflows of a county i adjusted by its population, suggests that migratory events are exponentially bounded at the county level since ζi is well described by a Gaussian distribution. The distribution of ζi is computed here for all the counties with at least 50.000 inhabitants. We also show the result of the two-sided KS test under the null hypothesis that ζi follows a Gaussian distribution.Full size imageHeterogeneity of international inflowsThe highest share of international inflows is concentrated in large cities. About 40% of the international inflows are destined to the top 10 (~2.6%) largest metro areas of the U.S. New York is the first with 8.5% of international inflows, followed by Los Angeles and Miami with 5.4% and 5.0%, respectively. Indeed, international inflows Yk scale superlinearly with the population Sk of the metro area k (Fig. 7A), thus larger cities have more immigrants per capita than smaller cities.Fig. 7: International inflow scales superlinearly with city size.Panel (A) shows the number of international immigrants as a function of the city size S for the 384 U.S. metro areas. The performance of the model Y = Y0Sθ, in which θ = 1.19 (95% CI [1.13, 1.24]) and Y0 = 4.10−4 is a normalization constant, is assessed by the coefficient of determination R2. Note that the spread of empirical data around the model narrows as the size of the city increases. Panel (B) shows the rank of the metro areas and the residues, which captures the deviation from the null model thus highligthing cities receiving more/less than expected international inflows. Names of the cities are followed by two-letter state abbreviations.Full size imageInterestingly, this gain with scale is also observed in socioeconomic city metrics as crime, GDP, innovation and wealth creation due to the manifestation of nonlinear agglomeration phenomena38,39,40. Using Y = Y0Sθ as a null model, we can compute deviations from the average behavior by means of residuals given by (log ({Y}_{k}/{Y}_{0}{S}_{k}^{theta }))38. The rank of the residues (Fig. 7B) shows that college towns are among the top metro areas receiving more international inflows than expected, while large cities as Los Angeles, New York, Atlanta, and Chicago are among the metro areas receiving less international inflows than expected.The spatial distribution of international inflows within cities is shown in Supplementary Fig. 10. The highest share of inflows is concentrated at core counties, and the percentage of inflows decreases dramatically with the distance from the core county. This result suggests that inflow of international migrants is an important component of population growth, particularly at the core regions of large cities.Robustness of our findingsPatterns of population redistribution change from time to time in the U.S., and are affected by several factors. For instance, in the 1960s non-metropolitan counties lost about 3 million people due to outflows to metropolitan counties, while the reverse trend was observed in the 1970s when non-metropolitan counties experienced net inflows of about 2.6 million people41. Wardwell and Brown in41 indicate that three factors might be among the main reasons of such change, namely economic decentralization, preference for rural living, and modernization of rural life. The temporal influence of factors as socioeconomic conditions, transportation infrastructure, natural amenities, and land use and development on population growth in rural and suburban areas is explored in42. Changes in rural migration patterns are also studied in43, where age-specific rural migration patterns from 1950 to 1995 are analyzed. In44, the authors explore redistribution trends across U.S. counties from 1980 to 1995 split into three five year periods (1980–1985, 1985–1990, 1990–1995), and45 analyzes changes in age-specific nationwide migration patterns from 1950 to 2010.The spatial structure of migration patterns may indeed change from time to time; our results correspond to the current intra- and inter-city redistribution trends, based on the most recent ACS migration flow files. We present a thorough empirical and statistical analysis of domestic migration flows among U.S. cities ans counties. Our study also introduces a framework that can be used for analyzing and comparing internal redistribution of people across different time periods. Indeed, we extended our analysis for two other time periods, 2005–2009 and 2010–2014. With respect to the spatial distribution of intra- and inter-city flows, similar trends are observed in both periods (Supplementary Figs. 12, 13), namely inter-city flows are responsible for the highest share of inflows to core counties, and intra-city flows are responsible for the highest share of inflows to external counties. We also explored the role of population density in driving netflows between counties within the same metro area in 2005–2009 and 2010–2014. The results (Supplementary Figs. 14 and 15) indicate that 95.7% of cities were dominated by intra-city moves to lower density counties in 2005–2009, and this percentage dropped to 76.1% in 2010–2014. Our findings indicate that the trends we report here are taking place since 2005 but with different intensities.The robustness of our findings is checked with additional migration data from the Internal Revenue Service (IRS), which reports the year-to-year address changes on individual tax returns filled with the IRS46. The results obtained with the analysis of IRS datasets from periods 2015–2016, 2016–2017, 2017–2018, 2018–2019 (Supplementary Figs. 16, 17, 18, 19), reveal similar trends to those we found using ACS data. Particularly, we observe that, for all periods considered, the correlation between intra-city netflow/S and distance to core county is stronger than we found with ACS data, thus highlighting the role of intra-city flows in driving population to external regions of cities. The main difference between both datasets is in the percentage of intra- and inter-city inflows and outflows: while ACS data indicates that both flows have approximately the same contribution to the total flows, the IRS data indicates that, besides the core county, intra-city flows are responsible for about 80% of inflows and outflows of metro areas. More

Abascal F, Zardoya R, Telford MJ (2010) TranslatorX: Multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res 38:7–13Article 

Google Scholar 
Akashi HD, Cádiz Díaz A, Shigenobu S, Makino T, Kawata M (2016) Differentially expressed genes associated with adaptation to different thermal environments in three sympatric Cuban Anolis lizards. Mol Ecol 25:2273–2285CAS 
PubMed 
Article 

Google Scholar 
Allen JA (1877) The influence of Physical conditions in the genesis of species. Radic Rev 1:108–140
Google Scholar 
Barghi N, Hermisson J, Schlötterer C (2020) Polygenic adaptation: A unifying framework to understand positive selection. Nat Rev Genet 21:769–781CAS 
PubMed 
Article 

Google Scholar 
Blem CR (1990) Avian energy storage. Curr Ornithol 7:59–113
Google Scholar 
Blix AS (2016) Adaptations to polar life in mammals and birds. J Exp Biol 219:1093–1105PubMed 
Article 

Google Scholar 
Boyle EA, Li YI, Pritchard JK (2017) An expanded view of complex traits: From polygenic to omnigenic. Cell 169:1177–1186CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Cannon B, Nedergaard J (2010) Thyroid hormones: Igniting brown fat via the brain. Nat Med 16:965–967CAS 
PubMed 
Article 

Google Scholar 
Castruita JAS, Westbury MV, Lorenzen ED (2020) Analyses of key genes involved in Arctic adaptation in polar bears suggest selection on both standing variation and de novo mutations played an important role. BMC Genom 21:1–8
Google Scholar 
Cherel Y, Gilles J, Handrich Y, Le Maho Y (1994) Nutrient reserve dynamics and energetics during long-term fasting in the king penguin (Aptenodytes patagonicus). J Zool 234:1–12Article 

Google Scholar 
Colles A, Liow LH, Prinzing A (2009) Are specialists at risk under environmental change? Neoecological, paleoecological and phylogenetic approaches. Ecol Lett 12:849–863PubMed 
PubMed Central 
Article 

Google Scholar 
Cristofari R, Bertorelle G, Ancel A, Benazzo A, Le Maho Y, Ponganis PJ et al. (2016) Full circumpolar migration ensures evolutionary utility in the Emperor penguin. Nat Commun 7:1–9Article 

Google Scholar 
Cristofari R, Liu X, Bonadonna F, Cherel Y, Pistorius P, Le Maho Y et al. (2018) Climate-driven range shifts of the king penguin in a fragmented ecosystem. Nat Clim Change 8:245–251Article 

Google Scholar 
Descamps S, Aars J, Fuglei E, Kovacs KM, Lydersen C, Pavlova O et al. (2017) Climate change impacts on wildlife in a High Arctic archipelago-Svalbard, Norway. Glob Change Biol 23:490–502Article 

Google Scholar 
Díaz-Franulic I, Raddatz N, Castillo K, González-Nilo FD, Latorre R (2020) A folding reaction at the C-terminal domain drives temperature sensing in TRPM8 channels. PNAS 117:20298–20304PubMed 
PubMed Central 
Article 

Google Scholar 
Duchamp C, Marmonier F, Denjean F, Lachuer J, Eldershaw TPD, Rouanet JL et al. (1999) Regulatory, cellular and molecular aspects of avian muscle non-shivering thermogenesis. Ornis Fennica 76:151–165
Google Scholar 
Elliott KH, Welcker J, Gaston AJ, Hatch SA, Palace V, Hare JF et al. (2013) Thyroid hormones correlate with resting metabolic rate, not daily energy expenditure, in two charadriiform seabirds. Biol Open 2:580–586CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Frost PGH, Siegfried WR, Greenwood PJ (1975) Arterio-venous heat exchange systems in the Jackass penguin Spheniscus demersus. J Zool 175:231–241Article 

Google Scholar 
Fu J, Gaetani S, Oveisi F, Lo Verme J, Serrano A, Rodríguez de Fonseca F et al. (2003) Oleylethanolamide regulates feeding and body weight through activation of the nuclear receptor PPAR-alpha. Nature 425:90–93CAS 
PubMed 
Article 

Google Scholar 
Gavryushkina A, Heath TA, Ksepka DT, Stadler T, Welch D, Drummond AJ (2017) Bayesian total-evidence dating reveals the recent crown radiation of penguins. Syst Biol 66:57–73PubMed 

Google Scholar 
Geering K, Kraehenbuhl JP, Rossier BC (1987) Maturation of the catalytic alpha unit of Na, K-ATPase during intracellular transport. J Cell Biol 105:2613–2619CAS 
PubMed 
Article 

Google Scholar 
Gilg O, Kovacs KM, Aars J, Fort J, Gauthier G, Grémillet D et al. (2012) Climate change and the ecology and evolution of Arctic vertebrates. Ann N. Y Acad Sci 1249:166–190PubMed 
Article 

Google Scholar 
Goldsmith R, Sladen WJ (1961) Temperature regulation of some Antarctic penguins. J Physiol 157:251–262CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Groscolas R (1990) Metabolic adaptations to fasting in emperor and king penguins. In Penguin Biology (LS Davis & JT Darby ed.), 269-296 San Diego: Academic PressGroscolas R, Robin JP (2001) Long-term fasting and re-feeding in penguins. Comp Biochem Physiol Part A Mol Integr Physiol 128:645–655CAS 

Google Scholar 
Halaas JL, Gajiwala KS, Maffei M, Cohen SL, Chait BT, Rabinowitz D et al. (1995) Weight-reducing effects of the plasma protein encoded by the obese gene. Science 269:543–546CAS 
PubMed 
Article 

Google Scholar 
Han MV, Demuth JP, McGrath CL, Casola C, Hahn MW (2009) Adaptive evolution of young gene duplicates in mammals. Genome Res 19:859–867CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hao Y, Xiong Y, Cheng Y, Song G, Jia C, Qu Y et al. (2019) Comparative transcriptomics of 3 high-altitude passerine birds and their low-altitude relatives. PNAS 116:11851–11856CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Hunt BG, Ometto L, Wurm Y, Shoemaker D, Soojin VY, Keller L et al. (2011) Relaxed selection is a precursor to the evolution of phenotypic plasticity. PNAS 108:15936–15941CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Iannello S, Milazzo P, Belfiore F (2007) Animal and human tissue Na,K-ATPase in normal and insulin-resistant states: regulation, behaviour and interpretative hypothesis on NEFA effects. Obes Rev 8:231–251CAS 
PubMed 
Article 

Google Scholar 
Ishii S, Amano I, Koibuchi N (2021) The role of thyroid hormone in the regulation of cerebellar development. Endocrinol Metab 36:703–716CAS 
Article 

Google Scholar 
Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C et al. (2014) Whole genome analyses resolve early branches in the tree of life of modern birds. Science 346:1320–1331CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kajimura S, Spiegelman BM, Seale P (2015) Brown and beige fat: Physiological roles beyond heat generation. Cell Metab 22:546–559CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kooyman GL, Gentry RL, Bergman WP, Hammel HT (1976) Heat loss in penguins during immersion and compression. Comp Bioch Physiol Part A 54:75–80CAS 
Article 

Google Scholar 
Kumar V, Kutschera VE, Nilsson MA (2015) Genetic signatures of adaptation revealed from transcriptome sequencing of Arctic and red foxes. BMC Genom 16:1–13Article 

Google Scholar 
Kumar S, Stecher G, Tamura K (2016) MEGA7: Molecular evolutionary genetics analysis version 70 for bigger datasets. Mol Biol Evol 33:1870–1874CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lahti DC, Johnson NA, Ajie BC, Otto SP, Hendry AP, Blumstein DT et al. (2009) Relaxed selection in the wild. Trends Ecol Evol 24:487–496PubMed 
Article 

Google Scholar 
Lemberger T, Saladin R, Vazquez M, Assimacopoulos F, Staels B, Desvergne B et al. (1996) Expression of the peroxisome proliferator-activated receptor alpha gene is stimulated by stress and follows a diurnal rhythm. J Biol Chem 271:1764–1769CAS 
PubMed 
Article 

Google Scholar 
Li FG, Li H (2019) A time-dependent genome-wide SNP-SNP interaction analysis of chicken body weight. BMC Genom 20:1–9Article 

Google Scholar 
Li C, Zhang Y, Li J, Kong L, Hu H, Pan H et al. (2014) Two Antarctic penguin genomes reveal insights into their evolutionary history and molecular changes related to the Antarctic environment. Gigascience 3:2047–2217XArticle 

Google Scholar 
Lin Z, Chen L, Chen X, Zhong Y, Yang Y, Xia W et al. (2019) Biological adaptations in the Arctic cervid, the reindeer (Rangifer tarandus). Science 364:eaav6312CAS 
PubMed 
Article 

Google Scholar 
Liu S, Lorenzen ED, Fumagalli M, Li B, Harris K, Xiong Z et al. (2014) Population genomics reveal recent speciation and rapid evolutionary adaptation in polar bears. Cell 157:785–794CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Liu S, Westbury MV, Dussex N, Mitchell KJ, Sinding MH, Heintzman PD et al. (2021) Ancient and modern genomes unravel the evolutionary history of the rhinoceros family. Cell 184:4874–4885CAS 
PubMed 
Article 

Google Scholar 
Lowell BB, Spiegelman BM (2000) Towards a molecular understanding of adaptive thermogenesis. Nature 404:652–660CAS 
PubMed 
Article 

Google Scholar 
Löytynoja A (2013) Phylogeny-aware alignment with PRANK. Methods Mol Biol 1079:155–170Article 

Google Scholar 
Lynch VJ, Bedoya-Reina OC, Ratan A, Sulak M, Drautz-Moses DI, Perry GH et al. (2015) Elephantid genomes reveal the molecular bases of woolly mammoth adaptations to the Arctic. Cell Rep 12:217–228CAS 
PubMed 
Article 

Google Scholar 
Madeira F, Park YM, Lee J, Buso N, Gur T, Madhusoodanan N et al. (2019) The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res 47:636–641Article 

Google Scholar 
Matos-Cruz V, Schneider ER, Mastrotto M, Merriman DK, Bagriantsev SN, Gracheva EO (2017) Molecular prerequisites for diminished cold sensitivity in ground squirrels and hamsters. Cell Rep 21:3329–3337CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Monera OD, Sereda TJ, Zhou NE, Kay CM, Hodges RS (1995) Relationship of sidechain hydrophobicity and α-helical propensity on the stability of the single-stranded amphipathic α-helix. J Pept Sci 1:319–329CAS 
PubMed 
Article 

Google Scholar 
Myers BR, Sigal YM, Julius D (2009) Evolution of Thermal Response Properties in a Cold-Activated TRP Channel. PloS one 4:e5741PubMed 
PubMed Central 
Article 

Google Scholar 
Ng CS, Chen CK, Fan WL, Wu P, Wu SM, Chen JJ et al. (2015) Transcriptomic analyses of regenerating adult feathers in chicken. BMC Genom 16:1–16CAS 
Article 

Google Scholar 
Ohno H, Shinoda K, Spiegelman BM, Kajimura S (2012) PPARg agonists induce a white-to-brown fat conversion through stabilisation of PRDM16 protein. Cell Metab 15:395–404CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Oldham MC, Horvath S, Geschwind DH (2006) Conservation and evolution of gene coexpression networks in human and chimpanzee brains. PNAS 103:17973–17978CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pan H, Cole TL, Bi X, Fang M, Zhou C, Yang Z et al. (2019) High-coverage genomes to elucidate the evolution of penguins. GigaScience 8:1–17Article 

Google Scholar 
Pelleymounter MA, Cullen MJ, Baker MB, Hecht R, Winters D, Boone T et al. (1995) Effects of the obese gene product on body weight regulation in ob/ob mice. Science 269:540–543CAS 
PubMed 
Article 

Google Scholar 
Poirier H, Niot I, Monnot MC, Braissant O, Meunier-Durmort C, Costet P et al. (2001) Differential involvement of peroxi-some-proliferator-activated receptors alpha and delta in fibrate and fatty-acid-mediated inductions of the gene encoding liver fatty-acid-binding protein in the liver and the small intestine. Biochem J 355:481–488CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pond SK, Frost S, Muse SV (2005) HyPhy: hypothesis testing using phylogenies. Bioinform 21:676–679CAS 
Article 

Google Scholar 
R Core Team (2013) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria
Google Scholar 
Ramasamy S, Ometto L, Crava CM, Revadi S, Kaur R, Horner DS et al. (2016) The evolution of olfactory gene families in Drosophila and the genomic basis of chemical-ecological adaptation in Drosophila suzukii. Genome Biol Evol 8:2297–2311PubMed 
PubMed Central 
Article 

Google Scholar 
Ramos B, González-Acuña D, Loyola DE, Johnson WE, Parker PG, Massaro M et al. (2018) Landscape genomics: natural selection drives the evolution of mitogenome in penguins. BMC Genom 19:1–17Article 

Google Scholar 
Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H et al. (2019) g:Profiler: a web server for functional enrichment analysis and conversions of gene lists. Nucleic Acids Res 47:191–198Article 

Google Scholar 
Rey B, Roussel D, Romestaing C, Belouze M, Rouanet JL, Desplanches D et al. (2010) Up-regulation of avian uncoupling protein in cold-acclimated and hyperthyroid ducklings prevents reactive oxygen species production by skeletal muscle mitochondria. BMC Physiol 10:1–12Article 

Google Scholar 
Roussel D, Le Coadic M, Rouanet JL, Duchamp C (2020) Skeletal muscle metabolism in sea-acclimatised king penguins I Thermogenic mechanisms. J Exp Biol 223:pjeb233668Article 

Google Scholar 
Rowland LA, Bal NC, Periasamy M (2015) The role of skeletal‐muscle‐based thermogenic mechanisms in vertebrate endothermy. Biol Rev 90:1279–1297PubMed 
Article 

Google Scholar 
Savini G, Scolari F, Ometto L, Rota-Stabelli O, Carraretto D, Gomulski LM et al. (2021) Viviparity and habitat restrictions may influence the evolution of male reproductive genes in tsetse fly (Glossina) species. BMC Biol 19:1–13Article 

Google Scholar 
Scholander PF (1955) Evolution of climatic adaptation in homeotherms. Evolution 9:15–26Article 

Google Scholar 
Schymkowitz J, Borg J, Stricher F, Nys R, Rousseau F, Serrano L (2005) The FoldX web server: an online force field. Nuc Acids Res 33:382–388Article 

Google Scholar 
Seale P, Kajimura S, Yang W, Chin S, Rohas LM, Uldry M et al. (2007) Transcriptional control of brown fat determination by PRDM16. Cell Metab 6:38–54CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Seale P, Conroe HM, Estall J, Kajimura S, Frontini A, Ishibashi J et al. (2011) Prdm16 determines the thermogenic program of subcutaneous white adipose tissue in mice. J Clin Investig 121:96–105CAS 
PubMed 
Article 

Google Scholar 
Smith MD, Wertheim JO, Weaver S (2015) Less is more: An adaptive branch-site random effects model for efficient detection of episodic diversifying selection. Mol Biol Evol 32:1342–1353CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sonna LA, Fujita J, Gaffin SL, Lilly CM (2002) Effects of heat and cold stress on mammalian gene expression. J Appl Physiol 92:1725–1742CAS 
PubMed 
Article 

Google Scholar 
Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S et al. (2016) The genecards suite: From gene data mining to disease genome sequence analysis. Curr Protoc Bioinform 54:1–30Article 

Google Scholar 
Storey KB, Storey JM (1992) Natural freeze tolerance in ectothermic vertebrates. Annu Rev Physiol 54:619–637CAS 
PubMed 
Article 

Google Scholar 
Storey JD, Bass AJ, Dabney A, Robinson D (2017) qvalue: Q-value estimation for false discovery rate control R package version 2150Supek F, Bošnjak M, Škunca N, Šmuc T (2011) REVIGO summarises and visualises long lists of gene ontology terms. PloS one 6:e21800CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Talbot DA, Duchamp C, Rey B, Hanuise N, Rouanet JL, Sibille B et al. (2004) Uncoupling protein and ATP/ADP carrier increase mitochondrial proton conductance after cold adaptation of king penguins. J Physiol 558:123–135CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tattersall GJ, Sinclair BJ, Withers PC, Fields PA, Seebacher F, Cooper CE et al. (2012) Coping with thermal challenges: Physiological adaptations to environmental temperatures. Compr Physiol 2:2151–2202PubMed 
Article 

Google Scholar 
Teulier L, Rouanet JL, Letexier D, Romestaing C (2010) Cold-acclimation-induced non-shivering thermogenesis in birds is associated with upregulation of avian UCP but not with innate uncoupling or altered ATP efficiency. J Exp Biol 213:2476–2482CAS 
PubMed 
Article 

Google Scholar 
Thomas DB, Fordyce RE (2008) The heterothermic loophole exploited by penguins. Aust J Zool 55:317–321Article 

Google Scholar 
Tigano A, Reiertsen TK, Walters JR, Friesen VL (2018) A complex copy number variant underlies differences in both colour plumage and cold adaptation in a dimorphic seabird. BioRxiv 507384. https://doi.org/10.1101/507384Toyomizu M, Ueda M, Sato S, Seki Y, Sato K, Akiba Y (2002) Cold-induced mitochondrial uncoupling and expression of chicken UCP and ANT mRNA in chicken skeletal muscle. FEBS Lett 529:313–318CAS 
PubMed 
Article 

Google Scholar 
Trucchi E, Gratton P, Whittington JD, Cristofari R, Le Maho Y, Stenseth NC et al. (2014) King penguin demography since the last glaciation inferred from genome-wide data. Proc R Soc B 281:20140528PubMed 
PubMed Central 
Article 

Google Scholar 
Trucchi E, Cristofari R, Le Bohec C (2019) Reply to: The role of ocean dynamics in king penguin range estimation. Nat Clim Change 9:122–122Article 

Google Scholar 
Vermillion KL, Anderson KJ, Hampton M, Andrews MT (2015) Gene expression changes controlling distinct adaptations in the heart and skeletal muscle of a hibernating mammal. Physiol Genom 47:58–74CAS 
Article 

Google Scholar 
Vézina F, Gustowska A, Jalvingh KM, Chastel O, Piersma T (2015) Hormonal correlates and thermoregulatory consequences of moulting on metabolic rate in a northerly wintering shorebird. Physiol Biochem Zool 82:129–142Article 

Google Scholar 
Vianna JA, Fernandes FA, Frugone MJ, Figueiró HV, Pertierra LR, Noll D et al. (2020) Genome-wide analyses reveal drivers of penguin diversification. PNAS 117:22303–22310CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wang Q, Tan X, Jiao S (2014) Analyzing cold tolerance mechanism in transgenic Zebrafish (Danio rerio). PloS one 9:e102492PubMed 
PubMed Central 
Article 

Google Scholar 
Wertheim JO, Murrell B, Smith MD, Kosakovsky Pond SL, Scheffler K (2015) RELAX: Detecting relaxed selection in a phylogenetic framework. Mol Biol Evol 32:820–832CAS 
PubMed 
Article 

Google Scholar 
Wollenberg Valero KC, Pathak R, Prajapati I, Bankston S, Thompson A, Usher J et al. (2014) A candidate multimodal functional genetic network for thermal adaptation. PeerJ 2:e578PubMed 
PubMed Central 
Article 

Google Scholar 
Yang Z (2007) PAML 4: A program package for phylogenetic analysis by maximum likelihood. Mol Biol Evol 24:1586–1591CAS 
PubMed 
Article 

Google Scholar 
Yang J, Bromage TG, Zhao Q, Xu BH, Gao WL, Tian HF et al. (2011) Functional evolution of leptin of Ochotona curzoniae in adaptive thermogenesis driven by cold environmental stress. PloS one 6:e19833CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Yang S, Lu X, Wang Y, Xu L, Chen X, Yang F et al. (2020) A paradigm of thermal adaptation in penguins and elephants by tuning cold activation in TRPM8. PNAS 117:8633–8638CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Yin Y, Wu M, Zubcevic L, Borschel WF, Lander GC, Lee SY (2018) Structure of the cold- and menthol-sensing ion channel TRPM8. Science 359:237–241CAS 
PubMed 
Article 

Google Scholar 
Yudin NS, Larkin DM, Ignatieva EV (2017) A compendium and functional characterization of mammalian genes involved in adaptation to Arctic or Antarctic environments. BMC Genet 18:33–43Article 

Google Scholar 
Zelcer N, Sharpe LJ, Loregger A, Kristiana I, Cook EC, Phan L et al. (2014) The E3 ubiquitin ligase MARCH6 degrades squalene monooxygenase and affects 3-hydroxy-3-methyl-glutaryl coenzyme A reductase and the cholesterol synthesis pathway. Mol Cell Biol 34:1262–1270PubMed 
PubMed Central 
Article 

Google Scholar 
Zhang Y, Proenca R, Maffei M, Barone M, Leopold L, Friedman JM (1994) Positional cloning of the mouse obese gene and its human homologue. Nature 372:425–432CAS 
PubMed 
Article 

Google Scholar 
Zhang G, Li C, Li Q, Li B, Larkin DM, Lee C et al. (2014) Comparative genomics reveals insights into avian genome evolution and adaptation. Science 346:1311–1320CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

We analyzed a total of 2196 morbidity reports received in Israel between November 2018 and October 2021 (Fig. 1A). The majority of these were for adult bats 1,783 (81.2%), of which 1432 (80.3%) were from urban areas (settlements populated by 30,000 or more people) and the remaining 351 (19.7%) were from rural areas, including small villages, nature reserves, and army bases. Out of all these adult cases, the animal’s sex was identified in only 295 cases (16.6%), with 171 (58%) being females and 124 (42%) being males. 413 (18.8% of all cases) were pups, up to 4 months old, from throughout the country, without sex identification.We found a dramatic and significant increase in adult bat morbidity during winter (see Fig. 1B) (p  More

Data collectionBetween June 27 and August 25, 2016, 6600 km of simultaneous visual and passive acoustic line transect surveys were completed on the National Oceanic and Atmospheric Administration (NOAA) ship Henry B. Bigelow5. Survey effort was distributed along saw tooth track lines spanning the continental slope from Virginia (US) to the southern tip of Nova Scotia (Canada) (36–42 N) and on several larger track lines over the abyssal plain. Two teams of visual observers independently recorded sightings of marine mammals using high-powered Fujinon binoculars (25 × 150; Fujifilm, Valhalla, NY) as well as environmental conditions (e.g. sea state) every 30 min.The speed of sound in water was collected three times each day (morning, noon, evening) by measuring conductivity, temperature, and depth (CTD) at specific intervals in the water column. The sound speed closest to the depth of the towed hydrophone array was extracted. On alternating survey days, Simrad EK60 single beam scientific echosounders operating at frequencies of 18, 38, 70, 120 and 200 kHz were used to collect active acoustic data.When possible during daylight hours (06:00–18:00 ET), passive acoustic data were collected continuously using a custom-built linear array composed of eight hydrophone elements and a depth sensor (Keller America Inc. PA7FLE, Newport News, VA) within two oil-filled modular sections separated by 30 m of cable (Fig. 1). The array was towed 300 m behind the vessel at approximately 5–10 m depth while the vessel was in waters more than 100 m deep and underway at speeds of 16–20 km/h. For more details see DeAngelis et al.31, with the only change being that two APC hydrophones and one Reson hydrophone in the aft section were replaced with HTI-96-Min hydrophones (High Tech, Inc., Long Beach, MS). The HTI’s had a flat frequency response from 1 to 30 kHz (− 167 dB re V/uPa ± 1.5 dB). Recordings were made using the acoustical software PAMGuard (v.1.15.02)34. This analysis used the data recorded by the last two 192 kHz sampled hydrophones in the array (MF5 and MF6).Figure 1The linear towed array included eight hydrophone elements and a depth sensor within two oil-filled modular sections separated by 30 m of cable. Six hydrophones sampled at 192 kHz (MF1–MF6) and two sampled at 500 kHz. The hydrophones were connected to two National Instruments sound cards (NI-USB-6356). A high pass filter of 1 kHz was applied by the recording system to reduce the amount of vessel noise in the recordings. This analysis used the passive acoustic data from MF5 and MF6. The schematic is not to scale.Full size imageClick detection and 2D event localizationThe passive acoustic data were filtered using a Butterworth band pass filter (4th order) between 2 and 20 kHz and decimated to 96 kHz to improve sperm whale click resolution. Clicks were automatically detected using the PAMGuard (v.2.01.03) general sperm whale click detector with a trigger threshold of 12 dB.Using PAMGuard’s bearing time display, all detections were reviewed to classify click types and mark click trains as “events” based on consistent changes in bearing, audible sound, ICI and spectral characteristics. Each event was marked to an individual level, tracking a whale from the first to the last detected click15,35. All events containing usual clicks were localized with PAMGuard’s Target Motion Analysis (TMA) module’s 2D simplex optimization algorithm. For further analysis, events were truncated at a slant range of 6500 m (Supplementary Fig. S1).Echosounder analysisA regression analysis was run using the R package MASS36. To account for overdispersion, a negative binomial generalized linear model (GLM) with a log link function was applied to a dataset of the daily acoustic detections33. Echosounder state (active versus passive), month (June, July, August), and habitat type (slope or abyssal) were included as covariates, with the total number of daily detections as the response variable. The track line distance covered per day was used as an offset for effort. The best fitting model was selected based on backwards stepwise selection using Akaike’s information criterion (AIC) and the single-term deletion method using Chi-squared goodness-of-fit tests37.3D localizationExtracting a .wav clip for each click and attributing metadataAn automated process was developed using the R package PAMpal38 (v. 0.14.0) to extract the time of each click in the marked events from PAMGuard databases, generate a .wav clip for each click, and attribute all metadata (e.g., event 2D localization, array depth, radial distance, sea state, and sound speed) necessary for estimating the click depth.Slant delayUsing the methods established by DeAngelis et al.31 and custom Matlab R2021a (MathWorks Inc., Natick, NA) scripts, the multipath arrival of clicks and surface reflected echoes were used to mathematically convert the linear array into a 2D planar array and estimate 3D localizations. Using the .wav clips exported from PAMPal, the time delay between the click and the corresponding surface reflected echo, known as the slant delay, was measured via autocorrelation. Within the autocorrelation solution’s envelope of correlation values, the optimal slant delay was measured using the peak with the highest correlation value above a threshold of 0.02 and within an expected time window after the direct click of 0.0005–0.015 s. Although theoretically a surface reflected echo could have arrived less than a millisecond ( 5 min) were categorized as U shaped, and as shallow ( 1600 m) based on the maximum click depth (Fig. 3).Figure 3Example of click depths (m) over time (min) for events categorized as (a) U shaped and shallow ( 1600 m).Full size imageClick depths were then binned at 400 m intervals to account for an animal’s unknown horizontal movement over time as well as uncertainty in the estimated click depths, and the total time an animal spent within each depth bin was calculated. For each event with a U shaped click depth pattern, the depth bin in which the bottom phase occurred6 was determined. Finally, to assess if a whale was diving in the water column or close to the seafloor, the depth bin in which the 90th percentile of the click depths was recorded was compared to the bin including the seafloor depth. If the whale was more than 400 m above the seafloor, it was determined to be diving in the water column.Distance samplingDepth-corrected average horizontal perpendicular distancesFor each event, a depth-corrected average horizontal perpendicular distance was calculated using the TMA derived perpendicular slant range and the average depth or an assumed depth in the Pythagorean theorem19,31. The weighted mean, first quartile, and third quartile of the average depths were tested as assumed depths for events excluded from 3D localization. If depth was greater than or equal to the slant range the perpendicular distance was coerced to 0, indicating the whale was diving directly below the track line. The resulting distribution of depth-corrected perpendicular distances that aligned most with distance sampling theory was used in the final distance analysis.Acoustic density and abundance estimationThe R package Distance41,42 (v.1.0.4) was used to estimate two separate detection functions based on the uncorrected slant ranges and the depth-corrected perpendicular distances. Half-normal, uniform, and hazard rate key functions were tested with cosine, simple polynomial, and Hermite polynomial adjustment terms. The best fitting models were selected based on the AIC, the Kolmogorov–Smirnov (K–S) test, the Cramer-von Mises (CvM) test, quantile–quantile plots, and visual review of the fitted models43. The probability of detection, abundance, and effective strip half width (ESW) were then estimated for foraging sperm whales.Permitting authorityData used in this manuscript were collected during surveys that were completed under U.S. Marine Mammal Protection Act permit numbers 17355 and 21371 issued to the Northeast Fisheries Science Center. More

Cavicchioli, R. et al. Scientists’ warning to humanity: microorganisms and climate change. Nature Reviews Microbiology 17, 569–586 (2019).CAS 
Article 

Google Scholar 
Azam, F. et al. The ecological role of water-column microbes in the sea. Marine Ecology Progress Series 10, 257–263 (1983).ADS 
Article 

Google Scholar 
Jiao, N. et al. Microbial production of recalcitrant dissolved organic matter: long-term carbon storage in the global ocean. Nature Reviews Microbiology 8, 593–599 (2010).CAS 
Article 

Google Scholar 
Zhang, C. et al. Evolving paradigms in biological carbon cycling in the ocean. National Science Review 5, 481–499 (2018).CAS 
Article 

Google Scholar 
Mestre, M. et al. Sinking particles promote vertical connectivity in the ocean microbiome. Proceedings of the National Academy of Sciences 115, E6799–E6807 (2018).ADS 
CAS 
Article 

Google Scholar 
Baumas, C. M. J. et al. Mesopelagic microbial carbon production correlates with diversity across different marine particle fractions. The ISME Journal 15, 1695–1708 (2021).CAS 
Article 

Google Scholar 
Ortega-Retuerta, E., Joux, F., Jeffrey, W. H. & Ghiglione, J. F. Spatial variability of particle-attached and free-living bacterial diversity in surface waters from the Mackenzie River to the beaufort sea (canadian arctic). Biogeosciences 10, 2747–2759 (2013). BG.ADS 
Article 

Google Scholar 
Ganesh, S., Parris, D. J., DeLong, E. F. & Stewart, F. J. Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone. The ISME Journal 8, 187–211 (2014).CAS 
Article 

Google Scholar 
Chen, S. et al. Interactions between marine group ii archaea and phytoplankton revealed by population correlations in the northern coast of south china sea. Frontiers in Microbiology 12 (2022).Eloe, E. A. et al. Compositional differences in particle-associated and free-living microbial assemblages from an extreme deep-ocean environment. Environmental Microbiology Reports 3, 449–458 (2011).Article 

Google Scholar 
Salazar, G. et al. Particle-association lifestyle is a phylogenetically conserved trait in bathypelagic prokaryotes. Mol Ecol 24, 5692–706 (2015).Article 

Google Scholar 
Karner, M. & Herndl, G. J. Extracellular enzymatic activity and secondary production in free-living and marine-snow-associated bacteria. Marine Biology 113, 341–347 (1992).CAS 
Article 

Google Scholar 
Grossart, H.-P., Tang, K. W., Kiørboe, T. & Ploug, H. Comparison of cell-specific activity between free-living and attached bacteria using isolates and natural assemblages. FEMS Microbiology Letters 266, 194–200 (2007).CAS 
Article 

Google Scholar 
Fierer, N. et al. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. The ISME Journal 6, 1007–1017 (2012).CAS 
Article 

Google Scholar 
Leff, J. W. et al. Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe. Proceedings of the National Academy of Sciences 112, 10967–10972 (2015).ADS 
CAS 
Article 

Google Scholar 
Chen, Y. et al. Large amounts of easily decomposable carbon stored in subtropical forest subsoil are associated with r-strategy-dominated soil microbes. Soil Biology and Biochemistry 95, 233–242 (2016).CAS 
Article 

Google Scholar 
Hou, S. et al. Benefit from decline: the primary transcriptome of Alteromonas macleodii str. Te101 during Trichodesmium demise. The ISME Journal 12, 981–996 (2018).CAS 
Article 

Google Scholar 
Cleveland, C. C., Nemergut, D. R., Schmidt, S. K. & Townsend, A. R. Increases in soil respiration following labile carbon additions linked to rapid shifts in soil microbial community composition. Biogeochemistry 82, 229–240 (2007).CAS 
Article 

Google Scholar 
Ho, A., Di Lonardo, D. P. & Bodelier, P. L. E. Revisiting life strategy concepts in environmental microbial ecology. FEMS Microbiology Ecology 93 (2017).Xing, J. et al. Fluxes, seasonal patterns and sources of various nutrient species (nitrogen, phosphorus and silicon) in atmospheric wet deposition and their ecological effects on Jiaozhou Bay, North China. Sci Total Environ 576, 617–627 (2017).ADS 
CAS 
Article 

Google Scholar 
Zhang, L., Xiong, L., Li, J. & Huang, X. Long-term changes of nutrients and biocenoses indicating the anthropogenic influences on ecosystem in Jiaozhou Bay and Daya Bay, China. Mar Pollut Bull 168, 112406 (2021).CAS 
Article 

Google Scholar 
Zhang, X. et al. Effects of organic nitrogen components from terrestrial input on the phytoplankton community in Jiaozhou Bay. Marine Pollution Bulletin 174, 113316 (2022).CAS 
Article 

Google Scholar 
Sharp, J. et al. Final dissolved organic carbon broad community intercalibration and preliminary use of DOC reference materials. Marine Chemistry 77 (2002).Walters, W. et al. Improved bacterial 16S rrna gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. mSystems 1 (2016).Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nature Biotechnology 37, 852–857 (2019).CAS 
Article 

Google Scholar 
Li, D. et al. MEGAHIT v1.0: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods 102, 3–11 (2016).CAS 
Article 

Google Scholar 
Mikheenko, A., Saveliev, V. & Gurevich, A. MetaQUAST: evaluation of metagenome assemblies. Bioinformatics 32, 1088–90 (2016).CAS 
Article 

Google Scholar 
Yu, K. et al. Recovery of high-qualitied genomes from a deep-inland salt lake using BASALT. bioRxiv https://doi.org/10.1101/2021.03.05.434042 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Kang, D. D. et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 7, e7359 (2019).Article 

Google Scholar 
Wu, Y. W., Simmons, B. A. & Singer, S. W. MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics 32, 605–7 (2016).CAS 
Article 

Google Scholar 
Alneberg, J. et al. Binning metagenomic contigs by coverage and composition. Nat Methods 11, 1144–6 (2014).CAS 
Article 

Google Scholar 
Nayfach, S. et al. New insights from uncultivated genomes of the global human gut microbiome. Nature 568 (2019).Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res 25, 1043–55 (2015).CAS 
Article 

Google Scholar 
Olm, M. R, Brown, C. T. & Banfield, J. F. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. The ISME Journal 5 (2017).Albanese, D. & Donati, C. Large-scale quality assessment of prokaryotic genomes with metashot/prok-quality. F1000Research 10 (2021).Bowers, R. M. et al. Minimum information about a single amplified genome (misag) and a metagenome-assembled genome (mimag) of bacteria and archaea. Nature Biotechnology 35, 725–731 (2017).CAS 
Article 

Google Scholar 
Parks, D. H. et al. A complete domain-to-species taxonomy for bacteria and archaea. Nat Biotechnol 38, 1079–1086 (2020).CAS 
Article 

Google Scholar 
Chaumeil, P. A., Mussig, A. J., Hugenholtz, P. & Parks, D. H. GTDB-Tk: a toolkit to classify genomes with the genome taxonomy database. Bioinformatics (2019).Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32, 1792–7 (2004).CAS 
Article 

Google Scholar 
Capella-Gutierrez, S., Silla-Martinez, J. M. & Gabaldon, T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972–3 (2009).CAS 
Article 

Google Scholar 
Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2–approximately maximum-likelihood trees for large alignments. PloS one 25, e9490–e9490 (2010).ADS 
Article 

Google Scholar 
NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRP367774 (2022).NCBI Sequence Read Archive https://identifiers.org/ncbi/insdc.sra:SRP367809 (2022).Tao, J. Jiaozhou bay 16S rDNA & metagenome dataset. figshare https://doi.org/10.6084/m9.figshare.19690459.v6 (2022). More

Amundsen, T. In Animal Signals: Signalling and Signal Design in Animal Communication (eds. Espmark, Y., Amundsen, T. & Rosenqvist, G.) 133–154 (Tapir Academic Press, 2000).Amundsen, T. Why are female birds ornamented? Trends Ecol. Evol. 15, 149–155 (2000).CAS 
PubMed 
Article 

Google Scholar 
Lande, R. Sexual dimorphism, sexual selection and adaptation in polygenic characters. Evolution 34, 292–305 (1980).PubMed 
Article 

Google Scholar 
Poissant, J., Wilson, A. J. & Coltman, D. W. Sex-specific genetic variance and the evolution of sexual size dimorphism: a systematic review of cross-sex genetic correlations. Evolution 64, 97–107 (2009).PubMed 
Article 

Google Scholar 
Nordeide, J. T., Kekäläinen, J., Janhunen, M. & Kortet, R. Female ornaments revisited—are they correlated with offspring quality? J. Anim. Ecol. 82, 26–38 (2013).PubMed 
Article 

Google Scholar 
Prum, R. O. The Evolution of Beauty: How Darwin’s Forgotten Theory of Mate Choice Shapes the Animal World and Us (Doubleday, 2017).Clark, C. J. & Rankin, D. Subtle, pervasive genetic correlation between the sexes in the evolution of dimorphic hummingbird tail ornaments. Evolution 74, 528–543 (2020).PubMed 
Article 

Google Scholar 
LeBas, N. R. Female finery is not for males. Trends Ecol. Evol. 21, 170–173 (2006).PubMed 
Article 

Google Scholar 
Kraaijeveld, K., Kraaijeveld-Smit, F. J. L. & Komdeur, J. The evolution of mutual ornamentation. Anim. Behav. 74, 657–677 (2007).Article 

Google Scholar 
Tobias, J. A., Montgomerie, R. & Lyon, B. E. The evolution of female ornaments and weaponry: social selection, sexual selection and ecological competition. Philos. Trans. R. Soc. B 367, 2274–2293 (2012).Article 

Google Scholar 
Hare, R. M. & Simmons, L. W. Sexual selection and its evolutionary consequences in female animals. Biol. Rev. 94, 1464–7931 (2019).Article 

Google Scholar 
Hernández, A., Martínez-Gómez, M., Beamonte-Barrientos, R. & Montoya, B. Colourful traits in female birds relate to individual condition, reproductive performance and male-mate preferences: a meta-analytic approach. Biol. Lett. 17, 20210283 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
Tsuboi, M., Gonzalez-Voyer, A., Höglund, J. & Kolm, N. Ecology and mating competition influence sexual dimorphism in Tanganyikan cichlids. Evol. Ecol. 26, 171–185 (2012).Article 

Google Scholar 
Andersson, M. Sexual Selection (Princeton Univ. Press, 1994).Doutrelant, C., Fargevieille, A. & Grégoire, A. Evolution of female coloration: what have we learned from birds in general and blue tits in particular. Adv. Study Behav. 52, 123–202 (2020).Article 

Google Scholar 
Dunn, P. O., Armenta, J. K. & Whittingham, L. A. Natural and sexual selection act on different axes of variation in avian plumage color. Sci. Adv. 1, e1400155 (2015).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Cotton, S., Fowler, K. & Pomiankowski, A. Do sexual ornaments demonstrate heightened condition-dependent expression as predicted by the handicap hypothesis? Proc. Biol. Sci. 271, 771–783 (2004).PubMed 
PubMed Central 
Article 

Google Scholar 
Bonduriansky, R. & Rowe, L. Sexual selection, genetic architecture, and the condition dependence of body shape in the sexually dimorphic fly Prochyliza xanthostoma (Piophilidae). Evolution 59, 138–151 (2005).PubMed 
Article 

Google Scholar 
Johnstone, R. A., Rands, S. A. & Evans, M. R. Sexual selection and condition-dependence. J. Evol. Biol. 22, 2387–2394 (2009).CAS 
PubMed 
Article 

Google Scholar 
Cotton, S., Fowler, K. & Pomiankowski, A. Heightened condition dependence is not a general feature of male eyespan in stalk-eyed flies (Diptera: Diopsidae). J. Evol. Biol. 17, 1310–1316 (2004).CAS 
PubMed 
Article 

Google Scholar 
David, P. et al. Male sexual ornament size but not asymmetry reflects condition in stalk-eyed flies. Proc. R. Soc. Lond. B 265, 2211–2216 (1998).Article 

Google Scholar 
Bolund, E., Schielzeth, H. & Forstmeier, W. No heightened condition dependence of zebra finch ornaments—a quantitative genetic approach. J. Evol. Biol. 23, 586–597 (2010).CAS 
PubMed 
Article 

Google Scholar 
Zahavi, A. Mate selection-a selection for a handicap. J. Theor. Biol. 53, 205–214 (1975).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Meunier, J., Figueiredo Pinto, S., Burri, R. & Roulin, A. Eumelanin-based coloration and fitness parameters in birds: a meta-analysis. Behav. Ecol. Sociobiol. 65, 559–567 (2011).Article 

Google Scholar 
Weaver, R. J., Santos, E. S. A., Tucker, A. M., Wilson, A. E. & Hill, G. E. Carotenoid metabolism strengthens the link between feather coloration and individual quality. Nat. Commun. 9, 73 (2018).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
White, T. E. Structural colours reflect individual quality: a meta-analysis. Biol. Lett. 16, 20200001 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Cohen, J. Statistical Power Analysis for the Behavioral Sciences (Taylor & Francis Inc., 1988)Andersson, M. Sexual selection, natural selection and quality advertisement. Biol. J. Linn. Soc. 17, 375–393 (1982).Article 

Google Scholar 
Walther, B. A. & Clayton, D. H. Elaborate ornaments are costly to maintain: evidence for high maintenance handicaps. Behav. Ecol. 16, 89–95 (2005).Article 

Google Scholar 
Folstad, I. & Karter, A. K. Parasites, bright males and the immunocompetence handicap. Am. Nat. 139, 603–622 (1992).Article 

Google Scholar 
Alonso-Alvarez, C., Bertrand, S., Faivre, B., Chastel, O. & Sorci, G. Testosterone and oxidative stress: the oxidation handicap hypothesis. Proc. R. Soc. Lond. B 274, 819–825 (2007).CAS 

Google Scholar 
Weaver, R. J., Koch, R. E. & Hill, G. E. What maintains signal honesty in animal colour displays used in mate choice? Philos. Trans. R. Soc. B 372, 20160343 (2017).Article 

Google Scholar 
Emlen, D. J., Warren, I. A., Johns, A., Dworkin, I. & Lavine, L. C. A mechanism of extreme growth and reliable signaling in sexually selected ornaments and weapons. Science 337, 860–864 (2012).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Huhta, E. Plumage brightness of prey increases predation risk: an among-species comparison. Ecology 84, 1793–1799 (2003).Article 

Google Scholar 
Tibbetts, E. A. & Dale, J. A socially enforced signal of quality in a paper wasp. Nature 432, 18–222 (2004).Article 

Google Scholar 
Webster, M. S., Ligon, R. A. & Leighton, G. M. Social costs are an underappreciated force for honest signalling in animal aggregations. Anim. Behav. 143, 167–176 (2018).Article 

Google Scholar 
Sheldon, B. C. Differential allocation: tests, mechanisms and implications. Trends Ecol. Evol. 15, 397–402 (2000).CAS 
PubMed 
Article 

Google Scholar 
Johnstone, R. A., Reynolds, J. D. & Deutsch, J. C. Mutual mate choice and sex differences in choosiness. Evolution 50, 1382–1391 (1996).PubMed 
Article 

Google Scholar 
Promislow, D. E. L., Montgomerie, R. & Martin, T. E. Mortality costs of sexual dimorphism in birds. Proc. R. Soc. Lond. B 250, 143–150 (1992).ADS 
Article 

Google Scholar 
Guindre-Parker, S. & Love, O. P. Revisiting the condition-dependence of melanin-based plumage. J. Avian Biol. 45, 29–33 (2014).Article 

Google Scholar 
Roulin, A. & Dijkstra, C. Genetic and environmental components of variation in eumelanin and phaeomelanin sex-traits in the barn owl. Heredity 90, 359–364 (2003).CAS 
PubMed 
Article 

Google Scholar 
Jawor, J. M. & Breitwisch, R. Melanin ornaments, honesty, and sexual selection. Auk 120, 249–265 (2003).Article 

Google Scholar 
Gunderson, A. R., Frame, A. M., Swaddle, J. P. & Forsyth, M. H. Resistance of melanized feathers to bacterial degradation: is it really so black and white? J. Avian Biol. 39, 539–545 (2008).Article 

Google Scholar 
Ruiz-de-Castañeda, R., Burtt, E. H. Jr., González-Braojos, S. & Moreno, J. Bacterial degradability of an intrafeather unmelanized ornament: a role for feather-degrading bacteria in sexual selection? Biol. J. Linn. Soc. 105, 409–419 (2012).Article 

Google Scholar 
Tazzyman, S. J., Iwasa, Y. & Pomiankowski, A. Signaling efficacy drives the evolution of larger sexual ornaments by sexual selection. Evolution 68, 216–229 (2014).PubMed 
Article 

Google Scholar 
Dale, J., Dey, C. J., Delhey, K., Kempenaers, B. & Valcu, M. The effects of life history and sexual selection on male and female plumage colouration. Nature 527, 367–370 (2015).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Guilford, T. & Dawkins, M. S. Receiver psychology and the evolution of animal signals. Anim. Behav. 42, 1–14 (1991).Article 

Google Scholar 
Tazzyman, S. J., Iwasa, Y. & Pomiankowski, A. The handicap process favors exaggerated, rather than reduced, sexual ornaments. Evolution 68, 2534–2549 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Peters, J. L. et al. Assessing publication bias in meta-analyses in the presence of between-study heterogeneity. J. R. Stat. Soc. Ser. A. 173, 575–591 (2010).MathSciNet 
Article 

Google Scholar 
Dumbacher, J. P. & Fleischer, R. C. Phylogenetic evidence for colour pattern convergence in toxic pitohuis: Müllerian mimicry in birds? Proc. Biol. Sci. 268, 1971–1976 (2001).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Jønsson, K. A., Delhey, K., Sangster, G., Ericson, P. G. P. & Irestedt, M. The evolution of mimicry of friarbirds by orioles (Aves: Passeriformes) in Australo-Pacific archipelagos. Proc. R. Soc. B Biol. Sci. B 283, 20160409 (2016).Article 

Google Scholar 
Ord, T. J. & Stuart-Fox, D. Ornament evolution in dragon lizards: multiple gains and widespread losses reveal a complex history of evolutionary change. J. Evol. Biol. 19, 797–808 (2006).CAS 
PubMed 
Article 

Google Scholar 
Moher, D., Liberati, A., Tetzlaff, J., Altman, D. G. & PRISMA Group. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. PLoS Med. 6, e1000097 (2009).O’Dea, R. E. et al. Preferred reporting items for systematic reviews and meta-analyses in ecology and evolutionary biology: a PRISMA extension. Biol. Rev. 96, 1695–1722 (2021).PubMed 
Article 

Google Scholar 
LeBas, N. R., Hockham, L. R. & Ritchie, M. G. Nonlinear and correlational sexual selection on ‘honest’ female ornamentation. Proc. R. Soc. Lond. B 270, 2159–2165 (2003).Article 

Google Scholar 
Ouzzani, M., Hammady, H., Fedorowicz, Z. & Elmagarmid, A. Rayyan—a web and mobile app for systematic reviews. Syst. Rev. 5, 210 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Rohatgi, A. WebPlotDigitizer. Software version 4.5. https://automeris.io/WebPlotDigitizer (2000).Sidney, S. Nonparametric Statistics for the Behavioral Sciences (McGraw-Hill,1956).Friedman, H. Simplified determination of statistical power, magnitude of effect and research sample sizes. Educ. Psychol. Meas. 42, 521–526 (1982).Article 

Google Scholar 
Nakagawa, S. & Cuthill, I. C. Effect size, confidence interval and statistical significance: a practical guide for biologists. Biol. Rev. 82, 591–605 (2007).PubMed 
Article 

Google Scholar 
Verhulst, S. & Nilsson, J. A. The timing of birds’ breeding seasons: a review of experiments that manipulated timing of breeding. Philos. Trans. R. Soc. Lond. B 363, 399–410 (2008).Article 

Google Scholar 
Brown, M. E. In Current Ornithology (eds. Nolan, V. & Ketterson, E. D.) 67–135 (Plenum Press, 1996).Labocha, M. K. & Hayes, J. P. Morphometric indices of body condition in birds: a review. J. Ornithol. 153, 1–22 (2012).Article 

Google Scholar 
Sánchez, C. A. et al. On the relationship between body condition and parasite infection in wildlife: a review and meta-analysis. Ecol. Lett. 20, 1869–1884 (2018).Article 

Google Scholar 
Arnholt, A. T. & Evans, B. BSDA: Basic statistics and data analysis. R package version 1.2.0. https://cran.r-project.org/package=BSDA (2017).Jackson, D., White, I. R., Price, M., Copas, J. & Riley, R. D. Borrowing of strength and study weights in multivariate and network meta-analysis. Stat. Methods Med. Res. 26, 2853–2868 (2017).MathSciNet 
PubMed 
Article 

Google Scholar 
Hadfield, J. D. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R Package. J. Stat. Softw. 33, 1–22 (2010).Article 

Google Scholar 
Jetz, W., Thomas, G. H., Joy, J. B., Hartmann, K. & Mooers, A. O. The global diversity of birds in space and time. Nature 491, 444–448 (2012).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Nakagawa, S. & De Villemereuil, P. A general method for simultaneously accounting for phylogenetic and species sampling uncertainty via Rubin’s rules in comparative analysis. Syst. Biol. 68, 632–641 (2019).PubMed 
Article 

Google Scholar 
Cinar, O., Nakagawa, S. & Viechtbauer, W. Phylogenetic multilevel meta-analysis: a simulation study on the importance of modeling the phylogeny. Methods Ecol. Evol. 13, 383–395 (2022).Article 

Google Scholar 
Viechtbauer, W. Conducting meta-analyses in R with the metafor package. J. Stat. Softw. 36, 1–48 (2010).Article 

Google Scholar 
Egger, M., Davey Smith, G., Schneider, M. & Minder, C. Bias in meta-analysis detected by a simple, graphical test. BMJ 315, 629–634 (1997).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Duval, S. & Tweedie, R. Trim and fill: a simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics 56, 455–463 (2000).CAS 
PubMed 
MATH 
Article 

Google Scholar 
Nakagawa, S. et al. Methods for testing publication bias in ecological and evolutionary meta-analyses. Methods Ecol. Evol. 13, 4–21 (2022).Article 

Google Scholar 
Nakagawa, S. & Santos, E. S. A. Methodological issues and advances in biological meta-analysis. Evol. Ecol. 26, 1253–1274 (2012).Article 

Google Scholar 
Billerman, S. M., Keeney, B. K., Rodewald, P. G. & Schulenberg, T. S. Birds of the World (Cornell Laboratory of Ornithology, 2000). More