Introduction

Plants rely on chemical cues to identify and mount appropriate responses to herbivores [1,2]. Among these cues are herbivore-induced plant volatiles (HIPVs) [3,4]. When perceived by undamaged plants, HIPVs can enhance plant defenses, and thus increase plant resistance to future herbivory [59]. HIPVs can directly induce plant defenses [7,8,10,11]. HIPVs can also prime plant defenses. In this case, HIPV-exposed plants respond more strongly and/or more rapidly to a secondary stimulus such as herbivore attack [5,7,8,12].

Plant responses to environmental cues depend on temporal patterns. For example, artificial damage that mimics spatiotemporal patterns of chewing herbivores induces a similar response to actual herbivore attack, whereas damage that ignores these patterns often does not [13]. The emission of HIPVs also follows discrete temporal dynamics [5,14], which, again, is likely to influence defense responses in receiver plants. Green leaf volatiles (GLVs), including (Z)-3-hexenyl acetate, (Z)-3-hexenal and (Z)-3-hexen-1-ol are catabolic products of the lipoxygenase/hydroperoxide lyase pathway and are emitted within minutes following damage [15,16]. Indole and terpenes are emitted within hours of herbivore attack [2,5]. Indole can mediate interactions between plants as well as between plants and other organisms [1719]. Terpenes are a very diverse group of volatiles (> 80,000 known structures) that are also important as both direct and indirect plant defenses [6,2022]. Indole and terpenes take longer to produce in comparison to GLVs and are less transiently emitted [2,5,15]. In maize, GLVs have been implicated in direct induction [8,10], and both GLVs and indole are known to prime defenses, and are thus considered as ‘bioactive’ in information transfer between plants [5,7].

Although it is clear that plants can either be primed or directly induced following HIPV exposure, if and how these two phenomena integrate in the context of a natural HIPV exposure sequence is unclear. Studies have yet to monitor induction and priming dynamics in real time and during continuous exposure to HIPVs. It is possible that, through time, continuous HIPV exposure may result in a self-enforced positive feedback loop, whereby priming enhances induction and induction serves as a priming event for subsequent induction.

Further, temporal trends might be complicated by dynamic environmental conditions, such as light fluctuations; over a 24 hr period, plants can experience both full sun and darkness. In darkness, the emission of many volatiles, especially those released via stomata, is severely hindered [2326].

In order to understand the kinetics of information transfer between plants and the importance of dynamic temporal emission patterns, we monitored the emission kinetics of terpenes as a marker of defense activation in maize plants that were exposed to HIPVs. We leveraged a highly temporally-resolved volatile multiplexing system that allowed us to track continuous volatile emissions in sender and receiver plants for the first time in real time. We further measured foliar terpene pools, stress hormones and defense gene expression in receiver plants, and determined the impact of day : night light transitions in the observed defense induction kinetics. Finally, we measured defenses upon sequences of short-term HIPV and GLV exposure to unravel the relative importance of direct induction and priming. Together, our experiments reveal how the natural kinetics of herbivore attack result in a clocked defense response in neighboring plants via volatile information transfer.

Results

Delayed burst in induced volatile emissions in plants exposed to volatiles of an herbivore-attacked neighbor

When maize plants are exposed to volatiles from herbivore-attacked neighbors, they start to release terpenes [7,11]. To understand the temporal dynamics of this phenomenon, we conducted a detailed time-course analysis of terpene release by receiver plants exposed to volatiles from sender plants under attack by Spodoptera exigua caterpillars (Fig 1). Herbivore feeding began around 12:00 and within 2 hr, we detected a small induction of sesquiterpenes [C H +, m/z = 205.20], monoterpenes [C H +, m/z = 137.13], 4,8-dimethylnona-1,3,7- triene (DMNT) [C11H19+, m/z = 151.15] and 4,8,12-trimethyltrideca-1,3,7,11-tetraene (TMTT) [C16H27+, m/z = 219.21] in undamaged receiver plants. Volatile emission was severely impaired during the night, regardless of treatment. Interestingly, as soon as the light was restored, we observed a strong burst of terpene release from receiver plants, which was 2-5 times higher than the release the day before (Fig 1). Thus, exposure to volatiles from an herbivore-attacked plant triggers a delayed burst in terpene emission in neighboring plants.

