Hyperspectral imaging to characterize plant–plant communication in response to insect herbivory
Plant Methodsvolume 14, Article number: 54 (2018)
In studies of plant stress signaling, a major challenge is the lack of non-invasive methods to detect physiological plant responses and to characterize plant–plant communication over time and space.
We acquired time series of phytocompound and hyperspectral imaging data from maize plants from the following treatments: (1) individual non-infested plants, (2) individual plants experimentally subjected to herbivory by green belly stink bug (no visible symptoms of insect herbivory), (3) one plant subjected to insect herbivory and one control plant in a separate pot but inside the same cage, and (4) one plant subjected to insect herbivory and one control plant together in the same pot. Individual phytocompounds (except indole-3acetic acid) or spectral bands were not reliable indicators of neither insect herbivory nor plant–plant communication. However, using a linear discrimination classification method based on combinations of either phytocompounds or spectral bands, we found clear evidence of maize plant responses.
We have provided initial evidence of how hyperspectral imaging may be considered a powerful non-invasive method to increase our current understanding of both direct plant responses to biotic stressors but also to the multiple ways plant communities are able to communicate. We are unaware of any published studies, in which comprehensive phytocompound data have been shown to correlate with leaf reflectance. In addition, we are unaware of published studies, in which plant–plant communication was studied based on leaf reflectance.
Plant defenses to insect herbivory include a wide range of co-evolutionary adaptations, and they are intensively studied from evolutionary, ecological and crop management perspectives. The leaf tissue damage by herbivores elicit complex cascades of physiological responses in surrounding cells (locally responses) and expression of hundreds of genes, and they induce whole-plant hormonal responses [1,2,3]. Plant defenses to insect herbivory involves triggering of so-called “herbivore-associated molecular patterns” (HAMPs) [3, 4], which lead to synthesis and emission of volatile organic compounds (VOCs) [1,2,3]. The synthesis and emission of VOCs in response to herbivory may enable faster systemic plant defense responses than hormonal signaling via the vascular tissue . Many VOCs are known to be involved in the regulation of plant responses to herbivory, including jasmonic acid (JA) [3, 5,6,7], salicylic acid (SA) [3, 8,9,10], and reactive oxygen species [11,12,13]. For instance, it was demonstrated experimentally that mechanically damaged sagebrush plants released a pulse of methyl jasmonate that induced direct resistance in wild tobacco . Furthermore, wild tobacco plants with clipped sagebrush neighbors had increased levels of putative defensive oxidative enzymes (polyphenol oxidase) and showed reduced levels of leaf damage by grasshoppers and cutworms .
According to Maffei et al. , reactive oxygen species and intracellular calcium signatures belong to early events of plant defenses in response to biotic stressors, and they are responsible for most of the cascading biochemical reactions. The oxidative metabolism, particularly hydrogen peroxide synthesis, is related to a wide variety of reactions and signaling cascades, which are necessary for all aspects of plant growth and defense against biotic and abiotic stresses . In addition, hydrogen peroxide is involved in: development of individual root hairs xylem differentiation and lignification, wall loosening and cross-linking, root/shoot coordination and stomatal control, and hypersensitivity reactions [13, 17]. In addition to aerial plant–plant communication, there is a growing body of research describing below-ground plant–plant communication through the roots [18,19,20,21,22].
To date, most studies of plant defenses to insect herbivory and plant–plant communication have been characterized and quantified by means of either headspace collection of volatiles in combination with gas chromatography/mass spectrometry analysis [23,24,25] or through invasive sampling of leaf tissue to quantify up/down regulations of phytocompounds [26,27,28]. In addition, experimental infestation of plants with arthropods is commonly used to confirm induction of plant defenses by a wide range of stressors and to confirm plant–plant communication [9, 29,30,31,32,33,34,35]. However, such arthropod bioassays have the potential to further alter plant traits and do not enabled non-invasive monitoring of plant response over time. Availability of a non-invasive technique to quantify plant responses to biotic stressors, such as, a herbivore, would greatly improve the ability to characterize plant responses to biotic stressors over time. Moreover, such techniques would enable more precise detection of when (based on time and intensity of stressor) significant physiological changes occur and therefore when invasive tissue sampling should be performed in order to optimize the likelihood of detecting physiological and molecular responses. A non-invasive method to detect and quantify plant–plant communication would be of considerable importance in ecological studies of insect herbivory and adaptive/evolutionary plant responses. Finally, reflectance based detection of spectral bands responding to plant–plant communication may be of considerable relevance to the application of remote sensing technologies in agricultural systems, as it may greatly enhance the ability to detect emerging insect infestation hotspots.
The research field of plant phenomics is largely driven by deployment of a range of molecular and imaging technologies [36, 37]. The interest in plant imaging as part of high throughput phenotyping and modern crop breeding, reflects the growing understanding of and appreciation for the many important and complex ways that growing plants respond to their surrounding environment and stressors, including physiological responses to herbivory by insects. For almost three decades, it is known that general leaf compounds, such as lipids, oils, protein, nitrogen, lignin, starch, cellulose, sugars, and chlorophyll can be quantified based on spectral information in specific portions of the solar incident spectrum [38, 39]. As part of applications of imaging technologies to studies of plant responses to stressors , a considerable body of research describes detection and characterization of chlorophyll content in growing plants [40,41,42] and more broadly detection of leaf reflectance responses to biotic plant stressors [43, 44]. Despite a large body of research into applications of imaging technologies to detect and characterize plant stress signals, there are only a limited number of studies, in which changes in reflectance features have been associated with phytocompounds other than leaf pigments, such as, chlorophyll [42, 45, 46]. An exception is a number of studies, in which the authors demonstrated association between specific reflectance features and leaf potassium content [47,48,49]. Use of leaf reflectance data to detect and characterize plant responses to abiotic and biotic stressors is based on the assumption that induced stress interferes with photosynthesis, chemical composition, and physical structure of the plant and affects the absorption of light energy and thus alters the reflectance spectrum of the plants .
In this study, the main hypothesis was that insect herbivory causes changes in leaf phytocompound levels, and these physiological defense responses are associated with detectable changes in phytocompound levels and in certain spectral bands of leaf reflectance profiles. As a secondary hypothesis, we predict that plant–plant communication (from plant with herbivory to an adjacent control plant without herbivory) will elicit both a change in phytocompound composition of leaves and also cause a corresponding change in leaf reflectance. To address these hypotheses, we acquired time series hyperspectral imaging data [before (baseline) and after 6, 12, 24, and 48 h of herbivory] from maize plants (Zea mays L.) inside cages subjected to the following four treatments (Fig. 1): (T1) individual non-infested plants (control). (T2) individual plants experimentally subjected to herbivory by green belly stink bug Dichelops melacanthus (Dallas) (Hemiptera: Pentatomidae). (T3) one plant subjected to insect herbivory (T3A) and one control plant (T3B) in separate pots but inside the same cage, and therefore the possibility of plant–plant communication via air space. (T4) one plant subjected to insect herbivory (T4A) and one control plant (T4B) together in the same pot, and therefore possibility of plant–plant communication via roots and air space. It is important to emphasize that this study was based on maize plants only being exposed to a low level of induced herbivory with no visible symptoms of insect herbivory. Using a separate set of maize plants under the same treatments, we obtained phytocompound data [jasmonic acid (JA), 12-oxo-phytodienoic acid (OPDA), auxins (in particular indole-3acetic acid (IAA)), salicylic acid (SA), hydrogen peroxide, carotenoids, and chlorophyll a and b] from all treatments and most time points. We are unaware of any published studies, in which comprehensive phytocompound data have been shown to correlate with leaf reflectance. In addition, we are unaware of published studies in which plant–plant communication was studied based on leaf reflectance. Thus, this study provides considerable novelty to the fields of insect–plant interactions and plant phenomics.
Test plants and insects
Potted (1 L) maize plants [Zea mays L. (Poaceae)] of hybrid Pioneer P2530 (non Bacillus thuringiensis) with four leaves [stage V3 ] were used in this study. Green belly stink bugs [Dichelops melacanthus (Dallas)] were obtained from a colony maintained in a climate-controlled room [temperature 25 ± 2 °C, relative humidity (RH) 60 ± 10% and a photoperiod of 14:10-h L:D] for at least 10 generations. Soybean [Glycine max (L.) Merrill (Fabaceae)] grains and plants were used as feeding and oviposition substrates, respectively. In Brazil, the green belly stink bug is a serious pest when feeding on young stems of maize plants . Infestation of maize plants consisted of transferring two adult couples (four individuals) to a tulle fabric meshed bag (15 cm × 10 cm) , which was secured around the stem and lowest leaf of each maize plant (Fig. 1). Each experimental unit (either a single, treatments 1 or 2) or two, treatments 3 or 4 maize plant) were placed inside a wooden cage (30 × 25 × 25 cm), in which all six sides were covered with 100 µm plastic mesh. This choice of cages was to minimize the volatile exchange among experimental units, and at the same time enable some air to circulate to minimize humidity build-up and to ensure similar temperature conditions inside experimental units. All experimental units were kept in the same climate-controlled room (temperature 25 ± 2 °C, 60 ± 10% RH and 14 L: 10 D h photoperiod), with all experimental units exposed to the same lighting regime (illuminance not measured). All treatments were replicated five times for hyperspectral imaging and three times for phytocompound data collections. As cages remained closed at all times when inside the climate-controlled room and only opened for a few minutes to perform hyperspectral imaging, the risk of plant–plant communication was kept as low as possible.
Hyperspectral imaging of maize leaves
We acquired time series hyperspectral imaging data: before (baseline) and after 6, 12, 24, and 48 h of herbivory. After baseline hyperspectral imaging, individual maize plants were randomly assigned to treatments. Hyperspectral imaging was performed in a room immediately adjacent to the climate-controlled room, where the experimental units were maintained. At each hyperspectral imaging event, a single cage at a time was transferred from the climate-controlled room to the hyperspectral imaging room, and hyperspectral imaging was completed within 1–2 min per plant. The order of maize plants and experimental units being imaged was randomized each time.
Hyperspectral imaging conditions and settings are similar to those described in a range of other hyperspectral imaging studies [49, 53,54,55,56,57,58]. All hyperspectral images were collected with artificial lighting from 15 W, 12 V light bulbs mounted in 2 angled rows, one on either side of the lens, with 3 bulbs in each row. A voltage stabilizer (Tripp-Lite, PR-7b, www.radioreference.com) powered the lighting. Ambient climate conditions were between 24–26 °C and 50–60% relative humidity. A piece of white Teflon was used for white calibration, and “relative reflectance” refers to proportional reflectance compared to that obtained from Teflon. Consequently, relative reflectance values ranged from 0 to 1. For each combination of maize plant and time point, we imaged an area corresponding to about 2 cm2 (equivalent to 64,000 pixels) on the second leaf after it had been placed in a horizontal position. Each hyperspectral image from a given combination of plant and time point was divided into two halves. A hyperspectral push broom camera (PIKA II, Resonon Inc., Bozeman, MT) was used. The objective lens had a 35 mm focal length (maximum aperture of F1.4) and was optimized for the visible and NIR spectra. The main specifications of the spectral camera are as follows: interface, Firewire (IEEE 1394b), output, digital (12 bit), 160 bands (spectral) by 640 pixels (spatial), angular field of view, 7 degrees, and spectral resolution of < 3 nm. We acquired reflectance in 240 spectral bands in the range from 380 to 1031 nm. However, to avoid spectral bands with low signal/noise ratio, we only included 230 spectral bands from 405 to 1012 nm, and they were spectrally 3X-binned (averaged) into 75 spectral bands with a spectral resolution of about 9 nm. We calculated average leaf reflectance in each of the spectral bands for all combinations of maize plant and time point and used these data for statistical analyses. With five time points (0, 6, 12, 24, and 48 h), six plant treatments (Fig. 1), five replications, and two hyperspectral images, a total of 300 average leaf reflectance profiles (observations) were included in the statistical analyses. In other words, the number of observations was 4 times the number of spectral bands (N = 75). This ratio between observations and spectral bands is very important for the robustness of the data classification based on average leaf reflectance. Moreover, model over-fitting due to the Hughes phenomenon or violation of the principle of parsimony  is a major concern when the number of explanatory variables is similar or exceeds the number of observations [44, 58, 60, 61].
Phytocompound analyses were carried out using the second leaf of maize plants from all treatment groups (Fig. 1) before infestation (baseline) and after 12, 24, and 48 h. Three replicates were used for each combination of treatment and exposure time.
A relative quantification of four phyto-signaling compounds (IAA, JA, SA, and OPDA) was performed using a method adapted from Trapp et al. . After collection leaves were immediately frozen in liquid nitrogen and storage at − 80 °C until the extraction.
The leaves were ground in a crucible, and 100 mg of ground tissue was transferred to 1.5 mL tubes and extracted twice with 1 mL of methanol (70%, v v−1). After the extraction, the solutions were dried in a SpeedVac and re-suspended in 100 µL of methanol containing 20 ng mL−1 of each internal standard (d5-IAA, d5-JA, d4-SA). These samples were analyzed by HPLC–MS/MS on an Agilent 1100 HPLC system (Agilent Technologies, Böblingen, Germany) connected to a LTQ Iontrap mass spectrometer (Thermo Scientific, Bremen, Germany) equipped with an Electrospray ionization source.
The chromatographic separation was performed in a Luna Phenyl-Hexyl column (150 × 4.6 mm, 5 μm; Phenomenex, Aschaffenburg, Germany), using gradient elution with formic acid (0.05%, v v−1) and methanol (with 0.05% of formic acid, v v−1) as mobile phases A and B, respectively. The elution profile was: 0–10 min, 42–55% B in A; 10–13 min, 55–100% B; 13–15 min 100% B; 15–15.1 min 100–42% B in A; and 15.1–20 min 42% B in A. The mobile phase flow rate was 1.1 mL min−1.
The ionization parameters were optimized for each phyto-signaling compound separately. Except IAA and d5-IAA (which were analyzed in positive ionization mode), all the others phyto-signaling compounds were analyzed in negative mode and monitored by the following SRM channels: IAA (176 > 130); d5-IAA (181 > 135); JA (209 > 59); d5-JA (214 > 62); OPDA (291 > 165); SA (137 > 93); d4-SA (141 > 97). Precursors and fragment ions were selected with an isolation window of 2 days.
Photosynthetic pigments (chlorophyll a and b, carotenoid) were measured based on fresh material from second leaves of maize plants (same as area being imaged with hyperspectral camera). Leaves of maize (fresh weight, 0.1 g) were frozen in liquid nitrogen and ground. Subsequently, ground leaf samples were transferred to a glass vial containing 5 mL of aqueous acetone solution (80%, v v−1) for 24 h, filtered and the absorbance was measured in an Analytik Jena spectrophotometer (Specord 200 Plus, Analytic Jena AG, Jena, Germany) at 663.6, 646.8 and 470 nm. Carotenoid and chlorophyll concentrations were calculated from the absorbance of extract at using the adapted formula by Lichtenthaler HK  as below:
To quantify hydrogen peroxide content, fresh tissue was collected from second leaf on maize plants (about 50 mg), and each sample was ground to a fine powder in liquid nitrogen and were homogenized in ice bath with 5 mL of trichloroacetic acid (0.1%, w v−1). The homogenate was centrifuged at 20,000 g for 30 min at 4 °C. The supernatant (0.5 mL) was added to 0.5 mL 10 mM potassium phosphate buffer (pH 7.0) and 1 mL of a 1 M L−1 potassium iodide solution. The absorbance of supernatant was checked at 390 nm . The hydrogen peroxide content was calculated based on a standard curve prepared with hydrogen peroxide solutions (0–500 µM L−1) obtained by dilution of hydrogen peroxide (30%, w w−1).
Experimental design and data analysis
Data processing and analyses were conducted in PC-SAS 9.4 (SAS Institute, NC). Initially, we conducted an analyses of variance (proc anova option = repeated and option = tukey) for repeated measures. An initial analyses of variance for repeated measures showed a highly significant effect of hours of insect herbivory (P < 0.001), and also a highly significant interaction between hours of insect herbivory and plant treatment (P = 0.001). The key focus of this study was to examine possible differences in maize plant responses to treatments, so effect of hours of insect herbivory was eliminated by rank-transformation of phytocompound data [ranked within each time point (0, 12, 24, and 48 h)] before further analysis. Using baseline data only (prior to insect herbivory), we conducted analysis of variance (proc anova) to examine the effect of the experimental design by comparing phytocompound levels and leaf reflectance from maize plants: (1) alone, (2) two plants in separate pots but in the same cage (T3A and T3B), and (3) two plants in the same pot (T4A and T4B). To analyze responses to insect herbivory and plant–plant communication, the maize plants were assigned to one of four groups:
Control: baseline data from T1 and T2 plants, and all additional data acquired from T2 plants,
Herbivory: all data acquired from T1, T3A, and T4A plants after onset of herbivory.
Air only: all data acquired from T3B plants after onset of herbivory.
Air and root: all data acquired from T4B plants after onset of herbivory.
We conducted analyses of variance (proc anova option = repeated and option = tukey) for repeated measures (with either rank-transformed phytocompound data or selected spectral bands) to compare average maize plant responses by these four groups.
The main data analysis was conducted in two separate parts, and, in separate analyses, we used the same experimental design in studies of phytocompound data and leaf reflectance data. Firstly, the purpose was to characterize maize plant responses to insect herbivory; secondly, the purpose was to characterize plant–plant communication. In order to characterize maize plant responses to insect herbivory, we examined data from the Control and Herbivory groups described above. In other words, maize plants that were either alone and not infested or plants that were directly infested, and a discriminant classification model (proc discrim)  was developed. Initially, stepwise linear discriminant analysis (proc stepwise) was used to only select the explanatory variables (phytocompound data or spectral bands) with significant contribution to each of the discriminant classification models. The selected subset of explanatory variables was used to generate the discriminant classification models, and their accuracies were quantified on the basis of 80% of the data being randomly selected as training data set and the remaining 20% of the data used for independent validation. This validation procedure was repeated 10 times to calculate average classification accuracies.
The second part of the main data analysis consisted of characterizing plant–plant communication, and this was conducted based on plants assigned to the “Air only” and “Air and root” groups described above. Moreover, the discriminant classification models (based on phytocompound data or leaf reflectance data from Control and Herbivory groups of plants) were used to interpret data from T3B and T4B plants. Under the assumption of plant–plant communication, T3B and T4B plants were predicted to be classified as subjected to herbivory (and not control).
Phytocompound responses to insect herbivory
At baseline, none of the eight phytocompounds varied significantly among the three groups of plants (P > 0.05). Thus, the experimental design (presence/absence of a second maize plant in the same cage) did not appear to significantly affect the phytocompound levels. Figure 2 shows average phytocompound levels after 12–48 h, after maize plants had been divided into four groups. Based on tukey comparisons, three phytocompounds (IAA, JA, and hydrogen peroxide) varied significantly among the four groups (P < 0.05), but only IAA showed a significant difference (up-regulation) in response to insect herbivory (Control vs. Herbivory).
Based on stepwise forward selection, we determined that IAA, hydrogen peroxide and chlorophyll a and b contributed significantly to the discriminant classification of Control and Herbivory plants. The validation of the discriminant classification model showed that the four phytocompounds could be used to classify Control and Herbivory plants with 80.7% accuracy (Control = 88.4% and Herbivory = 72.9%). Thus, although only IAA showed a significant up-regulation in response to insect herbivory (Fig. 2), a combination of phytocompounds provided strong indication of a maize plant response to insect herbivory. In following, the discriminant classification model was used to interpret phytocompound from T3B and T4B plants. We found that 67% of T3B plants were classified as Herbivory plants, which suggested that their phytocompound composition was more similar to that of plants subjected to insect herbivory than to control plants. We also found that 89% of T4B plants were classified as Herbivory plants. Thus, we demonstrated plant–plant communication, and this was particularly evident when root communication was possible.
Reflectance detection of insect herbivory and plant–plant communication
Figure 3a shows average baseline reflectance profiles of maize plants divided into three experimental groups. Average reflectance curves were also presented as their difference from that of Control plants (Fig. 3b). As the curves were generally < 1, it is seen that presence of a second plant inside cages generally caused a decrease in average leaf reflectance, especially in T4B plants in the spectral range from 405 to 700 nm. As indicated by black dots in Fig. 3b, eight spectral bands (near 500, 630, and 700 nm) showed significant difference in response to the experimental design.
Figure 4 shows the average leaf reflectance profiles of maize plants, when divided into four treatment groups. It is important to highlight that hyperspectral imaging data were acquired under controlled experimental conditions and that only a low and brief level of herbivory-induced stress level was imposed. Using stepwise forward selection, we found that 10 spectral bands (black dots in Fig. 4) contributed significantly to the discriminant classification model of Control and Herbivory plants. Based on analysis of variance for repeated measures, none of the 10 selected spectral bands showed significant difference in response to plant treatments (not shown). In other words, none of the spectral bands could be used individually as indicators of insect herbivory. However, we used the 10 spectral bands as explanatory variables in a discriminant classification model to classify Control and Herbivory plants, and this classification was associated with an overall accuracy of 79.3% (Control = 79.4% and Herbivory = 79.2%). Thus, despite only subtle between-class variation in average leaf reflectance data, it was possible to identify a reliable leaf reflectance response to low level of insect herbivory. In the second part of the hyperspectral imaging data analysis, the discriminant classification model (derived from plants assigned to Control and Herbivory) was used to interpret leaf reflectance data acquired from T3B and T4B plants. Interestingly, 76% (T3B plants) and 77% (T4B plants) were classified as subjected to insect herbivory. This result strongly suggested that plant–plant communication from infested plants (T3A and T4B) to T3B and T4B plants triggered a detectable leaf reflectance response, so that the average reflectance profiles of T3B and T4B plants were more similar to those of Herbivory plants.
This study was conducted to determine whether phytocompound data and hyperspectral imaging data could be used to detect maize plant response to both direct insect herbivory and also to plant–plant communication from an infested to a non-infested plant. We demonstrated that individual phytocompounds (except IAA) as well as reflectance in individual spectral bands were not reliable indicators of neither insect herbivory nor–plant communication. However, using a linear disciminination classification method based on combinations of either phytocompounds or spectral bands, we found clear evidence of maize plant responses. The potential of using changes in average reflectance data as indicator biotic stress is of particular interest to a broader community of researchers studying plant–insect interactions.
In a large body of literature, arthropods have been used as indicators of plant responses to a range of stressors and to plant–plant communication [9, 29,30,31,32,33,34,35]. As part of elucidating the mechanisms responsible for such bio-responses, it is well-established that herbivore-induced stress elicits plant physiological responses, which lead to changes in phytocompound levels (including OPDA, JA, IAA, and hydrogen peroxide) [1, 2, 65,66,67], and this was partially confirmed in this study. It is also well-described (and corroborated in this study) that herbivore-induced stress in one plant can elicit a phytocompound response in adjacent control plants, either induced by VOCs [1, 2, 14, 15, 68,69,70] and/or through root-based signaling [18,19,20,21,22]. Finally, there is ample evidence of how biotic stress by herbivores induces detectable changes in leaf reflectance [43, 44]. However, this is the first report describing herbivore-induced stress in plants subjected to insect herbivory causing a detectable and statistically significant change in leaf reflectance in an adjacent control plant. Thus, we are providing the first report of plant–plant communication of biotic stress being detected and quantified through the means of non-invasive method, namely analyses of leaf reflectance.
Plants have developed various strategies to defend themselves against herbivores and pathogens either alone or in “community” tactics . Although some of these strategies are constitutive, most are inducible in response to stressors, so that plants enter a “primed state” [71, 72]. In general, induced stress interferes with photosynthesis, chemical composition, and physical structure of the plant and affects the absorption of light energy and thus alters the reflectance spectrum of the plants under stress . In plants, immediate responses to herbivory are typically followed by a rapid and transient generation of reactive oxygen species in response to damage in the plant tissue [11,12,13, 73]. In a concentration-dependent manner, the hydrogen peroxide is a reactive oxygen species responsible by important functions in plant development and metabolism such as homeostasis, acclimatization, and plant defense processes . In addition, hydrogen peroxide acts as a signaling molecule in plants to a variety of others signaling molecules and plant hormones [14, 75, 76], since it can induce the synthesis or activation of transcription factors that are associated with the induction of several enzymes of the antioxidant system [77, 78].
Literature reviews of hyperspectral vegetation indices have identified spectral bands in the 900–940 nm region to correlate well with wet biomass and leaf area index across a series of crops [79, 80]. It is also been determined that spectral bands near 948, 975 and 1004 nm are meaningful when developing spectral band indices to detect iron deficiency in crops . In addition, the relationship reflectance at 970 nm/900 nm was proposed as a meaningful water index of crops . However, the vast majority of hyperspectral indices are based on spectral bands from 450 to 800 nm [79, 81,82,83,84,85]. In the current study, eight spectral bands (near 500, 630, and 700 nm) showed significant difference in response to the experimental design (Fig. 3b). This result implies that some level of plant–plant communication was already occurring prior to the onset of insect herbivory. We also demonstrated that a combination of 10 spectral bands (see black dots in Fig. 4), especially between 700 and 750 nm, could be used to classify maize plants with/without insect herbivory. Most studies of reflectance responses to insect-induced stress have shown that herbivory causes an increase in reflectance . Sensitivity analyses across wide spectral ranges (between 400 and 2500 nm) have shown that increased reflectance near 700 nm wavelength is a consistent indication of stress [86, 87]. Moreover, an increase in reflectance near 700 nm appears to be linked to reduced absorption due to lower chlorophyll concentration [42, 46, 88]. However especially regarding piercing-sucking herbivores, there are important exceptions suggesting that feeding behavior may be important. As an example, greenbug [Schizaphis graminum (Rondani) (Hemiptera: Aphididae)] infestations in wheat [Triticum aestivum L. (Poaceae)] caused a decrease in reflectance in the UV-light portion of the spectrum . In addition, herbivory by bird cherry-oat aphids [Rhopalosiphum padi (L.) (Hemiptera: Aphididae)] also caused a decrease in leaf reflectance . In a study of Russian wheat aphids [Diuraphis noxia (Mordvilko) (Hemiptera: Aphididae)] and greenbugs (S. graminum), both sucking herbivorous insects, it was determined that the strongest reflectance responses were detected in spectral bands between 625–635 and 680–695 nm . The results presented in this study underscore the importance of disentangling reflectance responses associated with inherent and more general plant–plant communication from those being potentially providing specific information about signaling of biotic stress imposed by herbivorous insects. As this is the first study of plant–plant communication and leaf reflectance data, we do not have comparative data.
An important detail in this study was the structure and handling of cages used to hold the maize plants. That is, we decided to keep all cages in the same room, as it was considered critical that they were all under the same abiotic conditions (i.e. light, temperature and relative humidity). In addition due to concerns about mold and build-up of humidity (especially when two plants were in a cage compared to cages with one plant), it was considered necessary to allow some air movement in and out of cages. Also, use of hermetically closed cages could also affect the feeding behavior and survival of the stink bugs. However at the same time, we obviously had concerns about possible volatile-based communication among plants in different cages. It is possible that more significant plant responses would have been detected, if we had been able to control and maintain air flow inside each of the cages, but that was not possible due to logistical limitations.
In conclusion, we demonstrated plant–plant communication from an infested to a non-infested plant when based on phytocompound data as well as on average leaf reflectance data from a selected combination of spectral bands. Thus, we have provided initial evidence of how remote sensing may be considered a powerful non-invasive method to increase our current understanding of both direct plant responses to biotic stressors but also to the multiple ways plant communities are able to communicate. Due to plant–plant communication, our results suggest that in laboratory studies thorough caging and/or separation of replicated plants, and of plants subjected to different treatments, is of major importance when conducting experimental laboratory studies. We are aware of how growth chamber space at universities is often shared by students and faculty from multiple groups. In addition, multiple plants may be grown in flats or trays to save space. The results presented here suggest that interference via plant–plant communication may adversely affect the quality and/or consistency of the data being collected, unless individual plants are effectively “protected” from plant–plant communication. In an applied agricultural context, the results from this study are encouraging as they suggest that “hotspot” detection of emerging insect infestations in agricultural crops may be more likely if based on analyses of reflectance values in spectral bands. That is, if a small cluster of crop plants are infested by insects and emit a stress signal to adjacent plants, which in turn also respond to the stressor, then it may be easier to spatially localize the emerging infestation hotspot in crop fields. Thus, confirmation of leaf reflectance based detection of plant–plant communication is highly relevant to development of accurate and reliable methods in precision agriculture to determine when and where crop plants are under stress.
Karban R. The language of plant communication (and how it compares to animal communication). In: Gagliano M, Ryan JC, Vieira P, editors. The language of plants. Minneapolis: University of Minnesota Press; 2017. p. 3–26.
Raguso RA, Kessler A. Speaking in chemical tongues: decoding the language of plant volatiles. In: Gagliano M, Ryan JC, Vieira P, editors. The language of plants. Minneapolis: University of Minnesota Press; 2017. p. 27–61.
Stahl E, Hilfiker O, Reymond P. Plant–arthropod interactions: Who is the winner? Plant J. 2017;93:703–28.
Mithöfer A, Boland W. Recognition of herbivory-associated molecular patterns. Plant Physiol. 2008;146:825–31.
Santino A, Taurino M, Domenico S, Bonsegna S, Poltronieri P, Pastor V, Flors V. Jasmonate signaling in plant development and defense response to multiple (a)biotic stresses. Plant Cell Rep. 2013;32:1085–98.
Bosch M, Wright LP, Gershenzon J, Wasternack C, Hause B, Schaller A, Stintzi A. Jasmonic acid and its precursor 12-oxophytodienoic acid control different aspects of constitutive and induced herbivore defenses in tomato. Plant Physiol. 2014;166:396–410.
Grinberg-Yaari M, Alagarmalai J, Lewinsohn E, Perl-Treves R, Soroker V. Role of jasmonic acid signaling in tomato defense against broad mite, Polyphagotarsonemus latus (Acari: Tarsonemidae). Arthropod–Plant Interact. 2015;9:361–72.
Uehara T, Sugiyama S, Matsuura H, Arie T, Masuta C. Resistant and susceptible responses in tomato to cyst nematode are differentially regulated by salicylic acid. Plant Cell Physiol. 2010;51:1524–36.
Erb M, Meldau S, Howe GA. Role of phytohormones in insect-specific plant reactions. Trends Plant Sci. 2012;17:250–9.
Kanno H, Hasegawa M, Kodama O. Accumulation of salicylic acid, jasmonic acid and phytoalexins in rice, Oryza sativa, infested by the white-backed planthopper, Sogatella furcifera (Hemiptera: Delphacidae). Appl Entomol Zool. 2011;47:27–34.
Cheeseman JM. Hydrogen peroxide and plant stress: a challenging relationship. Plant Stress. 2007;1:4–15.
Chaouch S, Queval G, Vanderauwera S, Mhamdi A, Vandorpe M, Langlois-Meurinne M, Van Breusegem F, Saindrenan P, Noctor G. Peroxisomal hydrogen peroxide is coupled to biotic defense responses by isochorismate synthase1 in a daylength-related manner. Plant Physiol. 2010;153:1692–705.
Baxter A, Mittler R, Suzuki N. ROS as key players in plant stress signalling. J Exp Bot. 2014;65:1229–40.
Karban R. The ecology and evolution of induced resistance against herbivores. Functional Ecology. 2011;25:339–47.
Karban R, Baldwin IT, Baxter KJ, Laue G, Felton GW. Communication between plants: induced resistance in wild tobacco plants following clipping of neighboring sagebrush. Oecologia. 2000;125:66–71.
Maffei ME, Mithofer A, Boland W. Insects feeding on plants: rapid signals and responses preceding the induction of phytochemical release. Phytochemistry. 2007;68:2946–59.
Kuzniak E, Urbanek H. The involvement of hydrogen peroxide in plant responses to stresses. Acta Physiol Plant. 2000;22:195–203.
Estabrook EM, Yoder JI. Plant–plant communications: rhizosphere signaling between parasitic angiosperms and their hosts. Plant Physiol. 1998;116:1–7.
Dicke M, Bruin J. Chemical information transfer between plants: back to the future. Biochem Syst Ecol. 2001;29:981–94.
Hiltpold I, Turlings TC. Belowground chemical signaling in maize: when simplicity rhymes with efficiency. J Chem Ecol. 2008;34:628–35.
Song YY, Zeng RS, Xu JF, Li J, Shen X, Yihdego WG. Interplant communication of tomato plants through underground common mycorrhizal networks. PLoS ONE. 2010;5:e13324.
Babikova Z, Gilbert L, Bruce TJ, Birkett M, Caulfield JC, Woodcock C, Pickett JA, Johnson D. Underground signals carried through common mycelial networks warn neighbouring plants of aphid attack. Ecol Lett. 2013;16:835–43.
Bicchi C, Cordero C, Iori C, Rubiolo P, Sandra P. Headspace sorptive extraction (HSSE) in the headspace analysis of aromatic and medicinal plants. J Sep Sci. 2000;23:539–46.
Dewulf J, Langenhove HV. Analysis of volatile organic compounds using gas chromatography. Trends Anal Chem. 2002;21:637–42.
Ciccioli P, Brancaleoni E, Frattoni M, Maris C. Sampling of atmospheric volatile organic compounds (VOCs) with sorbent tubes and their analysis by GC-MS. In: Burden FR, Foerstner U, McKelvie ID, Guenther A, editors. Environmental monitoring handbook. New York: McGraw-Hill Publisher; 2002.
Balcke GU, Handrick V, Bergau N, Fichtner M, Henning H, Stellmach H, Tissier H, Hause B, Frolov A. An UPLC–MS/MS method for highly sensitive high-throughput analysis of phytohormones in plant tissues. Plant Methods. 2012;8:1–11.
Liu H, Li X, Xiao J, Wang S. A convenient method for simultaneous quantification of multiple phytohormones and metabolites: application in study of rice-bacterium interaction. Plant Methods. 2012;8:1–12.
Trapp MA, Souza GD, Rodrigues-Filho E, Boland W, Mithofer A. Validated method for phytohormone quantification in plants. Front Plant Sci. 2014;5:417.
Aljbory Z, Chen M-S. Indirect plant defense against insect herbivores: a review. Insect Sci. 2018;25:2–23.
Erb M, Ton J, Degenhardt J, Turlings TCJ. Interactions between arthropod-induced aboveground and belowground defenses in plants. Plant Physiol. 2008;146:867–74.
Frost CJ, Mescher MC, Carlson JE. Moraes CMd: Plant defense priming against herbivores: getting ready for a different battle. Plant Physiol. 2008;146:818–24.
Karban R, Yang LH, Edwards KF. Volatile communication between plants that affects herbivory: a meta-analysis. Ecol Lett. 2014;17:44–52.
Moreira X, Nell CS, Katsanis A, Rasmann S, Mooney KA. Herbivore specificity and the chemical basis of plant-plant communication in Baccharis salicifolia (Asteraceae). New Phytologist. 2016. https://doi.org/10.1111/nph.14164.
Morrell K, Kessler A. Plant communication in a widespread goldenrod: keeping herbivores on the move. Funct Ecol. 2017;31:1049–61.
War AR, Sharma HC, Paulraj MG, War MY, Ignacimuthu S. Herbivore induced plant volatiles: their role in plant defense for pest management. Plant Signal Behav. 2011;6:1973–8.
Furbank RT. Plant phenomics: from gene to form and function. Funct Plant Biol. 2009;36:v–vi.
Finkel E. With ‘phenomics’, plant scientists hope to shift breeding Into overdrive. Science. 2009;325:380–1.
Curran PJ. Remote sensing of foliar chemistry. Remote Sens Environ. 1989;30:271–8.
Kokaly RF, Asner GP, Ollinger SV, Martin ME, Wessman CA. Characterizing canopy biochemistry from imaging spectroscopy and its application to ecosystem studies. Remote Sens Environ. 2009;113:S78–91.
Siebke K, Ball MC. Non-destructive measurement of chlorophyll b:a ratios and identification of photosynthetic pathways in grasses by reflectance spectroscopy. Funct Plant Biol. 2009;36:857–66.
Scholes JD, Rolfe SA. Chlorophyll fluorescence imaging as tool for understanding the impact of fungal diseases on plant performance: a phenomics perspective. Funct Plant Biol. 2009;36:880–92.
Carter GA, Knapp AK. Leaf optical properties in higher plants: linking spectral characteristics to stress and chlorophyll concentration. Am J Bot. 2001;88:677–84.
Prabhakar M, Prasad YG, Rao MN. Remote sensing of biotic stress in crop plants and its applications for pest management. In: Venkateswarlu B, Shanker AK, Shanker C, Maheswari M, editors. Crop stress and its management: perspectives and strategies. New York: Springer; 2012. p. 517–49.
Nansen C, Elliot NC. Remote sensing and reflectance profiling in entomology. Annu Rev Entomol. 2016;61:139–58.
Carter GA, Miller RL. Early detection of plant stress by digital imaging within narrow stress-sensitive wavebands. Remote Sens Environ. 1994;50:295–302.
Carter GA. Responses of leaf spectral reflectance to plant stress. Am J Bot. 1993;80:239–43.
Severtson D, Callow N, Flower K, Neuhaus A, Olejnik M, Nansen C. Unmanned aerial vehicle canopy reflectance data detects potassium deficiency and green peach aphid susceptibility in canola. Precision Agric. 2016;17:659–77. https://doi.org/10.1007/s11119-016-9442-0.
Lacoste C, Nansen C, Thompson S, Moir-Barnetson L, Mian A, McNee M, Flower KC. Increased susceptibility to aphids of flowering wheat plants exposed to low temperatures. Environ Entomol. 2015;44(3):610–8. https://doi.org/10.1093/ee/nvv024.
Nansen C, Sidumo AJ, Martini X, Stefanova K, Roberts JD. Reflectance-based assessment of spider mite “bio-response” to maize leaves and plant potassium content in different irrigation regimes. Comput Electron Agric. 2013;97:21–6.
Ritchie SW, Hanway JJ. How a corn plant develops, vol. 48. Ames: IOWA State University; 1986. p. 21.
Bortolotto OC, Mikami AY, Bueno AF, Silva GV, Queiroz AP. Aspectos biológicos de Dichelops melacanthus em três temperaturas, alimentados com grãos imaturos de milho 2B688Hx e 2B688. Ciência Rural. 2016;46:254–9.
Souza CM, Baldin ELL, Ribeiro LP, Silva IF, Morando R, Bicalho KU, Vendramim JD, Fernandes JB. Lethal and growth inhibitory activities of Neotropical Annonaceae-derived extracts, commercial formulation, and an isolated acetogenin against Helicoverpa armigera. J Pest Sci. 2016;90:701–9.
Wang Y, Nansen C, Zhang Y. Integrative insect taxonomy based on morphology, mitochondrial DNA, and hyperspectral reflectance profiling. Zool J Linnean Soc. 2016. https://doi.org/10.1111/zoj.12367.
Nansen C, Singh K, Mian A, Allison BJ, Simmons CW. Using hyperspectral imaging to characterize consistency of coffee brands and their respective roasting classes. J Food Eng. 2016;190:34–9.
Voss SC, Magni P, Dadour I, Nansen C. Reflectance-based determination of age and species of blowfly puparia. Int J Legal Med. 2016;131(1):263–74. https://doi.org/10.1007/s00414-016-1458-5.
Luo C, Wei C, Nansen C. How do “mute” cicadas produce their calling songs? PLoS ONE. 2015;10:e0118554.
Nansen C, Coelho AJ, Mendes JV, Parra JRP. Reflectance-based identification of parasitized host eggs and adult Trichogramma specimens. J Exp Biol. 2014;217:1187–92.
Nansen C, Geremias LD, Xue Y, Huang F, Parra JR. Agricultural case studies of classification accuracy, spectral resolution, and model over-fitting. Appl Spectrosc. 2013;67:1332–8.
Hawkins DM. The problem of overfitting. J Chem Inf Comput Sci. 2004;1:1–12.
Defernez M, Kemsley EK. The use and misuse of chemometrics for treating classification problems. Trends Anal Chem. 1997;16:216–21.
Kemsley EK. Discriminant analysis of high-dimensional data: a comparison of principal components analysis and partial least squares data reduction methods. Chemomet Intell Lab Syst. 1996;33:47–61.
Lichtenthaler HK. Chlorophylls and carotenoids: pigments of photosynthetic biomembranes. Methods Enzymol. 1987;148:350–82.
Jessup W, Dean RT, Gebicki JM. Iodometric determination of hydroperoxides in lipid and proteins. Methods Enzymol. 1994;233:289–303. https://doi.org/10.1016/s0076-6879(94)33032-8.
Fisher RA. The use of multiple measurements in taxonomic problems. Ann Eugen. 1936;7:179–88.
Wasternack C, Hause B. Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. Ann Bot. 2007;2013(111):1021–58.
Paudel JR, Bede JC. Ethylene signaling modulates herbivore-induced defense responses in the model legume Medicago truncatula. Mol Plant–Microbe Interact. 2015;28:569–79.
Berens ML, Berry HM, Mine A, Argueso CT, Tsuda K. Evolution of hormone signaling networks in plant defense. Ann Rev Phytopathol. 2017;55:401–25. https://doi.org/10.1146/annurev-phyto-080516-035544.
Karban R. Plant sensing and communication. Chicago: University of Chicago Press; 2015.
Farmer EE, Ryan CA. Interplant communication: airborne methyl jasmonate induces synthesis of proteinase inhibitors in plant leaves. Proc Natl Acad Sci USA. 1990;87:7713–6.
Ankala A, Kelley RY, Rowe DE, Williams WP, Luthe DS. Foliar herbivory triggers local and long distance defense responses in maize. Plant Sci. 2013;199–200:103–12.
Paré PW, Tumlinson JH. Plant volatiles as a defense against insect herbivores. Plant Physiol. 1999;121:325–31.
Kessler A, Halitschke R, Diezel C, Baldwin IT. Priming of plant defense responses in nature by airborne signaling between Artemisia tridentata and Nicotiana attenuata. Oecologia. 2006;148:280–92.
Cheeseman JM. Hydrogen peroxide functions as a stress signal in plants. Bot Bull Acad Sin. 2005;46:1–10.
Ślesak I, Libik M, Karpinska B, Karpinski S, Miszalski Z. The role of hydrogen peroxide in regulation of plant metabolism and cellular signalling in response to environmental stresses. Acta Bioquimica Polonica. 2007;54:39–50.
Breusegem FV, Vranová E, Dat JF, Inzé D. The role of active oxygen species in plant signal transduction. Plant Sci. 2001;161:405–14.
Petrov VD, Van Breusegem F. Hydrogen peroxide—a central hub for information flow in plant cells. AoB Plants. 2012;2012:pls014.
Agarwal S, Sairam RK, Srivastava GC, Tyagi A, Meena RC. Role of ABA, salicylic acid, calcium and hydrogen peroxide on antioxidant enzymes induction in wheat seedlings. Plant Sci. 2005;169:559–70.
Quan LJ, Zhang B, Shi WW, Li HY. Hydrogen peroxide in plants: a versatile molecule of the reactive oxygen species network. J Integr Plant Biol. 2008;50:2–18.
Thenkabail PS, Smith RB. Pauw Ed: hyperspectral vegetation indices and their relationships with agricultural crop characteristics. Remote Sens Environ. 2000;71:158–82.
Thenkabail PS. Optimal hyperspectral narrowbands for discriminating agricultural crops. Remote Sens Rev. 2001;20:257–91.
Delalieux S, Aardt JANV, Zarco-Tejada PJ, Kempeneers P, Verstraeten WW, Coppin P. Development of robust hyperspectral indices for the detection of deviations of normal plant state. EARSeL e Proc. 2007;6:82–93.
Peñuelas J, Isla R, Filella I, Araus JL. Visible and near-infrared reflectance assessment of salinity effects on barley. Crop Sci. 1997;37:198–202.
Carroll MW, Glaser JA, Hellmich RL, Hunt TE, Sappington TW, Calvin D, Copenhaver K, Fridgen J. Use of spectral vegetation indices derived from airborne hyperspectral imagery for detection of European corn borer infestation in Iowa corn plots. J Econ Entomol. 2008;101:1614–23.
Coops NC, Stone C, Culvenor DS, Chisholm L. Assessment of crown condition in eucalypt vegetation by remotely sensed optical indices. J Environ Qual. 2004;33:956–64.
Zhu Y, Yao X, Tian Y, Liu X, Cao W. Analysis of common canopy vegetation indices for indicating leaf nitrogen accumulations in wheat and rice. Int J Appl Earth Observ Geoinform. 2008;10:1–10.
Cibula WG, Carter GA. Identification of a far-red reflectance response to ectomycorrhizae in slash pine. Int J Remote Sens. 1992;13:925–32.
Carter GA. Primary and secondary effects of water come on the spectral reflectance of leaves. Am J Bot. 1991;78:916–24.
Carter GA. Ratios of leaf reflectance in narrow wavebands as indicators of plant stress. Int J Remote Sens. 1994;15:697–703.
Riedell WE, Blackmer TM. Leaf reflectance spectra of cereal aphid-damaged wheat. Crop Sci. 1999;39:1835–40.
Lacoste C, Nansen C, Thompson S, Moir-Barnetson L, Mian A, McNee M, Flower KC. Increased susceptibility to aphids of flowering wheat plants exposed to low temperatures. Environ Entomol. 2015;44:610–8.
Conceived and design the experiments: LPR and CN. Performed the experiments: LPR, CN, ALSK, JAWF, MAT. Performed the chemical analysis: LPR, ALSK, MT, and AM. Analyzed the data: CN. Wrote the paper: CN, LPR, JAWF, MAT, MT, and AM. All authors read and approved the final manuscript.
The authors thank Carmen Cella dos Santos, Magda Alana Pompelli, Marcelo Rigo, and Zelinda Meneguzzi (Epagri/Cepaf) for technical assistance and insect rearing. We thank the Brazilian National Counsel of Technological and Scientific Development (CNPq, Process Number 445518/2014-6) for partial financial support. Marilia Almeida Trapp gratefully acknowledges financial support by a Capes-Humboldt Research Fellowship. We wish to thank two anonymous reviewers for very important and constructive feedback.
The authors declare that they have no competing interests.
Availability of data and materials
The data generated or analysed as part of this study are included in this published article or are available upon request from the corresponding author (especially raw hyperspectral images).
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.