Ecology

1.Cunningham, A. A., Daszak, P. & Wood, J. L. N. One health, emerging infectious diseases and wildlife: two decades of progress? Philos. Trans. R. Soc. B Biol. Sci. 372, 20160167 (2017).2.Suzan, G., Esponda, F., Carrasco-Hernández, R. & Aguirre, A. A. in New Directions in Conservation Medicine: Applied Cases of Ecological Health (eds. Aguirre, A. A., Ostfeld, R. & Daszak, P.). 135–150 (Oxford University Press USA, 2012).3.Hussain, S., Ram, M. S., Kumar, A., Shivaji, S. & Umapathy, G. Human presence increases parasitic load in endangered lion-tailed macaques (Macaca silenus) in its fragmented rainforest habitats in Southern India. PLoS ONE 8, 1–8 (2013).
Google Scholar 
4.Junge, R. E., Barrett, M. A. & Yoder, A. D. Effects of anthropogenic disturbance on indri (Indri indri) health in Madagascar. Am. J. Primatol. 73, 632–642 (2011).PubMed 
Article 

Google Scholar 
5.Friggens, M. M. & Beier, P. Anthropogenic disturbance and the risk of flea-borne disease transmission. Oecologia 164, 809–820 (2010).PubMed 
Article 

Google Scholar 
6.Woodroffe, R. et al. Contact with domestic dogs increases pathogen exposure in endangered African wild dogs (Lycaon pictus). PLoS ONE 7, e30099 (2012).7.Crowl, T. A., Crist, T. O., Parmenter, R. R., Belovsky, G. & Lugo, A. E. The spread of invasive species and infectious disease as drivers of ecosystem change. Front. Ecol. Environ. 6, 238–246 (2008).Article 

Google Scholar 
8.Keesing, F., Holt, R. D. & Ostfeld, R. S. Effects of species diversity on disease risk. Ecol. Lett. 9, 485–498 (2006).CAS 
PubMed 
Article 

Google Scholar 
9.Gámez-Virués, S. et al. Landscape simplification filters species traits and drives biotic homogenization. Nat. Commun. 6, 8568 (2015).10.Alberdi, A., Aizpurua, O., Bohmann, K., Zepeda-Mendoza, M. L. & Gilbert, M. T. P. Do vertebrate gut metagenomes confer rapid ecological adaptation? Trends Ecol. Evol. 31, 689–699 (2016).PubMed 
Article 
PubMed Central 

Google Scholar 
11.Hooper, L. V., Littman, D. R. & Macpherson, A. J. Interactions between the microbiota and the immune system. Science 336, 1268–1273 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
12.Shapira, M. Gut microbiotas and host evolution: scaling up symbiosis. Trends Ecol. Evol. 31, 539–549 (2016).PubMed 
Article 
PubMed Central 

Google Scholar 
13.Brugman, S. et al. A comparative review on microbiota manipulation: lessons from fish, plants, livestock, and human research. Front. Nutr. 5, 1–15 (2018).Article 
CAS 

Google Scholar 
14.Wasimuddin et al. Astrovirus infections induce age-dependent dysbiosis in gut microbiomes of bats. ISME J. 12, 2883–2893 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
15.Wasimuddin et al. Adenovirus infection is associated with altered gut microbial communities in a non-human primate. Sci. Rep. 9, 1–12 (2019).CAS 
Article 

Google Scholar 
16.Wilkins, L. J., Monga, M. & Miller, A. W. Defining dysbiosis for a cluster of chronic diseases. Sci. Rep. 9, 1–10 (2019).
Google Scholar 
17.Brüssow, H. Problems with the concept of gut microbiota dysbiosis. Microb. Biotechnol. 13, 423–434 (2020).PubMed 
Article 
PubMed Central 

Google Scholar 
18.Otto, S. P. Adaptation, speciation and extinction in the Anthropocene. Proc. R. Soc. B Biol. Sci. 285, 20182047 (2018).19.Amato, K. R. et al. Habitat degradation impacts black howler monkey (Alouatta pigra) gastrointestinal microbiomes. ISME J. 7, 1344–1353 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
20.Ingala, M. R., Becker, D. J., Bak Holm, J., Kristiansen, K. & Simmons, N. B. Habitat fragmentation is associated with dietary shifts and microbiota variability in common vampire bats. Ecol. Evol. https://doi.org/10.1002/ece3.5228 (2019)21.Juan, P. A. S., Hendershot, J. N., Daily, G. C. & Fukami, T. Land-use change has host-specificinfluenc on avian gut microbiomes. ISME J. https://doi.org/10.1038/s41396-019-0535-4 (2019)22.Barelli, C. et al. Habitat fragmentation is associated to gut microbiota diversity of an endangered primate: implications for conservation. Sci. Rep. 5, 14862 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
23.de Juan, S., Thrush, S. F. & Hewitt, J. E. Counting on β-diversity to safeguard the resilience of estuaries. PLoS ONE 8, 1–11 (2013).
Google Scholar 
24.Burnham, K. P. & Anderson, D. R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach (Springer, 2002).25.Oksanen, J. et al. vegan: Community Ecology Package. R package version 2.5-5. https://github.com/vegandevs/vegan (2019).26.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 
PubMed Central 

Google Scholar 
27.Cohen, J. Statistical Power Analysis for the Behavioral Sciences (Lawrence Erlbaum Associates, 1988).28.Gillingham, M. A. F. et al. Offspring microbiomes differ across breeding sites in a panmictic species. Front. Microbiol. 10, 35 (2019).29.Anderson, M. J. Distance-based tests for homogeneity of multivariate dispersions. Biometrics 62, 245–253 (2006).PubMed 
Article 
PubMed Central 

Google Scholar 
30.Mandal, S. et al. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb. Ecol. Heal. Dis. 26, 1–7 (2015).
Google Scholar 
31.Douglas, G. M. et al. PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol. 38, 669–673 (2020).Article 
CAS 

Google Scholar 
32.Louca, S. & Doebeli, M. Efficient comparative phylogenetics on large trees. Bioinformatics 34, 1053–1055 (2018).CAS 
PubMed 
Article 

Google Scholar 
33.Barbera, P. et al. EPA-ng: massively parallel evolutionary placement of genetic sequences. Syst. Biol. 68, 365–369 (2019).PubMed 
Article 

Google Scholar 
34.Czech, L., Barbera, P. & Stamatakis, A. Genesis and Gappa: processing, analyzing and visualizing phylogenetic (placement) data. Bioinformatics 36, 3263–3265 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
35.Nyhus, P. J. Human—wildlife conflict and coexistence. Annu. Rev. Environ. Resour. 41, 143–171 (2016).Article 

Google Scholar 
36.Foden, W. B. et al. Climate change vulnerability assessment of species. WIREs Clim. Chang. 10, 1–36 (2019).Article 

Google Scholar 
37.Beck, J. M. et al. Multicenter comparison of lung and oral microbiomes of HIV-infected and HIV-uninfected individuals. Am. J. Respir. Crit. Care Med. 192, 1335–1344 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
38.Pita, L., Rix, L., Slaby, B. M., Franke, A. & Hentschel, U. The sponge holobiont in a changing ocean: from microbes to ecosystems. Microbiome 6, 46 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
39.Rosado, P. M. et al. Marine probiotics: increasing coral resistance to bleaching through microbiome manipulation. ISME J. 13, 921–936 (2019).CAS 
PubMed 
Article 

Google Scholar 
40.Wang, L. et al. Corals and their microbiomes are differentially affected by exposure to elevated nutrients and a natural thermal anomaly. Front. Mar. Sci. 5, 1–16 (2018).Article 

Google Scholar 
41.Zaneveld, J. R., McMinds, R. & Thurber, R. V. Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat. Microbiol. 2, 17121 (2017).42.Rocca, J. D. et al. The Microbiome Stress Project: toward a global meta-analysis of environmental stressors and their effects on microbial communities. Front. Microbiol. 10, 3272 (2019).43.Gillingham, M. A. F. et al. Bioaccumulation of trace elements affects chick body condition and gut microbiome in greater flamingos. Sci. Total Environ. 761, 143250 (2020).44.Chase, J. M. Stochastic community assembly causes higher biodiversity in more productive environments. Science 328, 1388–1392 (2010).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
45.Jiménez, R. R., Alvarado, G., Estrella, J. & Sommer, S. Moving beyond the host: unraveling the skin microbiome of endangered Costa Rican amphibians. Front. Microbiol. 10, 1–18 (2019).Article 

Google Scholar 
46.Wang, J. et al. Phylogenetic beta diversity in bacterial assemblages across ecosystems: deterministic versus stochastic processes. ISME J. 7, 1310–1321 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
47.Chase, J. M. & Myers, J. A. Disentangling the importance of ecological niches from stochastic processes across scales. Philos. Trans. R. Soc. B Biol. Sci. 366, 2351–2363 (2011).Article 

Google Scholar 
48.Pound, K. L., Lawrence, G. B. & Passy, S. I. Beta diversity response to stress severity and heterogeneity in sensitive versus tolerant stream diatoms. Divers. Distrib. 25, 374–384 (2019).Article 

Google Scholar 
49.Zhou, J. & Ning, D. Stochastic Community Assembly: does it matter in microbial ecology? Microbiol. Mol. Biol. Rev. 81, 1–32 (2017).Article 

Google Scholar 
50.Nicholas, R. A. J. & Ayling, R. D. Mycoplasma bovis: disease, diagnosis, and control. Res. Vet. Sci. 74, 105–112 (2003).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
51.Ley, D. H. in Diseases of Poultry (eds. et al.) (Blackwell Publishing, 2008).52.Groebel, K., Hoelzle, K., Wittenbrink, M. M., Ziegler, U. & Hoelzle, L. E. Mycoplasma suis invades porcine erythrocytes. Infect. Immun. 77, 576–584 (2009).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
53.do Nascimento, N. C., Santos, A. P., Guimaraes, A. M. S., Sanmiguel, P. J. & Messick, J. B. Mycoplasma haemocanis—the canine hemoplasma and its feline counterpart in the genomic era. Vet. Res. 43, 66 (2012).54.Hardham, J. M. et al. Transfer of Bacteroides splanchnicus to Odoribacter gen. nov. as Odoribacter splanchnicus comb. nov., and description of Odoribacter denticanis sp. nov., isolated from the crevicular spaces of canine periodontitis patients. Int. J. Syst. Evol. Microbiol. 58, 103–109 (2008).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
55.Kaakoush, N. O. Insights into the role of Erysipelotrichaceae in the human host. Front. Cell. Infect. Microbiol. 5, 1–4 (2015).Article 
CAS 

Google Scholar 
56.Ormerod, K. L. et al. Genomic characterization of the uncultured Bacteroidales family S24-7 inhabiting the guts of homeothermic animals. Microbiome 4, 1–17 (2016).Article 

Google Scholar 
57.Herrmann, E. et al. RNA-based stable isotope probing suggests Allobaculum spp. as particularly active glucose assimilators in a complex murine microbiota cultured in vitro. Biomed Res. Int. 2017, 1829685 (2017).58.Greetham, H. L. et al. Allobaculum stercoricanis gen. nov., sp. nov., isolated from canine feces. Anaerobe 10, 301–307 (2004).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
59.Silva, Y. P., Bernardi, A. & Frozza, R. L. The role of short-chain fatty acids from gut microbiota in gut-brain communication. Front. Endocrinol. 11, 1–14 (2020).60.Wiegel, J., Tanner, R. & Rainey, F. A. in The Prokaryotes: Volume 4: Bacteria: Firmicutes, Cyanobacteria (eds. Dworkin, M., Falkow, S., Rosenberg, E., Schleifer, K.-H. & Stackebrandt, E.) 654–678 (Springer US, 2006).61.Tamanai-Shacoori, Z. et al. Roseburia spp.: a marker of health? Future Microbiol 12, 157–170 (2017).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
62.Freier, T. A., Beitz, D. C., Li, L. & Hartman, P. A. Characterization of Eubacterium coprostanoligenes sp. nov., a Cholesterol-Reducing Anaerobe. Int. J. Syst. Bacteriol. 44, 137–142 (1994).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
63.Venegas, D. P. et al. Short chain fatty acids (SCFAs)-mediated gut epithelial and immune regulation and its relevance for inflammatory bowel diseases. Front. Immunol. 10, 277 (2019).64.MetaCyc. MetaCyc Pathway: pyrimidine deoxyribonucleotides biosynthesis from CTP. https://biocyc.org/META/NEW-IMAGE?type=PATHWAY&object=PWY-7210&show-citations=T (2020).65.Caspi, R. et al. The MetaCyc database of metabolic pathways and enzymes. Nucleic Acids Res. 46, D633–D639 (2018).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
66.MetaCyc. MetaCyc Pathway: poly(glycerol phosphate) wall teichoic acid biosynthesis. https://biocyc.org/META/NEW-IMAGE?type=PATHWAY&object=TEICHOICACID-PWY (2020).67.Brown, S., Santa Maria, J. P. & Walker, S. Wall teichoic acids of gram-positive bacteria. Annu. Rev. Microbiol. 67, 313–336 (2013).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
68.MetaCyc. MetaCyc Pathway: L-lysine biosynthesis II. https://metacyc.org/META/NEW-IMAGE?type=PATHWAY&object=PWY-2941 (2020).69.Hutton, C. A., Perugini, M. A. & Gerrard, J. A. Inhibition of lysine biosynthesis: an evolving antibiotic strategy. Mol. Biosyst. 3, 458–465 (2007).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
70.Wanner, S. et al. Wall teichoic acids mediate increased virulence in Staphylococcus aureus. Nat. Microbiol. 2, 1–12 (2017).
Google Scholar 
71.MetaCyc. MetaCyc Pathway: formaldehyde assimilation II (assimilatory RuMP Cycle). https://biocyc.org/META/NEW-IMAGE?type=PATHWAY&object=PWY-1861 (2020).72.Chen, N. H., Djoko, K. Y., Veyrier, F. J. & McEwan, A. G. Formaldehyde stress responses in bacterial pathogens. Front. Microbiol. 7, 1–17 (2016).
Google Scholar 
73.Tauseef, S. M., Premalatha, M., Abbasi, T. & Abbasi, S. A. Methane capture from livestock manure. J. Environ. Manag. 117, 187–207 (2013).CAS 
Article 

Google Scholar 
74.Dale, V. H., Brown, S., Calderón, M. O., Montoya, A. S. & Martínez, R. E. Estimating baseline carbon emissions for the eastern Panama Canal watershed. Mitig. Adapt. Strateg. Glob. Chang 8, 323–348 (2003).Article 

Google Scholar 
75.Schmid, J. et al. Ecological drivers of Hepacivirus infection in a neotropical rodent inhabiting landscapes with various degrees of human environmental change. Oecologia https://doi.org/10.1007/s00442-018-4210-7 (2018)76.Adler, G. H. & Beatty, R. P. Changing reproductive rates in a neotropical forest rodent, Proechimys semispinosus. J. Anim. Ecol. 66, 472 (1997).Article 

Google Scholar 
77.Adler, G. H. Fruit and seed exploitation by Central American spiny rats, Proechimys semispinosus. Stud. Neotrop. Fauna Environ. 30, 237–244 (1995).78.Hoch, G. A. & Adler, G. H. Removal of black palm (Astrocaryum standleyanum) seeds by spiny rats (Proechimys semispinosus). J. Trop. Ecol. 13, 51–58 (1997).Article 

Google Scholar 
79.Endries, M. J. & Adler, G. H. Spacing patterns of a tropical forest rodent, the spiny rat (Proechimys semispinosus), in Panama. J. Zool. 265, 147–155 (2005).Article 

Google Scholar 
80.Adler, G. H. The island syndrome in isolated populations of a tropical forest rodent. Oecologia 108, 694–700 (1996).PubMed 
Article 
PubMed Central 

Google Scholar 
81.Caporaso, J. G. et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. PNAS 108, 4516–4522 (2011).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
82.Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
83.Menke, S. et al. Oligotyping reveals differences between gut microbiomes of free-ranging sympatric Namibian carnivores (Acinonyx jubatus, Canis mesomelas) on a bacterial species-like level. Front. Microbiol. 5, 526 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
84.Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
85.Callahan, B. J. et al. DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
86.Callahan, B. J., Mcmurdie, P. J. & Holmes, S. P. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 11, 2639–2643 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
87.Quast, C. et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 41, 590–596 (2013).Article 
CAS 

Google Scholar 
88.Yilmaz, P. et al. The SILVA and ‘All-species Living Tree Project (LTP)’ taxonomic frameworks. Nucleic Acids Res. 42, 643–648 (2014).Article 
CAS 

Google Scholar 
89.Glöckner, F. O. et al. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J. Biotechnol. 261, 169–176 (2017).PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
90.Pedregosa, F. et al. Scikit-learn: machine learning in python. J. Mach. Learn. Res. 12, 2825–2830 (2011).91.Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2—āpproximately maximum-likelihood trees for large alignments. PLoS ONE 5, e9490 (2010).92.Huson, D. H. & Scornavacca, C. Dendroscope 3: an interactive tool for rooted phylogenetic trees and networks. Syst. Biol. 61, 1061–1067 (2012).PubMed 
Article 
PubMed Central 

Google Scholar 
93.R. Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. http://www.r-project.org/index.html (2017).94.McMurdie, P. J. & Holmes, S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 8, e61217 (2013).95.Davis, N. M., Proctor, Di. M., Holmes, S. P., Relman, D. A. & Callahan, B. J. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome 6, 1–14 (2018).Article 

Google Scholar 
96.Shannon, C. E. A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423 (1948).Article 

Google Scholar 
97.Faith, D. P. Conservation evaluation and phylogenetic diversity. Biol. Conserv. 61, 1–10 (1992).Article 

Google Scholar 
98.Bates, D., Mächler, M., Bolker, B. M. & Walker, S. C. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).99.Mcmurdie, P. J. & Holmes, S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput. Biol. 10, e1003531 (2014).100.Kim, Y. S., Unno, T., Kim, B.-Y. & Park, M. Sex differences in gut microbiota. World J. Mens. Health 38, 48–60 (2020).PubMed 
Article 
PubMed Central 

Google Scholar 
101.Kolodny, O. et al. Coordinated change at the colony level in fruit bat fur microbiomes through time. Nat. Ecol. Evol. 3, 116–124 (2019).PubMed 
Article 
PubMed Central 

Google Scholar 
102.Kartzinel, T. R., Hsing, J. C., Musili, P. M., Brown, B. R. P. & Pringle, R. M. Covariation of diet and gut microbiome in African megafauna. Proc. Natl Acad. Sci. USA 116, 23588–23593 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
103.Zuur, A., Ieno, E. N., Walker, N., Saveliev, A. A. & Smith, G. M. Mixed Effects Models and Extensions in Ecology with R (Springer Science & Business Media, 2009).104.Nakagawa, S. & Schielzeth, H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 4, 133–142 (2013).Article 

Google Scholar 
105.Brooks, M. E. et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R. J. 9, 378–400 (2017).Article 

Google Scholar 
106.Lozupone, C. & Knight, R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 71, 8228–8235 (2005).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
107.Lozupone, C. A., Hamady, M., Kelley, S. T. & Knight, R. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl. Environ. Microbiol. 73, 1576–1585 (2007).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
108.Anderson, M. J. Permutational Multivariate Analysis of Variance (PERMANOVA). https://doi.org/10.1002/9781118445112.stat07841. (2017)109.Anderson, M. J. & Walsh, D. C. I. PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: what null hypothesis are you testing? Ecol. Monogr. 83, 557–574 (2013).Article 

Google Scholar 
110.Li, H. et al. Pika population density is associated with the composition and diversity of gut microbiota. Front. Microbiol. 7, 1–9 (2016).
Google Scholar 
111.Weiss, S. et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome 5, 1–18 (2017).Article 

Google Scholar 
112.Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300 (1995).
Google Scholar 
113.Fackelmann, G. gfackelmann/human-encroachment-into-wildlife-gut-microbiomes: Release 1.0.0. https://doi.org/10.5281/zenodo.4725220. (2021) More