Delayed burst in induced volatile emissions in plants exposed to volatiles of a herbivore-attacked neighbor. Emission kinetics of herbivore-induced plant volatile (HIPV)-induced terpenes in undamaged receiver plants are shown. Dark green points represent mean emission of herbivore-damaged sender plants connected to undamaged receiver plants, with the emissions from damaged sender plants only subtracted. Black points represent undamaged sender plants connected to undamaged receiver plants, with the emissions from undamaged sender plants only subtracted. Blue rectangles represent the dark phase. Abbreviations. DMNT, 4,8-dimethylnona-1,3,7-triene; TMTT, 4,8,12-trimethyltrideca-1,3,7,11-tetraene. Error bars = SE. n = 8-10. Compounds were identified based on their molecular weight + 1, as all compounds were protonated. Sesquiterpenes: m/z = 205.20, Monoterpenes: m/z = 137.13, DMNT: m/z = 151.15, TMTT: m/z = 219.21.

The delayed burst in terpene emission is not explained by overaccumulation during the night

Why do receiver plants show a delayed burst in volatile-induced terpenoid release? During the night, plants close their stomata to limit water loss [27], which may also impair terpene release, thus leading to an accumulation of terpenes in the leaves and a burst once stomata open again [23,28]. To determine if terpenes accumulate above daytime levels in maize leaves of HIPV-exposed plants during the night, we measured internal foliar pools of the major sesquiterpenes, β-caryophyllene, β-farnesene, and α-bergamotene, as well as the homoterpene, TMTT over time. Internal pools of β-caryophyllene, β-farnesene and TMTT were marginally induced after 3 hours, but only became significantly induced after 8 hr of exposure to HIPVs (Fig 2A, 2B and 2D, Table 1). For all three, accumulation remained higher during the night (16.75 hr) as well as the following day. α-bergamotene only began accumulating in leaves on the second day (Fig 2C, Table 1). Thus, internal terpene pools remain comparable between night and day, suggesting that terpenes do not continue accumulating during the night and that, even when terpenes are emitted in large amounts on day two, internal pools do not drop below nighttime levels and that volatiles. Nevertheless, the volatiles that are present during the night, may be released suddenly the second day.

The delayed burst in terpene emission is not explained by terpene overaccumulation during the night. Accumulation of terpenes and induction of terpene biosynthesis genes in receiver plants was measured over time. A-D) Internal foliar pools of terpenes in receiver plants. E-H) Expression of terpene biosynthesis genes in receiver plants. Dark green bars represent receiver plants exposed to herbivore-induced plant volatiles (HIPVs) and light grey bars represent receiver plants connected to undamaged sender plants. Blue rectangles represent the dark phase. Abbreviations: TMTT, 4,8,12-trimethyltrideca-1,3,7,11-tetraene; FPPS3, Farnesyl pyrophosphate synthase 3; TPS2, Terpene synthase 2; TPS10, Terpene synthase 10; CYP92C5, Dimethylnonatriene/trimethyltetradecatetraene ynthase. . = p < 0.1, * = p < 0.05, ** = p < 0.01, *** = p < 0.001 as determined by a Welch’s two-sample t-test. Bars = mean ± SE. n = 4-6.

Welch’s two-sample t-test results comparing foliar terpene pools, biosynthesis genes and phytohormone levels between HIPV-exposed and control receiver plants. Bold values: p < 0.05, underlined values: p < 0.1. Abbreviations: β-car, β-caryophyllene; β-farn, β-farnesene; α-berg, α-bergamotene; TMTT, 4,8,12- trimethyltrideca-1,3,7,11-tetraene; FPPS3, Farnesyl pyrophosphate synthase 3; TPS2, Terpene synthase 2; TPS10, Terpene synthase 10; CYP92C5, Dimethylnonatriene/trimethyltetradecatetraene synthase; OPDA = 12- oxophytodienoic acid; JA = jasmonic acid; JA-Ile = jasmonic acid-isoleucine; OPR7 = oxo-phytodienoate reductase 7.

Some control plants accumulated measurable amounts of β-caryophyllene at night, likely due to biological and experimental variation. This type of variation was not visible in total sesquiterpene emissions measured by PTR-MS (Fig 1), likely because β-caryophyllene is a minor sesquiterpene in B73 [29] and closed stomata block sesquiterpene emission at night [23]. To get insight into terpene biosynthesis we measured the expression of terpene synthases in receiver plants. In maize, farnesene pyrophosphate synthase 3 (FPPS3), terpene synthase 2 and 10 (TPS2 and TPS10) as well as dimethylnonatriene/trimethyltetradecatetraene synthase (CYP92C5) are rate limiting for terpene production [3032]. All genes were induced by HIPV exposure during daytime, but not during the night. FPPS3 was slightly induced after 3 hr and significantly induced after 8 hr of HIPV exposure (Fig 2E, Table 1). TPS2 was slightly induced 3, 8 and 22 hr after the onset of HIPV exposure (Fig 2F, Table 1). TPS10 was significantly induced after 3 and 22 hr of HIPV exposure (Fig 2G, Table 1). CYP92C5 was only significantly induced after 3 hr of HIPV exposure (Fig 2H, Table 1). The patterns of volatile terpene biosynthesis do not support a scenario where the terpene burst on the second day is due to continued biosynthesis but lack of emission during the night.

