Ecology

Global biodiversity has been dramatically declining over the last decades1,2,3,4. The current biodiversity crisis is primarily driven by human-induced factors, the most serious of which are land-use change, habitat fragmentation, and climate change5. While global public awareness of climate change matters is high6,7, public recognition of biodiversity loss has, historically, been low8. The understanding of biodiversity concepts highly varies among countries and social groups9,10,11: in Nigeria, the biodiversity concept was known of 20.5% of non-professional Nigerians (with basic education or no formal training) while among 88.8% of professionals with tertiary education, it reached 88.8%; 60% of participants in a study in Switzerland had never heard the term biodiversity and Chinese farmers in another pilot study have never heard about biodiversity. In the European Union, the global leader of the environmental movement on both the political and discursive levels12,13, in 2018, 71% of EU citizens had heard of biodiversity, but only around 41% of these knew what biodiversity meant14. This illiteracy is a significant constraint for conservation strategies because the development and success of actions to halt and reverse biodiversity loss strongly rely on public support15.If general awareness of biodiversity loss is low, knowledge about plant diversity is even lower16. Plants have traditionally been overlooked, and expressions such as “plant blindness”, defined as a human tendency to ignore plant species17, perfectly illustrate the situation in terms of plant conservation. And yet, current estimates suggest that two out of five plant species are threatened with extinction18. Moreover, plants play a crucial role in the world ecosystems by providing habitat, shelter, oxygen, and food, including for humans19. Local community support boosts the effectiveness of biodiversity conservation actions20,21,22. However, how biodiversity is perceived and the benefits it provides to local populations have a significant influence on this support23. Therefore, stopping the loss of plant biodiversity and the impact it has on ecosystem health and human well-being must also strive to raise public awareness on the importance of plant conservation24.A big challenge, however, is to engage people with conservation. Nowadays, in a world where a large part of the human population lives in urban areas, the contact of people with nature is declining. This is a trend that will be even more accentuated in the future25. Perhaps society’s interest in plants is decreasing because of limited exposure to plants in daily lives, schools, and work. However, by critically examining our roles as plant scientists and educators, we realize that there are probably things we could, and should, do differently. New strategies to connect people to nature are required to spark people’s interest in and knowledge of plants. Citizen science programs and mobile applications (apps) are noteworthy initiatives that are helping to achieve this goal.Citizen science is defined as the general public involvement in scientific research activities and currently is a mainstream approach to collect information and data on a wide range of scientific subjects26,27. The development of mobile technologies and the widespread use of smartphones have boosted citizen science and enabled the development of mobile apps, which are digital tools that integrate, in real-time, data from multiple sources28.The goal of this article is to show how citizen science and mobile apps can be used as educational tools to raise awareness about plant biodiversity and conservation among the general public. We focused on formal education activities, at the Bachelor of Science (BSc) level, that were designed to collect data on various aspects of plant community and functional ecology. We also present the outcomes of two informal education initiatives that used citizen science to gather data on the distribution of plant diversity. We discuss these activities and results in light of their potential to engage the public into biodiversity conservation, and as educational and outreach tools.Formal education: UniversityDuring the COVID-19 pandemic (2021), Ecology practical classes of the Bologna Bachelor Degree in Biology (Faculty of Sciences of the University of Lisbon) had to be adapted to remote learning. Fortunately, during the States of Emergency imposed by the Portuguese Government, citizens were allowed to take brief walks. Taking advantage of citizen’s ability to briefly travel outdoors, we created three activities for students, as alternatives to those typically carried out in the classroom/campus, which we describe below.Activity 1—Analysis of the impact of disturbance on plant diversity in grasslandsThe objective of this activity was for students to explore the impact of disturbance and site attributes (such as soil type) on the diversity of the herbaceous plant community and its associated pollinators. This was undertaken in grasslands located near their homes, within walking distance (due to COVID lockdown movement restrictions). To achieve this goal, we developed a comprehensive sampling protocol that included methods for (i) selecting and characterizing sampling sites based on the level of human perturbation, (ii) soil characterization, (iii) sampling, identifying, and registering plants using the iNaturalist/Biodiversity4All platform and Flora-on web (Box 1), and (iv) pollinator sampling (Supplementary Data 1). To ensure accurate plant and pollinators identification, all observations were verified by professors responsible for each topic.First, each student chose one sampling site and teachers, using photographs, classified all sites regarding their perturbation level (low, medium, and high). Then, using the sampling protocol, students were invited to study different aspects of their sampling site, in loco or at their homes. Soil samples were analysed using simple methods and available household instruments (such as plastic cups, kitchen scale, and oven). Students were introduced to soil biodiversity as well as soil parameters (humidity, texture, structure, infiltration and draining) during the remote classes. Plants were sampled using a home-made 1 m2 quadrat. All species within were counted and identified to the lowest taxonomic level possible, using the mentioned apps and website. Before plant sampling, students were also asked to count and identify pollinators within their quadrats (broad taxonomic groups, bees, butterflies, flies, beetles) for 5 min, again using the apps to aid identification.Following field sampling, students were asked to calculate two taxonomic indices of plant communities. These included species richness, which measures the number of different species that occur in a sample, and the Simpson Diversity Index, which evaluates the probability that two individuals randomly selected from a sample will belong to the same species. Students also calculated functional diversity indices such as Functional Richness and Functional Dissimilarity, since functional diversity explores functional differences between species and how these differences reflect and affect the interactions with the environment and with other species29. Then, students assessed the relation between these indices and perturbation level. They analysed several functional traits of plants that are likely to respond to local perturbation (e.g., height, leaf size). Finally, they attempted to relate plant indices with the occurrence of pollinators.Overall, students sampled 147 grasslands that were affected by low (n = 17); medium (n = 86) and high (n = 40) levels of perturbation, scattered across mainland Portugal (Fig. 1a). In total, 3015 observations corresponding to 543 species of plant and 88 of insects (Fig. 1b) were registered in the iNaturalist/Biodiversity4All project Ecologia2_FCUL, created specifically to record all of the diversity data associated with this activity. Other registered taxa included six species of molluscs and 13 of arachnids, and other occasional soil macrofauna.Fig. 1: Analysis of the impact of disturbance on plant diversity in grasslands.a Location of grasslands sampled; b Banner and overview of main results of the project created in the platform iNaturalist/Biodiversity4All to register the sampled species; c Boxplots include data of the taxonomic diversity indices (plant species richness and Simpson Diversity Index) of sampled grasslands at three different perturbation levels: low, medium and high. Central lines represent median values, box limits indicate the upper and lower quartiles, whiskers correspond to 1.5 × the interquartile range above and below the upper and lower quartiles and points are the outliers. Boxplots with different letters indicate statistically significant differences among perturbation levels based on multiple pairwise comparisons.Full size imageThe results showed that the number of species (richness) decreased consistently with the level of perturbation. Simpson Diversity Index values increased, indicating low diversity values in highly perturbed herbaceous plant communities (Fig. 1c). Results revealed a trend towards an increase in the proportion of species with lower stature as perturbation increased. However, with no clear relationship with either biodiversity or perturbation. Finally, results indicated no clear relation of pollinator abundance or richness with plant richness and diversity, although field records relate a lower number of pollinators as wind intensity increased. In fact, pollinator sampling is extremely weather sensitive, which may have contributed to the lack of consistent relationships between pollinator diversity and perturbation.Box 1 Citizen science platforms and apps used for formal and informal educational activitiesiNaturalist (https://www.inaturalist.org/home): is a social network of naturalists, citizen scientists, and biologists that is based on mapping and sharing biodiversity observations. They describe themselves as “an online social network of people sharing biodiversity information in order to help each other learn about nature”. iNaturalist may be accessed via website or mobile app. Records are validated by the iNaturalist community. Observations reached approximately 110 million as of July 2022. This app allows the development of both open-access and registration-restricted projects. BioDiversity4All (https://www.biodiversity4all.org/) is a Portuguese biodiversity citizen science platform created by the Biodiversity for All Association. This platform was founded in 2010 and is currently linked to the “iNaturalist” network43. All the projects presented in this article were developed on the Biodiversity4All platform.Flora-on (https://flora-on.pt/): this portal contains occurrence data of vascular plants from the Portuguese flora collected by project collaborators (over 575,000 records as of July 2022). Flora-on was created by the Botanical Society of Portugal (SPBotânica), a Portuguese association devoted to the promotion and study of botany in Portugal. Botanists and naturalists provide most of the data, but occasional contributors are welcomed. Records are supervised by the portal editors, ensuring the dataset’s quality level. The portal includes stunning images of leaves, flowers, fruits, and other plant parts for 2299 of the 3300 taxa occurring in Portugal44. Additionally, the portal includes a powerful search engine that allows geographical, morphological, and taxonomical searches.LeafBite (https://zoegp.science/leafbyte): is a free, open-source iPhone app that measures total leaf area as well as consumed leaf area when herbivory is present45.Leaf-IT is a free and simple Android app created for scientific purposes. It was designed to measure leaf area under challenging field conditions. It has simple features for area calculation and data output, and can be used for ecological research and education46.Activity 2—Leaf trait assessment of shrub and tree speciesStudents were asked to assess three leaf traits Specific leaf area (SLA), Specific leaf mass (LMA), and Leaf Water Content (LWC) of two or three shrub or tree species. Each species should ideally fall into one of three functional groups known for their water adaptations, namely Hydrophytes, Mesophytes and Xerophytes. Students were challenged to choose charismatic Mediterranean species that grew nearby, such as Olea europaea, Nerium oleander or Phillyrea angustifolia. Alternatively, they could take the “Quercus challenge”, which involved ranking the Portuguese oak species based on their drought tolerance. A detailed protocol was developed to assist students for this purpose (Supplementary Data 2). In this protocol was demonstrated how to calculate the leaf area using the LeafBite and Leaf-IT apps (Box 1).The students calculated the SLA, LMA, and LWC of a total of 104 species (Supplementary Data 3) belonging to the main functional groups under study. Regarding the “Quercus challenge”, they were able to classify the six most representative oak species in Portugal and confirm the relationship among these indices and their tolerance to drought (Fig. 2).Fig. 2: Leaf trait assessment of shrub and tree species: Quercus challenge.Classification of Portuguese oak species regarding their drought tolerance (higher tolerance, left-up, lower tolerance right-down).Full size imageOne of the students, accomplished to present his own learning experience related to these activities at the XXIII Conference of the Environmental Research Network of Portuguese-speaking Nations – REALP, under the title “Plant Ecology during Confinement – A Digital Approach”.Activity 3—Evaluating the impact on the biodiversity of lawn management at the University of Lisbon campusAlthough, after the lockdown, practical classes returned to the laboratories and the field in 2021/22, we continued to use the iNaturalist/Biodiversity4All platform and the Flora-on website for biodiversity registering and identification, because of the success of the activities, as evidenced by the positive comments we received from students.The goal of this activity was to study the impact of lawn management on plant diversity and pollination on the University of Lisbon campus. To accomplish this, the students described the herbaceous communities and pollinators on four lawns (named C8, RL, RR, and TT) that had different management practices (mowing and irrigation). A comprehensive document with sampling guidelines was developed (Supplementary Data 4).The project Ecologia 2 Relvados 2022 registered 100 plant and 17 pollinator species (Fig. 3a). Given that the sampling took place during a cold and rainy week, which limited pollinator activity, the low number of pollinators registered was expectable (Lawson and Rands 2019). Following these analyses, the TT lawn (Fig. 3b), which had low levels of mowing and no watering, showed a significantly higher value of diversity, indicating it had the best management strategy for these systems (Fig. 3c), if the goal is to increase biodiversity.Fig. 3: Evaluating the impact on the biodiversity of lawn management at the University of Lisbon campus.a Banner and overview of main results of the project Ecologia 2 Relvados created in the platform iNaturalist/BioDiversity4All to register the sampled species; b Location of the lawns sampled in the Campus of the University of Lisbon; c Boxplots include data of the taxonomic diversity indices (plant species richness and Simpson Diversity Index) of sampled grasslands. Central lines represent median values, box limits indicate the upper and lower quartiles, whiskers correspond to 1.5 × the interquartile range above and below the upper and lower quartiles and points are the outliers. Boxplots with different letters indicate statistically significant differences among lawns based on multiple pairwise comparisons.Full size imageInformal education: BioBlitzesIntense biological surveys known as “BioBlitz” are carried out to record all organisms found in certain locations, such as cities, protected areas, or even entire countries. They are being used all over the world to collect and share georeferenced biodiversity data30. We developed two Plant Bioblitzes based on the BioDiversity4All/iNaturalist and Flora-on platforms. Social media, such as Facebook, Instagram, and Twitter, were used to promote these events and engage citizens (Fig. 4). The BioBlitzes were developed by SPBotânica in collaboration with BioDiversity4All.Fig. 4: Bioblitz I & II – Flora of Portugal.Posters created for the promotion of the two Flora of Portugal Bioblitzes.Full size imageBioblitz I & II – Flora of PortugalThe celebration of Fascination of Plants Day (18th of May) served as the backdrop for the organization of two-weekend Bioblitzes: Bioblitz Flora of Portugal I and Bioblitz Flora of Portugal II.In 2021, the Bioblitz was solely focused on project members, which meant that only those who had voluntarily joined the initiative could participate. In total, the 119 project members registered 4234 observations of 890 plant species. In contrast, the 2022 Bioblitz was an open project (no registration required). In total, the 323 observers made 6547 records of 1198 species. To evaluate the impact of the Bioblitz events, we compared the data registered in BioDiverstiy4All during the weekends of both events (2021 and 2022) with (i) the data registered in the platform during the equivalent weekends of 2019 and 2020 and (ii) also during the weekends before both Bioblitzes. The number of species, observations, and observers increased significantly from 2019 to 2020, 2021, and 2022, but, when comparing values from 2020 with 2021 and 2022, this rise was only verified during the Bioblitz weekends, proving the importance of Bioblitzes in this increase (Fig. 5).Fig. 5: Number of observations, species and observers registered on the BioDiversity4All/iNaturalist platform over equivalent weekends in 2019, 2020, 2021, and 2022.Numbers for 2021 and 2022 correspond to the weekends in which Bioblitzes I & II – Flora of Portugal were conducted, as well as previous ones.Full size image More

Bell, A. M. Behavioural differences between individuals and two populations of stickleback (Gasterosteus aculeatus). J. Evol. Biol. 18, 464–473 (2005).Article 
CAS 
PubMed 

Google Scholar 
Dochtermann, N. A. & Jenkins, S. H. Behavioural syndromes in Merriam’s kangaroo rats (Dipodomys merriami): A test of competing hypotheses. Proc. R. Soc. Lond. B 274, 2343–2349 (2007).
Google Scholar 
Tremmel, M. & Müller, C. Insect personality depends on environmental conditions. Behav. Ecol. 24, 386–392 (2013).Article 

Google Scholar 
Zidar, J. et al. A comparison of animal personality and coping styles in the red junglefowl. Anim. Behav. 130, 209–220 (2017).Article 

Google Scholar 
Sih, A., Bell, A. & Johnson, J. C. Behavioral syndromes: An ecological and evolutionary overview. Trends Ecol. Evol. 19, 372–378 (2004).Article 
PubMed 

Google Scholar 
Réale, D. & Dingemanse, N. J. Personality and individual social specialization. In Social behaviour: Genes, ecology and evolution (eds Székely, T. et al.) 417–441 (Cambridge University Press, 2010).Chapter 

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

Google Scholar 
Dingemanse, N. J., Kazem, A. J. N., Reale, D. & Wright, J. Behavioural reaction norms: Animal personality meets individual plasticity. Trends Ecol. Evol. 25, 81–89 (2010).Article 
PubMed 

Google Scholar 
Wolf, M., van Doorn, G. S. & Weissing, F. J. Evolutionary emergence of responsive and unresponsive personalities. Proc. Natl. Acad. Sci. USA 105, 15825–15830 (2008).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ólafsdóttir, G. Á. & Magellan, K. Interactions between boldness, foraging performance and behavioural plasticity across social contexts. Behav. Ecol. Sociobiol. 70, 1879–1889 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Mathot, K. J., Wright, J., Kempenaers, B. & Dingemanse, N. J. Adaptive strategies for managing uncertainty may explain personality-related differences in behavioural plasticity. Oikos 121(7), 1009–1020 (2012).Article 

Google Scholar 
Réale, D., Reader, S. M., Sol, D., McDougall, P. T. & Dingemanse, N. J. Integrating animal temperament within ecology and evolution. Biol. Rev. 82, 291–318 (2007).Article 
PubMed 

Google Scholar 
Coppens, C. M., de Boer, S. F. & Koolhaas, J. M. Coping styles and behavioural flexibility: Towards underlying mechanisms. Philos. Trans. R. Soc. B Biol. Sci. 365, 4021–4028 (2010).Article 

Google Scholar 
Benus, R. F., Daas, S. D., Koolhaas, J. M. & van Oortmerssen, G. A. Routine formation and flexibility in social and non-social behaviour of aggressive and non-aggressive male mice. Behaviour 112, 176–193 (1990).Article 

Google Scholar 
Dall, S. R., Houston, A. I. & McNamara, J. M. The behavioural ecology of personality: Consistent individual differences from an adaptive perspective. Ecol. Lett. 7, 734–739 (2004).Article 

Google Scholar 
Mitchell, D. J. & Biro, P. A. Is behavioural plasticity consistent across different environmental gradients and through time?. Proc. R. Soc. B. 284(1860), 20170893 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Stamps, J. A. Individual differences in behavioural plasticities. Biol. Rev. 91, 534–567 (2016).Article 
PubMed 

Google Scholar 
Stamps, J. A. & Biro, P. A. Personality and individual differences in plasticity. Curr. Opin. Behav. Sci. 12, 18–23 (2016).Article 

Google Scholar 
Dingemanse, N. J., Both, C., Drent, P. J. & Tinbergen, J. M. Fitness consequences of avian personalities in a fluctuating environment. Proc. R. Soc. Lond. B 271, 847 (2004).Article 

Google Scholar 
Smith, B. R. & Blumstein, D. T. Fitness consequences of personality: a meta-analysis. Behav. Ecol. 19, 448–455 (2008).Article 

Google Scholar 
Dingemanse, N. J. & Réale, D. Natural selection and animal personality. Behaviour 142, 1159–1184 (2005).Article 

Google Scholar 
Duque-Wilckens, N., Trainor, B. C. & Marler, C. A. Aggression and territoriality. In Encyclopedia of animal behavior (ed. Choe, J. C.) 539–546 (Elsevier, 2019).Chapter 

Google Scholar 
AmphibiaWeb. AmphibiaWeb: Information on amphibian biology and conservation. Available at https://amphibiaweb.org (2022).Ringler, M. et al. Acoustic ranging in poison frogs—It is not about signal amplitude alone. Behav. Ecol. Sociobiol. 71, 114 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Ringler, M., Ursprung, E. & Hödl, W. Site fidelity and patterns of short- and long-term movement in the brilliant-thighed poison frog Allobates femoralis (Aromobatidae). Behav. Ecol. Sociobiol. 63, 1281–1293 (2009).Article 

Google Scholar 
Ringler, M., Ringler, E., Magaña Mendoza, D. & Hödl, W. Intrusion experiments to measure territory size: Development of the method, tests through simulations, and application in the frog Allobates femoralis. PLoS ONE 6, e25844 (2011).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ringler, E., Ringler, M., Jehle, R. & Hödl, W. The female perspective of mating in A. femoralis, a territorial frog with paternal care—A spatial and genetic analysis. PLoS ONE 7, e40237 (2012).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ursprung, E., Ringler, M., Jehle, R. & Hödl, W. Strong male/male competition allows for nonchoosy females: High levels of polygynandry in a territorial frog with paternal care. Mol. Ecol. 20, 1759–1771 (2011).Article 
PubMed 

Google Scholar 
Pröhl, H. Territorial behavior in dendrobatid frogs. J Herpetol 39, 354–365 (2005).Article 

Google Scholar 
Peignier, M. et al. Exploring links between personality traits and their social and non-social environments in wild poison frogs. Behav. Ecol. Sociobiol. 76, 93 (2022).Article 

Google Scholar 
Chaloupka, S. et al. Repeatable territorial aggression in a Neotropical poison frog. Front. Ecol. Evol. 10, 398 (2022).Article 

Google Scholar 
Amézquita Torres, A. et al. Masking interference and the evolution of the acoustic communication system in the Amazonian dendrobatid frog Allobates femoralis. Evolution 60, 1874–1887 (2006).
Google Scholar 
Rodríguez López, C., Amézquita Torres, A., Ringler, M., Pašukonis, A. & Hödl, W. Calling amplitude flexibility and acoustic spacing in the territorial frog Allobates femoralis. Behav. Ecol. Sociobiol. 74, 1–10 (2020).
Google Scholar 
Asab. Guidelines for the treatment of animals in behavioural research and teaching. Anim. Behav. 159, 1–11 (2020).
Google Scholar 
Du Percie Sert, N. et al. The ARRIVE guidelines 2.0: Updated guidelines for reporting animal research. PLoS Biol 18, e3000410 (2020).Article 

Google Scholar 
Ringler, E., Mangione, R. & Ringler, M. Where have all the tadpoles gone? Individual genetic tracking of amphibian larvae until adulthood. Mol. Ecol. Resour. 15, 737–746 (2015).Article 
PubMed 

Google Scholar 
Ringler, M. et al. High-resolution forest mapping for behavioural studies in the Nature Reserve ‘Les Nouragues’, French Guiana. J. Maps 12, 26–32 (2016).Article 

Google Scholar 
Keith, D. A. et al. A function-based typology for Earth’s ecosystems. Nature 610, 513–518 (2022).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Kaefer, I. L., Montanarin, A., da Costa, R. S. & Lima, P. A. Temporal patterns of reproductive activity and site attachment of the brilliant-thighed frog Allobates femoralis from central Amazonia. J. Herpetol. 46, 549–554 (2012).Article 

Google Scholar 
Rasband, W. S. ImageJ (U. S. National Institutes of Health, 1997–2021).Bolger, D. T., Morrison, T. A., Vance, B., Lee, D. & Farid, H. A computer-assisted system for photographic mark–recapture analysis. Methods Ecol. Evol. 3, 813–822 (2012).Article 

Google Scholar 
Narins, P. M., Hödl, W. & Grabul, D. S. Bimodal signal requisite for agonistic behavior in a dart-poison frog, Epipedobates femoralis. Proc. Natl. Acad. Sci. USA 100, 577–580 (2003).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Gasser, H., Amézquita Torres, A. & Hödl, W. Who is calling? Intraspecific call variation in the aromobatid frog Allobates femoralis. Ethology 115, 596–607 (2009).Article 

Google Scholar 
Hödl, W. Dendrobates femoralis (Dendrobatidae): a handy fellow for frog bioacoustics in Proceedings of the 4th Ordinary General meeting of the Societas Europaea Herpetologica, (ed.van Gelder, J. J., Strijbosch, H. & Bergers, P.) (1987).Ursprung, E., Ringler, M. & Hödl, W. Phonotactic approach pattern in the neotropical frog Allobates femoralis: A spatial and temporal analysis. Behaviour 146, 153–170 (2009).Article 

Google Scholar 
Sonnleitner, R., Ringler, M., Loretto, M.-C. & Ringler, E. Experience shapes accuracy in territorial decision-making in a poison frog. Biol. Lett. 16, 20200094 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Hödl, W. Phyllobates femoralis (Dendrobatidae): Rufverhalten und akustische Orientierung der Männchen (Freilandaufnahmen) in Bundesstaatliche Hauptstelle für Wissenschaftliche Kinematographie (1983).Tumulty, J. P. et al. Brilliant-thighed poison frogs do not use acoustic identity information to treat territorial neighbours as dear enemies. Anim. Behav. 141, 203–220 (2018).Article 

Google Scholar 
Fernandes, I. Y. et al. Unlinking the speciation steps: Geographical factors drive changes in sexual signals of an Amazonian Nurse-Frog through body size variation. Evol. Biol. 48, 81–93 (2021).Article 

Google Scholar 
Garcia, M. J. et al. Dueling frogs: do male green tree frogs (Hyla cinerea) eavesdrop on and assess nearby calling competitors?. Behav. Ecol. Sociobiol. 73(2), 1041 (2019).Article 

Google Scholar 
Gingras, B., Böckle, M., Herbst, C. T. & Fitch, W. T. Call acoustics reflect body size across four clades of anurans. J Zool 289(2), 143–150 (2013).Article 

Google Scholar 
Stoffel, M. A., Nakagawa, S. & Schielzeth, H. rptR: Repeatability estimation and variance decomposition by generalized linear mixed-effects models. Methods Ecol. Evol. 8, 1639–1644 (2017).Article 

Google Scholar 
Fox, J. et al. Package ‘sem’: Structural Equation Models. https://CRAN.R-project.org/package=sem (2022).Burnham, K. P. & Anderson, D. R. Model selection and multimodel inference: A practical information-theoretic approach 2nd edn. (Springer, 2002).MATH 

Google Scholar 
Hertel, A. G., Niemelä, P. T., Dingemanse, N. J. & Mueller, T. A guide for studying among-individual behavioral variation from movement data in the wild. Mov. Ecol. 8, 30 (2020).Article 
PubMed 
PubMed Central 

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

Google Scholar 
Bürkner, P.-C. brms: An R package for Bayesian multilevel models using stan. J. Stat. Softw. 80, 1–28 (2017).Article 

Google Scholar 
Gelman, A., Hwang, J. & Vehtari, A. Understanding predictive information criteria for Bayesian models. Stat. Comput. 24, 997–1016 (2014).Article 
MathSciNet 
MATH 

Google Scholar 
Whalen, A. & Hoppitt, W. J. E. Bayesian model selection with network based diffusion analysis. Front. Psychol. 7, 409 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
R Core Team. R: A language and environment for statistical computing (R Foundation for Statistical Computing, 2019).
Google Scholar 
Ryan, M. J., Bartholomew, G. A. & Rand, A. S. Energetics of reproduction in a neotropical frog, Physalaemus pustulosus. Ecology 64, 1456–1462 (1983).Article 

Google Scholar 
Taigen, T. L. & Wells, K. D. Energetics of vocalization by an anuran amphibian (Hyla versicolor). J. Comp. Physiol. 155, 163–170 (1985).Article 

Google Scholar 
Pough, F. H. & Taigen, T. L. Metabolic correlates of the foraging and social behaviour of dart-poison frogs. Anim. Behav. 39, 145–155 (1990).Article 

Google Scholar 
Bell, A. M., Hankison, S. J. & Laskowski, K. L. The repeatability of behaviour: A meta-analysis. Anim. Behav. 77, 771–783 (2009).Article 
PubMed 
PubMed Central 

Google Scholar 
Kelleher, S. R., Silla, A. J. & Byrne, P. G. Animal personality and behavioral syndromes in amphibians: A review of the evidence, experimental approaches, and implications for conservation. Behav. Ecol. Sociobiol. 72, 10539 (2018).Article 

Google Scholar 
Moser-Purdy, C., MacDougall-Shackleton, E. A. & Mennill, D. J. Enemies are not always dear: Male song sparrows adjust dear enemy effect expression in response to female fertility. Anim. Behav. 126, 17–22 (2017).Article 

Google Scholar  More

Data sourcesMODIS on the Aqua satellite provides a global coverage within 1–2 days. All images acquired by this satellite mission from January 2003 to December 2020 were used in our study to detect global coastal phytoplankton blooms, with a total of 0.76 million images. MODIS Level-1A images were downloaded from the Ocean Biology Distributed Active Archive Center (OB.DAAC) at NASA Goddard Space Flight Center (GSFC), and were subsequently processed with SeaDAS software (version 7.5) to obtain Rayleigh-corrected reflectance (Rrc (dimensionless), which was converted using the rhos (in sr−1) product (rhos × π) from SeaDAS)41, remote sensing reflectance (Rrs (sr−1)) and quality control flags (l2_flags). If a pixel was flagged by any of the following, it was then removed from phytoplankton bloom detection: straylight, cloud, land, high sunglint, high solar zenith angle and high sensor zenith angle (https://oceancolor.gsfc.nasa.gov/atbd/ocl2flags/). MODIS level-3 product for aerosol optical thicknesses (AOT) at 869 nm was also obtained from OB.DAAC NASA GSFC (version R2018.0), which was used to examine the impacts of aerosols on bloom trends.We examined the algal blooms in the EEZs of 153 ocean-bordering countries (excluding the EEZs in the Caspian Sea or around the Antarctic), 126 of which were found with at least one bloom in the past two decades. The EEZ dataset is available at https://www.marineregions.org/download_file.php?name=World_EEZ_v11_20191118.zip. The EEZs are up to 200 nautical miles (or 370 km) away from coastlines, which include all continental shelf areas and offer the majority of marine resources available for human use. Regional statistics of algal blooms were also performed for LMEs. LMEs encompass global coastal oceans and outer edges of coastal currents areas, which are defined by various distinct features of the oceans, including hydrology, productivity, bathymetry and trophically dependent populations42. Of the 66 LMEs identified globally, we excluded the Arctic and Antarctic regions and examined 54 LMEs. The boundaries of LMEs were obtained from https://www.sciencebase.gov/catalog/item/55c77722e4b08400b1fd8244.We used HAEDAT to validate our satellite-detected phytoplankton blooms in terms of presence or absence. The HAEDAT dataset (http://haedat.iode.org) is a collection of records of HAB events, maintained under the UNESCO Intergovernmental Oceanographic Commission and with data archives since 1985. For each HAB event, the HAEDAT records its bloom period (ranging from days to months) and geolocation. We merged duplicate entries when both the recorded locations and times of the HAEDAT events were very similar to one another, and a total number of 2,609 HAEDAT events were ultimately selected between 2003 and 2020.We used the ¼° resolution National Oceanic and Atmospheric Administration Optimum Interpolated SST (v. 2.1) data to examine the potential simulating effects of warming on the global phytoplankton trends. We also estimated the SST gradients following the method of Martínez-Moreno33. As detailed in ref. 33, the SST gradient can be used as a proxy for the magnitude of oceanic mesoscale currents (EKE). We used the SST gradient to explore the effects of ocean circulation dynamics on algal blooms.Fertilizer uses and aquaculture production for different countries was used to examine the potential effects of nutrient enrichment from humans on global phytoplankton bloom trends. Annual data between 2003 and 2019 on synthetic fertilizer use, including nitrogen and phosphorus, are available from https://ourworldindata.org/fertilizers. Annual aquaculture production includes cultivated fish and crustaceans in marine and inland waters, and sea tanks, and the data between 2003 and 2018 are available from https://ourworldindata.org/grapher/aquaculture-farmed-fish-production.The MEI, which combines various oceanic and atmospheric variables36, was used to examine the connections between El Niño–Southern Oscillation activities and marine phytoplankton blooms. The dataset is available from https://psl.noaa.gov/enso/mei/.Development of an automated bloom detection methodA recent study by the UNESCO Intergovernmental Oceanographic Commission revealed that globally reported HAB events have increased6. However, such an overall increasing trend was found to be highly correlated with recently intensified sampling efforts6. Once this potential bias was accounted for by examining the ratio between HAB events to the number of samplings5, there was no significant global trend in HAB incidence, though there were increases in certain regions. With synoptic, frequent, and large-scale observations, satellite remote sensing has been extensively used to monitor algal blooms in oceanic environments17,18,19. For example, chlorophyll a (Chla) concentrations, a proxy for phytoplankton biomass, has been provided as a standard product by NASA since the proof-of-concept Coastal Zone Color Scanner (1978–1986) era43,44. The current default algorithm used to retrieve Chla products is based on the high absorption of Chla at the blue band45,46, which often shows high accuracy in the clear open oceans but high uncertainties in coastal waters. This is because, in productive and dynamic coastal oceans, the absorption of Chla in the blue band can be obscured by the presence of suspended sediments and/or coloured dissolved organic matter (CDOM)47. To address this problem, various regionalized Chla algorithms have been developed48. Unfortunately, the concentrations of the water constituents (CDOM, sediment and Chla) can vary substantially across different coastal oceans. As a result, a universal Chla algorithm that can accurately estimate Chla concentrations in global coastal oceans is not currently available.Alternatively, many spectral indices have been developed to identify phytoplankton blooms instead of quantifying their bloom biomass, including the normalized fluorescence line height21 (nFLH), red tide index49 (RI), algal bloom index47 (ABI), red–blue difference (RBD)50, Karenia brevis bloom index50 (KBBI) and red tide detection index51 (RDI). In practice, the most important task for these index-based algorithms is to determine their optimal thresholds for bloom classification. However, such optimal thresholds can be regional-or image-specific20, due to the complexity of optical features in coastal waters and/or the contamination of unfavourable observational conditions (such as thick aerosols, thin clouds, and so on), making it difficult to apply spectral-index-based algorithms at a global scale.To circumvent the difficulty in determining unified thresholds for various spectral indices across global coastal oceans, an approach from a recent study to classify algal blooms in freshwater lakes52 was adopted and modified here. In that study, the remotely sensed reflectance data in three visible bands (red, green and blue) were converted into two-dimensional colour space created by the Commission Internationale del’éclairage (CIE), in which the position on the CIE chromaticity diagram represented the colour perceived by human eyes (Extended Data Fig. 1a). As the algal blooms in freshwater lakes were manifested as greenish colours, the reflectance of bloom-containing pixels was expected to be distributed in the green gamut of the CIE chromaticity diagram; the stronger the bloom, the closer the distance to the upper border of the diagram (the greener the water).Here, the colour of phytoplankton blooms in the coastal oceans can be greenish, yellowish, brownish, or even reddish53, owing to the compositions of bloom species (diatoms or dinoflagellates) and the concentrations of different water constituents. Furthermore, the Chla concentrations of the coastal blooms are typically lower than those in inland waters, thus demanding more accurate classification algorithms. Thus, the algorithm proposed by Hou et al.52 was modified when using the CIE chromaticity space for bloom detection in marine environments. Specifically, we used the following coordinate conversion formulas to obtain the xy coordinate values in the CIE colour space:$$begin{array}{c}x=X/(X+Y+Z)\ y=Y/(X+Y+Z)\ X=2.7689R+1.7517G+1.1302B\ Y=1.0000R+4.5907G+0.0601B\ Z=0.0000R+0.0565G+5.5943Bend{array}$$
(1)
where R, G and B represent the Rrc at 748 nm, 678 nm (fluorescence band) and 667 nm in the MODIS Aqua data, respectively. By contrast, the R, G and B channels used in Hou et al.52 were the red, green and blue bands. We used the fluorescence band for the G channel because, for a given region, the 678 nm signal increases monotonically with the Chla concentration for blooms of moderate intensity21, which is similar to the response of greenness to freshwater algal blooms. Thus, the converted y value in the CIE coordinate system represents the strength of the fluorescence. In practice, for pixels with phytoplankton blooms, the converted colours in the chromaticity diagram will be located within the green, yellow or orange–red gamut (see Extended Data Fig. 1a); the stronger the fluorescence signal is, the closer the distance to the upper border of the CIE diagram (larger y value). By contrast, for bloom-free pixels without a fluorescence signal, their converted xy coordinates will be located in the blue or purple gamut. Therefore, we can determine a lower boundary in the CIE two-dimensional coordinate system to separate bloom and non-bloom pixels, similar to the method proposed by Hou et al.52.We selected 53,820 bloom-containing pixels from the MODIS Rrc data as training samples to determine the boundary of the CIE colour space. These sample points were selected from nearshore waters worldwide where frequent phytoplankton blooms have been reported (Extended Data Fig. 2); the algal species included various species of dinoflagellates and diatoms20. A total of 80 images was used, which were acquired from different seasons and across various bloom magnitudes, to ensure that the samples used could almost exhaustively represent the different bloom conditions in the coastal oceans.We combined the MODIS FLHRrc (fluorescence line height based on Rrc) and enhanced red–green–blue composite (ERGB) to delineate bloom pixels manually. The FLHRrc image was calculated as:$$begin{array}{c}{{rm{FLH}}}_{{rm{Rrc}}}={R}_{{rm{rc}}678}times {F}_{678}-[{R}_{{rm{rc}}667}times {F}_{667}+({R}_{{rm{rc}}748}times {F}_{748}\ ,,-,{R}_{{rm{rc}}667}times {F}_{667})times (678-667)/(748-667)]end{array}$$
(2)
where Rrc667, Rrc678 and Rrc748 are the Rrc at 667, 678 and 748 nm, respectively, and F667, F678 and F748 are the corresponding extraterrestrial solar irradiance. ERGB composite images were generated using Rrc of three bands at 555 (R), 488 (G) and 443 nm (B). Although phytoplankton-rich and sediment-rich waters have high FLHRrc values, they appear as darkish and bright features in the ERGB images (Extended Data Fig. 3), respectively21. In fact, visual examination with fluorescence signals and ERGB has been widely accepted as a practical way to delineate coastal algal blooms on a limited number of images21,54,55. Note that the FLHRrc here was slightly different from the NASA standard nFLH product56, as the latter is generated using Rrs (corrected for both Rayleigh and aerosol scattering) instead of Rrc (with residual effects of aerosols). However, when using the NASA standard algorithm to further perform aerosol scattering correction over Rrc, 20.7% of our selected bloom-containing pixels failed to obtain valid Rrs (without retrievals or flagged as low quality), especially for those with strong blooms (see examples in Extended Data Fig. 4). Likewise, we also found various nearshore regions with invalid Rrs retrievals. By contrast, Rrc had valid data for all selected samples and showed more coverage in nearshore coastal waters. The differences between Rrs and Rrc were because the assumptions for the standard atmospheric correction algorithm do not hold for bloom pixels or nearshore waters with complex optical properties57. In fact, Rrc has been used as an alternative to Rrs in various applications in complex waters58,59.We converted the Rrc data of 53,820 selected sample pixels into the xy coordinates in the CIE colour space (Extended Data Fig. 1a). As expected, these samples of bloom-containing pixels were located in the upper half of the chromaticity diagram (the green, yellow and orange–red gamut) (Extended Data Fig. 1a). We determined the lower boundary of these sample points in the chromaticity diagram, which represents the lightest colour and thus the weakest phytoplankton blooms; any point that falls above this boundary represents stronger blooms. The method to determine the boundary was similar to Hou et al.52: we first binned the sample points according to the x value in the chromaticity diagram and estimated the 1st percentile (Q1%) of the corresponding Y for each bin; then, we fit the Q1% using two-order polynomial regression. Sensitivity analysis with Q0.3% (the three-sigma value) resulted in minor changes ( 1/3 AND y  > y2), it is classified as a ‘bloom’ pixel.Depending on the local region and application purpose, the meaning of ‘phytoplankton bloom’ may differ. Here, for a global application, the pixelwise bloom classification is based on the relationship (represented using the CIE colour space) between Rrc in the 667-, 678- and 754-nm bands derived from visual interpretation of the 80 pairs of FLHRrc and ERGB imagery. Instead of a simple threshold, we used a lower boundary of the sample points in the chromaticity diagram to define a bloom. In simple words, a pixel is classified as a bloom if its fluorescence signal is detectable (the associated xy coordinate in the CIE colour space located above the lower boundary). Histogram of the nFLH values from the 53,820 training pixels demonstrated the minimum value of ~0.02 mW cm−2 μm−1 (Extended Data Fig. 1a), which is in line with the lower-bound signal of K. brevis blooms on the West Florida shelf21,47. Note that, such a minimum nFLH is determined from the global training pixels, and it does not necessarily represent a unified lower bound for phytoplankton blooms across the entire globe, especially considering that fluorescence efficiency may be a large variable across different regions. Different regions may have different lower bounds of nFLH to define a bloom, and such variability is represented by the predefined boundary in the CIE chromaticity diagram in our study. Correspondingly, although the accuracy of Chla retrievals may have large uncertainties in coastal waters, the histogram of the 53,820 training pixels shows a lower bound of ~1 mg m−3 (Extended Data Fig. 1a). Similarly to nFLH, such a lower bound may not be applicable to all coastal regions, as different regions may have different lower bounds of Chla for bloom definition.Although the MODIS cloud (generated by SeaDAS with Rrc869 0.12) and Index2 ( More

Chazdon, R. L. et al. Carbon sequestration potential of second-growth forest regeneration in the Latin American tropics. Sci. Adv. 2, e1501639 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
IKI. The Bonn Challenge. https://www.bonnchallenge.org/ (2022).UNCCD. Land Degradation Neutrality. https://www.unccd.int/land-and-life/land-degradation-neutrality/overview (2022).Bastin, J.-F. et al. The global tree restoration potential. Science 365, 76–79 (2019).Article 
CAS 
PubMed 

Google Scholar 
Brancalion, P. H. S. et al. Global restoration opportunities in tropical rainforest landscapes. Sci. Adv. 5, eaav3223 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Strassburg, B. B. N. et al. Global priority areas for ecosystem restoration. Nature 586, 724–729 (2020).Article 
CAS 
PubMed 

Google Scholar 
Erbaugh, J. T. et al. Global forest restoration and the importance of prioritizing local communities. Nat. Ecol. Evol. 4, 1472–1476 (2020).Article 
CAS 
PubMed 

Google Scholar 
Fleischman, F. et al. Restoration prioritization must be informed by marginalized people. Nature 607, E5–E6 (2022).Article 
CAS 
PubMed 

Google Scholar 
Chaturvedi, R. et al. Restoration Opportunities Atlas of India. www.india.restorationatlas.org/methodology (2022).McLain, R., Lawry, S., Guariguata, M. R. & Reed, J. Toward a tenure-responsive approach to forest landscape restoration: a proposed tenure diagnostic for assessing restoration opportunities. Land Use Policy 104, 103748 (2021).Article 

Google Scholar 
Binod, B., Bhattarcharjee, A. & Ishwar, N. M. Bonn Challenge and India: Progress on Restoration Efforts Across States and Landscapes (IUCN, 2018).Government of India. Aspirational Districts Phase 1 (vikaspedia, 2018).Government of India. Census of India. https://censusindia.gov.in/2011census/dchb/DCHB.html (2011).DeFries, R. et al. Land management can contribute to net zero. Science 376, 1163–1165 (2022).Article 
CAS 
PubMed 

Google Scholar 
Borah, B., Bhattacharya, A. & Ishwar, N. M. Bonn Challenge and India. Progress On Restoration Efforts Across States and Landscapes. https://www.bonnchallenge.org/pledges/india (2018).Gopalakrishna, T. et al. Existing land uses constrain climate change mitigation potential of forest restoration in India. Conserv. Lett. https://doi.org/10.1111/conl.12867 (2022).Dhyani, S. et al. Agroforestry to achieve global climate adaptation and mitigation targets: are South Asian countries sufficiently prepared? Forests 12, 303 (2021).Article 

Google Scholar 
Nerlekar, A. N. et al. Removal or utilization? Testing alternative approaches to the management of an invasive woody legume in an arid Indian grassland. Restor. Ecol. https://doi.org/10.1111/rec.13477 (2022).Coleman, E. A. et al. Limited effects of tree planting on forest canopy cover and rural livelihoods in Northern India. Nat Sustain 4, 997–1004 (2021).Article 

Google Scholar 
Ramprasad, V., Joglekar, A. & Fleischman, F. Plantations and pastoralists: afforestation activities make pastoralists in the Indian Himalaya vulnerable. Ecol. Soci. https://doi.org/10.5751/ES-11810-250401 (2020).DeFries, R. et al. Improved household living standards can restore dry tropical forests. Biotropica https://doi.org/10.1111/btp.12978 (2021).Lele, S., Khare, A. & Mokashi, S. Estimating and Mapping CFR Potential (ATREE, 2020).Agarwala, M. et al. Impact of biogas interventions on forest biomass and regeneration in southern India. Global Ecol. Conservation 11, 213–223 (2017).Article 

Google Scholar 
Menon, A. & Schmidt-Vogt, D. Effects of the COVID-19 pandemic on farmers and their responses: a study of three farming systems in Kerala. South India. Land 11, 144 (2022).
Google Scholar 
Fremout, T. et al. Diversity for Restoration (D4R): Guiding the selection of tree species and seed sources for climate‐resilient restoration of tropical forest landscapes. J. Appl. Ecol. 59, 664–679 (2022).Article 

Google Scholar 
Hughes, K. A. et al. Can restoration of the commons reduce rural vulnerability? A Quasi-experimental comparison of COVID-19 livelihood-based coping strategies among rural households in three Indian States. Int. J. Common. 16, 189 (2022).Article 

Google Scholar 
Madhusudan, M. D. & Vanak, A. Mapping the Distribution and Extent of India’s Semi-arid Open Natural Ecosystems. https://doi.org/10.1002/essoar.10507612.1 (2021).Vanak, A. T., Hiremath, A. J., Ganesh, T. & Rai, N. D. Filling in the (Forest) Blanks: the Past, Present and Future of India’s Savanna Grasslands (ATREE, 2017).Oxford Poverty & Human Development Initiative. Global Multidimensional Poverty Index 2018. The Most Detailed Picture to Date of the World’s Poorest People. https://ophi.org.uk/wp-content/uploads/G-MPI_2018_2ed_web.pdf (2018).Hijmans, R. J. et al. raster: Geographic Data Analysis and Modeling. https://rspatial.org/raster (2023).Bivand, R. et al. rgdal: Bindings for the ‘Geospatial’ Data Abstraction Library. https://cran.r-project.org/web/packages/rgdal/index.html (2023).QGIS Development Team. QGIS Geographic Information System (Open Source Geospatial Foundation Project, 2022). More

Wilkinson, D. At cross purposes. Nature 412, 485 (2001).Article 
ADS 
CAS 
PubMed 

Google Scholar 
de Bary, H. A. Über Symbiose [On Symbiosis]. Tageblatt für die Versammlung Dtsch. Naturforscher und Aerzte (in Cassel) [Daily J. Conf. Ger. Sci. Phys.] (in Ger. 51, 121–126 (1878).Lücking, R., Leavitt, S. D. & Hawksworth, D. L. Species in lichen-forming fungi: balancing between conceptual and practical considerations, and between phenotype and phylogenomics. Fungal Div.109, 99–154 (Springer, Netherlands, 2021).de Vries, J. & Archibald, J. M. Plant evolution: Landmarks on the path to terrestrial life. New Phytol. 217, 1428–1434 (2018).Article 
PubMed 

Google Scholar 
Ahmadjian, V. The Lichen Symbiosis (John Wiley & Sons, 1993).
Google Scholar 
Lücking, R., Hodkinson, B. P. & Leavitt, S. D. The 2016 classification of lichenized fungi in the Ascomycota and Basidiomycota-approaching one thousand genera. Bryologist 119, 361–416 (2016).Article 

Google Scholar 
Schneider, K., Resl, P. & Spribille, T. Escape from the cryptic species trap: lichen evolution on both sides of a cyanobacterial acquisition event. Mol. Ecol. 25, 3453–3468 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Wedin, M., Döring, H. & Gilenstam, G. Saprotrophy and lichenization as options for the same fungal species on different substrata: Environmental plasticity and fungal lifestyles in the Stictis-Conotrema complex. New Phytol. 164, 459–465 (2004).Article 

Google Scholar 
Muggia, L., Baloch, E., Stabentheiner, E., Grube, M. & Wedin, M. Photobiont association and genetic diversity of the optionally lichenized fungus Schizoxylon albescens. FEMS Microbiol. Ecol. 75, 255–272 (2011).Article 
CAS 
PubMed 

Google Scholar 
Sanders, W. B., Moe, R. L. & Ascaso, C. Ultrastructural study of the brown alga Petroderma maculiforme (Phaeophyceae) in the free-living state and in lichen symbiosis with the intertidal marine fungus Verrucaria tavaresiae (Ascomycotina). Eur. J. Phycol. 40, 353–361 (2005).Article 
CAS 

Google Scholar 
Vondrák, J. et al. From Cinderella to Princess. Preslia 94, 143–181 (2022).Article 

Google Scholar 
Hawksworth, D. L. The variety of fungal-algal symbioses, their evolutionary significance, and the nature of lichens. Bot. J. Linn. Soc. 96, 3–20 (1988).Article 

Google Scholar 
Larsson, K. H. & Ryvarden, L. Corticioid fungi of Europe 1. Acanthobasidium–Gyrodontium. Synop. Fungorum 43, 1–266 (2021).
Google Scholar 
Albertini, J. B., von Schweinitz, L. D. Conspectus fungorum in Lusatiae Superioris agro Niskiensi crescentium, e methodo Persooniana. (DE: Sumtibus Kummerianis, Lipsiae 1805) https://doi.org/10.5962/bhl.title.3601.Poelt, J. & Jülich, W. Über die Beziehungen zweier corticioider Basidiomyceten zu Algen. Österr. Bot. Zeitschrift 116, 400–410 (1969).Article 

Google Scholar 
Voytsekhovich, A., Ordynets, O. & Akimov, Y. Optionally lichenized fungi of Hyphodontia (Agaricomycetes, Schizoporaceae) and their photobiont composition. Aктyaльнi Пpoблeми Бoтaнiки Ta Eкoлoгiї. Maтepiaли Miжнapoднoї Кoнфepeнцiї Moлoдиx Учeниx 65 (2013).Voytsekhovich, A., Mikhailyuk, T., Akimov, Y., Ordynets, A., Gustavs, L. Optionally lichenized fungi of Hyphodontia (Agaricomycetes, Schizoporaceae). 8th Congress of the International Symbiosis Society, Lisbon, 12–18 July 2015. Lisbon, PT:, 217 (Conf. abstract) (2015).Gustavs L, Schiefelbein U, Darienko T, P. T. Symbioses of the green algal genera Coccomyxa and Elliptochloris (Trebouxiophyceae, Chlorophyta). in Algal and Cyanobacteria Symbioses (ed. Grube M, Seckbach J) 169–208 (2017).Darienko, T., Gustavs, L., Eggert, A., Wolf, W. & Pröschold, T. Evaluating the species boundaries of green microalgae (Coccomyxa, Trebouxiophyceae, Chlorophyta) using integrative taxonomy and DNA barcoding with further implications for the species identification in environmental samples. PLoS ONE 10, 1–31 (2015).Article 

Google Scholar 
Malavasi, V. et al. DNA-based taxonomy in ecologically versatile microalgae: A re-evaluation of the species concept within the coccoid green algal genus Coccomyxa (Trebouxiophyceae, Chlorophyta). PLoS ONE 11, e0151137 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Green, T. G. A., Nash, T. H. Lichen Biology. In Lichen Biology, Second Edition 152–181 (Cambridge University Press, Cambridge, 2008) https://doi.org/10.1017/CBO9780511790478.Lindgren, H. et al. Cophylogenetic patterns in algal symbionts correlate with repeated symbiont switches during diversification and geographic expansion of lichen-forming fungi in the genus Sticta (Ascomycota, Peltigeraceae). Mol. Phylogenet. Evol. 150, 106860 (2020).Article 
PubMed 

Google Scholar 
Kulichová, J., Škaloud, P. & Neustupa, J. Molecular diversity of green corticolous microalgae from two sub-mediterranean European localities. Eur. J. Phycol. 49, 345–355 (2014).Article 

Google Scholar 
Pröschold, T. & Darienko, T. The green puzzle Stichococcus (Trebouxiophyceae, Chlorophyta): New generic and species concept among this widely distributed genus. Phytotaxa 441, 113–142 (2020).Article 

Google Scholar 
Meier, F. A., Scherrer, S. & Honegger, R. Faecal pellets of lichenivorous mites contain viable cells of the lichen-forming ascomycete Xanthoria parietina and its green algal photobiont. Trebouxia arboricola. Biol. J. Linn. Soc. 76, 259–268 (2002).Article 

Google Scholar 
Bernicchia, A. & Gorjón, S. P. Corticiaceae s.l. 1008 (2010), ISBN: 9788890105791.Parmasto, E. Descriptiones taxorum novorum. Combinationes novae. Proc. Acad. Sci. Est. SSR. Biol. 16, 377–394 (1967).
Google Scholar 
Hjortstam, K., Larsson, K., Ryvarden, L. & Eriksson, J. The Corticiaceae of North Europe. (Oslo: Fungiflora, 1988).Jaag, O. Coccomyxa schmidle Monographie einer algengattung. Beitr. Kryptogamenflora Schweiz 8, 1–132 (1933).
Google Scholar 
Oberwinkler, F. Die gattungen der Basidiolichenen. Vorträge aus dem Gesamtgebiet der Botanik. Herausgegeb. v. d. Deutsch. bot. Ges. Neue Folge 4, 139–169 (1970).
Google Scholar 
Poelt, J. Basidienflechten, eine in den Alpen lange übersehene Pflanzengruppe. Jahrb. Vereins Schutze Alpenpfl. Tiere 40, 81–92 (1975).
Google Scholar 
Eriksson, J., Hjortstam, K. The Corticiaceae of North Europe. Vol. 6. (Grønlands Eskefabrikk, 1981).Oberwinkler, F. Basidiolichens. In Fungal Association 211–225 (Springer, Berlin Heidelberg, Berlin, 2001). https://doi.org/10.1007/978-3-662-07334-6_12.Chapter 

Google Scholar 
Jülich, W. A new lichenized Athelia from Florida. Persoonia 10, 149–151 (1978).
Google Scholar 
Zavada, M. S. & Simoes, P. The possible demi-lichenization of the basidiocarps of Trametes Versicolor (L.:Fries) pilat (polyporaceae). Northeast. Nat. 8, 101–112 (2001).
Google Scholar 
Neustroeva, N., Mukhin, V., Novakovskaya, I. & Patova, E. Biodiversity of symbiotic algae of wood decay Basidimycetes in the Central Urals. III Russ. Natl. Conf. “Information Technol. Biodivers. Res. 1, 83–92 (2020).
Google Scholar 
Zavada, M. S., DiMichele, L. & Toth, C. R. The possible demi-lichenization of Trametes versicolor (L.: Fries) Pilát (Polyporaceae): The transfer of fixed 14CO2 from epiphytic algae to T. versicolor. Northeast. Nat. 11, 33–40 (2004).Article 

Google Scholar 
Mukhin, V. A., Patova, E. N., Kiseleva, I. S., Neustroeva, N. V. & Novakovskaya, I. V. Mycetobiont symbiotic algae of wood-decomposing fungi. Russ. J. Ecol. 47, 133–137 (2016).Article 
CAS 

Google Scholar 
Sanders, W. B. & Masumoto, H. Lichen algae: The photosynthetic partners in lichen symbioses. Lichenologist 53, 347–393 (2021).Article 

Google Scholar 
Krause, G. & Weis, E. Chlorophyll fluorescence and photosynthesis: the basics. Annu. Rev. Plant Biol. 42(1), 313–349 (1991).Article 
CAS 

Google Scholar 
Lüttge, U. & Büdel, B. Resurrection kinetics of photosynthesis in desiccation-tolerant terrestrial green algae (Chlorophyta) on tree bark. Plant Biol. 12, 437–444 (2010).Article 
PubMed 

Google Scholar 
Lange, O. L. Moisture content and CO2 exchange of lichens: I. Influence of temperature on moisture-dependent net photosynthesis and dark respiration in Ramalina maciformis. Oecologia 45, 82–87 (1980).Article 
ADS 
PubMed 

Google Scholar 
Palmqvist, K. & Sundberg, B. Light use efficiency of dry matter gain in five macrolichens: Relative impact of microclimate conditions and species-specific traits. Plant Cell Environ. 23, 1–14 (2000).Article 

Google Scholar 
Vondrak, J. & Kubásek, J. Algal stacks and fungal stacks as adaptations to high light in lichens. Lichenol. 45(1), 115 (2013).Article 

Google Scholar 
Smith, N. G. & Dukes, J. S. Plant respiration and photosynthesis in global-scale models: Incorporating acclimation to temperature and CO2. Glob. Chang. Biol. 19, 45–63 (2013).Article 
ADS 
PubMed 

Google Scholar 
Medeiros, P. M. & Simoneit, B. R. T. Analysis of sugars in environmental samples by gas chromatography-mass spectrometry. J. Chromatogr. A 1141, 271–278 (2007).Article 
CAS 
PubMed 

Google Scholar 
Honegger, R. Functional aspects of the lichen symbiosis. Annu. Rev. Plant Physiol. Plant Mol. Biol. 42, 553–578 (1991).Article 
CAS 

Google Scholar 
Honegger, R. The lichen symbiosis—What is so spectacular about it?. Lichenologist 30, 193–212 (1998).Article 

Google Scholar 
Kirk, P. M. et al. (eds) Dictionary of the Fungi 10th edn. (CABI, Netherlands, 2008).
Google Scholar 
Ahmadjian, V. The lichen alga Trebouxia: Does it occur free-living?. Plant Syst. Evol. 158, 243–247 (1988).Article 

Google Scholar 
Sanders, W. B. Complete life cycle of the lichen fungus Calopadia puiggarii (Pilocarpaceae, Ascomycetes) documented in situ: Propagule dispersal, establishment of symbiosis, thallus development, and formation of sexual and asexual reproductive structures. Am. J. Bot. 101, 1836–1848 (2014).Article 
PubMed 

Google Scholar 
Rindi, F. & Guiry, M. Composition and spatial variability of terrestrial algal assemblages occurring at the bases of urban walls in Europe. Phycologia 43, 225–235 (2004).Article 

Google Scholar 
Stonyeva, M. P., Uzunov, B. A. & Gärtner, G. Aerophytic green algae, epimycotic on Fomes fomentarius (L. ex Fr.) Kickx. Annu. Sofia Univ “St. Kliment Ohridski”. Fac. Biol. 99, 19–25 (2015).
Google Scholar 
Aras, S. & Cansaran, D. Isolation of DNA for sequence analysis from herbarium material of some lichen specimens. Turk. J. Bot. 30, 449–453 (2006).
Google Scholar 
Hall, T. BioEdit: A userfriendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 41, 95–98 (1999).CAS 

Google Scholar 
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Posada, D. jModelTest: Phylogenetic model averaging. Mol. Biol. Evol. 25, 1253–1256 (2008).Article 
CAS 
PubMed 

Google Scholar 
Huelsenbeck, J. P. & Ronquist, F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17, 754–755 (2001).Article 
CAS 
PubMed 

Google Scholar 
Vondrák, J. & Kubásek, J. Algal stacks and fungal stacks as adaptations to high light in lichens. Lichenol. 45, 115–124 (2013).Article 

Google Scholar 
Kubásek, J., Hájek, T. & Glime, J. M. Bryophyte photosynthesis in sunflecks: Greater relative induction rate than in tracheophytes. J. Bryol. 36, 110–117 (2014).Article 

Google Scholar 
Kubásek, J. et al. Moss stomata do not respond to light and CO2 concentration but facilitate carbon uptake by sporophytes: A gas exchange, stomatal aperture, and C-13-labelling study. New Phytol. 230, 1815–1828 (2021).Article 
PubMed 

Google Scholar 
Feige, G. & Kremer, B. Unusual carbohydrate pattern in Trentepohlia species. Phytochemistry 19, 1844–1845 (1980).Article 
CAS 

Google Scholar 
Tonon, T., Li, Y. & McQueen-Mason, S. Mannitol biosynthesis in algae: More widespread and diverse than previously thought. New Phytol. 213, 1573–1579 (2017).Article 
PubMed 

Google Scholar 
Gustavs, L., Görs, M. & Karsten, U. Polyol patterns in biofilm-forming aeroterrestrial green algae (Trebouxiophyceae, Chlorophyta). J. Phycol. 47, 533–537 (2011).Article 
PubMed 

Google Scholar  More

Cui, Y. Z. et al. Rapid growth in nitrogen dioxide pollution over Western China, 2005–2013. Atmos. Chem. Phys. 16, 6207–6221. https://doi.org/10.5194/acp-16-6207-2016 (2016).Article 
ADS 
CAS 

Google Scholar 
Gu, J. B. et al. Ground-Level NO2 concentrations over China inferred from the Satellite OMI and CMAQ model simulations. Remote Sens. 9, 519. https://doi.org/10.3390/rs9060519 (2017).Article 
ADS 

Google Scholar 
Cui, Y. Z. et al. Spatio-Temporal heterogeneous impacts of the drivers of NO2 pollution in Chinese cities: Based on satellite observation data. Remote Sens. 14, 3487. https://doi.org/10.3390/rs14143487 (2022).Article 
ADS 

Google Scholar 
Huang, Z. Y., Xu, X. K., Ma, M. G. & Shen, J. W. Assessment of NO2 population exposure from 2005 to 2020 in China. Environ. Sci. Pollut. Res. https://doi.org/10.1007/s11356-022-21420-6 (2022).Article 

Google Scholar 
Zheng, Z. H., Yang, Z. W., Wu, Z. F. & Marinello, F. Spatial variation of NO2 and its impact factors in China: An application of sentinel-5P products. Remote Sens. 11, 1939. https://doi.org/10.3390/rs11161939 (2019).Article 
ADS 

Google Scholar 
Bignal, K. L., Ashmore, M. R., Headley, A. D., Stewart, K. & Weigert, K. Ecological impacts of air pollution from road transport on local vegetation. Appl. Geochem. 22, 1265–1271. https://doi.org/10.1016/j.apgeochem.2007.03.017 (2007).Article 
ADS 
CAS 

Google Scholar 
Zhu, Y. J. et al. Spatiotemporally mapping of the relationship between NO2 pollution and urbanization for a megacity in Southwest China during 2005–2016. Chemosphere 220, 155–162. https://doi.org/10.1016/j.chemosphere.2018.12.095 (2019).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Stieb, D. M. et al. A national study of the association between traffic-related air pollution and adverse pregnancy outcomes in Canada, 1999–2008. Environ. Res. 148, 513–526. https://doi.org/10.1016/j.envres.2016.04.025 (2016).Article 
CAS 
PubMed 

Google Scholar 
Hu, Y. et al. Associations between total mortality and personal exposure to outdoor-originated NO2 in 271 Chinese cities. Atmos. Environ. https://doi.org/10.1016/j.atmosenv.2020.118170 (2021).Article 

Google Scholar 
Han, K. M. Temporal analysis of OMI-Observed tropospheric NO2 columns over east Asia during 2006–2015. Atmosphere 10, 658 (2019).Article 
ADS 
CAS 

Google Scholar 
EEA. Air quality in Europe—2016 report. European Environment Agency EEA Report No 28/2016. Retrieved 2 Dec 2016 from: http://www.eea.europa.eu/publications/air-quality-in-europe-2016Ahmad, A. et al. A comparative study on capability of different tree species in accumulating heavy metals from soil and ambient air. Chemosphere 172, 459–467. https://doi.org/10.1016/j.chemosphere.2017.01.045 (2017).Article 
ADS 
CAS 

Google Scholar 
Erin, R. D., Bryan, K. P., Amy, X. L. & Ronald, C. C. Laboratory measurements of stomatal NO2 deposition to native California trees and the role of forests in the NOx cycle. Atmos. Chem. Phys. 22, 14023–14041. https://doi.org/10.5194/acp-20-14023-2020 (2020).Article 
CAS 

Google Scholar 
Takahashi, M. et al. Differential assimilation of nitrogen dioxide by 70 taxa of roadside trees at an urban pollution level. Chemosphere 61, 633–639. https://doi.org/10.1016/j.chemosphere.2005.03.033 (2005).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Guo, L. L., Li, B. F. & Chen, H. A. A review of urban Micro-climate research on block scale in China. Urban Dev. Stud. 24, 75–81. https://doi.org/10.3969/j.issn.10063862.2017.01.010 (2017).Article 

Google Scholar 
Jung, S. & Yoon, S. Analysis of the effects of floor area ratio change in urban street canyons on microclimate and particulate matter. Energies 14, 714. https://doi.org/10.3390/en14030714 (2021).Article 
CAS 

Google Scholar 
Yin, S. et al. Quantifying air pollution attenuation within urban parks: An experimental approach in Shanghai, China. Environ. Pollut. 159, 2155–2163. https://doi.org/10.1016/j.envpol.2011.03.009 (2011).Article 
CAS 
PubMed 

Google Scholar 
Lin, C., Feng, X. F. & Heal, M. R. Temporal persistence of intra-urban spatial contrasts in ambient NO2, O3 and Ox in Edinburgh, UK. Atmos. Pollut. Res. 7, 734–741. https://doi.org/10.1016/j.apr.2016.03.008 (2016).Article 

Google Scholar 
Brantley, H. L., Hagler, G. S. W., Deshmukh, P. J. & Baldauf, R. W. Field assessment of the effects of roadside vegetation on near-road black carbon and particulate matter. Sci. Total Environ. 468, 120–129. https://doi.org/10.1016/j.scitotenv.2013.08.001 (2014).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Irga, P. J., Burchett, M. D. & Torpy, F. R. Does urban forestry have a quantitative effect on ambient air quality in an urban environment?. Atmos. Environ. 120, 173–181. https://doi.org/10.1016/j.atmosenv.2015.08.050 (2015).Article 
ADS 
CAS 

Google Scholar 
Tong, Z. M., Baldauf, R. W., Isakov, V., Deshmunk, P. & Zhang, K. M. Roadside vegetation barrier design to mitigate near-road air pollution impacts. Sci. Total Environ. 541, 920–927. https://doi.org/10.1016/j.scitotenv.2015.09.067 (2016).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Setälä, H., Viippola, V., Rantalainen, A. L., Pennanen, A. & Yli-Pelkonen, V. Does urban vegetation mitigate air pollution in northern conditions?. Environ. Pollut. 183, 104–112. https://doi.org/10.1016/j.envpol.2012.11.010 (2013).Article 
CAS 
PubMed 

Google Scholar 
Xing, Y. & Brimblecombe, P. Role of vegetation in deposition and dispersion of air pollution in urban parks. Atmos. Environ. 201, 73–83. https://doi.org/10.1016/j.atmosenv.2018.12.027 (2019).Article 
ADS 
CAS 

Google Scholar 
Xu, C., Wang, Y. P. & Li, L. L. Study on spatiotemporal distribution of the tropospheric NO2 column concentration in China and its relationship to energy consumption based on the time-series data from 2005 to 2013. Energy Sources Part A 42, 2130–2144. https://doi.org/10.1080/15567036.2019.1607931 (2020).Article 
CAS 

Google Scholar 
Xu, J. H., Lindqvist, H., Liu, Q. F., Wang, K. & Wang, L. Estimating the spatial and temporal variability of the ground-level NO2 concentration in China during 2005–2019 based on satellite remote sensing. Atmos. Pollut. Res. 12, 57–67. https://doi.org/10.1016/j.apr.2020.10.008 (2021).Article 
CAS 

Google Scholar 
Daniel, L. G. et al. TROPOMI NO2 in the United States: A detailed look at the annual averages, weekly cycles, effects of temperature, and correlation with surface NO2 concentrations. Earths Feature 9, 4. https://doi.org/10.1029/2020EF001665 (2021).Article 
CAS 

Google Scholar 
Mavroidis, I. & Chaloulakou, A. Long-term trends of primary and secondary NO2 production in the Athens area. Variation of the NO2/NOx ratio. Atmos. Environ. 45, 6872–6879. https://doi.org/10.1016/j.atmosenv.2010.11.006 (2011).Article 
ADS 
CAS 

Google Scholar 
Van der, A. R. J. et al. Detection of the trend and seasonal variation in tropospheric NO2 over China. J. Geophys. Res. Atmos. https://doi.org/10.1029/2005JD006594 (2006).Article 

Google Scholar 
Salama, D. S. et al. Satellite observations for monitoring atmospheric NO2 in correlation with the existing pollution sources under arid environment. Model. Earth Syst. Environ. 8, 4103–4121. https://doi.org/10.1007/s40808-022-01352-3 (2022).Article 
PubMed 

Google Scholar 
Ahmad, S. S. & Aziz, N. Spatial and temporal analysis of ground level ozone and nitrogen dioxide concentration across the twin cities of Pakistan. Environ. Monit. Assess. 185, 3133–3147. https://doi.org/10.1007/s10661-012-2778-7 (2013).Article 
CAS 
PubMed 

Google Scholar 
Khaled, G., Abdulaziz, A., Watheq, A. & Mumin, A. Analysis of NOx, NO and NO2 ambient levels in Dhahran, Saudi Arabia. Urban Clim. 21, 232–242. https://doi.org/10.2495/AIR170081 (2017).Article 

Google Scholar 
Casquero-Vera, J. A. et al. Impact of primary NO2 emissions at different urban sites exceeding the European NO2 standard limit. Sci. Total Environ. 646, 1117–1125 (2019).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Desyana, R. D., Sulistyantara, B., Nasrullah, N. & Fatimah, I. S. Study of the effectiveness of several tree canopy types on roadside green belt in influencing the distribution of NO2 gas emitted from transportation. EES https://doi.org/10.1088/1755-1315/58/1/012045 (2017).Article 

Google Scholar 
Rotach, M. W. Profiles of turbulence statistics in and above an urban street canyon. Atmos. Environ. 29, 1473–1486. https://doi.org/10.1016/1352-2310(95)00084-C (1995).Article 
ADS 
CAS 

Google Scholar 
Luo, M. Study on Air Pollutants Removal Effects of Green Space with Different Community Structures (Huazhong Agricultural University, 2013).
Google Scholar 
Rao, M., George, L. A., Rosenstiel, T. N., Shandas, V. & Dinno, A. Assessing the relationship among urban trees, nitrogen dioxide, and respiratory health. Environ. Pollut. 194, 96–104. https://doi.org/10.1016/j.envpol.2014.07.011 (2014).Article 
CAS 
PubMed 

Google Scholar 
Yli-Pelkonen, V., Viippola, V., Kotze, D. J. & Setala, H. Greenbelts do not reduce NO2 concentrations in near-road environments. Urban Clim. 21, 306–317. https://doi.org/10.1016/j.uclim.2017.08.005 (2017).Article 

Google Scholar 
Fantozzi, F., Monaci, F., Blanusa, T. & Bargagli, R. Spatio-temporal variations of ozone and nitrogen dioxide concentrations under urban trees and in a nearby open area. Urban Clim. 12, 119–127. https://doi.org/10.1016/j.uclim.2015.02.001 (2015).Article 

Google Scholar 
Nie, L., Deng, Z. H. & Chen, Q. B. SO2 and NOx purify-cation ability of forest in Kunming City. J. West China For. Sci. 44, 116–120 (2015).
Google Scholar 
Baldauf, R. Roadside vegetation design characteristics that can improve local, near-road air quality. Transp. Res. Part D 52, 354–361. https://doi.org/10.1016/j.trd.2017.03.013 (2017).Article 

Google Scholar 
Lai, D. Y., Liu, Y. Q., Liao, M. C. & Yu, B. Q. Effects of different tree layouts on outdoor thermal comfort of green space in summer Shanghai. Urban Clim. 47, 101398 (2023).Article 

Google Scholar 
Lai, D., Liu, W., Gan, T., Liu, K. & Chen, Q. A review of mitigating strategies to improve the thermal environment and thermal comfort in urban outdoor spaces. Sci. Total Environ. 661, 337–353 (2019).Article 
ADS 
CAS 
PubMed 

Google Scholar  More

In co-culture with the bacterivorous flagellate Poteriospumella lacustris, the prey bacterium Pseudomonas putida exhibited a characteristic succession of predation defenses. The initial and the final defense differed substantially from one another with regard to their mechanism and their population-level benefits to the bacteria.Our results strongly indicate that the initial bacterial defense falls into the category of chemical defense, and is regulated by phenotypic plasticity. This would require P. putida to be able to sense predator density and to regulate the excretion of inhibitory substances accordingly. Because a considerable proportion of the P. putida genome is known to be involved in regulation and signal transduction allowing for very flexible responses to environmental triggers [41] both conditions are likely to be met. The filtrate exposure tests (Fig. 3) provide specific evidence for the ability of P. putida KT2440 to up- and downregulate the excretion of compounds inhibiting flagellate growth in response to grazing pressure. Previous research [25] corroborated the ability of P. putida to escape grazing from bacterivorous flagellates through induced responses like aggregation or biofilm formation.To provide a possible characterization for the apparent bacterial toxin, the whole-genome sequences of P. putida KT2440 obtained here were aligned against the antiSMASH [42] database. The output suggests the existence of non-ribosomal peptide synthetase clusters mediating the production of pyoverdines, a particular class of siderophores. The latter are molecules released by bacteria into the environment, which enhance the uptake of essential metals like, e.g., iron under deficient conditions. Specific pyoverdines associated with P. putida KT2440 have previously been identified [43]. Recent findings have shown that the benefits from siderophore production are not limited to competitive advantages gained from enhanced resource exploitation [44]. Pyoverdines were also demonstrated to determine the virulence of Pseudomonads via the damage of mitochondria in colonized hosts [45]. Moreover, pyoverdines were shown to be involved in the inducible defense of P. putida against predatory myxobacteria [46]. Such multiple functions have been reported for a number of bacterial metabolites, especially in Pseudomonads [47], and the particular combination of pyoverdin effects would explain the observed simultaneous flagellate inhibition and promoted bacterial growth.In contrast to the initial chemical defense of P. putida, the subsequent filamentation clearly provides an example of rapid evolution. Although the responsible mutation(s) could only be pinpointed in a few isolates so far (Table S1), there is no doubt about the genetic manifestation and heritability of the filamentous phenotype due to its demonstrated non-reversible nature.Only recently, similar observations were made by long-term co-cultivation of Pseudomonas fluorescence with the amoeboid predator Neaglena grubei [48]. In that system, protective adaptations like enhanced biofilm formation and altered motility were traced down to mutations in two particular genes (wspF, amrZ).From the perspective of the bacterial population, filamentation appears to be a much less efficient defense mechanism than toxin production. This is clearly reflected by the ratio of prey to predator biomass, which differed by two orders of magnitude between the initial and final defense (Table 6). It raises the question of why bacteria would abandon a highly effective form of defense in favor of a much less effective one. As demonstrated experimentally, adaptation of predators to the toxin can be excluded as a cause (Fig. 4). Moreover, it was not instantly evident how the small-sized flagellate was ultimately able to persist in large numbers given a very high proportion of completely inedible prey individuals (Fig. 1D and Fig. S2).Table 6 Average abundance of predator and prey during the temporary steady state following the initial bacterial defense (day 13–16) and during the final steady state (beyond day 30).Full size tableTo develop a comprehensive understanding of the system addressing the questions raised above, we set up a semi-continuous differential equation model to simulate the dynamics of predator and prey phenotypes. The model considers seven state variables (carbon, densities of four bacterial phenotypes, flagellate density, and toxin concentration) whose dynamics are controlled by nine processes (Table 3, Fig. 2). In addition to microbial growth and grazing, the model implements a phenotypically plastic predation defense (toxin production) as well as a genetic defense (filamentation) which arises via mutation. The particular assumptions implemented in the model are as follows:Dual effect of bacterial metabolitesIn line with the above discussion on siderophore-like compounds, secondary metabolites excreted by P. putida were assumed to exhibit a dual function, both inhibiting the growth of flagellates and allowing for a more efficient exploitation of the resources by bacteria. The inhibition of predators was demonstrated directly (Figs. 3 and 4) while enhanced resource exploitation was inferred from bacterial abundances in co-cultures exceeding the carrying capacity observed in predator-free controls (Fig. 1A, day 11–18).Metabolite production is costlyThe production of bacterial metabolites was assumed to be associated with a slight fitness cost [49] since resources are diverted from reproduction, thus resulting in a lowered growth rate of toxin-producing bacteria. The assumed fitness cost of 11% (parameter cBx in Table 5) allowed for the best agreement between simulated and observed data and is in agreement with data on the cost of pyoverdine production by P. aeruginosa [50]. The cost only manifests when toxin production is upregulated.Predator recognition and quorum sensing interactIn the model, the production of bacterial metabolites is upregulated when the two conditions of high flagellate abundance and high bacterial abundance coincide. That is, the expression of the toxin-based bacterial defense is assumed to be jointly controlled by predator recognition and quorum sensing (QS). Examples for such joint control of bacterial defenses have been reported previously [8, 26, 51]. The involvement of QS in chemical defense strategies is particularly likely as effective toxin concentrations can only be reached when producers are highly abundant. While multiple QS systems have been described for other Pseudomonads, only a single system has been identified in P. putida KT2440 so far [52, 53].Mutation rates are conditional on stressThe emergence of mutations resulting in the filamentation of P. putida was assumed to be conditional on a high ambient concentration of bacterial metabolites. The latter was considered as a proxy for bacterial stress which can affect mutagenesis either directly or indirectly by a variety of mechanisms [54,55,56]. Without this assumption, the almost synchronous appearance of filaments in all replicates at a late point in time would be very difficult to explain. Specifically, if mutation frequencies were high, filaments would become the predominant phenotype early (Fig. S3) which contradicts observations. On the other hand, if frequencies were low but unconditional, the timing of filament appearance should vary between replicates, which is in contrast to observations either (Fig. 1B).Filamentation is associated with a fitness costMeasurements of growth rate constants revealed a significant fitness disadvantage of filamentous isolates in comparison to single-celled, undefended isolates (p  More

An, Z., Guo, F., Chen, Y., Bai, G. & Chen, Z. Rhizosphere bacterial and fungal communities during the growth of Angelica sinensis seedlings cultivated in an Alpine uncultivated meadow soil. PeerJ 8, e8541. https://doi.org/10.7717/peerj.8541 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Munkholm, L. J., Heck, R. J. & Deen, B. Long-term rotation and tillage effects on soil structure and crop yield. Soil Tillage Res. 127, 85–91. https://doi.org/10.1016/j.still.2012.02.007 (2013).Article 

Google Scholar 
Jiao, X. L. et al. Effects of maize rotation on the physicochemical properties and microbial communities of American ginseng cultivated soil. Sci. Rep. 9, 8615. https://doi.org/10.1038/s41598-019-44530-7 (2019).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Wang, X., Chen, Y., Guo, F., Yuan, H. & Guo, Y. Effects of medicinal crop stubbles on physiological and biochemical characteristics of Angelica sinensis seedings. J. Chin. Med. Mater. 40, 2002–2006 (2017).
Google Scholar 
Jin, Y. et al. Effect of various crop residues on growth and disease resisitance of Angelica sinensis seedlings in Min County. Acta Pratacul. Sin. 27, 69–78 (2018).MathSciNet 

Google Scholar 
Bai, G., Guo, F., Chen, Y., Yuan, H. & Xiao, W. Differences in physiological resistance traits of Angelica sinensis seedlings from uncultivated and cultivated fields in Min County. Acta Pratacul. Sin. 28, 86–95 (2019).
Google Scholar 
Bai, G. et al. Regulated effects of preceding crop on soil property and cultivating seedlings for Angelica sinensis on cultivated farmland. Chin. J. Eco-Agric. 28, 701–712. https://doi.org/10.13930/j.cnki.cjea.190719 (2020).Article 
CAS 

Google Scholar 
Mendes, R., Garbeva, P. & Raaijmakers, J. M. The rhizosphere microbiome: Significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 37, 634–663 (2013).Article 
CAS 
PubMed 

Google Scholar 
Tkacz, A., Cheema, J., Chandra, G., Grant, A. & Poole, P. S. Stability and succession of the rhizosphere microbiota depends upon plant type and soil composition. Int. Soc. Microb. Ecol. 9, 2349–2359. https://doi.org/10.1038/ismej.2015.41 (2015).Article 
CAS 

Google Scholar 
Berg, G. et al. Microbiome definition re-visited: Old concepts and new challenges. Microbiome 8, 103. https://doi.org/10.1186/s40168-020-00875-0 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Chaparro, J. M., Badri, D. V. & Vivanco, J. M. Rhizosphere microbiome assemblage is affected by plant development. ISME J. 8, 790–803. https://doi.org/10.1038/ismej.2013.196 (2014).Article 
CAS 
PubMed 

Google Scholar 
Uroz, S. et al. Specific impacts of beech and Norway spruce on the structure and diversity of the rhizosphere and soil microbial communities. Sci. Rep. 6, 27756. https://doi.org/10.1038/srep27756 (2016).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Chamberlain, L. A. et al. Crop rotation, but not cover crops, influenced soil bacterial community composition in a corn-soybean system in southern Wisconsin. Appl. Soil Ecol. 154, 103603. https://doi.org/10.1016/j.apsoil.2020.103603 (2020).Article 

Google Scholar 
Classen, A. T. et al. Direct and indirect effects of climate change on soil microbial and soil microbial-plant interactions: What lies ahead?. Ecosphere 6, 130. https://doi.org/10.1890/es15-00217.1 (2015).Article 

Google Scholar 
Tiemann, L. K. et al. Crop rotational diversity enhances belowground communities and functions in an agroecosystem. Ecol. Lett. 18, 761–771. https://doi.org/10.1111/ele.12453 (2015).Article 
CAS 
PubMed 

Google Scholar 
Maldonado, S. et al. Enhanced crop productivity and sustainability by using native phosphate solubilizing rhizobacteria in the agriculture of arid zones. Front. Sustain. Food Syst. 4, 607355. https://doi.org/10.3389/fsufs.2020.607355 (2020).Article 

Google Scholar 
Gómez Expósito, R., de Bruijn, I., Postma, J. & Raaijmakers, J. M. Current insights into the role of rhizosphere bacteria in disease suppressive soils. Front. Microbiol.y 8, 2529. https://doi.org/10.3389/fmicb.2017.02529 (2017).Article 

Google Scholar 
Li, X., Rui, J., Mao, Y., Yannarell, A. & Mackie, R. Dynamics of the bacterial community structure in the rhizosphere of a maize cultivar. Soil Biol. Biochem. 68, 392–401. https://doi.org/10.1016/j.soilbio.2013.10.017 (2014).Article 
CAS 

Google Scholar 
Fierer, N. et al. Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. Int. Soc. Microb. Ecol. 6, 1007–1017. https://doi.org/10.1038/ismej.2011.159 (2012).Article 
CAS 

Google Scholar 
Kuffner, M. et al. Culturable bacteria from Zn- and Cd-accumulating Salix caprea with differential effects on plant growth and heavy metal availability. J. Appl. Microbiol. 108, 1471–1484. https://doi.org/10.1111/j.1365-2672.2010.04670.x (2010).Article 
CAS 
PubMed 

Google Scholar 
De Corato, U. Disease-suppressive compost enhances natural soil suppressiveness against soil-borne plant pathogens: A critical review. Rhizosphere 13, 100192. https://doi.org/10.1016/j.rhisph.2020.100192 (2020).Article 

Google Scholar 
Brookes, P. C., Landman, A., Pruden, G. & Jenkinson, D. S. Chloroform fumigation chloroform fumigation and the release of soil nitrogen: a rapid direct extraction method to measure microbial biomass nitrogen in soil. Soil Biol. Biochem. 17, 837–842 (1985).Article 
CAS 

Google Scholar 
Arnebrant, K. & Schnürer, J. Changes in atp content during and after chloroform fumigation. Soil Biol. Biochem. 22, 875–877 (1990).Article 
CAS 

Google Scholar 
Toju, H. et al. Community composition of root-associated fungi in a Quercus-dominated temperate forest: “codominance” of mycorrhizal and root-endophytic fungi. Ecol. Evol. 3, 1281–1293. https://doi.org/10.1002/ece3.546 (2013).Article 
PubMed 
PubMed Central 

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

Google Scholar 
Bokulich, N. A. et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 10, 57–59. https://doi.org/10.1038/nmeth.2276 (2013).Article 
CAS 
PubMed 

Google Scholar 
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200 (2011).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Haas, B. J. et al. Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21, 494–504. https://doi.org/10.1101/gr.112730.110 (2011).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Edgar, R. C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–998. https://doi.org/10.1038/nmeth.2604 (2013).Article 
CAS 
PubMed 

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

Google Scholar 
Quast, C. et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. https://doi.org/10.1093/nar/gks1219 (2013).Article 
CAS 
PubMed 

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

Google Scholar 
Louca, S., Parfrey, L. W. & Doebeli, M. Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272 (2016).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Nguyen, N. H. et al. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. https://doi.org/10.1016/j.funeco.2015.06.006 (2016).Article 

Google Scholar 
Sisk-Hackworth, L., Ortiz-Velez, A., Reed, M. B. & Kelley, S. T. Compositional data analysis of periodontal disease microbial communities. Front. Microbiol. 12, 617949. https://doi.org/10.3389/fmicb.2021.617949 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Khan, M. A. W. et al. Deforestation impacts network co-occurrence patterns of microbial communities in Amazon soils. FEMS Microbiol. Ecol. 95, fiy230. https://doi.org/10.1093/femsec/fiy230 (2019).Article 
CAS 
PubMed 

Google Scholar 
Zhang, B., Zhang, J., Liu, Y., Shi, P. & Wei, G. Co-occurrence patterns of soybean rhizosphere microbiome at a continental scale. Soil Biol. Biochem. 118, 178–186. https://doi.org/10.1016/j.soilbio.2017.12.011 (2018).Article 
CAS 

Google Scholar 
Huang, M., Jiang, L., Zou, Y., Xu, S. & Deng, G. Changes in soil microbial properties with no-tillage in Chinese cropping systems. Biol. Fertil. Soils 49, 373–377. https://doi.org/10.1007/s00374-013-0778-6 (2013).Article 

Google Scholar 
Unger, P. W. & Cassel, D. K. Tillage implement disturbance effects on soil properties related to soil and water conservation: A literature review. Soil Tillage Res. 19, 363–382 (1991).Article 

Google Scholar 
Alvarez, R. & Steinbach, H. S. A review of the effects of tillage systems on some soil physical properties, water content, nitrate availability and crops yield in the Argentine Pampas. Soil Tillage Res. 104, 1–15. https://doi.org/10.1016/j.still.2009.02.005 (2009).Article 

Google Scholar 
Essel, E. et al. Bacterial and fungal diversity in rhizosphere and bulk soil under different long-term tillage and cereal/legume rotation. Soil Tillage Res. 194, 104302. https://doi.org/10.1016/j.still.2019.104302 (2019).Article 

Google Scholar 
Zhu, Q., Wang, N., Duan, B., Wang, Q. & Wang, Y. Rhizosphere bacterial and fungal communities succession patterns related to growth of poplar fine roots. Sci. Total Environ. 756, 143839. https://doi.org/10.1016/j.scitotenv.2020.143839 (2021).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Guseva, K. et al. From diversity to complexity: Microbial networks in soils. Soil Biol. Biochem. 169, 108604. https://doi.org/10.1016/j.soilbio.2022.108604 (2022).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Jiang, B. et al. Analysis of microbial community structure and diversity in surrounding rock soil of different waste dump sites in fushun western opencast mine. Chemosphere 269, 128777. https://doi.org/10.1016/j.chemosphere.2020.128777 (2020).Article 
ADS 
MathSciNet 
CAS 
PubMed 

Google Scholar 
Liu, J. et al. Pecan plantation age influences the structures, ecological networks, and functions of soil microbial communities. Land Degrad. Dev. 33, 3294–3309. https://doi.org/10.1002/ldr.4389 (2022).Article 

Google Scholar 
Lv, X. et al. Strengthening insights in microbial ecological networks from theory to applications. mSystems 4, e00124-19. https://doi.org/10.1128/mSystems.00124-19 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Toju, H., Kishida, O., Katayama, N. & Takagi, K. Networks depicting the fine-scale co-occurrences of fungi in soil Horizons. PLoS ONE 11, e0165987. https://doi.org/10.1371/journal.pone.0165987 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Chun, S. J., Cui, Y., Baek, S. H., Ahn, C. Y. & Oh, H. M. Seasonal succession of microbes in different size-fractions and their modular structures determined by both macro- and micro-environmental filtering in dynamic coastal waters. Sci. Total Environ. 784, 147046. https://doi.org/10.1016/j.scitotenv.2021.147046 (2021).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Cardinale, M., Grube, M., Erlacher, A., Quehenberger, J. & Berg, G. Bacterial networks and co-occurrence relationships in the lettuce root microbiota. Environ. Microbiol. 17, 239–252. https://doi.org/10.1111/1462-2920.12686 (2015).Article 
CAS 
PubMed 

Google Scholar 
Zhou, Z. et al. Increases in bacterial community network complexity induced by biochar-based fertilizer amendments to karst calcareous soil. Geoderma 337, 691–700. https://doi.org/10.1016/j.geoderma.2018.10.013 (2019).Article 
ADS 
CAS 

Google Scholar 
Olesen, J. M., Bascompte, J., Dupont, Y. L. & Jordano, P. The modularity of pollination networks. Proc. Natl. Acad. Sci. U.S.A. 104, 19891–19896. https://doi.org/10.1073/pnas.0706375104 (2007).Article 
ADS 
PubMed 
PubMed Central 
MATH 

Google Scholar 
Eisenhauer, N. et al. Root biomass and exudates link plant diversity with soil bacterial and fungal biomass. Sci. Rep. 7, 44641. https://doi.org/10.1038/srep44641 (2017).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Hassan, M. K., McInroy, J. A. & Kloepper, J. W. The interactions of rhizodeposits with plant growth-promoting Rhizobacteria in the rhizosphere: A review. Agriculture 9, 142. https://doi.org/10.3390/agriculture9070142 (2019).Article 
CAS 

Google Scholar 
Sasse, J., Martinoia, E. & Northen, T. Feed your friends: Do plant exudates shape the root microbiome?. Trends Plant Sci. 23, 25–41. https://doi.org/10.1016/j.tplants.2017.09.003 (2018).Article 
CAS 
PubMed 

Google Scholar 
Zhang, F., Xu, X., Wang, G., Wu, B. & Xiao, Y. Medicago sativa and soil microbiome responses to Trichoderma as a biofertilizer in alkaline-saline soils. Appl. Soil Ecol. 153, 103573. https://doi.org/10.1016/j.apsoil.2020.103573 (2020).Article 

Google Scholar 
Woźniak, A. Chemical properties and enzyme activity of soil as affected by tillage system and previous crop. Agriculture 9, 262. https://doi.org/10.3390/agriculture9120262 (2019).Article 
CAS 

Google Scholar 
Choudhary, M. et al. Changes in soil biology under conservation agriculture based sustainable intensification of cereal systems in Indo-Gangetic Plains. Geoderma 313, 193–204. https://doi.org/10.1016/j.geoderma.2017.10.041 (2018).Article 
ADS 
CAS 

Google Scholar 
Ai, C. et al. Distinct responses of soil bacterial and fungal communities to changes in fertilization regime and crop rotation. Geoderma 319, 156–166. https://doi.org/10.1016/j.geoderma.2018.01.010 (2018).Article 
ADS 
CAS 

Google Scholar 
Gałązka, A., Gawyjołek, K., Perzyński, A., Gałązka, R. & Jerzy, K. Changes in enzymatic activities and microbial communities in soil under long-term maize monoculture and crop rotation. Pol. J. Environ. Stud. 26, 39–46. https://doi.org/10.15244/pjoes/64745 (2017).Article 
CAS 

Google Scholar 
Tremblay, C., Deslauriers, A., Lafond, J., Lajeunesse, J. & Paré, M. Effects of soil pH and fertilizers on haskap (Lonicera caerulea L) vegetative growth. Agriculture 9, 56. https://doi.org/10.3390/agriculture9030056 (2019).Article 
CAS 

Google Scholar 
Sirisuntornlak, N. et al. Interactive effects of silicon and soil pH on growth, yield and nutrient uptake of maize. SILICON 13, 289–299. https://doi.org/10.1007/s12633-020-00427-z (2021).Article 
CAS 

Google Scholar 
Xu, Y., Ge, Y., Song, J. & Rensing, C. Assembly of root-associated microbial community of typical rice cultivars in different soil types. Biol. Fertil. Soils 56, 249–260. https://doi.org/10.1007/s00374-019-01406-2 (2019).Article 
CAS 

Google Scholar 
Putranta, H., Permatasari, A. K., Sukma, T. A. & Dwandaru, W. S. B. The effect of pH, electrical conductivity, and nitrogen (N) in the soil at yogyakarta special region on tomato plant growth. TEM J.-Technol. Educ. Manag. Inform. 8, 860–865. https://doi.org/10.18421/TEM83-24 (2019).Article 

Google Scholar 
Wang, J. et al. Effects of alternate partial root-zone irrigation on soil microorganism and maize growth. Plant Soil 302, 45–52. https://doi.org/10.1007/s11104-007-9453-8 (2007).Article 
CAS 

Google Scholar 
Yang, X., Zhu, K., Loik, M. E. & Sun, W. Differential responses of soil bacteria and fungi to altered precipitation in a meadow steppe. Geoderma 384, 114812. https://doi.org/10.1016/j.geoderma.2020.114812 (2021).Article 
ADS 
CAS 

Google Scholar 
Balota, E. L., Colozzi Filho, A., Andrade, D. S. & Dick, R. P. Long-term tillage and crop rotation effects on microbial biomass and C and N mineralization in a Brazilian Oxisol. Soil Tillage Res. 77, 137–145. https://doi.org/10.1016/j.still.2003.12.003 (2004).Article 

Google Scholar 
Franchini, J., Crispino, C., Souza, R., Torres, E. & Hungria, M. Microbiological parameters as indicators of soil quality under various soil management and crop rotation systems in southern Brazil. Soil Tillage Res. 92, 18–29. https://doi.org/10.1016/j.still.2005.12.010 (2007).Article 

Google Scholar 
Li, X., Wang, T., Chang, S. X., Jiang, X. & Song, Y. Biochar increases soil microbial biomass but has variable effects on microbial diversity: A meta-analysis. Sci. Total Environ. 749, 141593. https://doi.org/10.1016/j.scitotenv.2020.141593 (2020).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Lynch, J. M. & Panting, L. M. Effects of season, cultivation and nitrogen fertiliser on the size of the soil microbial biomass. J. Sci. Food Agric. 33, 249–252 (1982).Article 
CAS 

Google Scholar 
Tan, G. et al. Effects of biochar application with fertilizer on soil microbial biomass and greenhouse gas emissions in a peanut cropping system. Environ. Technol. 42, 9–19. https://doi.org/10.1080/09593330.2019.1620344 (2021).Article 
CAS 
PubMed 

Google Scholar 
Liu, C. et al. Linkages between nutrient ratio and the microbial community in rhizosphere soil following fertilizer management. Environ. Res. 184, 109261. https://doi.org/10.1016/j.envres.2020.109261 (2020).Article 
CAS 
PubMed 

Google Scholar 
Li, H. et al. Film mulching, residue retention and N fertilization affect ammonia volatilization through soil labile N and C pools. Agric. Ecosyst. Environ. 308, 107272. https://doi.org/10.1016/j.agee.2020.107272 (2021).Article 
CAS 

Google Scholar 
Jiao, P. et al. Bacteria are more sensitive than fungi to moisture in eroded soil by natural grass vegetation restoration on the Loess Plateau. Sci. Total Environ. 756, 143899. https://doi.org/10.1016/j.scitotenv.2020.143899 (2021).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Sommer, J. et al. The tree species matters: Belowground carbon input and utilization in the myco-rhizosphere. Eur. J. Soil Biol. 81, 100–107. https://doi.org/10.1016/j.ejsobi.2017.07.001 (2017).Article 
CAS 

Google Scholar 
Yu, K., Pieterse, C. M. J., Bakker, P. A. H. M. & Berendsen, R. L. Beneficial microbes going underground of root immunity. Plant Cell Environ. 42, 2860–2870. https://doi.org/10.1111/pce.13632 (2019).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
de Varennes, A. & Goss, M. J. The tripartite symbiosis between legumes, rhizobia and indigenous mycorrhizal fungi is more efficient in undisturbed soil. Soil Biol. Biochem. 39, 2603–2607. https://doi.org/10.1016/j.soilbio.2007.05.007 (2007).Article 
CAS 

Google Scholar 
Wang, X. et al. Mycorrhizal symbiosis modulates the rhizosphere microbiota to promote rhizobia-legume symbiosis. Mol. Plant 14, 503–516. https://doi.org/10.1016/j.molp.2020.12.002 (2020).Article 
CAS 
PubMed 

Google Scholar 
Zhang, R., Vivanco, J. M. & Shen, Q. The unseen rhizosphere root-soil-microbe interactions for crop production. Curr. Opin. Microbiol. 37, 8–14. https://doi.org/10.1016/j.mib.2017.03.008 (2017).Article 
PubMed 

Google Scholar 
Berendsen, R. L., Pieterse, C. M. & Bakker, P. A. The rhizosphere microbiome and plant health. Trends Plant Sci. 17, 478–486. https://doi.org/10.1016/j.tplants.2012.04.001 (2012).Article 
CAS 
PubMed 

Google Scholar  More