Ecology

1.
Webster, M. S., Marra, P. P., Haig, S. M., Bensch, S. & Holmes, R. T. Links between worlds: unraveling migratory connectivity. Trends Ecol. Evol. 17, 76–83 (2002).
Google Scholar 
2.
Newton, I. The migration ecology of birds (Elsevier-Academic Press, London, 2008).
Google Scholar 

3.
Schaub, M., Kania, W. & Köppen, U. Variation of primary production during winter induces synchrony in survival rates in migratory white storks Ciconia ciconia. J. Anim. Ecol. 74, 656–666 (2005).
Google Scholar 

4.
Higuchi, H. et al. Migration of Honey-buzzards Pernis apivorus based on satellite tracking. Ornithol. Sci. 4, 109–115 (2005).
Google Scholar 

5.
Takekawa, J. et al. Geographic variation in Bar-headed Geese Anser indicus: connectivity of wintering areas and breeding grounds across a broad front. Wildfowl 59, 100–123 (2009).
Google Scholar 

6.
Batbayar, N. & Lee, H. Steppe eagle migration from Mongolia to India. In Bird migration across the Himalayas: wetland functioning amidst mountains and glaciers (eds Prins, H. H. T. & Namgali, T.) 117–127 (Cambridge University Press, Cambridge, 2017).
Google Scholar 

7.
Dixon, A., Rahman, L., Sokolov, A. & Sokolov, V. A. Peregrine falcons crossing the ‘roof of the world.’ In Bird migration across the Himalayas: wetland functioning amidst mountains and glaciers (eds Prins, H. H. & Namgali, T.) 53–67 (Cambridge University Press, Cambridge, 2017).
Google Scholar 

8.
Zalles, J. I. & Bildstein, K. L. Raptor watch: a global directory of raptor migration sites (Hawk Mountain Sanctuary, Kempton, 2000).
Google Scholar 

9.
Den Besten, J. W. Migration of Steppe Eagles Aquila nipalensis and other raptors along the Himalayas past Dharamsala, India, in autumn 2001 and spring 2002. Forktail 20, 9–13 (2004).
Google Scholar 

10.
Juhant, M. A. & Bildstein, K. L. Raptor migration across and around the Himalayas. In Bird migration across the Himalayas: wetland functioning amidst mountains and glaciers, pp 98–116 (eds Prins, H. H. & Namgali, T.) (Cambridge University Press, Cambridge, 2017).
Google Scholar 

11.
Clark, N. E., Boakes, E. H., Mcgowan, P. J. K., Mace, G. M. & Fuller, R. A. Protected areas in South Asia have not prevented habitat loss: a study using historical models of land-use change. PLoS ONE 8, e65298 (2013).
ADS  PubMed  PubMed Central  CAS  Google Scholar 

12.
Malakoff, D., Wigginton, N. S., Fahrenkamp-Uppenbrink, J. & Wible, B. Rise of the urban planet. Science 80, 272 (2016).
Google Scholar 

13.
Yasue, M., Feare, C. J., Bennun, L. & Fiedler, W. The epidemiology of H5N1 Avian influenza in wild birds: why we need better ecological data. Bioscience 56, 923–929 (2006).
Google Scholar 

14.
Yanjie, Xu., Gong, P., Wielstra, B. & Si, Y. Southward autumn migration of waterfowl facilitates cross-continental transmission of the highly pathogenic avian influenza H5N1 virus. Sci. Rep. 6, 30262 (2016).
Google Scholar 

15.
Palm, E. C. et al. Mapping migratory flyways in Asia using dynamic Brownian bridge movement models. Mov. Ecol. 3, 3 (2015).
PubMed  PubMed Central  Google Scholar 

16.
Parr, N. et al. High altitude flights by ruddy shelduck Tadorna ferruginea during trans-Himalayan migrations. J. Avian Biol. 48, 1310–1315 (2017).
Google Scholar 

17.
Galushin, V. M. A huge urban population of birds of prey in Delhi India. Ibis (Lond. 1859) 113, 522 (1971).
Google Scholar 

18.
Kumar, N., Jhala, Y. V., Qureshi, Q., Gosler, A. G. & Sergio, F. Human-attacks by an urban raptor are tied to human subsidies and religious practices. Sci. Rep. 9, 2545 (2019).
ADS  PubMed  PubMed Central  Google Scholar 

19.
Kumar, N. et al. The population density of an urban raptor is inextricably tied to human cultural practices. Proc. R. Soc. B Biol. Sci. 286, 20182932 (2019).
Google Scholar 

20.
Naoroji, R. Birds of prey of the Indian subcontinent (Christopher Helm, London, 2006).
Google Scholar 

21.
Ferguson-Lees, J. & Christie, D. A. Raptors of the world. (2001).

22.
Choudhury, A. Migration of Black-eared Kite or Large Indian Kite Milvus migrans lineatus(Gray) from Mongolia to North-Eastern India. J. Bombay Nat. Hist. Soc. 102, 229–230 (2003).
Google Scholar 

23.
Forsman, D. Identification of black-eared kite. Bird. World 16, 56–60 (2003).
Google Scholar 

24.
DeCandido, R., Subedi, T., Siponen, M., Sutasha, K. & Pierce, A. Flight identification of Milvus migrans lineatus ‘Black-eared’Kite and Milvus migrans govinda ‘Pariah’Kite in Nepal and Thailand. Bird. ASIA 20, 32–36 (2013).
Google Scholar 

25.
Scott, G. R. Elevated performance: the unique physiology of birds that fly at high altitudes. J. Exp. Biol. 214, 2455 (2011).
PubMed  CAS  Google Scholar 

26.
Sergio, F. et al. Individual improvements and selective mortality shape lifelong migratory performance. Nature 515, 410–413 (2014).
ADS  PubMed  CAS  Google Scholar 

27.
Sergio, F. et al. Migration by breeders and floaters of a long-lived raptor: implications for recruitment and territory quality. Anim. Behav. 131, 59–72 (2017).
Google Scholar 

28.
Panuccio, M., Agostini, N., Mellone, U. & Bogliani, G. Circannual variation in movement patterns of the Black Kite (Milvus migrans migrans): a review. Ethol. Ecol. Evol. 26, 1–18 (2014).
Google Scholar 

29.
Kokko, H. Competition for early arrival in migratory birds. J. Anim. Ecol. 68, 940–950 (1999).
Google Scholar 

30.
Sergio, F., Blas, J., Forero, M. G., Donazar, J. A. & Hiraldo, F. Sequential settlement and site dependence in a migratory raptor. Behav. Ecol. 18, 811–821 (2007).
Google Scholar 

31.
Bildstein, K. L. Migrating raptors of the world: their ecology & conservation (Cornell University Press, Cornell, 2006).
Google Scholar 

32.
Flack, A. et al. Costs of migratory decisions: a comparison across eight white stork populations. Science Advances 2, e1500931 (2016).
ADS  PubMed  PubMed Central  Google Scholar 

33.
Board, C. P. C. Solid waste management in slaughterhouses (Ministry of Environment and Forests, Government of India, 2004).
Google Scholar 

34.
Kumar, N. et al. Habitat selection by an avian top predator in the tropical megacity of Delhi: human activities and socio-religious practices as prey-facilitating tools. Urban Ecosyst. 21, 339–349 (2018).
Google Scholar 

35.
Meyburg, B.-U. & Meyburg, C. GPS-Satelliten-Telemetrie bei einem adulten Schwarzmilan (Milvus migrans): Aufenthaltsraum während der Brutzeit, Zug und Überwinterung. Popul. Greifvogel und Eulenarten 6, 311–352 (2009).
Google Scholar 

36.
Blanco, G. et al. Integrating population connectivity into pollution assessment: overwintering mixing reveals flame retardant contamination in breeding areas in a migratory raptor. Environ. Res. 166, 553–561 (2018).
PubMed  CAS  Google Scholar 

37.
Sergio, F. et al. No effect of satellite tagging on survival, recruitment, longevity, productivity and social dominance of a raptor, and the provisioning and condition of its offspring. J. Appl. Ecol. 52, 1665–1675 (2015).
Google Scholar 

38.
Tanferna, A., López-Jiménez, L., Blas, J., Hiraldo, F. & Sergio, F. Different location sampling frequencies by satellite tags yield different estimates of migration performance: pooling data requires a common protocol (migration estimates by satellite tracking). PLoS ONE 7, e49659 (2012).
ADS  PubMed  PubMed Central  CAS  Google Scholar 

39.
Seaman, D. E. & Powell, R. A. An evaluation of the accuracy of Kernel density estimators for home range analysis. Ecology 77, 2075–2085 (1996).
Google Scholar 

40.
Terraube, J. et al. Broad wintering range and intercontinental migratory divide within a core population of the near-threatened pallid harrier. Divers. Distrib. 18, 401–409 (2012).
Google Scholar 

41.
DeCandido, R., Gurung, S., Subedi, T. & Allen, D. The east–west migration of Steppe Eagle Aquila nipalensis and other raptors in Nepal and India. Bird ASIA 19, 18–25 (2013).
Google Scholar 

42.
Subedi, T. R. et al. Population structure and annual migration pattern of Steppe Eagles at Thoolakharka Watch Site, Nepal, 2012–2014. J. Raptor Res. 51, 165–171 (2017).
Google Scholar  More

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript. More

1.
Ockendon, N. et al. Mechanisms underpinning climatic impacts on natural populations: altered species interactions are more important than direct effects. Glob. Change Biol. 20, 2221–2229 (2014).
Article  Google Scholar 
2.
Cahill, A. E. et al. How does climate change cause extinction? Proc. R. Soc. B 280, 20121890 (2013).
Article  Google Scholar 

3.
Penuelas, J., Filella, I. & Comas, P. Changed plant and animal life cycles from 1952 to 2000 in the Mediterranean region. Glob. Change Biol. 8, 531–544 (2002).
Article  Google Scholar 

4.
Parmesan, C. et al. Poleward shifts in geographical ranges of butterfly species associated with regional warming. Nature 399, 579–583 (1999).
CAS  Article  Google Scholar 

5.
Freeman, B. G., Yaw, J. A. L., Sunday, J. M. & Hargreaves, A. L. Expanding, shifting and shrinking: the impact of global warming on species’ elevational distributions. Glob. Ecol. Biogeogr. 27, 1268–1276 (2018).
Article  Google Scholar 

6.
Draper, A. M. & Weissburg, M. J. Impacts of global warming and elevated CO2 on sensory behavior in predator–prey interactions: a review and synthesis. Front. Ecol. Evol. 7, 72 (2019).
Article  Google Scholar 

