Philippa Kaur

1.Sogin, M. L. et al. Microbial diversity in the deep sea and the underexplored “rare biosphere”. P Natl Acad. Sci. USA 103, 12115–12120 (2006).CAS 
ADS 

Google Scholar 
2.Pedros-Alio, C. The rare bacterial biosphere. Ann. Rev. Mar. Sci. 4, 449–466 (2012).PubMed 

Google Scholar 
3.Lynch, M. D. J. & Neufeld, J. D. Ecology and exploration of the rare biosphere. Nat. Rev. Microbiol. 13, 217–229 (2015).CAS 
PubMed 

Google Scholar 
4.Campbell, B. J., Yu, L. Y., Heidelberg, J. F. & Kirchman, D. L. Activity of abundant and rare bacteria in a coastal ocean. P Natl Acad. Sci. USA 108, 12776–12781 (2011).CAS 
ADS 

Google Scholar 
5.Gobet, A. et al. Diversity and dynamics of rare and of resident bacterial populations in coastal sands. Isme J. 6, 542–553 (2012).PubMed 

Google Scholar 
6.Wilhelm, L. et al. Rare but active taxa contribute to community dynamics of benthic biofilms in glacier-fed streams. Environ. Microbiol. 16, 2514–2524 (2014).CAS 
PubMed 

Google Scholar 
7.Lawson, C. E. et al. Rare taxa have potential to make metabolic contributions in enhanced biological phosphorus removal ecosystems. Environ. Microbiol. 17, 4979–4993 (2015).CAS 
PubMed 

Google Scholar 
8.Newton, R. J. & Shade, A. Lifestyles of rarity: understanding heterotrophic strategies to inform the ecology of the microbial rare biosphere. Aquat. Micro. Ecol. 78, 51–63 (2016).
Google Scholar 
9.Lauro, F. M. et al. The genomic basis of trophic strategy in marine bacteria. P Natl Acad. Sci. USA 106, 15527–15533 (2009).CAS 
ADS 

Google Scholar 
10.Klappenbach, J. A., Dunbar, J. M. & Schmidt, T. M. RRNA operon copy number reflects ecological strategies of bacteria. Appl Environ. Micro. 66, 1328–1333 (2000).CAS 
ADS 

Google Scholar 
11.Roller, B. R. K., Stoddard, S. F. & Schmidt, T. M. Exploiting rRNA operon copy number to investigate bacterial reproductive strategies. Nat. Microbiol. 1, 16160 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
12.Polz, M. F. & Cordero, O. X. Bacterial evolution: genomics of metabolic trade-offs. Nat. Microbiol. 1, 16181 (2016).CAS 
PubMed 

Google Scholar 
13.Giovannoni, S. J. SAR11 bacteria: the most abundant plankton in the oceans. Ann. Rev. Mar. Sci. 9, 231–255 (2017).PubMed 

Google Scholar 
14.Elser, J. J. et al. Biological stoichiometry from genes to ecosystems. Ecol. Lett. 3, 540–550 (2000).
Google Scholar 
15.Hessen, D. O., Elser, J. J., Sterner, R. W. & Urabe, J. Ecological stoichiometry: an elementary approach using basic principles. Limnol. Oceanogr. 58, 2219–2236 (2013).CAS 
ADS 

Google Scholar 
16.Acharya, K., Kyle, M. & Elser, J. J. Biological stoichiometry of Daphnia growth: an ecophysiological test of the growth rate hypothesis. Limnol. Oceanogr. 49, 656–665 (2004).CAS 
ADS 

Google Scholar 
17.Hendrixson, H. A., Sterner, R. W. & Kay, A. D. Elemental stoichiometry of freshwater fishes in relation to phylogeny, allometry and ecology. J. Fish. Biol. 70, 121–140 (2007).
Google Scholar 
18.Matzek, V. & Vitousek, P. M. N: P stoichiometry and protein: RNA ratios in vascular plants: an evaluation of the growth-rate hypothesis. Ecol. Lett. 12, 765–771 (2009).PubMed 

Google Scholar 
19.Ghoul, M. & Mitri, S. The ecology and evolution of microbial competition. Trends Microbiol 24, 833–845 (2016).CAS 
PubMed 

Google Scholar 
20.Laland, K., Matthews, B. & Feldman, M. W. An introduction to niche construction theory. Evol. Ecol. 30, 191–202 (2016).PubMed 
PubMed Central 

Google Scholar 
21.Ratzke, C., Barrere, J. & Gore, J. Strength of species interactions determines biodiversity and stability in microbial communities. Nat. Ecol. Evol. 4, 376–383 (2020).PubMed 

Google Scholar 
22.Větrovský, T. & Baldrian, P. The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses. Plos ONE 8, e57923 (2013).PubMed 
PubMed Central 
ADS 

Google Scholar 
23.Molenaar, D., van Berlo, R., de Ridder, D. & Teusink, B. Shifts in growth strategies reflect tradeoffs in cellular economics. Mol. Syst. Biol. 5, 323–323 (2009).PubMed 
PubMed Central 

Google Scholar 
24.Nemergut, D. R. et al. Decreases in average bacterial community rRNA operon copy number during succession. Isme J. 10, 1147–1156 (2016).CAS 
PubMed 

Google Scholar 
25.Wu, L. W. et al. Microbial functional trait of rRNA operon copy numbers increases with organic levels in anaerobic digesters. Isme J. 11, 2874–2878 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
26.Wu, L. W. et al. Global diversity and biogeography of bacterial communities in wastewater treatment plants. Nat. Microbiol. 4, 2579–2579 (2019).PubMed 

Google Scholar 
27.Stoddard, S. F., Smith, B. J., Hein, R., Roller, B. R. K. & Schmidt, T. M. rrnDB: improved tools for interpreting rRNA gene abundance in bacteria and archaea and a new foundation for future development. Nucleic Acids Res. 43, D593–D598 (2015).CAS 
PubMed 

Google Scholar 
28.Dai, T. et al. Dynamics of coastal bacterial community average ribosomal RNA operon copy number reflect its response and sensitivity to ammonium and phosphate. Environ. Pollut. 260, 113971 (2020).CAS 
PubMed 

Google Scholar 
29.Zhou, J., Deng, Y., Luo, F., He, Z. & Yang, Y. Phylogenetic molecular ecological network of soil microbial communities in response to elevated CO2. Mbio 2, e00122–00111 (2011).PubMed 
PubMed Central 

Google Scholar 
30.Vellend, M. Conceptual Synthesis in Community Ecology. Q Rev. Biol. 85, 183–206 (2010).PubMed 

Google Scholar 
31.Frey, E. Evolutionary game theory: Theoretical concepts and applications to microbial communities. Phys. A 389, 4265–4298 (2010).MathSciNet 
CAS 
MATH 

Google Scholar 
32.Boeuf, D. et al. Biological composition and microbial dynamics of sinking particulate organic matter at abyssal depths in the oligotrophic open ocean. P Natl Acad. Sci. USA 116, 11824 (2019).CAS 

Google Scholar 
33.Thompson, L. R. et al. A communal catalogue reveals Earth’s multiscale microbial diversity. Nature 551, 457–463 (2017).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
34.Gralka, M., Szabo, R., Stocker, R. & Cordero, O. X. Trophic Interactions and the Drivers of Microbial Community Assembly. Curr. Biol. 30, R1176–R1188 (2020).CAS 
PubMed 

Google Scholar 
35.Gorter, F. A., Manhart, M. & Ackermann, M. Understanding the evolution of interspecies interactions in microbial communities. Philos. Trans. R. Soc. Lond. B Biol. Sci. 375, 20190256 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
36.Gandhi, S. R., Korolev, K. S. & Gore, J. Cooperation mitigates diversity loss in a spatially expanding microbial population. P Natl Acad. Sci. USA 116, 23582–23587 (2019).CAS 

Google Scholar 
37.Calatayud, J. et al. Positive associations among rare species and their persistence in ecological assemblages. Nat. Ecol. Evol. 4, 40–45 (2020).PubMed 

Google Scholar 
38.Furman, O. et al. Stochasticity constrained by deterministic effects of diet and age drive rumen microbiome assembly dynamics. Nat. Commun. 11, 1904 (2020).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
39.Tardy, V. et al. Stability of soil microbial structure and activity depends on microbial diversity. Env. Microbiol. Rep. 6, 173–183 (2014).CAS 

Google Scholar 
40.Coyte, K. Z., Schluter, J. & Foster, K. R. The ecology of the microbiome: networks, competition, and stability. Science 350, 663–666 (2015).CAS 
PubMed 
ADS 

Google Scholar 
41.Blanchet, F. G., Cazelles, K. & Gravel, D. Co-occurrence is not evidence of ecological interactions. Ecol. Lett. 23, 1050–1063 (2020).PubMed 

Google Scholar 
42.Blasche, S. et al. Metabolic cooperation and spatiotemporal niche partitioning in a kefir microbial community. Nat. Microbiol. 6, 196–208 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
43.Chatzinikolaou, E. et al. Spatio-temporal benthic biodiversity patterns and pollution pressure in three Mediterranean touristic ports. Sci. Total Environ. 624, 648–660 (2018).CAS 
PubMed 
ADS 

Google Scholar 
44.Filippini, G. et al. Sediment bacterial communities associated with environmental factors in Intermittently Closed and Open Lakes and Lagoons (ICOLLs). Sci. Total Environ. 693, 133462 (2019).CAS 
PubMed 
ADS 

Google Scholar 
45.Pesant, S. et al. Open science resources for the discovery and analysis of Tara Oceans data. Sci. Data 2, 150023 (2015).CAS 
PubMed 
PubMed Central 

Google Scholar 
46.Huse, S. M. et al. Exploring microbial diversity and taxonomy using SSU rRNA hypervariable tag sequencing. Plos Genet. 4, e1000255 (2008).PubMed 
PubMed Central 

Google Scholar 
47.Salas-González, I. et al. Coordination between microbiota and root endodermis supports plant mineral nutrient homeostasis. Science 371, eabd0695 (2021).PubMed 

Google Scholar 
48.Caporaso, J. G. et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. P Natl Acad. Sci. USA 108, 4516 (2011).CAS 
ADS 

