Skip to main content

Species-specific interference exerted by the shrub Cistus clusii Dunal in a semi-arid Mediterranean gypsum plant community



The gypsovag shrub Cistus clusii is locally dominant in semi-arid gypsum plant communities of North-Eastern Spain. This species commonly grows in species-poor patches even though it has nurse potential, suggesting interference on neighbouring species. Other Cistus species exert a chemically mediated interference on plant communities, suggesting that it might be a common phenomenon in this genus. This study aimed investigating whether C. clusii exerts chemically mediated interference on neighbouring species in gypsum plant communities. We tested in a greenhouse whether aqueous extracts from C. clusii leaves (L), roots (R) and a mixture of both (RL) affected germination, seedling survival, and growth of nine native species of gypsum communities, including C. clusii itself. We further assessed in the field richness and abundance of plants under the canopy of C. clusii compared to Gypsophila struthium (shrub with a similar architecture having a nurse role) and in open patches. Finally, we specifically assessed in the field the influence of C. clusii on the presence of the species tested in the greenhouse experiment.


Aqueous extracts from C. clusii (R and RL) negatively affected either germination or survival in four of nine species. In the field, richness and abundance of plants were lower under the canopy of C. clusii than under G. struthium, but higher than in open patches. Specifically, five of nine species were less frequent than expected under the canopy of C. clusii.


Cistus clusii shows species-specific interference with neighbouring species in the community, which may be at least partially attributable to its phytotoxic activity. To our knowledge, this is the first report of species-specific interference by C. clusii.


Both interference and facilitation influence the composition and structure of plant communities [1, 2]. The interplay between interference and facilitation may be particularly important in arid and semi-arid environments, where abiotic conditions make facilitation an important process that affects local species composition [3]. Shrubs can act as nurse plants by providing under their canopy favourable microhabitats for plants and thus can harbour highly diverse microcosms [4, 5]. In early developmental stages, plants are highly vulnerable to abiotic stress and often need facilitation by nurse plants to establish [6, 7]. Once the facilitated seedlings become adults, they may exert an adverse effect on nurse plants by competitive exclusion for water, nutrients, light or space [8,9,10,11].

Some plant species produce chemical compounds that are released to the local environment through volatilisation, leaf leachates, root exudates or leaf litter decomposition [12,13,14]. Although these compounds can have positive effects by promoting plant growth [15] or increasing species richness [16], they are usually phytotoxic and act as selective agents that affect the performance of other species negatively [17]. The most common phytotoxic compounds found in plants are terpenes and phenolic compounds [18, 19] that may inhibit or reduce germination capacity, may cause a delay in germination time and may hamper root elongation or nutrient absorption thereby reducing plant survival and growth [20,21,22,23]. The release of phytotoxic compounds may have ecological implications because non-competitive species can take advantage of the affected plants since the latter become less competitive [24]. Therefore, chemically-mediated interference might be a way for nurse plants to gain a competitive advantage over neighbouring plants that may be potentially competitive [17].

In arid and semi-arid areas of the Mediterranean region, chemical interference is likely to be a common phenomenon, given the abundance of aromatic plants that produce potential phytotoxic compounds [25,26,27]. This production may be promoted by stressful conditions as lack of water or nutrients, salinization and high solar radiation [28]. Furthermore, in arid and semi-arid environments, the effect of those compounds may increase because of their relative accumulation in the soil [29] and an intensification of plant sensitivity [28]. The accumulation of phytotoxic compounds in the soil may have incremental effects along plant life-span [29], being the longer-lived plants more affected by phytotoxic effects due to more prolonged exposure.

Mediterranean gypsum plant communities are mainly composed of well-adapted gypsophytes (i.e., gypsum soil specialists). However, in those communities the gypsovag (i.e., non-specialist) rosemary-leaved rockrose Cistus clusii Dunal often forms locally dominant populations that are associated with species-poor plant communities. Cistus clusii, a multi-branched perennial shrub (0.5–1.0 m tall), is distributed throughout the western Mediterranean region on alkaline soils, including gypsum, marls and limestones [30], and is very tolerant to dry environments [31]. This shrub may have nurse potential because it provides shade under its canopy due to its multi-branched architecture, creating favourable microenvironments in which other species can establish [32]. Nevertheless, C. clusii commonly grows in patches isolated from other species (personal observation) which suggests that it may exert interference on other plant species in the community.

Several studies have shown that other Cistus species have phytotoxic effects, either inhibiting the germination and growth of hetero-specific seedlings via foliar exudates (e.g., Cistus ladanifer; [33,34,35]) or immobilising nutrients in the soil via root exudates (e.g., Cistus albidus; [36]). Based on the species-poor patches of C. clusii observed in gypsum plant communities and the phytotoxicity of other Cistus species, we postulated that the gypsovag C. clusii could exert phytotoxic effects on other plants beneath its canopy. To date, no study has demonstrated phytotoxicity in C. clusii, even though it is known to produce phenolic compounds and terpenes [37,38,39].

The aim of this study was to test whether or not C. clusii interferes with neighbouring species in gypsum plant communities through chemical mechanisms of interference. A controlled seeding experiment was performed in a greenhouse to identify potential phytotoxic effects of aqueous extracts from C. clusii leaves, roots and a mixture of both plant tissues on the development of neighbouring species in the community. We hypothesized that (1) aqueous extracts from C. clusii would affect those species negatively and more specifically, we expected that (2) the effect of leaf extracts would be more intense than the effect of root extracts, as previously shown for other species [40]. Moreover, we expected the effect of the mixture of leaf and root extracts to be the most significant through a combined activity of both extracts. We predicted that (3) affections by aqueous extracts from C. clusii would be exhibited on seed germination inhibition or delay, early survival decline and plant growth reduction [20,21,22,23]. Complementary to the greenhouse experiment, plant spatial associations were evaluated in the field in the local vicinity of C. clusii, compared to a shrub species of similar architecture in the community. Plant spatial associations were assessed as an indicator for biotic interactions to disentangle the potential interference exerted by this shrub on neighbouring plants in gypsum plant communities. We hypothesized that (4) interference exerted by C. clusii would imply less plant–plant associations, resulting in an impoverishment of species in its local vicinity compared to the other shrub. Since the effects of phytotoxic compounds can increase along plant life-span [29], we expected 5) a more evident depletion of perennial plants (especially at adult stage) than of annual plants in the local vicinity of C. clusii.