7.
Romero, G. Q. et al. Global predation pressure redistribution under future climate change. Nat. Clim. Change 8, 1087–1091 (2018).
Article  Google Scholar 

8.
Hudson, P. J., Cattadori, I. M., Boag, B. & Dobson, A. P. Climate disruption and parasite–host dynamics: patterns and processes associated with warming and the frequency of extreme climatic events. J. Helminthol. 80, 175–182 (2006).
CAS  Article  Google Scholar 

9.
Robinson, C. J., Bohannan, B. J. M. & Young, V. B. From structure to function: the ecology of host-associated microbial communities. Microbiol. Mol. Biol. Rev. 74, 453–476 (2010).
CAS  Article  Google Scholar 

10.
Longo, A. V. & Zamudio, K. R. Temperature variation, bacterial diversity and fungal infection dynamics in the amphibian skin. Mol. Ecol. 26, 4787–4797 (2017).
Article  Google Scholar 

11.
Longo, A. V., Savage, A. E., Hewson, I. & Zamudio, K. R. Seasonal and ontogenetic variation of skin microbial communities and relationships to natural disease dynamics in declining amphibians. R. Soc. Open Sci. 2, 140377 (2015).
Article  CAS  Google Scholar 

12.
Kohl, K. D. & Yahn, J. Effects of environmental temperature on the gut microbial communities of tadpoles. Environ. Microbiol. 18, 1561–1565 (2016).
Article  Google Scholar 

13.
Fontaine, S. S., Novarro, A. J. & Kohl, K. D. Environmental temperature alters the digestive performance and gut microbiota of a terrestrial amphibian. J. Exp. Biol. 221, 187559 (2018).
Article  Google Scholar 

14.
Woodhams, D. C. et al. Interacting symbionts and immunity in the amphibian skin mucosome predict disease risk and probiotic effectiveness. PLoS ONE 9, e96375 (2014).
Article  CAS  Google Scholar 

15.
Muletz-Wolz, C. R. et al. Inhibition of fungal pathogens across genotypes and temperatures by amphibian skin bacteria. Front. Microbiol. 8, 1551 (2017).
Article  Google Scholar 

16.
Bestion, E. et al. Climate warming reduces gut microbiota diversity in a vertebrate ectotherm. Nat. Ecol. Evol. 1, 0161 (2017).
Article  Google Scholar 

17.
Meyer, E. A., Cramp, R. L., Bernal, M. H. & Franklin, C. E. Changes in cutaneous microbial abundance with sloughing: possible implications for infection and disease in amphibians. Dis. Aquat. Organ. 101, 235–242 (2012).
Article  Google Scholar 

18.
Flury, S. & Gessner, M. O. Experimentally simulated global warming and nitrogen enrichment effects on microbial litter decomposers in a marsh. Appl. Environ. Microbiol. 77, 803–809 (2011).
CAS  Article  Google Scholar 

19.
Belden, L. K. & Harris, R. N. Infectious diseases in wildlife: the community ecology context. Front. Ecol. Environ. 5, 533–539 (2007).
Article  Google Scholar 

20.
Bernabé, T. N. et al. Warming weakens facilitative interactions between decomposers and detritivores, and modifies freshwater ecosystem functioning. Glob. Change Biol. 24, 3170–3186 (2018).
Article  Google Scholar 

21.
Hoekman, D. Turning up the heat: temperature influences the relative importance of top-down and bottom-up effects. Ecology 91, 2819–2825 (2010).
Article  Google Scholar 

22.
Becker, C. G. et al. Low-load pathogen spillover predicts shifts in skin microbiome and survival of a terrestrial-breeding amphibian. Proc. R. Soc. B 286, 20191114 (2019).
Article  Google Scholar 

23.
Jiménez, R. R. & Sommer, S. The amphibian microbiome: natural range of variation, pathogenic dysbiosis, and role in conservation. Biodivers. Conserv. 26, 763–786 (2017).
Article  Google Scholar 

24.
Greenspan, S. E. et al. Arthropod–bacteria interactions influence assembly of aquatic host microbiome and pathogen defense. Proc. R. Soc. B 286, 20190924 (2019).
CAS  Article  Google Scholar 

25.
Walker, W. A. in The Microbiota in Gastrointestinal Pathophysiology: Implications for Human Health, Prebiotics, Probiotics, and Dysbiosis (eds Floch, M. H. et al.) 227–232 (Academic Press, 2016).

26.
Zaneveld, J. R., McMinds, R. & Vega Thurber, R. Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat. Microbiol. 2, 17121 (2017).
CAS  Article  Google Scholar 

27.
Jin Song, S. et al. Engineering the microbiome for animal health and conservation. Exp. Biol. Med. 244, 494–504 (2019).
CAS  Article  Google Scholar 

28.
Kriss, M., Hazleton, K. Z., Nusbacher, N. M., Martin, C. G. & Lozupone, C. A. Low diversity gut microbiota dysbiosis: drivers, functional implications and recovery. Curr. Opin. Microbiol. 44, 34–40 (2018).
Article  Google Scholar 

29.
Collins, M. et al. in Climate Change 2013: The Physical Science Basis (eds. Stocker, T. F. et al.) 1029–1136 (IPCC, Cambridge Univ. Press, 2013).

30.
Sabagh, L. T., Ferreira, G. L., Branco, C. W. C., Rocha, C. F. D. & Dias, N. Y. N. Larval diet in bromeliad pools: a case study of tadpoles of two species in the genus Scinax (Hylidae). Copeia 2012, 683–689 (2012).
Article  Google Scholar 

31.
Angilletta, M. J., Niewiarowski, P. H. & Navas, C. A. The evolution of thermal physiology in ectotherms. J. Therm. Biol. 27, 249–268 (2002).
Article  Google Scholar 

32.
Becker, C. G., Longo, A. V., Haddad, C. F. B. & Zamudio, K. R. Land cover and forest connectivity alter the interactions among host, pathogen and skin microbiome. Proc. R. Soc. B 284, 20170582 (2017).
Article  Google Scholar 

33.
Fukami, T. & Nakajima, M. Community assembly: alternative stable states or alternative transient states? Ecol. Lett. 14, 973–984 (2011).
Article  Google Scholar 

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

35.
Longo, A. V. & Zamudio, K. R. Environmental fluctuations and host skin bacteria shift survival advantage between frogs and their fungal pathogen. ISME J. 11, 349–361 (2017).
Article  Google Scholar 

36.
Friedman, B. A. & Dugan, P. R. Identification of Zoogloea species and the relationship to zoogloeal matrix and floc formation. J. Bacteriol. 95, 1903–1909 (1968).
CAS  Article  Google Scholar 

37.
Gao, N. et al. Both widespread PEP-CTERM proteins and exopolysaccharides are required for floc formation of Zoogloea resiniphila and other activated sludge bacteria. Environ. Microbiol. 20, 1677–1692 (2018).
CAS  Article  Google Scholar 

38.
Merritt, R. W., Dadd, R. H. & Walker, E. D. Feeding behavior, natural food, and nutritional relationships of larval mosquitos. Annu. Rev. Entomol. 37, 349–376 (1992).
CAS  Article  Google Scholar 

39.
Moghadam, F. S. & Zimmer, M. Effects of warming and nutrient enrichment on how grazing pressure affects leaf litter–colonizing bacteria. J. Environ. Qual. 43, 851–858 (2014).
Article  CAS  Google Scholar 

40.
Zander, A., Bersier, L. & Gray, S. M. Effects of temperature variability on community structure in a natural microbial food web. Glob. Change Biol. 23, 56–67 (2017).
Article  Google Scholar 

41.
Aguirre, A. A. & Tabor, G. M. Global factors driving emerging infectious diseases: Impact on wildlife populations. Ann. NY Acad. Sci. 1149, 1–3 (2008).
Article  Google Scholar 

42.
Greenspan, S. E. et al. Infection increases vulnerability to climate change via effects on host thermal tolerance. Sci. Rep. 7, 9349 (2017).
Article  CAS  Google Scholar 

43.
Neely, W. J. et al. Synergistic effects of warming and disease linked to high mortality in cool-adapted terrestrial frogs. Biol. Conserv. 245, 108521 (2020).
Article  Google Scholar 

44.
Raffel, T. R. et al. Disease and thermal acclimation in a more variable and unpredictable climate. Nat. Clim. Change 3, 146–151 (2013).
Article  Google Scholar 

45.
Raffel, T. R., Halstead, N. T., Mcmahon, T. A., Davis, A. K. & Rohr, J. R. Temperature variability and moisture synergistically interact to exacerbate an epizootic disease. Proc. R. Soc. B 282, 20142039 (2015).
Article  Google Scholar 

46.
Greenspan, S. E. et al. White blood cell profiles in amphibians help to explain disease susceptibility following temperature shifts. Dev. Comp. Immunol. 77, 280–286 (2017).
CAS  Article  Google Scholar 

47.
Dézerald, O. et al. Food-web structure in relation to environmental gradients and predator–prey ratios in tank-bromeliad ecosystems. PLoS ONE 8, e71735 (2013).
Article  CAS  Google Scholar 

48.
Kitching, R. L. Food Webs and Container Habitats: The Natural History and Ecology of Phytotelmata (Cambridge Univ. Press, 2000).

49.
Richardson, B. A. The bromeliad microcosm and the assessment of faunal diversity in a Neotropical forest. Biotropica 31, 321–336 (1999).
Article  Google Scholar 

50.
Leroy, C. et al. What drives detrital decomposition in Neotropical tank bromeliads? Hydrobiologia 802, 85–95 (2017).
Article  Google Scholar 

51.
Giongo, A. et al. Seasonal physiological parameters and phytotelmata bacterial diversity of two bromeliad species (Aechmea gamosepala and Vriesea platynema) from the Atlantic Forest of Southern Brazil. Diversity 11, 111 (2019).
CAS  Article  Google Scholar 

52.
Frank, J. H. & Lounibos, L. P. Insects and allies associated with bromeliads: a review. Terr. Arthropod Rev. 1, 125–153 (2009).
CAS  Article  Google Scholar 

53.
Ruano-Fajardo, G., Toledo, L. F. & Mott, T. Jumping into a trap: high prevalence of chytrid fungus in the preferred microhabitats of a bromeliad-specialist frog. Dis. Aquat. Organ. 121, 223–232 (2016).
Article  Google Scholar 

54.
Gomez-Hoyos, D. A. et al. Phytotelmata selection by anurans and implications for their conservation at Las Tablas Protected Zone, Costa Rica. Alytes 35, 1–11 (2018).
Google Scholar 

