Aboveground Herbivory Influences Belowground Defense Responses in Maize

The European corn borer (ECB; Ostrinia nubilalis) is an economically damaging insect pest of maize (Zea mays L.), an important cereal crop widely grown globally. Among inbred lines, the maize genotype Mp708 has shown resistance to diverse herbivorous insects, although several aspects of the defense mechanisms of Mp708 plants are yet to be explored. Here, the changes in root physiology arising from short-term feeding by ECB on the shoot tissues of Mp708 plants was evaluated directly using transcriptomics, and indirectly by monitoring changes in growth of western corn rootworm (WCR; Diabrotica virgifera virgifera) larvae. Mp708 defense responses negatively impacted both ECB and WCR larval weights, providing evidence for changes in root physiology in response to ECB feeding on shoot tissues. There was a significant downregulation of genes in the root tissues following short-term ECB feeding, including genes needed for direct defense (e.g., proteinase inhibitors and chitinases). Our transcriptomic analysis also revealed specific regulation of the genes involved in hormonal and metabolite pathways in the roots of Mp708 plants subjected to ECB herbivory. These data provide support for the long-distance signaling-mediated defense in Mp708 plants and suggest that altered metabolite profiles of roots in response to ECB feeding of shoots likely negatively impacted WCR growth.


INTRODUCTION
Plants have evolved various strategies to defend themselves against insect attack, mostly through direct and indirect defenses (War et al., 2018;Nalam et al., 2019). While feeding on host plants, insects secrete herbivore-associated molecular patterns (HAMPs), which are analogous to pathogen-associated molecular patterns (PAMPs), or microbe-associated molecular patterns in pathogens or microbes. Among HAMPs, oral secretions (OS), saliva, frass, and herbivoreassociated endosymbionts from chewing insects have been identified in influencing plant defenses (Basu et al., 2018). Plants presumably use surface immune receptors to recognize these insect HAMPs (Steinbrenner et al., 2020), and further induce a multitude of downstream defenses to curb the insect growth and/or colonization on host plants (Basu et al., 2018;Acevedo et al., 2019;Zogli et al., 2020).
The European corn borer (ECB; Ostrinia nubilalis), known as the "billion-dollar bug, " is a damaging pest of maize (Zea mays L.) and estimated to cost a billion dollars annually for plant treatment and lost crop yield (Hutchison et al., 2010). ECB management using transgenic maize plants was highly successful for several years, however, recently there are concerns regarding field-evolved ECB resistance to transgenic maize plants . ECB feeds on the aboveground parts of the maize plants and uses OS or saliva to modulate the plant defense responses. It was suggested that the indole-3-acetic acid (IAA) present in the OS of ECB may function as a HAMP in modulating maize defense responses (Dafoe et al., 2013). Another study showed that ECB saliva strongly induced direct defenses in maize (Louis et al., 2013). ECB saliva induced defense-related transcripts including lipoxygenase (LOX) and 12-phytodienoic acid (OPR) genes, which are involved in the biosynthesis of jasmonic acid (JA; Louis et al., 2013). It was also shown that ECB caterpillars utilized host-specific salivary factors in modulating plant defenses (Louis et al., 2013). In addition to ECB OS/saliva, ECB frass was also reported to modulate maize defense responses (Ray et al., 2016b).
Maize inbred line Mp708, which exhibits broad-based resistance to insect pests, was developed by crossing the insect resistant Mp704 and susceptible Tx601 plants (Williams et al., 1990). Several studies demonstrated that Mp708 provides resistance to both above and belowground feeding caterpillars (Pechan et al., 2000;Gill et al., 2011;Castano-Duque et al., 2017) as well as to piercing-sucking insects (Louis et al., 2015;Varsani et al., 2016Varsani et al., , 2019Pingault et al., 2021). Mp708's resistance to insect pests is mediated by the defense protein Maize Insect Resistance1-Cysteine Protease (Mir1-CP), which damages the peritrophic matrix of caterpillars and significantly reduces caterpillar growth (Pechan et al., 2000(Pechan et al., , 2002. Other defense-related proteins such as proteases, protease inhibitors, and peroxidases were also induced by caterpillar feeding on maize plants, which could also negatively impact caterpillar growth on Mp708 plants (Castano-Duque and Luthe, 2018). In addition, several other factors are also involved in providing maize resistance to chewing insects. For example, resistance to western corn rootworm (WCR; Diabrotica virgifera virgifera) can be influenced by root architecture (Branson et al., 1982), biomechanical strength (Meihls et al., 2012), or biochemical composition (Van Dam, 2009) of host plants. Among direct defense, plants can adapt their morphology by increasing trichome densities that create a physical barrier for insects (Tian et al., 2012). At the molecular level, chitinases, protease inhibitors, and related compounds are synthesized to counter insect invasion (Jamal et al., 2013;War et al., 2018).
Phytohormones also play a key role in modulating maize defense against chewing insects. Among them, phytohormones such as JA, salicylic acid (SA), and ethylene (ET) are the critical players in providing resistance to chewing pests (Erb et al., 2012). ECB feeding on maize foliage resulted in an increased accumulation of phytoalexins such as benzoxianoids and kauralexins (Dafoe et al., 2011). Interestingly, JA acts as a regulator of benzoxianoid and kauralexin biosynthesis (Oikawa et al., 2002), when acting in concert with ET. Whereas, JA can act synergistically with ET, SA and indole acetic acid can negatively regulate JA levels and actions (Koornneef and Pieterse, 2008). SA, JA, and abscisic acid have also been found in the fall armyworm (FAW; Spodoptera frugiperda) saliva (Acevedo et al., 2019). Application of these phytohormones on plants with the same concentrations detected in the caterpillar saliva induced the expression of mir1 in maize and lowered herbivore-induced defenses in tomato (Acevedo et al., 2019).
Aboveground to belowground communication is vital for plants to activate appropriate defenses against various stresses. When attacked by insect pests, plants activate defense responses both locally and systemically (Balmer et al., 2013). For example, several studies indicate that aboveground feeding by insect pests trigger belowground root-based defenses and vice-versa (Erb et al., 2009(Erb et al., , 2015Louis et al., 2015;Varsani et al., 2016). Similarly, caterpillar feeding on maize foliage can modulate long distance defense-related gene expression in roots (Ankala et al., 2013). Based on previous findings, we hypothesize that the aboveground feeding by ECB will trigger belowground defenses and modulate transcriptomic signatures in the belowground root tissues of the insect resistant Mp708 genotype. In this study, we used RNA-seq to investigate the dynamics of maize root transcriptome after short-term (24 h) aboveground feeding by ECB on maize foliage. We also evaluated how changes in root physiology induced by ECB herbivory affect larval growth of WCR, a major root-feeding insect pest of maize.

Plant and Insect Materials
Maize inbred lines, Mp708 and Tx601, used in this study were grown in soil mixed with vermiculite and perlite (PRO-MIX BX BIOFUNGICIDE + MYCORRHIZAE, Premier Tech Horticulture Ltd., Canada) in growth chambers with 14:10 (L:D) h photoperiod, 160 µE light m −2 s −1 , 25 • C, and 50-60% relative humidity. All plants for the experiments were used at the V2-V3 developmental stage (∼2 weeks; Ritchie et al., 1998). These plants were grown in 3.8 cm × 21.0 cm plastic Cone-tainers (Hummert International, Earth City, MO, United States). ECB and WCR were obtained from Benzon Research Inc., PA, and Crop Characteristics Inc., MN, respectively. ECB eggs obtained from Benzon Research Inc. were reared on artificial diet in growth chamber with 14:10 (L:D) h photoperiod at 23 • C as described previously (Louis et al., 2013). Non-diapausing WCR eggs were supplied by Crop Characteristics Inc. in 2015. WCRs were originally collected from susceptible non-Bt maize fields in Minnesota, United States, and this colony was maintained on non-Bt maize plants at Crop Characteristics Inc. for more than 80 generations (Scott Sargent, Crop Characteristics Inc. personal communication). Thus, these WCR populations were never exposed to selection for Bt resistance, and therefore, remain susceptible to Bt plants. WCR eggs obtained were kept in growth chamber with 14:10 (L:D) h photoperiod at 23 • C.
Larvae were hatched at 23 • C in growth chambers and were used for experiments.

Insect Bioassays
For insect bioassays, two newly molted 2nd instar ECB each weighing approximately 8-12 mg were introduced into the whorl region of maize plants. Individual plants were caged and after 5 days ECB larvae were collected. Only one ECB larva was able to be recovered from each maize plant due to ECB cannibalism, and the final weight individual larva was recorded. Similarly, to understand the interaction between ECB and WCR, two newly molted 5th instar ECB were introduced into the whorl region of Mp708 plants. After 24 h, ECB larvae were removed from the whorl region, and two newly hatched neonates of WCR were used to infest each plant. A fine-tipped paint brush (size O) was used to transfer neonate WCR larvae to maize plants as described previously (Jakka et al., 2016). Briefly, newly hatched WCR larvae were placed on the exposed nodal roots of maize plants and a thin layer of soil was used to cover the exposed roots. Bioassay plants were kept in a growth chamber with 14:10 (L:D) h photoperiod, 25 • C, and 50-60% relative humidity. These plants were watered as needed. After 12 days, the soil content of Cone-tainers were combined and the WCR were collected as described previously (Gill et al., 2011). Briefly, the soil content of all pants within a treatment was combined and sifted by hand until the collection of WCR was occasional. Subsequently, the soil was sifted through two mesh screens of different sizes (2.000 and 0.833 mm) to find the remaining WCR larvae in the soil. Finally, we ceased to monitor the presence of WCR in the soil content until none were obtained after 15-20 min of incessant searching. Mp708 plants that were uninfested with ECB on the whorl region were used as the control. Both bioassays were replicated twice with similar results and data from both independent experiments were combined and presented in Figure 1. The insect bioassay data were analyzed using mixed-model analysis of variance with experiment replication as the random effect. Analysis of variance (ANOVA) was performed using PROC GLIMMIX (SAS 9.3, SAS Institute, Cary, NC, United States), and means, when significant, were separated using the least significant difference procedure (P < 0.05).

Western-Blot Analysis
The hallmark of Mp708's defense is the accumulation of Mir1-CP protein both at the site of insect feeding or in systemic leaves or tissues away from the site of insect feeding (i.e., distal roots; Pechan et al., 2000;Louis et al., 2015;Varsani et al., 2016). To determine the presence of the Mir1-CP protein in Mp708 roots after aboveground feeding by ECB, we performed a western-blot analysis. Protein extraction and western-blot analysis were performed as described previously (Ankala et al., 2009;Varsani et al., 2016). Protein concentration was determined with Qubit reagent (Q33211) and Qubit 4 Fluorometer. A total of 60 µg of protein was used for electrophoresis. Upon gel electrophoresis, proteins were transferred to PVDF membrane using the iBlot 2 gel transfer FIGURE 1 | Mp708 provides enhanced resistance to ECB and WCR. (A) Larval weight of ECB feeding on Tx601 and Mp708 plants. Two newly molted 2nd instar ECB larvae each were introduced into the whorl region of maize plants. Individual plants were caged and after 5 days ECB larvae were collected, and final weight were recorded. (B) Aboveground feeding by ECB negatively impacts the growth of belowground feeding WCR. Two newly molted 5th instar ECB (+ECB) were introduced into the whorl region of Mp708 plants. After 24 h, ECB larvae were removed from the whorl region, and two newly hatched neonates of WCR were introduced. 12 days after the introduction of WCR, the fresh weight of WCR was measured. Mp708 plants that were uninfested with ECB (-ECB) on the whorl region were used as the control. For both (A,B), different letters above the bars indicate values that are significantly different from each other (P < 0.05). Error bars represent ± SE.
system. Images were finally captured with iBright FL1000 Imaging System.

RNA Extraction, RNA-seq Libraries Construction and Sequencing
At the V2-V3 stage of Mp708 plants, two newly molted 5th instar ECB were introduced into the whorl region. Individual plants were caged, and root samples were collected 24 h after infestation. ECB uninfested Mp708 plants were used as the control. In total, three biological replicates were produced for each treatment. Each biological replicate consisted of three individual maize tissues (root) pooled together. For each replicate, three ECBinfested and uninfested root tissues were collected and flashfrozen in liquid nitrogen. Maize Mp708 root tissues (80-100 mg) were ground using 2010 Geno/Grinder (SPEX SamplePrep, NJ, United States) for 40 s at 1,400 strokes min −1 . Total RNA was extracted from the homogenized tissue using Qiagen RNeasy Plant Mini Kit. Extracted total RNA was quantified through Nanodrop 2000c Spectrophotometer (Thermo Scientific TM). Then, RNA-seq libraries were constructed by using the mRNAseq standard TruSeq protocol from Illumina. RNA-seq libraries were sequenced in 50 bp paired-end with 20 million reads on average per library.

Analysis of RNA-Seq Libraries
The quality check of the RNA-seq libraries was performed with FASTQC (Andrews, 2010) and reads with a Phred score lower than 20 and length below 45 base pairs were removed with Trimmomatic v0.39 (Bolger et al., 2014). Then, trimmed reads were mapped on the maize reference genome v4 1 with Tophat2 (Kim et al., 2013) using the following parameters: 1 mismatch (-N 1), 0 splicing mismatch (-m 0), unique mapped reads (-g 1 -M). The transcripts reconstruction was performed with Cufflinks v2.2.1 with the following parameters: quantification against the reference annotation only (-G), multi-read-correct (-u) and fragbiais-correct (-b). The differentially expressed gene analysis was performed with Cuffdiff 2.2.1. Differentially expressed genes (DEGs) were identified with the following parameters: multipletesting corrected P-values ≤ 5% and expression ratios of | log 2 (Infested/Control)| ≥ log 2 (2). Gene ontology (GO) were analyzed with MaizeMine 2 by using the reference annotation as a template.

Functional Annotation
The GOBU package was used for enrichment calculations (Lin et al., 2006). The full set of maize gene annotation was used as the reference comparison set against down or upregulated DEGs. The P-values were calculated using Fisher's exact test and corrected for multiple testing with FDR method using the R module called "P-adjust."

Mp708 Provided Enhanced Resistance to European Corn Borer and Aboveground Feeding by European Corn Borer Negatively Impacted the Growth of Belowground Feeding Western Corn Rootworm
Previously, it was shown that Mp708 provides enhanced resistance to different feeding guilds of insect pests (Pechan et al., 2000;Lopez et al., 2007;Gill et al., 2011;Louis et al., 2015;Varsani et al., 2016;Pingault et al., 2021). Here, we examined whether Mp708 provides enhanced resistance to stem boring insect, ECB. Two independent experiments with similar experimental set up were performed. In the first and second bioassays,10 and 14 plants each of Tx601 and Mp708, respectively, were used and the combined data is presented in Figure 1A. Our results show that ECB growth was lower on the Mp708 plants compared to the maize inbred line Tx601 ( Figure 1A). As mentioned before, the hallmark of Mp708's defense is the accumulation of Mir1-CP at the site (i.e., leaf) and distal parts (i.e., roots) of the plant after insect infestation (Louis et al., 2015;Varsani et al., 2016). To determine if Mir1-CP accumulates in the roots after foliar feeding by ECB, we monitored the Mir1-CP levels in the roots of Mp708 plants using immunoblot analysis. Our results suggest an increased accumulation of Mir1-CP levels in the roots after ECB feeding in the aboveground parts of the plant for 24 h (Supplementary Figure 1). Next, we monitored whether aboveground feeding by ECB on Mp708 plants can induce resistance to subsequent herbivory by belowground feeding of WCR. In the first bioassay, out of 16 WCR larvae introduced on eight Mp708 plants, we recovered 12 and 13 WCR larvae on Mp708 plants that were preinfested with ECB (+ECB) and Mp708 plants that were uninfested with ECB (-ECB), respectively. Similarly, in the second bioassay, out of 24 WCR larvae introduced on 12 Mp708 plants, we recovered 17 and 19 WCR larvae on Mp708 plants that were preinfested with ECB (+ECB) and Mp708 plants that were uninfested with ECB (-ECB), respectively, and the fresh weight of each WCR was recorded. Here, the datasets from two independent experiments with similar experimental set up were combined pooled and the combined data is presented in Figure 1B. Our insect bioassay results indicate that aboveground feeding by ECB for 24 h provided enhanced resistance to subsequent herbivory by belowground feeding of WCR ( Figure 1B). Collectively, these data suggest that aboveground feeding by ECB can trigger belowground defenses in maize.

European Corn Borer Feeding Activated Hormone Pathways and Suppressed Biosynthesis of Secondary Metabolites in the Roots of Mp708 Plants
After 24 h of ECB feeding on the whorl region of the Mp708 plants, we identified 4,793 genes that were differentially expressed in Mp708 roots. Among these, 1,629 genes were upregulated, while 3,164 genes were downregulated (Figure 2 and Supplementary Table 1). The gene ontology enrichment was performed for each up and downregulated gene sets. The top ten gene ontology molecular functions terms enriched in upregulated DEGs were related to "protein kinase activity" (141 genes), "carbohydrate binding" (48 genes), "protein serine/threonine kinase activity" (97 genes), "Oxidoreductase activity, acting on paired donors, with incorporation, or reduction of molecular oxygen" (51 genes; Figure 2 and Supplementary Table 2). Molecular function terms enriched in downregulated genes included "peroxidase activity" (69 genes), "oxidoreductase activity" (294 genes), "antioxidant activity" (72 genes; Figure 2 and Supplementary Table 2).
The gene ontology biological process enrichment analysis revealed that downregulated genes are related with functions enriched in "response to oxidative stress, " "oxidation-reduction process, " or "carbohydrate metabolic process" while upregulated genes are related to functions as: "protein phosphorylation, " "regulation of nucleic acid-templated transcription, " "regulation of transcription, DNA-templated" (Supplementary Table 2).

Transcription Factors Involved in Phytohormone Pathways Were Activated in Roots in Response to European Corn Borer Feeding on Whorl Region of Mp708 Plants
Among the DEGs, 670 were transcription factors (TFs) and 18 families were represented by >1% of the TFs (Figure 3). Among the 18 families, 13 were composed of at least 50% of downregulated genes (Figure 3 Table 3). The large number of upregulated ERF transcription factors in Mp708 plants further confirm our previous observation that ET acts as a central node in mediating maize defense against various feeding guilds of insect pests (Louis et al., 2015).

Genes Related to Phytohormone Pathways Were Activated in Roots Upon European Corn Borer Feeding on Mp708 Shoots
To gain a better understanding for the role of phytohormones in defense against root feeders in response to aboveground herbivory, we investigated the DEGs involved in phytohormonerelated pathways. Thirty-one DEGs were found and divided between ET (11 DEGs), JA (12 DEGs), and SA (8 DEGs) pathways (Figure 4 and Supplementary Table 4). Among the 11 DEGs related to the ET pathway, six genes were upregulated: Ethylene-insensitive protein 2 (EIN2; Zm00001d039341), Ethylene-Insensitive-Like 8 (EIL8; Zm00001d016924), and four members of the 1-Aminocyclopropane-1-Carboxylic Acid Oxidase (ACO) family (ZmACO6: Zm00001d000163, ZmACO4: Zm00001d024852, ZmACO4: Zm00001d024853, and ZmACO20: Zm00001d052136). Interestingly, two gene members of the Acetyl-CoA synthetase (ACS) family were both downregulated (ACS1: Zm00001d039487 and ACS6: Zm00001d033862; Figure 4 and Supplementary Table 4). In contrast to ET, a majority of the DEGs involved in JA pathway were downregulated (11 genes) while only two genes were upregulated: JASMONATE RESISTANT 1, (JAR1; Zm00001d039345) and OPR5 (Zm00001d003584; Figure 4 and Supplementary Table 4). The three remaining genes of the OPR family, along with the six genes of the OPDA family and two genes of the ACS family, were downregulated (Figure 4 and Supplementary Table 4). All the DEGs involved in the SA pathway were downregulated in the roots of Mp708 plants after ECB feeding on the shoot region for 24 h (Figure 4 and Supplementary Table 4).

Altered Expression of Genes Involved in Plant Direct Defenses in Mp708 Roots After European Corn Borer Feeding on the Whorl Region
Mir1-CP, encoded by mir1, is characterized as a key component for the resistance to insect pests in maize (Pechan et al., 2000;Louis et al., 2015). Our results indicate that mir1 (Zm00001d036542) was upregulated in Mp708 roots at 24 hpi of ECB feeding on the whorl region with a high fold-change (FC; log 2 FC = 2.075). Protease inhibitors have also been identified in playing a role in direct plant defense mechanisms in response to insect attack (Koiwa et al., 1997). Here, one gene encoding a serpin, Zm00001d049148, was downregulated. Zm00001d002960, a member of Kunitz type inhibitor (KTI) family, was also downregulated. However, another KTI (smart serine protease inhibitor), WIN2 (Zm00001d015279), was upregulated. Three genes encoding Bowman-Birk type inhibitors were differentially expressed, among them two (Zm00001d003304 and Zm00001d040027) were downregulated and one, Zm00001d012939, was upregulated (Supplementary Table 1). Chitinases have also been identified as playing a role in plant defense and suppressing herbivore defenses (Ray et al., 2016a). Among the DEGs, 11 genes were functionally annotated as "chitinase family related, " "Chinase A, " "basic chitinase, " or "homolog of carrot EP3-3 chitinase." Among the 11 genes encoding chitinases, five were upregulated (Zm00001d036370, Zm00001d014840, Zm00001d043988, Zm00001d003190, and Zm00001d027525) in the roots of Mp708 plants after ECB feeding for 24 h on the shoot region, with a log 2 FC from 1.05 (Zm00001d027525) to 2.38 (Zm00001d036370), and six genes were downregulated (Zm00001d010911, Zm00001d037656, Zm00001d036366, Zm00001d048903, Zm00001d017152, and Zm00001d042769) with a log 2 FC from -1.39 to -3.06 (Supplementary Table 1).

Expression Pattern of Genes Involved in Metabolite Pathways That Modulate Defense in Maize Roots Against Root Feeding Herbivores
Root feeding insects are often attracted to metabolites released from the host plant, such as flavonoids, sugars, amino acids, or terpenes (Erb et al., 2013). Genes involved in these pathways were mined from the maize annotation using KEGG ID (Supplementary Table 5) and the expression profile is provided in Figure 5. One and ten genes were upregulated in the phenylpropanoid and flavonoid pathways, respectively, (Figure 5). Among the upregulated genes, one gene encoded a flavanone 3-hydroxylase (F3H, Zm00001d001960), two genes encoded chalcone and stilbene synthase family proteins (CHS, Zm00001d052673 and Zm00001d007403), and three genes encoded 2-oxoglutarate (2OG) and Fe(II)dependent oxygenase superfamily proteins (Zm00001d037487, Zm00001d042192, and Zm00001d051896), including DMR6 (DOWNY MILDEW RESISTANT 6, Zm00001d051896). DMR6 was proven to be involved in plant immunity for Arabidopsis (Van Damme et al., 2008). Zm00001d002287 (FLS2, Leucine-rich receptor-like protein kinase family protein) was downregulated with a log 2 FC (-6.04). FLS2 was described as a pattern recognition receptor that perceives the PAMP (Chinchilla et al., 2006). Three genes encoding 2-oxoglutarate (2OG) and Fe(II)-dependent oxygenase superfamily proteins were also downregulated (Zm00001d002564, Zm00001d012456, and Zm00001d047744). Three genes encoding UDP-glucosyl transferases were downregulated (Figure 5). Among the DEGs involved in amino acid metabolism, two genes were upregulated (Zm00001d028303, Zm00001d037455), and the remaining 13 genes were downregulated, as well as the seven and five DEGs involved in sugar and terpene pathways, respectively, (Figure 5). These data indicated a complex change in root metabolism within 24 h of aboveground feeding by ECB, with potential changes in the profiles and abundances of metabolites secreted by roots.

DISCUSSION
Plant response to insect herbivory varies from one genotype to another and from one pest to another. The maize genotype Mp708 has been investigated in the past for its response to FAW or WCR (Pechan et al., 2000;Castano-Duque et al., 2017). Here, we used RNA-seq to investigate the dynamics of Mp708 root transcriptome before and after ECB feeding (24 h) on the aboveground part of the plant. Previously, it was reported that short-term feeding by ECB did not impact ECB larval weight, suggesting a physiological adaptation of ECB to maize defense mechanisms (Dafoe et al., 2011). Here, the ECB larval weight was reduced after feeding on Mp708 compared to Tx601, which is susceptible to ECB, suggesting that Mp708 had resistance mechanisms not present in Tx601. Consistent with our earlier studies with corn leaf aphids (Varsani et al., 2016), our current work also documented an increased accumulation of the Mir1-CP protein in Mp708 roots following ECB herbivory on the whorl region of the plants.
At the transcriptional level, aboveground feeding by ECB for 24 h resulted in the downregulation of a majority of the DEGs in the roots of Mp708 plants. This suggests that ECB feeding resulted in the shutdown of plant functional mechanisms after ECB feeding for 24 h, including the genes involved in reactive oxygen species (ROS) metabolism. ROS and ROS metabolism play an important role in cell growth and development, as well as in the responses to biotic and abiotic stresses, and are well known modulators of plant defense (Koch et al., 2016). Both ROS and antioxidants respond to insect salivary compounds (Shinya et al., 2016). In this study, several genes encoding peroxidases were downregulated, suggesting that modulation of ROS was part of the early response of Mp708 plants to ECB herbivory. There was a significant positive enrichment of genes involved in protein phosphorylation, which might be important for plant defense signaling in Mp708 plants (Park et al., 2012;Jagodzik et al., 2018).
Compounds released by plant tissues can act as signals for insects to locate or avoid a host plant for herbivory (Erb et al., 2013). CO 2 emitted by plant root as well as sugars, terpenes, and amino acids have been characterized as attractants for herbivory but other compounds such as flavonoids can act as repellents for herbivory (Erb et al., 2013). Several genes encoding proteins associated with flavonoid biosynthesis were upregulated in roots of Mp708 plants challenged with ECB, and genes associated with sugar, amino acid, and terpene metabolism were downregulated. These data indicated that even short-term herbivory of Mp708 plants by ECB changed root physiology, possibly strengthening defense mechanisms against root attacking pests, such as the WCR. The negative effects of aboveground ECB challenge on WCR growth would support this contention.
Ethylene was also described as an attractant to WCR (Erb et al., 2013). In our transcriptomic data set, several genes encoding proteins of the ET pathway (Wang et al., 2002;Lacey and Binder, 2014) were upregulated. These included ET biosynthetic enzymes ACO and an ET receptor (EIN). Similarly, JA is another important hormone required for plant defense, and several proteins are required for its biosynthesis and perception of action. Additionally, JA regulates the transcription of many genes (Staswick and Tiryaki, 2004;Ankala et al., 2009;Afrin et al., 2015;Van Moerkercke et al., 2019;Yuan et al., 2019). ECB feeding on aboveground tissues downregulated a large majority of genes involved in JA pathway in the roots. However, JAR1 (Zm00001d039345; Staswick and Tiryaki, 2004) and COI (Zm00001d012338) were upregulated after 24 h of ECB feeding, suggesting that the JA-Ile formation was induced by ECBfeeding. It is known that JA-Ile binding to COI1-JAZ induces JA responsive genes (Pauwels et al., 2010;Zhang et al., 2017), such as mir1. Here, Mp708 showed an accumulation of mir1 transcripts in response to ECB feeding. Furthermore, it was previously shown in Mp708 that Mir1-CP accumulation was positively regulated by JA (Ankala et al., 2009).
Insect herbivores use diverse strategies to manipulate the JA-SA antagonism in plants to improve fitness outcomes. In many cases, herbivores use their salivary secretion to activate SA pathway and suppress JA-regulated defenses (Weech et al., 2008;Diezel et al., 2009). In Mp708 plants, ECB herbivory downregulated the genes in the SA pathway after 24 h of feeding, indicating that the defense response of Mp708 plants was not shifted by ECB feeding to increase SA responses. Potentially, this could be part of the resistance mechanism of Mp708 plants relative to other genetically related but susceptible maize lines. It has also been shown that ET not only had a synergistic interaction with JA (Penninckx et al., 1998) but was also a modulator of JA (Bostock, 2005;Von Dahl and Baldwin, 2007;Pieterse et al., 2009). ET can influence SA levels, but SA can antagonistically interact with JA. Thus, the crosstalk between ET, JA, and SA is a finely regulated process (Yang et al., 2019). In response to ECB feeding, early defense mechanisms were mostly associated with the increase of ET pathway genes, but the apparent antagonism of JA and SA was not induced by ECB feeding for 24 h. Future research will address how chewing insects influence ET, JA, and SA regulated resistance mechanisms in Mp708 plants compared to related susceptible lines.

CONCLUSION
European corn borer herbivory of aerial tissues in maize genotype Mp708 triggered large scale transcriptomic modulation in root tissues after 24 h of ECB introduction. Many genes associated with ROS mitigation were downregulated. However, several genes encoding protein kinases were upregulated, suggesting signaling cascades responsive to long-distance signals from the shoots. Transcriptomic analysis showed that genes encoding proteins associated with phytohormones ET and JA were activated after short-term ECB feeding. A similar upregulation of genes related to SA were not observed, indicating that ECB were not manipulating the JA-SA relationship in the root system. These transcriptomic data suggests that resistance mechanisms in Mp708 plants to chewing insects might be based on finely tuned regulation of plant hormonal pathways, and subsequent changes in plant metabolism.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI (accession: PRJNA684328).

AUTHOR CONTRIBUTIONS
JL conceived and designed the research. LP performed the computational analysis. SB and PZ performed the experiments. GS and NP contributed to methods development, reagents, and data analysis. WPW developed the maize genotype used in the study. LP, GS, and JL wrote the article with contributions from all authors. All authors have read and agreed to the final version of the manuscript.

FUNDING
This work was partially supported by the Nebraska Agricultural Experiment Station with funding from the Hatch Act (Accession #1007272) through the USDA National Institute of Food and Agriculture and funds from the University of Nebraska-Lincoln to JL. GS and NP were supported by the USDA-ARS CRIS project 3042-21000-034-00D.