Study area

The study was conducted in Sierra de Alcubierre (41°41′N 0°32′W, municipality of Leciñena), in the Middle Ebro Valley, Zaragoza (NE Spain), one of the largest gypsum outcrops in Europe [41]. This area has a semi-arid Mediterranean climate with high continental influence. Average precipitation is 367 mm year−1, and average annual temperature is 14.5 °C (Zuera ‘Aspasa’ meteorological station, 1973–2012 period; source: Gobierno de Aragón, The landscape is characterised by low hills (480 m a.s.l. average) with mainly gypsiferous lithology, and flat valleys, most of which have been cultivated. In the gypsiferous hills, plant communities are composed predominantly of highly specialised flora (the gypsophytes Helianthemum squamatum (L.) Pers., Gypsophila struthium Loefl. ssp. hispanica (Willk.) G. López, Ononis tridentata L. and Lepidium subulatum L.) and some widespread Mediterranean shrub species, e.g., Rosmarinus officinalis L., Thymus vulgaris L. and C. clusii [42]. The vegetation structure is a scattered scrubland comprising large open areas interspersed with patches of vegetation. This unique habitat (Gypsophiletalia) has a high ecological value and is listed as a conservation priority in international directives [43].

Greenhouse experiment

A greenhouse seeding experiment was performed under controlled conditions to identify potential phytotoxic effects of C. clusii leaf and root aqueous extracts on the germination, early survival and growth of nine species. The choice of species included the most abundant perennials co-occurring with C. clusii at the study area [44]. Selected species were G. struthium ssp. hispanica (hereafter G. struthium), H. squamatum, Helianthemum syriacum (Jacq.) Dum. Cours., R. officinalis, T. vulgaris, Helichrysum stoechas (L.) Moench, Linum suffruticosum Orteg. ex Planch. and Stipa lagascae Roem. & Schult. Cistus clusii was also included in the experiment to test autotoxicity. Ripe fruits were collected at the study site from ten similar-sized individuals per test species. Most species seeds were collected in June 2015, except for G. struthium seeds, which were collected in September 2014 and R. officinalis seeds, which were collected in February 2015, matching the fructification peak respectively. Seeds were separated from the fruits, discarding any malformed seeds.

Aqueous solutions were prepared from leaves and roots of C. clusii and used as watering treatments in the experiment. Water was used as a solvent to simulate the leaching of phytotoxic compounds by rainfall. Solutions were prepared by ‘cool extraction’, soaking fresh plant material in distilled water for 24 h at room temperature in darkness [45]. For fresh material, we used leaves recently collected from natural communities and roots from plants grown for 3 months in a nursery. It was unfeasible to collect a sufficient amount of roots from natural communities given the difficulties encountered due to the deep taproots of C. clusii [46]. The seeds used to grow C. clusii in the nursery were collected from the same population as the collected leaves. The treatments were aqueous leaf extracts 1.5 g l−1 concentration (L), aqueous root extracts 0.025 g l−1 concentration (R), a mixture of both extracts 1.525 g l−1 concentration (RL) and water as the control (C). The water: leaf and water: root volumetric ratios were equivalent, and were within the range that occurs in natural conditions (see Additional file 1).

The seeding experiment was performed in July 2015 in a greenhouse maintained at 25 °C during the day and 15 °C during the night. Trays (60 × 40 × 20 cm) were filled with a peat-based substrate in which a known seed mixture was sown (9 species × 15 seeds per species in each tray). Seeding density was 0.06 seeds cm−2. Each of the four extract treatments (L, R, RL and C) had five replicates (trays). To assure germination, hard seeds were pretreated to break coat-imposed dormancy [47]; specifically, H. squamatum, H. syriacum and T. vulgaris seeds were mechanically scarified using sandpaper [48] and C. clusii seeds received a dry-heat shock at 100 °C for 5 min [49]. Before sowing, all seeds were soaked in distilled water for 20 h to stimulate germination. Extract treatments (L, R, RL and C) were applied twice a week by watering trays with one litre of the specific aqueous solution. To record potential effects of the extract treatments on the delay of germination and possible cumulative effects on seedling survival, the experiment was monitored once per week. Germinated seedlings were labeled with a toothpick indicating the date of emergence, and seedling survival was recorded throughout the experiment. To avoid any position effects, trays were randomized once a week. After 10 weeks, living seedlings were harvested and washed, and the below-ground and above-ground parts of each plant were separated and kept in individual paper bags. Plants were dried in an oven at 70 °C for 48 h and weighed using a 0.01 mg precision balance. Total dry biomass and the ratio of below-ground/above-ground biomass were used as growth estimators.

Vegetation survey

To test the potential interference exerted by C. clusii on neighbouring plants under natural conditions, a vegetation survey was conducted in May 2014, at the peak of vegetation growth. We surveyed the plants growing under the canopy of C. clusii, and also under the canopy of G. struthium for comparison purpose. The latter is a gypsophyte shrub which has a proved nurse role [50] and, to our knowledge, without any phytotoxic effects. Both shrubs have similar architecture, providing similar soil temperature and surface compaction under their canopies, which were improved, compared to open patches (see Additional file 2). To obtain comparable samples from the surrounding open patches, we surveyed the vegetation in paired areas placed in a random direction ≥ 50 cm away from each sampled target plant. Sampled areas were defined by circles matching the size of the area under the canopy of the paired target plant [51]. For both focal species, 25 sets of paired plant-open patches were sampled (n = 100 circles). All plants growing within the circles were recorded and identified to the species level. To assess the potentially differential phytotoxic effects along plant life-span, plants were categorised as either annual (short-lived) or perennial (long-lived), and within perennial, as either seedling or adult. For each category, we estimated richness (number of species present) and abundance (number of individuals present of all plant species) at each microsite: in open patches, under the canopy of C. clusii, and under the canopy of G. struthium.

Data analyses

In the greenhouse experiment, the effects of the extract treatments either on total germination, as well as germination delay or on seedling survival gradual decline were evaluated considering germination and survival rates through time. Differences among extract treatments in seed germination rate and seedling survival rate of the nine test species were evaluated using Cox proportional hazard models and, for data visualisation, Kaplan–Meier curves [52]. Pairwise comparisons among extract treatments were performed with Tukey’s post hoc tests. For each test species, differences in total biomass and the ratio of below-ground/above-ground biomass among extract treatments were tested using linear mixed models (LMMs) with the tray as a random factor. Time since germination and size of seedlings were strongly correlated; therefore, the number of weeks from germination to the end of the experiment was included as a covariate. To attain the assumption of normality, the dependent variables were log-transformed.

Differences in richness and abundance among microsites were analysed by fitting generalized linear models (GLMs) with the assumption of a Poisson error distribution and log link function. The size of the sampled area (circle area) was included as a continuous covariate because it might have influenced the number of plants recorded. When a significant effect of the microsite was found, Tukey’s post hoc tests were applied for pairwise comparisons. In addition, for each of the nine test species used in the greenhouse experiment, G-tests (log likelihood ratio tests) were implemented to compare the observed frequencies with the haphazardly expected frequencies at each microsite (in open patches, under the canopy of C. clusii, and under the canopy of G. struthium). The expected frequencies were estimated as the total observed frequency of each species multiplied by the proportion of the area occupied by each microsite.

All statistical analyses were performed using R software [53]. To fit Cox models and construct Kaplan–Meier curves for germination and survival, the ‘survival’ package was used [54]. To fit LMMs for growth, the ‘nlme’ package was used [55]. To fit GLMs for richness and abundance, the ‘stats’ package was used [53]. All pairwise comparisons were performed using the ‘multcomp’ package [56].


Greenhouse experiment

Cistus clusii aqueous extracts had a negative effect on seed germination rates of three of the nine test species (H. squamatum, H. stoechas and C. clusii; Fig. 1). Seed germination of H. squamatum was lower in the trays subjected to extracts than it was in the control trays and did not differ significantly among extract treatments. Seed germination of H. stoechas was lower in all extract treatments compared to the control treatment, being the lowest in the RL treatment. Seed germination of C. clusii was significantly lower in the R and the RL treatments than it was in the control treatment; however, germination rates did not differ significantly between the trays watered with leaf extracts (L) and the control trays. At the end of the experiment, < 10% of R. officinalis seeds had germinated (Fig. 1); therefore, this species was excluded from the survival and growth analyses.

Fig. 1

Kaplan-Meier curves representing seed germination rate over time (weeks from the experiment start) under each extract treatment: C (control), L (leaf extracts), R (root extracts) and RL (root and leaf extracts mixture) for each test species separately. Different letters represent statistically significant differences between extract treatments after Tukey’s post hoc tests (P < 0.05)

Survival rates differed significantly among extract treatments in two of the test species (H. syriacum and L. suffruticosum). Seedling survival of L. suffruticosum was lower in the R and the RL treatments compared to both the L treatment and the control treatment. Seedling survival of H. syriacum was the lowest in the control treatment and the RL treatment and the highest in the L treatment (Fig. 2).

Fig. 2

Kaplan-Meier curves representing seedling survival rate over time (weeks since germination) under each extract treatment: C (control), L (leaf extracts), R (root extracts) and RL (root and leaf extracts mixture) for each test species separately. Different letters represent statistically significant differences between extract treatments after Tukey’s post hoc tests (P < 0.05)

Neither total biomass nor the ratio of below-ground/above-ground biomass differed significantly among extract treatments for any of the test species (see Additional files 3, 4).

Vegetation survey

Microsite had a significant effect on richness and abundance of annual and perennial plants (Fig. 3; see Additional file 5). In all cases, richness and abundance were lower in open patches than under the canopies of both shrubs. Significantly fewer species of perennial adults were found under the canopy of C. clusii than under the canopy of G. struthium; however, richness under the shrubs did not differ significantly for seedlings of perennial species and annual plants. The abundances of annual plants and perennial adult plants were lower under the canopy of C. clusii than under the canopy of G. struthium; however, the abundance of perennial seedlings did not differ significantly between the two shrubs (Fig. 3; see Additional file 5).

Fig. 3

Mean richness and abundance of annuals, perennial seedlings and perennial adults among microsites: in open patches, under the canopy of C. clusii and under the canopy of G. struthium. Microsite effect was significant in GLMs in all cases (P < 0.001). Different letters indicate statistically significant differences between microsites after Tukey’s post hoc tests (P < 0.05)

Except for G. struthium, H. stoechas and S. lagascae, the seedlings of the test species were significantly less frequent in open patches than they were under the canopy of C. clusii. Among adults, H. syriacum and H. squamatum were less frequent in open patches than they were under the canopy of C. clusii. The other test species showed no difference in frequency between these two microsites (Table 1; see Additional file 6). Most test species were significantly less frequent in open patches than they were under the canopy of G. struthium, except seedlings of S. lagascae and adults of H. squamatum and R. officinalis (Table 1; see Additional file 6). Five of the test species were significantly less frequent under the canopy of C. clusii than they were under the canopy of G. struthium; specifically, seedlings of H. syriacum, H. stoechas and C. clusii, and adults of T. vulgaris, H. stoechas, S. lagascae and C. clusii. The observed frequency of C. clusii seedlings under the canopy of C. clusii was significantly higher than expected (Table 1; see Additional file 6).

Table 1 Summary of the pairwise comparisons of the G-test, indicating significant differences (P < 0.05) between observed frequencies (f) and expected frequencies (fe) of the nine test species at each microsite


We combined a controlled experiment with a field survey aiming to disentangle the potential chemical mechanisms of interference exerted by C. clusii on neighbouring plant species in gypsum plant communities of the Middle Ebro Valley. While the controlled experiment allowed us to isolate the phytotoxic effect of C. clusii root and leaf aqueous extracts on the early establishment of the test species, the field survey showed a more complex picture. In the field, chemical interference influences the net plant–plant interactions outcome together with facilitation and competition for resources [57]. Thus, a complementary assessment considering experimental and field effects of C. clusii on our test species can help unravel the relative relevance of its potential chemical interference compared to other types of interference, i.e., competition for resources (Table 2).

Table 2 Comparison of the potential effects exerted by C. clusii under experimental and field conditions

The greenhouse experiment confirmed our hypothesis that aqueous extracts of C. clusii affect the development of some species from gypsum plant communities. Supplementary chemical analyses of C. clusii tissues confirmed the presence of water-soluble terpenes and phenolic compounds with potential phytotoxic activity (see Additional file 7). We hypothesised that aqueous extracts from C. clusii would affect germination rates, survival rates and growth of the test species negatively. Even though it was not visible for all test species, our hypothesis was supported by experimental results for germination and survival rates. However, it was not evident for seedling biomass, manifesting that C. clusii aqueous extracts do not affect seedling growth, at least at the short-term. Chemically inhibition or retardation of germination and seedling survival decline may have ecological implications in the community. This may result in an advantage of low-competitive species over the affected species at early life stages [21, 24], likely causing a species shift in the community. Germination inhibition by phytotoxic compounds is a phenomenon widely reported by other studies in semi-arid communities [21, 58].

Diverse effects of aqueous extracts were found on germination and survival rates depending on the plant tissue tested in the experiment. Other studies found that leaves of allelopathic plants contained more water-soluble phytotoxic compounds than roots [59], likely resulting in stronger phytotoxic effects, as observed by Dorning and Cipollini [40] in an invasive shrub. Our additional chemical analyses confirmed that leaves from C. clusii contain more water-soluble potential phytotoxic compounds than its roots. We expected the leaf extracts to exert a stronger negative effect on the test species than root extracts. On the contrary, despite containing fewer compounds than did its leaves, solutions containing root extracts more often had negative effects on germination and seedling survival compared to pure leaf extracts. As predicted, roots combined with leaves was the most inhibiting treatment likely due to a synergic effect of the compounds contained in both plant tissues [60].

The outcomes of the field survey and the experiment denoted that chemical interference could explain why some species are less frequent in the local vicinity of C. clusii. This fact was especially evident for H. stoechas because the low number of individuals found under the canopy of C. clusii compared to those under the canopy of G. struthium paralleled the low germination of seeds treated with C. clusii extracts. These results confirm our hypothesis of the lessening of plant–plant associations due to the chemical interference exerted by C. clusii on neighbouring plants, deriving to an impoverishment of species around this shrub. It has already been evidenced in semi-arid plant communities that phytotoxic effects on neighbouring plants result in species-poor islands around the allelopathic plant [61]. Other negative interactions between C. clusii and neighbouring species in the field could not be corroborated in our seeding experiment because other factors beyond the releasing of phytotoxic compounds may influence plant establishment under natural conditions [62]. For example, in the field, fewer than the expected number of individuals of H. syriacum were found under the canopy of C. clusii; however, in the experiment, C. clusii aqueous extracts had a positive effect on H. syriacum seedling survival. Similarly, T. vulgaris and S. lagascae were found in low abundance under the canopy of C. clusii; however, the extract treatments did not significantly affect their performance. Negative interactions between those species and C. clusii are thus more likely caused by competition for space and resources rather than chemical interference [57].

The limitations of the greenhouse experiment could have led to an underestimation of the phytotoxic effects of C. clusii that may occur under natural conditions. For example, the effects of hydrophobic compounds present in C. clusii (e.g., β-pinene; [39]) were not tested in the experiment because they were not extracted in the aqueous solutions (see Additional file 4). Those compounds can be released to the environment by volatilisation and are potentially phytotoxic [58, 63]. Also, soil microorganisms are known to transform chemical compounds [64], which can increase phytotoxicity under natural conditions [65], dissimilar to the controlled conditions of the experiment. On the other hand, it should be noted that phytotoxic compounds can have a cumulative effect over a plant life-span [29] and our seeding experiment did not last long enough to detect potential long-term negative effects of C. clusii. Nevertheless, accordingly to our expectations, long-term negative effects were found in the field survey since C. clusii consistently harboured fewer species than G. struthium, mainly perennials at adult stages.

Both in the greenhouse experiment and in the field C. clusii did not show a negative effect on G. struthium, suggesting that this species tolerates C. clusii. In the field, R. officinalis also appeared to exhibit tolerance to C. clusii; however, this could not be confirmed experimentally because R. officinalis exhibited an overall very low germination rate. Those species may have adapted to the potential phytotoxic compounds of C. clusii because they frequently co-occur with this species in gypsum plant communities of the Middle Ebro Valley. Tolerance to ‘chemical neighbours’ is a well-studied co-evolutionary phenomenon that allows a species to coexist with phytotoxic plant species [61, 66, 67]. Moreover, although C. clusii extracts reduced seedling survival of H. squamatum and L. suffruticosum in the greenhouse experiment, this was not evident in natural conditions. Possibly, in nature, the net interaction outcome between C. clusii and those species tends to be neutral because of positive interactions. Cistus clusii has the potential of behaving as a nurse plant since the area under its canopy provides microsites similar to known nurse plants in this habitat (e.g., G. struthium; [50]) and may suit the establishment of C. clusii-tolerant species [4]. Indeed, there was a positive net effect of C. clusii when richness and plant abundance were compared to open patches. Positive effects of C. clusii on the establishment of other species compared to open patches had been documented in semi-arid plant communities before [32].

Few of the studies that investigated phytotoxicity in plants also evaluated autotoxicity [33]. In the greenhouse experiment, C. clusii aqueous extracts inhibited the germination of its own seeds, indicating a phytotoxic potential against itself. In natural conditions, more seedlings of C. clusii than expected were found under its canopy, probably because of high seed accumulation [68]; however, fewer than the expected numbers of adult C. clusii were found under its canopy, in agreement with an auto-inhibiting effect. The low establishment of adult C. clusii in the vicinity of C. clusii shrubs could have important implications for C. clusii population dynamics as auto-inhibition could lead to a strong reliance on other nurse species to establish under the highly restrictive conditions that occur in gypsum environments.

Despite ameliorating micro-environmental conditions under its canopy (see Additional file 2), and having a nurse role compared to open patches [32], C. clusii did not present such a positive role as the other shrub with similar architecture, suggesting interference with neighbouring plants in the community. We found potential phytotoxic compounds in C. clusii leaves and roots, and the associated vegetation showed a species-specific sensitivity to C. clusii. Among test species, there were possible C. clusii-vulnerable species (H. stoechas, H. syriacum, S. lagascae, and T. vulgaris), but also potential C. clusii-tolerant species (G. struthium, H. squamatum, L. suffruticosum, and R. officinalis). Species-specific phytotoxicity has been previously described [27] and might have important ecological implications for the dynamics of plant communities, by affecting the recruitment of some species and thereby their abundance in the community [69]. Besides C. clusii, other gypsovags that are very common in gypsum environments (e.g., T. vulgaris and R. officinalis) leach chemical compounds with known phytotoxic activity to their local environments [70, 71]. Phytotoxicity might be a mechanism that allows gypsovags to succeed in competition for resources with neighbouring plant species that may be better adapted to the harsh conditions in gypsum soils [72]. For example, phosphorous is scarce in gypsum hills [73, 74], and excluding other plants from the local vicinity might be a means of minimising local phosphorous depletion.


This study provides novel results of species-specific interference of C. clusii on other plant species. Phytotoxicity of C. clusii at least partly affects species richness in its local vicinity in gypsum plant communities in the Middle Ebro Valley. The importance of the role of the phytotoxicity of C. clusii in plant–plant interaction outcomes at the community level should be investigated in other Mediterranean plant communities.



aqueous leaf extract


aqueous root extract


aqueous root and leaf extract


control (water)


linear mixed models


generalized linear models


  1. 1.

    Gross N, Kunstler G, Liancourt P, De Bello F, Suding KN, Lavorel S. Linking individual response to biotic interactions with community structure: a trait-based framework. Funct Ecol. 2009;23:1167–78.

    Article  Google Scholar 

  2. 2.

    Cáceres Y, Llambí LD, Rada F. Shrubs as foundation species in a high tropical alpine ecosystem: a multi-scale analysis of plant spatial interactions. Plant Ecol Diver. 2015;8:147–61.

    Article  Google Scholar 

  3. 3.

    Soliveres S, Maestre FT. Plant–plant interactions, environmental gradients and plant diversity: a global synthesis of community-level studies. Persp Plant Ecol Evol Syst. 2014;16:154–63.

    Article  Google Scholar 

  4. 4.

    Callaway RM. Direct mechanisms for facilitation. Positive interactions and interdependence in plant communities. Netherlands: Springer; 2007. p. 15–116. Accessed 9 Feb 2016.

  5. 5.

    Cavieres LA, Badano EI. Do facilitative interactions increase species richness at the entire community level? J Ecol. 2009;97:1181–91.

    Article  Google Scholar 

  6. 6.

    Escudero A, Romao RL, de la Cruz M, Maestre FT. Spatial pattern and neighbour effects on Helianthemum squamatum seedlings in a Mediterranean gypsum community. J Veg Sci. 2005;16:383–90.

    Google Scholar 

  7. 7.

    Barberá GG, Navarro-Cano JA, Castillo VM. Seedling recruitment in a semi-arid steppe: the role of microsite and post-dispersal seed predation. J Arid Environ. 2006;67:701–14.

    Article  Google Scholar 

  8. 8.

    Miriti MN. Ontogenetic shift from facilitation to competition in a desert shrub. J Ecol. 2006;94:973–9.

    Article  Google Scholar 

  9. 9.

    Valiente-Banuet A, Verdú M. Temporal shifts from facilitation to competition occur between closely related taxa. J Ecol. 2008;96:489–94.

    Article  Google Scholar 

  10. 10.

    Schöb C, Prieto I, Armas C, Pugnaire FI. Consequences of facilitation: one plant’s benefit is another plant’s cost. Funct Ecol. 2014;28:500–8.

    Article  Google Scholar 

  11. 11.

    Armas C, Pugnaire FI. Ontogenetic shifts in interactions of two dominant shrub species in a semi-arid coastal sand dune system. J Veg Sci. 2009;20:535–46.

    Article  Google Scholar 

  12. 12.

    Inderjit Duke SO. Ecophysiological aspects of allelopathy. Planta. 2003;217:529–39.

    CAS  Article  Google Scholar 

  13. 13.

    Bertin C, Yang X, Weston LA. The role of root exudates and allelochemicals in the rhizosphere. Plant Soil. 2003;256:67–83.

    CAS  Article  Google Scholar 

  14. 14.

    Zhang C, Fu S. Allelopathic effects of leaf litter and live roots exudates of Eucalyptus species on crops. Allelop J. 2010;26:91–100.

    Google Scholar 

  15. 15.

    Tsubo M, Nishihara E, Nakamatsu K, Cheng Y, Shinoda M. Plant volatiles inhibit restoration of plant species communities in dry grassland. Basic Appl Ecol. 2012;13:76–84.

    Article  Google Scholar 

  16. 16.

    Ehlers BK, Charpentier A, Grøndahl E. An allelopathic plant facilitates species richness in the Mediterranean garrigue. J Ecol. 2014;102:176–85.

    Article  Google Scholar 

  17. 17.

    Callaway RM, Ridenour WM. Novel weapons: invasive success and the evolution of increased competitive ability. Front Ecol Environ. 2004;2:436–43.

    Article  Google Scholar 

  18. 18.

    Langenheim JH. Higher plant terpenoids: a phytocentric overview of their ecological roles. J Chem Ecol. 1994;20:1223–80.

    CAS  Article  Google Scholar 

  19. 19.

    Inderjit. Plant phenolics in allelopathy. Bot Rev. 1996;1:186–202.

    Article  Google Scholar 

  20. 20.

    Rice EL. Allelopathy. 2nd ed. New York: Academic Press; 1984.

    Google Scholar 

  21. 21.

    Escudero A, Albert MJ, Pita JM, Pérez-García F. Inhibitory effects of Artemisia herba-alba on the germination of the gypsophyte Helianthemum squamatum. Plant Ecol. 2000;148:71–80.

    Article  Google Scholar 

  22. 22.

    Ridenour W, Callaway R. The relative importance of allelopathy in interference: the effects of an invasive weed on a native bunchgrass. Oecologia. 2001;126:444–50.

    Article  Google Scholar 

  23. 23.

    Nielsen JA, Frew RD, Whigam PA, Callaway RM, Dickinson KJM. Germination and growth responses of co-occurring grass species to soil from under invasive Thymus vulgaris. Allelop J. 2015;35:139–52.

    Google Scholar 

  24. 24.

    Verdú M, Traveset A. Early emergence enhances plant fitness: a phylogenetically controlled meta-analysis. Ecology. 2005;86:1385–94.

    Article  Google Scholar 

  25. 25.

    De Almeida LFR, Frei F, Mancini E, De Martino L, De Feo V. Phytotoxic activities of mediterranean essential oils. Molecules. 2010;15:4309–23.

    Article  Google Scholar 

  26. 26.

    Dhima K, Vasilakoglou I, Garane V, Ritzoulis C, Lianopoulou V, Panou-Philotheou E. Competitiveness and essential oil phytotoxicity of seven annual aromatic plants. Weed Sci. 2010;58:457–65.

    CAS  Article  Google Scholar 

  27. 27.

    Linhart YB, Gauthier P, Keefover-Ring K, Thompson JD. Variable phytotoxic effects of Thymus vulgaris (Lamiaceae) terpenes on associated species. Int J Plant Sci. 2015;176:20–30.

    Article  Google Scholar 

  28. 28.

    Pedrol N, González L, Reigosa MJ. Allelopathy and abiotic stress. Berlin: Springer; 2006. p. 171–209.

    Google Scholar 

  29. 29.

    Sosa T, Valares C, Alías JC, Lobón NC. Persistence of flavonoids in Cistus ladanifer soils. Plant Soil. 2010;337:51–63.

    CAS  Article  Google Scholar 

  30. 30.

    Demoly J-P, Montserrat P. Cistus L. Castroviejo, S (eds) Flora Ibérica. Madrid: Real Jardín Botánico, CSIC; 1990. p. 319–37.

  31. 31.

    Munné-Bosch S, Jubany-Marí T, Alegre L. Enhanced photo- and antioxidative protection, and hydrogen peroxide accumulation in drought-stressed Cistus clusii and Cistus albidus plants. Tree Physiol. 2003;23:1–12.

    Article  Google Scholar 

  32. 32.

    Armas C, Pugnaire FI. Plant interactions govern population dynamics in a semi-arid plant community. J Ecol. 2005;93:978–89.

    Article  Google Scholar 

  33. 33.

    Alías JC, Sosa T, Escudero JC, Chaves N. Autotoxicity against germination and seedling emergence in Cistus ladanifer L. Plant Soil. 2006;282:327–32.

    Article  Google Scholar 

  34. 34.

    Chaves N, Sosa T, Alias JC, Escudero JC. Germination inhibition of herbs in Cistus ladanifer L. soils: possible involvement of allelochemicals. Allelop J. 2003;11:31–42.

    Google Scholar 

  35. 35.

    Chaves N, Escudero JC. Allelopathic effect of Cistus ladanifer on seed germination. Funct Ecol. 1997;11:432–40.

    Article  Google Scholar 

  36. 36.

    Castells E, Peñuelas J, Valentine DW. Are phenolic compounds released from the Mediterranean shrub Cistus albidus responsible for changes in N cycling in siliceous and calcareous soils? New Phytol. 2004;162:187–95.

    CAS  Article  Google Scholar 

  37. 37.

    Barrajón-Catalán E, Fernández-Arroyo S, Roldán C, Guillén E, Saura D, Segura-Carretero A, et al. A systematic study of the polyphenolic composition of aqueous extracts deriving from several Cistus genus species: evolutionary relationship. Phytochem Anal. 2011;22:303–12.

    Article  Google Scholar 

  38. 38.

    Hernández I, Alegre L, Munné-Bosch S. Drought-induced changes in flavonoids and other low molecular weight antioxidants in Cistus clusii grown under Mediterranean field conditions. Tree Physiol. 2004;24:1303–11.

    Article  Google Scholar 

  39. 39.

    Morales-Soto A, Oruna-Concha MJ, Elmore JS, Barrajón-Catalán E, Micol V, Roldán C, et al. Volatile profile of Spanish Cistus plants as sources of antimicrobials for industrial applications. Ind Crops Prod. 2015;74:425–33.

    CAS  Article  Google Scholar 

  40. 40.

    Dorning M, Cipollini D. Leaf and root extracts of the invasive shrub Lonicera maackii inhibit seed germination of three herbs with no autotoxic effects. Plant Ecol. 2006;184:287–96.

    Article  Google Scholar 

  41. 41.

    Machı́n J, Navas A. Spatial analysis of gypsiferous soils in the Zaragoza province (Spain), using GIS as an aid to conservation. Geoderma. 1998;87:57–66.

    Article  Google Scholar 

  42. 42.

    Braun-Blanquet J, Bolòs O. Les groupements végétaux du bassin moyen de l’Ebre et leur dynamisme. Anales de la Estación Experimental de Aula Dei. 1958;.

  43. 43.

    Anonymous. Council directive 92/43/EEC of 21 May 1992 on the conservation of natural habitats and of wild fauna and flora. European Community, Brussels, Belgium. 1992.

  44. 44.

    Pueyo Y, Alados CL, Barrantes O, Komac B, Rietkerk M. Differences in gypsum plant communities associated with habitat fragmentation and livestock grazing. Ecol Appl. 2008;18:954–64.

    CAS  Article  Google Scholar 

  45. 45.

    Herranz JM, Ferrandis P, Copete MA, Duro EM, Zalacaín A. Effect of allelopathic compounds produced by Cistus ladanifer on germination of 20 Mediterranean taxa. Plant Ecol. 2005;184:259–72.

    Article  Google Scholar 

  46. 46.

    Guerrero-Campo J. Respuestas de la vegetación y de la morfología de las plantas a la erosión del suelo: Valle del Ebro y Prepirineo aragonés. Consejo de Protección de la Naturaleza de Aragón; 1998.

  47. 47.

    Thanos CA, Georghiou K, Kadis C, Pantazi C. Cistaceae: a plant family with hard seeds. Isr J Bot. 1992;41:251–63.

    Google Scholar 

  48. 48.

    Pérez-García F, González-Benito ME. Seed germination of five Helianthemum species: effect of temperature and presowing treatments. J Arid Environ. 2006;65:688–93.

    Article  Google Scholar 

  49. 49.

    Castro J, Romero-García AT. Dormancy and germination in Cistus clusii (Cistaceae): effect of biotic and abiotic factors. Plant Ecol. 1999;144(1):103–14.

    Article  Google Scholar 

  50. 50.

    Navarro-Cano JA, Ferrer-Gallego PP, Laguna E, Ferrando I, Goberna M, Valiente-Banuet A, et al. Restoring phylogenetic diversity through facilitation. Restor Ecol. 2016;24:449–55.

    Article  Google Scholar 

  51. 51.

    Cavieres LA, Badano EI, Sierra-Almeida A, Gómez-González S, Molina-Montenegro MA. Positive interactions between alpine plant species and the nurse cushion plant Laretia acaulis do not increase with elevation in the Andes of central Chile. New Phytol. 2006;169:59–69.

    Article  Google Scholar 

  52. 52.

    Bewick V, Cheek L, Ball J. Statistics review 12: survival analysis. Crit Care. 2004;8:389–94.

    Article  Google Scholar 

  53. 53.

    R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria: ISBN 3-900051-07-0; 2014.

  54. 54.

    Therneau TM. A package for survival analysis in S. R Pack Ver. 2014;2:37.

    Google Scholar 

  55. 55.

    Pinheiro J, Bates D, DebRoy S, Sarkar D. Linear and nonlinear mixed effects models. R Pack Ver. 2007;3:57.

    Google Scholar 

  56. 56.

    Hothorn T, Bretz F, Westfall P. Simultaneous inference in general parametric models. Biom J. 2008;50:346–63.

    Article  Google Scholar 

  57. 57.

    Inderjit, Del Moral R. Is separating resource competition from allelopathy realistic? Bot Rev. 1997;63:221–30.

    Article  Google Scholar 

  58. 58.

    Arroyo AI, Pueyo Y, Pellissier F, Ramos J, Espinosa-Ruiz A, Millery A, et al. Phytotoxic effects of volatile and water soluble chemicals of Artemisia herba-alba. J Arid Environ. 2018;151:1–8.

    Article  Google Scholar 

  59. 59.

    Whitehead DC, Dibb H, Hartley RD. Bound phenolic compounds in water extracts of soils, plant roots and leaf litter. Soil Biol Biochem. 1983;15:133–6.

    CAS  Article  Google Scholar 

  60. 60.

    Kruse M, Strandberg M, Strandberg B. Ecological effects of allelopathic plants—a review. NERI Technical Report, No. 315; 2000. p. 11–64.

  61. 61.

    Arroyo AI, Pueyo Y, Saiz H, Alados CL. Plant–plant interactions as a mechanism structuring plant diversity in a Mediterranean semi-arid ecosystem. Ecol Evol. 2015;5:5305–17.

    Article  Google Scholar 

  62. 62.

    Inderjit Callaway RM. Experimental designs for the study of allelopathy. Plant Soil. 2003;256:1–11.

    CAS  Article  Google Scholar 

  63. 63.

    Chowhan N, Singh HP, Batish DR, Kohli RK. Phytotoxic effects of β-pinene on early growth and associated biochemical changes in rice. Acta Physiol Plant. 2011;33:2369–76.

    CAS  Article  Google Scholar 

  64. 64.

    Jilani G, Mahmood S, Chaudhry AN, Hassan I, Akram M. Allelochemicals: sources, toxicity and microbial transformation in soil—a review. Ann Microbiol. 2008;58:351–7.

    CAS  Article  Google Scholar 

  65. 65.

    Gagliardo RW, Chilton WS. Soil transformation of 2(3H)-benzoxazolone of rye into phytotoxic 2-amino-3H-phenoxazin-3-one. J Chem Ecol. 1992;18:1683–91.

    CAS  Article  Google Scholar 

  66. 66.

    Grøndahl E, Ehlers BK. Local adaptation to biotic factors: reciprocal transplants of four species associated with aromatic Thymus pulegioides and T. serpyllum. J Ecol. 2008;96:981–92.

    Article  Google Scholar 

  67. 67.

    Jensen CG, Ehlers BK. Genetic variation for sensitivity to a thyme monoterpene in associated plant species. Oecologia. 2010;162:1017–25.

    Article  Google Scholar 

  68. 68.

    Callaway RM. Indirect mechanisms for facilitation. Positive interactions and interdependence in plant communities. Netherlands :Springer; 2007. p. 117–77. Available from: Accessed 27 Apr 2017.

  69. 69.

    Schenk HJ. Root competition: beyond resource depletion. J Ecol. 2006;94:725–39.

    Article  Google Scholar 

  70. 70.

    Thompson JD, Chalchat JC, Michet A, Linhart YB, Ehlers B. Qualitative and quantitative variation in monoterpene co-occurrence and composition in the essential oil of Thymus vulgaris chemotypes. J Chem Ecol. 2003;29:859–80.

    CAS  Article  Google Scholar 

  71. 71.

    Vokou D, Vareltzidou S, Katinakis P. Effects of aromatic plants on potato storage: sprout suppression and antimicrobial activity. Agr Ecosyst Environ. 1993;47:223–35.

    Article  Google Scholar 

  72. 72.

    Escudero A, Palacio S, Maestre FT, Luzuriaga AL. Plant life on gypsum: a review of its multiple facets. Biol Rev. 2015;90:1–18.

    Article  Google Scholar 

  73. 73.

    Guerrero-Campo J, Alberto F, Hodgson J, Garcı́a-Ruiz JM, Montserrat-Martı́ G. Plant community patterns in a gypsum area of NE Spain. I. Interactions with topographic factors and soil erosion. J Arid Environ. 1999;41:401–10.

    Article  Google Scholar 

  74. 74.

    Pueyo Y, Alados CL, Maestro M, Komac B. Gypsophile vegetation patterns under a range of soil properties induced by topographical position. Plant Ecol. 2007;189:301–11.

    Article  Google Scholar 

Download references

Authors’ contributions

All authors conceived the study; AF and BKE designed the greenhouse experiment; AF, YP and CLA designed the field surveys; AF performed the experiment, the chemical analyses and the field surveys; AF, BKE and YP designed the manuscript; AF analyzed the data and wrote the draft; BKE, YP and CLA commented on the manuscript. All authors read and approved the final manuscript.

Authors’ information

AF is a Ph.D. candidate in ecology, studying plant–plant interactions in semi-arid gypsum plant communities. She is experienced in the ecological restoration of gypsum environments degraded by opencast mining. BKE is an evolutionary ecologist with research interests in understanding how genetic variation within dominant plant species affect the ecology and evolution of associated species and how the local environments and intraspecific variation affect competition, facilitation and co-existence. CLA is an ecologist interested in the effects of environmental stress on ecosystem conservation, with special focus on Mediterranean and mountain grazing ecosystems. YP is a plant ecologist with research interests in the functioning of semiarid ecosystems, including plant–plant and plant-soil interactions and eco-hydrology, focusing on understanding the effects of environmental and human influences on plant communities for ecosystem conservation and restoration.


We thank L. Lauridsen and S. Slotsbo for their help with the greenhouse experiment and the supplementary chemical analyses, M. L. Giner, P. Sánchez, P. Bravo, J. Rodríguez and C. Capistrós for their help in the fieldwork, M. Ballesteros and A. I. Arroyo for the statistical advice and B. MacWhirter for the English proofreading.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its additional files (see Additional files 4, 5 and 6).

Consent for publication

Not applicable.

Ethics approval and consent to participate

Collections of specimens used in the study complied with regional and national guidelines. No permissions were requested because the area where samples were taken is not under any protective figure.


This work was supported by the Spanish Ministry of Economy, Industry and Competitiveness (CGL2012-37508, CGL2016-80783-R). Funding provided by CGL2012-37508, CGL2016-80783-R research projects was used to pay the expenses derived from the field surveys and the laboratory analyses. This ministry also supported A. Foronda with a Ph.D. Grant (BES-2013-063852), favouring the data analyses and interpretation and the writing of the manuscript. A. Foronda also received a Grant from this ministry for a short research stay in Denmark (EEBB-I-15-10298), which allowed the collaboration with B. K. Ehlers and the performance of the greenhouse experiment and the laboratory analyses.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Ana Foronda.

Additional files

Additional file 1.

Methods and results for the determination of aqueous extracts concentration.

Additional file 2.

Methods and results for the micro-environmental conditions under the canopies of C. clusii and G. struthium.

Additional file 3.

Effects of the extract treatments on seedling biomass.

Additional file 4.

Dataset including the data about seed germination, seedling survival and seedling biomass from the greenhouse experiment.

Additional file 5.

Dataset including the data about richness and abundance from the field surveys.

Additional file 6.

Dataset including the data about the observed and expected frequencies of the test species under field conditions.

Additional file 7.

Methods and results for the chemical analyses of leaf and root extracts.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Foronda, A., Ehlers, B.K., Alados, C.L. et al. Species-specific interference exerted by the shrub Cistus clusii Dunal in a semi-arid Mediterranean gypsum plant community. BMC Ecol 18, 49 (2018).

Download citation


  • Aqueous extracts
  • Germination
  • Interference
  • Nurse plant
  • Survival
  • Richness