55.
Haddad, C. F. B. et al. Guia dos Anfíbios da Mata Atlântica: Diversidade e Biologia (Anolis Books, 2013).

56.
Gosner, K. L. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica 16, 183–190 (1960).
Google Scholar 

57.
Oksanen, J. et al. vegan: community ecology package. R package version 2.4-1 https://CRAN.R-project.org/package=vegan (2016).

58.
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2018).

59.
Caporaso, J. G. et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624 (2012).
CAS  Article  Google Scholar 

60.
Caporaso, J. G. et al. EMP 16S Illumina Amplicon Protocol (protocols.io, 2018).

61.
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K. & Schloss, P. D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120 (2013).
CAS  Article  Google Scholar 

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

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

64.
Amir, A. et al. Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems 2, e00191–16 (2017).
Google Scholar 

65.
Katoh, K., Misawa, K., Kuma, K. & Miyata, T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 30, 3059–3066 (2002).
CAS  Article  Google Scholar 

66.
Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS ONE 5, e9490 (2010).
Article  CAS  Google Scholar 

67.
Bokulich, N. A. et al. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 6, 90 (2018).
Article  Google Scholar 

68.
DeSantis, T. Z. et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072 (2006).
CAS  Article  Google Scholar 

69.
McDonald, D. et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 6, 610–618 (2012).
CAS  Article  Google Scholar 

70.
Bokulich, N. A. et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 10, 57–59 (2013).
CAS  Article  Google Scholar 

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

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

73.
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  Article  Google Scholar 

74.
Lozupone, C., Lladser, M. E., Knights, D., Stombaugh, J. & Knight, R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 5, 169–172 (2011).
Article  Google Scholar 

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

76.
Deegan, J. On the occurrence of standardized regression coefficients greater than one. Educ. Psychol. Meas. 38, 873–888 (1978).
Article  Google Scholar 

77.
JMP v.14.0.0 (SAS Institute, 2019).

78.
Warnes, G. et al. gplots: various R programming tools for plotting data. R package version 3.0.1.1 https://CRAN.R-project.org/package=gplots (2019). More

1.
Paustian, K. et al. Climate-smart soils. Nature 532, 49–57. https://doi.org/10.1038/nature17174 (2016).
ADS  CAS  Article  PubMed  Google Scholar 
2.
Knicker, H. Soil organic N-An under-rated player for C sequestration in soils?. Soil Biol. Biochem. 43, 1118–1129 (2011).
CAS  Article  Google Scholar 

3.
Schmidt, M. W. I. et al. Persistence of soil organic matter as an ecosystem property. Nature 478, 49–56 (2011).
ADS  CAS  Article  Google Scholar 

4.
Lehmann, J. & Kleber, M. The contentious nature of soil organic matter. Nature 528, 60–68. https://doi.org/10.1038/nature16069 (2015).
ADS  CAS  Article  PubMed  Google Scholar 

5.
Cotrufo, M. F. et al. Formation of soil organic matter via biochemical and physical pathways of litter mass loss. Nat. Geosci. 8, 776. https://doi.org/10.1038/ngeo2520 (2015).
ADS  CAS  Article  Google Scholar 

6.
Liang, C., Schimel, J. P. & Jastrow, J. D. The importance of anabolism in microbial control over soil carbon storage. Nat. Microb. https://doi.org/10.1038/nmicrobiol.2017.105 (2017).
Article  Google Scholar 

7.
Czaban, W., Rasmussen, J., Nicolaisen, M. & Fomsgaard, I. S. Dissipation kinetics of asparagine in soil measured by compound-specific analysis with metabolite tracking. Biol. Fertil. Soils 52, 911–916. https://doi.org/10.1007/s00374-016-1132-6 (2016).
CAS  Article  Google Scholar 

8.
Farrell, M. et al. Rapid peptide metabolism: a major component of soil nitrogen cycling?. Global Biogeochem. Cycles 25, 107760 (2011).
Article  Google Scholar 

9.
Apostel, C., Dippold, M. A., Bore, E. & Kuzyakov, Y. Sorption of Alanine changes microbial metabolism in addition to availability. Geoderma 292, 128–134. https://doi.org/10.1016/j.geoderma2017.01.016 (2017).
ADS  CAS  Article  Google Scholar 

10.
Khan, K. S., Mack, R., Castillo, X., Kaiser, M. & Joergensen, R. G. Microbial biomass, fungal and bacterial residues, and their relationships to the soil organic matter C/N/P/S ratios. Geoderma 271, 115–123. https://doi.org/10.1016/j.geoderma.2016.02.019 (2016).
ADS  CAS  Article  Google Scholar 

11
Miltner, A., Bombach, P., Schmidt-Brücken, B. & Kästner, M. SOM genesis: microbial biomass as a significant source. Biogeochemistry 111, 1–15 (2012).
Article  Google Scholar 

12.
Lavallee, J. M., Soong, J. L. & Cotrufo, M. F. Conceptualizing soil organic matter into particulate and mineral-associated forms to address global change in the 21st century. Glob. Change Biol. 26, 261–273. https://doi.org/10.1111/gcb.14859 (2020).
ADS  Article  Google Scholar 

13
Cotrufo, M. F., Ranalli, M. G., Haddix, M. L., Six, J. & Lugato, E. Soil carbon storage informed by particulate and mineral-associated organic matter. Nat. Geosci. 12, 989. https://doi.org/10.1038/s41561-019-0484-6 (2019).
ADS  CAS  Article  Google Scholar 

14.
Nannipieri, P. & Eldor, P. The chemical and functional characterization of soil N and its biotic components. Soil Biol. Biochem. 41, 2357–2369 (2009).
CAS  Article  Google Scholar 

15.
Schimel, J. P. & Bennett, J. Nitrogen mineralization: challenges of a changing paradigm. Ecology 85, 591–602 (2004).
Article  Google Scholar 

16.
Jan, M. T., Roberts, P., Tonheim, S. K. & Jones, D. L. Protein breakdown represents a major bottleneck in nitrogen cycling in grassland soils. Soil Biol. Biochem. 41, 2272–2282. https://doi.org/10.1016/j.soilbio.2009.08.013 (2009).
CAS  Article  Google Scholar 

17.
Jones, D. L. & Kielland, K. Amino acid, peptide and protein mineralization dynamics in a taiga forest soil. Soil Biol. Biochem. 55, 60–69 (2012).
CAS  Article  Google Scholar 

18.
Mariano, E., Jones, D. L., Hill, P. W. & Trivelin, P. C. O. Mineralisation and sorption of dissolved organic nitrogen compounds in litter and soil from sugarcane fields. Soil Biol. Biochem. 103, 522–532. https://doi.org/10.1016/j.soilbio.2016.10.004 (2016).
CAS  Article  Google Scholar 

19.
Enggrob, K. L., Larsen, T., Larsen, M., Elsgaard, L. & Rasmussen, J. The influence of hydrolysis and derivatization on the determination of amino acid content and isotopic ratios in dual labeled (13C, 15N) white clover. Rapid Commun Mass Spectrom https://doi.org/10.1002/rcm.8300 (2018).
Article  Google Scholar 

20.
Rubaek, G. H. In Long-term field experiments—a unique research platform (eds Christensen, B. T. et al.) (Aarhus University, Denmark, 2008).
Google Scholar 

21.
Christensen, B. T., Petersen, J. & Trentemoller, U. M. The Askov long-term experiments on animal manure and mineral fertilizers: The Lermarken site 1894–2004. Report No. 121 (2006).

22.
Dumont, M. G. & Murrell, J. C. Stable isotope probing—linking microbial identity to function. Nat. Rev. Microbiol. 3, 499–504. https://doi.org/10.1038/nrmicro1162 (2005).
CAS  Article  PubMed  Google Scholar 

23.
Jones, D. L. & Hodge, A. Biodegradation kinetics and sorption reactions of three differently charged amino acids in soil and their effects on plant organic nitrogen availability. Soil Biol. Biochem. 31, 1331–1342 (1999).
CAS  Article  Google Scholar 

24.
Gonod, L. V., Jones, D. L. & Chenu, C. Sorption regulates the fate of the amino acids lysine and leucine in soil aggregates. Eur. J. Soil Sci. 57, 320–329 (2006).
CAS  Article  Google Scholar 

25.
Kogel-Knabner, I. & Rumpel, C. In Advances in Agronomy, vol 149 (ed. Sparks, D. L.) 1–48 (Elsevier, Amsterdam, 2018).
Google Scholar 

26.
Kirchman, D. L., Newell, S. Y. & Hodson, R. E. Incorporation versus biosynthesis of leucine – implications for measuring rates of protein-synthesis and biomass production by bacteria in marine systems. Mar. Ecol. Prog. Ser. 32, 47–59. https://doi.org/10.3354/meps032047 (1986).
ADS  CAS  Article  Google Scholar 

27.
Dippold, M. A. & Kuzyakov, Y. Biogeochemical transformations of amino acids in soil assessed by position-specific labelling. Plant Soil 373, 385–401. https://doi.org/10.1007/s11104-013-1764-3 (2013).
CAS  Article  Google Scholar 

28.
Simelyte, E., Rimpilainen, M., Zhang, X. & Toivanen, P. Role of peptidoglycan subtypes in the pathogenesis of bacterial cell wall arthritis. Ann. Rheum. Dis. 62, 976–982. https://doi.org/10.1136/ard.62.10.976 (2003).
CAS  Article  PubMed  PubMed Central  Google Scholar 

29.
Schneewind, O. & Missiakas, D. M. Protein secretion and surface display in Gram-positive bacteria. Philos. Trans. R. Soc. B Biol. Sci. 367, 1123–1139. https://doi.org/10.1098/rstb.2011.0210 (2012).
CAS  Article  Google Scholar 

30.
Vollmer, W., Blanot, D. & de Pedro, M. A. Peptidoglycan structure and architecture. FEMS Microbiol. Rev. 32, 149–167. https://doi.org/10.1111/j.1574-6976.2007.00094.x (2008).
CAS  Article  PubMed  Google Scholar 

31.
Tamez-Hidalgo, P., Christensen, B. T., Lever, M. A., Elsgaard, L. & Lomstein, B. A. Endospores, prokaryotes, and microbial indicators in arable soils from three long-term experiments. Biol. Fertil. Soils 52, 101–112. https://doi.org/10.1007/s00374-015-1057-5 (2016).
CAS  Article  Google Scholar 

32.
Hemkemeyer, M., Christensen, B. T., Martens, R. & Tebbe, C. C. Soil particle size fractions harbour distinct microbial communities and differ in potential for microbial mineralisation of organic pollutants. Soil Biol. Biochem. 90, 255–265. https://doi.org/10.1016/j.soilbio.2015.08.018 (2015).
CAS  Article  Google Scholar 

