Genetic variation for sensitivity to a thyme monoterpene in associated plant species.

Recent studies have shown that plant allelochemicals can have profound effects on the performance of associated species, such that plants with a history of co-existence with "chemical neighbour" plants perform better in their presence compared to naïve plants. This has cast new light on the complexity of plant-plant interactions and plant communities and has led to debates on whether plant communities are more co-evolved than traditionally thought. In order to determine whether plants may indeed evolve in response to other plants' allelochemicals it is crucial to determine the presence of genetic variation for performance under the influence of specific allelochemicals and show that natural selection indeed operates on this variation. We studied the effect of the monoterpene carvacrol-a dominant compound in the essential oil of Thymus pulegioides-on three associated plant species originating from sites where thyme is either present or absent. We found the presence of genetic variation in both naïve and experienced populations for performance under the influence of the allelochemical but the response varied among naïve and experienced plant. Plants from experienced populations performed better than naïve plants on carvacrol soil and contained significantly more seed families with an adaptive response to carvacrol than naïve populations. This suggests that the presence of T. pulegioides can act as a selective agent on associated species, by favouring genotypes which perform best in the presence of its allelochemicals. The response to the thyme allelochemical varied from negative to neutral to positive among the species. The different responses within a species suggest that plant-plant interactions can evolve; this has implications for community dynamics and stability.

Introduction

Many plant species produce secondary compounds with allelopathic effects (reviews in Rice 1979; Langenheim 1994). Some of these plants have the ability, via leaching of allelochemicals, to modify their local environment in a way that puts associated species under selective pressure to adapt to the direct as well as indirect effects of the allelochemicals. In trees, conifers are known for their high production and diversity of terpenes and phenolics which significantly modify the soil conditions (White 1991; Pavolainen et al. 1998; Prescott et al. 2000) and affect the diversity of the understory vegetation (Iason et al. 2005). Some invasive species can nearly eradicate entire native plant communities via the modifications that come about from the leaching of their allelochemicals. For instance, Centaurea maculosa and Centaurea diffusa are invasive species in North America. Both of these species leach specific allelochemicals which negatively impact the performance of other plant species and modify the local soil biota creating a positive feedback on themselves. The negative effects of the allelochemicals are much greater in North American communities compared to native Eurasian ones (e.g. Callaway and Aschehoug 2000; Callaway et al. 2004; see also Blair et al. 2008 for discussions on detection and stability of specific allelochemicals in soil). Prati and Bossdorf (2004) showed similar responses in their study on the allelopathic effects of the invasive plant, Alliaria petiolata. The high production of phenolic compounds in Vaccinium myrtillus differentially affects the germination and growth of two spruce species, having a more negative effect on the introduced rather than native spruce species (Mallik and Pellissier 2000). These studies demonstrate that sensitivity to specific plant allelochemicals can vary with the history of co-existence of the interacting plants.

Terpenes are the largest class of secondary compounds in plants (Langenheim 1994; Gershenzon and Dudareva 2007) and their ecological role is a key component in both antagonistic and mutualistic interactions. Their functions vary from having either inhibitory or stimulating effects on a number of associated organisms such as bacteria, fungi, other plants and herbivores (Langenheim 1994; Gershenzon and Dudareva 2007). In a community context, the role of terpenes in soil nutrient cycling may be particularly important. It is known that terpenes in the soil can alter the nitrogen cycle by inhibiting the nitrification process, thereby decreasing the amount of available nitrate in the soil, and increasing that of ammonium (White 1991; Pavolainen et al. 1998; Hättenschwiler and Vitousek 2000; Castells et al. 2003). Presence of terpenes in the soil may thereby affect the competitive interactions between plants in the community by favouring plants that are good at adjusting or adapting to the local changes in nitrogen availability.

One plant genus famous for its production of terpenes is the aromatic Thymus. The common Mediterranean species Thymus vulgaris has a genetically determined polymorphism for the production of monoterpenes constituting the main component of its essential oil. A total of at least seven different so-called chemotypes, each dominated by a different monoterpene, are known from this species (Vernet et al. 1986). Thyme populations differ with respect to chemotypes produced, and individual populations are usually dominated by one to two different types (Thompson et al. 2003). The grass Bromus erectus shows local adaptation to its home non-phenolic chemotype by germinating and growing better on soil collected under its “home” thyme chemotype compared to soil from sites where T. vulgaris produces a phenolic chemotype (Ehlers and Thompson 2004). Thyme species found in northern Europe also show variation in their terpene production (Sáez and Stahl-Biskup 2002). Grøndahl and Ehlers (2008) studied four different plant species originating from sites dominated by either Thymus pulegioides (producing a phenolic monoterpene), or T. serpyllum (producing sesquiterpenes) and showed that plants germinated and grew better on soil treated with their “home” terpene. Plants that showed an adaptive response to their local thyme terpene originated from sites where associated plants and thyme showed a significantly positive spatial association.