Study siteThe experiment was conducted at Madingley Wood, Cambridgeshire, UK (0◦3.2´E, 52◦12.9´N) during summer 2018. Madingley Wood is an established field site with an ongoing long-term study of the blue tit and great tit populations. During the autumn and winter birds are caught from feeding stations using mist nets and they are fitted with British Trust of Ornithology (BTO) ID rings. Since 2012, blue tits and great tits have been fitted with RFID tags (BTO Special Methods permit to HMR), which enables collecting data remotely about their foraging behavior and social relationships. The study site has 90 nest boxes that are monitored annually during the breeding season. In 2018 chicks (n = 325) fledged successfully from 45 nest boxes (blue tits = 21, great tits = 24) and they were all ringed and fitted with RFID tags when they were approximately 10 days old. Because new juvenile flocks were arriving at our study site throughout the summer, we also conducted several mist-netting and ringing sessions in July and August to maintain a high proportion of blue tits and great tits ringed and RFID tagged for the experiments (on average 89%, see below). The study protocol was approved by the Animal Users Committee at the Department of Zoology, University of Cambridge.Food itemsWe investigated birds’ foraging choices by offering them colored almond flakes at bird feeders that were distributed throughout the wood. Before beginning the experiments, we allowed the birds to become familiar with the food items by providing plain ‘control’ almonds (plain and not colored) in paired feeders (1.5 m apart) at three locations (approximately 170 m from each other). The feeders were surrounded with metal cages to exclude larger birds, and we placed plastic buckets under the feeders to collect any spilled almonds and minimize birds’ opportunities to forage from the ground. We introduced the feeders at the beginning of June when the nestlings had fledged and were beginning to forage independently, and continued to provide these plain almonds in between our learning experiments (Fig. 3c).In the learning experiments, almond flakes were dyed with non-toxic food dye (Classikool Concentrated Droplet Food Colouring). We used three different color pairs: green (Leaf Green) and red (Bright Red), purple (Lavender Purple) and blue (Royal Blue), and orange (Satsuma Orange) and yellow (Dandelion Yellow). Almond flakes were dyed by soaking them for approximately 20 min in a solution of 900 ml of water and 30 ml of food dye and then left to air dry for 48 h. In the avoidance learning experiments, we made half of the almond flakes unpalatable by soaking them for one hour in 67% solution of chloroquine phosphate, following previously established methods from avoidance learning studies with birds in captivity14,23,24,25. The food dye was added to the solution during the last 20 min.Red and green are common colors used by aposematic, or cryptic prey, respectively9. Therefore, we investigated whether blue tits and great tits had initial color biases towards these colors before starting the main experiment. Because we did not want the birds in our study population to have any experience of the colors before the main experiment, this pilot study was conducted in Newbury, which is 130 km from our main study site. Birds were simultaneously presented with two feeders containing red and green almonds (both palatable) for 30 min and the number of almonds of each color taken by blue tits or great tits was recorded using binoculars. The position of the feeders was switched after 15 min to control for any preferences for feeder location, and the test was repeated on 9 different days. We did not find any evidence that birds had initial color preferences (t-test: t = 0, df = 15.69, p = 1). For the other two learning experiments, we chose color pairs that were available as a food dye and as different from red and green in the visible spectrum as possible to avoid generalization across experiments. These color pairs (blue/purple and yellow/orange) had similar contrast ratios as green and red, based on their RGB values (measured from photographs, see Supplementary Information). Although avian and human vision is different, the discriminability of colors is likely to be similar51, and rapid avoidance learning in each experiment shows that all colors were easily distinguishable. This was the main requirement for testing social information use, and subtle differences in color pair discriminability should only introduce noise to our data but not influence our conclusions.Learning experiments with colored almondsWe conducted three avoidance learning experiments with different color pairs throughout the summer: red/green, blue/purple, and yellow/orange (unpalatable/palatable). In addition, we conducted a reversal-learning experiment with the blue/purple color pair by making both colors palatable after birds had acquired avoidance to blue almonds. Each experiment followed a similar protocol, in which birds were presented with colored almonds at the same three feeding stations where they were previously offered plain almonds. Each feeding station had two feeders, where one contained the palatable color and the other contained the unpalatable color (except in the reversal learning test when both colors were palatable). We switched the side of the feeders every day to make sure that birds learned to associate palatability with an almond color and not a feeder position. The feeders were filled at least once a day (or more often if necessary) to make sure that birds always had access to both colors. We continued each avoidance learning experiment until >90% of all recorded visits were to the feeder with palatable almonds, indicating that most birds in the population had learned to discriminate the colors. This took 7 days in the red/green experiment and 8 days in the other two color pairs (blue/purple and yellow/orange). The reversal learning experiment was finished after 9 days when 50% of the visits were to the previously unpalatable color (blue), indicating that most birds had reversed their learned avoidance towards it.Recording visits to feedersWe monitored visits to all feeders using RFID antennas and data loggers (Francis Scientific Instruments, Ltd) that scanned birds’ unique RFID tag codes when they landed on a perch attached to the feeder. During the learning experiments, each day we also recorded videos from all three feeding stations (using Go Pro Hero Action Camera and Canon Legria HF R66 Camcorder). From the videos, we monitored the proportion of blue tits and great tits that did not have RFID tags and were therefore not recorded when visiting the feeders. We calculated the estimated RFID tag coverage for each day of the experiments by watching at least 100 visits to the feeders from the videos (divided equally among the three feeding stations) and recording whether blue tits and great tits had an RFID tag or not. We realized that the number of untagged individuals was very high (approximately 50% of all visiting birds) when we started the experiment with the first color pair (red/green; see Supplementary Fig. 3). We, therefore, stopped the experiment after two days and caught birds from the feeding stations with mist nets to fit RFID tags to new individuals. To maintain a high number of individuals RFID tagged for the other color pairs, we conducted a mist netting session a day before starting each experiment, as well as 4–5 days after it. We always switched the feeders back to containing plain almonds during mist-netting sessions to ensure that this would not interfere with the learning experiments. Apart from the first two days of the red/green experiment, the RFID tag coverage was on average 89% throughout the experiments (varying between 80 and 95%, Supplementary Fig. 3).Birds were recorded every time that they visited the feeders, i.e., landed on the RFID antenna. However, it is possible that birds did not take the almond during every visit. To get an estimate of how often birds landed on the antenna without taking the almond, and whether this differed between palatable and unpalatable colors, we analyzed the visits to the feeders from the video recordings. We watched videos from the five first days of each experiment (i.e., different color pairs) and analyzed 60 visits to each color (divided approximately equally among the three feeding stations). We recorded whether the feeding event happened (birds ate the almond at the feeder or flew away with it) or whether birds left the feeder without sampling the almond. Because the number of visits to the unpalatable feeder was low during the last days of the avoidance learning experiments, we decided not to analyze avoidance learning videos after day five (but recorded visits from all days of the reversal learning experiment). We found that in avoidance learning experiments birds started to ‘reject’ unpalatable almonds after two days, i.e., they sometimes landed on the feeder but flew away without taking the almond (see Supplementary Fig. 4a). This change was not observed at palatable feeders where birds continued to consume almonds at a similar rate as at the beginning of the experiment (Supplementary Fig. 4a). In reversal learning, the proportion of visits that did not include a feeding event did not differ between purple and blue almonds: birds showed similar hesitation towards both colors at the beginning of the experiment, but this wariness decreased when the experiment progressed, with birds taking the almond during most of their visits (Supplementary Fig. 4b).Statistical analyses and model validationForaging choices in learning experimentsWe first analyzed how birds’ foraging choices changed during the learning experiments using generalized linear mixed-effects models with a binomial error distribution. The number of times an individual visited each feeder on each day of the experiment was used as a bounded response variable, and this was explained by species (blue tit/great tit), individuals’ age (juvenile/adult), and day of the experiment (continuous variable), as well as bird identity as a random effect. When analyzing avoidance learning, initial exploration of data suggested that results were similar across all three experiments, so we combined the experiments in the same model. To investigate whether learning curves differed between the species or age groups, the day of the experiment was included as a second-order polynomial term, and we started model selections with models that included a three-way interaction between species, age, and day2. Best-fitting models were selected based on Akaike’s information criterion (see Supplementary Tables 1 and 2).Social networkTo investigate if birds used social information in their foraging choices, we first constructed a social network of the bird population based on their visits to feeders outside of the learning experiments, i.e., when birds were presented with plain almonds (in total 92 days, see Supplementary Information for the robustness of analysis to exclusion of network data before or after the experiment). We used only these data as individuals were likely to vary in their hesitation to visit novel colored almonds. We used a Gaussian mixture model to detect the clusters of visits (‘gathering events’) at the feeders52 and then calculated association strengths between individuals based on how often they were observed in the same group (gambit of the group approach). These associations (network edges) were calculated using the simple ratio index, SRI35.$$frac{x}{x+{y}_{{mathrm{A}}}+{y}_{{mathrm{B}}}+{y}_{{{mathrm{AB}},}}}$$
(2)
where x is the number of samples where individuals A and B co-occurred in the same group, yA is the number of samples where only individual A was seen, yB is the number of samples where only individual B was seen, and yAB is the number of samples where both A and B were observed in the same sample but not together. Network associations, therefore, estimated the probability that two individuals were in the same group at a given time, with the values scaled between 0 (never observed in the same group) and 1 (always observed in the same group).Social information use during avoidance learning: model descriptionIf social avoidance learning was occurring, then the more birds observed negative responses of others feeding on the unpalatable feeder, the less likely they would be to choose the unpalatable feeder themselves. Thus, we expected the probability of bird j choosing the unpalatable option at time t to decrease with ({R}_{-,j}left(tright)) (the real number of negative feeding events observed by j prior to time t). Likewise, if appetitive social learning was occurring, then the more birds observed positive responses of others feeding on the palatable feeder, the more likely they would be to choose the palatable feeder themselves (rather than the unpalatable feeder). So, we also expected the probability of j choosing the unpalatable option at time t to decrease as ({R}_{+,j}left(tright),)(the real number of positive events observed by j prior to time t) increased.However, we could not test for an effect of ({R}_{-,j}left(tright)) and ({R}_{+,j}left(tright)) directly, since birds often ate the almond away from the feeder, and therefore the real number of observed feeding events could not be measured. Instead, we aimed to test for a pattern following the social network that is consistent with these social learning processes. We reasoned that the probability that one individual i, observes a specific feeding event by another individual j, was proportional to the network connection between them, aij (probability they are in the same feeding group at a given time). Therefore, in each avoidance learning experiment (i.e., different color pair), we calculated the expected number of negative feeding events observed, prior to each choice (occurring at time t) as$${O}_{-,i}left(tright)={sum }_{j}{N}_{-,j}left(tright){a}_{{ij}},$$
(3)
where ({N}_{-,j}left(tright)) was the number of times j had visited unpalatable almonds prior to time t (i ≠ j), and summation is across all birds in the network, and likewise for the expected number of positive feeding events:$${O}_{+,i}left(tright)={sum }_{j}{N}_{+,j}left(tright){a}_{{ij}},$$
(4)
where ({N}_{+,j}left(tright)) was the number of times j had visited palatable almonds prior to time t (i ≠ j).We analyzed whether the expected observations of positive and/or negative feeding events of others influenced the foraging choices in the avoidance learning experiments using generalized linear mixed-effects models with a binomial error distribution. We used each choice (i.e., visit a feeder) as a binary response variable (1 = unpalatable chosen, 0 = palatable chosen), with the probability that unpalatable feeder is chosen on feeding event E given by ({p}_{E}={p}_{-,{i}(E)}left({t}_{E}right)), where i(E) is the individual that fed during event E and ({t}_{E}) is the time at which event E occurred. We then modeled the probability of i choosing the unpalatable option at time t as:$${p}_{-,i}left(tright)={rm{logit}}left(alpha +{beta }_{{rm{asoc}}+}{N}_{+,i}left(tright)+{beta }_{{rm{asoc}}-}{N}_{-,i}left(tright)+{beta }_{{rm{soc}}+}{O}_{+,i}left(tright)+{beta }_{s{rm{oc}}-}{O}_{-,i}left(tright)+{{{{rm{B}}}}}_{i}right),$$
(5)
where ({N}_{+,i}left(tright)) is the number of times a choosing individual had visited the palatable feeder (positive personal information), ({N}_{-,i}left(tright)) is the number of times a choosing individual had visited the unpalatable feeder (negative personal information), ({O}_{+,i}left(tright)) is the expected number of observed positive (positive social information) and ({O}_{-,i}left(tright)) observed negative feeding events (negative social information). Bird identity was included as a random effect, ({{rm{{B}}}}_{i}) (age and species were later added as variables, see below). Parameters ({beta }_{{rm{asoc}}+}) and ({beta }_{{rm{asoc}}-}) are the effects of asocial learning about the palatable and unpalatable foods, ({beta }_{{rm{soc}}+}) is the effect of social learning about the palatable food, and ({beta }_{{rm{soc}}-})is the effect of social avoidance learning about the unpalatable food. Estimation of these parameters, with associated Wald tests and confidence intervals, allowed us to make inferences about which effects were operating and the size of these effects. To aid model fitting we standardized all predictor variables and then back-transformed the effects to the original scale (see Supplementary Tables 3–5 for the model outputs). To assess the importance of asocial and social effects, we also ran separate models that excluded either asocial or social parameters and compared them to the initial model in Eq. (5) using Akaike’s information criterion (see Supplementary Table 6). However, in most cases, this reduced model fit significantly, and we, therefore, kept all parameters in the final models.Our approach took ({O}_{-,{i}}left(tright)) as a measure of ({R}_{-,j}left(tright)), and ({O}_{+,{i}}left(tright)) as a measure of ({R}_{+,j}left(tright))-, which we termed the ‘expected’ number of observations of each type. Strictly speaking, ({O}_{-,{i}}left(tright)) and ({O}_{+,{i}}left(tright)) were upper limits on the expected number of observations, assuming that birds observed all feeding events in the groups in which they were present, whereas only an unknown proportion of such events (({p}_{o})) was observed. Therefore, the real expected number of negative/positive observations would be (Eleft({R}_{-,j}left(tright)right)={p}_{o}{O}_{-,{i}}left(tright)) and (Eleft({R}_{+,j}left(tright)right)={p}_{o}{O}_{+,{i}}left(tright)) respectively. Thus, the coefficient, ({beta }_{{rm{soc}}-}), for the effect of ({O}_{-,{i}}left(tright)) could be interpreted as ({beta }_{s{rm{oc}}-}={p}_{o}acute{{beta }_{{rm{soc}}-}}) where (acute{{beta }_{s{rm{oc}}-}}) is the effect per observation. Note that since ({p}_{o}le 1), and(,{beta }_{{rm{soc}}-}=acute{{beta }_{s{rm{oc}}-}}{p}_{o}), ({beta }_{s{rm{oc}}-}) is more likely to underestimate than overestimate the effect per observation of a negative feeding event. An analogous argument applies to the coefficient, ({beta }_{{rm{soc}}+}), for the effect of ({O}_{+,{i}}left(tright)).Social information use during avoidance learning: extension to test for species effectsAfter fitting the initial model shown in Eq. (5), we further broke down the model to test whether individuals were more likely to learn socially by observing conspecifics than heterospecifics. This was done by splitting the expected number of observed positive and negative feeding events to observations of conspecifics (({O}_{+{rm{C}},i}left(tright)), ({O}_{-{rm{C}},i}left(tright))) and heterospecifics (({O}_{+{rm{H}},i}left(tright)), ({O}_{-{rm{H}},i}left(tright))), and including these in the model as separate explanatory variables thus:$${p}_{-,i}(t)={rm{logit}}left(begin{array}{c}alpha +{beta }_{{rm{asoc}}+}{N}_{+,i}(t)+{beta }_{{rm{asoc}}-}{N}_{-,i}(t)\ +{beta }_{{rm{soc}},+{rm{H}}}{O}_{+{rm{H}},i}(t)+{beta }_{{rm{soc}},-{rm{H}}}{O}_{-{rm{H}},i}(t)\ +{beta }_{{rm{soc}},+{rm{C}}}{O}_{+{rm{C}},i}(t)+{beta }_{{rm{soc}},-{rm{C}}}{O}_{-{rm{C}},i}(t) \ +{{{{rm{B}}}}}_{i}end{array}right)$$
(6a)
with ({beta }_{{rm{soc}},-{rm{H}}}) and ({beta }_{{rm{soc}},-{rm{C}}}) giving the effect of a negative observation of a heterospecific and conspecific, respectively, whereas ({beta }_{{rm{soc}},+{rm{H}}}) and ({beta }_{{rm{soc}},+{rm{C}}}) give the effect of positive observation of a heterospecific and conspecific, respectively. In general –/+ subscripts refer to negative/positive feeding events and C/H subscripts to feeding events by conspecifics/heterospecifics. By re-parameterizing the model thus:$${p}_{-,i}left(tright)={rm{logit}}left(begin{array}{c}alpha +{beta }_{{rm{asoc}}+}{N}_{+,i}left(tright)+{beta }_{{rm{asoc}}-}{N}_{-,i}left(tright) \ +{beta }_{{rm{soc}},{rm{H}}+}{O}_{+,i}left(tright)+{beta }_{{rm{soc}},{rm{H}}-}{O}_{-,i}left(tright)\ +left({beta }_{{rm{soc}},{rm{C}}+}-{beta }_{{rm{soc}},{rm{H}}+}right){O}_{+{rm{C}},i}left(tright)+left({beta }_{{rm{soc}},{rm{C}}-}-{beta }_{{rm{soc}},{rm{H}}-}right){O}_{-{rm{C}},i}left(tright)\ +{{{{rm{B}}}}}_{i}end{array}right)$$
(6b)
we were able to test for a difference between observations of negative feeds by conspecifics and heterospecifics (left({beta }_{{rm{soc}},{rm{C}}-}-{beta }_{{rm{soc}},{rm{H}}-}right)) and between observations of positive feeds by conspecifics and heterospecifics (left({beta }_{{rm{soc}},{rm{C}}+}-{beta }_{{rm{soc}},{rm{H}}+}right)).For all experiments there was no evidence for a difference between ({beta }_{{rm{soc}},-{rm{H}}}) and ({beta }_{{rm{soc}},-{rm{C}}}) (yellow/orange: Z = 0.803, p = 0.42; red/green: Z = 0.065, p = 0.95; blue/purple: Z = 1.113, p = 0.27). However, there was some evidence of a difference between ({beta }_{{rm{soc}},+{rm{H}}}) and ({beta }_{{rm{soc}},+{rm{C}}}) in two of the three experiments (yellow/orange: Z = 1.359, p = 0.17; red/green: Z = 1.417, p = 0.16; blue/purple: Z = 0.729, p = 0.47). Consequently, we reduced the model down to:$${p}_{-,i}left(tright)={rm{logit}}left(begin{array}{c}alpha +{beta }_{{rm{asoc}}+}{N}_{+,i}left(tright)+{beta }_{{rm{asoc}}-}{N}_{-,i}left(tright)\ +{beta }_{{rm{soc}},-}{O}_{-,i}left(tright)+{beta }_{{rm{soc}},+{rm{H}}}{O}_{+{rm{H}},i}left(tright)+{beta }_{{rm{soc}},+{rm{C}}}{O}_{+{rm{C}},i}left(tright)\ +{{{{rm{B}}}}}_{i}end{array}right)$$
(7)
for further analysis, i.e., with different effects for observations of conspecific/heterospecific positive feeds, but not of negative feeds. We did this for all color combinations (including blue/purple) to allow comparison across experiments (see Table 1). The R code used to run these models can be found in Supplementary data53 in ‘GLMM models Orange Yellow final.r’.Social information use during avoidance learning: simulations to test for a network effectNext, we tested whether the social effects we detected followed the social network. When using a network-based diffusion analysis (NBDA43), researchers can compare a network model with one in which the network has homogeneous connections among all individuals, but we found this to be unreliable for our model. Instead, we used a simulation approach to generate a null distribution for the null hypothesis of homogeneous social effects, taking the size of the social effects from the fitted models. We ran 1000 simulations (using the same procedure described above) for all social effects that were found to be significant in each avoidance learning model (each color pair; see Table 1). The total number of expected observations was kept equal, but we homogenized the observation effect across all birds by replacing the probability of bird i observing a feed by bird j, previously ({a}_{{ij}}), with ({sum }_{i}{a}_{{ij}}/n), where n is the number of birds in the experiment, (i.e., all birds had the same probability of observing each feeding event). The model was fitted to the simulated data each time to extract the Z value (Wald test statistic) of the social effect of interest. The distribution of these values was then used as a null distribution to test whether our observed social effect differed from the effects that did not follow the social network. To this end, we calculated the proportion of simulations that yielded a Z value as extreme or more extreme than that observed (judged by distance in either direction from the mean of the null distribution). The R code used to run these simulations can be found in Supplementary data53 in ‘Simulations to test if network effects follow network Orange Yellow.r’.Social information use during avoidance learning: extension to test for age effectsWe then aimed to test whether each of the three social effects detected differed based on the age class of the observed individual (adult versus juveniles). We, therefore, split the negative expected observations ({O}_{-,i}left(tright)) into the expected observations of adults ({O}_{-{rm{A}},i}left(tright)) and juveniles ({O}_{-{rm{J}},i}left(tright)), each with its associated coefficient in the model ({beta }_{{rm{soc}},-{rm{A}}}) and ({beta }_{{rm{soc}},-{rm{J}}}). Likewise, we split positive observations of conspecifics as ({O}_{+{rm{CA}},i}left(tright)) and ({O}_{+{rm{CJ}},i}left(tright)) and positive observations of heterospecifics as ({O}_{+{rm{HA}},i}left(tright)) and ({O}_{+{rm{HJ}},i}left(tright)) to give the model:$${p}_{-,i}left(tright)={rm{logit}}left(begin{array}{c}alpha +{beta }_{{rm{asoc}}+}{N}_{+,i}left(tright)+{beta }_{{rm{asoc}}-}{N}_{-,i}left(tright)\ +{beta }_{{rm{soc}},-{rm{A}}}{O}_{-{rm{A}},i}left(tright)+{beta }_{{rm{soc}},-{rm{J}}}{O}_{-{rm{J}},i}left(tright)\ +{beta }_{{rm{soc}},+{rm{HA}}}{O}_{+{rm{HA}},i}left(tright)+{beta }_{{rm{soc}},+{rm{HJ}}}{O}_{+{rm{HJ}},i}left(tright)\ +{beta }_{{rm{soc}},+{rm{CA}}}{O}_{+{rm{CA}},i}left(tright)+{beta }_{{rm{soc}},+{rm{CJ}}}{O}_{+{rm{CJ}},i}left(tright)\ +{{{{rm{B}}}}}_{i}end{array}right)$$
(8a)
As before, –/+ subscripts refer to negative/positive feeding events, C/H subscripts to feeding events by conspecifics/heterospecifics, and A/J subscripts to feeding events by adults/juveniles. We also fitted a re-parameterized version allowing us to test for a difference between expected observations of adults and observations of juveniles for each of the three social effects:$${p}_{-,i}left(tright)={rm{logit}}left(begin{array}{c}alpha +{beta }_{{rm{asoc}}+}{N}_{+,i}left(tright)+{beta }_{{rm{asoc}}-}{N}_{-,i}left(tright)\ +{beta }_{{rm{soc}},-{rm{J}}}{O}_{-,i}left(tright)+left({{beta }_{{rm{soc}},-{rm{A}}}-beta }_{{rm{soc}},-{rm{J}}}right){O}_{-{rm{A}},i}left(tright)\ +{beta }_{{rm{soc}},+{rm{H}}}{O}_{+{rm{H}},i}left(tright)+left({{beta }_{{rm{soc}},+{rm{HA}}}-beta }_{{rm{soc}},+{rm{HJ}}}right){O}_{+{rm{JA}},i}left(tright)\ +{beta }_{{rm{soc}},+{rm{C}}}{O}_{+{rm{C}},i}left(tright)+left({{beta }_{{rm{soc}},+{rm{CA}}}-beta }_{{rm{soc}},+{rm{CJ}}}right){O}_{+{rm{CA}},i}left(tright)\ +{{{{rm{B}}}}}_{i}end{array}right)$$
(8b)
The R code used to run these models can be found in Supplementary data53 in ‘GLMM models Orange Yellow final.r’. The main results of each model are presented in Table 2 and full model outputs in Supplementary Tables 3–5.Social information use during reversal learningTo investigate social information use during reversal learning, we used the order of acquisition diffusion analysis (OADA), a variant of NBDA43, which explores the order in which individuals acquire a behavioral trait44. The rate of social transmission between two individuals is assumed to be linearly proportional to their network connection, and the spread of trait acquisition is therefore predicted to follow the network patterns if individuals are using social information. We used NBDA to investigate whether the order of individuals’ first visit to the previously unpalatable blue almonds (mimics) followed the network. We fitted several different models that included (i) only asocial learning, (ii) social transmission of information following a homogeneous network (equal associations among all individuals), or (iii) social transmission of information following our observed network. Models that included social transmission were further divided into models with equal or different transmission rates from adults and juveniles, and from conspecifics and heterospecifics, by constructing separate networks for each adult/juvenile and conspecific/heterospecific combination. To investigate whether asocial or social learning rates differed between blue tits and great tits, we included species as an individual-level variable. We then compared different social transmission models that assumed that species differed in both asocial and social learning rates, only in asocial or only in social learning rates, or that they did not differ in either (see Table 3). The best-supported model was selected using a model-averaging approach with Akaike’s information criterion corrected for small sample sizes. All analyses were conducted with the software R.3.6.154, using lme455, asnipe56, and NBDA57 packages.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

1.IPCC Climate Change 2014: Synthesis report. In Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (eds Core Writing Team et al.) (IPCC, 2014).
Google Scholar 
2.Keogan, K. et al. Global phenological insensitivity to shifting ocean temperatures among seabirds. Nat. Clim. Chang. 8, 313–318 (2018).ADS 
Article 

Google Scholar 
3.Chambers, L. E. et al. Phenological changes in the southern hemisphere. PLoS ONE 8, E75514 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
4.Thackeray, S. J. et al. Phenological sensitivity to climate across taxa and trophic levels. Nature 535, 241–245 (2016).ADS 
CAS 
PubMed 
Article 