33.
Suarez-Tapia, A., Thomsen, I. K., Rasmussen, J. & Christensen, B. T. Residual N effect of long-term applications of cattle slurry using winter wheat as test crop. Field Crops Res. 221, 257–264. https://doi.org/10.1016/j.fcr.2017.10.013 (2018).
Article  Google Scholar 

34.
Castellano, M. J., Mueller, K. E., Olk, D. C., Sawyer, J. E. & Six, J. Integrating plant litter quality, soil organic matter stabilization, and the carbon saturation concept. Glob. Change Biol. 21, 3200–3209. https://doi.org/10.1111/gcb.12982 (2015).
ADS  Article  Google Scholar 

35.
Jensen, J. L., Schjonning, P., Christensen, B. T. & Munkholm, L. J. Suboptimal fertilisation compromises soil physical properties of a hard-setting sandy loam. Soil Res. 55, 332–340. https://doi.org/10.1071/sr16218 (2017).
Article  Google Scholar 

36.
Madigan, M. T., Bender, K. S., Buckley, D. H., Sattley, W. M. & Stahl, D. A. Brock Biology of Microorganisms 15th edn. (Pearson Education Limited, London, 2018).
Google Scholar 

37.
Rousk, J. et al. Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 4, 1340–1351 (2010).
Article  Google Scholar 

38.
Burns, R. G. In Soil Enzymes (ed. Burns, R. G.) (Academic Press, Cambridge, 1978).
Google Scholar 

39.
Detmers, F. J. M., Kunji, E. R. S., Lanfermeijer, F. C., Poolman, B. & Konings, W. N. Kinetics and specificity of peptide uptake by the oligopeptide transport system of Lactococcus lactis. Biochemistry 37, 16671–16679. https://doi.org/10.1021/bi981712t (1998).
CAS  Article  PubMed  Google Scholar 

40.
Scherrer, R., Beaman, T. C. & Gerhardt, P. Macromolecular sieving by dormant spore of bacillus-cereus. J. Bacteriol. 108, 868–870 (1971).
CAS  Article  Google Scholar 

41.
Nannipieri, P., Trasar-Cepeda, C. & Dick, R. P. Soil enzyme activity: a brief history and biochemistry as a basis for appropriate interpretations and meta-analysis. Biol. Fertil. Soils 54, 11–19. https://doi.org/10.1007/s00374-017-1245-6 (2018).
CAS  Article  Google Scholar 

42.
Simpson, J., Warren, C. & Adams, P. Potential protease activity and organic nitrogen concentration are rapid tests and accurate indicators of N-availability in Tasmanian Eucalyptus nitens plantations. Soil Biol. Biochem. 115, 152–160. https://doi.org/10.1016/j.soilbio.2017.08.015 (2017).
CAS  Article  Google Scholar 

43.
Lipson, D. & Nasholm, T. The unexpected versatility of plants: organic nitrogen use and availability in terrestrial ecosystems. Oecologia 128, 305–316 (2001).
ADS  Article  Google Scholar 

44.
Mortensen, R. M. In Oversigt over Landsforsøgene 2016 (ed. Pedersen, J. B.) 255–257 (Aarhus, SEGES, 2016).
Google Scholar 

45.
Rasmussen, J., Gjettermann, B., Eriksen, J., Jensen, E. S. & Hogh-Jensen, H. Fate of N-15 and C-14 from labelled plant material: Recovery in perennial ryegrass-clover mixtures and in pore water of the sward. Soil Biol. Biochem. 40, 3031–3039 (2008).
CAS  Article  Google Scholar 

46.
Colas, D., Doumeng, C., Pontalier, P. Y. & Rigal, L. Green crop fractionation by twin-screw extrusion: Influence of the screw profile on alfalfa (Medicago sativa) dehydration and protein extraction. Chem. Eng. Process. 72, 1–9. https://doi.org/10.1016/j.cep.2013.05.017 (2013).
CAS  Article  Google Scholar 

47.
Petersen, S. O. & Klug, M. J. Effects of sieving, storage, and incubation-temperature on the phospholipid fatty-acid profile of a soil microbial community. Appl. Environ. Microbiol. 60, 2421–2430 (1994).
CAS  Article  Google Scholar 

48.
Petersen, S. O., Frohne, P. S. & Kennedy, A. C. Dynamics of a soil microbial community under spring wheat. Soil Sci. Soc. Am. J. 66, 826–833 (2002).
ADS  CAS  Article  Google Scholar 

49.
Fierer, N., Schimel, J. P. & Holden, P. A. Variations in microbial community composition through two soil depth profiles. Soil Biol. Biochem. 35, 167–176. https://doi.org/10.1016/s0038-0717(02)00251-1 (2003).
CAS  Article  Google Scholar 

50.
Stromberger, M. E., Keith, A. M. & Schmidt, O. Distinct microbial and faunal communities and translocated carbon in Lumbricus terrestris drilospheres. Soil Biol. Biochem. 46, 155–162. https://doi.org/10.1016/j.soilbio.2011.11.024 (2012).
CAS  Article  Google Scholar 

51.
Frostegard, A., Baath, E. & Tunlid, A. Shifts in the structure of soil microbial communities in limed forests as revealed by phospholipid fatty-acid analysis. Soil Biol. Biochem. 25, 723–730. https://doi.org/10.1016/0038-0717(93)90113-p (1993).
Article  Google Scholar  More

1.
Ramakrishna, W., Yadav, R. & Li, K. Plant growth promoting bacteria in agriculture: two sides of a coin. Appl. Soil Ecol. 138, 10–18 (2019).
Google Scholar 
2.
Yang, J. et al. Comparison of the rhizosphere soil microbial community structure and diversity between powdery mildew-infected and noninfected strawberry plants in a greenhouse by high-throughput sequencing technology. Curr. Microbiol. 77, 1724–1736 (2020).
CAS  PubMed  Google Scholar 

3.
Nimnoi, P., Pongsilp, N. & Ruanpanun, P. Monitoring the efficiency of Streptomyces galilaeus strain KPS-C004 against root knot disease and the promotion of plant growth in the plant-parasitic nematode infested soils. Biol. Control 114, 158–166 (2017).
Google Scholar 

4.
Wei, X., Wang, X., Cao, P., Gao Z., Chen, A. J. & Han, J. Microbial community changes in the rhizosphere soil of healthy and rusty Panax ginseng and discovery of pivotal fungal genera associated with rusty roots. Biomed Res. Int. 2020, 8018525 (2020).
PubMed  PubMed Central  Google Scholar 

5.
Karlidag, H., Esitken, A., Turan, M. & Sahin, F. Effects of root inoculation of plant growth promoting rhizobacteria (PGPR) on yield, growth and nutrient element contents of leaves of apple. Sci. Hortic. 114, 16–20 (2007).
CAS  Google Scholar 

6.
Turan, M. et al. Plant growth-promoting rhizobacteria improved growth, nutrient, and hormone content of cabbage (Brassica oleracea) seedlings. Turk. J. Agric. For. 38, 327–333 (2014).
ADS  CAS  Google Scholar 

7.
Mantelin, S. & Touraine, B. Plant growth-promoting bacteria and nitrate availability: impacts on root development and nitrate uptake. J. Exp. Bot. 55, 27–34 (2004).
CAS  PubMed  Google Scholar 

8.
Sivakumar, P. K., Parthasarthi, R. & Lakshmipyriya, V. P. Encapsulation of plant growth promoting inoculant in bacterial alginate beads enriched with humic acid. Int. J. Curr. Microbiol. Appl. Sci. 3, 415–422 (2014).
Google Scholar 

9.
Nimnoi, P., Pongsilp, N. & Lumyong, S. The establishment and activity of bacterial inoculant immobilized in agar compared with those of alginate and perlite after being introduced into soil. Chiang Mai J. Sci. 44, 751–767 (2017).
CAS  Google Scholar 

10.
Piromyou, P. et al. Indigenous microbial community structure in rhizosphere of Chinese kale as affected by plant growth-promoting rhizobacteria inoculation. Pedosphere 23, 577–592 (2013).
CAS  Google Scholar 

11.
Li, T. et al. Changes in soil bacterial community structure as a result of incorporation of Brassica plants compared with continuous planting eggplant and chemical disinfection in greenhouses. PLoS ONE 12, e0173923 (2017).
PubMed  PubMed Central  Google Scholar 

12.
Bouyoucos, G. J. Hydrometer method improved for making particle size analyses of soils. Agron. J. 54, 464–465 (1962).
Google Scholar 

13.
Walkley, A. & Black, I. A. An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil Sci. 37, 29–38 (1934).
ADS  CAS  Google Scholar 

14.
Bray, R. H. & Kurtz, L. T. Determination of total, organic, and available forms of phosphorus in soils. Soil Sci. 59, 39–46 (1945).
ADS  CAS  Google Scholar 

15.
Jackson, M. L. Soil Chemical Analysis (Prentice Hall, New Delhi, India, 1967).
Google Scholar 

16.
Bergersen, F. J. Measurement of nitrogen fixation by direct means. In Methods for Evaluating Biological Nitrogen Fixation (ed. Bergersen, F. J.) 65–110 (Wiley, New York, 1980).
Google Scholar 

17.
Bradstreet, R. B. Kjeldahl method for organic nitrogen. Anal. Chem. 26, 185–187 (1954).
CAS  Google Scholar 

18.
Cavell, A. J. The colorimetric determination of phosphorus in plant materials. J. Sci. Food Agric. 6, 479–480 (1955).
CAS  Google Scholar 

19.
Walinga, I., van Vark, W. & Houba, V. J. G. Plant Analysis Procedures (Wageningen Agricultural University Academic Publisher, Wageningen, 1989).
Google Scholar 

20.
Helmke, P. A. & Sparks, D. L. Lithium, sodium, potassium, rubidium and cesium, in Methods of Soil Analysis. Part 3, Chemical Methods (ed. Sparks, D. L.) 551–574 (Soil Science Society of America, Inc. and American Society of Agronomy, Inc., Madison, Wisconsin, USA, 1996).
Google Scholar 

21.
Apprill, A., McNally, S., Parsons, R. & Weber, L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 75, 129–137 (2015).
Google Scholar 

22.
Parada, A. E., Needham, D. M. & Fuhrman, J. A. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414 (2016).
CAS  PubMed  Google Scholar 

23.
Magoč, T. & Salzberg, S. L. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963 (2011).
PubMed  PubMed Central  Google Scholar 

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

25.
Haas, B. J. et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21, 494–504 (2011).
CAS  PubMed  PubMed Central  Google Scholar 