The above studies demonstrate that plant allelochemicals can act as a strong selective force driving evolution in plant–plant interactions. Fundamental for local adaptation to evolve is the presence of genetic variation for an adaptive response upon which selection can operate. So far studies of associated plant species’ response to the presence of a “chemical neighbour” have focused on a “population-level” response by comparing the performance of species from different origins and history with the plant producing the allelochemical. Variation in performance and sensitivity among individual plants to a plant allelochemical has been documented (e.g. Callaway et al. 2005a). However, as far as we know, no one has yet examined the presence of genetic variation in associated species for sensitivity to specific plant allelochemicals, and asked if such variation is in accordance with selection operating on this.

In this study we looked for the presence of genetic variation for sensitivity to the thyme monoterpene carvacrol, a dominant component of the essential oil in T. pulegioides. We use maternal seed families from three different species originating from sites where they either grow together with T. pulegioides (experienced) or from sites where T. pulegioides is absent (naïve). All three study species are common in plant communities where T. pulegioides dominates the ground vegetation. We ask the following questions:

Do associated plant species show genetic variation for sensitivity to the thyme monoterpene?

Do experienced and naïve plants differ in their response to carvacrol?

Is there variation in sensitivity to the thyme monoterpene between experienced and naïve plants, which may then suggest that the presence of T. pulegioides in a plant community can act as a selective agent on the genetic variation in associated species?

Materials and methods

Study sites and study species

Five sites in Denmark, ranging from Eastern Jutland to the western part of Funen were chosen for this study. At two of the sites Thymus pulegioides dominates the ground cover of the vegetation and at the other sites T. pulegioides is absent. At sites with T. pulegioides, the essential oil of the thyme plants is dominated by the phenolic monoterpene carvacrol (Grøndahl et al. 2008). All sites are dry grassland habitat, and located in the same geographic region around the Eastern Sea.

We chose three plant species, Plantago lanceolata L. (Plantaginaceae), Campanula rotundifolia L. (Campanulaceae) and Agrostis capillaris L. (Poaceae), common at the T. pulegioides sites. All three plant species are short lived, outcrossing perennials, mainly wind pollinated (Agrostis and Plantago) or insect pollinated (Campanula), with seeds dispersed by wind (Agrostis and Campanula) or animals (Plantago). Seeds were collected from sites where thyme was either a dominant component of the ground vegetation or from sites where thyme was absent. We refer to these as “experienced” and “naïve”, respectively, with reference to thyme chemistry. We collected seeds from two naïve and two experienced populations of P. lanceolata, and from one naïve and one experienced population of both C. rotundifolia and A. capillaris. Seeds were collected from up to 20 different maternal plants per species per population at fruit ripening the year prior to the experiment and kept in paper bags until sowing. At the sites where T. pulegioides is present, seeds were collected from plants growing in the middle of thyme tufts. For both experienced and naïve populations, seeds were collected from maternal plants with an interplant distance of at least 2 m to avoid sampling of seed families from genetically identical plants.

Soil preparation

Two types of soil were prepared; soil treated with the thyme terpene carvacrol and control soil. The soil used was standard greenhouse soil (Pindstrup no. 1) with a low fertilisation level (approximately 0.078 g nitrogen, 0.039 g phosphorus and 0.13 g potassium l−1 soil). Soil with a low fertilisation level was chosen to mimic the general low fertilisation level in dry grassland habitats from which the study species originated. Soil treated with carvacrol contained 50 μl carvacrol 100 g−1 soil dry weight. This corresponds to approximately 500 μg g−1 and is within the range of concentration plants can experience in natural sites [230–1,500 μg g−1 dry weight, see Grøndahl and Ehlers (2008) for details].