Google Scholar 
5.Forrest, J. & Miller-Rushing, A. J. Toward a synthetic understanding of the role of phenology in ecology and evolution. Philos. Trans. R. Soc. Lond. Series B Biol. Sci. 365, 3101–3112 (2010).Article 

Google Scholar 
6.Møller, A. P., Rubolini, D. & Lehikoinen, E. Populations of migratory bird species that did not show a phenological response to climate change are declining. Proc. Natl. Acad. Sci. 105, 16195–16200 (2008).ADS 
PubMed 
PubMed Central 
Article 

Google Scholar 
7.Hipfner, J. M. Matches and mismatches: Ocean climate, prey phenology and breeding success in a zooplanktivorous seabird. Mar. Ecol. Prog. Ser. 368, 295–304 (2008).ADS 
Article 

Google Scholar 
8.Moyes, K. et al. Advancing breeding phenology in response to environmental change in a wild red deer population. Glob. Change Biol. 17, 2455–2469 (2011).ADS 
Article 

Google Scholar 
9.Poloczanska, E. S. et al. Global imprint of climate change on marine life. Nat. Clim. Chang. 3, 919–925 (2013).ADS 
Article 

Google Scholar 
10.Hauser, D. D. et al. Decadal shifts in autumn migration timing by Pacific Arctic beluga whales are related to delayed annual sea ice formation. Glob. Change Biol. 23, 2206–2217 (2017).ADS 
Article 

Google Scholar 
11.Beltran, R. S., Kirkham, A. L., Breed, G. A., Ward Testa, J. & Burns, J. M. Reproductive success delays moult phenology in a polar mammal. Sci. Rep. 9, 1–12 (2019).CAS 
Article 

Google Scholar 
12.Oosthuizen, W. C. et al. Individual heterogeneity in life-history trade-offs with age at first reproduction in capital breeding elephant seals. Popul. Ecol. 61, 421–435 (2019).Article 

Google Scholar 
13.Ling, J. K. The skin and hair of the southern elephant seal, Mirounga leonina (Linn). IV. Annual cycle of pelage follicle activity and moult. Aust. J. Zool. 60, 259–271 (2012).Article 

Google Scholar 
14.Johanos, T. C., Becker, B. L. & Ragen, T. J. Annual reproductive cycle of the female Hawaiian monk seal (Monachus schauinslandi). Mar. Mamm. Sci. 10, 13–30 (1994).Article 

Google Scholar 
15.Beltran, R. S., Burns, J. M. & Breed, G. A. Convergence of biannual moulting strategies across birds and mammals. Proc. R. Soc. B Biol. Sci. 285, 20180318 (2018).Article 

Google Scholar 
16.Boyd, I., Arnbom, T. & Fedak, M. Water flux, body composition, and metabolic rate during moult in female southern elephant seals (Mirounga leonina). Physiol. Zool. 66, 43–60 (1993).Article 

Google Scholar 
17.Postma, M., Bester, M. N. & De Bruyn, P. J. N. Spatial variation in female southern elephant seal mass change assessed by an accurate non-invasive photogrammetry method. Antarct. Sci. 25, 731–740 (2013).ADS 
Article 

Google Scholar 
18.Paterson, W. et al. Seals like it hot: Changes in surface temperature of harbour seals (Phoca vitulina) from late pregnancy to moult. J. Therm. Biol. 37, 454–461 (2012).Article 

Google Scholar 
19.Slip, D. J. & Burton, H. R. Population status and seasonal haulout patterns of the southern elephant seal (Mirounga leonina) at Heard Island. Antarct. Sci. 11, 38–47 (1999).ADS 
Article 

Google Scholar 
20.Kirkman, S. et al. Variation in the timing of moult in southern elephant seals at Marion Island. S. Afr. J. Wildlife Res. 33, 79–84 (2003).
Google Scholar 
21.Ling, J. & Bryden, M. Southern elephant seal Mirounga leonina (Linnaeus), 1758. In Handbook of Marine Mammals Vol. 2 (eds Ridgway, S. H. & Harrison, R. J.) 297–327 (Academic Press, 1981).
Google Scholar 
22.Constable, A. J. et al. Climate change and southern ocean ecosystems I: how changes in physical habitats directly affect marine biota. Glob. Change Biol. 20, 3004–3025 (2014).ADS 
Article 

Google Scholar 
23.Rogers, A. et al. Antarctic futures: An assessment of climate-driven changes in ecosystem structure, function, and service provisioning in the Southern Ocean. Ann. Rev. Mar. Sci. 12, 87–120 (2020).CAS 
PubMed 
Article 

Google Scholar 
24.Hindell, M. A. & Burton, H. R. Seasonal haul-out patterns of the southern elephant seal (Mirounga leonina L.), at Macquarie Island. J. Mammal. 69, 81–88 (1988).Article 

Google Scholar 
25.Postma, M., Bester, M. N. & De Bruyn, P. J. N. Age-related reproductive variation in a wild marine mammal population. Polar Biol. 36, 719–729 (2013).Article 

Google Scholar 
26.Le Boeuf, B. J. & Laws, R. M. Elephant seals: An introduction to the genus. In Elephant Seals: Population Ecology, Behavior, and Physiology (eds Le Beouf, B. J. & Laws, R. M.) 1–28 (University Of California Press, 1994).Chapter 

Google Scholar 
27.De Bruyn, P. J. N. et al. Sex at sea: Alternative mating system in an extremely polygynous mammal. Anim. Behav. 82, 445–451 (2011).Article 

Google Scholar 
28.Oosthuizen, W. C., Pradel, R., Bester, M. N. & De Bruyn, P. J. N. Making use of multiple surveys: Estimating breeding probability using a multievent-robust design capture-recapture model. Ecol. Evolut. 9, 836–848 (2019).Article 

Google Scholar 
29.Bester, M. N. et al. The marine mammal programme at the Prince Edward Islands: 38 years of research. Afr. J. Mar. Sci. 33, 511–521 (2011).Article 

Google Scholar 
30.Pistorius, P. A., De Bruyn, P. J. N. & Bester, M. N. Population dynamics of southern elephant seals: A synthesis of three decades of demographic research at Marion Island. Afr. J. Mar. Sci. 33, 523–534 (2011).Article 

Google Scholar 
31.Pistorius, P. A., Bester, M. N. & Kirkman, S. P. Survivorship of a declining population of southern elephant seals, Mirounga leonina, in relation to age, sex and cohort. Oecologia 121, 201–211 (1999).ADS 
CAS 
PubMed 
Article 

Google Scholar 
32.Oosthuizen, W. C., Bester, M. N., Altwegg, R., McIntyre, T. & De Bruyn, P. J. N. Decomposing the variance in southern elephant seal weaning mass: Partitioning environmental signals and maternal effects. Ecosphere 6, 139 (2015).Article 

Google Scholar 
33.Daniel, R. G., Jemison, L. A., Pendleton, G. W. & Crowley, S. M. Molting phenology of harbor seals on Tugidak Island, Alaska. Mar. Mamm. Sci. 19, 128–140 (2003).Article 

Google Scholar 
34.Chaise, L. L. et al. Local weather and body condition influence habitat use and movements on land of molting female southern elephant seals (Mirounga leonina). Ecol. Evol. 8, 6081–6090 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
35.Dingemanse, N. J. & Dochtermann, N. A. Quantifying individual variation in behaviour: Mixed-effect modelling approaches. J. Anim. Ecol. 82, 39–54 (2013).PubMed 
Article 

Google Scholar 
36.Burnham, K. P., Anderson, D. R. & Huyvaert, K. P. AIC model selection and multimodel inference in behavioral ecology: Some background, observations, and comparisons. Behav. Ecol. Sociobiol. 65, 23–35 (2011).Article 

Google Scholar 
37.Zuur, A., Ieno, E. N., Walker, N., Saveliev, A. A. & Smith, G. M. Mixed Effects Models and Extensions in Ecology with R (Springer Science & Business Media, 2009).MATH 
Book 

Google Scholar 
38.Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).Article 

Google Scholar 
39.R Core Team. R: A language and environment for statistical computing. (R Foundation for Statistical Computing, Vienna, Austria, 2020). https://www.R-project.org/.40.Nakagawa, S. & Schielzeth, H. Repeatability for Gaussian and non-Gaussian data: A practical guide for biologists. Biol. Rev. 85, 935–956 (2010).PubMed 

Google Scholar 
41.Nakagawa, S. & Schielzeth, H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 4, 133–142 (2013).Article 

Google Scholar 
42.Grosbois, V. et al. Assessing the impact of climate variation on survival in vertebrate populations. Biol. Rev. 83, 357–399 (2008).CAS 
PubMed 
Article 

Google Scholar 
43.Forcada, J., Trathan, P. N. & Murphy, E. J. Life history buffering in Antarctic mammals and birds against changing patterns of climate and environmental variation. Glob. Change Biol. 14, 2473–2488 (2008).ADS 
Article 

Google Scholar 
44.Bradshaw, C. J., Hindell, M. A., Sumner, M. D. & Michael, K. J. Loyalty pays: Potential life-history consequences of fidelity to marine foraging regions by southern elephant seals. Anim. Behav. 68, 1349–1360 (2004).Article 

Google Scholar 
45.Cotté, C., Park, Y.-H., Guinet, C. & Bost, C.-A. Movements of foraging king penguins through marine mesoscale eddies. Proc. R. Soc. B Biol. Sci. 274, 2385–2391 (2007).Article 

Google Scholar 
46.Barbraud, C. & Weimerskirch, H. Antarctic birds breed later in response to climate change. Proc. Natl. Acad. Sci. 103, 6248–6251 (2006).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
47.Hindell, M. A. et al. Long-term breeding phenology shift in royal penguins. Ecol. Evol. 2, 1563–1571 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
48.Zimova, M. et al. Function and underlying mechanisms of seasonal colour moulting in mammals and birds: What keeps them changing in a warming world?. Biol. Rev. 93, 1478–1498 (2018).PubMed 
Article 

Google Scholar 
49.Thompson, P. & Rothery, P. Age and sex differences in the timing of moult in the common seal, Phoca vitulina. J. Zool. 212, 597–603 (1987).Article 

Google Scholar 
50.Cronin, M., Gregory, S. & Rogan, E. Moulting phenology of the harbour seal in south-west Ireland. J. Mar. Biol. Assoc. UK 94, 1079–1086 (2014).Article 

Google Scholar 
51.Kirkman, S. et al. Variation in the timing of the breeding haulout of female southern elephant seals at Marion Island. Aust. J. Zool. 52, 379–388 (2004).Article 

Google Scholar 
52.Hindell, M. A., Slip, D. J. & Burton, H. R. Body mass loss of moulting female southern elephant seals, Mirounga leonina, at Macquarie Island. Polar Biol. 14, 275–278 (1994).Article 

Google Scholar 
53.Carlini, A., Marquez, M., Daneri, G. & Poljak, S. Mass changes during their annual cycle in females of southern elephant seals at King George Island. Polar Biol. 21, 234–239 (1999).Article 

Google Scholar 
54.Worthy, G. A. J., Morris, P. A., Costa, D. P. & Le Boeuf, B. J. Moult energetics of the northern elephant seal (Mirounga angustirostris). J. Zool. 227, 257–265 (1992).Article 

Google Scholar 
55.Arnbom, T., Fedak, M. & Rothery, P. Offspring sex ratio in relation to female size in southern elephant seals, Mirounga leonina. Behav. Ecol. Sociobiol. 35, 373–378 (1994).Article 

Google Scholar 
56.Proffitt, K. M., Garrott, R. A., Rotella, J. J. & Wheatley, K. E. Environmental and senescent related variations in Weddell seal body mass: Implications for age-specific reproductive performance. Oikos 116, 1683–1690 (2007).Article 

Google Scholar 
57.Bérubé, C. H., Festa-Bianchet, M. & Jorgenson, J. T. Individual differences, longevity, and reproductive senescence in bighorn ewes. Ecology 80, 2555–2565 (1999).Article 

Google Scholar 
58.Oosthuizen, W. C., Péron, G., Pradel, R., Bester, M. N. & De Bruyn, P. J. N. Positive early-late life-history trait correlations in elephant seals. Ecology 102, e03288. https://doi.org/10.1002/ecy.3288 (2021).59.Kirby, V. L. & Ortiz, C. L. Hormones and fuel regulation in fasting elephant seals. In Elephant Seals: Population Ecology, Behavior, and Physiology (eds Le Beouf, B. J. & Laws, R. M.) 374–386 (University of California Press, 1994).Chapter 

Google Scholar 
60.Yochem, P. K. & Stewart, B. S. Hair and fur. In Encyclopedia of Marine Mammals 3rd edn (eds Würsig, B. et al.) 447–448 (Academic Press, 2018).Chapter 

Google Scholar 
61.Ashwell-Erickson, S., Fay, F., Elsner, R. & Wartzok, D. Metabolic correlates of molting and regeneration of pelage in Alaskan harbour and spotted seals (Phoca vitulina and Phoca largha). Can. J. Zool. 64, 1086–1092 (1986).CAS 
Article 

Google Scholar 
62.Laws, R. M. & Sinha, A. A. Reproduction. In Antarctic Seals: Research Methods and Techniques (ed. Laws, R. M.) 228–267 (Cambridge University Press, 1993).Chapter 

Google Scholar 
63.Smout, S., King, R. & Pomeroy, P. Environment-sensitive mass changes influence breeding frequency in a capital breeding marine top predator. J. Anim. Ecol. 89, 384–396 (2020).PubMed 
Article 

Google Scholar 
64.Clutton-Brock, T. & Sheldon, B. C. Individuals and populations: the role of long-term, individual-based studies of animals in ecology and evolutionary biology. Trends Ecol. Evol. 25, 562–573 (2010).PubMed 
Article 

Google Scholar 
65.Cordes, L. S. & Thompson, P. M. Variation in breeding phenology provides insights into drivers of long-term population change in harbour seals. Proc. R. Soc. B: Biol. Sci. 280, 20130847 (2013).Article 

Google Scholar 
66.Barbraud, C. et al. Contrasted demographic responses facing future climate change in Southern Ocean seabirds. J. Anim. Ecol. 80, 89–100 (2011).PubMed 
Article 

Google Scholar 
67.Atkinson, S., DeMaster, D. P. & Calkins, D. G. Anthropogenic causes of the western Steller sea lion Eumetopias jubatus population decline and their threat to recovery. Mamm. Rev. 38, 1–18 (2008).Article 

Google Scholar 
68.Desprez, M., McMahon, C. R., Hindell, M. A., Harcourt, R. & Gimenez, O. Known unknowns in an imperfect world: incorporating uncertainty in recruitment estimates using multi-event capture-recapture models. Ecol. Evol. 3, 4658–4668 (2013).PubMed 
PubMed Central 
Article 

Google Scholar 
69.Griffen, B. D. Reproductive skipping as an optimal life history strategy in the southern elephant seal, Mirounga leonina. Ecol. Evol. 8, 9158–9170 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
70.Rohwer, S., Viggiano, A. & Marzluff, J. M. Reciprocal tradeoffs between moult and breeding in albatrosses. Condor 113, 61–70 (2011).Article 

Google Scholar 
71.Grissot, A., Graham, I. M., Quinn, L., Bråthen, V. S. & Thompson, P. M. Breeding status influences timing but not duration of moult in the northern fulmar Fulmarus glacialis. Ibis 162, 446–459 (2019).Article 

Google Scholar 
72.Condit, R., Beltran, R.S., Robinson, P.W., Crocker, D.E. & Costa, D.P. Birth timing after the long feeding migration in elephant seals. https://doi.org/10.1101/2020.10.02.324210 (in review). More

1.Cámara, A. & Santero-Sánchez, R. Economic, social, and environmental impact of a sustainable fishereis model in Spain. Sustainability 11, 6311 (2019).Article 

Google Scholar 
2.Department of Fisheries Malaysia. Annual Fisheries statistics 2010–2019. https://www.dof.gov.my/index.php/pages/view/82 (2020).3.Department of Fisheries, Thailand. The annual marine fisheries statistics (2010–2019) based on the sample survey. https://elibonline.fisheries.go.th/elib/cgi-bin/opacexe.exe?op=dsp&bid=10498&lang=0&db=Main&pat=&cat=sub&skin=s&lpp=20&catop=edit&scid=zzz (2020).4.Jha, S., Deepti, V., Ravali, V. & Sujatha, K. Studies on some aspects of biology of Uranoscopus cognatus Cantor, 1849 (Pisces: Uranoscopidae) off Visakhapatnam, central eastern coast of India. Indian J. Mar. Sci. 48, 85–92 (2019).
Google Scholar 
5.Clark, M. R. et al. The impacts of deep-sea fisheries on benthic communities: A review. ICES J. Mar. Sci. 73(1), 51–69 (2016).Article 

Google Scholar 
6.Van Denderen, P. D. et al. Evaluating impacts of bottom trawling and hypoxia on benthic communities at the local, habitat, and regional scale using a modelling approach. ICES J. Mar. Sci. 77(1), 578–589 (2019).
Google Scholar 
7.Erdoğan Sağlam, N. & Sağlam, C. Population parameters of stargazer (Uranoscopus scaber Linnaeus, 1758) in the southeastern Black Sea region during the 2011–2012 fishing season. J. Appl. Ichthyol. 29, 1313–1317 (2013).Article 

Google Scholar 
8.Matsunuma, M. et al. Fishes of Terengganu: East Coast of Malay Peninsula, Malaysia (National Museum of Nature and Science, 2011).
Google Scholar 
9.Vilasri, V. Family Uranoscopidae. In Fishes of Southern Taiwan (eds Koeda, K. & Ho, H. S.) 1097–1105 (National Museum of Marine Biology & Aquarium, 2019).
Google Scholar 
10.Starks, E. C. The Osteology and Relationships of the Uranoscopoid Fishes (Stanford University Press, 1923).
Google Scholar 
11.Pietsch, T. W. Phylogenetic relationships of trachinoid fishes of the family Uranoscopidae. Copeia 1989, 253–303 (1989).Article 

Google Scholar 
12.Kishimoto, H. Uranoscopidae. In FAO Species Identification Guide for Fisheries Purposes (eds Carpenter, K. E. & Niem, V. H.) 3519–3531 (FAO, 2001).
Google Scholar 
13.Randall, J. E. & Arnold, R. J. Uranoscopus rosette, a new species of stargazer (Uranoscopidae: Trachinoidei) from the Red Sea. Aqua. Int. J. Ichthyol. 18, 209–219 (2012).
Google Scholar 
14.Jung-chen, H. & Hin-Kiu, M. Stargazers (Uranoscopidae) have exceptionally more bile. Kuroshio Sci. 9–1, 17–26 (2015).
Google Scholar 
15.Vilasri, V. Comparative anatomy and phylogenetic systematics of the family Uranoscopidae (Actinopterygii: Perciformes). Mem. Fac. Fish. Hokkaido Univ. 55, 1–106 (2013).
Google Scholar 
16.Fricke, R., Eschmeyer, W. N. & Van der Laan, R. (eds). Eschmeyer’s catalog of fishes: genera, species, references. http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp (2020).17.Froese, R. & Pauly, D. Uranoscopidae. Fishbase https://www.fishbase.se/Summary/FamilySummary.php?ID=378 (2019).18.Fricke, R. Two new species of stargazers of the genus Uranoscopus (Teleostei: Uranoscopidae) from the western Pacific Ocean. Zootaxa 4476, 157–167 (2018).PubMed 
Article 

Google Scholar 
19.Fricke, R., Jawad, L. A., Al-Kharusi, L. H. & Al-Mamry, J. M. New record and redescription of Uranoscopus crassiceps Alcock, 1890 (Uranoscopidae) From Oman, Arabian Sea, Northwestern Indian Ocean, based on adult specimens. Cybium 37, 143–147 (2013).
Google Scholar 
20.Department of Fisheries Malaysia. Valid Local Name of Malaysian Marine Fishes (Department of Fisheries Malaysia, 2009).
Google Scholar 
21.Spalding, M. D. et al. Marine ecoregions of the world: A bioregionalization of coastal and shelf areas. Bioscience 57, 573–583 (2007).Article 

Google Scholar 
22.Voris, H. K. Maps of Pleistocene sea levels in Southeast Asia: Shorelines, river systems and time durations. J. Biogeogr. 27, 1153–1167 (2000).Article 

Google Scholar 
23.Rohfritsch, A. & Borsa, P. Genetic structure of Indian scad mackerel Decapterus russelli: Pleistocene vicariance and secondary contact in the Central Indo-West Pacific Seas. Heredity 95, 315 (2005).CAS 
PubMed 
Article 

Google Scholar 
24.Lohman, D. J. et al. Biogeography of the Indo-Australian archipelago. Annu. Rev. Ecol. Evol. Syst. 42, 205–226 (2011).Article 

Google Scholar 
25.Crandall, E. D. et al. The molecular biogeography of the Indo-Pacific: Testing hypotheses with multispecies genetic patterns. Glob. Ecol. Biogeogr. 28, 943–960 (2019).Article 

Google Scholar 
26.Reece, J. S., Bowen, B. W., Joshi, K., Goz, V. & Larson, A. Phylogeography of two moray eels indicates high dispersal throughout the Indo-Pacific. J. Hered. 101, 391–402 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
27.Akib, N. A. M. et al. High connectivity in Rastrelliger kanagurta: influence of historical signatures and migratory behaviour inferred from mtDNA cytochrome b. PLoS ONE 10, e0119749 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
28.Jamaludin, N. A. et al. Phylogeography of the Japanese scad, Decapterus maruadsi (Teleostei; Carangidae) across the Central Indo-West Pacific: evidence of strong regional structure and cryptic diversity. Mitochondrial DNA A 2, 1–13 (2020).
Google Scholar 
29.Gaither, M. R., Toonen, R. J., Robertson, D. R., Planes, S. & Bowen, B. W. Genetic evaluation of marine biogeographical barriers: perspectives from two widespread Indo-Pacific snappers (Lutjanus kasmira and Lutjanus fulvus). J. Biogeogr. 37, 133–147 (2010).Article 

Google Scholar 
30.Gaither, M. R. et al. Phylogeography of the reef fish Cephalopholis argus (Epinephelidae) indicates Pleistocene isolation across the Indo-Pacific Barrier with contemporary overlap in the Coral Triangle. BMC Evol. Biol. 11, 189 (2011).PubMed 
PubMed Central 
Article 

Google Scholar 
31.Timm, J. & Kochzius, M. Geological history and oceanography of the Indo-Malay Archipelago shape the genetic population structure in the false clown anemonefish (Amphiprion ocellaris). Mol. Ecol. 17, 3999–4014 (2008).PubMed 
Article 

Google Scholar 
32.Otwoma, L. M. & Kochzius, M. Genetic population structure of the coral reef sea star Linckia laevigata in the Western Indian Ocean and Indo-West Pacific. PLoS ONE 11, 10 (2016).Article 
CAS 

Google Scholar 
33.Williams, S. T., Jara, J., Gomez, E. & Knowlton, N. The marine Indo-West Pacific break: Contrasting the resolving power of mitochondrial and nuclear genes. Integr. Comp. Biol. 42, 941–952 (2002).CAS 
PubMed 
Article 