26.
Edgar, R. C. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998 (2013).
CAS  PubMed  PubMed Central  Google Scholar 

27.
Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541 (2009).
CAS  PubMed  PubMed Central  Google Scholar 

28.
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2013).
CAS  PubMed  Google Scholar 

29.
Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004).
CAS  PubMed  PubMed Central  Google Scholar 

30.
Muyzer, G., de Waal, E. C. & Uitterlinden, A. G. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl. Environ. Microbiol. 59, 695–700 (1993).
CAS  PubMed  PubMed Central  Google Scholar 

31.
R Core Team. R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna, 2013).
Google Scholar 

32.
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, New York, 2009).
Google Scholar 

33.
Nimnoi, P. & Pongsilp, N. Marine bacterial communities in the upper gulf of Thailand assessed by lllumina next-generation sequencing platform. BMC Microbiol. 20, 19 (2020).
CAS  PubMed  PubMed Central  Google Scholar 

34.
Gu, Y. et al. The effect of microbial inoculant origin on the rhizosphere bacterial community composition and plant growth-promotion. Plant Soil 452, 105–117 (2020).
CAS  Google Scholar 

35.
Schoebitz, M., López, M. D. & Roldán, A. Bioencapsulation of microbial inoculants for better soil–plant fertilization. A review. Agron. Sustain. Dev. 33, 751–765 (2013).
CAS  Google Scholar 

36.
Wu, Z., Guo, L., Qin, S. & Li, C. Encapsulation of R. planticola Rs-2 from alginate–starch–bentonite and its controlled release and swelling behavior under simulated soil conditions. J. Ind. Microbiol. Biotechnol. 39, 317–327 (2012).
CAS  PubMed  Google Scholar 

37.
García, J. A. L., Probanza, A., Ramos, B., Barriuso, J. & Mañero, F. J. G. Effects of inoculation with plant growth promoting rhizobacteria (PGPRs) and Sinorhizobium fredii on biological nitrogen fixation, nodulation and growth of Glycine max cv. Osumi. Plant Soil 267, 143–153 (2004).
Google Scholar 

38.
Glick, B. R. Plant growth-promoting bacteria: mechanisms and application. Scientifica 2012, 963401 (2012).
PubMed  PubMed Central  Google Scholar 

39.
Jiménez-Guerrero, I. et al. The Sinorhizobium (Ensifer) fredii HH103 type 3 secretion system suppresses early defense responses to effectively nodulate soybean. Mol. Plant Microbe Interact. 28, 790–799 (2015).
PubMed  Google Scholar 

40.
Temprano-Vera, F. et al. Sinorhizobium fredii strains HH103 and NGR234 form nitrogen fixing nodules with diverse wild soybeans (Glycine soja) from central China but are ineffective on northern China accessions. Front. Microbiol. 9, 2843 (2018).
PubMed  PubMed Central  Google Scholar 

41.
Rodríguez-Navarro, D. N., Temprano, F. & Orive, R. Survival of Rhizobium sp. (Hedysarum coronarium L.) on peat-based inoculants and inoculated seeds. Soil Biol. Biochem. 23, 375–379 (1991).
Google Scholar 

42.
Khavazi, K., Rejali, F., Seguin, P. & Miransari, M. Effects of carrier, sterilisation method, and incubation on survival of Bradyrhizobium japonicum in soybean (Glycine max L.) inoculants. Enzyme Microb. Technol. 41, 780–784 (2007).
CAS  Google Scholar 

43.
Daza, A. et al. Perlite as a carrier for bacterial inoculants. Soil Biol. Biochem. 32, 567–572 (2000).
CAS  Google Scholar 

44.
Bashan, Y., Hernandez, J. P., Leyva, L. A. & Bacilio, M. Alginate microbeads as inoculant carriers for plant growth-promoting bacteria. Biol. Fertil. Soils 35, 359–368 (2002).
Google Scholar 

45.
Kiran, K. K., Koteswaraiah, P. & Chandra, T. S. Production of halophilic α-amylase by immobilized cells of moderately halophilic Bacillus sp. strain TSCVKK. Br. Microbiol. Res. J. 2, 146–157 (2012).
CAS  Google Scholar 

46.
Schoebitz, M., Ceballos, C. & Ciampi, L. Effect of immobilized phosphate solubilizing bacteria on wheat growth and phosphate uptake. J. Soil Sci. Plant. Nutr. 13, 1–10 (2013).
Google Scholar 

47.
Wu, Z., He, Y., Chen, L., Han, Y. & Li, C. Characterization of Raoultella planticola Rs-2 microcapsule prepared with a blend of alginate and starch and its release behavior. Carbohydr. Polym. 110, 259–267 (2014).
CAS  PubMed  Google Scholar 

48.
He, Y. et al. Viability evaluation of alginate-encapsulated Pseudomonas putida Rs-198 under simulated salt-stress conditions and its effect on cotton growth. Eur. J. Soil Biol. 75, 135–141 (2016).
Google Scholar 

49.
Guo, L., Wu, Z., Rasool, A. & Li, C. Effects of free and encapsulated co-culture bacteria on cotton growth and soil bacterial communities. Eur. J. Soil Biol. 53, 16–22 (2012).
CAS  Google Scholar 

50.
Wu, Z., Zhao, Y., Kaleem, I. & Li, C. Preparation of calcium-alginate microcapsuled microbial fertilizer coating Klebsiella oxytoca Rs-5 and its performance under salinity stress. Eur. J. Soil Biol. 47, 152–159 (2011).
CAS  Google Scholar 

51.
Minaxi, J. S. Efficacy of rhizobacterial strains encapsulated in nontoxic biodegradable gel matrices to promote growth and yield of wheat plants. Appl. Soil Ecol. 48, 301–308 (2011).
Google Scholar 

52.
Jain, R., Saxena, J. & Sharma, V. The evaluation of free and encapsulated Aspergillus awamori for phosphate solubilization in fermentation and soil–plant system. Appl. Soil Ecol 46, 90–94 (2010).
Google Scholar 

53.
Jain, R., Saxena, J. & Sharma, V. Differential effects of immobilized and free forms of phosphate-solubilizing fungal strains on the growth and phosphorus uptake of mung bean plants. Ann. Microbiol. 64, 1523–1534 (2014).
CAS  Google Scholar 

54.
Sankaralingam, S., Shankar, T., Sendeshkannan, K., Ramasubburayan, R. & Prakash, S. Production of protease from Pseudomonas sp. by immobilization approach on different matrices. Eur. J. Appl. Sci. 4, 146–156 (2012).
CAS  Google Scholar 

55.
Ding, J. et al. Effect of 35 years inorganic fertilizer and manure amendment on structure of bacterial and archaeal communities in black soil of northeast China. Appl. Soil Ecol. 105, 187–195 (2016).
Google Scholar 

56.
Li, F., Chen, L., Zhang, J., Yin, J. & Huang, S. Bacterial community structure after long-term organic and inorganic fertilization reveals important associations between soil nutrients and specific taxa involved in nutrient transformations. Front. Microbiol. 8, 187 (2017).
PubMed  PubMed Central  Google Scholar 

57.
Ghosh, P. K. & Maiti, T. K. Structure of extracellular polysaccharides (EPS) produced by rhizobia and their functions in legume–bacteria symbiosis: a review. Achiev. Life Sci. 10, 136–143 (2016).
Google Scholar 

58.
O’Brien, F. J. M., Dumont, M. G., Webb, J. S. & Poppy, G. M. Rhizosphere bacterial communities differ according to fertilizer regimes and cabbage (Brassica oleracea var. capitata L.) harvest time, but not aphid herbivory. Front. Microbiol. 9, 1620 (2018).
PubMed  PubMed Central  Google Scholar 

59.
Zhang, X. et al. Rhizosphere microbial community structure is selected by habitat but not plant species in two tropical seagrass beds. Front. Microbiol. 11, 161 (2020).
ADS  PubMed  PubMed Central  Google Scholar 

60.
Wang, M. et al. Soil chemical property changes in eggplant/garlic relay intercropping systems under continuous cropping. PLoS ONE 9, e111040 (2014).
ADS  PubMed  PubMed Central  Google Scholar  More

1.
Siuda, K. Kleszcze Polski (Acari: Ixodida). Część II. Systematyka i Rozmieszczenie (Polskie Towarzystwo Parazytologiczne, Warsaw, 1993).
Google Scholar 
2.
Nowak-Chmura, M. & Siuda, K. Ticks of Poland. Review of contemporary issues and latest research. Ann. Parasitol. 58(3), 125–155 (2012).
PubMed  Google Scholar 

3.
Stańczak, J., Cieniuch, S., Lass, A., Biernat, B. & Racewicz, M. Detection and quantification of Anaplasma phagocytophilum and Babesia spp. in Ixodes ricinus ticks from urban and rural environment, northern Poland, by real-time polymerase chain reaction. Exp. Appl. Acarol. 66(1), 63–81. https://doi.org/10.1007/s10493-015-9887-2 (2015).
PubMed  PubMed Central  Article  Google Scholar 

4.
Filippova, N.A. Ixodid ticks of the subfamily Ixodinae. (Izdatielstwo Nauka (Fauna SSSR. Paukoobraznyje) 4, 1977).

5.
Nijhof, A. M. et al. Ticks and associated pathogens collected from domestic animals in the Netherlands. Vector. Borne. Zoonot. Dis. 7, 585–595. https://doi.org/10.1089/vbz.2007.0130 (2007).
Article  Google Scholar 

6.
Siuda, K. & Nowak-Chmura, M. Fauna of ticks (Ixodida) in the Beskid Wyspowy. In Arthropods: Human and Animal Parasites (eds Buczek, A. & Błaszak, C. Z.) 13–24 (Akapit, New York, 2011).
Google Scholar 

7.
Nowak-Chmura, M. Fauna kleszczy (Ixodida) Europy Środkowej (Wydawnictwo Naukowe Uniwersytetu Pedagogicznego, Kraków, 2013).
Google Scholar 

8.
Król, N., Obiegala, A., Pfeffer, M., Lonc, E. & Kiewra, D. Detection of selected pathogens in ticks collected from cats and dogs in the Wrocław Agglomeration, South-West Poland. Parasit. Vectors 9(1), 351. https://doi.org/10.1186/s13071-016-1632-0 (2016).
CAS  PubMed  PubMed Central  Article  Google Scholar 

9.
Siuda, K., Nowak, M., Gierczak, M., Wierzbowska, I. & Faber, M. Kleszcze (Acari: Ixodida) pasożytujące na psach i kotach domowych w Polsce. Wiad. Parazytol. 53, 155 (2007).
Google Scholar 