The delayed burst in terpenoid emission is associated with clocked jasmonate production

Volatile release in maize is regulated by jasmonates [33,34]. To understand whether the delayed terpene burst is associated with jasmonate signaling, we measured the levels of 12- oxo-phytodienoic acid (OPDA), jasmonic acid (JA) and jasmonic acid-isoleucine (JA-Ile), as well as the expression of oxo-phytodienoate reductase 7 (OPR7), which is critical for JA biosynthesis [35]. We found a significant induction of OPDA, JA-Ile production and OPR7 expression at the beginning of the second day, 22 hr after the onset of volatile exposure (Fig 3, Table 1). Thus, jasmonate production is temporally aligned with the delayed terpene burst at the onset of the second day.

The delayed burst in terpene emission is associated with clocked jasmonate production. Foliar jasmonate concentrations (A-C) and jasmonate biosynthesis (OPR7; D) in receiver plants over time are shown. Dark green bars represent receiver plants exposed to herbivore-induced plant volatiles (HIPVs) and light grey bars represent receiver plants connected to undamaged sender plants. Blue rectangles represent the dark phase. Abbreviations. OPDA = 12-oxophytodienoic acid, JA = jasmonic acid, JA-Ile = jasmonic acid-isoleucine, OPR7 = oxo-phytodienoate reductase 7. . = p < 0.1, * = p < 0.05, ** = p < 0.01, *** = p < 0.001 as determined by a Welch’s two-sample t-test. Bars = mean ± SE. n = 3-6.

The delayed volatile burst is conserved under continuous light

To test whether the delayed burst in terpene volatiles is linked to light : dark transitions, we exposed maize plants to HIPVs under continuous light. Similar to exposure during a normal light regime, the largest burst of terpene emission occurred ca. 12-18 hr after the onset of volatile exposure (Fig 4), suggesting that the temporal delay in volatile emission in receiver plants is not dependent on light fluctuations, but is otherwise clocked (Fig 4).

The delayed volatile burst is conserved under continuous light. Emission kinetics of herbivore-induced plant volatiles (HIPV)-induced terpenes in undamaged receiver plants under continuous light are shown. Plants were grown under normal light conditions. Lights were left on continuously following the start of the treatment. Dark green points represent the mean emission of herbivore damaged sender plants connected to undamaged receiver plants with the emissions from damaged sender plants only subtracted. Black points represent undamaged sender plants connected to undamaged receiver plants, with the emissions from undamaged sender plants only subtracted. Yellow rectangles represent when the lights would typically be turned off. Abbreviations: DMNT, 4,8- dimethylnona-1,3,7-triene; TMTT, 4,8,12-trimethyltrideca-1,3,7,11- tetraene. Error bars = SE. n = 8-10. Compounds were identified based on their molecular weight + 1, as all compounds were protonated. Sesquiterpenes: m/z = 205.20, Monoterpenes: m/z = 137.13, DMNT: m/z = 151.15, TMTT: m/z = 219.21.

In order to control for potential differences in induction capacity at different times of day, we designed a second continuous light experiment. However, instead of starting the herbivory treatment at ca. 12:00 hr, we started it around 20:00 hr. Similarly, we observed an initial bump in terpene emission followed by a larger burst 10-12 hr after HIPV exposure (Supplemental Fig 1). In both continuous light experiments the ‘second day’ burst occurred somewhat earlier than under normal light conditions (Fig 4 and Fig 1). It is possible that stomata, which are closed in the dark, delay the onset of the terpene burst under dark : light transitions.

The delayed volatile burst cannot be fully explained by the emission kinetics of bioactive herbivory-induced volatiles

A simple explanation for the delayed volatile burst may be that the bioactive volatiles are more strongly emitted from sender plants at the onset of day two, thus triggering a stronger volatile response in the receiver plants at this time. To test this hypothesis, we analyzed volatile emission kinetics of S. exigua-infested plants over time. We focused our analysis on GLVs and indole, which are known to induce and/or prime volatile release in neighboring maize plants [7,8,10,11]. Terpenes are not known to prime or induce volatile release in maize, and were thus not included in the sender plant analysis [11]. We analyzed GLV and indole emission data of the sender plants from the experiments depicted in Fig 1 and Fig 4 (Fig 5), and then correlated their emission with the induction of terpenes in receiver plants (Fig. 6). We found strong positive correlations between hexenyl acetate, hexenal, hexen-1-ol and indole emissions in sender plants and terpene responses in receiver plants when plants had a dark period (Fig 6). Interestingly, this was not the case under continuous light, where the apex of bioactive volatile emission in sender plants was not temporally aligned with the apex of terpene responses in receiver plants (Fig 6). Sender plant emissions even showed slight negative correlations with terpene induction in receiver plants under continuous light. Thus, the emission kinetics of bioactive volatiles from sender plants cannot fully explain the delayed terpene burst.