Google Scholar 
34.Supmee, V., Sangthong, P., Songrak, A. & Suppapan, J. Population genetic structure of Asiatic Hard Clam (Meretrix meretrix) in Thailand based on Cytochrome Oxidase subunit I gene sequence. Biodiversitas 21, 2702–2709 (2020).Article 

Google Scholar 
35.Hui, M. et al. Comparative genetic population structure of three endangered giant clams (Cardiidae: Tridacna species) throughout the Indo-West Pacific: Implications for divergence, connectivity and conservation. J. Molluscan Stud. 82, 403–414 (2016).Article 

Google Scholar 
36.Panithanarak, T., Karuwancharoen, R., Na-Nakorn, U. & Nguyen, T. T. Population genetics of the spotted seahorse (Hippocampus kuda) in Thai waters: Implications for conservation. Zool. Stud. 49, 564–576 (2010).CAS 

Google Scholar 
37.Kasim, N. S. et al. Recent population expansion of longtail tuna Thunnus tonggol (Bleeker, 1851) inferred from the mitochondrial DNA markers. PeerJ 8, 9679 (2020).Article 

Google Scholar 
38.Canales-Aguirre, C. B., Ferrada-Fuentes, S., Galleguillos, R., Oyarzun, F. X. & Hernández, C. E. Population genetic structure of Patagonian toothfish (Dissostichus eleginoides) in the Southeast Pacific and Southwest Atlantic Ocean. PeerJ 6, e4173 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
39.Sato, M. et al. Genetic structure and demographic connectivity of marbled flounder (Pseudopleuronectes yokohamae) populations of Tokyo Bay. J. Sea Res. 142, 79–90 (2018).ADS 
Article 

Google Scholar 
40.Borsa, P. Genetic structure of round scad mackerel Decapterus macrosoma (Carangidae) in the Indo-Malay archipelago. Mar. Biol. 142, 575–581 (2003).CAS 
Article 

Google Scholar 
41.Eytan, R. I. & Hellberg, M. E. Nuclear and mitochondrial sequence data reveal and conceal different demographic histories and population genetic processes in Caribbean reef fishes. Evolution 64, 3380–3397 (2010).CAS 
PubMed 
Article 

Google Scholar 
42.Tan, M. P., Jamsari, A. F. J. & Siti Azizah, M. N. Genotyping of microsatellite markers to study genetic structure of the wild striped snakehead Channa striata in Malaysia. J. Fish. Biol. 88, 1932–1948 (2016).CAS 
PubMed 
Article 

Google Scholar 
43.Tan, M. P., Jamsari, A. F. J. & Siti Azizah, M. N. Phylogeographic pattern of the striped snakehead, Channa striata in Sundaland: Ancient river connectivity, geographical and anthropogenic signatures. PLoS ONE 7, 1–11 (2012).
Google Scholar 
44.Tan, M. P., Jamsari, A. F. J., Muhlisin, Z. A. & Siti Azizah, M. N. Mitochondrial genetic variation and population structure of the striped snakehead, Channa striata in Malaysia and Sumatra. Indonesia. Biochem. Syst. Ecol. 60, 99–105 (2015).CAS 
Article 

Google Scholar 
45.Haponski, A. E. & Stepien, C. A. Phylogenetic and biogeographical relationships of the Sander pikeperches (Percidae: Perciformes): Patterns across North America and Eurasia. Biol. J. Linn. Soc. Lond. 110, 156–179 (2013).Article 

Google Scholar 
46.Milá, B., Van Tassell, J. L., Calderón, J. A., Rüber, L. & Zardoya, R. Cryptic lineage divergence in marine environments: Genetic differentiation at multiple spatial and temporal scales in the widespread intertidal goby Gobiosoma bosc. Ecol. Evol. 7, 5514–5523 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
47.Piganeau, G., Gardner, M. & Eyre-Walker, A. A broad survey of recombination in animal mitochondria. Mol. Biol. Evol. 21, 2319–2325 (2004).CAS 
PubMed 
Article 

Google Scholar 
48.Avise, J. C. Molecular Markers, Natural History, and Evolution (Sinauer Associates Inc, 2004).
Google Scholar 
49.De Mandal, S., Chhakchhuak, L., Gurusubramanian, G. & Kumar, N. S. Mitochondrial markers for identification and phylogenetic studies in insects–A Review. DNA Barcodes 2, 1–9 (2014).CAS 
Article 

Google Scholar 
50.Simon, C. et al. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann. Entomol. Soc. Am. 87, 651–701 (1994).CAS 
Article 

Google Scholar 
51.Hoofer, S. R., Reeder, S. A., Hansen, E. W. & Van Den Bussche, R. A. Molecular phylogenetics and taxonomic review of noctilionoid and vespertilionoid bats (Chiroptera: Yangochiroptera). J. Mammal. 84, 809–821 (2003).Article 

Google Scholar 
52.Hewitt, G. M. Speciation, hybrid zones and phylogeography or seeing genes in space and time. Mol. Ecol. 10, 537–549 (2001).CAS 
PubMed 
Article 

Google Scholar 
53.Surya, S. et al. Morphometry and length-weight relationship of Uranoscopus marmoratus Cuvier, 1829 (Family: Uranoscopidae) from Palk Bay, India. Int. J. Biol. Sci. 5, 1–10 (2016).ADS 

Google Scholar 
54.Narejo, N. T. Morphometric characters and their relationships in Gudusia chapra (Hamilton) from Keenjhar Lake (Distt: Thatta), Sindh, Pakistan. Pak. J. Zool. 42, 101–104 (2010).
Google Scholar 
55.Gan, H. M., Nur Ilham Syahadah, M. Y., Vilasri, V., Tun Nurul Aimi, M. J. & Tan, M. P. Four whole mitogenome sequences of yellowtail stargazers (Uranoscopus cognatus cantor 1849) from East Peninsular Malaysia and West Coast of Thailand. Mitochondrial DNA B 4, 256–258 (2019).Article 

Google Scholar 
56.Panjarat, S. Sustainable fisheries in the Andaman Sea coast of Thailand. Division for Ocean Affairs and the Law of the Sea Office of Legal Affairs. (The United Nations, 2008).57.Derrick, B., Noranarttragoon, P., Zeller, D., Teh, L. C. & Pauly, D. Thailand’s missing marine fisheries catch (1950–2014). Front. Mar. Sci. 4, 402 (2017).Article 

Google Scholar 
58.Sampantamit, T., Ho, L., Van Echelpoel, W., Lachat, C. & Goethals, P. Links and trade-offs between fisheries and environmental protection in relation to the sustainable development goals in Thailand. Water 12, 399 (2020).Article 

Google Scholar 
59.Chong, V., Lee, P. & Lau, C. Diversity, extinction risk and conservation of Malaysian fishes. J. Fish Biol. 76, 2009–2066 (2010).CAS 
PubMed 
Article 

Google Scholar 
60.Lim, H. C., Ahmad, A. T., Nuruddin, A. A. & Mohd Nor, S. A. Cytochrome b gene reveals panmixia among Japanese Threadfin Bream, Nemipterus japonicus (Bloch, 1791) populations along the coasts of Peninsular Malaysia and provides evidence of a cryptic species. Mitochondrial DNA A 27, 575–584 (2016).CAS 
Article 

Google Scholar 
61.Nabilsyafiq, M. H. et al. ND5 gene marker reveals recent population expansion of wild pearse’s mudskipper (Periophthalmus novemradiatus Hamilton) inhabits Setiu wetlands in east Peninsular Malaysia. Malays. Appl. Biol. 48, 87–93 (2019).
Google Scholar 
62.Zhou, Y. et al. Importance of incomplete lineage sorting and introgression in the origin of shared genetic variation between two closely related pines with overlapping distributions. Heredity 118, 211–220 (2017).CAS 
PubMed 
Article 

Google Scholar 
63.Lessios, H. A. The great American schism: divergence of marine organisms after the rise of the Central American Isthmus. Annu. Rev. Ecol. Evol. Syst. 39, 63–91 (2008).Article 

Google Scholar 
64.Avise, J. Molecular Markers, Natural History and Evolution (Chapman and Hall, 1994).Book 

Google Scholar 
65.Nelson, J. S., Grande, T. C. & Wilson, M. V. Fishes of the World (John Wiley and Sons, 2016).Book 

Google Scholar 
66.Young, J. Z. Memoirs: On the autonomic nervous system of the Teleostean Fish Uranoscopus scaber. J. Cell Sci. 2, 491–536 (1931).Article 

Google Scholar 
67.Day, J., Clark, J. A., Williamson, J. E., Brown, C. & Gillings, M. Population genetic analyses reveal female reproductive philopatry in the oviparous Port Jackson shark. Mar. Freshw. Res. 70, 986–994 (2019).Article 

Google Scholar 
68.Roycroft, E. J., Le Port, A. & Lavery, S. D. Population structure and male-biased dispersal in the short-tail stingray Bathytoshia brevicaudata (Myliobatoidei: Dasyatidae). Conserv. Genet. 20, 717–728 (2019).CAS 
Article 

Google Scholar 
69.King, T. L., Eackles, M. S., Spidle, A. P. & Brockmann, H. J. Regional differentiation and sex-biased dispersal among populations of the horseshoe crab Limulus polyphemus. Trans. Am. Fish. Soc. 134, 441–465 (2005).Article 

Google Scholar 
70.Lane, A. & Shine, R. Intraspecific variation in the direction and degree of sex-biased dispersal among sea-snake populations. Mol. Ecol. 20, 1870–1876 (2011).PubMed 
Article 

Google Scholar 
71.Casale, P., Laurent, L., Gerosa, G. & Argano, R. Molecular evidence of male-biased dispersal in loggerhead turtle juveniles. J. Exp. Mar. Biol. Ecol. 267, 139–145 (2002).CAS 
Article 

Google Scholar 
72.Wyrtki, K. Physical Oceanography of the Southeast Asian Waters (University of California, 1961).
Google Scholar 
73.Barber, P., Palumbi, S., Erdmann, M. & Moosa, M. Sharp genetic breaks among populations of Haptosquilla pulchella (Stomatopoda) indicate limits to larval transport: patterns, causes, and consequences. Mol. Ecol. 11, 659–674 (2002).CAS 
PubMed 
Article 

Google Scholar 
74.Kamarudin, K. R. & Esa, Y. Phylogeny and phylogeography of Barbonymus schwanenfeldii (Cyprinidae) from Malaysia inferred using partial cytochrome b mtDNA gene. J. Trop. Biol. Conserv. 5, 1–13 (2009).
Google Scholar 
75.Tan, M. P. et al. Genetic diversity of the Pearse’s Mudskipper Periophthalmus novemradiatus (Perciformes: Gobiidae) and characterization of its complete mitochondrial genome. Thalassas 36, 103–113 (2020).Article 

Google Scholar 
76.Roberts, T. R. & Khaironizam, M. Z. Trophic polymorphism in the Malaysian fish Neolissochilus soroides and other old world barbs (Teleostei, Cyprinidae). Nat. Hist. Bull. Siam Soc. 56, 25–53 (2008).
Google Scholar 
77.Forsman, A. Rethinking phenotypic plasticity and its consequences for individuals, populations and species. Heredity 115, 276 (2015).CAS 
PubMed 
Article 

Google Scholar 
78.Vasileva, E. Morphokaryological variability and divergence of stargazers (Uranoscopus, perciformes) from the Mediterranean Sea basin: I. Divergence and taxonomic state of the Black Sea Stargazer. J. Ichthyol. 52, 476–484 (2012).Article 

Google Scholar 
79.Aljanabi, S. M. & Martinez, I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Ress 25, 4692–4693 (1997).CAS 
Article 

Google Scholar 
80.Ward, R. D., Zemlak, T. S., Innes, B. H., Last, P. R. & Hebert, P. D. DNA barcoding Australia’s fish species. Philos. Trans. R. Soc. Lond., B Biol. Sci. 360, 1847–1857 (2005).CAS 
Article 

Google Scholar 
81.López, J. A., Chen, W. J. & Ortí, G. Esociform phylogeny. Copeia 2004, 449–464 (2004).Article 

Google Scholar 
82.Mat Jaafar, T. N., Taylor, M. I., Mohd Nor, S. A., Bruyn, M. D. & Carvalho, G. R. Comparative genetic stock structure in three species of commercially exploited Indo-Malay Carangidae (Teleosteii, Perciformes). J. Fish Biol. 96, 337–349 (2020).CAS 
PubMed 
Article 

Google Scholar 
83.Tamura, K., Stecher, G., Peterson, D., Filipski, A. & Kumar, S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
84.Farris, J. S., Källersjö, M., Kluge, A. G. & Bult, C. Testing significance of incongruence. Cladistics 10, 315–319 (1994).Article 

Google Scholar 
85.Swofford, D. L. PAUP: Phylogenetic Analysis Using Parsimony (and Other Methods), Version 4.0 Beta 10 (Sinauer Associates, 2002).
Google Scholar 
86.Librado, P. & Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452 (2009).CAS 
PubMed 
Article 

Google Scholar 
87.Excoffier, L. & Lischer, H. E. L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567 (2010).PubMed 
Article 

Google Scholar 
88.Hasegawa, M., Kishino, H. & Yano, T. A. Dating of the huma-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol. 22, 160–174 (1985).CAS 
PubMed 
Article 

Google Scholar 
89.Kimura, M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucletide sequences. J. Mol. Evol. 16, 111–120 (1980).ADS 
CAS 
Article 

Google Scholar 
90.Felsenstein, J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 39, 783–791 (1985).PubMed 
Article 

Google Scholar 
91.Bandelt, H. J., Forster, P. & Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (1999).CAS 
PubMed 
Article 

Google Scholar 
92.Dupanloup, I., Schneider, S. & Excoffier, L. A simulated annealing approach to define the genetic structure of populations. Mol. Ecol. 11, 2571–2581 (2002).CAS 
PubMed 
Article 

Google Scholar 
93.Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Stat. Methodol. 57, 289–300 (1995).MathSciNet 
MATH 

Google Scholar 
94.Nei, M. Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. U.S.A. 70, 3321–3323 (1973).ADS 
CAS 
PubMed 
PubMed Central 
MATH 
Article 

Google Scholar 
95.Lynch, M. & Crease, T. The analysis of population survey data on DNA sequence variation. Mol. Biol. Evol. 7, 377–394 (1990).CAS 
PubMed 

Google Scholar 
96.Hudson, R. R., Slatkin, M. & Maddison, W. P. Estimation of levels of gene flow from DNA sequence data. Genetics 132, 583–589 (1992).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
97.Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123, 585–595 (1989).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
98.Fu, Y.-X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147, 915–925 (1997).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
99.Harpending, H. Infertility and forager demography. Am. J. Phys. Anthropol. 93, 385–390 (1994).CAS 
PubMed 
Article 

Google Scholar 
100.Slatkin, M. & Hudson, R. R. Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics 129, 555–562 (1991).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
101.Rogers, A. R. & Harpending, H. Population growth makes waves in the distribution of pairwise genetic differences. Mol. Biol. Evol. 9, 552–569 (1992).CAS 
PubMed 

Google Scholar 
102.Schneider, S. & Excoffier, L. Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: Application to human mitochondrial DNA. Genetics 152, 1079–1089 (1999).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
103.Yildirim, Y. Genetic structure of Pleurobranchaea maculata in New Zealand (Massey University, 2016).
Google Scholar 
104.Hubbs, C. & Lagler, K. The fishes of the Great Lakes region 213 (The University of Michigan Press, 1958).
Google Scholar 
105.Kishimoto, H. A new stargazer, Uranoscopus flavipinnis, from Japan and Taiwan with redescription and neotype designation of U. japonicus. Japan. J. Ichthyol. 34, 1–14 (1987).
Google Scholar 
106.Kishimoto, H. Redescription and lectotype designation of the stargazer Uranoscopus kaianus Günther. Copeia 1984, 1009–1011 (1984).Article 

Google Scholar 
107.Gomon, M. F. & Johnson, J. A new fringed stargazer (Uranoscopidae: Ichthyscopus) with descriptions of the other Australian species of the genus. Mem. Queensl. Mus. 43, 597–619 (1999).
Google Scholar 
108.Rainboth, W. J. Fishes of the Cambodian Mekong (Food and Agriculture Org, 1996).
Google Scholar 
109.Imamura, H. & Matsuura, K. Redefinition and phylogenetic relationships of the family Pinguipedidae (Teleostei: Perciformes). Ichthyol. Res. 50, 259–269 (2003).Article 

Google Scholar 
110.Hammer, Ø., Harper, D. A. & Ryan, P. D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4, 9 (2001).
Google Scholar 
111.Seah, Y. G., Nabilsyafiq, M. & Mazlan, A. G. Preliminary study on the morphology and biology of coexist Nemipterus furcosus and Nemipterus tambuloides from Terengganu Waters Peninsular Malaysia. J. Fish. Aquat. Sci. 11, 418–424 (2016).Article 

Google Scholar 
112.Johnson, R. A. & Wichern, D. W. Multivariate statistical analysis (Prentice Hall Upper Saddle River, 1998).MATH 

Google Scholar  More