10.
Day, M. J. One health: The importance of companion animal vector-borne diseases. Parasit. Vectors 4, 49. https://doi.org/10.1186/1756-3305-4-492011 (2011).
PubMed  PubMed Central  Article  Google Scholar 

11.
Rijpkema, S., Golubić, D., Molkenboer, M., Verbeek-De Kruif, N. & Schellekens, J. Identification of four genomic groups of Borrelia burgdorferi sensu lato in Ixodes ricinus ticks collected in a Lyme borreliosis endemic region of northern Croatia. Exp. Appl. Acarol. 20, 23–30 (1996).
CAS  PubMed  Article  Google Scholar 

12.
Wójcik-Fatla, A., Szymańska, J., Wdowiak, L., Buczek, A. & Dutkiewicz, J. Coincidence of three pathogens (Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum and Babesia microti) in Ixodes ricinus ticks in the Lublin makroregion. Ann. Agric. Environ. Med. 16(1), 151–158 (2009).
PubMed  Google Scholar 

13.
Massung, R. F. et al. Nested PCR assay for detection of granulocytic ehrlichiae. J. Clin. Microbiol. 36(4), 1090–1095 (1998).
CAS  PubMed  PubMed Central  Article  Google Scholar 

14.
Persing, D. H. et al. Detection of Babesia microti by polymerase chain reaction. J. Clin. Microbiol. 30, 2097–2103 (1992).
CAS  PubMed  PubMed Central  Article  Google Scholar 

15.
Sroka, J., Szymańska, J. & Wójcik-Fatla, A. The occurence of Toxoplasma gondii and Borrelia burgdorferi sensu lato in Ixodes ricinus ticks from east Poland with the use of PCR. Ann. Agric. Environ. Med. 16(2), 313–319 (2009).
PubMed  Google Scholar 

16.
Claerebout, E. et al. Ticks and associated pathogens collected from dogs and cats in Belgium. Parasit. Vectors. 6, 183. https://doi.org/10.1186/1756-3305-6-183 (2013).
PubMed  PubMed Central  Article  Google Scholar 

17.
Schreiber, C. et al. Pathogens in ticks collected from dogs in Berlin/Brandenburg, Germany. Parasit. Vectors. 7, 535. https://doi.org/10.1186/s13071-014-0535-1 (2014).
PubMed  PubMed Central  Article  Google Scholar 

18.
Davies, S. et al. Prevalence of ticks and tick-borne pathogens: Babesia and Borrelia species in ticks infesting cats of Great Britain. Vet. Parastiol. 15(244), 129–135. https://doi.org/10.1016/j.vetpar.2017.07.033 (2017).
Article  Google Scholar 

19.
Geurden, T. et al. Detection of tick-borne pathogens in ticks from dogs and cats in different European countries. Ticks. Tick. Borne. Dis. 9(6), 1431–1436. https://doi.org/10.1016/j.ttbdis.2018.06.013 (2018).
PubMed  Article  Google Scholar 

20.
Zajkowska, P. Ticks (Acari:Ixodida) attacking domestic dogs in the Malopolska voivodeship, Poland in Arthropods: In the contemporary world (eds. Buczek, A. & Błaszak Cz.). 87–99 (Koliber, 2015).

21.
Zygner, W., Jaros, S. & Wędrychowicz, H. Prevalence of Babesia canis, Borrelia afzelii, and Anaplasma phagocytophilum infection in hard ticks removed from dogs in Warsaw (central Poland). Vet. Parasitol. 153(1–2), 139–142. https://doi.org/10.1016/j.vetpar.2008.01.036 (2008).
PubMed  Article  Google Scholar 

22.
Welc-Falęciak, R., Rodo, A., Siński, E. & Bajer, A. Babesia canis and other tick-borne infections in dogs in Central Poland. Vet. Parasitol. 166(3–4), 191–198. https://doi.org/10.1016/j.vetpar.2009.09.038 (2009).
PubMed  Article  Google Scholar 

23.
Kocoń, A., Nowak-Chmura, M., Kłyś, M. & Siuda K. Ticks (Acari: Ixodida) attacking domestic cats (Felis catus L.) in southern Poland in Arthropods: In Urban and Suburban Environments (eds. Buczek, A. & Błaszak Cz.). 51–61 (Koliber, 2017).

24.
Roczeń-Karczmarz, M. et al. Comparison of the occurrence of tick-borne diseases in ticks collected from vegetation and animals in the same area. Med. Weter. 74(8), 484–488. https://doi.org/10.21521/mw.6107 (2018).
Article  Google Scholar 

25.
Matysek, M. Kleszcze w Tatrach. 2(44), 46–57 (Tatry, 2013).

26.
Haitlinger, R. Przyczynek do znajomości Acarina drobnych ssaków subalpejskiej strefy Tatr. Wiad Parazytol. 26(6), 715 (1980).
Google Scholar 

27.
Siuda, K., et al. Ryzyko atakowania przez kleszcze Ixodes ricinus ludzi w atrakcyjnych turystycznie obszarach województwa małopolskiego. Streszczenia Materiałów II Konferencji “Medycyna Podróży”. Białystok 12–14.05.2005, 44 (2005).

28.
Siuda, K., Lenčáková, D., Stanko, M., Nowak, M. & Pet’ko, B. Ixodes ricinus ticks as vectors of Lyme Borreliosis in the Carpathians Region of Poland. Konferencia Slovenska Parazitologicka Spoločnosi, Zoonózy, ich pôvodcovia a vectory, Smolenice, Slovakia, 4–6.04.2005. Zbornik abstraktov, 52 (2005).

29.
Siuda, K. & Nowak, M. Tick fauna (Ixodida) of the Polish Carpathian Mountains. In 12th International Congress of Acarology, Amsterdam, The Netherlands, 21–26.08.2006. Abstract Book Bruin J., 192 (2006).

30.
Danielová, V. et al. Extension of Ixodes ricinus ticks and agents of tick-borne diseases to mountain areas in the Czech Republic. Int. J. Med. Microbiol. 296, 48–53. https://doi.org/10.1016/j.ijmm.2006.02.007 (2006).
PubMed  Article  Google Scholar 

31.
Majláth, I. & Majláthová, V. Atlas kliešťov Slovenska. 41–47 (2015).

32.
Tarageľová, V. R., Mahríková, L., Selyemová, D., Václav, R. & Derdáková, M. Natural foci of Borrelia lusitaniae in a mountain region of Central Europe. Ticks. Tick. Borne. Dis. 7(2), 350–356. https://doi.org/10.1016/j.ttbdis.2015.12.006 (2016).
PubMed  Article  Google Scholar 

33.
Asman, M. et al. Detection of protozoans Babesia microti and Toxoplasma gondii and their co-existence in ticks (Acari: Ixodida) collected in Tarnogórski district (Upper Silesia, Poland). Ann. Agric. Environ. Med. 22(1), 80–83. https://doi.org/10.5604/12321966.1141373 (2015).
PubMed  Article  Google Scholar 

34.
Pet’ko, B., et al. Epizootiological aspect of Lyme borreliosis in the city agglomerations of the Carpathian regions of Slovakia and Poland and their peripheral part in Stawonogi. Pasożyty i nosicele (eds. Buczek, A. & Błaszak Cz.). 157–164, (KGM, 2001).

35.
Asman, M. et al. Anaplasma phagocytophilum, Babesia microti, Borrelia burgdorferi sensu lato, and Toxoplasma gondii in Ixodes ricinus (Acari, Ixodida) ticks collected from Slowinski National Park (Northern Poland). J. Vector. Ecol. 42(1), 200–202. https://doi.org/10.1111/jvec.12258 (2017).
PubMed  Article  Google Scholar 

36.
Asman, M. et al. The occurrence of three tick-borne pathogens in Ixodes ricinus ticks collected from the area of the Kraków-Częstochowa Upland (Southern Poland). Acarologia. 58(4), 969–975. https://doi.org/10.24349/acarologia/20184301 (2018).
Article  Google Scholar 

37.
Asman, M., Witecka, J., Solarz, K., Zwonik, A. & Szilman, P. Occurrence of Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum and Babesia microti in Ixodes ricinus ticks collected from selected areas of Opolskie Province in south-west Poland. Ann. Agric. Environ. Med. 26(4), 544–547. https://doi.org/10.26444/aaem/110214 (2019).
CAS  PubMed  Article  Google Scholar  More

1.
Kasischke, E. S. et al. Can. J. Res. 40, 1313–1324 (2010).
Article  Google Scholar 
2.
Walker, X. J. et al. Nature 572, 520–523 (2019).
Article  Google Scholar 

3.
Smolyanitsky, V. November 2019–April 2020 Arctic Seasonal Review (WMO, 2020); https://go.nature.com/2FdYUyX

4.
Hu, Y., Fernandez-Anez, N., Smith, T. E. & Rein, G. Int. J. Wildland Fire 27, 293–312 (2018).
Article  Google Scholar 

5.
Alaska Division of Forestry AK Fire Info: Alaska Wildland Fire Information https://go.nature.com/3hal6HB (2020).

6.
Prosperi, P. et al. Clim. Change 161, 415–432 (2020).
Article  Google Scholar 

7.
Kirdyanov, A. V. et al. Environ. Res. Lett. 15, 034061 (2020).
Article  Google Scholar 

8.
Olefeldt, D. et al. Nat. Commun. 7, 13043 (2016).
Article  Google Scholar 

9.
Gibson, C. M. et al. Nat. Commun. 9, 3041 (2018).
Article  Google Scholar 

10.
Usup, A., Hashimoto, Y., Takahashi, H. & Hayasaka, H. Tropics 14, 1–19 (2004).
Article  Google Scholar 

11.
Turetsky, M. R. et al. Nat. Geosci. 13, 138–143 (2020).
Article  Google Scholar 

12.
Xu, W. et al. Remote Sens. Environ. 193, 138–149 (2017).
Article  Google Scholar 

13.
Wooster, M. J., Xu, W. & Nightingale, T. Remote Sens. Environ. 120, 236–254 (2012).
Article  Google Scholar 

14.
Kirillina, K., Shvetsov, E. G., Protopopova, V. V., Thiesmeyer, L. & Yan, W. Environ. Res. Lett. 15, 035009 (2020).
Article  Google Scholar  More