Green leaf volatile (GLV) and indole emissions in Spodoptera exigua-damaged plants. Light green points represent the mean emission of herbivore-damaged sender plants. Grey points represent mean emissions of undamaged sender plants. Blue rectangles represent the dark phase. Yellow rectangles represent when the lights would typically be turned off. Continuous light-exposed plants were grown under normal light conditions, however lights were left on continuously following the start of the treatment. Error bars = SE. n = 8-10. Compounds were identified based on their molecular weight + 1, as all compounds were protonated. Sesquiterpenes: m/z = 205.20, Monoterpenes: m/z = 137.13, DMNT: m/z = 151.15, TMTT: m/z = 219.21.

Receiver plant terpene emissions are tightly correlated with bioactive sender plant signals under normal light conditions. The left panels depict scatter plot correlation matrices of bioactive volatile emissions from damaged sender plants and terpene emissions from (HIPV)-exposed receiver plants. Only data from the first measurement following the addition of herbivores to sender plants are included. Upper scatter plot (in blue box) shows correlations under normal light conditions and lower scatter plot (yellow box) depicts correlations under continuous light. Continuous light-exposed plants were grown under normal light conditions, however, lights were left on continuously following the start of the treatment. Each black point represents the mean value of all individuals at a given time point after herbivory began. Panels on the right hand side depict heat maps based on the value of Pearson’s correlation coefficient between two given compounds. Numbers in the center of each square are Pearson’s correlation coefficient. Correlation coefficients contained in a pink rectangles indicate a significant correlation (p < 0.05). Hexenyl acetate (HAC), Hexenal (Hexa), Hexen-1-ol (Hexo) and indole (Ind) were from sender plants, and sesquiterpenes (SQT), monoterpenes (MNT), 4,8-dimethylnona-1,3,7-triene (DMNT) and 4,8,12-trimethyltrideca-1,3,7,11-tetraene (TMTT) were from receiver plants.

The combination of direct induction and priming can explain the delayed terpene burst in receiver plants

Based on the above observations, we reasoned that, under natural conditions, the strong delayed volatile burst may be due to the interaction of an initial volatile burst that primes plants for higher volatile release at the onset of the second day, and a secondary trigger in the form of bioactive volatiles released from sender plants at the onset of the second day. To test this hypothesis, we exposed receiver plants to volatiles from control plants or volatiles from herbivory-induced plants for 1.25 hr. We then disconnected the receiver plants and exposed them to clean air for 17 hr, until the beginning of the next day. Half of the plants were then exposed to (Z)-3-hexenyl acetate (HAC) as a secondary trigger. HAC was selected as an inducer as it showed the strongest correlation with all measured terpenes (Fig 6). Following the short exposure to HIPVs, we detected a generally small, but significant, induction of sesquiterpenes, monoterpenes, DMNT and TMTT emissions in receiver plants (Fig 7, Table 2). During the night, no differences between treatments were detected any more. At the beginning of the next day, we found a slight induction in terpene emissions, most apparent for sesquiterpenes and TMTT, in plants that had been exposed to herbivory induced volatiles 18 hr prior (Fig 7, Table 2). (Z)-3-hexenyl acetate exposure at the beginning of the next day also induced terpene emissions. Furthermore, (Z)-3-hexenyl acetate induced a stronger release of all terpenes in plants that had been exposed to HIPVs the previous day. In an independent experiment, we also tested whether HAC sensitivity per se may be higher at the start of the light period. We exposed maize seedlings to HAC dispensers either 30 min or 4 hr 30 min after lights came on. We found induction in terpene emission to be similar regardless of when induction began (Supplemental Fig 2). Thus, exposure to HIPVs prompts maize plants to release more terpenes the next day and also primes maize plants to respond more strongly to a secondary HIPV trigger. Together, these two phenomena result in a pronounced terpene burst.