Google Scholar 
49.Wear, E. K., Wilbanks, E. G., Nelson, C. E. & Carlson, C. A. Primer selection impacts specific population abundances but not community dynamics in a monthly time-series 16S rRNA gene amplicon analysis of coastal marine bacterioplankton. Environ. Microbiol. 20, 2709–2726 (2018).CAS 
PubMed 
PubMed Central 

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

Google Scholar 
51.Webb, C. O., Ackerly, D. D., McPeek, M. A. & Donoghue, M. J. Phylogenies and community ecology. Annu Rev. Ecol. Syst. 33, 475–505 (2002).
Google Scholar 
52.Wu, L. W. et al. Long-term successional dynamics of microbial association networks in anaerobic digestion processes. Water Res. 104, 1–10 (2016).PubMed 

Google Scholar 
53.Galand, P. E., Casamayor, E. O., Kirchman, D. L. & Lovejoy, C. Ecology of the rare microbial biosphere of the Arctic Ocean. P Natl Acad. Sci. USA 106, 22427–22432 (2009).CAS 
ADS 

Google Scholar 
54.Ju, F. & Zhang, T. Bacterial assembly and temporal dynamics in activated sludge of a full-scale municipal wastewater treatment plant. Isme J. 9, 683–695 (2015).CAS 
PubMed 
ADS 

Google Scholar  More