Regional observations of volatile organic compounds
Measurements were conducted in summer 2017 and spring 2018 in the area of Pisa, Italy, as part of the HOMING project (Hunting Organic Molecules In NaviGation). Measurements consisted of: (i) a pilot study (summer 2017) to investigate the volatile organic compounds (VOCs) emitted by three representative local ecosystems surrounding the bird’s home aviary; (ii) a 2-month intensive field campaign (spring 2018) at the bird’s home aviary to monitor VOCs along with meteorological variables; (iii) three flights on board of a Cessna aircraft to sample at ca. 180 m, within the birds typical flight altitude (spring 2018). Specifically, (i) was conducted to identify the chemical composition of surrounding forest sites and test whether they can be smelled distinctly from the aviary. (ii) Was performed to derive the olfactory maps developed by birds based on the assumption that birds and our analytical equipment have comparable detection thresholds. (iii) Was needed to examine any regional scale gradient used by pigeons when flying. The home aviary is operated by the University of Pisa and is located at the rural site Arnino (43°39′25.7″N 10°18′14.7″E), outside the city of Pisa (11 km North–East), close to the Tyrrenian sea (1.8 km West) and the mouth of river Arno (3.2 km North) (Fig. 1).
Figure 1

Map of the sampling area. The point designated as “home” refers to the birds aviary of Arnino (43°39′25.7″N 10°18′14.7″E) located in Tuscany in the area of Pisa-Livorno. Points on the map show the three sites used for releasing birds during a navigational experiment and the main biogenic and anthropogenic sites in the area. Photos show an aerial view of the Arnino field site during the airborne sampling (left) and a homing pigeon equipped with a GPS-tag before being released for testing its homing performance (right). Figure drawn with Igor WaveMetrics.

Full size image

Arnino houses hundreds homing pigeons (Columba livia) used for navigational experiments. Volatile organic compound mixing ratios were determined with state-of-the-art on-line and off-line analytical techniques (PTR-MS and GC–MS) and speciated in their isomeric and enantiomeric forms.
Navigational experiments entailed releasing 174 pigeons equipped with GPS loggers (N. 80/2013-A), who had never left their home before, from unfamiliar sites 50–70 km from their home (Fig. 1). For logistical reasons (power availability and shelter for technological equipment at release points), no simultaneous measurement of air composition and bird release was possible. Therefore, we used results from the most recent available flight experiments (summer 2016 and summer 2017) to determine general homing performance indicators and examined the latter in conjunction with available meteorological information, as modelled air masses trajectories.
The daily variation of volatile organic compounds at the aviary site
Meteorological parameters including, temperature, relative humidity, wind direction and wind speed were measured in May–June 2018 during the ground-based campaign at the home aviary (Supplementary Fig. 1). Air temperature and relative humidity ranged between 12–25 °C and 60–85%, respectively (Fig. 2 and Supplementary Fig. 1), showing the expected mirrored diel cycle, whereby temperature is highest and RH lowest by day. Wind direction and speed exhibited a repeating daily pattern during the 2 months of measurements (Fig. 2 and Supplementary Fig. 1).
Figure 2

Diel cycle of meteorological parameters (a,b) and measured VOCs (c–f). The volumetric mixing ratio (VMR) of the volatile organic compounds are reported as their measured mass fragment by PTR-MS: m/z 63 (c), m/z 69 (d), m/z 81 (e) and m/z 107 (f); identified respectively as dimethyl sulphide (c), isoprene (d), sum of monoterpenes (e) and sum of xylenes (f). Meteorological parameters are plotted as mean campaign values (marker), median (line) and interquartile range (shaded area). VOC box plots report mean campaign values (marker), median (line), interquartile range (box) and 10th and 90th percentiles (whiskers).

Full size image

The proximity to the coast exposes the measurement site to influence by the local sea breeze system. Air masses reaching the aviary came from inland (mean wind angle 150°) every day until midday, and thereafter between 12:00 and 20:00–21:00 (local time) the air masses were advected from the sea (mean angle 270°). The wind speed, like the wind direction, also showed a diurnal profile. Figure 2 shows the diurnal profile of some compounds measured by proton transfer reaction mass spectrometry (PTR-MS). They have contrasting diel cycles and hereafter investigated further: dimethyl sulphide (DMS, m/z 63), isoprene (m/z 69), sum of monoterpenes (MT, m/z 81) and sum of xylenes (m/z 107). The full mass list and identities are reported in Supplementary Table 4 with other compound diel variations in Supplementary Fig. 2. We refer to the abbreviations included at the beginning of the manuscript for a reader not familiar with the terms used herein.
Dimethyl sulphide is a biogenic compound emitted by phytoplankton activity in seawater21,22,23; it is insoluble in water, and so readily transfers to the air at a rate dependent on the temperature and wind speed. Our measurements show that DMS mixing ratios at the site remained low and constant until 12:00, possibly reflecting weak soil emissions from the land24, to increase to a maximum at 20:00 due to continued marine emission, until the reversal of the wind back to offshore. The maximum hourly averaged mixing ratio of DMS is 0.17 ± 0.8 ppbv, comparable to 0.12 ppbv measured by Derstroff et al.25 at a coastal site on Cyprus.
Isoprene mixing ratios showed a bimodal distribution with two maxima at 11:00–12:00, and 18:00–20:00; similar to methyl ethyl ketone (MEK) and isoprene oxidation products methyl vinyl ketone, methacrolein and isoprene peroxides (MVK + MACR + ISOPOOH) and acetic acid (Supplementary Fig. 2). At both morning and afternoon maxima, the isoprene mixing ratio reached 0.7 ± 0.6 ppbv (maximum hourly average ± 1σ standard deviation). This is 10 times higher than the value observed during wintertime at a similar rural Mediterranean site in Spain26, but lower than the value measured in summertime by Steinbacher and coauthors27 at a rural site in the Po valley, Italy. Isoprene emission from terrestrial vegetation is mainly light dependent28, therefore near-source measurements typically exhibit a broad maximum around noon when irradiation is highest29,30. The two isoprene maxima observed here are the result of the sea breeze system, which brings isoprene poor marine influenced air to the site around midday.
Monoterpenes (C10H16) are emitted from vegetation, and are to a large part responsible for the scent of forests. Their sum at the site reached 1.6 ± 1 ppbv as a maximum hourly average at 01:00. Despite monoterpene emissions being also a function of temperature they do not show the same profile as isoprene. Monoterpene mixing ratios were higher at night, decreasing slowly from 03:00 and then sharply at 08:00, and remaining low through the day before rising back to early morning levels in the late evening. While isoprene is formed in the plant and released directly, mainly in response to light, many plant species emit monoterpenes day and night from resin duct storage pools contained in the leaves or needles31. The rates of monoterpene emission therefore depend strongly on temperature, and although lower in the night, continued emission into a shallower boundary layer can generate significant mixing ratios, as has been observed elsewhere (for example in the boreal forest32). As expected, air advected from the more vegetated inland had higher monoterpene mixing ratios than air from the sea. Although, marine monoterpene and isoprene emissions have been observed previously33, in this region terrestrial vegetation emissions clearly dominate. A similar profile (higher nighttime daytime ambient mixing ratio) was reported by Davison et al.34 from the Mediterranean forest of Castelporziano, however, the average daytime and nighttime concentrations were 3 times lower than reported here.
Xylenes (Fig. 2), benzene, toluene, and sum of trimethylbenzenes (Supplementary Fig. 2) are aromatic compounds predominantly emitted from anthropogenic sources, for example fossil fuel use. Mixing ratios increased in the morning at 8:00, probably due to traffic emissions, and decreased steadily between 9:00–12:00, as the wind direction changed. They remained low in the afternoon as the wind came from the sea, reflecting the expected absence of significant sources from this sector. The mixing ratios of all compounds increased again after 17:00, slowly for benzene and toluene and faster for xylenes and trimethylbenzene, suggesting different anthropogenic sources of such species (i.e. road traffic and industries). Although evidence for biogenic aromatic species has been documented35,36, here no evidence of natural sources of aromatic compounds was found. The maximum mixing ratios for benzene, toluene, sum of xylenes and trimethylbenzene, were 0.16 ± 0.08 ppbv, 0.32 ± 0.7 ppbv, 0.22 ± 0.09 ppbv and 0.11 ± 0.9 (hourly average maximum ± 1σ standard deviation), respectively. Similar results were found for benzene and toluene at a rural Mediterranean site in Spain (0.19 and 0.41 ppbv, respectively26).
Chiral monoterpenes
Many monoterpenes exist in two chiral forms, meaning they exist in nature as two non-superimposable mirror image forms (enantiomers), often with different biological activities. Insects can perceive each enantiomer differently and each may act as entirely distinct chemical signals37. We investigated the abundance of the prevalent chiral monoterpenes at the home aviary and at three ecosystems surrounding the aviary North (lake), East (mixed forest) and South (pine forest). This was to test whether unique mixtures could be ascertained at each site which is a condition of the mosaic hypothesis. At all sites, including the aviary, the (−) configuration dominated over the (+) for the monoterpene species α-pinene, β-pinene and limonene. This has been also observed in the tropical rainforest for α-pinene38. Relative ratios for the same monoterpene enantiomeric pair do not change significantly across the investigated sites, however, their absolute concentrations and relative abundance to the other terpene species do differ (see Fig. 3).
Figure 3

Volume mixing ratio (VMR) distribution of enantiomers from three ecosystems surrounding the birds’ aviary and at the birds’ aviary. Category plots show mean measured values with their standard deviation from three days of sampling during summer 2017. Box plots show mean campaign values (marker), median (line), interquartile range (box) and 10th and 90th percentiles (whiskers) from 2 months intensive field campaign at the birds’ aviary in spring 2018 (home) and from three flights conducted at 180 m altitude over the region during the airborne campaign. The x axes shows the retention time (min) of the chiral molecules in the gas chromatography mass spectrometer (GC–MS) and shows the good resolution achieved with the method. The left y axes indicate the volume mixing ratio of all molecules, except for (−)limonene measured at the aviary, whose VMR is indicated by the right y axes.

Full size image