The combination of volatile priming and sender emission kinetics can explain the delayed terpene burst in receiver plants. Sender plants were connected to receiver plants 30 min prior to herbivore exposure on sender plants and left connected for 1.25 hr following exposure (time between perforated vertical lines). After 1.25 hr, chambers were disconnected and measurements were collected from receiver plant chambers only. The following day, after light was restored, plants were treated with (Z)-3-hexenyl acetate (HAC) dispensers to simulate bioactive signals (indicated by red solid vertical line). * = p < 0.05, ** = p < 0.01, *** = p < 0.001 as determined by aligned rank transformed nonparametric factorial repeated measures ANOVA. Abbreviations: HIPV, herbivore-induced plant volatile; DMNT, 4,8-dimethylnona-1,3,7-triene ; TMTT, 4,8,12-trimethyltrideca-1,3,7,11-tetraene. Colored points represent mean emissions standardized by fresh weight (fw). Error bars = SE. n = 12-16. Compounds were identified based on their molecular weight + 1, as all compounds were protonated. Sesquiterpenes: m/z = 205.20, Monoterpenes: m/z = 137.13, DMNT: m/z = 151.15, TMTT: m/z = 219.21.

Aligned rank transformed nonparametric factorial repeated measures ANOVA results from data presented in Fig 7. Bold values: p < 0.05 and underlined values: p < 0.1. Abbreviations: HAC, (Z)-3-hexenyl acetate; SQT, sesquiterpenes; MNT, monoterpenes; DMNT, 4,8-dimethylnona-1,3,7-triene ; TMTT, 4,8,12-trimethyltrideca-1,3,7,11-tetraene.

Discussion

Herbivore-induced plant volatiles (HIPVs) play an important role in mediating interactions between damaged and undamaged plants [2]. However, the kinetics and temporal dynamics of this information transfer remain poorly understood. We show that, somewhat surprisingly, upon continuous exposure to HIPVs, receiver plants show a pronounced activation of defenses at the onset of the second day. We find that this is the result of volatile-mediated priming on the first day that yields a clocked response to the volatiles that are perceived on the next day. Here, we discuss the mechanisms and biological implications of this phenomenon.

Plants are well known to respond to herbivore-induced volatiles by increasing their own defenses. Defense activation can happen directly, with plants increasing their defense hormone production and volatile release [8,10], and/or via priming effects, with plants increasing their defenses more strongly upon a secondary stimulus [5,7]. Here, we show that these two phenomena operate together to trigger strong defense activation in receiver plants with predictable temporal kinetics. So far, studies in maize and other plants showing priming effects have often exposed receiver plants to HIPVs for several hours, overnight, or even for multiple days, and induced them, e.g. by simulated herbivory, on the following day in the absence of HIPVs [5,8,3638]. Indeed, such a setup would reveal clear priming, with minor or no direct induction by the volatile treatment. Our findings demonstrate that more natural continuous exposure to volatiles results in the same pattern without the need for another stimulus. We thus conclude that HIPVs are sufficient to trigger robust, clocked defense activation in neighboring plants.

What is the mechanism that results in a strong activation of defenses on the onset of the second day of exposure to HIPVs? In maize, green leaf volatiles (GLVs) directly induce and prime jasmonate production and terpene release [7,8,10]. Indole can prime, but not directly induce these two responses [5]. GLVs and indole again can interact to increase defense priming upon a secondary stimulus [7]. Given these considerations, together with the observed volatile release, response kinetics, and our manipulative experiments, we can draw up the following scenario: On day one, maize plants are exposed to GLVs, which trigger a small burst of terpene release in receiver plants. At the same time, receiver plants are primed to respond more strongly on the next day, once light is restored. At the onset of the second day, four things happen simultaneously. First, the emission of GLVs from the neighboring plants increases. Second, the defense priming mechanism responding to these cues kicks into gear. Third, terpene production is activated. Fourth, stomata open and enable volatile diffusion. Together, these elements result in a strong terpene burst. The orchestration of these elements is noteworthy and results in a predictable response pattern under variable HIPV exposure. Further experiments will reveal whether similar patterns are observed in other plants, and how they may operate on a mechanistic level. Additionally, considering we observed an induction in jasmonates coinciding with terpene bursts. In the future, it will be interesting to explore the kinetics of non-volatile jasmonate-dependent defenses, as they are likely to follow similar patterns. Our experiments reveal the importance of tracking plant defense responses in real time as they perceive the dynamic natural volatile blend of herbivore attacked neighbors.

What is the ecological relevance of the clocked temporal kinetics of defense activation in neighboring plants? We propose several hypotheses. First, responding most strongly to volatile cues on the second day may avoid unnecessary energy expenditure [1,13]. GLVs are emitted upon both herbivore attack and mechanical damage [39]. However, when damage is not sustained, GLVs will dissipate rapidly [15]. Thus, responding most strongly to repeated GLV exposure could avoid false negatives and allow plants to respond more robustly to the presence of actual herbivory. Second, sending out volatiles on the second day may maximize indirect defenses. Attracting natural enemies too early to herbivore-free plants could have fitness costs for both the plant and herbivore natural enemies. Responding on the second day could be advantageous, as by this time, the natural enemies would have located and interacted with the herbivores on sender plants, and would thus be ready to move over to the next plant in anticipation of herbivore arrival. As fitness outcomes are likely highly nuanced, especially considering HIPVs also attract herbivores [40], investigation into information transfer on the multi-plant scale, namely in the context of multi-trophic interactions, will be critical. Considerable work will be needed to understand whether the observed defense activation pattern has any adaptive benefit.