Current and past scenariosThe Arbutus unedo distribution presence database and eight environmental variables were used to predict the species potential distribution under past, current, and future climatic conditions. The current fitted model predicted the A. unedo suitable habitat with good performance, when evaluated using AUC (AUC = 0.92)31.The results suggest that the species distribution was mainly determined by the following attributes: (1) mean temperature of the driest quarter (BIO9), (2) annual mean temperature (BIO1), (3) slope and (4) temperature annual range (BIO7), summing up to 78.6% of the model’s total contributions (Table 1). In a recent study for the same species in Portugal, the authors considered that the variables precipitation seasonality (BIO15) and slope had a net significant influence on the species’ habitat suitability16. The variable slope was expected to influence the A. unedo distribution, since this tree is ecologically adapted to rocky slopes, which are edaphically unsuitable for woodland, where the vegetation is limited to shrub communities. The Arbutus unedo is particularly intolerant to shade, and this high light requirement is particularly needed for fruit production, restricting this species to settlement with open habitats21. This aspect was confirmed in a study about a woodland disturbance, where the A. unedo widespread initially, and, afterward, the species declined as canopy woodland re-developed35.The Arbutus unedo was adapted to a wide range of climatic conditions considering the response curves analysis (Fig. 1). According to Santiso24, this species developed a conservative strategy in the use of nutrients and water, when they are scarce. The Arbutus unedo plasticity and evolvability explained the current species persistence throughout its distribution range, which will be decisive in future response to climate change. The species distribution prediction for the present is in agreement with the species distribution according to Caudullo et al.38 (Fig. 6). Indeed, the A. unedo is widespread in Portugal, Galicia, and southwest of Spain, occupying the coastal belt from Tunisia to Morocco along the North of Africa, and from the Iberian Peninsula to Turkey across southern Europe, yet occurring in the northern Iberian Peninsula, western France, and south-western Ireland. The results showed that MaxEnt fitted model suitability prediction was high in those regions, which was further supported by the high observed AUC. Furthermore, it is important to stress the current potential distribution meaningful accuracy, and, therefore, the ability to extrapolate past and future prediction scenarios for A. unedo distribution in the Mediterranean Basin. Not surprisingly, the highest predicted suitability for the species occurs mostly in its native range, in regions where the mean January temperature is above 4ºC, a limit required for the species’ survival23. Unfortunately, this variable is only available in the WorldClim database for the present, but not for the future or past periods. Nevertheless, the species’ current spatial distribution was overlapped to the areas with mean January temperature higher than 4ºC, and a high overlap is reckoned, thus confirming this factor ecological importance17,23. The analysis of the species’ predicted suitability show that high-suitability areas [0.8–1.0] are more frequent in current conditions, compared to other scenarios, and thus confirming the suitability of climate current conditions for the species’ development (Fig. 2d).Figure 6Arbutus unedo L. distribution area and occurrences. Black dots represent the species’ occurrences (Supplementary Table S1), and the green area represents the species native distribution range for A. unedo downloaded from the data sets in38. Equal-area projection EPSG:3035. The map was made using the sf 1.0–3 and the terra 1.4–14 (https://rspatial.org/terra/) packages in R 4.1.032,33.Full size imageThe predicted habitat suitability for the LGM and MH scenarios differed, as expected, from the obtained for the current climatic conditions. The climatic conditions in the MH, were more suitable than during the LGM, due to climate warming, thus allowing the species expansion. According to the LGM projections, the species occurred in the Mediterranean islands (Sicily, Sardinia, and Balearic), North of Italy, south and eastern coast of Spain, Catalonia, North Africa, and spots in Portugal, Greece, and Turkey (Fig. 2a). These were suitable putative cryptic refugia areas that remained during the LGM. Keppel et al.39 defined refugia as locations to which species retreated during periods of adverse climate, and could potentially expand from, when environmental conditions turn out to be more suitable for the species. These locations/habitats were responsible for the species’ survival under changing environmental conditions for millennia. Identifying and characterizing climate refugia provides an important context for understanding the modern species distribution development, traits, and local adaptation40. The Iberian, Italian and Balkan peninsulas, which remained relatively ice-free during the ice ages were identified as the main glacial tree refugia areas in Europe41,42.The Betic mountains, in the Iberian Peninsula (IP), were sought to be a species’ refuge from the last ice age, nearby the sea (Valencia region, Spain), where fossil pollen pieces of evidence were found in the Canal de Navarrés peat deposit4, with 34 ka, and in the Siles lake (ca. 19 ka) located in the Segura mountains of southern Spain43 (Fig. 5, numbers 1 and 4, respectively). These mountains were ice-free during the Late Pleistocene, and the persistence of a mild climate could explain that other thermophilic species have prevailed in this particular region during the Full Glacial, proved by high genetic diversity levels44, and with the pollen spectra suggesting the region as a glacial refugium for temperate and Mediterranean trees9,13,43.The Iberian Peninsula border could have been a refuge for the species, due to the importance of the sea as a temperature regulator. Indeed, pollen evidence was found in the Basque mountains (northern IP, Fig. 5, number 3: 21 ka)45. Despite the information given by the MaxEnt modelling, the A. unedo reconstruction for the LGM was observed in other areas, since paleo-evidences were found in unexpected northern sites (Donatella Magri, personal communication). Certainly, fossil charcoal pieces of evidence were found in the central region of Portugal (Serra do Sicó), nearby the sea, supporting the presence of thermophilous taxa, including A. unedo, during the Full Glacial (24 ka, Fig. 5, number 2), and the presence of other Mediterranean taxa, such as olive tree and Pistacia lentiscus L., suggested that this may have been refuge zones for thermophilous plants46. Corroboration of the species expansion later during the warm and humid period of the Holocene (see Fig. 3b; green area) was confirmed by pollen remains (Fig. 5), particularly inland and northwardly, in the IP, and, also, in a molecular population genetics study made in Portugal, with northward haplotype migration47.Santiso et al.17 in a study about the A. unedo phylogeography, concluded that this species had two clades, separated during the late Pleistocene, before the LGM, suggesting that it may have coincided with the hardest glaciations recorded in the Quaternary. One clade occupied the Atlantic Iberia and, possibly, North Africa, while the other occurred in the western Mediterranean Basin in Spain. Besides, the A. unedo possibly persisted in the late Quaternary in the western Mediterranean, based on chloroplast DNA observations, and the results from the current study supported this interpretation, and by the fossil evidence (Supplementary Table 3). The same strong genetic differentiation between the western and eastern Mediterranean Basin was also found in olive lineages, and this pattern is congruent in other Mediterranean shrubs and tree species9,13,15. In another species (maritime pine), Bucci et al.44 confirmed the existence of a genetic divide between eastern and western lineages, also previously described by Burban and Petit48 based on mitochondrial DNA. Additionally, maximum haplotypic diversity for this species was found in south-eastern and central Spain, which, therefore, may be considered a biodiversity hotspot and a strong signal for refuge, as long-term populations tend to harbour more diversity than recently expanding ones41. In their study with A. unedo in Portugal using cpSSR, Ribeiro et al.47 found signals of two putative refugia in southern and central littoral in the country, also supported by macrofossil and pollen remains. Furthermore, according to Médail and Diadema41, several regions in the Western Mediterranean (large Mediterranean islands, North Africa and Catalonia), could have played a role in the case of A. unedo, and sought as refugia locations.The Arbutus unedo suitability maps obtained for past conditions support the claim that MH climatic conditions were more suitable for the species expansion than the LGM ones (Fig. 2a,b, and d). Consequently, spatial distribution changes analysis between past and current scenarios (Fig. 3) showed that A. unedo potential distribution expanded extensively from the LGM-current and the MH-current periods, indicating that more suitable areas were available for the species at present. Pollen records during the Early-Middle Holocene supported the Mediterranean presence of this species, in particular, nearby the sea (Fig. 5), which is under the MaxEnt predicted eastward dispersion (Fig. 3a, b). In another species bird-dispersed, the Myrtus communis L. (myrtle), a spread from west to east was verified, following genetic differentiation during the Pleistocene15, since the western region of the Mediterranean Basin had a milder climate, compared to the eastern one, during the LGM.The Arbutus unedo had, probably, a considerable ability to disperse, migrating over thousands of kilometres and even crossing sea stretches, allowing the species expansion during the MH, when the climatic conditions became more favourable17. In the future, A. unedo migration possibility will depend on climatic change pace and, also, on seed dispersal, particularly long-distance dispersal events49, since long-distance migration fitness will be useful in future change scenarios, allowing species to progress to newly suitable areas, as happened in the past. Nevertheless, from past evidence in Ireland, Britain, and across Europe, trees migrated faster than would be expected, during the warming period at the beginning of the Holocene, and that biotic and/or abiotic vectors must have been involved in dispersal50. This is the case with A. unedo, the seeds are dispersed by different type of birds that eat the soft fruit and, some of them are migratory species, as the European robin (Erithacus rubecula L.)51.Future scenariosThe impact of climate change on the A. unedo potential distribution was assessed under two representative concentration pathways (RCP 4.5 and RCP8.5) for the years 2050 and 2070. The results showed that climate change, under both moderate and high emission scenarios, will affect the species distribution range. A decrease in the predicted suitable areas was, generally, observed, since the climate becomes less favourable for the species in the future.According to the results, under the RCP 4.5 scenario, the potential distribution area will increase up to 2050, and decrease afterward according to the prediction for 2070 (Fig. 2d, e, f). Considering the RCP 8.5 scenario, the suitable area will exhibit a net decrease from the present up to 2050 and continue this trend until 2070. It is also expected that he medium–high suitable areas will gradually decrease from the present to the future, especially for scenario RCP 8.5. These results are in agreement with those obtained by other authors, confirming that Mediterranean species (e.g. Quercus sp.) distribution areas will be negatively affected by future climate change52,53.Several studies concluded that global warming will influence species distributions by causing expansions, contractions, or shifts in the species ranges2. As expected, the results from the current study showed that suitable areas will contract under future climate scenarios when compared to the current conditions, though suitable areas will emerge. These effects’ impact will depend on the climate change scenario severity. Despite contraction and expansion effects, the presence of the species will gradually decrease from the 2050s to the 2070s. Moreover, a species’ shift toward the North will be verified, because of suitable areas emergence, observed mainly in the RCP 8.5 scenario. Those areas will mostly emerge in the North of France, South of the United Kingdom, and Ireland, implying species’ latitudinal migration. The species’ northward displacement is consistent with climate change studies’ results obtained by several authors, including A. unedo16,20,54,55. Nevertheless, according to Gerassis et al.20, under climate change, the expected habitat disruption and fragmentation could lead to very adverse conditions for A. unedo survival in the future, which could undoubtedly conduct to a possible species’ presence decline in most of the current distribution area. Moreover, distribution models that predict climate‐induced range shifts do not account for spatial dispersal variation56, but adaptive dispersal evolution always reduced neutral genetic diversity across the species’ range. This means that the species’ genetic pool might be erased, depending on the climate change velocity, amongst other conditions, like landscape fragmentation and competition with other species/crops16,57. Additionally, the species ability to migrate mainly through seeds, dispersed by migratory frugivorous birds, bird abundance, and the velocity of climate change, are key issues for future species survival23,58.These results suggested that future changes in environmental conditions may lead to suitable habitat loss in areas where the species had persisted and with a possible range shift towards the North. These findings also revealed that with continuous future climate warmth, the current potential distribution A. unedo areas will become unsuitable or contract, leading to significant changes in the species’ current distribution pattern and putative presence loss. The possibility of species’ migration will ultimately depend on its capability to keep pace with the changing conditions and the velocity to adapt to environmental changes, such as those presented by habitat and climate modifications56.The Mediterranean Basin is one of the most vulnerable climate change hotspots in the world, thus understanding how future climate changes will disturb Mediterranean plant species distribution will be key for tree management planning and conservation design. Nevertheless, further investigation is needed for species well adapted in this region to assess the impacts of climate change in their current and future potential distributions. Those studies including past climate impact on species distribution should be complemented with phylogeographic methods and paleoclimate reconstructions to locate refuges. Other species distribution models, besides MaxEnt, could be tested, although the lack of absence records data limits considerably the modelling approach. Additionally, the human influence magnitude on the predicted ecological niche should be further studied in detail, including urbanization, industrialization and putative tourism pressure, particularly in coastal areas. More

Forest fires are disasters that cause extensive damage to the entire world in economic, ecological, and environmental ways. These fires can be caused by natural reasons, such as high temperatures that can create spontaneous combustion of dry fuel such as sawdust, leaves, lightning, etc., or by human activities, such as unextinguished campfires, arson, inappropriately burned debris, etc1. According to research, 90% of the world’s forest fire incidents have occurred as a result of the abovementioned human carelessness1. The increase in carbon dioxide levels in the atmosphere due to forest fires contributes to the greenhouse effect and climate change. Additionally, ash destroys much of the nutrients in the soil and can cause erosion, which may result in floods and landslides.
At earlier times, forest fires were detected using watchtowers, which were not efficient because they were based on human observations. In recent history and even the present day, several forest fire detection methods have been implemented, such as watchtowers, satellite image processing methods, optical sensors, and digital camera-based methods2, although there are many drawbacks, such as inefficiency, power consumption, latency, accuracy and implementation costs. To address these drawbacks, a forest fire detection system using wireless sensor networks is proposed in this paper.Wireless sensor networks (WSNs) are self-configured and infrastructure-free wireless networks that help monitor physical or environmental conditions and pass these data through the network to a designated location or sink where the data can be observed and analyzed3. Efficiency and low power consumption are the major advantages of a WSN. In the proposed detection system, wireless sensor nodes are deployed according to cellular architecture to cover the entire area with sensors to monitor temperature, relative humidity, light intensity level, and carbon monoxide (CO) level using a microcontroller, transceiver module, and power components. The power supply to the sensor node is provided using batteries as the primary power supply, and solar panels are used as the secondary power supply. These sensor nodes are specially designed with a spherical shape to withstand damage caused by environmental conditions as well as animals.The sensor readings for each parameter are checked with a preset threshold ratio and a ratio that is calculated continuously in the node in real time, and only the ratios that exceed the preset ratio are sent from the sensor node to the base station for further analytical processing. The network utilized for this transmission is in the architecture of tree topology considering facts such as low power consumption, reduced latency, less complexity, etc. Cluster heads are used in this network to gather data from several sensor nodes and pass them on to the base station or the gateway node. The gateway node is an interface that connects the network with the secondary analysis process.For the analysis process, a machine learning regression model was used along with threshold ratio analysis to enhance detection accuracy. For the training and testing process of the model, data were collected during the fire and no fire situations in different areas and under different climatic conditions. During the data collection process, 7000 data samples were collected, where a data sample included temperature, relative humidity, light intensity level, and CO level at a particular time. Eighty percent of the collected data were randomly used as training data for the model, and the remaining 20% were used as test data.If the outcome of the machine learning model indicates a fire in a specific area, a text message will be sent to the mobile phone numbers of the authorized officers in responsible units. As this process is designed with a minimum delay, the fire can be detected within the initial stage, and the responsible parties can take necessary actions in a shorter period, which will minimize the damage.Related workForest fire detection has been a focus of many researchers for the last decade because of increased forest fire case reports from all over the world due to severe damage to society and the environment. Many methods have been proposed to detect forest fires, such as camera-based systems, WSN-based systems, and machine learning application-based systems, with both positive and negative aspects and performance figures of detection. Due to the higher probability of accurate and early detection due to the use of multiple sensor sources and deployment of sensor nodes in areas not visible to satellites, wireless sensor networks have a more positive outlook, and they have become the more applicable technology in many fields4.Many researchers have focused on environmental parameters, such as air temperature, relative humidity, barometric pressure, sound, light intensity, soil moisture, and wind speed and direction, along with gases, such as CO, CO2, methane, H2, and hydrocarbons apart from smoke, to detect forest fire conditions by considering the variations in these parameters during a fire5,6, and sensors have been selected according to the range, sensitivity, power consumption, and cost7,8.As supplying power to a sensor node is a challenging task in forested areas, utilizing only battery options is difficult because they do not last long, and distributing power using a wire would require a higher cost to deploy in a large forest. Therefore, many researchers have proposed solar-powered systems as secondary power sources along with rechargeable batteries as the main power source4,6, while some researchers have proposed solar batteries because they last longer9. To reduce the power consumption of sensor nodes, techniques such as keeping selected components active while others are deactivated have been proposed10,11,12.Most WSN-based detection systems are centered around a base station due to the memory and processing limitations of the nodes. Important and partially processed data are transferred to the base station through wireless media for processing and enabling relevant actions, while the base station also acts as the gateway between the sensor nodes and the system user4,9.When constructing a WSN, communicating data among the relevant entities is the main objective, and star topology and mesh topology-based networks have been proposed in many papers because of the different attributes in their systems. A mesh topology was chosen over a star topology because of its ability to self-organize, self-configure, and automatically establish among nodes in a network13. As a smaller number of nodes involved for transmission results in minimum energy consumption, concepts based on cluster heads have been used14. To minimize the loss of energy and data packages during transfer, a cluster-tree network topology structure was proposed15. Considering the sensing range of a node, fault tolerance, and energy consumption, a paper has proposed applying the on-demand k-coverage technique that provides event detection using static nodes with variable sensing ranges. This technique utilizes the maximum detection performance with the minimum power consumption for an event16. A survey on rare event detection has mentioned many event detection strategies that deliver maximized detection capability, minimized detection delay and low energy consumption, such as duty cycle, component deactivation, overpopulation/node redundancy, collaboration, and energy harvesting17.To reduce deployment cost and power consumption, a paper proposes a novel localization scheme that divides the whole forest area into different grids and allocates them to respective zones with another 8 neighboring grids. One centroid node from those grids, which is called the initiator node, predicts whether the zone is highly active (HA), medium active (MA), and low active (LA). Here, HA zones send data continuously to the base node through the interior node, MA zones send data periodically, and LA zones do not transfer data in the status that manages power consumption effectively18. Another author proposed obtaining data from the sensors every 2 min if there is the potential of a forest fire or obtaining data every 15 min otherwise to reduce the energy wastage19.To place sensor nodes in the most effective configuration to detect fire conditions, a sensor node was proposed at three different heights to perform different parameter measurements, while some authors have suggested covering sensor nodes to avoid direct sunlight exposure and minimize the false alarm rate4,10. As the network connectivity of service providers in forest areas is not robust, communication techniques that use dedicated network paths such as LoRa, ZigBee, and XBee have been used as the communication infrastructure. When considering attributes such as transmission range, high security, low power consumption (LPWAN protocol), and other relevant configurations, most papers have suggested using the LoRa module for transmission13,19.Most papers have suggested having threshold value-based fire detection on a sensor node, and if the exceeded threshold has remained the same, then a sink determines the location and will send an alarm to the fire department11. Because of the environmental parameter variations according to the place and time, threshold values are configured by the user considering geographic situations, climatic changes, seasonal changes, etc. after sensors obtain the data from the surrounding10.A fusion information process was proposed, where information from multiple sources is considered in making the final decision, which is better than using those sources individually, and two algorithms based on the threshold ratio method and Dempster-Shafer theory were used4. To enhance the detection accuracy, machine learning applications have been proposed in many papers based on different machine learning approaches, such as support vector machine (SVM) classification14 and regression techniques, such as logistic regression. However, applying machine learning techniques to fire detection systems has many limitations, such as the limited amount of energy, the energy required for data processing, the short range of communication and limited computations, the complexity of ML algorithms when executing on sensor nodes, and the difficulty of being distributed on every sensor node20,21. More

1.Thompson, R. S., Parish, J. C., Maidment, S. C. R. & Barrett, P. M. Phylogeny of the ankylosaurian dinosaurs (Ornithischia: Thyreophora). J. Syst. Palaeontol. 10(2), 301–312 (2012).Article 

Google Scholar 
2.Arbour, V. M. & Currie, P. J. Ankylosaurid dinosaur tail clubs evolved through stepwise acquisition of key features. J. Anat. 227, 514–523. https://doi.org/10.1111/joa.12363 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
3.Brown, C. M. An exceptionally preserved armored dinosaur reveals the morphology and allometry of osteoderms and their horny epidermal coverings. PeerJ 5, e4066. https://doi.org/10.7717/peerj.4066 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
4.Butler, R. J. & Barrett, P. M. Palaeoenvironmental controls on the distribution of Cretaceous herbivorous dinosaurs. Sci. Nat. 95, 1027–1032 (2008).CAS 
Article 

Google Scholar 
5.Brown, C. M. et al. Dietary palaeoecology of an Early Cretaceous armoured dinosaur (Ornithischia; Nodosauridae) based on floral analysis of stomach contents. R. Soc. Open Sci. 7, 200305. https://doi.org/10.1098/rsos.200305 (2020).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
6.Bourke, J. M., Porter, W. R. & Witmer, L. M. Convoluted nasal passages function as efficient heat exchangers in ankylosaurs (Dinosauria: Ornithischia: Thyreophora). PLoS ONE 13, e0207381. https://doi.org/10.1371/journal.pone.0207381 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
7.Mallon, J. C. & Anderson, J. S. The functional and palaeoecological implications of tooth morphology and wear for the megaherbivorous dinosaurs from the dinosaur park formation (Upper Campanian) of Alberta, Canada. PLoS ONE 9(6), e98605. https://doi.org/10.1371/journal.pone.0098605 (2014).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
8.Ősi, A., Prondvai, E., Mallon, J. & Bodor, E. R. Diversity and convergences in the evolution of feeding adaptations in ankylosaurs (Dinosauria: Ornithischia). Hist. Biol. 29(4), 539–570. https://doi.org/10.1080/08912963.2016.1208194 (2017).Article 

Google Scholar 
9.Hayashi, S., Carpenter, K., Scheyer, T. M., Watanabe, M. & Suzuki, D. Function and evolution of ankylosaur dermal armor. Acta Palaeontol. Pol. 55(2), 213–228 (2010).Article 

Google Scholar 
10.Norman, D. B. Scelidosaurus harrisonii from the Early Jurassic of Dorset, England: Cranial anatomy. Zool. J. Linn. Soc. 188(1), 1–81. https://doi.org/10.1093/zoolinnean/zlz074 (2019).Article 

Google Scholar 
11.Galton, P. M. Endocranial casts of the plated dinosaur Stegosaurus (Upper Jurassic, Western USA): A complete undistorted cast and the original specimens of Othniel Charles Marsh. In The Armored Dinosaurs (ed. Carpenter, K.) 103–129 (Indiana University Press, 2001).
Google Scholar 
12.Galton, P. M. Skull bones and endocranial casts of stegosaurian dinosaur Kentrosaurus Hennig, 1915 from Upper Jurassic of Tanzania, East Africa. Geol. Palaeontol. 22, 123–143 (1988).
Google Scholar 
13.Kuzmin, I. et al. The braincase of Bissektipelta archibaldi—new insights into endocranial osteology, vasculature, and paleoneurobiology of ankylosaurian dinosaurs. Biol. Commun. 65(2), 85–156. https://doi.org/10.21638/spbu03.2020.201 (2020).Article 

Google Scholar 
14.Leahey, L. G., Molnar, R. E., Carpenter, K., Witmer, L. M. & Salisbury, S. W. Cranial osteology of the ankylosaurian dinosaur formerly known as Minmi sp. (Ornithischia: Thyreophora) from the Lower Cretaceous Allaru Mudstone of Richmond, Queensland, Australia. PeerJ 3, e1475. https://doi.org/10.7717/peerj.1475 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
15.Paulina-Carabajal, A., Lee, Y. N., Kobayashi, Y., Lee, H. J. & Currie, P. J. Neuroanatomy of the ankylosaurid dinosaurs Tarchia teresae and Talarurus plicatospineus from the Upper Cretaceous of Mongolia, with comments on endocranial variability among ankylosaurs. Palaeogeogr. Palaeoclimatol. Palaeoecol. 494, 135–146. https://doi.org/10.1016/j.palaeo.2017.11.030 (2018).Article 

Google Scholar 
16.Pereda-Suberbiola, X. & Galton, P. M. On the taxonomic status of the dinosaur Struthiosaurus austriacus Bunzel from the Late Cretaceous of Austria. C. R. Acad. Sci. Paris II 315, 1275–1280 (1992).
Google Scholar 
17.Pereda-Suberbiola, X. & Galton, P. M. Revision of the cranial features of the dinosaur Struthiosaurus austriacus Bunzel (Ornithischia: Ankylosauria) from the Late Cretaceous of Europe. Neues Jahrb. Geol. Palaontol. Abh. 191, 173–200 (1994).
Google Scholar 
18.Pereda-Suberbiola, X. & Galton, P. M. Reappraisal of the nodosaurid ankylosaur Struthiosaurus austriacus Bunzel from the Upper Cretaceous Gosau Beds of Austria. In The Armored Dinosaurs (ed. Carpenter, K.) 173–210 (Indiana University Press, 2001).
Google Scholar 
19.Pereda-Suberbiola, X. A revised census of European Late Cretaceous nodosaurids (Ornithischia: Ankylosauria): Last occurrence and possible extinction scenarios. Terra Nova 4, 641–648 (1992).ADS 
Article 

Google Scholar 
20.Pereda-Suberbiola, X. Ankylosaurian dinosaur remains from the Upper Cretaceous of Laño (Iberian Peninsula). Est. Mus. Cienc. Nat. de Álava 14(Número especial 1), 273–288 (1999).
Google Scholar 
21.Pereda-Suberbiola, X., Astibia, H. & Buffetaut, E. New remains of the armoured dinosaur Struthiosaurus from the Late Cretaceous of the Iberian peninsula (Lafio locality, Basque-Cantabric Basin). Bull. Soc. Géol. Fr. 166, 207–211 (1995).
Google Scholar 
22.Garcia, G. & Pereda-Suberbiola, X. A new species of Struthiosaurus (Dinosauria: Ankylosauria) from the upper cretaceous of Villeveyrac (Southern France). J. Vertbr. Paleontol 23(1), 156–165 (2003).Article 

Google Scholar 
23.Nopcsa, F. Dinosaurierreste aus Siebenbürgen V. Geologica Hungarica, Series Palaeontologica 4, 1–76 (1929).
Google Scholar 
24.Codrea, V. et al. More than just Nopcsa’s Transylvanian dinosaurs: A look outside the Haţeg Basin. Palaeogeogr. Palaeoclimatol. Palaeoecol. 293, 391–405. https://doi.org/10.1016/j.palaeo.2009.10.027 (2010).Article 

Google Scholar 
25.Ősi, A. & Prondvai, E. Sympatry of two ankylosaurs (Hungarosaurus and cf. Struthiosaurus) in the Santonian of Hungary. Cretac. Res. 44, 30–38. https://doi.org/10.1016/j.cretres.2013.03.006 (2013).Article 

Google Scholar 
26.Bunzel, E. Die Reptilfauna der Gosaformation in der Neuen Welt bei Wiener-Neustadt. Abhandlungen der Kaiserlich-Königlichen Geologischen Reichsanstalt. 5, 1–18 (1871).
Google Scholar 
27.Franzosa, J. & Rowe, T. Cranial endocast of the Cretaceous theropod dinosaur Acrocanthosaurus atokensis. J. Vertebr. Paleontol. 25, 859–864 (2005).Article 

Google Scholar 
28.Witmer, L. M. & Ridgely, R. C. New insights into the brain, braincase, and ear region of tyrannosaurs (Dinosauria, Theropoda), with implications for sensory organization and behavior. Anat. Rec. 292(9), 1266–1296. https://doi.org/10.1002/ar.20983 (2009).Article 

Google Scholar 
29.Schade, M., Rauhut, O. W. M. & Evers, S. W. Neuroanatomy of the spinosaurid Irritator challengeri (Dinosauria: Theropoda) indicates potential adaptations for piscivory. Sci. Rep. 10(9259), 1613–1616. https://doi.org/10.1038/s41598-020-66261 (2020).Article 

Google Scholar 
30.Witmer, L. M., Ridgely, R. C., Dufeau, D. L. & Semones, M. C. Using CT to peer into the past: 3D visualization of the brain and ear regions of birds, crocodiles, and nonavian dinosaurs. In Anatomical imaging: towards a new morphology (eds Endo, H. & Frey, R.) 67–87 (Springer, 2008).Chapter 

Google Scholar 
31.Evers, S. W. et al. Neurovascular anatomy of the protostegid turtle Rhinochelys pulchriceps and comparisons of membranous and endosseous labyrinth shape in an extant turtle. Zool. J. Linn. Soc. 187, 800–828. https://doi.org/10.1093/zoolinnean/zlz063 (2019).Article 

Google Scholar 
32.Ősi, A., Pereda Suberbiola, X. & Földes, T. Partial skull and endocranial cast of the ankylosaurian dinosaur Hungarosaurus from the Late Cretaceous of Hungary: Implications for locomotion. Palaeontol. Electro. 17(1), 18p (2014).
Google Scholar 
33.Sampson, S. D. & Witmer, L. M. Craniofacial anatomy of Majungasaurus crenatissimus (Theropoda: Abelisauridae) from the Late Cretaceous of Madagascar. J. Vertebr. Paleontol. 27, 32–104. https://doi.org/10.1671/0272-4634(2007)27[32:CAOMCT]2.0.CO;2 (2007).Article 

Google Scholar 
34.Walsh, S. A., Barrett, P. M., Milner, A. C., Manley, G. & Witmer, L. M. Inner ear anatomy is a proxy for deducing auditory capability and behaviour in reptiles and birds. Proc. R. Soc. B 276, 1355–1360 (2009).Article 

Google Scholar 
35.Porter, W. R. & Witmer, L. M. Vascular patterns in the heads of dinosaurs: Evidence for blood vessels, sites of thermal exchange, and their role in physiological thermoregulatory strategies. Anat. Rec. 303, 1075–1103. https://doi.org/10.1002/ar.24234 (2019).Article 

Google Scholar 
36.Benson, R. B. J., Starmer-Jones, E., Close, R. A. & Walsh, S. A. Comparative analysis of vestibular ecomorphology in birds. J. Anat. 231, 990–1018. https://doi.org/10.1111/joa.12726 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
37.Bronzati, M. et al. Deep evolutionary diversification of semicircular canals in archosaurs. Curr. Biol. 31(12), 2520-2529.e6. https://doi.org/10.1016/j.cub.2021.03.086 (2021).CAS 
Article 
PubMed 

Google Scholar 
38.Hanson, M., Hoffman, E. A., Norell, M. A. & Bhullar, B. S. The early origin of a birdlike inner ear and the evolution of dinosaurian movement and vocalization. Science 372(6542), 601–609. https://doi.org/10.1126/science.abb4305 (2021).ADS 
CAS 
Article 
PubMed 

Google Scholar 
39.Miyashita, T., Arbour, V. M., Witmer, L. M. & Currie, P. J. The internal cranial morphology of an armoured dinosaur Euoplocephalus corroborated by X-ray computed tomographic reconstruction. J. Anat. 219, 661–675. https://doi.org/10.1111/j.1469-7580.2011.01427.x (2011).Article 
PubMed 
PubMed Central 

Google Scholar 
40.Paulina-Carabajal, A., Lee, Y. N. & Jacobs, L. L. Endocranial morphology of the primitive nodosaurid dinosaur Pawpawsaurus campbelli from the Early Cretaceous of North America. PLoS ONE 11, e0150845. https://doi.org/10.1371/journal.pone.0150845 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
41.Spoor, F. & Thewissen, J. G. M. Comparative and functional anatomy of balance in aquatic mammals. In Sensory Evolution on the Threshold, Adaptations in Secondarily Aquatic Vertebrates (eds Thewissen, J. G. M. & Nummela, S.) 65–81 (University of California Press, 2008).
Google Scholar 
42.Sobral, G. & Müller, J. Archosaurs and their kin: The ruling reptiles. In Evolution of the Vertebrate Ear (eds Clack, J. A. et al.) 285–326 (Springer, 2016).Chapter 

Google Scholar 
43.Arbour, V. M. Estimating impact forces of tail club strikes by ankylosaurid dinosaurs. PLoS ONE 4(8), e6738. https://doi.org/10.1371/journal.pone.0006738 (2009).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
44.Park, J. Y. et al. A new ankylosaurid skeleton from the Upper Cretaceous Baruungoyot Formation of Mongolia: Its implications for ankylosaurid postcranial evolution. Sci. Rep. 11, 4101. https://doi.org/10.1038/s41598-021-83568-4 (2021).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
45.Witmer, L. M., Chatterjee, S., Franzosa, J. & Rowe, T. Neuroanatomy of flying reptiles and implications for flight, posture and behaviour. Nature 425, 950–953. https://doi.org/10.1038/nature02048 (2003).ADS 
CAS 
Article 
PubMed 

Google Scholar 
46.Walsh, S. A. et al. Avian cerebellar floccular fossa size is not a proxy for flying ability in birds. PLoS ONE 8(6), e67176. https://doi.org/10.1371/journal.pone.0067176 (2013).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
47.Ferreira-Cardoso, S. et al. Floccular fossa size is not a reliable proxy of ecology and behaviour in vertebrates. Sci. Rep. 7(1), 2017. https://doi.org/10.1038/s41598-017-01981-0 (2005).CAS 
Article 

Google Scholar 
48.Ezcurra, M. D. et al. Enigmatic dinosaur precursors bridge the gap to the origin of Pterosauria. Nature 588(7838), 445–449. https://doi.org/10.1038/s41586-020-3011-4 (2020) (Epub 2020 Dec 9).ADS 
CAS 
Article 
PubMed 

Google Scholar 
49.Carpenter, K., Sanders, F., McWhinney, L. A. & Wood, L. Evidence for predator–prey relationships: examples for Allosaurus and Stegosaurus. In The Carnivorous Dinosaurs (ed. Carpenter, K.) 325–350 (Indiana University Press, 2005).
Google Scholar 
50.Mallison, H. CAD assessment of the posture and range of motion of Kentrosaurus aethiopicus HENNIG 1915. Swiss. J. Geosci. 103, 211–233. https://doi.org/10.1007/s00015-010-0024-2 (2010).Article 

Google Scholar 
51.Arbour, V. M., Lech-Hernes, N. L., Guldberg, T. E., Hurum, J. H. & Currie, P. J. An ankylosaurid dinosaur from Mongolia with in situ armour and keratinous scale impressions. Acta Palaeontol. Pol. 58(1), 55–64. https://doi.org/10.4202/app.2011.0081 (2013).Article 

Google Scholar 
52.Scheyer, T. M. & Sander, P. M. Histology of ankylosaur osteoderms: Implications for systematics and function. J. Vertebr. Paleontol. 24, 874–893 (2004).Article 

Google Scholar 
53.Carpenter, K. Ankylosauria. In The Complete Dinosaur 2nd edn (eds Brett-Surman, M. K. et al.) 505–526 (Indiana University Press, 2012).
Google Scholar 
54.Jerison, H. J. Evolution of the brain and intelligence 482 (Academic Press, 1973).
Google Scholar 
55.Marugan-Lobon, J., Chiappe, L. M. & Farke, A. A. The variability of inner ear orientation in saurischian dinosaurs: Testing the use of semicircular canals as a reference system for comparative anatomy. PeerJ 1, e124 (2013).Article 

Google Scholar 
56.Benoit, J. et al. A test of the lateral semicircular canal correlation to head posture, diet and other biological traits in “ungulate” mammals. Sci. Rep. 10, 19602. https://doi.org/10.1038/s41598-020-76757-0 (2020).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
57.Coombs, W. P. Jr. Forelimb muscles of the Ankylosauria (Reptilia, Ornithischia). J. Paleontol. 52, 642–657 (1978).
Google Scholar 
58.Maidment, S. C. R. & Barrett, P. M. Osteological correlates for quadrupedality in ornithischian dinosaurs. Acta Palaeontol. Pol. 59(1), 53–70 (2014).
Google Scholar  More

Symbiont diversity and distributionThe present study provides the first evidence of symbiosis in R. alisae, a species of nudibranchs. This is the most multiple symbiosis that have ever been recorded for marine invertebrates. While many organisms establish an exclusively one-on-one relationship with a single microbial species or microbes belonging to the same functional group5,12, there are also organisms that harbor multiple microbial species, in which symbiont–symbiont and host–symbiont interactions occur. Six phylotypes of chemoautotrophic bacteria were reported for mussel Idas sp. from a cold seep area11 and five extracellular symbionts for the gutless oligochaete worm Olavius algarvensis34. However, in these cases, symbioses involving bacteria and marine invertebrates are either endosymbiotic microbes co-occurring inside the host bacteriocytes5,11 or ectosymbiotic microbes associated with the external surfaces of the animals3,4,9,15,34, with the exception of scaly-foot snail from hydrothermal vents having partnerships simultaneously with epi- and endosymbiontic microbes35.Bacterial symbionts in R. alisae have appeared to be more diverse than was previously known for marine invertebrates. It is evident that the detected symbiont phylotypes differ greatly from all other known symbionts found in marine invertebrates. Labrenzia (Rodobacteriales) and Maritalea (Rhizobiales) have not been recorded as forming symbiotic associations with invertebrates or plants so far, although other members of the families Rodobacteriales and Rhizobiales are well known symbionts14. Strains of Bradyrhizobium, Burkholderia, Achromobacter, and Stenotrophomonas are reported as symbionts of plants, interacting with a vast majority of nodulating legume species and efficient in biological nitrogen fixation36. This may be important when considering the nature of these symbionts in the nudibranch. Symbioses between cyanobacteria and marine organisms are commonly found among marine plants, fungi, sponges, ascidians, corals, and protists37,38. Synechococcus, identified as dominant symbiont clones of R. alisae (Table S2), is a unicellular cyanobacterium common in the marine environment, providing a range of beneficial functions including photosynthesis, nitrogen fixation, UV protection, and production of defensive toxins8,9,37. Symbiotic interactions between actinobacteria and their host have been observed in insects, human, animals, and plants, where the bacteria provide the host with protection against pathogens and produce essential nutrients39. However, none of the members of the clade Actinobacteria recorded in R. alisae are known to live symbiotically.Arrangement of symbiotic associationDespite the high diversity of bacteria, they are well organized in the host. Dense clusters of rod-shaped bacteria, Labrenzia, Maritalea, Bradyrhizobium, Burcholderia, Achromobacter, and Stenotrophomonas, were found within host-derived vacuoles, referred to as bacteriocytes, inside epithelial cells of R. alisae (Fig. 3). Although such arrangement differs from that typical of bacteriocytes, which are usually considered as specialized cells of the hosts for harboring bacteria, it resembles that reported for scaly-food snail from hydrothermal vents, which harbor symbionts in the esophageal gland35. Bacteriocytes in the gastropod Lurifax vitreus found near hydrothermal vents also constitute a portion of the mantle epithelium; they have large vacuoles containing many live and dividing bacteria40. Each bacteriocyte was densely packed with certain symbionts, and the bacteriocytes were randomly distributed within the epithelium cells. A distinctly regular distribution pattern was observed in the gill epithelium of the mussel Bathymodiolus sp.: the thiotrophic symbionts occupy the apical region, and the methanotrophic symbionts are more abundant in the basal region of bacteriocytes4. In the mussel Idas sp., however, there is no spatial pattern of the six distinct bacterial phylotypes, and the symbionts are mixed within bacteriocytes11.Synechococcus dominated the cytoplasm of intestinal epithelium and, more rarely, epidermis cells, mainly as specialized cell type referred to as nitrogen-fixing heterocysts. They are visually similar to cyanobacteria from corals and sponges8,37.The phylogenetic diversity and the spatial organization of the symbiotic community in R. alisae were determined by the 16S rRNA analysis, which was consistent with the results of FISH and TEM. Unlike most symbioses of marine invertebrates when bacteria house specialized host cells5,11 or cover epidermis7,15, symbiotic association of R. alisae exhibited spatial partitioning between symbionts, which were unevenly distributed between the tissues (Table S2). It has been established that different members of the microbial community can complement each other in acquisition of various restrictive nutrients, confirming the importance of the functional diversity of symbionts41. Thus, Stenotrophomonas rhizophila and Bradyrhizobium build a beneficial association in the rhizosphere and can act synergistically on promoting growth and nutrient uptake of soybean36. Cyanobacteria can interact synergistically with beneficial members from the endophytic microbiome of rice seedlings42. The location of bacterium in the organism of R. alisae may, in fact, depend on the specific metabolic and ecological roles that the symbionts play, and also on the interaction with bacterium belonging to different physiological groups.Nature of symbiosisSymbiotic associations between microbes and invertebrates are acquired mainly in a nutrient-depleted environment where symbionts usually provide their hosts with essential nutrients and high-energy compounds1. In contrast to known symbioses between microbes and gutless invertebrates, which obtain nutrients exclusively from the bacteria, R. alisae, like most nudibranch species, is a sponge-eating predator. However, due to the lack of adipose tissue, sponges are distinguished by a low lipid content (0.4 to 1.5% of wet weight)43 and also by specific proteinaceous spongin fibers and chitin, a polysaccharide similar to cellulose that can be indigestible for some predators, which together indicate their low nutritional value. Furthermore, R. alisae feeds exclusively on the sponge O. pennata; therefore, in habitats with low prey availability, this nudibranch has to survive starvation while searching for sponge assemblages. We suppose that symbiotic bacteria of R. alisae contribute to the utilization of low-quality food, similarly to symbiotic bacteria from the genera Rhodobacter, Burkholderia, and Aeromonas associated with the detritivorous isopod Asellus aquaticus44.A fatty acid analysis, as a useful approach to clarifying the nature of symbiosis5,20,32, has confirmed the trophic interaction between symbionts and the nudibranch host (Table S2). Among the fatty acids of symbiotic bacteria in R. alisae, OBFA are a major acyl constituent of membranes in Stenotrophomonas45 and also in Actinobacteria, Arthrobacter, Iamia, Ilumatobacter, and Kocuria46. Cis-vaccenic acid is a major component of Maritalea30. Omega-cyclohexyl tridecanoic acid (cyclo19:0) is specific to Bradyrhizobium47, Burkholderia, and Achromobacter48. Linoleic acid is produced by cyanobacteria including marine species of Synecoccocus33; in nudibranch, it obviously serves as a precursor in the synthesis of arachidonic acid (20:4n-6), thus, providing additional evidence for the transfer of fatty acids from symbionts to the host. Mollusks are capable of converting linoleic acid to arachidonic acid, since they have enzymes required for its synthesis21. The presence of these bacteria-specific markers and the abundance of arachidonic acid confirm the metabolic role of symbionts in the nudibranch host.Among nutrients, biologically available nitrogen can be considered a restrictive nutrient for the sponge-eating R. alisae, which can be acquired with the help of nitrogen-fixing symbionts, also referred to as diazotrophs. R. alisae harbors Bradyrhizobium, Burkholderia, Achromobacter, and Stenotrophomonas that are efficient in biological nitrogen fixation previously found to be associated with nodulating legume species36. Symbiotic nitrogen fixers are known to be associated with a variety of marine invertebrates such as wood-boring bivalves, corals, sponges, sea urchins, tunicates, and polychaetes7,8,37. Moreover, the protection of the enzyme nitrogenase that catalyzes N2 fixation against oxygen is an important physiological requirement in bacteria such as symbiotic Bradyrhizobium, Burkholderia, Achromobacter, and Stenotrophomonas that are located in bacteriocytes and provide this protection. Synechococcus is known as a nitrogen-fixer37,49. It performs N2 fixation in heterocysts where nitrogenase is restricted under oxic conditions. Indeed, heterocysts of Synechococcus are abundant in the intestine cells of R. alisae (Fig. 5B–D).Nitrate assimilation is one of the major processes of nitrogen acquisition by many heterotrophic bacteria and cyanobacteria50,51. Symbionts of R. alisae can play an important role in the process of nitrate utilization through denitrification, dissimilatory nitrate reduction, and assimilatory nitrate reduction as a nitrogen source and synthesize it into organic nitrogen. The nitrate reducers, Labrenzia52, Stenotrophomonas53, Maritalea30, and Rhodobacteraceae29 are widely represented in R. alisae. Synechococcus also utilizes nitrate, nitrite, or ammonium for growth50. Thus, symbiotic bacteria may play a significant role in the N-budget of the nudibranch mollusk.The symbiotic bacteria of R. alisae, including Bradyrhizobium, Maritalea, Labrenzia, Burkholderia, Achromobacter, Stenotrophomonas, Arthrobacter, Iamia, Ilumatobacter, and Kocuria, are known as carboxydotrophic or carbon monoxide (CO) oxidizers54,55. Despite the toxicity of CO for multicellular organisms, numerous aerobic and anaerobic microorganisms can use CO as a source of energy and/or carbon for cell growth56. The marine worm Olavius algarvensis establishes symbiosis with chemosynthetic bacteria using CO, a substrate previously not known to play a role in symbiotic associations with marine invertebrates, as an energy source57. We do not rule out that the R. alisae symbionts also might exploit CO as carbon and energy source. Despite this, assumption may seem impossible taking in account the CO toxicity, but, since many invertebrates (mollusks, tube worm, etc.) use toxic sulfate, thiosulfate, and methane as an energy source1,15, this hypothesis is worth to be addressed.An important component of skeleton in marine sponges of the family Microcionidae, including O. pennata, is the structural polysaccharide chitin58. Some bacteria are capable of hydrolyzing chitin via the activity of chitinolytic enzymes and can utilize chitin as a source of carbon, nitrogen, and/or energy59. Chitinase activity was documented for strains of Labrenzia60, Burkholderia61, Arthrobacter62, Achromobacter63, Stenotrophomonas64, Alcaligenes65, and actinobacteria59 associated with R. alisae. Thus, these bacteria can work synergistically to digest chitin and spongin, contributing to feeding success of the host nudibranch which depends solely on low-quality, nitrogen- and carbon-deficient food available.Furthermore, direct evidence has confirmed that many bioactive compounds from invertebrates are produced by symbiotic microorganisms66,67. Many biologically active compounds including toxic and deterrent secretions have been identified in nudibranchs of the family Discodorididae68. Symbiotic bacteria may exhibit toxic activity to provide the host nudibranch with chemical defense against predators and environment. Bacteria, especially actinobacteria, living in a symbiotic relationship with R. alisae may help the host in defense, since nudibranch lack a shell, and secondary metabolites of bacteria can provide them with chemical defense against predators and environment, as has been reported for some marine invertebrates2,9,10.In complex associations, the integration and coexistence of symbionts depend on supplementary partnerships and mutual contribution to the host’s metabolism41. The most intensively studied cases are highly specialized associations, where both partners can only exist in close relationship with one another. The relatively high diversity of microbes in R. alisae complicates understanding the complex pattern of molecular and cellular interactions between the host and its symbionts. More

1.Barrington, D. J. & Ghadouani, A. Application of hydrogen peroxide for the removal of toxic cyanobcteria and other phytoplankton from waste water. Environ. Sci. Technol. 4(23), 8916–8921 (2008).ADS 

Google Scholar 
2.Lurling, M., Meng, D. & Fassen, E. L. Effects of hydrogen peroxide and ultrasound on biomass reduction and toxin release in cyanobacterium, Microcytis aeruginosa. Toxins 6(12), 3260–3281 (2014).PubMed 
PubMed Central 

Google Scholar 
3.Ghime, D. & Ghosh, P. Advanced oxidation processes: A powerful treatment option for the removal of recalcitrant organic compounds. In Advanced Oxidation Processes-Applications, Trends, and Prospects (IntechOpen, 2020).4.Rahdar, S., Igwegbe, C. A., Ghasem, M. & Ahmadi, S. Degradation of aniline by the combined process of ultrasound and hydrogen peroxide (US/H2O2). MethodsX 6, 492–499 (2019).PubMed 
PubMed Central 

Google Scholar 
5.Derakhshan, Z. et al. Evaluation of kenaf fibers as moving bed biofilm carriers in algal membrane photobioreactor. Ecotoxicol. Environ. Saf. 152, 1–7 (2018).PubMed 
CAS 

Google Scholar 
6.Shekoohiyan, S. et al. Performance evaluation of cyanobacteria removal from water reservoirs by biological method. Afr. J. Microbiol. Res. 7(17), 1729–1734 (2013).CAS 

Google Scholar 
7.Cooper, W. J., Zika, R., Petasne, R. G. & Plane, J. M. Photochemical formation of hydrogen peroxide in natural waters exposed to sunlight. Environ. Sci. Technol. 22, 1156–1160. https://doi.org/10.1021/es00175a004 (1988).ADS 
Article 
PubMed 
CAS 

Google Scholar 
8.Cooper, W. J., Lean, D. R. S. & Carey, J. H. Spatial and temporal patterns of hydrogen peroxide in lake waters. Can. J. Fish. Aquat. Sci. 46, 1227–1231. https://doi.org/10.1139/f89-158 (1989).Article 
CAS 

Google Scholar 
9.Cory, R. M. et al. Seasonal dynamics in dissolved organic matter, hydrogen peroxide, and cyanobacterial blooms in Lake Erie. Front. Mar. Sci. https://doi.org/10.3389/fmars.2016.00054 (2016).Article 

Google Scholar 
10.Caverzan, A. et al. Plant responses to stresses: Role of ascorbate peroxidase in the antioxidant protection. Genet. Mol. Biol. 35(4), 1011–1019 (2012).PubMed 
PubMed Central 
CAS 

Google Scholar 
11.Sharma, P., Jha, A. B., Dubey, R. S. & Pessarakli, M. Reactive oxygen species, oxidative damage, and antioxidative defense mechanism in plants under stressful conditions. J. Bot. 2012, 1–26 (2012).
Google Scholar 
12.Ugya, A. Y., Imam, T. S., Li, A., Ma, J. & Hua, X. Antioxidant response mechanism of freshwater microalgae species to reactive oxygen species production: A mini review. J. Chem. Ecol. 36(2), 174–193 (2020).CAS 

Google Scholar 
13.Rastogi, R. P., Singh, S. P., Häder, D.-P. & Sinha, R. P. Detection of reactive oxygen species (ROS) by the oxidant-sensing probe 2′,7′-dichlorodihydrofluorescein diacetate in the cyanobacterium Anabaena variabilis PCC 7937. Biochem. Biophys. Res. Commun. 397(3), 603–607 (2010).PubMed 
CAS 

Google Scholar 
14.Foyer, C. H. Reactive oxygen species, oxidative signaling and the regulation of photosynthesis. Environ. Exp. Bot. 154, 134–142 (2018).PubMed 
PubMed Central 
CAS 

Google Scholar 
15.Gill, S. S. & Tuteja, N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 48(12), 909–930 (2010).PubMed 
CAS 

Google Scholar 
16.Ma, Z. & Gao, K. Spiral breakage and photoinhibition of Arthrospira platensis (Cyanophyta) caused by accumulation of reactive oxygen species under solar radiation. Environ. Exp. Bot. 68(2), 208–213 (2010).CAS 

Google Scholar 
17.Welkie, D. G. et al. A hard day’s night: Cyanobacteria in diel cycles. Trends Microbiol. 27(3), 231–242 (2019).PubMed 
CAS 

Google Scholar 
18.Latifi, A., Ruiz, M. & Zhang, C. C. Oxidative stress in cyanobacteria. FEMS Microbiol. Rev. 33(2), 258–278 (2009).PubMed 
CAS 

Google Scholar 
19.Lea-Smith, D. J., Bombelli, P., Vasudevan, R. & Howe, C. J. Photosynthetic, respiratory and extracellular electron transport pathways in cyanobacteria. Biochim. Biophys. Acta (BBA) Bioenerg. 1857(3), 247–255 (2016).CAS 

Google Scholar 
20.Raja, V., Majeed, U., Kang, H., Andrabi, K. I. & John, R. Abiotic stress: Interplay between ROS, hormones and MAPKs. Environ. Exp. Bot. 137, 142–157 (2017).CAS 

Google Scholar 
21.Asada, S., Fukuda, K., Oh, M., Hamanishi, C. & Tanaka, S. Effect of hydrogen peroxide on the metabolism of articular chondrocytes. Inflamm. Res. 48(7), 399–403 (1999).PubMed 
CAS 

Google Scholar 
22.Nishiyama, Y. & Murata, N. Revised scheme for the mechanisms of photoinhibition and its application to enhance the abiotic stress tolerance of the photosynthetic machinery. Appl. Microbiol. Biotechnol. 98(21), 8777–8796 (2014).PubMed 
CAS 

Google Scholar 
23.Mikula, P., Zezulka, S., Jancula, D. & Marsalek, B. Metabolic activity and membrane integrity changes in Microcystis aeruginosa—New findings on hydrogen peroxide toxicity in cyanobacteria. Eur. J. Phycol. 47(3), 195–206 (2012).CAS 

Google Scholar 
24.Huisman, J. & Hulot, F. D. Population dynamics of harmful cyanobacteria. In Harmful Cyanobacteria, 143–176 (Springer, 2005).25.Bergström, A. K. The use of TN:TP and DIN:TP ratios as indicators for phytoplankton nutrient limitation in oligotrophic lakes affected by N deposition. Aquat. Sci. 72(3), 277–281 (2010).
Google Scholar 
26.Downing, J. A. & McCauley, E. The nitrogen: Phosphorus relationship in lakes. Limnol. Oceanogr. 37(5), 936–945 (1992).ADS 
CAS 

Google Scholar 
27.Horne, A. J. & Goldman, C. R. Limnology Vol. 2 (McGraw-Hill, 1994).
Google Scholar 
28.Mittler, R. Abiotic stress, the field environment and stress combination. Trends Plant Sci. 11(1), 15–19. https://doi.org/10.1016/j.tplants.2005.11.002 (2006).Article 
PubMed 
CAS 

Google Scholar 
29.Saints, M., Diaz, P., Monza, J. & Borsani, O. Heat stress results in loss of chloroplast Cu/Zn superoxide dismutase and increased damage to Photosystem II in combined drought-heat stressed Lotus japonicus. Physiol. Plant 140(1), 46–56. https://doi.org/10.1111/j.1399-3054.2010.01383.x (2010).Article 
CAS 

Google Scholar 
30.Suzuki, N., Rivero, R. M., Shulaev, V., Blumwald, E. & Mittler, R. Abiotic and biotic stress combinations. New Phytol. 203(1), 3–43. https://doi.org/10.1111/nph.12797 (2014).Article 

Google Scholar 
31.Asaeda, T. & Barnuevo, A. Oxidative stress as an indicator of niche-width preference of mangrove Rhizophora stylosa. For. Ecol. Manag. 432, 73–82 (2019).
Google Scholar 
32.Asaeda, T., Senavirathna, M. D. H. J., Vamsi Krishna, L. & Yoshida, N. Impact of regulated water levels on willows (Salix subfragilis) at a flood-control dam, and the use of hydrogen peroxide as an indicator of environmenal stress. Ecol. Eng. 127, 96–102 (2019).
Google Scholar 
33.Asaeda, T., Senavirathna, M. D. H. J. & Vamsi Krishna, L. Evaluation of habitat preferance of invasive macrophyte Egeria densa in different channel slopes using hydrogen peroxide as an indicator. Front. Plant Sci. 11, 422. https://doi.org/10.3389/fpls.2020.00422 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
34.Diaz, J. & Plummer, S. Production of extracellular reactive oxygen species by phytoplankton: Past and future directions. J. Plankton Res. 40(6), 655–666 (2018).PubMed 
PubMed Central 
CAS 

Google Scholar 
35.Drábková, M., Admiraal, W. & Maršálek, B. Combined exposure to hydrogen peroxide and PAR selective effects on cyanobacteria, green algae, and diatoms. Environ. Sci. Technol. 41(1), 309–314 (2007).ADS 
PubMed 

Google Scholar 
36.Bouchard, J. N. & Purdie, D. A. Effect of elevated temperature, darkness and hydrogen peroxide treatment on oxidative stress and cell death in the bloom-forming toxic cyanobacterium Microcystis aeruginosa. J. Phycol. 47(6), 1316–1325 (2011).PubMed 
CAS 

Google Scholar 
37.Leunert, F., Eckert, W., Paul, A., Gerhardt, V. & Grossart, H. P. Phytoplanktonic response to UV-generated hydrogen peroxide from natural organic matter. J. Plankton Res. 36(1), 185–197. https://doi.org/10.1093/plankt/fbt096 (2014).Article 
CAS 

Google Scholar 
38.Wang, B. et al. Optimization method for Microcystis bloom mitigation by hydrogen peroxide and its stimulative effects on growth of chlorophytes. Chemosphere 228, 503–512 (2019).ADS 
PubMed 
CAS 

Google Scholar 
39.Foo, S. C., Chapman, I. J., Hartnell, D. M., Turner, A. D. & Franklin, D. J. Effects of H2O2 on growth, metabolic activity and membrane integrity in three strains of Microcystis aeruginosa. Environ. Sci. Pollut. Res. 27(31), 38916–38927 (2020).CAS 

Google Scholar 
40.Barrington, D. J., Reichwaldt, E. S. & Ghadouani, A. The use of hydrogen peroxide to remove cyanobacteria and microcystins from waste stabilization ponds and hypereutrophic systems. Ecol. Eng. 50, 86–94 (2013).
Google Scholar 
41.Drábková, M., Matthijs, H., Admiraal, W. & Maršálek, B. Selective effects of H2O2 on cyanobacterial photosynthesis. Photosynthetica 45(3), 363–369 (2007).
Google Scholar 
42.Marsac, N. T. D. Occurrence and nature of chromatic adaptation in cyanobacteria. J. Bacteriol. 130(1), 82–91 (1977).
Google Scholar 
43.Garcia, P. E., Queimalinos, C. & Dieguez, M. C. Natural levels and photo-production rates of hydrogen peroxide (H2O2) in Andean Patagonian aquatic sysyems: Influence of the dissolved organic matter pool. Chemosphere 217, 550–557 (2019).ADS 
PubMed 
CAS 

Google Scholar 
44.Herrmann, R. The daily changing pattern of hydrogen peroxide in New Zealand surface waters. Environ. Toxicol. Chem. 15(5), 652–662 (1996).CAS 

Google Scholar 
45.Spoof, L. et al. Elimination of cyanobacteria and microcystins in irrigation water—Effects of hydrogen peroxide treatment. Environ. Sci. Pollut. Res. 27(8), 8638–8652. https://doi.org/10.1007/s11356-019-07476-x (2020).Article 
CAS 

Google Scholar 
46.Lopez, C. V. G. et al. Protein measuremements of microalgae and cyanobacterial biomass. Bioresour. Technol. 101(19), 7587–7591 (2010).PubMed 

Google Scholar 
47.Vesterkvist, P. S. M., Misiorek, J. O., Spoof, L. E. M., Toivola, D. M. & Meriluoto, J. A. O. Comparative cellular toxicity of hydrophilic and hydrophobic microcystins on Caco-2 cells. Toxins 4(11), 1008–1023 (2012).PubMed 
PubMed Central 
CAS 

Google Scholar 
48.Preece, E. P., Hardy, F. J., Moore, B. C. & Bryan, M. A review of microcystin detections in estuarine and marine waters: Environmental implications and human health risk. Harmful Algae 61, 31–45 (2017).CAS 

Google Scholar 
49.Pham, T.-L. & Utsumi, M. An overview of the accumulation of microcystins in aquatic ecosystems. J. Environ. Manag. 213, 520–529 (2018).CAS 

Google Scholar 
50.Goldman, J. C., McCarthy, J. J. & Peavey, D. G. Growth rate influence on the chemical composition of phytoplankton in oceanic waters. Nature 279(5710), 210–215 (1979).ADS 
CAS 

Google Scholar 
51.Paerl, H. W., Fulton, R. S. 3rd., Moisander, P. H. & Dyble, J. Harmful freshwater algal blooms, with an emphasis on cyanobacteria. Sci. World. J. 1, 76–113 (2001).CAS 

Google Scholar 
52.Xie, L., Xie, P., Li, S., Tang, H. & Liu, H. The low TN:TP ratio, a case or result of Microcystis blooms?. Water Res. 37(9), 2073–2080 (2003).PubMed 
CAS 

Google Scholar 
53.Asaeda, T., Rashid, M. H. & Schoelynck, J. Tissue hydrogen peroxide concentration can explain the invasiveness of aquatic macrophytes: A modeling perspective. Front. Environ. Sci. 8, 292 (2021).ADS 

Google Scholar 
54.Hesse, K., Dittman, E. & Borner, T. Consequences of impaired microcystin production for light-dependent growth and pigmentation of Microcystis aeruginosa PCC 7806. FEMS Microbiol. Ecol. 37(1), 39–43 (2001).CAS 

Google Scholar 
55.Tilzer, M. M. Light‐dependence of photosynthesis and growth in cyanobacteria: Implications for their dominance in eutrophic lakes. N. Z. J. Mar. Freshwater Res. 21(3), 401-412 (1987).Article 
CAS 

Google Scholar 
56.Iwase, S. & Abe, Y. Identification and change in concentration of musty-odor compounds during growth in blue–green algae. J. Mar. Sci. Technol. 8(1), 27–33 (2010).
Google Scholar 
57.Abeynayaka, H. D. L., Asaeda, T. & Kaneko, Y. Buoyancy limitation of filamentous cyanobacteria under prolonged pressure due to the gas vesicle collapse. Environ. Manag. 60(2), 293–303 (2017).ADS 

Google Scholar 
58.Rippka, R., Deruelles, J., Waterbury, J. B., Herdman, M. & Stanier, R. Y. Generic assignments, strain histories and properties of pure cultures of cyanobacteria. Microbiology 111(1), 1–61 (1979).
Google Scholar 
59.Jana, S. & Choudhuri, M. A. Glycolate metabolism of three submersed aquatic angiosperms during ageing. Aquat. Bot. 12, 345–354 (1982).CAS 

Google Scholar 
60.Veljovic-Jovanovic, S., Noctor, G. & Foer, C. H. Are leaf hydrogen peroxide concentrations commonly overestimated? The potential influence of artefactual interference by tissue phenolics and ascorbate. Plant Physiol. Biochem. 40, 501–507 (2002).CAS 

Google Scholar 
61.Cheeseman, J. M. Hydrogen peroxide concentrations in leaves under natular conditions. J. Exp. Bot. 57(10), 2435–2444 (2006).PubMed 
CAS 

Google Scholar 
62.Queval, G., Hager, J., Gakiere, B. & Noctor, G. Why are literature data for H2O2 contents so variable? A discussion of potential difficulties in the quantitative assay of leaf extracts. J. Exp. Bot. 59(2), 135–146. https://doi.org/10.1093/jxb/erm193 (2008).Article 
PubMed 
CAS 

Google Scholar 
63.Aebi, H. Catalase in vitro. Methods Enzymol. 105, 121–126 (1984).PubMed 
CAS 

Google Scholar 
64.Nakano, Y. & Asada, K. Hydrogen peroxide is scavenged by ascorbate-specific peroxidase in spinach chloroplasts. Plant Cell Physiol. 22(5), 867–880 (1981).CAS 

Google Scholar 
65.Ahmad, P., Jaleel, C. A., Salem, M. A., Nabi, G. & Sharma, S. Roles of enzymatic and non enzymatic antioxidants in plants during abiotic stress. Crit. Rev. Biotechnol. 30(3), 161–175 (2010).PubMed 
CAS 

Google Scholar  More

Description of study areasIGR plant (33° 42/N, 46° 13/E) is located along the edge of the mountains of Zagros forests and 25 km from Ilam city. Its main activity, to supply gas to the western provinces of Iran, started in 2007. It converts sour gas to sweet gas and also produces various products such as pastil sulfur, ethane, and liquefied gas. The refinery has two chimneys, which release waste gases into the atmosphere. Oak trees are the main tree species of the Zagros forests around the refinery; these are exposed to various air pollutants and different elements from this source. Based on random analysis of exhaust emissions, sulfur dioxide and sulfide hydrogen are the major pollutants emitted from the flare gases of this refinery plant34. The sampling points have an average altitude of about 1000–1250 m and a slope of less than 20%. The climate of the region is semiarid and influenced by Mediterranean winds. The predominant wind direction was west and southwest. The highest and lowest air temperatures were 41.4 °C and − 11.3 °C, respectively. The average annual rainfall was 71.94 mm (http://www.amarilam.ir).Samples collection and analysesAll methods were carried out in accordance with the relevant institutional, national, and international guidelines and legislation. Besides they were discussed and approved by the Research Ethics Committee of Tarbiat Modares University. The formal identification of the Quercus brantii Lindl. was performed by H. Dadkhah-Aghdash based on colorful Flora of Iran35. The permissions or licenses to collect Brant oak (Quercus brantii Lindl.) trees in Zagros forests were obtained. A voucher specimen of Brant oaks were collected and deposited at the Herbarium of department of Plant Biology of Tarbiat Modares University.We studied different distances (1000, 1500, 2000, 2500, and 10,000 m [control]) in an easterly direction from the gas refinery. The map of study area was drawn by software of ArcGIS version of 10.5, https://desktop.arcgis.com (Fig. 5). At each distance, three soil samples taken from the depth of 0–20 cm with a plastic gardening shovel, 30 healthy and mature leaves were collected from a certain height (nearly the middle of the canopy) and the outer canopy of three Brant oak trees in the late spring, summer, and autumn of 2019. These trees with average height and diameter at breast height of 5.5 m and 45 cm were selected randomly. The leaf and soil samples were put into polyethylene bags and transported to the laboratory for analysis36.Figure 5Locations of collection sites of soil samples and Brant oak leaves at five different distances (1000, 1500, 2000, 2500 and 10,000 m [control]) from the gas refinery (drawn by H. Dadkhah-Aghdash using software of ArcGIS Desktop. version of 10.5. ESRI, California, US. https://desktop.arcgis.com).Full size imageIn the lab, firstly the leaves were categorized into two types: unwashed leaves and leaves washed with ethylenediaminetetraacetic acid (EDTA) solution to remove some atmospheric dusts and particles deposition. The leaf and soil samples were dried for 10 days until they reached a constant weight at lab temperature. The leaves were grinded and homogenized, soils were sieved with ASTM mesh (DAMAVAND, Iran) with a diameter of 2 mm and homogenized.To determine the pH and electrical conductivity (EC) of soils, 2 g of the soil samples were shaken in 10 ml of double-distilled water with a ratio of 1:5; after 1 h, the pH and electrical conductivity (EC) of the solution were measured by a digital pH meter (Fan Azma Gostar Company, Iran) and EC meter (Sartorius, PT-20, USA). The analysis of the particle sizes of the soil was carried out using the hydrometer method and texture class was determined with a soil texture triangle37.According to different U.S.EPA protocols that were modified by following references, the soil and leaf samples were prepared and dissolved. The digestion of soil samples was conducted with a mixture of concentrated HF–HClO4–HNO338. Approximately 0.5 g of dry soil sample was digested with 10 mL of HCl on a hot plate at ~ 180 °C until the solution was reduced to 3 mL. Approximately 5 mL of HF (40%, w/w), 5 mL of HNO3 (63%, w/w), and 3 mL of HClO4 (70%, w/w) were then added and the solution was digested. This process was continued with adding 3 mL of HNO3, 3 mL of HF, and 1 mL of HClO4 until the silicate minerals had fully disappeared. This solution was transferred to a 25 mL volumetric tube, and 1% HNO3 was added to bring the sample up to a constant volume for the element’s determinations. After filtering the digested samples, the concentrations of sulfur (S), arsenic (As), chromium (Cr), copper (Cu), lead (Pb), zinc (Zn), manganese (Mn), and nickel (Ni) were measured via inductively coupled plasma mass spectrometry (ICP-MS,7500 CS, Agilent, US). The procedures of quality assurance and quality control (QA/QC) were performed.To quantify element contents from soil samples, external standards with calibration levels were used. The precision and the repeatability of the analysis were tested on the instrument by analyzing three replicate samples.According to Liang et al.39 leaf samples were acid digested and sieved powder samples were placed in the acid-washed tubes and 10 mL of 65% nitric acid was added to it. The solution was placed at room temperature overnight (12 h) after that, it was placed for 4 h at 100 °C and then 4 h at 140 °C until the solution color was clear. After cooling, the solution was diluted by deionized water to 50 mL and then passed through Whatman filter paper until 25 mL of the filtrate volume was provided. Each sample was digested three times and the average of measurements is reported. Total plant elements were measured by using the ICP-MS (7500 CS, Agilent, US). A control sample was also used beside each sample to determine the background pollution during digestion. To confirm the accuracy of the methodology and to ensure the extraction of trace elements from the leaf samples, the standard solution of each studied elements was used.Measuring of pollution levels of different elements in soils and leavesFor assessment of contamination levels (concentration) of different elements in soils and trees, common indices of pollution including geoaccumulation index (Igeo), pollution index (PI), pollution load index (PLI), enrichment factor of plants (EFplant), bioconcentration factor (BCF), air originated metals (AOM ), metal accumulation index (MAI) were used.Igeo was calculated using the following (Eq. 1):$${text{I}}_{{{text{geo}}}} = log_{2} left[ {{text{C}}_{{text{n}}} / 1.5{text{ B}}_{{text{n}}} } right]$$
(1)
where Cn is the measured concentration of the element n, Bn is the geoaccumulation background for this element and 1.5 is a constant coefficient used to eliminate potential variations in the baseline data40. The Igeo classifies samples into seven grades:  5 for extremely polluted30.The first PI is expressed as (Eq. 2):$${text{PI }} = {text{ C}}_{{text{i}}} /{text{S}}_{{text{i}}}$$
(2)
where Ci is the concentration of element i in the soil (mg kg−1) and Si is the soil quality standard or reference value for element i (mg kg−1). The PLI for different elements is calculated via the (Eq. 3):$${text{PLI}} = left( {{text{PI}}_{{1}} times {text{ PI}}_{{2}} times {text{ PI}}_{{3}} times cdots times {text{PI}}_{{text{n}}} } right)^{{{1}/{text{n}}}}$$
(3)
The PLI of soils is classified as follows: PLI  More

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