AnimalsTissue samples for expression analysis were collected on August 24 and September 9 2016 from a population of spring-spawning herring kept in captivity at University of Bergen; the rearing of captive herring was approved by the Norwegian national animal ethics committee (Forsøksdyrutvalget FOTS ID-5072). The tissue samples used for ATAC-sequencing were collected at Hästskär on June 19 2019 from wild spring-spawning Atlantic herring in the Baltic Sea, which do not require ethical permission.Genome scan and genetic diversity at the TSHR locusWe used a 2 × 2 contingency X2 test to estimate the extent of SNP allele frequency differences between seven spring- and seven autumn-spawning herring populations from the Northeast Atlantic (Supplementary Table 2), and thus, identify the major genomic regions associated with seasonal reproduction. The SNP allele frequencies were generated in a previous study16 and were derived from Pool-seq data. For the X2 test, we formed two groups, spring and autumn spawners, and summed the reads supporting the reference and the alternative alleles for the pools included in each group.To characterize genetic diversity at the TSHR locus, we calculated nucleotide diversity (π) and Tajima’s D for the same seven spring- and seven autumn-spawning Atlantic herring populations used in the genome scan (n = ~41–100 per pool) (Supplementary Table 2). The whole-genome re-sequence data of these pools were previously reported by Han et al.16 (for details of the pools used here see Supplementary Table 2). Unbiased nucleotide diversity π and Tajima’s D were calculated for each pool using the program PoPoolation 1.2.236, which accounts for the truncated allele frequency spectrum of pooled data. In brief, a pileup file of chromosome 15 was generated from each pool BAM file using samtools v.1.1037,38. Indels and 5 bp around indels were removed to exclude spurious SNPs due to misalignments around indels. To minimize biases in the π and Tajima’s D calculations, which are sensitive to sequencing errors and coverage fluctuations39, the coverage of each pileup file was subsampled without replacement to a uniform value based on a per-pool coverage distribution (the target coverage corresponded to the 5th percentile of the coverage frequency distribution, which was used as the minimum coverage allowed for an SNP to be included in the analysis) (Supplementary Table 2). We also calculated the diversity parameters but skipping the coverage subsampling step and obtained very similar results with both approaches (Supplementary Fig. 5), thus, we decided to keep working with the subsampled datasets as coverage subsampling is recommended by the software developers36. To exclude spurious SNPs associated with repetitive sequences and copy number variants, we applied a maximum coverage filter equivalent to the per-pool 99th percentile of the coverage frequency distribution (Supplementary Table 2). A minimum base quality of 20, a minimum mapping quality of 20, and a minor allele count of 2 were required to retain high quality SNPs for further analysis. Both, π and Tajima’s D statistics were calculated using a sliding window approach with a window size of 10 kb and a step size of 2 kb (the selected window-step combination offered a good genomic resolution while reducing the noise from single SNPs, after testing windows of 5, 10, 20, 40, 50, and 100 kb for non-overlapping and overlapping windows with a step size equivalent to 20% of the window size). Only windows with a coverage fraction ≥ 0.5 were included in the computations. In addition, we estimated the effective allele frequency difference, or delta allele frequency (dAF), between spring and autumn spawning groups at the TSHR locus using the formula dAF = abs(mean(spring pools)−mean(autumn pools)). For each of the diversity parameters, we assessed whether the mean differences between sets of SNPs within chr 15: 8.85–8.95 Mb (215 SNPs) and outside (214 635 SNPs) the TSHR region were statistically significant among spring- and autumn-spawning groups using a Wilcoxon test. Data postprocessing, statistical tests, and plotting were performed in the R environment40 (for specific parameters used in PoPoolation, see the associated code to this publication).Identification of the 5.2 kb structural variantSequences spanning the entire TSHR gene plus 10 kb upstream and downstream from two reference assemblies, ASM96633v115 and Ch_v2.0.218, were aligned using BLAST41 and the output was subsequently processed with a custom R script40. Repeats were annotated by CENSOR21 for the region harboring the 5.2 kb structural variant. To validate the structural variant, long-range PCR was performed with genomic DNA from two spring- and autumn-spawning Atlantic herring in a 20 μL reaction containing 0.8 mM dNTPs, 0.3 μM each of the forward and reverse 5.2kb-confirm primers (Supplementary Table 1) and 0.75 U PrimeSTAR GXL DNA Polymerase (TaKaRa) following the program: 95 °C for 3 min, 35 cycles of 98 °C for 10 s, 58 °C for 20 s and 68 °C for 2 min 30 s, and a final extension of 10 min at 68 °C.ATAC-seq analysisBSH and brain without BSH were dissected from two spring-spawning herring caught in the Baltic Sea and transported to the lab on dry ice, then kept at –80 °C before nuclei isolation. ATAC-seq libraries were constructed according to the Omni-ATAC protocol with minor modifications42. Briefly, tissue was homogenized in 2 ml homogenization buffer (5 mM CaCl2, 3 mM Mg(Ac)2, 10 mM Tris-HCl (pH = 7.8), 0.017 mM PMSF, 0.17 mM ß-mercaptoethanol, 320 mM Sucrose, 0.1 mM EDTA and 0.1% NP-40) with a Dounce homogenizer on ice. 400 μl suspension was transferred to a 2 ml tube for the density gradient centrifugation with different concentrations of Iodixanol solution. After centrifugation, a 200 μl fraction containing the nuclei band was collected, stained with Trypan blue and counted with a Countess II Automated Cell Counter (Thermo Fisher Scientific). An aliquot of 100,000–200,000 nuclei was used as input in a 50 μl transposition reaction containing 2 X TD buffer and 100 nM assembled Tn5 transposase for a 30-min incubation at 37 °C. Tagmented DNA was purified with a Zymo clean kit (Zymo Research). Purified DNA was used for an initial pre-amplification for 5 cycles, and the additional amplification cycle was determined by qPCR based on the “R vs Cycle Number” plot43. Amplified libraries were purified with a Zymo clean kit again, and library concentrations and qualities were evaluated using the 2200 TapeStation System (Agilent Technologies).ATAC-seq was performed with a MiniSeq High Output Kit (150 cycles) on a MiniSeq instrument (Illumina) and 7–9 million reads were generated for each ATAC-seq library. Quality control, trimming, mapping, and peak calling of the sequenced reads were conducted following the ENCODE ATAC-seq pipeline (https://www.encodeproject.org/atac-seq/). The trimmed reads were aligned to the Atlantic herring reference genome (Ch_v2.0.2)18 with Bowtie244 and the mapping rate was 85–95%. Duplicate reads, reads with low mapping quality and those aligned to the mitochondria genome were removed. The remaining reads (4–5 million) were subjected to peak calling by MACS245, where 22–32 K peaks were called. Sequenced library qualities were further evaluated by calculating the TSS enrichment score and checking the library complexity with the Non-Redundant Fraction (NRF) and PCR Bottlenecking Coefficients (PBC1 and 2). Finally, conserved peaks between two biological replicates were identified by evaluating the irreproducible discovery rate (IDR).Genotyping of six differentiated variants and haplotype analysisAll six genetic variants, including the 5.2 kb structural variant, two non-coding SNPs, two missense SNPs and the copy number variant of C-terminal 22aa repeat, were genotyped in 45 spring-, 67 autumn-spawning Atlantic herring and 13 Pacific herring. TaqMan Custom SNP assays were performed to genotype the four SNPs in 5 μl reactions with a template of 20 ng genomic DNA (ThermoFisher Scientific). Copy number of the C-terminal 22aa repeat was determined by the PCR product size generated with geno22aa primers. Genotyping of the 5.2 kb structural variant was performed in a PCR reaction containing two forward primers (geno5.2kb-1F and geno5.2kb-2F) and one reverse primer (geno5.2kb-R), which generated PCR products with different sizes between spring and autumn spawners. All the primers used for genotyping are listed in Supplementary Table 1.Tissue expression profiles by quantitative PCRTotal RNA was prepared from gonad, heart, spleen, kidney, gills, intestine, hypothalamus and saccus vasculosus (BSH), and brain without BSH (brain) of six adult spring-spawning Atlantic herring using RNeasy Mini Kit (Qiagen). RNA was then reverse transcribed into cDNA with a High-Capacity cDNA Reverse Transcription Kit (ThermoFisher Scientific). TaqMan Gene Expression assay (ThermoFisher Scientific) containing 0.3 μM primers and 0.25 μM TaqMan probe (Integrated DNA Technologies) was performed to compare the relative expression levels of TSHR among different tissues. qPCR with SYBR Green chemistry was used for TSHB and DIO2 in a 10 μl reaction of SYBR Green PCR Master Mix (ThermoFisher Scientific) and 0.3 μM primers, with a program composed of an initial denaturation for 10 min at 95 °C followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. Ct values were first normalized to the housekeeping gene ACTIN, then the average expression for each gene in the gonad was assumed to be 1 for the subsequent calculation of the relative expression in other tissues.Plasmid constructsThe coding sequence for the herring single-chain TSH (scTSH) was designed following a strategy previously used for mammalian gonadotropins46 that contained an in-frame fusion of the cDNA sequences (5′–3′) of herring TSH beta subunit (NCBI: XM_012836756.1) and alpha subunit (NCBI: XM_012822755.1) linked by six histidines and then the C-terminal peptide of the hCG beta subunit. The designed sequence should generate a protein with a size of 30.6 kDa. Both scTSH and spring herring TSHR cDNA sequences were synthesized in vitro and cloned in the expression vector pcDNA3.1 by Genscript (Leiden, Netherlands). pcDNA3.1 plasmid expressing human TSHR was kindly provided by Drs. Gilbert Vassart and Sabine Costagliola (Université libre de Bruxelles, Belgium). Then, the spring herring TSHR and human TSHR plasmids were used as templates for site-directed mutagenesis to generate constructs coding for different mutant herring or human TSHRs. Plasmids for the dual-luciferase assay, including pGL4.29[luc2P/CRE/Hygro] containing cAMP response elements (CREs) to drive the transcription of luciferase gene luc2P and pRL-TK monitoring the transfection efficiency, were purchased from Promega. Five ng of each plasmid was used to transform the XL1-Blue competent cells (Agilent), plasmid DNA was subsequently extracted from 200 ml overnight culture of a single transformant clone using an EndoFree plasmid Maxi Kit (Qiagen).Cell cultureChinese hamster ovary (CHO) (ATCC CCL-61) and human embryonic kidney 293 (HEK293) (ATCC CRL-1573) cells were maintained in DMEM supplemented with 5% (CHO) or 10% FBS (HEK293), 100 U/ml penicillin, 100 μg/ml streptomycin and 292 μg/ml l-Glutamine (ThermoFisher Scientific) at 37 °C with 5% CO2. Epithelioma Papulosum Cyprini (EPC) cells (ATCC CRL-2872) were cultured in EMEM (Sigma) supplemented with 10% FBS, 100 U/ml penicillin, 100 μg/ml streptomycin, 292 μg/ml l-Glutamine and 1 mM Sodium Pyruvate (ThermoFisher Scientific) at 26 °C with 5% CO2.Production of recombinant herring scTSHCHO cells were transfected with the scTSH expression plasmid using Lipofectamine 3000 (Invitrogen), stable clones were subsequently selected with 500 μg/ml G418 (Invitrogen) and screened for producing scTSH by western blot using a polyclonal antisera against the sea bass alpha subunit47. A positive clone was expanded in 225 cm2 cell culture flasks (Corning) in culture medium containing 5% FBS until confluence, then the cells were maintained in serum-free DMEM for hormone production for 7 days at 25 °C48. After 7 days, culture medium containing scTSH or without (negative control) was centrifuged at 15000 x g for 15 min and concentrated by ultrafiltration using Centricon Plus-70 / Ultracel PL-30 (Merck Millipore Ltd.). Then, western blotting was performed to confirm TSH production. Concentrated medium containing herring scTSH was denatured at 94 °C for 5 min in 0.1% SDS and 50 mM 2-mercaptoethanol, and then treated with 2.5 units of peptide-N-glycosidase F (Roche Diagnostics) at 37 °C for 2 h in 20 mM sodium phosphate with 0.5% Nonidet P-40, pH 7.5. All samples were run in 12% SDS-PAGE in the reducing condition and transferred to a PVDF membrane (Immobilon P; Millipore Corp.), then blocked overnight with 5% skimmed milk at 4 °C. After blocking, the membrane was incubated with polyclonal antisera against the sea bass alpha subunit (dilution 1:2000) for 90 min at room temperature, washed, and then further incubated with 1:25000 goat anti-rabbit immunoglobulin G (IgG) horseradish peroxidase conjugate (Bio-Rad Laboratories) for 60 min at room temperature. Immunodetection was performed by chemiluminescence with a Pierce ECL Plus Western Blotting Substrate kit (ThermoFisher Scientific).Cell surface expressionA Rhotag (MNGTEGPNFYVPFSNKTGVVYEE) was inserted at the N-terminus of herring TSHR for flow cytometry analysis of receptor cell surface expression. Anti-Rhotag polyclonal antibody was kindly provided by Drs. Gilbert Vassart and Sabine Costagliola (Université Libre de Bruxelles, Brussels, Belgium). PBS containing 1% BSA and 0.05% sodium azide was prepared as the flow cytometry (FCM) buffer for the washing and antibody incubation steps. 2.2 × 106 EPC cells were seeded in a 100 mm poly-d-Lysine-treated petri dish the day before transfection. Each dish was transfected with 10 μg TSHR or empty pcDNA3.1 expression plasmid using 20 μl jetPRIME transfection reagent in 500 μl jetPRIME transfection buffer (Polyplus transfection). Cells were harvested 24 h after transfection, then washed once in cold PBS and fixed in 2% PFA for 10 min at room temperature. After fixation, cells were washed three times with FCM buffer, then incubated with anti-Rhotag antibody or FCM buffer (negative control) for 1 h at room temperature. Cells were washed again with FCM buffer three times and stained with Alexa Fluor 488-labeled chicken anti-mouse IgG (H + L) antibody (1:200 dilution, ThermoFisher Scientific) or FCM buffer (negative control) for 45 min in the dark. After the fluorescent staining, cells were washed three times and resuspended in FCM buffer before analysis on a CytoFLEX instrument (Beckman Coulter). A minimum of 100,000 events was recorded for each sample, fluorescence intensities of negative control and cells transfected with empty pcDNA3.1 plasmid were used as the background for gating strategy. Cell surface expression was represented by the mean fluorescence intensity of the positively stained cell population.Dual-luciferase reporter assayEPC or HEK293 cells were plated in a 48-well plate at a density of 1 × 105 cells/well the day before transfection. A total of 250 ng plasmid mixture containing pGL4.29[luc2P/CRE/Hygro], TSHR expression plasmid (or empty pcDNA3.1) and pRL-TK with the ratio of 20:5:1 was prepared to transfect each well of cells using jetPRIME transfection reagent (Polyplus). Medium was replaced by fresh medium containing 10% FBS (TSH-induced condition) or serum-free medium (constitutive activity condition) 4 h after transfection. On day three, cells were treated with serum-free medium containing different dilutions of the concentrated scTSH medium for 4 h (TSH-induced condition) or directly subjected to the luminescence measurement without TSH induction (constitutive activity condition). Luminescence was measured using a Dual-Luciferase Reporter assay (Promega) on an Infinite M200 Microplate Reader (Tecan Group Ltd., Switzerland), and luciferase activity was represented as the ratio of firefly (pGL4.29[luc2P/CRE/Hygro]) to Renilla (pRL-TK) luminescence.5′-RACE to identify the herring DIO2 TSSTotal RNA was prepared from brain of a spring-spawning Atlantic herring using the RNeasy Mini Kit (Qiagen). Six μg of the isolated RNA was used for 5′-RACE with a FirstChoiceTM RLM-RACE Kit (ThermoFisher Scientific). One μl cDNA or Outer RACE PCR product was used as PCR template in a 20 μL reaction containing 0.8 mM dNTPs, 0.3 μM of each forward and reverse primer (Supplementary Table 1) and 0.75 U PrimeSTAR GXL DNA Polymerase (TaKaRa). Amplification was carried out with an initial denaturation of 3 min at 95 °C, followed by 35 cycles of 98 °C for 10 s, 58 °C for 20 s and 68 °C for 40 s, and a final extension of 10 min at 68 °C. The final 5′ RACE product was sequenced at Eurofins Genomics (Ebersberg, Germany).Sequence conservation analysisGenomic sequences covering the TSHR locus were extracted from Ensembl Genome Browser for Atlantic herring and 11 other fish species, including Amazon molly (Poecilia formosa), denticle herring (Denticeps clupeoides), goldfish (Carassius auratus), guppy (Poecilia reticulata), Neolamprologus brichardi, Japanese medaka (Oryzias latipes), northern pike (Esox lucius), orange clownfish (Amphiprion percula), spotted gar (Lepisosteus oculatus), three-spined stickleback (Gasterosteus aculeatus) and spotted green pufferfish (Tetraodon nigroviridis). The extracted sequences were firstly aligned using progressiveCactus49,50, and a subsequent alignment was generated using the hal2maf program from halTools51 with Atlantic herring assembly (Ch_v2.0.2)18 as the co-ordinate backbone. This alignment was used for the downstream phastCons score calculation by running phyloFit24 and phastCons25 from the PHAST package with default parameters. Peaks were called by grouping signals with a minimum phastCons score of 0.2 within 500 bp region.Structure modeling of human and herring TSHRsIn order to explore the possible interactions of the variant residues with other receptor interacting proteins and to study intramolecular interactions, we built a structural homology model for the herring TSHR (herrTSHR) complexed with herring TSH and Gs-protein. The TSHR hinge region that harbors the Q370H substitution and the C-terminus containing the 22aa repeat were excluded from the homology model due to the lack of structural templates for these regions. The homology model was constructed by using the following structural templates of evolutionarily related class A GPCRs: (i) the leucine-rich repeat domain (LRRD) complexed with hormone was modeled based on the solved FSHR LRRD – FSH complex structure (Protein Data Bank (PDB) ID: 4AY9)52,53, this part of model included herring TSHR Cys33 – Asn296 and fragments of the hinge region Gln297 – Thr312 and Ser393 – Ile421; (ii) the available structural complex of β2-adrenoreceptor with Gs-protein (PDB ID: 3SN6)54 was used as the template to model the seven-transmembrane helix domain (7TMD) of herring TSHR in the active conformation; (iii) the extracellular loop 2 (ECL2) was built by using the ECL2 of μ-opioid receptor (PDB ID: 6DDE)55. To prepare the template for herring TSHR modeling, the fused T4-lysozyme and bound ligand of β2-adrenoreceptor were deleted, the ECL1 and ECL3 loops were adjusted manually to the loop length of herring TSHR. Due to the lack of third intracellular loop (ICL3) in the β2-adrenoreceptor structure, amino acid residues of herring TSHR ICL3 were manually added to the template. Since herring TSHR does not have the TMH5 proline, which is highly conserved among all class A GPCRs and responsible for the helical kinks and bulges within this region56, we assumed a rather regular (stretched) helix conformation for the herring TSHR TMH5 and therefore replaced the kinked β2-adrenoreceptor TMH5 template with a regular α-helix. Moreover, the ECL2 template was substituted with μOR ECL2 structure because of its higher sequence similarity with herring TSHR in this region. Finally, amino acid residues of this chimeric 7TMD template and FSHR N-terminus were mutated to the corresponding spring herring TSHR residues and sequence of the heterodimeric FSH ligand was substituted by the herring TSH. All homology models were generated by using SYBYL-X 2.0 (Certara, NJ, US). The 7TMD structure was then fused with FSHR N-terminus at position 421. The assembled complex was subsequently optimized by the energy minimization under constrained backbone atoms (the AMBER F99 force field was used), followed by a 2 ns molecular dynamics simulation (MD) of the side chains. The entire TSHR complex was energetically minimized without any constraints until converging at a termination gradient of 0.05 kcal/mol*Å. Next, for autumn herrTSHR modeling, the spring TSHR sequence was substituted with autumn TSHR sequence. For humTSHR, the spring TSHR sequence was substituted with human TSHR, and the herring TSH ligand was replaced by the bovine TSH sequence. Both complex models were energetically minimized until converging at a termination gradient of 0.05 kcal/mol*Å.To investigate the microenvironment around the L471M mutation at TMH2 position 2.51, local short MD’s of 4 ns on Met4712.51 (spring herrTSHR), Leu471 (autumn herrTSHR) or Phe461 (humTSHR) and its surrounding amino acids were performed. During MD simulations, backbone atoms of the entire complexes as well as all side chains, except residues at positions 1.47, 1.51, 1.54, 2.48, 2.52, and 2.55 that form the hydrophobic patch around position 2.51, were constrained.Statistics and reproducibilityResults were presented as the mean + SD (standard deviation) calculated from at least four biological replicates for each experiment, and at least two independent experiments were conducted for each assay. Unpaired two-tailed Student’s t test was performed to calculate the P-values and means were judged as statistically significant when P ≤ 0.05.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

1.Clapham, P. J. in Encyclopedia of marine mammals 489–492 (Elsevier, 2018).2.Stevick, P. T. et al. A quarter of a world away: female humpback whale moves 10 000 km between breeding areas. Biol. Lett. 7, 299–302 (2010).PubMed 
PubMed Central 
Article 

Google Scholar 
3.Nicol, S. et al. Southern Ocean iron fertilization by baleen whales and Antarctic krill. Fish. Fish. 11, 203–209 (2010).Article 

Google Scholar 
4.Smetacek, V. & Nicol, S. Polar ocean ecosystems in a changing world. Nature 437, 362–368 (2005).CAS 
PubMed 
Article 

Google Scholar 
5.Dunlop, R. A. Potential motivational information encoded within humpback whale non-song vocal sounds. J. Acoustical Soc. Am. 141, 2204–2213 (2017).Article 

Google Scholar 
6.Stimpert, A. K., Au, W. W. L., Parks, S. E., Hurst, T. & Wiley, D. N. Common humpback whale (Megaptera novaeangliae) sound types for passive acoustic monitoring. J. Acoustical Soc. Am. 129, 476–482 (2011).Article 

Google Scholar 
7.Van Opzeeland, I., Van Parijs, S., Kindermann, L., Burkhardt, E. & Boebel, O. Calling in the cold: pervasive acoustic presence of humpback whales (Megaptera novaeangliae) in Antarctic coastal waters. PLoS ONE 8, 1–7 (2013).
Google Scholar 
8.Siegel, V. Biology and Ecology of Antarctic Krill. (Springer, 2016).9.Atkinson, A. et al. Krill (Euphausia superba) distribution contracts southward during rapid regional warming. Nat. Clim. Change 9, 142–147 (2019).Article 

Google Scholar 
10.Loeb, V. J., Hofmann, E. E., Klinck, J. M., Holm-Hansen, O. & White, W. B. ENSO and variability of the Antarctic Peninsula pelagic marine ecosystem. Antarctic Sci. https://doi.org/10.1017/s0954102008001636 (2009).11.Loeb, V. J. & Santora, J. A. Climate variability and spatiotemporal dynamics of five Southern Ocean krill species. Prog. Oceanogr. 134, 93–122 (2015).Article 

Google Scholar 
12.Bombosch, A. et al. Predictive habitat modelling of humpback (Megaptera novaeangliae) and Antarctic minke (Balaenoptera bonaerensis) whales in the Southern Ocean as a planning tool for seismic surveys. Deep-Sea Res. Part I: Oceanographic Res. Pap. 91, 101–114 (2014).Article 

Google Scholar 
13.Brierley, A. S. et al. Antarctic krill under sea ice: elevated abundance in a narrow band just south of ice edge. Science 295, 1890–1892 (2002).CAS 
PubMed 
Article 

Google Scholar 
14.Rettig, S. et al. in 1st International Conference and Exhibition on Underwater Acoustics. (eds Papadakis, J. & Bjorno, L.) 1669–1674 (2013).15.Garland, E. C. et al. Humpback whale song on the Southern Ocean feeding grounds: Implications for cultural transmission. PLoS ONE 8, e79422 (2013).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
16.Stimpert, A. K., Peavey, L. E., Friedlaender, A. S. & Nowacek, D. P. Humpback whale song and foraging behavior on an Antarctic feeding ground. PLoS ONE 7, e51214 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
17.Filun, D. et al. Frozen verses: Antarctic minke whales (Balaenoptera bonaerensis) call predominantly during austral winter. R. Soc. Open Sci. 7, 192112 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
18.Boebel, O. The Expedition PS89 of the Research Vessel POLARSTERN to the Weddell Sea in 2014/2015. (Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven, 2015).19.Burkhardt, E. Whale sightings during Polarstern cruise PS96 (ANT-XXXI/2), https://doi.org/10.1594/PANGAEA.923113 (2020).20.Herr, H., Viquerat, S. & Siebert, U. Aerial cetacean survey Southern Ocean 2014/2015, https://doi.org/10.1594/PANGAEA.894938 (2018).21.Herr, H., Viquerat, S. & Siebert, U. Ship based cetacean survey Southern Ocean 2014/2015, https://doi.org/10.1594/PANGAEA.894873 (2018).22.National Oceanic and Atmospheric Administration & Department of Commerce. Climate Prediction Centre (CPC) Oceanic Nino Index (2019).23.Thomisch, K. et al. Spatio-temporal patterns in acoustic presence and distribution of Antarctic blue whales Balaenoptera musculus intermedia in the Weddell Sea. Endanger. Species Res. 30, 239–253 (2016).Article 

Google Scholar 
24.Širović, A. et al. Seasonality of blue and fin whale calls and the influence of sea ice in the Western Antarctic Peninsula. Deep Sea Res. Part II: Topical Stud. Oceanogr. 51, 2327–2344 (2004).Article 

Google Scholar 
25.Schall, E. et al. Large-scale spatial variabilities in the humpback whale acoustic presence in the Atlantic sector of the Southern Ocean. R. Soc. Open Sci. 7, 201347 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
26.Loeb, V., Hofmann, E. E., Klinck, J. M. & Holm-Hansen, O. Hydrographic control of the marine ecosystem in the South Shetland-Elephant Island and Bransfield Strait region. Deep Sea Res. Part II: Topical Stud. Oceanogr. 57, 519–542 (2010).Article 

Google Scholar 
27.Sallée, J.-B., Speer, K. & Rintoul, S. Zonally asymmetric response of the Southern Ocean mixed-layer depth to the Southern Annular Mode. Nat. Geosci. 3, 273–279 (2010).Article 
CAS 

Google Scholar 
28.Kim, Y. S. & Orsi, A. H. On the variability of Antarctic Circumpolar Current fronts inferred from 1992–2011 altimetry. J. Phys. Oceanogr. 44, 3054–3071 (2014).Article 

Google Scholar 
29.Yuan, X. ENSO-related impacts on Antarctic sea ice: a synthesis of phenomenon and mechanisms. Antarct. Sci. 16, 415 (2004).Article 

Google Scholar 
30.Meredith, M. P., Murphy, E. J., Hawker, E. J., King, J. C. & Wallace, M. I. On the interannual variability of ocean temperatures around South Georgia, Southern Ocean: Forcing by El Niño/Southern Oscillation and the southern annular mode. Deep Sea Res. Part II: Topical Stud. Oceanogr. 55, 2007–2022 (2008).Article 

Google Scholar 
31.Lovenduski, N. S. & Gruber, N. Impact of the Southern Annular Mode on Southern Ocean circulation and biology. Geophys. Res. Lett. 32, 1–4 (2005).Article 

Google Scholar 
32.Craig, A. S., Herman, L. M., Gabriele, C. M. & Pack, A. A. Migratory timing of humpback whales (Megaptera novaeangliae) in the central north Pacific varies with age, sex and reproductive status. Behaviour 140, 981–1001 (2003).Article 

Google Scholar 
33.Hofmann, E. E., Klinck, J. M., Locarnini, R. A., Fach, B. & Murphy, E. Krill transport in the Scotia Sea and environs. Antarct. Sci. 10, 406–415 (1998).Article 

Google Scholar 
34.Barendse, J. et al. Migration redefined? Seasonality, movements and group composition of humpback whales Megaptera novaeangliae off the west coast of South Africa. Afr. J. Mar. Sci. 32, 1–22 (2010).Article 

Google Scholar 
35.Witteveen, B. H., Foy, R. J., Wynne, K. M. & Tremblay, Y. Investigation of foraging habits and prey selection by humpback whales (Megaptera novaeangliae) using acoustic tags and concurrent fish surveys. Mar. Mammal. Sci. 24, 516–534 (2008).Article 

Google Scholar 
36.Brown, M. R., Corkeron, P. J., Hale, P. T., Schultz, K. W. & Bryden, M. M. Evidence for a sex-segregated migration in the humpback whale (Megaptera novaeangliae). Proc. R. Soc. Lond. B 259, 229–234 (1995).CAS 
Article 

Google Scholar 
37.International Whaling Commission. Report on the workshop on the comprehensive assessment of Southern Hemisphere humpback whales. J. Cetacea. Res. Manag. Spec. Issue 3, 1–50 (2011).
Google Scholar 
38.Findlay, K. P. et al. Humpback whale “super-groups” – A novel low-latitude feeding behaviour of Southern Hemisphere humpback whales (Megaptera novaeangliae) in the Benguela Upwelling System. PLOS ONE 12, e0172002 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
39.Gridley, T., Silva, M., Wilkinson, C., Seakamela, S. & Elwen, S. H. Song recorded near a super-group of humpback whales on a mid-latitude feeding ground off South Africa. J. Acoustical Soc. Am. 143, EL298–EL304 (2018).CAS 
Article 

Google Scholar 
40.Ross-Marsh, E., Elwen, S., Prinsloo, A., James, B. & Gridley, T. Singing in South Africa: monitoring the occurrence of humpback whale (Megaptera novaeangliae) song near the Western Cape. Bioacoustics, 1–17, https://doi.org/10.1080/09524622.2019.1710254 (2020).41.Cai, W. et al. Increasing frequency of extreme El Niño events due to greenhouse warming. Nat. Clim. Change 4, 111–116 (2014).Article 

Google Scholar 
42.Bengtson Nash, S. M. et al. Signals from the south; humpback whales carry messages of Antarctic sea‐ice ecosystem variability. Glob. Change Biol. 24, 1500–1510 (2018).Article 

Google Scholar 
43.Baumgartner, M. F. & Mussoline, S. E. A generalized baleen whale call detection and classification system. J. Acoustical Soc. Am. 129, 2889–2902 (2011).Article 

Google Scholar 
44.Klinck, H. et al. Long-range underwater vocalizations of the crabeater seal (Lobodon carcinophaga). J. Acoustical Soc. Am. 128, 474–479 (2010).Article 

Google Scholar 
45.Risch, D. et al. Mysterious bio-duck sound attributed to the Antarctic minke whale (Balaenoptera bonaerensis). Biol. Lett. 10, 20140175 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
46.Schall, E. & Van Opzeeland, I. Calls produced by Ecotype C killer whales (Orcinus orca) off the Eckstrom iceshelf, Antarctica. Aquat. Mamm. 43, 117–126 (2017).Article 

Google Scholar 
47.Van Opzeeland, I. et al. Acoustic ecology of Antarctic pinnipeds. Mar. Ecol. Prog. Ser. 414, 267–291 (2010).Article 

Google Scholar 
48.Dunlop, R. A., Cato, D. H. & Noad, M. J. Non-song acoustic communication in migrating humpback whales (Megaptera novaeangliae). Mar. Mammal. Sci. 24, 613–629 (2008).Article 

Google Scholar 
49.Schall, E. & El-Gabbas, A. Humpback-whale-acoustic-detection-and-environmental-modelling, https://github.com/elenaschall/Humpback-whale-acoustic-detection-and-environmental-modelling (GitHub, GitHub, 2021).50.Bioacoustics, Research & Program. Raven Pro: Interactive Sound Analysis Software (Version 1.5) http://ravensoundsoftware.com/ (The Cornell Lab of Ornithology, Ithaca, NY, 2014).51.Cavalieri, D., Parkinson, C., Gloersen, P. & Zwally, H. Sea ice concentrations from Nimbus-7 SMMR and DMSP SSM/I-SSMIS passive microwave data, version 1. Boulder, Colorado USA, NASA National Snow and Ice Data Center Distributed Active Archive Center 10, https://doi.org/10.5067/8GQ8LZQVL0VL (1996).52.Greene, C. A. Daily Antarctic sea ice concentration (2020).53.Marshall, G. J. Trends in the southern annular mode from observations and reanalyses. J. Clim. 16, 4134–4143 (2003).Article 

Google Scholar 
54.Marshall, G. & National Center for Atmospheric Research Staff (Eds). The climate data guide: Marshall Southern Annular Mode (SAM) index (Station-based), (2019).55.R Core Team. R: A language and environment for statistical computing, https://www.R-project.org/ (R Foundation for Statistical Computing, Vienna, Austria, 2018)56.National Oceanic and Atmospheric Administration (NOAA) & Climate Prediction Centre (CPC). Oceanic Nino Index (2019).57.Wood, S. N. Generalized additive models: an introduction with R (CRC press, 2017).58.Pinheiro, J., Bates, D., DebRoy, S. & Sarkar, D. _nlme: Linear and nonlinear mixed effects models. R package version 3.1-145 (R CoreTeam, 2020).59.Hyndman, R. et al. forecast: Forecasting functions for time series and linear models_. R. package version 8.11, http://pkg.robjhyndman.com/forecast (2020)..60.Schall, E. et al. Humpback whale acoustic presence in the Atlantic sector of the Southern Ocean, https://doi.org/10.5061/dryad.ncjsxkss0 (2021).61.Wessel, P. & Smith, W. H. A global, self‐consistent, hierarchical, high‐resolution shoreline database. J. Geophys. Res.: Solid Earth 101, 8741–8743 (1996).Article 

Google Scholar 
62.Amante, C. & Eakins, B. W. ETOPO1 arc-minute global relief model: procedures, data sources and analysis, https://doi.org/10.7289/V5C8276M (2009).63.Spreen, G., Kaleschke, L. & Heygster, G. Sea ice remote sensing using AMSR-E 89-GHz channels. J. Geophys Res. Oceans 113, C02S03 (2008).Article 

Google Scholar  More

1.De’Ath G, Fabricius KE, Sweatman H, Puotinen M. The 27-year decline of coral cover on the Great Barrier Reef and its causes. Proc Natl Acad Sci U.S.A. 2012;109:17995–9.PubMed 
PubMed Central 
Article 

Google Scholar 
2.Randall CJ, van Woesik R. Contemporary white-band disease in Caribbean corals driven by climate change. Nat Clim Chang. 2015;5:375–9.Article 

Google Scholar 
3.Maynard J, van Hooidonk R, Eakin CM, Puotinen M, Garren M, Williams G, et al. Projections of climate conditions that increase coral disease susceptibility and pathogen abundance and virulence. Nat Clim Chang. 2015;5:688–95.Article 

Google Scholar 
4.Cziesielski MJ, Schmidt-Roach S, Aranda M. The past, present, and future of coral heat stress studies. Ecol Evol. 2019;9:10055–66.PubMed 
PubMed Central 
Article 

Google Scholar 
5.Bourne D, Iida Y, Uthicke S, Smith-Keune C. Changes in coral-associated microbial communities during a bleaching event. ISME J. 2008;2:350–63.CAS 
PubMed 
Article 

Google Scholar 
6.van de Water JAJM, Chaib De Mares M, Dixon GB, Raina JB, Willis BL, Bourne DG, et al. Antimicrobial and stress responses to increased temperature and bacterial pathogen challenge in the holobiont of a reef-building coral. Mol Ecol. 2018;27:1065–80.PubMed 
Article 
CAS 

Google Scholar 
7.Sussman M, Mieog JC, Doyle J, Victor S, Willis BL, Bourne DG. Vibrio zinc-metalloprotease causes photoinactivation of coral endosymbionts and coral tissue lesions. PLoS ONE. 2009;4:1–14.8.Ben-Haim Y, Zicherman-Keren M, Rosenberg E. Temperature-regulated bleaching and lysis of the coral Pocillopora damicornis by the novel pathogen Vibrio coralliilyticus. Appl Environ Microbiol. 2003;69:4236–41.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
9.Garren M, Son K, Raina J-B, Rusconi R, Menolascina F, Shapiro OH, et al. A bacterial pathogen uses dimethylsulfoniopropionate as a cue to target heat-stressed corals. ISME J. 2014;8:999–1007.CAS 
PubMed 
Article 

Google Scholar 
10.Garren M, Son K, Tout J, Seymour JR, Stocker R. Temperature-induced behavioral switches in a bacterial coral pathogen. ISME J. 2016;10:1363–72.CAS 
PubMed 
Article 

Google Scholar 
11.Barbara GM, Mitchell JG. Marine bacterial organisation around point-like sources of amino acids. FEMS Microbiol Ecol. 2003;43:99–109.CAS 
PubMed 
Article 

Google Scholar 
12.Seymour JR, Marcos, Stocker R. Resource patch formation and exploitation throughout the marine microbial food web. Am Nat. 2009;173:E15–29.CAS 
PubMed 
Article 

Google Scholar 
13.Son K, Menolascina F, Stocker R. Speed-dependent chemotactic precision in marine bacteria. Proc Natl Acad Sci U.S.A. 2016;113:8624–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
14.Meron D, Efrony R, Johnson WR, Schaefer AL, Morris PJ, Rosenberg E, et al. Role of Flagella in virulence of the coral pathogen Vibrio coralliilyticus. Appl Environ Microbiol. 2009;75:5704–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
15.Ushijima B, Häse CC. Influence of chemotaxis and swimming patterns on the virulence of the coral pathogen Vibrio coralliilyticus. J Bacteriol. 2018;200:1–16.Article 

Google Scholar 
16.Crossland CJ, Barnes DJ, Borowitzka MA. Diurnal lipid and mucus production in the staghorn coral Acropora acuminata. Mar Biol. 1980;60:81–90.17.Davies PS. The role of zooxanthellae in the nutritional energy requirements of Pocillopora eydouxi. Coral Reefs. 1984;2:181–6.18.Rix L, de Goeij JM, Mueller CE, Struck U, Middelburg JJ, van Duyl FC, et al. Coral mucus fuels the sponge loop in warm-and cold-water coral reef ecosystems. Sci Rep. 2016;6:1–11.Article 
CAS 

Google Scholar 
19.Naumann MS, Haas A, Struck U, Mayr C, El-Zibdah M, Wild C. Organic matter release by dominant hermatypic corals of the Northern Red Sea. Coral Reefs. 2010;29:649–59.Article 

Google Scholar 
20.Wild C, Huettel M, Klueter A, Kremb SG, Rasheed MYM, Jørgensen BB. Coral mucus functions as an energy carrier and particle trap in the reef ecosystem. Nature. 2004;428:66–70.CAS 
PubMed 
Article 

Google Scholar 
21.Bythell JC, Wild C. Biology and ecology of coral mucus release. J Exp Mar Bio Ecol. 2011;408:88–93.Article 

Google Scholar 
22.Bakshani CR, Morales-Garcia AL, Althaus M, Wilcox MD, Pearson JP, Bythell JC, et al. Evolutionary conservation of the antimicrobial function of mucus: a first defence against infection. NPJ Biofilms Microbiomes. 2018;14:1–12.
Google Scholar 
23.Gibbin E, Gavish A, Krueger T, Kramarsky-Winter E, Shapiro O, Guiet R, et al. Vibrio coralliilyticus infection triggers a behavioural response and perturbs nutritional exchange and tissue integrity in a symbiotic coral. ISME J. 2019;13:989–1003.24.Gavish AR, Shapiro OH, Kramarsky-Winter E, Vardi A. Microscale tracking of coral-vibrio interactions. ISME Communications. 2021;1:1–18.25.Shapiro OH, Fernandez VI, Garren M, Guasto JS, Debaillon-Vesque FP, Kramarsky-Winter E, et al. Vortical ciliary flows actively enhance mass transport in reef corals. Proc Natl Acad Sci U.S.A. 2014;111:13391–6.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
26.Seymour JR, Ahmed T, Stocker R. A microfluidic chemotaxis assay to study microbial behavior in diffusing nutrient patches. Limnol Oceanogr Methods. 2008;6:477–88.CAS 
Article 

Google Scholar 
27.Penn K, Wang J, Fernando SC, Thompson JR. Secondary metabolite gene expression and interplay of bacterial functions in a tropical freshwater cyanobacterial bloom. ISME J. 2014;8:1866–78.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
28.Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:1–21.Article 
CAS 

Google Scholar 
29.Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:1–12.30.Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U.S.A. 2005;102:15545–50.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
31.Mootha VK, Lindgren CM, Eriksson K-F, Subramanian A, Sihag S, Lehar J, et al. PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34:267–73.CAS 
PubMed 
Article 

Google Scholar 
32.Schneider WR, Doetsch RN. Effect of viscosity on bacterial motility. J Bacteriol. 1974;117:696–701.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
33.Martinez VA, Schwarz-Linek J, Reufer M, Wilson LG, Morozov AN, Poon WCK. Flagellated bacterial motility in polymer solutions. Proc Natl Acad Sci U.S.A. 2014;111:17771–6.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
34.Kimes NE, Grim CJ, Johnson WR, Hasan NA, Tall BD, Kothary MH, et al. Temperature regulation of virulence factors in the pathogen Vibrio coralliilyticus. ISME J. 2012;6:835–46.CAS 
PubMed 
Article 

Google Scholar 
35.Kojima S, Yamamoto K, Kawagishi I, Homma M. The polar flagellar motor of Vibrio cholerae is driven by an Na+ motive force. J Bacteriol. 1999;181:1927–30.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
36.Sowa Y, Hotta H, Homma M, Ishijima A. Torque-speed relationship of the Na+-driven flagellar motor of Vibrio alginolyticus. J Mol Biol. 2003;327:1043–51.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
37.Milo R, Phillips R. Cell biology by the numbers. 1st ed. New York, NY: Garland Science; 2016.38.Crossland CJ. In situ release of mucus and DOC-lipid from the corals Acropora variabilis and Stylophora pistillata in different light regimes. Coral Reefs. 1987;6:35–42.CAS 
Article 

Google Scholar 
39.Wild C, Woyt H, Huettel M. Influence of coral mucus on nutrient fluxes in carbonate sands. Mar Ecol Prog Ser. 2005;287:87–98.40.Ducklow HW, Mitchell R. Composition of mucus released by coral reef coelenterates. Limnol Oceanogr. 1979;24:706–14.CAS 
Article 

Google Scholar 
41.Meikle P, Richards GN, Yellowlees D. Structural determination of the oligosaccharide side chains from a glycoprotein isolated from the mucus of the coral Acropora formosa. J Biol Chem. 1987;262:16941–7.CAS 
PubMed 
Article 

Google Scholar 
42.Coddeville B, Maes E, Ferrier-Pagès C, Guerardel Y. Glycan profiling of gel forming mucus layer from the scleractinian symbiotic coral Oculina arbuscula. Biomacromolecules. 2011;12:2064–73.CAS 
PubMed 
Article 

Google Scholar 
43.Hasegawa H, Häse CC. TetR-type transcriptional regulator VtpR functions as a global regulator in Vibrio tubiashii. Appl Environ Microbiol. 2009;75:7602–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
44.Ball AS, Chaparian RR, van Kessel JC. Quorum sensing gene regulation by LuxR/HapR master regulators in Vibrios. J Bacteriol. 2017;199:1–13.45.Rutherford ST, Van Kessel JC, Shao Y, Bassler BL. AphA and LuxR/HapR reciprocally control quorum sensing in vibrios. Genes Dev. 2011;25:397–408.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
46.Hammer BK, Bassler BL. Quorum sensing controls biofilm formation in Vibrio cholerae. Mol Microbiol. 2003;50:101–4.CAS 
PubMed 
Article 

Google Scholar 
47.Waters CM, Lu W, Rabinowitz JD, Bassler BL. Quorum sensing controls biofilm formation in Vibrio cholerae through modulation of cyclic Di-GMP levels and repression of vpsT. J Bacteriol. 2008;190:2527–36.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
48.Burger AH. Quorum Sensing in the Hawai’ian Coral Pathogen Vibrio coralliilyticus strain OCN008. University of Hawaii at Manoa; 2017.49.Yildiz FH, Schoolnik GK. Vibrio cholerae O1 El Tor: identification of a gene cluster required for the rugose colony type, exopolysaccharide production, chlorine resistance, and biofilm formation. Proc Natl Acad Sci U.S.A. 1999;96:4028–33.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
50.Fong JCN, Syed KA, Klose KE, Yildiz FH. Role of Vibrio polysaccharide (vps) genes in VPS production, biofilm formation and Vibrio cholerae pathogenesis. Microbiology. 2010;156:2757–69.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
51.Fong JCN, Karplus K, Schoolnik GK, Yildiz FH. Identification and characterization of RbmA, a novel protein required for the development of rugose colony morphology and biofilm structure in Vibrio cholerae. J Bacteriol. 2006;188:1049–59.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
52.Fong JCN, Yildiz FH. The rbmBCDEF gene cluster modulates development of rugose colony morphology and biofilm formation in Vibrio cholerae. J Bacteriol. 2007;189:2319–30.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
53.DiRita VJ, Mekalanos JJ. Periplasmic interaction between two membrane regulatory proteins, ToxR and ToxS, results in signal transduction and transcriptional activation. Cell. 1991;64:29–37.CAS 
PubMed 
Article 

Google Scholar 
54.Almagro-Moreno S, Root MZ, Taylor RK. Role of ToxS in the proteolytic cascade of virulence regulator ToxR in Vibrio cholerae. Mol Microbiol. 2015;98:963–76.CAS 
PubMed 
Article 

Google Scholar 
55.Lee SE, Ryu PY, Kim SY, Kim YR, Koh JT, Kim OJ, et al. Production of Vibrio vulnificus hemolysin in vivo and its pathogenic significance. Biochem Biophys Res Commun. 2004;324:86–91.56.Senoh M, Okita Y, Shinoda S, Miyoshi S. The crucial amino acid residue related to inactivation of Vibrio vulnificus hemolysin. Micro Pathog. 2008;44:78–83.CAS 
Article 

Google Scholar 
57.Bröms JE, Ishikawa T, Wai SN, Sjöstedt A. A functional VipA-VipB interaction is required for the type VI secretion system activity of Vibrio cholerae O1 strain A1552. BMC Microbiol. 2013;13:1–12.Article 
CAS 

Google Scholar 
58.Vizcaino MI, Johnson WR, Kimes NE, Williams K, Torralba M, Nelson KE, et al. Antimicrobial resistance of the coral pathogen Vibrio coralliilyticus and Caribbean sister phylotypes isolated from a diseased octocoral. Micro Ecol. 2010;59:646–57.Article 

Google Scholar 
59.Ritchie KB. Regulation of microbial populations by coral surface mucus and mucus-associated bacteria. Mar Ecol Prog Ser. 2006;322:1–14.CAS 
Article 

Google Scholar 
60.Nissimov J, Rosenberg E, Munn CB. Antimicrobial properties of resident coral mucus bacteria of Oculina patagonica. FEMS Microbiol Lett. 2009;292:210–5.CAS 
PubMed 
Article 

Google Scholar 
61.Shnit-Orland M, Kushmaro A. Coral mucus-associated bacteria: a possible first line of defense. FEMS Microbiol Ecol. 2009;67:371–80.CAS 
PubMed 
Article 

Google Scholar 
62.Rypien KL, Ward JR, Azam F. Antagonistic interactions among coral-associated bacteria. Environ Microbiol. 2010;12:28–39.CAS 
PubMed 
Article 

Google Scholar 
63.Alagely A, Krediet CJ, Ritchie KB, Teplitski M. Signaling-mediated cross-talk modulates swarming and biofilm formation in a coral pathogen Serratia marcescens. ISME J. 2011;5:1609–20.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
64.Stocker R, Seymour JR, Samadani A, Hunt DE, Polz MF. Rapid chemotactic response enables marine bacteria to exploit ephemeral microscale nutrient patches. Proc Natl Acad Sci U.S.A. 2008;105:4209–14.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
65.Polz MF, Hunt DE, Preheim SP, Weinreich DM. Patterns and mechanisms of genetic and phenotypic differentiation in marine microbes. Philos Trans R Soc B Biol Sci. 2006;361:2009–21.Article 

Google Scholar 
66.Taylor JR, Stocker R. Trade-offs of chemotactic foraging in turbulent water. Science. 2012;338:675–9.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
67.Krediet CJ, Ritchie KB, Cohen M, Lipp EK, Patterson Sutherland K, Teplitski M. Utilization of mucus from the coral Acropora palmata by the pathogen Serratia marcescens and by environmental and coral commensal bacteria. Appl Environ Microbiol. 2009;75:3851–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
68.Krediet CJ, Ritchie KB, Alagely A, Teplitski M. Members of native coral microbiota inhibit glycosidases and thwart colonization of coral mucus by an opportunistic pathogen. ISME J. 2013;7:980–90.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
69.Packer HL, Armitage JP. The chemokinetic and chemotactic behavior of Rhodobacter sphaeroides: two independent responses. J Bacteriol. 1994;176:206–12.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
70.Deepika D, Karmakar R, Tirumkudulu MS, Venkatesh KV. Variation in swimming speed of Escherichia coli in response to attractant. Arch Microbiol. 2015;197:211–22.CAS 
PubMed 
Article 

Google Scholar 
71.Zhulin IB, Armitage JP. Motility, chemokinesis, and methylation-independent chemotaxis in Azospirillum brasilense. J Bacteriol. 1993;175:952–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
72.Ramos HC, Rumbo M, Sirard J-C. Bacterial flagellins: mediators of pathogenicity and host immune responses in mucosa. Trends Microbiol. 2004;12:509–17.CAS 
PubMed 
Article 

Google Scholar 
73.Reed KC, Muller EM, van Woesik R. Coral immunology and resistance to disease. Dis Aquat Organ. 2010;90:85–92.CAS 
PubMed 
Article 

Google Scholar 
74.Ushijima B, Videau P, Poscablo D, Stengel JW, Beurmann S, Burger AH, et al. Mutation of the toxR or mshA genes from Vibrio coralliilyticus strain OCN014 reduces infection of the coral Acropora cytherea. Environ Microbiol. 2016;18:4055–67.CAS 
PubMed 
Article 

Google Scholar 
75.Ushijima B, Richards GP, Watson MA, Schubiger CB, Häse CC. Factors affecting infection of corals and larval oysters by Vibrio coralliilyticus. PLoS ONE. 2018;13:e0199475.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
76.Peterson KM, Mekalanos JJ. Characterization of the Vibrio cholerae ToxR regulon: identification of novel genes involved in intestinal colonization. Infect Immun. 1988;56:2822–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
77.Provenzano D, Klose KE. Altered expression of the ToxR-regulated porins OmpU and OmpT diminishes Vibrio cholerae bile resistance, virulence factor expression, and intestinal colonization. Proc Natl Acad Sci U.S.A. 2000;97:10220–4.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
78.Waters CM, Bassler BL. The Vibrio harveyi quorum-sensing system uses shared regulatory components to discriminate between multiple autoinducers. Genes Dev. 2006;20:2754–67.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
79.Mukherjee S, Bassler BL. Bacterial quorum sensing in complex and dynamically changing environments. Nat Rev Microbiol. 2019;17:371–82.80.Sikora AE, Zielke RA, Lawrence DA, Andrews PC, Sandkvist M. Proteomic analysis of the Vibrio cholerae type II secretome reveals new proteins, including three related serine proteases. J Biol Chem. 2011;286:16555–66.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
81.Korotkov KV, Sandkvist M, Hol WGJ. The type II secretion system: biogenesis, molecular architecture and mechanism. Nat Rev Microbiol. 2012;10:336–51.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
82.Stathopoulos C, Hendrixson DR, Thanassi DG, Hultgren SJ, St. Geme III JW, Curtiss III R. Secretion of virulence determinants by the general secretory pathway in Gram-negative pathogens: an evolving story. Microbes Infect. 2000;2:1061–72.83.Hood RD, Singh P, Hsu FS, Güvener T, Carl MA, Trinidad RRS, et al. A Type VI secretion system of Pseudomonas aeruginosa targets a toxin to bacteria. Cell Host Microbe. 2010;7:25–37.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
84.Zheng J, Ho B, Mekalanos JJ. Genetic analysis of anti-amoebae and anti-bacterial activities of the Type VI secretion system in Vibrio cholerae. PLoS ONE. 2011;6:e23876.85.MacIntyre DL, Miyata ST, Kitaoka M, Pukatzki S. The Vibrio cholerae type VI secretion system displays antimicrobial properties. Proc Natl Acad Sci U.S.A. 2010;107:19520–4.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
86.Lee SH, Hava DL, Waldor MK, Camilli A. Regulation and temporal expression patterns of Vibrio cholerae virulence genes during infection. Cell. 1999;99:625–34.87.Pennetzdorfer N, Lembke M, Pressler K, Matson JS, Reidl J, Schild S. Regulated proteolysis in Vibrio cholerae allowing rapid adaptation to stress conditions. Front Cell Infect Microbiol. 2019;9:1–9.Article 
CAS 

Google Scholar 
88.Liu R, Chen H, Zhang R, Zhou Z, Hou Z, Gao D, et al. Comparative transcriptome analysis of Vibrio splendidus JZ6 reveals the mechanism of its pathogenicity at low temperatures. Appl Environ Microbiol. 2016;82:2050–61.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
89.Hughes TP, Anderson KD, Connolly SR, Heron SF, Kerry JT, Lough JM, et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science. 2018;359:80–3.90.Vezzulli L, Previati M, Pruzzo C, Marchese A, Bourne DG, Cerrano C, et al. Vibrio infections triggering mass mortality events in a warming Mediterranean Sea. Environ Microbiol. 2010;12:2007–19.91.Zaneveld JR, Burkepile DE, Shantz AA, Pritchard CE, McMinds R, Payet JP, et al. Overfishing and nutrient pollution interact with temperature to disrupt coral reefs down to microbial scales. Nat Commun. 2016;7:1–12.Article 
CAS 

Google Scholar  More

1.Erb, K. H. et al. Unexpectedly large impact of forest management and grazing on global vegetation biomass. Nature 553, 73–76 (2018).CAS 
PubMed 
Article 

Google Scholar 
2.Luyssaert, S. et al. Old-growth forests as global carbon sinks. Nature 455, 213–215 (2008).CAS 
PubMed 
Article 

Google Scholar 
3.Drake, J. B. et al. Above-ground biomass estimation in closed canopy Neotropical forests using lidar remote sensing: factors affecting the generality of relationships. Glob. Ecol. Biogeogr. 12, 147–159 (2003).Article 

Google Scholar 
4.Lefsky, M. A. et al. Lidar remote sensing of above-ground biomass in three biomes. Glob. Ecol. Biogeogr. 11, 393–399 (2002).Article 

Google Scholar 
5.Duncanson, L. et al. The importance of consistent global forest aboveground biomass product validation. Surv. Geophys. 40, 979–999 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
6.Spawn, S. A., Sullivan, C. C., Lark, T. J. & Gibbs, H. K. Harmonized global maps of above and belowground biomass carbon density in the year 2010. Sci. Data 7, 112 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
7.Ottaviani, G. et al. The neglected belowground dimension of plant dominance. Trends Ecol. Evol. 35, 763–766 (2020).PubMed 
Article 

Google Scholar 
8.Jackson, L. E., Burger, M. & Cavagnaro, T. R. Roots, nitrogen transformations, and ecosystem services. Annu. Rev. Plant Biol. 59, 341–363 (2008).CAS 
PubMed 
Article 

Google Scholar 
9.Gill, R. A. & Jackson, R. B. Global patterns of root turnover for terrestrial ecosystems. New Phytol. 147, 13–31 (2000).Article 

Google Scholar 
10.Robinson, D. Implications of a large global root biomass for carbon sink estimates and for soil carbon dynamics. Proc. R. Soc. Lond. B 274, 2753–2759 (2007).CAS 

Google Scholar 
11.Bardgett, R. D., Mommer, L. & De Vries, F. T. Going underground: root traits as drivers of ecosystem processes. Trends Ecol. Evol. 29, 692–699 (2014).PubMed 
Article 

Google Scholar 
12.Ribeiro, S. C. et al. Above- and belowground biomass in a Brazilian Cerrado. For. Ecol. Manage. 262, 491–499 (2011).Article 

Google Scholar 
13.Mokany, K., Raison, R. J. & Prokushkin, A. S. Critical analysis of root:shoot ratios in terrestrial biomes. Glob. Chang. Biol. 12, 84–96 (2006).Article 

Google Scholar 
14.Saatchi, S. S. et al. Benchmark map of forest carbon stocks in tropical regions across three continents. Proc. Natl Acad. Sci. USA 108, 9899–9904 (2011).CAS 
PubMed 
Article 

Google Scholar 
15.Ruesch, A. S. & Gibbs, H. H. K. New IPCC Tier-1 Global Biomass Carbon Map for the Year 2000 (Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, 2008).16.Chen, J. L. & Reynolds, J. F. A coordination model of whole-plant carbon allocation in relation to water stress. Ann. Bot. 80, 45–55 (1997).CAS 
Article 

Google Scholar 
17.Franklin, O. et al. Modeling carbon allocation in trees: a search for principles. Tree Physiol. 32, 648–666 (2012).CAS 
PubMed 
Article 

Google Scholar 
18.Bloom, A. J., Chapin, F. S. & Mooney, H. A. Resource limitation in plants—an economic analogy. Annu. Rev. Ecol. Syst. 16, 363–392 (1985).Article 

Google Scholar 
19.Poorter, H. et al. Biomass allocation to leaves, stems and roots: meta-analyses of interspecific variation and environmental control. New Phytol. 193, 30–50 (2012).CAS 
PubMed 
Article 

Google Scholar 
20.Reich, P. in Plant Roots: The Hidden Half (eds. Waisel, Y. et al.) 205–220 (Marcel Dekker, 2006).21.Ledo, A. et al. Tree size and climatic water deficit control root to shoot ratio in individual trees globally. New Phytol. 217, 8–11 (2018).PubMed 
Article 

Google Scholar 
22.Qi, Y., Wei, W., Chen, C. & Chen, L. Plant root-shoot biomass allocation over diverse biomes: a global synthesis. Glob. Ecol. Conserv. 18, e00606 (2019).Article 

Google Scholar 
23.Reich, P. B. et al. Temperature drives global patterns in forest biomass distribution in leaves, stems, and roots. Proc. Natl Acad. Sci. USA 111, 13721–13726 (2014).CAS 
PubMed 
Article 

Google Scholar 
24.De Frenne, P. et al. Latitudinal gradients as natural laboratories to infer species’ responses to temperature. J. Ecol. 101, 784–795 (2013).Article 

Google Scholar 
25.Luo, Y. Terrestrial carbon-cycle feedback to climate warming. Annu. Rev. Ecol. Evol. Syst. 38, 683–712 (2007).Article 

Google Scholar 
26.Jackson, R. B. et al. A global analysis of root distributions for terrestrial biomes. Oecologia 108, 389–411 (1996).CAS 
PubMed 
Article 

Google Scholar 
27.Malhi, Y., Doughty, C. & Galbraith, D. The allocation of ecosystem net primary productivity in tropical forests. Philos. Trans. R. Soc. Lond. B 366, 3225–3245 (2011).CAS 
Article 

Google Scholar 
28.Roberts, D. R. et al. Cross-validation strategies for data with temporal, spatial, hierarchical, or phylogenetic structure. Ecography 40, 913–929 (2017).Article 

Google Scholar 
29.Cairns, M. A., Brown, S., Helmer, E. H. & Baumgardner, G. A. Root biomass allocation in the world’s upland forests. Oecologia 111, 1–11 (1997).PubMed 
Article 

Google Scholar 
30.McCarthy, M. C. & Enquist, B. J. Consistency between an allometric approach and optimal partitioning theory in global patterns of plant biomass allocation. Funct. Ecol. 21, 713–720 (2007).Article 

Google Scholar 
31.Barton, C. V. M. & Montagu, K. D. Effect of spacing and water availability on root:shoot ratio in Eucalyptus camaldulensis. For. Ecol. Manage. 221, 52–62 (2006).Article 

Google Scholar 
32.Enquist, B. J. & Niklas, K. J. Global allocation rules for patterns of biomass partitioning in seed plants. Science 295, 1517–1520 (2002).CAS 
PubMed 
Article 

Google Scholar 
33.Goward, S. N., Tucker, C. J. & Dye, D. G. North American vegetation patterns observed with the NOAA-7 advanced very high resolution radiometer. Vegetatio 64, 3–14 (1985).Article 

Google Scholar 
34.Manzoni, S., Jackson, R. B., Trofymow, J. A. & Porporato, A. The global stoichiometry of litter nitrogen mineralization. Science 321, 684–686 (2008).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
35.Kaiser, C., Franklin, O., Dieckmann, U. & Richter, A. Microbial community dynamics alleviate stoichiometric constraints during litter decay. Ecol. Lett. 17, 680–690 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
36.Jiao, F., Shi, X. R., Han, F. P. & Yuan, Z. Y. Increasing aridity, temperature and soil pH induce soil C-N-P imbalance in grasslands. Sci. Rep. 6, 19601 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
37.Sitch, S. et al. Recent trends and drivers of regional sources and sinks of carbon dioxide. Biogeosciences 12, 653–679 (2015).Article 

Google Scholar 
38.De Deyn, G. B., Cornelissen, J. H. C. & Bardgett, R. D. Plant functional traits and soil carbon sequestration in contrasting biomes. Ecol. Lett. 11, 516–531 (2008).PubMed 
Article 
PubMed Central 

Google Scholar 
39.Tjoelker, M. G., Craine, J. M., Wedin, D., Reich, P. B. & Tilman, D. Linking leaf and root trait syndromes among 39 grassland and savannah species. New Phytol. 167, 493–508 (2005).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
40.Personeni, E. & Loiseau, P. How does the nature of living and dead roots affect the residence time of carbon in the root litter continuum? Plant Soil 267, 129–141 (2004).CAS 
Article 

Google Scholar 
41.Tuanmu, M. N. & Jetz, W. A global 1-km consensus land-cover product for biodiversity and ecosystem modelling. Glob. Ecol. Biogeogr. 23, 1031–1045 (2014).Article 

Google Scholar 
42.Pan, Y., Birdsey, R. A., Phillips, O. L. & Jackson, R. B. The structure, distribution, and biomass of the world’s forests. Annu. Rev. Ecol. Evol. Syst. 44, 593–622 (2013).Article 

Google Scholar 
43.Jackson, R. B., Mooney, H. A. & Schulze, E. D. A global budget for fine root biomass, surface area, and nutrient contents. Proc. Natl Acad. Sci. USA 94, 7362–7366 (1997).CAS 
PubMed 
Article 

Google Scholar 
44.Genet, H., Bréda, N. & Dufrêne, E. Age-related variation in carbon allocation at tree and stand scales in beech (Fagus sylvatica L.) and sessile oak (Quercus petraea (Matt.) Liebl.) using a chronosequence approach. Tree Physiol. 30, 177–192 (2009).PubMed 
Article 

Google Scholar 
45.De Castro, E. A. & Kauffman, J. B. Ecosystem structure in the Brazilian Cerrado: a vegetation gradient of aboveground biomass, root mass and consumption by fire. J. Trop. Ecol. 14, 263–283 (1998).Article 

Google Scholar 
46.Ding, B. & Sun, J. Study on biomass of Korean pine plantation in east mountain areas of northeast China. Bull. Bot. Res. 9, 149–157 (1989).
Google Scholar 
47.Ding, B., Liu, S. & Cai, T. Studies on biological productivity of artificial forests of Dahurian larches. Chin. J. Plant Ecol. 14, 226–236 (1990).
Google Scholar 
48.Ding, B. & Sun, J. Accumulation and distribution of productivity and nutrient element in natural Manchurian ash. J. Northeast For. Univ. 4, 1–9 (1989).
Google Scholar 
49.Dossa, E. L., Fernandes, E. C. M., Reid, W. S. & Ezui, K. Above- and belowground biomass, nutrient and carbon stocks contrasting an open-grown and a shaded coffee plantation. Agrofor. Syst. 72, 103–115 (2008).Article 

Google Scholar 
50.Epron, D. et al. Do changes in carbon allocation account for the growth response to potassium and sodium applications in tropical Eucalyptus plantations? Tree Physiol. 32, 667–679 (2012).CAS 
PubMed 
Article 

Google Scholar 
51.Fonseca, W., Rey Benayas, J. M. & Alice, F. E. Carbon accumulation in the biomass and soil of different aged secondary forests in the humid tropics of Costa Rica. For. Ecol. Manage. 262, 1400–1408 (2011).Article 

Google Scholar 
52.Goodman, R. C. et al. Amazon palm biomass and allometry. For. Ecol. Manage. 310, 994–1004 (2013).Article 

Google Scholar 
53.Greenland, D. J. & Kowal, J. M. L. Nutrient content of the moist tropical forest of Ghana. Plant Soil 12, 154–173 (1960).CAS 
Article 

Google Scholar 
54.He, Y. et al. Carbon storage capacity of monoculture and mixed-species plantations in subtropical China. For. Ecol. Manage. 295, 193–198 (2013).Article 

Google Scholar 
55.Aiba, M. & Nakashizuka, T. Variation in juvenile survival and related physiological traits among dipterocarp species co‐existing in a Bornean forest. J. Veg. Sci. 18, 379–388 (2007).Article 

Google Scholar 
56.Jha, K. K. Carbon storage and sequestration rate assessment and allometric model development in young teak plantations of tropical moist deciduous forest, India. J. For. Res. 26, 589–604 (2015).CAS 
Article 

Google Scholar 
57.Kalita, R. M., Das, A. K. & Nath, A. J. Allometric equations for estimating above- and belowground biomass in Tea (Camellia sinensis (L.) O. Kuntze) agroforestry system of Barak Valley, Assam, northeast India. Biomass Bioenergy 83, 42–49 (2015).Article 

Google Scholar 
58.Kenzo, T. et al. Development of allometric relationships for accurate estimation of above- and below-ground biomass in tropical secondary forests in Sarawak, Malaysia. J. Trop. Ecol. 25, 371–386 (2009).Article 

Google Scholar 
59.Kenzo, T. et al. Allometric equations for accurate estimation of above-ground biomass in logged-over tropical rainforests in Sarawak, Malaysia. J. For. Res. 14, 365–372 (2009).CAS 
Article 

Google Scholar 
60.Kraenzel, M., Castillo, A., Moore, T. & Potvin, C. Carbon storage of harvest-age teak (Tectona grandis) plantations, Panama. For. Ecol. Manage. 173, 213–225 (2003).Article 

Google Scholar 
61.Kuyah, S., Dietz, J., Muthuri, C., van Noordwijk, M. & Neufeldt, H. Allometry and partitioning of above- and below-ground biomass in farmed eucalyptus species dominant in Western Kenyan agricultural landscapes. Biomass Bioenergy 55, 276–284 (2013).Article 

Google Scholar 
62.Liu, S., Cai, Y. & Cai, T. in Long-term Research on Forest Ecosystems (ed. Zhou, X.) 419–427 (Northeast Forestry Univ. Press, 1991).63.Luo, T. et al. Root biomass along subtropical to alpine gradients: global implication from Tibetan transect studies. For. Ecol. Manage. 206, 349–363 (2005).Article 

Google Scholar 
64.Markesteijn, L. & Poorter, L. Seedling root morphology and biomass allocation of 62 tropical tree species in relation to drought- and shade-tolerance. J. Ecol. 97, 311–325 (2009).Article 

Google Scholar 
65.McNicol, I. M. et al. Development of allometric models for above and belowground biomass in swidden cultivation fallows of northern Laos. For. Ecol. Manage. 357, 104–116 (2015).Article 

Google Scholar 
66.Aiba, M. & Nakashizuka, T. Sapling structure and regeneration strategy in 18 Shorea species co-occurring in a tropical rainforest. Ann. Bot. 96, 313–321 (2005).PubMed 
PubMed Central 
Article 

Google Scholar 
67.Menaut, J. C. & Cesar, J. Structure and primary productivity of Lamto savannas, Ivory Coast. Ecology 60, 1197–1210 (1979).Article 

Google Scholar 
68.Morais, V. A. et al. Estoques de carbono e biomassa de um fragmento de cerradão em Minas Gerais, Brasil. Cerne 19, 237–245 (2013).Article 

Google Scholar 
69.Mugasha, W. A. et al. Allometric models for prediction of above- and belowground biomass of trees in the miombo woodlands of Tanzania. For. Ecol. Manage. 310, 87–101 (2013).Article 

Google Scholar 
70.Návar, J. Plasticity of biomass component allocation patterns in semiarid Tamaulipan thornscrub and dry temperate pine species of northeastern Mexico. Polibotánica 31, 121–141 (2011).
Google Scholar 
71.Njana, M. A., Eid, T., Zahabu, E. & Malimbwi, R. Procedures for quantification of belowground biomass of three mangrove tree species. Wetl. Ecol. Manage. 23, 749–764 (2015).Article 

Google Scholar 
72.Nogueira Junior, L. R., Engel, V. L., Parrotta, J. A., de Melo, A. C. G. & Ré, D. S. Equações alométricas para estimativa da biomassa arbórea em plantios mistos com espécies nativas na restauração da Mata Atlântica. Biota Neotrop. 14, 1–9 (2014).Article 

Google Scholar 
73.Peichl, M. & Arain, M. A. Above- and belowground ecosystem biomass and carbon pools in an age-sequence of temperate pine plantation forests. Agric. For. Meteorol. 140, e20130084 (2006).Article 

Google Scholar 
74.Battles, J. J. et al. Vegetation composition, structure, and biomass of two unpolluted watersheds in the Cordillera de Piuchué, Chiloé Island, Chile. Plant Ecol. 158, 5–19 (2002).Article 

Google Scholar 
75.Ryan, C. M., Williams, M. & Grace, J. Above- and belowground carbon stocks in a miombo woodland landscape of Mozambique. Biotropica 43, 423–432 (2011).Article 

Google Scholar 
76.Saint-André, L. et al. Age-related equations for above- and below-ground biomass of a Eucalyptus hybrid in Congo. For. Ecol. Manage. 205, 199–214 (2005).Article 

Google Scholar 
77.Aryal, D. R., De Jong, B. H. J., Ochoa-Gaona, S., Esparza-Olguin, L. & Mendoza-Vega, J. Carbon stocks and changes in tropical secondary forests of southern Mexico. Agric. Ecosyst. Environ. 195, 220–230 (2014).Article 

Google Scholar 
78.Schepaschenko, D. et al. A dataset of forest biomass structure for Eurasia. Sci. Data 4, 170070 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
79.Schroth, G., D’Angelo, S. A., Teixeira, W. G., Haag, D. & Lieberei, R. Conversion of secondary forest into agroforestry and monoculture plantations in Amazonia: consequences for biomass, litter and soil carbon stocks after 7 years. For. Ecol. Manage. 163, 131–150 (2002).Article 

Google Scholar 
80.Schulze, E. D. et al. Rooting depth, water availability, and vegetation cover along an aridity gradient in Patagonia. Oecologia 108, 503–511 (1996).Article 

Google Scholar 
81.Stolbovoi, V. & McCallum, I. Land resources of Russia [CD] (International Institute for Applied Systems Analysis and the Russian Academy of Science, 2002); http://www.iiasa.ac.at/Research/FOR/russia_cd/guide.htm82.Wang, L. et al. Biomass allocation patterns across China’s terrestrial biomes. PLoS ONE 9, e93566 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
83.Wauters, J. B., Coudert, S., Grallien, E., Jonard, M. & Ponette, Q. Carbon stock in rubber tree plantations in Western Ghana and Mato Grosso (Brazil). For. Ecol. Manage. 255, 2347–2361 (2008).Article 

Google Scholar 
84.Williams-Linera, G. Biomass and nutrient content in two successional stages of tropical wet forest in Uxpanapa, Mexico. Biotropica 15, 275–284 (1983).Article 

Google Scholar 
85.Xu, Y. et al. Improving allometry models to estimate the above- and belowground biomass of subtropical forest, China. Ecosphere 6, 289 (2015).Article 

Google Scholar 
86.Youkhana, A. H. & Idol, T. W. Allometric models for predicting above- and belowground biomass of Leucaena-KX2 in a shaded coffee agroecosystem in Hawaii. Agrofor. Syst. 83, 331–345 (2011).Article 

Google Scholar 
87.Zhang, H. et al. Biogeographical patterns of biomass allocation in leaves, stems, and roots in China’s forests. Sci. Rep. 5, 15997 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
88.Castellanos, J., Maass, M. & Kummerow, J. Root biomass of a dry deciduous tropical forest in Mexico. Plant Soil 131, 225–228 (1991).Article 

Google Scholar 
89.Zheng, Z., Feng, Z., Cao, M., Li, Z. & Zhang, J. Forest structure and biomass of a tropical seasonal rain forest in Xishuangbanna, southwest China. Biotropica 38, 318–327 (2006).Article 

Google Scholar 
90.Návar, J. Root stock biomass and productivity assessments of reforested pine stands in northern Mexico. For. Ecol. Manage. 338, 139–147 (2015).Article 

Google Scholar 
91.Wang, X., Fang, J. & Zhu, B. Forest biomass and root–shoot allocation in northeast China. For. Ecol. Manage. 255, 4007–4020 (2008).Article 

Google Scholar 
92.Chen, D. K., Zhou, X. F., Zhao, H. X., Wang, Y. H. & Jing, Y. Y. Study on the structure, function and succession of the four types in natural secondary forest. J. Northeast For. Univ. 2, 1–20 (1982).
Google Scholar 
93.Chidumayo, E. N. Estimating tree biomass and changes in root biomass following clear-cutting of Brachystegia-Julbernardia (miombo) woodland in central Zambia. Environ. Conserv. 41, 54–63 (2014).Article 

Google Scholar 
94.Coll, L., Potvin, C., Messier, C. & Delagrange, S. Root architecture and allocation patterns of eight native tropical species with different successional status used in open-grown mixed plantations in Panama. Trees 22, 585–596 (2008).Article 

Google Scholar 
95.Das, D. K. & Chaturvedi, O. P. Structure and function of Populus deltoides agroforestry systems in eastern India: 1. dry matter dynamics. Agrofor. Syst. 65, 215–221 (2005).Article 

Google Scholar 
96.Ni, J. Estimating net primary productivity of grasslands from field biomass measurements in temperate northern China. Plant Ecol. 174, 217–234 (2011).Article 

Google Scholar 
97.Olson, R. et al. NPP Multi-Biome: Summary Data from Intensive Studies at 125 Sites, 1936–2006 (ORNL DAAC, accessed 19 June 2019); https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=135298.Perez, C. A. & Frangi, J. L. Grassland biomass dynamics along an altitudinal gradient in the pampa. J. Range Manage. 53, 518–528 (2007).Article 

Google Scholar 
99.Perez-Quezada, J. F. F., Delpiano, C. A. A., Snyder, K. A. A., Johnson, D. A. A. & Franck, N. Carbon pools in an arid shrubland in Chile under natural and afforested conditions. J. Arid Environ. 75, 29–37 (2011).Article 

Google Scholar 
100.Pornon, A., Boutin, M. & Lamaze, T. Contribution of plant species to the high N retention capacity of a subalpine meadow undergoing elevated N deposition and warming. Environ. Pollut. 245, 235–242 (2019).CAS 
PubMed 
Article 

Google Scholar 
101.Ramakrishnan, P. S. & Ram, S. C. Vegetation, biomass and productivity of seral grasslands of Cherrapunji in north-east India. Vegetatio 74, 47–53 (1988).Article 

Google Scholar 
102.Shaver, G. R., Laundre, J. A., Giblin, A. E. & Nadelhoffer, K. J. Changes in live plant biomass, primary production, and species composition along a riverside toposequence in Arctic Alaska, USA. Arct. Alp. Res. 28, 363–379 (2006).Article 

Google Scholar 
103.Smith, J. M. B. & Klinger, L. F. Aboveground:belowground phytomass ratios in Venezuelan paramo vegetation and their significance. Arct. Alp. Res. 17, 189–198 (2006).Article 

Google Scholar 
104.Sun, J. et al. Effects of grazing regimes on plant traits and soil nutrients in an alpine steppe, northern Tibetan Plateau. PLoS ONE 9, e108821 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
105.Wang, P. et al. Belowground plant biomass allocation in tundra ecosystems and its relationship with temperature. Environ. Res. Lett. 11, 055003 (2016).Article 
CAS 

Google Scholar 
106.Yang, Y., Fang, J., Ji, C. & Han, W. Above- and belowground biomass allocation in Tibetan grasslands. J. Veg. Sci. 20, 177–184 (2009).Article 

Google Scholar 
107.Yang, Y., Fang, J., Ma, W., Guo, D. & Mohammat, A. Large-scale pattern of biomass partitioning across China’s grasslands. Glob. Ecol. Biogeogr. 19, 268–277 (2010).Article 

Google Scholar 
108.Geng, H. L., Wang, Y. H., Wang, F. Y. & Jia, B. R. The dynamics of root-shoot ratio and its environmental effective factors of recovering Leymus chinensis steppe vegetation in Inner Mongolia, China. Acta Ecol. Sin. 28, 4629–4634 (2008).Article 

Google Scholar 
109.Hui, D. & Jackson, R. B. Geographical and interannual variability in biomass partitioning in grassland ecosystems: a synthesis of field data. New Phytol. 169, 85–93 (2006).CAS 
PubMed 
Article 

Google Scholar 
110.Jouquet, P., Tavernier, V., Abbadie, L. & Lepage, M. Nests of subterranean fungus-growing termites (Isoptera, Macrotermitinae) as nutrient patches for grasses in savannah ecosystems. Afr. J. Ecol. 43, 191–196 (2005).Article 

Google Scholar 
111.Leonid, U. et al. Impact of climate and grazing on biomass components of eastern Russia typical steppe. J. Integr. Agric. 13, 1183–1192 (2014).Article 

Google Scholar 
112.Lucash, M. S., Farnsworth, B. & Winner, W. E. Response of sagebrush steppe species to elevated CO2 and soil temperature. West. N. Am. Nat. 65, 80–86 (2005).
Google Scholar 
113.Luo, W. et al. Patterns of plant biomass allocation in temperate grasslands across a 2500-km transect in northern China. PLoS ONE 8, e71749 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
114.Barbour, M. G. Desert dogma reexamined: root/shoot productivity and plant spacing. Am. Midl. Nat. 89, 41–57 (1973).Article 

Google Scholar 
115.Becker, P., Sharbini, N. & Yahya, R. Root architecture and root:shoot allocation of shrubs and saplings in two lowland tropical forests: implications for life-form composition. Biotropica 31, 93–101 (1999).
Google Scholar 
116.Becker, P. & Castillo, A. Root architecture of shrubs and saplings in the understory of a tropical moist forest in lowland Panama. Biotropica 22, 242–249 (1990).Article 

Google Scholar 
117.Beier, C. et al. Carbon and nitrogen balances for six shrublands across Europe. Glob. Biogeochem. Cycles 23, GB4008 (2009).Article 
CAS 

Google Scholar 
118.Bhatt, Y. D., Rawat, Y. S. & Singh, S. P. Changes in ecosystem functioning after replacement of forest by Lantana shrubland in Kumaun Himalaya. J. Veg. Sci. 5, 67–70 (1994).Article 

Google Scholar 
119.Caldwell, M. M., White, R. S., Moore, R. T. & Camp, L. B. Carbon balance, productivity, and water use of cold-winter desert shrub communities dominated by C3 and C4 species. Oecologia 29, 275–300 (1977).PubMed 
Article 

Google Scholar 
120.De Viñas, I. C. R. et al. Biomass of root and shoot systems of Quercus coccifera shrublands in eastern Spain. Ann. For. Sci. 57, 803–810 (2000).Article 

Google Scholar 
121.Caravaca, F., Figueroa, D., Alguacil, M. M. & Roldán, A. Application of composted urban residue enhanced the performance of afforested shrub species in a degraded semiarid land. Bioresour. Technol. 90, 65–70 (2003).CAS 
PubMed 
Article 

Google Scholar 
122.Caravaca, F., Figueroa, D., Azcón-Aguilar, C., Barea, J. M. & Roldán, A. Medium-term effects of mycorrhizal inoculation and composted municipal waste addition on the establishment of two Mediterranean shrub species under semiarid field conditions. Agric. Ecosyst. Environ. 97, 95–105 (2003).Article 

Google Scholar 
123.Carrasco, L., Azcón, R., Kohler, J., Roldán, A. & Caravaca, F. Comparative effects of native filamentous and arbuscular mycorrhizal fungi in the establishment of an autochthonous, leguminous shrub growing in a metal-contaminated soil. Sci. Total Environ. 409, 1205–1209 (2011).CAS 
PubMed 
Article 

Google Scholar 
124.Carrillo-Garcia, Á., Bashan, Y. & Bethlenfalvay, G. J. Resource-island soils and the survival of the giant cactus, cardon, of Baja California Sur. Plant Soil 218, 207–214 (2000).CAS 
Article 

Google Scholar 
125.Carrión-Prieto, P. et al. Mediterranean shrublands as carbon sinks for climate change mitigation: new root-to-shoot ratios. Carbon Manage. 8, 67–77 (2017).Article 
CAS 

Google Scholar 
126.Deng, L., Han, Q. S., Zhang, C., Tang, Z. S. & Shangguan, Z. P. Above-ground and below-ground ecosystem biomass accumulation and carbon sequestration with Caragana korshinskii Kom plantation development. Land Degrad. Dev. 28, 906–917 (2017).Article 

Google Scholar 
127.Perkins, S. R. & Owens, M. K. Growth and biomass allocation of shrub and grass seedlings in response to predicted changes in precipitation seasonality. Plant Ecol. 168, 107–120 (2003).Article 

Google Scholar 
128.Gargaglione, V., Peri, P. L. & Rubio, G. Allometric relations for biomass partitioning of Nothofagus antarctica trees of different crown classes over a site quality gradient. For. Ecol. Manage. 259, 1118–1126 (2010).Article 

Google Scholar 
129.Hao, H. M. et al. Effects of shrub patch size succession on plant diversity and soil water content in the water-wind erosion crisscross region on the Loess Plateau. Catena 144, 177–183 (2016).Article 

Google Scholar 
130.Herwitz, S. R. & Olsvig-Whittaker, L. Preferential upslope growth of Zygophyllum dumosum Boiss. (Zygophyllaceae) roots into bedrock fissures in the northern Negev desert. J. Biogeogr. 16, 457–460 (1989).Article 

Google Scholar 
131.Hoffmann, A. & Kummerow, J. Root studies in the Chilean matorral. Oecologia 32, 57–69 (1978).PubMed 
Article 

Google Scholar 
132.Holl, K. D. Effects of above- and below-ground competition of shrubs and grass on Calophyllum brasiliense (Camb.) seedling growth in abandoned tropical pasture. For. Ecol. Manage. 109, 187–195 (1998).Article 

Google Scholar 
133.Hollister, R. D. & Flaherty, K. J. Above- and below-ground plant biomass response to experimental warming in northern Alaska. Appl. Veg. Sci. 13, 378–387 (2010).
Google Scholar 
134.Kizito, F. et al. Seasonal soil water variation and root patterns between two semi-arid shrubs co-existing with pearl millet in Senegal, West Africa. J. Arid Environ. 67, 436–455 (2006).Article 

Google Scholar 
135.Kummerow, J., Krause, D. & Jow, W. Root systems of chaparral shrubs. Oecologia 29, 163–177 (1977).PubMed 
Article 

Google Scholar 
136.León, M. F., Squeo, F. A., Gutiérrez, J. R. & Holmgren, M. Rapid root extension during water pulses enhances establishment of shrub seedlings in the Atacama Desert. J. Veg. Sci. 22, 120–129 (2011).Article 

Google Scholar 
137.Li, C. P. & Xiao, C. W. Above- and belowground biomass of Artemisia ordosica communities in three contrasting habitats of the Mu Us Desert, northern China. J. Arid Environ. 70, 195–207 (2007).Article 

Google Scholar 
138.Liang, Y. M., Hazlett, D. L. & Lauenroth, W. K. Biomass dynamics and water use efficiencies of five plant communities in the shortgrass steppe. Oecologia 80, 148–153 (1989).CAS 
PubMed 
Article 

Google Scholar 
139.Zan, Q., Wang, Y., Liao, B. & Zheng, D. Biomass and net productivity of Sonneratia apetala, S. caseolaris mangrove man-made forest. Wuhan Bot. Res. 19, 391–396 (2001).
Google Scholar 
140.Liao, B., Zheng, D. & Zheng, S. Studies on the biomass of Sonneratia caseolaris stand. For. Res. 3, 47–54 (1990).
Google Scholar 
141.Lufafa, A. et al. Allometric relationships and peak-season community biomass stocks of native shrubs in Senegal’s Peanut Basin. J. Arid Environ. 73, 260–266 (2009).Article 

Google Scholar 
142.Lusk, C. H. Leaf area and growth of juvenile temperate evergreens in low light: species of contrasting shade tolerance change rank during ontogeny. Funct. Ecol. 18, 820–828 (2004).Article 

Google Scholar 
143.Marsh, A. S., Arnone, J. A., Bormann, B. T. & Gordon, J. C. The role of Equisetum in nutrient cycling in an Alaskan shrub wetland. J. Ecol. 88, 999–1011 (2000).Article 

Google Scholar 
144.Martínez, F. et al. Belowground structure and production in a Mediterranean sand dune shrub community. Plant Soil 201, 209–216 (1998).Article 

Google Scholar 
145.Marziliano, P. A. et al. Estimating belowground biomass and root/shoot ratio of Phillyrea latifolia L. in the Mediterranean forest landscapes. Ann. For. Sci. 72, 585–593 (2015).Article 

Google Scholar 
146.Mauchamp, A., Montaña, C., Lepart, J., Rambal, S. & Montana, C. Ecotone dependent recruitment of a desert shrub, Flourensia cernua, in vegetation stripes. Oikos 68, 107–116 (1993).Article 

Google Scholar 
147.Mendoza-Ponce, A. & Galicia, L. Aboveground and belowground biomass and carbon pools in highland temperate forest landscape in central Mexico. Forestry 83, 497–506 (2010).Article 

Google Scholar 
148.Miller, P. C. & Ng, E. Root:shoot biomass ratios in shrubs in southern California and central Chile. Madrono 24, 215–223 (1977).
Google Scholar 
149.Mooney, H. A. & Rundel, P. W. Nutrient relations of the evergreen shrub, Adenostoma fasciculatum, in the California chaparral. Bot. Gaz. 140, 109–113 (1979).CAS 
Article 

Google Scholar 
150.Moro, M. J., Pugnaire, F. I., Haase, P. & Puigdefábregas, J. Effect of the canopy of Retama sphaerocarpa on its understorey in a semiarid environment. Funct. Ecol. 11, 425–431 (1997).Article 

Google Scholar 
151.Negreiros, D., Fernandes, G. W., Silveira, F. A. O. & Chalub, C. Seedling growth and biomass allocation of endemic and threatened shrubs of rupestrian fields. Acta Oecol. 35, 301–310 (2009).Article 

Google Scholar 
152.Nie, X., Yang, Y., Yang, L. & Zhou, G. Above- and belowground biomass allocation in shrub biomes across the northeast Tibetan Plateau. PLoS ONE 11, e0154251 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
153.Nobel, P. S., Quero, E. & Linares, H. Root versus shoot biomass: responses to water, nitrogen, and phosphorus applications for Agave lechuguilla. Bot. Gaz. 150, 411–416 (1989).Article 

Google Scholar 
154.Pacaldo, R. S., Volk, T. A. & Briggs, R. D. Greenhouse gas potentials of shrub willow biomass crops based on below- and aboveground biomass inventory along a 19-year chronosequence. Bioenergy Res. 6, 252–262 (2013).CAS 
Article 

Google Scholar 
155.Padilla, F. M., Miranda, J. D., Jorquera, M. J. & Pugnaire, F. I. Variability in amount and frequency of water supply affects roots but not growth of arid shrubs. Plant Ecol. 204, 261–270 (2009).Article 

Google Scholar 
156.Portsmuth, A., Niinemets, Ü., Truus, L. & Pensa, M. Biomass allocation and growth rates in Pinus sylvestris are interactively modified by nitrogen and phosphorus availabilities and by tree size and age. Can. J. For. Res. 35, 2346–2359 (2005).CAS 
Article 

Google Scholar 
157.Roth, G. A., Whitford, W. G. & Steinberger, Y. Jackrabbit (Lepus californicus) herbivory changes dominance in desertified Chihuahuan Desert ecosystems. J. Arid Environ. 70, 418–426 (2007).Article 

Google Scholar 
158.Ruiz-Peinado, R., Moreno, G., Juarez, E., Montero, G. & Roig, S. The contribution of two common shrub species to aboveground and belowground carbon stock in Iberian dehesas. J. Arid Environ. 91, 22–30 (2013).Article 

Google Scholar 
159.Rundel, P. W. Biomass, productivity, and nutrient allocation in subalpine shrublands and meadows of the Emerald Lake Basin, Sequoia National Park, California. Arct. Antarct. Alp. Res. 47, 115–123 (2015).Article 

Google Scholar 
160.Millikin, C. S. & Bledsoe, C. S. Biomass and distribution of fine and coarse roots from blue oak (Quercus douglasii) trees in the northern Sierra Nevada foothills of California. Plant Soil 214, 27–38 (1999).CAS 
Article 

Google Scholar 
161.Saura-Mas, S. & Lloret, F. Adult root structure of Mediterranean shrubs: relationship with post-fire regenerative syndrome. Plant Biol. 16, 147–154 (2014).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
162.Schenk, H. J. & Mahall, B. E. Positive and negative plant interactions contribute to a north-south-patterned association between two desert shrub species. Oecologia 132, 402–410 (2002).PubMed 
Article 

Google Scholar 
163.Silva, J. S., Rego, F. C. & Martins-Loução, M. A. Belowground traits of Mediterranean woody plants in a Portuguese shrubland. Ecol. Mediterr. 28, 5–13 (2002).Article 

Google Scholar 
164.Simões, M. P., Madeira, M. & Gazarini, L. Biomass and nutrient dynamics in Mediterranean seasonal dimorphic shrubs: strategies to face environmental constraints. Plant Biosyst. 146, 500–510 (2012).
Google Scholar 
165.Tao, Y., Zhang, Y. M. & Downing, A. Similarity and difference in vegetation structure of three desert shrub communities under the same temperate climate but with different microhabitats. Bot. Stud. 54, 59 (2013).PubMed 
PubMed Central 
Article 

Google Scholar 
166.Toscano, S., Scuderi, D., Giuffrida, F. & Romano, D. Responses of Mediterranean ornamental shrubs to drought stress and recovery. Sci. Hortic. 178, 145–153 (2014).Article 

Google Scholar 
167.Trubat, R., Cortina, J. & Vilagrosa, A. Nutrient deprivation improves field performance of woody seedlings in a degraded semi-arid shrubland. Ecol. Eng. 37, 1164–1173 (2011).Article 

Google Scholar 
168.Van Wijk, M. T., Williams, M., Gough, L., Hobbie, S. E. & Shaver, G. R. Luxury consumption of soil nutrients: a possible competitive strategy in above-ground and below-ground biomass allocation and root morphology for slow-growing arctic vegetation? J. Ecol. 91, 664–676 (2003).Article 

Google Scholar 
169.Walker, L. R., Clarkson, B. D., Silvester, W. B. & Clarkson, B. R. Colonization dynamics and facilitative impacts of a nitrogen-fixing shrub in primary succession. J. Veg. Sci. 14, 277–290 (2003).Article 

Google Scholar 
170.Wang, B. & Yang, X. S. Comparison of biomass and species diversity of four typical zonal vegetations. J. Fujian Coll. For. 29, 345–350 (2009).
Google Scholar 
171.Wang, M. & Li, H. Quantitative study on the soil water dynamics of various forest plantations in the Loess Plateau region in northwestern Shanxi. Acta Ecol. Sin. 2, 178–184 (1995).
Google Scholar 
172.Wang, P. et al. Seasonal changes and vertical distribution of root standing biomass of graminoids and shrubs at a Siberian tundra site. Plant Soil 407, 55–65 (2016).CAS 
Article 

Google Scholar 
173.Whittaker, R. H. & Woodwell, G. M. Dimension and production relations of trees and shrubs in the Brookhaven Forest, New York. J. Ecol. 56, 1–25 (1968).Article 

Google Scholar 
174.Xu, H., Li, Y., Xu, G. & Zou, T. Ecophysiological response and morphological adjustment of two Central Asian desert shrubs towards variation in summer precipitation. Plant Cell Environ. 30, 399–409 (2007).CAS 
PubMed 
Article 

Google Scholar 
175.Yan, Z. Biomass and its allocation in a 28-year-old Castanopsis kawakamii plantation. J. Fujian Coll. For. 2, 114–118 (1996).
Google Scholar 
176.Gong, Y. et al. Carbon storage and vertical distribution in three shrubland communities in Gurbantünggüt Desert, Uygur Autonomous Region of Xinjiang, northwest China. Chin. Geogr. Sci. 22, 541–549 (2012).Article 

Google Scholar 
177.Yu, Y., Shi, D., Qiuyi, J., He, L. & Cheng, G. On the biomass of secondary Schima superba forest in Hangzhou. J. Zhejiang For. Coll. 2, 157–161 (1993).
Google Scholar 
178.Kato, T. et al. Carbon dioxide exchange between the atmosphere and an alpine meadow ecosystem on the Qinghai-Tibetan Plateau, China. Agric. Meteorol. 124, 121–134 (2004).Article 

Google Scholar 
179.Li, Z., Zhu, Q. & Li, J. A comparison of photosynthetic carbon sequestration of four shrubs in Ningxia. Pratacultural Sci. 29, 352–357 (2012).CAS 

Google Scholar 
180.Zhu, X., Shi, Q. & Li, Y. A preliminary study on the Qinghai’s treasure house of the forest biomass and shrubs. Sci. Technol. Qinghai Agric. For. 1, 15–20 (1993).
Google Scholar 
181.Liao, B. & Zheng, D. Study on the forest biomass and productivity of olive wood. For. Res. 4, 22–29 (1991).
Google Scholar 
182.Liu, B., Liu, Z., Lü, X., Maestre, F. T. & Wang, L. Sand burial compensates for the negative effects of erosion on the dune-building shrub Artemisia wudanica. Plant Soil 374, 263–273 (2014).CAS 
Article 

Google Scholar 
183.Alguacil, M. M., Hernández, J. A., Caravaca, F., Portillo, B. & Roldán, A. Antioxidant enzyme activities in shoots from three mycorrhizal shrub species afforested in a degraded semi-arid soil. Physiol. Plant. 118, 562–570 (2003).CAS 
Article 

Google Scholar 
184.Axe, M. S., Grange, I. D. & Conway, J. S. Carbon storage in hedge biomass—a case study of actively managed hedges in England. Agric. Ecosyst. Environ. 250, 81–88 (2017).Article 

Google Scholar 
185.van den Hoogen, J. et al. Soil nematode abundance and functional group composition at a global scale. Nature 572, 194–198 (2019).PubMed 
Article 
CAS 

Google Scholar 
186.Erin, L. et al. h2o: R Interface for the ‘H2O’ Scalable Machine Learning Platform. R package v.3.32.0.2 (2020); https://github.com/h2oai/h2o-3187.Sagi, O. & Rokach, L. Ensemble learning: a survey. WIREs Data Min. Knowl. Discov. 8, e1249 (2018).
Google Scholar 
188.R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019).189.Gorelick, N. et al. Google Earth Engine: planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 202, 18–27 (2017).Article 