Airborne information transfer has the potential to play a major role in plant communities, as it functions as a signaling viaduct between plants not physically connected to one another [4,7,41]. Although it has been established for decades that plants transmit airborne chemical information between individuals, we are only scratching the surface of the dynamic nature of this phenomenon. The kinetics of information transfer at the detailed temporal resolution of this study provide some insights regarding how HIPVs act across time and space. On this basis, a more comprehensive understanding of volatile-mediated information transfer between plants can be built.

Methods

Plant and insect growth

V2-stage Zea mays (maize, B73) was used throughout this study. At this stage, maize has 4 leaves: two are fully developed, one is expanding and one is emerging. Maize plants were grown in commercial potting soil (Selmaterra, BiglerSamen, Switzerland) in 180 ml pots. Plants were grown in a greenhouse supplemented with artificial lights (ca. 300 µmol m-2 s-1). In between light and dark phases, plants were supplemented with lower light (ca. 60 µmol m- 2 s-1) for 15 min to more accurately simulate day : night transitions. The greenhouse was maintained at 22 ± 2 °C, 40–60% relative humidity, with a 14 h : 10 h, light : dark cycle.

Two-week-old maize seedlings were used for all experiments. Spodoptera exigua (Frontier Agricultural Sciences, USA) were reared from eggs on artificial diet [42] and used for experiments when they reached the 4th instar stage.

Herbivore treatment and experimental setup

The herbivore treatment consisted of adding three 4th instar larvae onto plants (sender plants). For controls, sender plants were left undamaged. Two principal experimental setups were used to determine the kinetics of herbivore-induced plant volatile (HIPV) emissions in both sender plants (damaged by herbivores) and receiver plants (exposed to HIPVs from sender plants):

Setup 1

Single plants (senders) in transparent glass chambers (Ø×H 12 × 45 cm) were connected with PTFE tubing to a second transparent glass chamber that was either left empty or filled with a second plant (receiver). Chambers were sealed other than an airflow inlet on the first chamber and an air outlet on the second chamber, allowing HIPVs from sender plants to pass through the second chamber. Cumulative HIPV emissions from senders and receivers were measured as described in the section “Volatile emission sampling”. In order to determine emissions from receiver plants alone, HIPV emissions from sender plants connected to an empty chamber were subtracted from sender plants connected to receiver plants (Supplemental Fig 3).

Setup 2

To isolate emissions solely from receiver plants, chambers containing sender and receiver plants were initially kept separate for basal HIPV emission profiling. Chambers were then connected as in setup 1 for 30 min before herbivore treatment. After herbivores were added to sender plants, chambers remained connected for 1.25 hr, after which time sender and receiver plants were once again disconnected and rearranged so that receiver plants had both clean air flowing through the inlet and an outlet for volatile collection and profiling (Supplemental Fig 3).

Preparation of (Z)-3-hexenyl acetate dispensers

(Z)-3-hexenyl acetate dispensers were prepared as previously described with modifications [7,10]. In brief, 1.5-ml glass vials (Ø×H 11.6 × 32 mm) containing ca. 100 mg glass wool were filled with 200 µl (Z)-3-hexenyl acetate (HAC; >98%; Sigma-Aldrich, Buchs, Switzerland) diluted 50 fold in EtOH and sealed with screw caps containing a rubber septum. The caps were pierced with a 1-µl glass capillary and sealed with PTFE tape. Vials were then covered with aluminum foil and equilibrated for at least 5 days before using. The (Z)-3 isomers of all GLVs measured in this study (hexenyl acetate, hexenal, and hexen-1-ol) are known to be bioactive [8,43] and interconvertible, with hexenal and hexenyl acetate for instance being transformed into hexen-1-ol. As such, we designed the dispensers to emit a molar concentration of HAC comparable to the molar concentration of all (Z)-3 GLVs emitted by herbivore-infested plants. HAC dispensers emitted HAC at a rate of 4.68 (±0.40 SD) nmol/hr and maize seedlings infested with 3 4th-instar S. exigua larvae emitted all GLVs at a rate of 1.94 nmol/hr (±0.25 SD). Of note, GLV emissions induced by caterpillars vary over time, and can be more than 2-fold higher than the average during times of strong active feeding (Supplemental Fig 4). Thus, the release rate of the dispensers is within the plant’s physiological range.

Volatile sampling