The terpene was added to the soil as follows: liquid carvacrol (Sigma-Aldrich) was mixed in separate Petri dishes with filter paper; the filter paper was cut into small pieces of approximately 1 cm2 and sealed in plastic bags for 24 h after which all liquid had soaked into the paper; the filter papers were mixed carefully into the soil. Soil with the filter paper was then left for a further 24 h sealed in plastic bags to give a homogeneous terpene concentration and was afterwards transferred into individual germination pots.

Experimental setup

Seeds from each maternal plant were germinated in control and carvacrol soil. After an initial growth period the seedlings were transferred to individual pots with the respective soil treatment.

The number of seeds and seedlings available per plant varied. In order to test for any among-family variation in response to soil treatment we chose to have a minimum of four and a maximum of eight plants/maternal plant per soil treatment. This resulted in a final set up of 12 seed families per population from the two naïve and two experienced populations of P. lanceolata, 17 naïve and 20 experienced seed families from C. rotundifolia and 20 naïve and 18 experienced seed families from A. capillaris. Individual pots were kept in a greenhouse with day temperatures reaching up to 23°C and night temperatures were 15°C. Pots were randomised twice a week. After a growth period of 4 months all plants were harvested. The plants were rinsed and divided into above- and belowground parts and for Plantago and Campanula, reproductive stalks were separated from the vegetative parts (Agrostis did not flower during the course of the experiment). The number of reproductive stalks was counted and the vegetative plant material was oven dried at 50°C for 48 h and then weighed. Biomass was assessed from a total of 620 plants of P. lanceolata, 504 plants of C. rotundifolia and 550 plants of A. capillaris.

Statistical analysis

Mixed model ANOVA was used to analyse the effect of treatment, population, family (within population) and their interaction on biomass of Agrostis and Campanula plants, with treatment considered a fixed effect and the remaining predictor variables as random. For Plantago, as we had two experienced and two naïve populations, the mixed model ANOVA additionally included the effect of origin (either experienced or naïve with respect to thyme), population nested within origin, and family nested within population nested within origin. Residual error plots were inspected visually to verify normality. These analyses were performed using the program JMP (SAS Institute 1999).

To assess the reaction norm of the seed families from naïve and experienced populations, their respective response to treatment was analysed using a maximum likelihood method. We refer to a “positive reaction norm” when a seed family performed better on carvacrol soil than on control and a “negative reaction norm” for a seed family with the reverse response. Each population was characterised by the number of reaction norms in each category. Assuming that the proportion of positive reaction norms is binomially distributed, we tested for differences in the proportion of positive reaction norms between experienced and naïve populations using a likelihood ratio test (LRT). Significance was assessed by assuming that under the null hypothesis of equal proportions, the LRT statistic is X 2 distributed with 1 df.

Results

Effect of treatments and origin on plant performance

Our results show that the effect of carvacrol depends on the origin of plants. This was demonstrated by significant population × treatment interactions for total, aboveground and root biomass of Agrostis, for aboveground biomass of Campanula and by a significant origin × treatment interaction for Plantago (Table 1).

Table 1

Summary of ANOVA for biomass of Agrostis capillaris, Campanula rotundifolia, and Plantago lanceolata, and spike production of C. rotundifolia and P. lanceolata in response to origin (Ori), Thymus pulegioides present or absent, soil treatments (Trt), Population (Pop) and family (Fam)

Agrostis capillaris	Total biomass R 2 = 0.42			Aboveground biomass R 2 = 0.45		Root biomass R 2 = 0.34		Number of spikes
Source	df	SS	F	SS	F	SS	F
Trt	1	0.11	0.03	17.6	15.02***	14.9	9.31**	na
Pop	1	306.3	73.4***	104.3	89.04***	53.1	33.11***	na
Fam (Pop)	36	713.6	4.75***	215.0	5.10***	212.0	3.67***	na
Trt × Pop	1	56.4	13.52***	9.5	8.12**	19.6	12.21***	na
Trt × Fam (Pop)	36	342.7	2.28***	82.5	1.95**	120.6	2.09**	na
Error	462	1,927.9		541.1		741.5		na

* p < 0.05

** p < 0.01

*** p < 0.001

For Agrostis (Table 1; Fig. 1a,d,g), the biomass of experienced plants was overall higher than that of naïve plants, with the largest difference between naïve and experienced plants on carvacrol soil. Moreover, naïve Agrostis plants had the lowest total and aboveground biomass on carvacrol soil, whereas experienced plants increased their total and root biomass on carvacrol soil compared to control soil (Fig. 1a,g).