At the aviary, (−)limonene was the most abundant species (0.89 ± 1.7 ppbv and 2.5 ± 2.7 ppbv), followed by (−)α-pinene (0.42 ± 0.57 ppbv and 1 ± 1 ppbv), (+)limonene (0.26 ± 0.29 ppbv and 0.37 ± 0.42 ppbv) and (−)β-pinene (0.17 ± 0.45 ppbv and 0.27 ± 0.49 ppbv, values indicate mean daytime and nighttime values ± 1σ standard deviation, respectively). The lake site emitted the least monoterpenes, followed by the pine forest and the mixed forest, the latter showing the highest concentration of measured terpenes. Since (−)limonene dominates the terpene blend measured at the aviary, we can surmise that the surrounding pine forest had a strong impact on the air chemical composition at the site (Fig. 3). A similar suite of biogenic molecules was found at all locations, albeit with differing ratios. In other words, there were no clear unique chemical markers for a particular area. Therefore, the source of biogenic molecules near to the aviary has the potential to interfere with or mask, odours from more distant sites. These findings are not consistent with the original “mosaic” hypothesis in which each location in the region has a unique chemical signature that can be determined from the aviary.
Olfactory maps: spatial and temporal distribution of VOCs
Soon after fledging, the young pigeons are housed in a loft with an associated aviary, and from the end of May are allowed to perform spontaneous flights around the loft typically until October, when the season of the experiments usually ends. During the first months after fledging39 young pigeons are exposed to the changing chemical conditions at the home site and according to the Papi olfactory hypothesis, they learn to associate wind-borne odours with wind directions thereby generating an olfactory map. It is thought that this olfactory map is continually updated throughout their life depending on conditions experienced40,41. A visual representation of this spatio-temporal information can be expressed in the form of a bivariate plot of the measured volatile organic compounds. Bivariate polar plots represent the mean campaign mixing ratio of a given compound as a function of the wind direction (angle) and wind speed (radius), similarly to a wind rose they can highlight spatial gradients in the surrounding of a measurement site. Figure 4 shows that methanol (a biogenic compound and biomass burning marker) mixing ratios were higher for higher wind speed and easterly winds (from inland). However, mixing ratios were elevated at both high and low wind speeds; indicating some methanol sources were also local. Dimethyl sulphide mixing ratios were larger for higher wind speeds originating from the North–West and South–West sectors, spanning from the mouth of the river Arno to the sea. The largest mixing ratios for isoprene came from inland, in particularly from the West–North–East for higher wind speeds and the South–East for lower ones. Smaller mixing ratios are found from the sea for large wind speed (North–West), suggestive of a weaker marine isoprene emission source. Small marine isoprene emissions have been measured previously, especially in chlorophyll rich waters12,33,42. Monoterpene mixing ratios are higher at lower wind speeds, mostly for winds coming from South–East but also East–North–East and South–West. The polar graph highlights the multiple sources of monoterpenes, and that the levels at the site are also locally influenced. The anthropogenic aromatic species xylenes and trimethylbenzene were larger when the air was transported from inland, in particular from the South-East, where heavy-traffic roads such as highway E80, the road FI-PI-LI and via AureliaSS1 are located. Additional sources of aromatics appear South of the aviary for high wind conditions, especially for trimethylbenzene. Trimethylbenzene (C6H3(CH3)3), is an aromatic hydrocarbon characterized by a strong odour, which is generally isolated from the C9 fraction of aromatics during petroleum distillation. South of the aviary, in the Stagno industrial area, there is a large petroleum refinery (Fig. 1).
Figure 4

Spatial distributions of measured VOCs reported as protonated masses and identified as: methanol (m/z 33), DMS (m/z 63), isoprene (m/z 69), sum of monoterpenes (m/z 81), sum of xylenes (m/z 107), trimethylbenzene (m/z 121). Mean campaign values of VOCs (colored scale) are showed as a function of wind direction (angle) and wind speed (radius). Figures drawn with R.

Full size image

This air composition measured at the aviary is the air breathed and smelled by the birds during their first months after fledging. The birds are therefore likely to be aware of, mixed biogenic and anthropogenic sources to the South-East, a more biogenic dominated source (with differing composition) to the West-North–East, and a marine source to the West. Thus, the ground based measurement campaign mapped several chemical compounds that vary distinctly at the site according to the wind direction and time of day. Which VOCs are useful for navigation will also depend on their respective atmospheric lifetimes. Following emission to the atmosphere VOCs are oxidized mainly by OH radicals, with small contributions by O3, NO3 and Cl radicals. Methanol, the most abundant species measured, has a lifetime of 12 days (based on OH reactivity only, assuming OH radical concentration of 2 × 106 molecules cm–343). The corresponding lifetime of isoprene is 1.4 h, therefore with a typical windspeed of 2.2 m/s the concentration of isoprene will be reduced to 1/e of the initial value in 11 km, while limonene in 2 km, xylenes in 33 km and DMS in 88 km. Therefore, over distances of ca. 50–100 km, VOCs with moderate lifetimes such as DMS and the aromatic species appear capable of creating regional gradients. The more reactive VOCs are usually photochemically transformed into less reactive products which themselves, or their combination, may create regional gradients. Significantly, the three regional gradients highlighted in Fig. 4, namely DMS, the aromatics and the monoterpenes are not aligned, rather they slope in different directions. Theoretically then, by comparing the level of these compounds to that experienced at the aviary the pigeon may, with reference to the wind-odour experience at the aviary, orient homeward. For example if the DMS level is lower at the release site than at the aviary, and the bird “knows” from the olfactory map learning phase that DMS comes from the West, then the birds homeward direction will have a westerly component. Having multiple gradients available would enable the pigeon to triangulate a homeward direction.
Regional spatial gradients
To verify that the spatial gradients observed from the olfactory maps exist at a larger regional scale we measured the air composition at the altitudes and over distances typically flown by homing pigeons (10–300 m, 100 km). Results from a flight campaign conducted on 26/05/2018 and 27/05/2018, in conjunction with the ground based campaign, are depicted in Figs. 3, 5 and Supplementary Fig. 3. Figure 3 shows the speciation of chiral molecules from ground measurements at four distinct sites (including the aviary) and above the sampled region considering all the airborne samples taken with the three flights. Interestingly, the regionally predominant compound at 180 m is α-pinene, rather than limonene observed at the ground. This is because α-pinene has a longer atmospheric lifetime (2.6 h), in comparison to β-pinene (1.8 h) and limonene (49 min)43. Figure 5 shows the mixing ratio of a biogenic precursor compound (−)α-pinene and the common terpene photochemical oxidation product nopinone44 sampled at different locations above the home aviary, and the release sites. The terpene airborne concentrations compare well with the ground-based ones (higher above the mixed Mediterranean forest, followed by the pine forest, followed by the lake site, as reported in Fig. 3). However, as expected, the terpene airborne concentrations are smaller compared to ground concentrations (40–50 pptv and 80–208 pptv, respectively) consistent with the ground based emissions being progressively oxidized by OH radicals and ozone as well as being mixed with the relatively clean air in the free troposphere above. This means that the birds will encounter strong vertical gradients of primary emitted species during flight with lower concentrations aloft. Furthermore oxidation chemistry will generate entirely different species, complicating the mixture, weakening the primary chemical signal and possibly obfuscating the originally emitted olfactory signals. In Fig. 5 we observe during the afternoon flights (with onshore winds) the nopinone mixing ratios generally increase to the East. This demonstrates that monoterpenes emitted at ground level during the day and in onshore winds are being oxidized to (among other species) nopinone on the time and space scales of the pigeon release experiments. Therefore, spatial gradients of reactive VOC do exist for the primary emitted species and their oxidation products, at a spatial scale compatible with the distances used for navigation experiments.
Figure 5

Spatial distributions of measured (−)α-pinene (a) and nopinone (b) from airborne sampling during three flights conducted on 26/05/2018 and 27/05/2018. Colored scales indicate respectively volume mixing ratio and counts of chromatograms peak area. Figures drawn with Igor WaveMetrics.

Full size image

Pigeon flight tracks and air masses trajectories analysis
To test whether the intensity of atmospheric odour signals could potentially aid in the homing of pigeons, we provide here a first preliminary analysis based on real pigeon tracks and simulated air masses trajectories. We hypothesize that overall stronger gradients of indicative aerial chemicals enable better navigational performance in pigeons18. To do this we examine the GPS logged tracks of homing birds in relation to the general gradients of atmospheric chemicals determined in the previous sections, and to the origin of air at the release sites. Bird release experiments were conducted over six days in summer 2016 and summer 2017 from three release sites simultaneously (Fig. 1). Those experiments were chosen for being the most recent results available for indices analysis, and for being conducted on sunny days, when air temperature, humidity and atmospheric pressure were comparable to those encountered during VOCs measurements (Supplementary Fig. 4). Each release experiment involved 9 or 10 individual birds, released singly every 5–10 min, over approximately 1 h and involved 27–30 individuals released on the same day from three sites (N = 174). The total number of tracks obtained was 143, due to either the loss of the device by some birds or the misfunctioning of the GSM GPS logger. Bird flight tracks were used to determine: the pigeon initial flight direction, the homing capabilities en route with the homing efficiency index (HEI45) and for the whole track with a mean aggregate azimuth penalty (MAAP); see Methods and supplementary information. Due to lack of power and shelter at the release sites, chemical information of atmospheric composition was not available during those experiments, therefore a meteorological approach, based on modelled air trajectories, was developed in order to examine how air mass transport was related to the birds’ homing performance. Specifically, we generated forward and backward trajectories of the air masses for each release day, for the site and time of release. Forward trajectories analysis showed that the prevalent wind direction during the days and time of pigeons’ flights was from west (supplementary information). We cannot yet test the bouquet of potential chemical information, but concentrate here on one traceable chemical gradient: DMS (Dimethyl sulphide) emanating from the Tyrrhenian sea. DMS, a chemical compound of marine origin (see Supplementary Fig. 5), is here identified to be a suitable candidate for homing by olfaction. Figures 2 and 4 show that it was among the most abundant compounds measured during daytime at the aviary, it follows a regular pattern of emission, is known to be detectable by birds, and it is atmospherically stable enough to survive transport over longer distances, likely decreasing along the West–East direction. Air masses reaching the aviary diurnally during the ground campaign were found to spend considerable time ( > 20 h in most cases) in the marine boundary layer in the past 24 h (see table 1 supplementary material). For these reasons, we tested if more marine-influenced air masses can influence the pigeon initial orientation (Fig. 6) and the homing path en route (homing efficiency index, supplementary material). We examined where each air mass was located (during 24 h prior to the release) reported in latitude and longitude coordinates and determined the time the air masses spent in the marine boundary layer, and the time the air masses spent over non-marine areas. As shown in Fig. 4, pigeons at the home loft are exposed to DMS associated with westerly winds; easterly winds are unlikely to carry DMS. According to the olfactory navigation hypothesis5 it is expected that a low atmospheric level of DMS at the release site is likely to produce a bias towards west in the orientation of the birds, regardless of their ultimate home direction. To test this hypothesis, the deviation from west of the individual mean vectors computed on the initial part of each track (within 10 km from the release site) is tested against the ratio of the time the air masses spent over sea versus over land (ratio sea/land) and the opposite (ratio land/sea). The Spearman ranking test highlighted that the deviation from west is positively correlated to the sea/land ratio of the air masses trajectories (n = 135, S = 0.163, p  More