Abstract
In grazing conditions, unpalatable species may induce either associational defense or neighbor contrast susceptibility in neighboring communities. Using surveys from eight grasslands, we tested whether various unpalatable species have the same impacts on neighboring communities in response to grazing. The studied unpalatable species were: Phlomis cancellata (an unpalatable nonpoisonous plant), Euphorbia boissieriana, E. microsciadia (poisonous plants), and Seseli transcaucasicum (a highly poisonous plant). Our results showed that, in the ungrazed grasslands, communities containing P. cancellata had lower biodiversity than communities without it. In the moderately- and heavily grazed grasslands, P. cancellata induced associational defense in the neighboring communities. In heavily grazed grasslands, both Euphorbia species promoted neighbor contrast susceptibility in the neighboring communities. Similarly, S. transcaucasicum in a heavily grazed grassland, induced neighbor contrast susceptibility. Different responses of plant community vulnerability among the studied unpalatable plants might be due to herbivore different foraging decisions. Accordingly, grazers selectively choose from other patches when facing P. cancellata and other plant individuals when there is a poisonous plant in a patch. Our results suggested that grazing intensity may not substantially affect the foraging decisions of sheep and goats in response to unpalatable species. We recommend monitoring the abundance of poisonous species to maintain the sustainable use of grasslands.
Similar content being viewed by others
Introduction
Foraging decisions of large herbivores in order to select species with high nutrients and energy affect plant communities1. Plant species adapt to grazing by using morphological and chemical defensive traits2. Defensive traits of neighboring species modify the degree of protection of a focal species against large herbivores3,4,5. Thus, whether species gain defense or become more susceptible to grazing also depends on its neighbors5,6.
In grazing conditions, plant species in a community can promote four types of associational effects: (1) Associational defense, plants can reduce herbivore damage by growing closely to unpalatable neighbors. (2) Associational susceptibility, plants experience increased damage by growing together with more palatable species. (3) Neighbor contrast defense, plants will be less chosen by neighboring with more palatable species. (4) Neighbour contrast susceptibility, plants will be more susceptible to herbivory by neighboring with a less palatable species1. Therefore, unpalatable species can induce either associational defense or neighbor contrast susceptibility. For livestock grazing, associational defense is more likely to happen than neighbor contrast susceptibility5.
Previous studies on evaluating the effects of unpalatable species on neighboring communities have tended to focus on plants that use structural defensive traits (e.g., Refs.7,8,9). Consequently, species using chemical defensive traits have rarely been studied. Moreover, there is inconsistency among the reported findings on the impacts of poisonous species on neighboring communities. For example, both negative (e.g., Refs.10,11) and positive impacts—e.g., Refs.4,10,12—have been reported.
Most of the research have been limited to study a single focal species, and the effects of unpalatable species at the community level received lesser attention13,14. Although phylogenetic diversity can be altered by different plant–plant interactions15, no studies have compared this measure between communities containing or without a chemically unpalatable or poisonous species to the best knowledge of the authors. Generally, in the literature on associational effects of unpalatable species, effects of grazing intensity have been neglected. Nevertheless, it is well-established that interaction between plants is likely to be changed with different grazing intensities6,16.
Here, we use eight datasets that were collected from ungrazed, moderately-, and heavily grazed sites in northeastern Iran. We evaluated the species and phylogenetic diversities as well as the species composition of plant communities containing and without four unpalatable plants: a nonpoisonous chemically unpalatable, two poisonous, and a highly poisonous species. We hypothesized that chemically unpalatable plants—with different toxicity levels—do not show the same associational effect. The aim of this study was to quantify the impacts of these species on biodiversity and species composition of neighboring communities. We compare these communities with communities without unpalatable species in different levels of large herbivore grazing.
Materials and methods
Study areas and datasets
We used data from eight surveys (Table 1) that were sampled from various landscapes in the mountainous ranges of northeastern Iran. Herbaceous plants were the dominant species in the selected sites. Land use change and grazing by sheep and goats were the primary disturbances in the study sites17,18,19,20,21,22. Two sites were partially-managed: the Kelilagh no-hunting zone (hereafter is called the Kelilagh) and the Heydari Wildlife Refuge (hereafter is called the Heydari). The other sites were neither protected nor restored. Information about all sites is presented in Table 1. The selected sites represent moderately- (i.e., the Fereizi and Heydari sites), heavily grazed (i.e., the Darrud, Arabchah, Zharf, and Boghmech sites), and ungrazed areas (i.e., the Najafi and Kelilagh sites).
Selected species
We considered the literature data as well as field observations to select unpalatable plants. Four species, ranging from nonpoisonous unpalatable to a highly poisonous plant, were selected:
-
(A)
Phlomis cancellata (Lamiaceae) is a fragrant, chemically unpalatable species. This plant is used in the traditional medicine of Iran19,23. Phlomis cancellata is a dominant species in the grazed areas of northeastern Iran19. This species is not poisonous and has a hemicryptophyte life form24.
-
(B)
Euphorbia boissieriana and E. microsciadia (Euphorbiaceae) are two poisonous species with a similar physiognomy (i.e., hemicryptophyte)24. These plants remain intact after grazing. If they are accidentally ingested, these species will cause severe harm to livestock25.
-
(C)
Seseli transcaucasicum (Apiaceae) remains intact after grazing. This plant is considered a highly poisonous species that is lethal if ingested by livestock. S. transcaucasicum has a chamaephyte life form26.
Data preparation
In each dataset, random plots containing > 10% canopy cover of unpalatable species were selected. These plots hereafter are called Contain Focal Species (CFSs). Given the size of the studied species, we estimated that there was at least one full-grown individual in plots that the focal species has a canopy cover of > 10%. Additionally, the nearest plots with no unpalatable species were selected. These plots hereafter are called Without Focal Species (WFSs). Environmental factors (e.g., elevation, aspect, slope degree) were evaluated to be the same to eliminate undesired sources of variability between CFSs and WFSs. If environmental conditions were not the same for a plot pair, this pair was excluded from further analyses. Finally, data from 180 sampling units of 1 m × 1 m and 51 sampling units of 3 m × 3 m area were analyzed. The plots of the Fereizi had a 9 m2 area, and those of the other sites had a 1 m2 area. We only compared CFSs and WFSs of each site, and no comparisons were made among the sites. We considered the vascular plant species and their canopy cover in WFSs and CFSs.
For P. cancellata, data from six sites were analyzable. For Euphorbia species, two sites and for S. transcaucasicum, one site had analyzable data. A data was considered analyzable if it had at least five CFSs-WFSs pairs. The data limitation was due to a limited distribution range of unpalatable species, e.g., S. transcaucasicum has only been recorded in the Arabchah. Further, Euphorbia species were rarely recorded from the ungrazed and moderately grazed sites. The data is presented in Supplementary file.
Data analysis
Species diversity
We compared the species diversity of the CFSs and WFSs. We used Hill numbers27. We calculated the species richness (q = 0 in the Hill numbers) and the reciprocal of the Gini-Simpson index (q = 2 in the Hill numbers). These indices evaluate the species diversity of a community at the level of rare and dominant species27, respectively. A coverage-based rarefaction and extrapolation method was employed for biodiversity calculations. This method was used to eliminate the effects of unequal sampling effort (e.g., an unequal number of samples) on our inferences28. We calculated the base coverage of each dataset. We used the method of Chao and Jost28 to calculate this value. Then, Species diversity and phylogenetic diversity indices were estimated at the base coverage. A 95% confidence interval (CI) for each index was estimated by using a bootstrap method. All species diversity calculations were conducted using the iNEXT package29. For each site, the results were reported as the percentage change in species diversity in CFSs relative to the species diversity of the WFSs.
Phylogenetic diversity
We used an angiosperm phylogeny of all vascular species that were collected in the plots. The non-angiosperms were excluded to prevent adding long branches that lead to the outliers generating in the phylogenetic diversity calculations30. The phylogenetic trees were created using the V.Phylomaker package31. We used the Phylogenetic Hill diversity indices to compare CFSs and WFSs in each site. The phylogenetic diversity at the level of rare (q = 0) and dominant (q = 2) species was calculated. We used a coverage-based rarefaction and extrapolation approach. A 95% CI was also estimated for each index. The iNEXT‐PD package was used for phylogenetic diversity calculations29,32,33. The results were reported as the percentage change in phylogenetic diversity of CFSs relative to the phylogenetic diversity of the WFSs.
Species composition
We used a transformation-based Principal Component Analysis (tb-PCA) to visualize the variation in species composition of CFSs and WFSs in each site. First, we applied the Hellinger transformation on the species data to eliminate the effects of double zeros on the next analyses34. This transformation was performed by using the vegan package35. Then, a PCA was conducted on the transformed data by using the rda function in the vegan package. We conducted an analysis of similarity (ANOSIM), number of permutations = 999, to test whether the species composition of CFSs and WFSs was significantly different. We used the anosim function in the vegan package for this analysis. All analyses were performed in R version 3.636.
Results
Phlomis cancellata
In the Najafi—an ungrazed site, CFSs showed a 75% lower species richness (q = 0) and 63% lower the reciprocal of the Gini-Simpson index (q = 2) compared to WFSs. The decrease in the species richness was significant in this site (Fig. 1). In the other ungrazed site—the Kelilagh, CFSs had a 60% lower species richness and 13% lower the reciprocal of the Gini-Simpson index than WFSs. In the Kelilagh, the species richness decrease was significant (Fig. 1). On average in the ungrazed areas, communities with P. cancellata (CFSs) had a significantly 67% lower species richness than communities without this plant (WFSs). The decreases in the reciprocal of the Gini-Simpson index in CFSs was not significant and had a value of 38%.
The phylogenetic diversity results showed that, in the ungrazed sites, CFSs had a 56% (the Najafi) and 38% (the Kelilagh) lower phylogenetic richness (q = 0) significantly compared to WFSs (Fig. 2). Considering the phylogenetic diversity of the dominant species (i.e., q = 2), CFSs had lower phylogenetic diversity than WFSs, with values of 26% (the Najafi) and 14% (the Kelilagh). The decrease in the phylogenetic diversity of dominant plants was significant in the Najafi and nonsignificant in the Kelilagh. Thus, on average, CFSs had a significantly 47% lower phylogenetic richness than WFSs in ungrazed sites. The decrease in the phylogenetic diversity of dominant species in CFSs had an average value of 20%.
In the moderately grazed sites, CFSs had an 11% higher species richness in the Fereizi and 10% higher in the Heydari than WFSs. This increase in the species richness was not significant in these sites (Fig. 1). On average, the species richness of CFSs in moderately grazed sites had a non-significantly 11% higher species richness than WFSs. Considering the reciprocal of the Gini-Simpson index, CFSs of the Fereizi and Heydari had higher diversity than WFSs with a value of 23% and 45%, respectively. This finding was nonsignificant in the Fereizi but significant in the Heydari (Fig. 1). The average increase in the reciprocal of the Gini-Simpson index in CFSs was 34% compared to WFSs in the moderately grazed sites.
In the moderately grazed areas, CFSs had higher phylogenetic richness than WFSs, with a value of 14% in the Fereizi and 21% in the Heydari. This phylogenetic richness increase was significant in neither of the sites (Fig. 2). Thus, on average CFSs had 18% higher phylogenetic richness than WFSs in the moderately grazed sites. Considering the phylogenetic diversity at the level of dominant species, the increase in diversity of CFSs for the Fereizi and Heydari was 16% and 9%, respectively. This finding was significant in the Fereizi and not significant in the Heydari (Fig. 2). Therefore, the average phylogenetic diversity of dominant species in CFSs was 12% higher compared to WFSs in moderately grazed areas.
In the heavily grazed sites, in the Darrud, CFSs had 37% higher species richness and 8% higher the reciprocal of the Gini-Simpson index than WFSs. These increases were not significant (Fig. 1). In the Zharf, species richness (8%) and the reciprocal of the Gini-Simpson index (17%) were higher in CFSs than WFSs. In this site, the difference between CFSs and WFSs was not significant (Fig. 1). On average, CFSs had non-significantly 23% higher species richness than WFSs in the heavily grazed sites. For the reciprocal of the Gini-Simpson index, the difference between CFSs and WFSs was not significant and CFSs on average had 11% higher diversity.
In the heavily grazed sites, CFSs had 44% (the Darrud) and 10% (the Zharf) higher phylogenetic richness (q = 0) than WFSs. The increase in phylogenetic richness was significant in the Darrud but not significant in the Zharf (Fig. 2). Considering the phylogenetic diversity at the level of dominant species (q = 2), CFSs had higher diversity than WFSs. This difference was significant in the Darrud and not significant in the Zharf with 22% and 7% values, respectively. On average, CFSs had 27% higher phylogenetic richness than WFSs in the heavily grazed sites. The increase in the phylogenetic diversity of dominant species in CFSs had an average value of 15%.
No significant difference was detected among the species compositions of CFSs and WFSs in the six sites. Figure 3 presents the tb-PCA and ANOSIM results for each of the study areas.
Euphorbia boissieriana and E. microsciadia
For E. boissieriana, in the Darrud—a heavily grazed site, CFSs had 19% lower species richness and 11% lower the reciprocal of the Gini-Simpson index than WFSs. The difference in species diversity was significant at neither of the levels (Fig. 3). Considering the phylogenetic diversity indices, the phylogenetic diversity of CFSs was 20% and 18% lower than WFSs for q = 0 and 2, respectively. This result was not significant at any of the q levels (Fig. 4).
For E. microsciadia, in the Boghmech (a heavily grazed site), CFSs had 0.1% (species richness) and 13% (the reciprocal of the Gini-Simpson index) lower diversity compared to WFSs. The difference in species diversity was not significant at the both levels (Fig. 4). Considering the phylogenetic diversity, CFSs had lower phylogenetic richness than WFSs. This difference was not significant and had a value of 11%. On the contrary, a 19% non-significantly increase in phylogenetic diversity of dominant species was observed in CFSs compared to WFSs.
No significant difference was detected between the species compositions of CFSs and WFSs for E. boissieriana. However, there was a significant difference between the species composition of the two communities for E. microsciadia. Figure 4 shows the species and phylogenetic diversity results along with tb-PCA results for CFSs and WFSs.
Seseli transcaucasicum
In the Arabchah—a heavily grazed site, CFSs had 9% lower species richness and 30% lower the reciprocal of the Gini-Simpson index than WFSs. This difference in species diversity was only significant at the level of dominant species (q = 2) (Fig. 5a). The phylogenetic richness (q = 0) of CFSs was non-significantly 3% lower than WFSs. Phylogenetic diversity of dominant species in CFSs was significantly 27% lower than WFSs (Fig. 5b). There was a significant difference between the species compositions of CFSs and WFSs (Fig. 5c). Table 2 summarizes our results.
Discussion
It has previously been reported that unpalatable species protect their neighbors from grazing4, our results indicated the type of association effect depends on the unpalatable species.
Phlomis cancellata: associational defense
Our results revealed that P. cancellata had the same effect on neighboring plants in areas with the same grazing condition. Our results indicated P. cancellata helped neighbors to compensate for grazing. As a result, CFSs had higher biodiversity than WFSs. However, neighboring with this plant had negative impacts in the ungrazed areas. The canopy of P. cancellata creates a shady condition that might be unfavorable in ungrazed conditions. Similar results were reported for Urtica thunbergiana. This species was a facilitator in the grazing conditions but a competitor in ungrazed areas2. It was also reported that Juncus effusus increased the diversity of neighboring communities in a grazed area and competed with neighbors in an ungrazed site. However, J. effusus affected the species composition of neighboring communities in the grazed site37. The relation of P. cancellata with neighbors in ungrazed conditions might be a competitive interaction. In the presence of large herbivores, species benefit from associating with this plant. This finding suggested that biotic stresses affect the plant–plant interactions. A similar condition was also reported for Filipendula ulmaria38.
Euphorbia boissieriana and E. microsciadia: neighbor contrast susceptibility
Our results revealed that for both studied Euphorbia species, WFSs showed non-significantly higher diversity than CFSs. These results suggested a case of neighbor contrast susceptibility. The phylogenetic diversity results implied that not every species were able to benefit from neighboring E. boissieriana or E. microsciadia, except for those species from relative taxa showing a similar evolutionary history. These two species, through biotic interaction, do not allow every plant to grow nearby. Allelopathy, which is common among Euphorbia species39,40,41,42, might be that biotic interaction. The allelopathic effect of E. microsciadia was possibly stronger than that of E. boissieriana. The species composition of communities with E. microsciadia was significantly different from communities without this plant. A study on E. schickendatzii reported that species diversity was higher in the plant communities containing E. schickendatzii when compared with communities without this plant. Also, the species composition of the two communities was significantly different12.
Seseli transcaucasicum: neighbor contrast susceptibility
The species composition of CFSs was different from that of WFSs. Furthermore, communities with S. transcaucasicum showed a lower species diversity and phylogenetic diversity than communities without this plant. Similar to the studied Euphorbia species, S. transcaucasicum induced neighbor contrast susceptibility in neighboring communities. Similar to E. boissieriana and E. microsciadia but with a higher magnitude, few species from relative lineages could grow near S. transcaucasicum. Gao et al.11 suggested that range expansion of poisonous species can result in an altered community structure.
Comparing the four species: herbivore grazing hierarchy
Herbivores make foraging decisions at different spatial scales simultaneously43,44. The variation in associational effects among the studied species was due to herbivore different foraging decisions at different spatial scales (within- or between-patches). Our finding suggested that grazers selectively choose another patch when facing P. cancellata. On the other hand, they choose from plant individuals when there was a poisonous plant (i.e., E. boissieriana, E. microsciadia, and S. transcaucasicum) in a selected patch. Furthermore, comparing different grazing levels, it can be declared that grazing intensity may not have a strong effect on the foraging decision of sheep and goats in response to unpalatable species.
Limitations
This study is an observational and not experimental research. Considering this limitation is important when extending our results to the other areas. Regarding our data, the Fereizi plots were 3 m × 3 m quadrats. We did not compare the results among the sites to avoid any biases that come from this limitation. Also, we used another data from a moderately grazed site to strengthen our inferences. For S. transcaucasicum, data from one site was available. As a result, we could not account for the variation among the sites. All study areas have similar climatic conditions (i.e., semi-arid climate). Therefore, we could not test whether different climatic conditions affect the impacts of unpalatable species on neighboring communities.
Conclusions
We have observed that the unfavorable microhabitat of P. cancellata in ungrazed areas could become a possible shelter in grazing conditions. However, neighboring with a poisonous species could increase the herbivore damage on a plant community. We encourage to perform an experimental study on the effects of these studied species on neighboring communities. We have used a coverage-based approach for comparing biodiversity between CFSs and WFSs. We suggest using this method for future studies dealing with plant community comparisons.
Data availability
The data regarding this study is presented in Supplementary Information.
References
Huang, Y. et al. Foraging responses of sheep to plant spatial micro-patterns can cause diverse associational effects of focal plant at individual and population levels. J. Anim. Ecol. 87, 863–873 (2018).
Suzuki, R. O. & Suzuki, S. N. Facilitative and competitive effects of a large species with defensive traits on a grazing-adapted, small species in a long-term deer grazing habitat. Plant Ecol. 212, 343–351 (2011).
Courant, S. & Fortin, D. Foraging decisions of bison for rapid energy gains can explain the relative risk to neighboring plants in complex swards. Ecology 91, 1841–1849 (2010).
Callaway, R. M., Kikodze, D., Chiboshvili, M. & Khetsuriani, L. Unpalatable plants protect neighbors from grazing and increase plant community diversity. Ecology 86, 1856–1862 (2005).
Barbosa, P. et al. Associational resistance and associational susceptibility: Having right or wrong neighbors. Annu. Rev. Ecol. Evol. Syst. 40, 1–20 (2009).
Graff, P., Aguiar, M. R. & Chaneton, E. J. Shifts in positive and negative plant interactions along a grazing intensity gradient. Ecology 88, 188–199 (2007).
Fidelis, A., Overbeck, G. E., Pillar, V. D. & Pfadenhauer, J. The ecological value of Eryngium horridum in maintaining biodiversity in subtropical grasslands. Austral Ecol. 34, 558–566 (2009).
Uytvanck, J. V., Maes, D., Vandenhaute, D. & Hoffmann, M. Restoration of woodpasture on former agricultural land: The importance of safe sites and time gaps before grazing for tree seedlings. Biol. Conserv. 141, 78–88 (2008).
Yu, F., Krüsi, B., Schütz, M., Schneller, J. & Wildi, O. Is vegetation inside Carex sempervirens tussocks highly specific or an image of the surrounding vegetation?. J. Veg. Sci. 17, 567–576 (2006).
Cheng, W. et al. Unpalatable weed Stellera chamaejasme L. provides biotic refuge for neighboring species and conserves plant diversity in overgrazing alpine meadows on the Tibetan Plateau in China. J. Mt. Sci. 11, 746–754 (2014).
Gao, F. et al. The expansion process of a Stellera chamaejasme population in a degraded alpine meadow of Northwest China. Environ. Sci. Pollut. Res. 26, 20469–20474 (2019).
Hierro, J. L. & Cock, M. C. Herbivore-mediated facilitation alters composition and increases richness and diversity in ruderal communities. Plant Ecol. 214, 1287–1297 (2013).
Wang, L. et al. Spatially complex neighboring relationships among grassland plant species as an effective mechanism of defense against herbivory. Oecologia 164, 193–200 (2010).
Sotomayor, D. A. & Lortie, C. J. Indirect interactions in terrestrial plant communities: Emerging patterns and research gaps. Ecosphere 6, art103 (2015).
Castillo, J. P., Verdú, M. & Valiente-Banuet, A. Neighborhood phylodiversity affects plant performance. Ecology 91, 3656–3663 (2010).
Verwijmeren, M., Smit, C., Bautista, S., Wassen, M. J. & Rietkerk, M. Combined grazing and drought stress alter the outcome of nurse: Beneficiary interactions in a semi-arid ecosystem. Ecosystems 22, 1295–1307 (2019).
Memariani, F., Joharchi, M. R., Ejtehadi, H. & Emadzade, K. A contribution to the flora and vegetation of Binalood mountain range, NE Iran: Floristic and chorological studies in Fereizi region. Ferdowsi Univ. Int. J. Biol. Sci. J. Cell Mol. Res. 1, 1–17 (2009).
Kartoolinezhad, D. & Moshki, A. Changes in Juniperus polycarpos community in response to physiographical factors (Hezarmasjed Mountain, Iran). Austrian J. For. Sci. 4, 215–232 (2014).
Erfanian, M. B., Ejtehadi, H., Vaezi, J. & Moazzeni, H. Plant community responses to multiple disturbances in an arid region of northeast Iran. Land Degrad. Dev. 30, 1554–1563 (2019).
Memariani, F. & Joharchi, M. R. Iris ferdowsii (Iridaceae), a new species of section Regelia from northeast of Iran. Phytotaxa 291, 192 (2017).
Erfanian, M.B., Sagharyan, M., Memariani, F. & Ejtehadi, H. Predicting range shifts of three endangered endemic plants of the Khorassan-Kopet Dagh floristic province under global change. Sci Rep 11, 9159 (2021).
Erfanian, M. B. et al. Plant community responses to environmentally friendly piste management in northeast Iran. Ecol. Evol. 9, 8193–8200 (2019).
Morteza-Semnani, K., Moshiri, K. & Akbarzadeh, M. The essential oil composition of Phlomis cancellata Bunge. J. Essent. Oil Res. 18, 672–673 (2006).
Atashgahi, Z., Ejtehadi, H., Mesdaghi, M. & Ghasemzadeh, F. Plant diversity of the Heydari Wildlife Refuge in northeastern Iran, with a checklist of vascular plants. Phytotaxa 340, 101–127 (2018).
Arzani, H., Motamedi, J., Aghajanlu, F., Rashtvand, S. & Zareii, A. Forage quality of important rangeland species in mountainous rangelands of Qazvin and Badman Zanjan. J. Range Watershed Manag. 69, 805–818 (2017).
Hosseini, S. S. Study of Plant Biodiversity in Relation to Physiographic Factors in Hezar Masjed Summit, Khorasan Razavi Province, NE Iran. MSc Thesis (Ferdowsi University of Mashhad, 2016) (in Persian).
Chao, A. et al. Rarefaction and extrapolation with Hill numbers: A framework for sampling and estimation in species diversity studies. Ecol. Monogr. 84, 45–67 (2014).
Chao, A. & Jost, L. Coverage-based rarefaction and extrapolation: Standardising samples by completeness rather than size. Ecology 93, 2533–2547 (2012).
Hsieh, T. C., Ma, K. H. & Chao, A. iNEXT: An R package for rarefaction and extrapolation of species diversity (Hill numbers). Methods Ecol. Evol. 7, 1451–1456 (2016).
Barber, N. A. et al. Grassland restoration characteristics influence phylogenetic and taxonomic structure of plant communities and suggest assembly mechanisms. J. Ecol. 107, 2105–2120 (2019).
Jin, Y. & Qian, H. V. PhyloMaker: An R package that can generate very large phylogenies for vascular plants. Ecography https://doi.org/10.1111/ecog.04434 (2019).
Chao, A., Chiu, C.-H. & Jost, L. Phylogenetic diversity measures based on Hill numbers. Philos. Trans. R. Soc. B Biol. Sci. 365, 3599–3609 (2010).
Chao, A. et al. Rarefaction and extrapolation of phylogenetic diversity. Methods Ecol. Evol. 6, 380–388 (2015).
Legendre, P. & Legendre, L. F. J. Numerical Ecology (Elsevier, 2012).
Oksanen, J. et al. vegan: Community Ecology Package. R package. https://CRAN.R-project.org/package=vegan (2017).
R Core Team. R: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, 2019).
Boughton, E. H., Quintana-Ascencio, P. F. & Bohlen, P. J. Refuge effects of Juncus effusus in grazed, subtropical wetland plant communities. Plant Ecol. 212, 451–460 (2011).
Axelsson, E. P. & Stenberg, J. A. Associational resistance mediates interacting effects of herbivores and competitors on fireweed performance. Basic Appl. Ecol. 15, 10–17 (2014).
Liu, Q. et al. Allelochemicals in the rhizosphere soil of Euphorbia himalayensis. J. Agric. Food Chem. 62, 8555–8561 (2014).
Gholamalipour Alamdari, E., Seifolahi, B., Avarseji, Z. & Biabavi, A. Evaluation of allelopathic effect of Euphorbia maculata weed on traits of germination, chlorophyll and carotenoids pigments of wheat cultivars. Iran. J. Seed Res. 5, 71–85 (2018).
Steenhagen, D. A. & Zimdahl, R. L. Allelopathy of leafy spurge (Euphorbia esula). Weed Sci. 27, 1–3 (1979).
da Silva, U. P., Furlani, G. M., Demuner, A. J., da Silva, O. L. M. & Varejão, E. V. V. Allelopathic activity and chemical constituents of extracts from roots of Euphorbia heterophylla L. Nat. Prod. Res. 33, 2681–2684 (2019).
Miller, A. M., McArthur, C. & Smethurst, P. J. Spatial scale and opportunities for choice influence browsing and associational refuges of focal plants. J. Anim. Ecol. 78, 1134–1142 (2009).
Huang, Y., Wang, L., Wang, D., Zeng, D.-H. & Liu, C. How does the foraging behavior of large herbivores cause different associational plant defenses?. Sci. Rep. 6, 20561 (2016).
Acknowledgements
This study was financially supported by the Ferdowsi University of Mashhad. The authors wish to thank A. Basiri for his assistance in the field sampling. Special thanks to M.R. Joharchi and Y. Nasseh for the identification of certain plant species. Sincere thanks to F. Sharifi for her comments during the revision of this manuscript. We are also most grateful to four anonymous reviewers whose constructive comments helped improve the manuscript.
Author information
Authors and Affiliations
Contributions
M.B.E. and H.E. designed the study. H.E., F.M., Z.A., M.M., M.S., M.D., M.H., and S.H. contributed data. M.B.E. analysed the data, prepared the figures, and wrote the draft. The final manuscript was revised and approved by all authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Erfanian, M.B., Memariani, F., Atashgahi, Z. et al. Unpalatable plants induce a species-specific associational effect on neighboring communities. Sci Rep 11, 14380 (2021). https://doi.org/10.1038/s41598-021-93698-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-93698-4
This article is cited by
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.