Fig. 1

Total, aboveground, and root biomass of Agrostis capillaris, Campanula rotundifolia, and Plantago lanceolata and numbers of spikes produces by C. rotundifolia and P. lanceolata. Plants grew on control soil (Control) and on soil treated with the Thymus pulegioides monoterpene carvacrol (Carvacrol). White bars Plant originating from naïve populations, grey bars plant originating from experienced populations

Experienced Campanula plants also had significantly higher aboveground biomass than naïve plants on carvacrol soil (Fig. 1b,e), whereas no difference in biomass between experienced and naïve plants was found on control soil. While naive Campanula plants had the lowest biomass on carvacrol soil, the experienced plants performed equally well on both soil treatments. The number of spikes produced did not differ between the naïve and experienced population but varied significantly between treatments, always being lower on carvacrol soil (Table 1; Fig. 1j).

In contrast to Agrostis and Campanula plants, the biomass of both experienced and naïve Plantago was significantly lower on carvacrol soil compared to control soil (Fig. 1c,f,i). However, a significant treatment × origin effect for aboveground biomass showed that experienced Plantago performed better on carvacrol soil than naïve Plantago, whereas no difference in aboveground biomass was found between experienced and naïve plants on control soil (Table 1; Fig. 1f). The number of spikes produced varied both with origin and populations (nested within origin; Table 1). Plantago from both experienced populations produced, overall, more spikes than plants from naïve sites (Fig. 1k). In particular, one experienced population produced almost twice as many spikes on both soil treatments as those from the other sites (data not shown).

Effect of treatments and family on plant performance

The significant interactions between treatment and family (nested within population; Table 1) for the total, aboveground, and root biomass of Agrostis suggest the presence of genetic variation for sensitivity to carvacrol. Combined with the significant population × treatment interactions, these show that naïve and experienced seed families differ in their response to the thyme monoterpene. Indeed reaction norms of naïve and experienced seed families show that experienced populations harbour significantly more seed families with a positive growth response to carvacrol than naïve populations (Fig. 2; Table 2).

Fig. 2

Reaction norm of individual seed families from naïve (a–c) and experienced (d–f) populations of A. capillaris, C. rotundifolia, and P. lanceolata. Experienced populations had been collected from sites with carvacrol-producing T. pulegioides. Each line represents the mean total biomass of plants from a seed family grown on control soil (Control) and on carvacrol soil (Carvacrol). Solid lines Positive reaction norm, dotted lines negative reaction norms

Table 2

Summary of likelihood ratio test (LRT) for differences between experienced and naïve populations in proportion of seed families showing a positive growth response on carvacrol soil

Species	Naive (N)/Experienced (E)	n	Positive reaction norm	Negative reaction norm	LRTa (P-value)
Total biomass
Plantago	N	24	1	23	4.57
Plantago	E	24	6	18	(P = 0.032)
Campanula	N	17	5	12	1.63
Campanula	E	20	10	10	(P = 0.2)
Agrostis	N	20	3	17	3.41
Agrostis	E	18	9	9	(P = 0.0009)
Pooled	N	61	8	53	8.6
Pooled	E	62	22	40	(P = 0.003)
Aboveground biomass
Plantago	N	24	1	23	4.57
Plantago	E	24	6	18	(P = 0.032)
Campanula	N	17	4	13	3.88
Campanula	E	20	11	9	(P = 0.048)
Agrostis	N	20	2	18	3.19
Agrostis	E	18	6	12	(P = 0.07)
Pooled	N	61	7	54	11.41
Pooled	E	62	23	39	(P = 0.0007)
Root biomass
Plantago	N	24	1	23	0
Plantago	E	24	1	23	(P = 1)
Campanula	N	17	6	11	0.1
Campanula	E	20	8	12	(P = 0.77)
Agrostis	N	20	9	11	7.82
Agrostis	E	18	15	3	(P = 0.005)
Pooled	N	61	16	45	2.19
Pooled	E	62	24	38	(P = 0.14)

aLRT null hypothesis assumes equal proportions of positive reaction norms between populations. P-values for LRT are χ2 distributed with 1 df