Volatile emissions were sampled using setup 1 and 2. Volatile profiling was measured with a high-throughput platform comprised of a proton transfer reaction time-of-flight mass spectrometer (PTR-ToF-MS; Tofwerk, Switzerland) and an automated headspace sampling system (Abon Life Sciences, Switzerland) supplied with clean airflow (0.8 L min-1). An outlet on receiver plants was accessible to the autosampler/PTR-ToF-MS system. The PTR- ToF-MS system drew air at 0.1 L min-1. Between samples, a zero gas measurement was performed for 3 s to avoid contamination. At each time point volatiles were continuously measured for 15-25 s and averaged to a single mean per sample. Complete mass spectra (0- 500 m/z) were recorded in positive mode at ca. 10 Hz. The PTR was operated at 100 °C and an E/N of approximately 120 Td. The volatile data extraction and processing were conducted using Tofware software package v3.2.2 (Tofwerk, Switzerland). Protonated compounds were identified based on their molecular weight + 1. During volatile collection LED lights (DYNA, heliospectra) were placed ca. 80 cm above the glass cylinders and provided light at ca. 300 μmol m-2 s -1. Identical light : dark cycle timing as in the greenhouse for plant growth was used. In order to quantify GLVs emitted from herbivore-infested plants and HAC dispensers, we collected total volatile emissions for 1 hr using an identical volatile collection and analysis method described in Erb et al., 2015. In brief, volatiles were collected in Super-Q adsorbent traps, which were then eluted in dichloromethane and injected into a GC-MS (Agilent, USA). Compound identities were confirmed with mass spectroscopy analysis and similarity to library matches (NIST search 2.2 Mass Spectral Library, USA). Emission rates of the Z-3 isomers of GLVS were quantified using a standard curve with synthetic compounds.

Foliar terpene pools and gene expression

Setup 1 was used to determine terpene pools and gene expression. After 3, 8, 16.75 and 22 hr of exposure to HIPVs, the oldest developing leaf of receiver plants was harvested and flash frozen on liquid nitrogen. The oldest developing leaf was chosen as it is the largest and highly responsive to bioactive HIPVs [10]. A new set of plants was used for each time point. Analysis of foliar terpene pools was conducted using slightly modified, previously described, methods [44]. In brief, ca. 15 mg of ground fresh frozen leaf tissue was added to a 20-ml precision thread headspace glass vial sealed with a magnetic screw cap fitted with a silicone/PTFE septum (Gertel GmbH & Co. KG, Germany). Immediately after a vial was prepared, an SPME fiber (100 μm polydimethylsiloxane coating; Supelco, USA) was inserted into the vial and volatiles were collected for 40 min at 50 °C. After collection, volatiles were thermally desorbed for 3 min at 220 °C and analyzed using GC-MS (Agilent, USA). Helium was used as the carrier gas at a flow-rate of 1 ml min-1 with a temperature gradient of 5 °C/min from 60 °C (1 min hold) to 250 °C. Compound identification was based on similarity to library matches (NIST search 2.2 Mass Spectral Library, USA). For quantification of gene expression, total RNA was extracted and purified from ca. 80 mg ground fresh frozen tissue using the GeneJET plant RNA extraction kit following the manufacturer’s instructions.

Genomic DNA was removed from 1 μg purified RNA using gDNA Eraser (PrimeScript RT Reagent Kit, Perfect Real Time) following manufacturer’s instructions (Takara Bio Inc., Kusatsu, Japan). Reverse transcription and cDNA was synthesized using PrimeScript Reverse Transcriptase (TaKaRa Bio). Gene expression was determined with quantitative reverse transcription polymerase chain reaction (qRT-PCR) using ORA SEE qPCR Mix (highQu GmbH, Germany) on an Applied Biosystems QuantStudio 5 Real-Time PCR system. The normalized expression (NE) values were calculated as: NE=(1/(PE Cttarget))/(1/(PE Ctreference)) where PE = primer efficiency and Ct = cycle threshold [45]. Ubiquitin (UBI1) was used as the reference gene. Gene identifiers and primer sequences are listed in Supplemental Table 1.

Phytohormone quantification