Google Scholar 
190.Heiberger, R. M. HH: Statistical Analysis and Data Display: Heiberger and Holland (2020).191.Hothorn, T. & Zeileis, A. partykit: A modular toolkit for recursive partytioning in R. J. Mach. Learn. Res. 16, 3905–3909 (2015).
Google Scholar 
192.Borkovec, M. & Madin, N. ggparty: ‘ggplot’ visualizations for the ‘partykit’ package (2019).193.Dormann, C. F. Effects of incorporating spatial autocorrelation into the analysis of species distribution data. Glob. Ecol. Biogeogr. 16, 129–138 (2007).Article 

Google Scholar 
194.Hutchinson, M., Xu, T., Houlder, D., Nix, H. & McMahon, J. ANUCLIM 6.0 User’s Guide (Australian National Univ., 2009).195.Fick, S. E. & Hijmans, R. J. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).Article 

Google Scholar 
196.Global Aridity and PET database (CGIAR-CSI, accessed 15 May 2018); http://www.cgiarcsi.community/data/global-aridity-and-pet-database197.CIESIN Gridded Population of the World, version 4 (GPWv4): Population Density Adjusted to Match 2015 Revision UN WPP Country Totals (NASA SEDAC, 2018); https://doi.org/10.7927/H4HX19NJ198.Venter, O. et al. Global terrestrial human footprint maps for 1993 and 2009. Sci. Data 3, 160067 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
199.SoilGrids (ISRIC, accessed 15 May 2018); https://www.soilgrids.org200.Entekhabi, D. et al. The soil moisture active passive (SMAP) mission. Proc. IEEE 98, 704–716 (2010).Article 

Google Scholar 
201.Fan, Y., Li, H. & Miguez-Macho, G. Global patterns of groundwater table depth. Science 339, 940–943 (2013).CAS 
PubMed 
Article 

Google Scholar 
202.Batjes, N. H. Harmonized soil property values for broad-scale modelling (WISE30sec) with estimates of global soil carbon stocks. Geoderma 269, 61–68 (2016).CAS 
Article 

Google Scholar 
203.Schaaf, C. & Wang, Z. MCD43A1 MODIS/Terra+Aqua BRDF/Albedo Model Parameters Daily L3 Global – 500m V006 (NASA LP DAAC, 2015); https://doi.org/10.5067/MODIS/MCD43A1C.006204.Didan, K. MOD13Q1 MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid V006 (NASA LP DAAC, 2015).205.Crowther, T. W. et al. Mapping tree density at a global scale. Nature 525, 201–205 (2015).CAS 
PubMed 
Article 

Google Scholar  More