Plantago also showed significant treatment × family interactions for total, aboveground, and root biomass, and for the number of spikes produced, suggesting genetic variation for sensitivity to the thyme monoterpene. For aboveground biomass which also showed a treatment × origin effect, this suggests that genetic variation in sensitivity to carvacrol differed between naïve and experienced populations. The number of families with a positive reaction norm in experienced versus naïve populations shows that experienced Plantago populations have significantly more seed families with positive reaction norms for both total and aboveground biomass compared to naive seed families (Fig. 2; Table 2).

We did not detect any treatment × family interactions in Campanula (Table 1) but we found that the experienced population had more seed families with a positive reaction norm than naïve families (Fig. 2; Table 2), though this was only significant for the aboveground biomass (Table 2).

In no instance did the number of families with a positive reaction norm from naïve populations exceed that of the experienced ones.

Discussion

In this study we document genetic variation for sensitivity to the thyme monoterpene carvacrol in three different plant species. We show a clear tendency for plants experienced with T. pulegioides to be less affected, and in the case of Agrostis, even to be positively affected, by the thyme terpene compared to plants naïve to thyme. We find especially that Agrostis—and to a lesser extent Campanula and Plantago—show patterns of genetic variation in sensitivity to the thyme monoterpene, which suggests that selection is operating on this. For all response variables in Agrostis, for aboveground biomass in both Campanula and Plantago, and for the latter also total biomass, we find that populations experienced with thyme have significantly more seed families with a positive growth response to the thyme terpene compared to populations naïve to thyme. This finding, and the fact that we in no instance found that naïve populations exceeded experienced ones in number of seed families with a positive growth response to carvacrol, is suggestive of a “footprint” of selection in the experienced populations for adaptation to the effects of the monoterpene leached by their thyme neighbour.

The presence of genetic variation in response to environmental cues is necessary for the ability of populations to evolve adaptations to the local environment. While we demonstrate that associated species show genetic variation for sensitivity to a local thyme monoterpene and that this variation is suggestive of selection operating at sites where thyme is present, the sampling of only one to two experienced and naïve populations of each plant species limits our ability to generalise on the adaptation of associated plant species to thyme allelochemicals. The observed variation in response to the thyme terpene between naïve and experienced seed families could be caused by genetic differences among populations due to different site condition not related to the presence or absence of T. pulegioides. Even though all populations were sampled in dry grassland habitats and within the same climatic region, small-scale differences in environmental conditions and reduced migration among populations could create population genetic differentiation. However, several findings suggest that the presence of the thyme terpene carvacrol is the main reason for a differential response. First, the fact that the only difference among soil treatments is the presence/absence of carvacrol itself strongly points to the thyme terpene as the main reason for the differences found between naïve and experienced plants. Second, for both Campanula and Plantago we found no difference in biomass among naïve and experienced plants on control soil which we would expect if either maternal effects or genetic differences among populations irrespective of thyme were important for our results. Differences in biomass between experienced and naïve Campanula and Plantago were only found on carvacrol soil. Third, for Agrostis, experienced plants performed better on both types of soil treatment; however, the difference in response to carvacrol between naïve and experienced sites was strongly suggestive of an adaptation to the specific thyme terpene. Naïve Agrostis had the lowest total and aboveground biomass on carvacrol soil compared to control soil while experienced Agrostis had the highest biomass on carvacrol soil compared to control soil (Fig. 1a)—a “reciprocal” response pattern expected under local adaptation (e.g. Kawecki and Ebert 2004)

Experienced plants grow better on carvacrol soil than naïve pants; however, the response to the thyme terpene varies among species from negative (Plantago), to negative and neutral (Campanula naïve versus experienced), and from neutral to positive (naïve versus experienced Agrostis). This difference in response between naïve and experienced plants suggests that the outcome of an interaction between a plant allelochemical and its associated species can evolve.

An increasing number of studies show that plant interactions can act as drivers of local evolution, and this affects the way we view community structure and stability (e.g. review by Lortie et al. 2004 and Callaway et al. 2005b). Turkington and Harper (1979) and Martin and Harding (1981) demonstrated that competition among plant species varied among plant genotypes. The performance of individual plant species under inter-specific competition increased when competing with plant genotypes from a native site compared to competition with plant genotypes of the same species but originating from a distant site. Callaway and colleagues have shown how interactions among allelochemicals from invasive plants and associated plant species are strongly affected by the evolutionary history between plants and that the performance of plants is higher among experienced plants than plants naïve to the allelochemicals (e.g. Callaway and Aschehoug 2000; Callaway and Ridenour 2004; Callaway et al. 2005b). The present study supports these findings and highlights the importance of within-species variation in the ability to respond to biotic neighbours, for an understanding of how plant interactions evolve and how this in turn affects the stability and dynamics of plant communities.