Plant treatments and tissue collection for phytohormone were identical to foliar terpene, and gene expression analysis. The phytohormones jasmonic acid (JA), JA-isoleucine (JA-Ile) and 12-oxo-phytodienoic acid (OPDA) were extracted, analzyed and quantified using a slightly modified version of the method detailed in [46]. In brief, ca. 80 mg of finely ground fresh frozen leaf tissue was extracted with 1 ml ethylacetate:formic acid (99.5:0.5, v/v) spiked with isotopically labelled forms of the abovementioned phytohormones. d5-JA was acquired from CDN Isotopes (Canada), d5-OPDA was acquired from OLChemIm (Czechia), and 13C6-JA-Ile was synthesized at the University of Neuchâtel according to a previously described method [47]. Extracts were dried and re-suspended in 200 μl 50% MeOH using a sonicating bath for 30 min. Phytohormones were analyzed using UPLC-MS-MS fitted with an Acuity BEH C18 column (Waters, USA) using a flow-rate of 400 μl/min and an injection volume of 2 μl. The two mobile phases used were 0.05% formic acid in water (A) and 0.05% formic acid in acetonitrile (B) with the following gradient conditions: 5-50% B over 5 min, 60-100% B over 3 min, 100% B for 4 min and a final re-equilibration at 5% B for 4 min.

Statistical analyses

All statistical analyses were performed in R version 4.2.2 [48]. Volatile emissions were analyzed by aligned rank transformed nonparametric factorial repeated measures ANOVA using the package ARTool, as individual plant emission kinetics were measured repeatedly over time [49]. Internal terpene pools, gene expression and phytohormone levels were analyzed using Welch’s t-tests between control and HIPV-exposed receiver plants within each time point.

Data availability statement

All data used to generate figures and tables presented in this study are available as supplemental material (Appendix A).

Acknowledgements

We would like to thank Dr. Bernardus Schimmel for assistance with gene expression quantification. Additionally, we thank all members of the Biotic Interactions and Chemical Ecology groups at the University of Bern for helpful discussions. This work was supported by the Swiss National Science Foundation (Grants Nr. 210651 and 200355), the State Secretariat for Education, Research, and Innovation SERI (Project CANWAS), the Horizon 2020 Marie Skłodowska-Curie Actions (Grant Nr. 886651) and the University of Bern.

Author Contributions

JMW and ME conceived the project. JMW, TC, LW and ME designed experiments. JMW and TMC conducted experiments. JMW, GG and ME analyzed and interpreted data. JMW wrote the initial draft of the manuscript with critical revisions from ME and input from all authors.

Supplemental Tables

Gene identifiers and qRT-PCR primer sequences used in this study.

Supplemental Figures

The delayed volatile burst is conserved under continuous light, regardless of the time of day. Emission kinetics of herbivore-induced plant volatile (HIPV)-induced terpenes in undamaged receiver plants under continuous light are shown. Herbivory on sender plants ca. 1-hr before lights would typically be turned off (around 20:00 hr). Plants were grown under normal light conditions. Lights were left on continuously following the start of the treatment. Dark green points represent the mean emission of herbivore damaged sender plants connected to undamaged receiver plants with the emissions from damaged sender plants only subtracted. Black points represent undamaged sender plants connected to undamaged receiver plants, with the emissions from undamaged sender plants only subtracted. Yellow rectangles represent when the lights would typically be turned off. Abbreviations: DMNT, 4,8- dimethylnona-1,3,7-triene; TMTT, 4,8,12-trimethyltrideca-1,3,7,11- tetraene. Error bars = SE. n = 6. Compounds were identified based on their molecular weight + 1, as all compounds were protonated. Sesquiterpenes: m/z = 205.20, Monoterpenes: m/z = 137.13, DMNT: m/z = 151.15, TMTT: m/z = 219.21.

Maize seedlings respond consistently to green leaf volatiles (GLVs) throughout the day. Seedlings were exposed to (Z)-3-hexenyl acetate (HAC) either at 07:40 hr, ca. 30 min after lights came on (Immediate, yellow points) or at 11:40 hr, ca. 4 hr and 30 min after (4-hr delay, blue points). Control (untreated) emission levels were subtracted from both treatments. Yellow and blue perforated vertical lines represent the moment that HAC dispensers were added to the immediate and 4-hr delay plants, respectively. Abbreviations: DMNT, 4,8- dimethylnona-1,3,7-triene; TMTT, 4,8,12-trimethyltrideca-1,3,7,11- tetraene. Error bars = SE. n = 4. Compounds were identified based on their molecular weight + 1, as all compounds were protonated. Sesquiterpenes: m/z = 205.20, Monoterpenes: m/z = 137.13, DMNT: m/z = 151.15, TMTT: m/z = 219.21.

Experimental setup schemes for volatile profiling

Green leaf volatile (GLV) emission from herbivore-damaged maize seedlings is highly variable over time. Representative GLV emissions from a single attacked plant on the second day of herbivory (18-24 hr of exposure to HIPVs) are shown in black. The perforated red horizontal line represents the average across all time points presented. Compounds were identified based on their molecular weight + 1, as all compounds were protonated. Sesquiterpenes: m/z = 205.20, Monoterpenes: m/z = 137.13, DMNT: m/z = 151.15, TMTT: m/z = 219.21.