How may thyme and its monoterpenes affect associated plants?

Thyme can be a dominant component of the vegetation in local plant communities and its effect on associated species is thus a combination of direct interaction between thyme plants and their neighbours as well as indirect interaction via the leaching of terpenes into the local environment. In this study we specifically address the effect of a thyme neighbour by studying the impact of its monoterpene on the performance of associated plant species. Monoterpenes are known to act directly as growth inhibitors of other plant species (Langenheim 1994; Vokou et al. 2003). Indirectly they may affect the performance of other plants via their effect on soil micro-organisms. Studies on aliphatic monoterpenes from pine have shown that they can affect both the amount and availability of nitrogen in soil. The amount of available nitrogen can be reduced by the effects of monoterpenes on heterotrophic organisms increasing the immobilisation of nitrogen and reducing nitrogen mineralization (White 1994). The availability of nitrogen is affected by monoterpenes inhibiting the nitrification process creating an increase in ammonium and a decrease in nitrate (White 1994; Pavolainen et al. 1998). However, the monoterpene concentrations at which these processes are affected vary. At high concentrations (>1,000 µg g−1 soil) net immobilization of ammonium occurred, whereas inhibition of the nitrification process occurred at monoterpene concentrations as low as 10-125 µg g−1 (White 1991, 1994). Similar effects may also be caused by phenolic monoterpenes in thyme. Analysis of soil collected under T. pulegioides plants and soil collected from the same site but where distance to a T. pulegiodes patch was at least 10 m, showed that the former had a significantly higher amount of ammonium and a higher ratio of ammonium:nitrate than soil collected away from thyme (B. K. Ehlers, unpublished result). In this study we used a monoterpene concentration of 500 µg g−1, i.e. within the range of what associated plants may experience under natural conditions (White 1994; Grøndahl and Ehlers 2008) and where the most likely indirect effects are the inhibition of nitrification. As nitrogen is the mineral which plants require in the largest amount (Bloom et al. 2003) and is also the nutrient which is most often limiting for plant growth in terrestrial ecosystems (Vitousek et al. 1997), the effects of monoterpenes on nitrogen cycling may be particularly important for associated plant species. Plants can absorb different sources of nitrogen but use mainly nitrate and ammonium. The assimilation of ammonium is energetically more favourable for the plant than nitrate but, unlike nitrate, ammonium cannot be stored in the plant and must be assimilated immediately (Crawford and Glass 1998). Moreover, high uptake of ammonium results in excretion of protons that acidify the rhizosphere and reduce the uptake of vital cations such as magnesium, potassium, and calcium. This effect (ammonium toxicity) varies among species. For instance, de Graaf et al. (1998) found that Arnica montana and Cirsium dissectum performed best using nitrate as their nitrogen source and had reduced survival and growth at increasing concentrations of ammonium, likely due to a decrease in cation uptake. In contrast, Calluna vulgaris grew best in soil that only contained ammonium compared to mixes of ammonium and nitrate or pure nitrate (de Graaf et al. 1998). Bloom et al. (2003) showed that Lycopersicon esculentum produced significantly heavier roots when growing with ammonium compared with nitrate as a nitrogen source, but without reduction in shoot mass, indicating that plant growth was not nutrient limited. This is similar to the response of the experienced Agrostis in the present study and could indicate a positive growth response to an increased ammonium concentration in carvacrol soil.

Terpenes of thyme may affect thyme-associated species differently depending on two things: their ability to respond and adapt to an alteration in soil nutrients, and their response and ability to adapt to the inhibiting and even toxic effects of monoterpenes on enzyme activity and biochemical processes of cells (e.g. White 1994) irrespective of their indirect effects on soil nutrients.

There is increasing interest in plant interactions as drivers of evolution in plant communities. This study adds to this field in two ways: we show that a dominant thyme monoterpene can act as a potential selective agent on associated plant species; and we demonstrate that where plants co-exist with thyme, associated plant species show genetic variation in sensitivity to this allelochemical in a way which is suggestive of selection operating.