Next Article in Journal
Copulation Duration and Sperm Precedence with Reference to Larval Diapause Induction in Monochamus alternatus Hope (Coleoptera: Cerambycidae)
Next Article in Special Issue
Resistance to Cassava Whitefly (Bemisia tabaci) among Eastern and Southern African Elite Cassava Genotypes
Previous Article in Journal
Mitogenomic Analysis and Phylogenetic Implications for the Deltocephaline Tribe Chiasmini (Hemiptera: Cicadellidae: Deltocephalinae)
Previous Article in Special Issue
Sublethal Effects of Spirotetramat, Cyantraniliprole, and Pymetrozine on Aphis gossypii (Hemiptera: Aphididae)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Omics Analysis Reveals Key Pathways in Cotton Defense against Mirid Bug (Adelphocoris suturalis Jakovlev) Feeding

1
National Key Laboratory of Cotton Bio-Breeding and Integrated Utilization, Chinese Academy of Agricultural Sciences, No. 38, Huanghe Road, Anyang 455000, China
2
Hubei Insect Resources Utilization and Sustainable Pest Management Key Laboratory, College of Plant, Science and Technology, Huazhong Agricultural University, Wuhan 430070, China
3
Green Agricultural Products Safety and Warning Laboratory, Research Center of Soil Resource Comprehensive Utilization and Ecological Environment in Western Inner Mongolia, Hetao College, Bayannur 015000, China
4
Henan Institute of Science and Technology, College of Life Science, Hualan St. 90, Xinxiang 453003, China
5
Anhui Provincial Center for Disease Control and Prevention, Hefei 230601, China
*
Authors to whom correspondence should be addressed.
Insects 2024, 15(4), 254; https://doi.org/10.3390/insects15040254
Submission received: 14 February 2024 / Revised: 5 April 2024 / Accepted: 6 April 2024 / Published: 8 April 2024
(This article belongs to the Collection Integrated Pest Management of Crop)

Abstract

:

Simple Summary

This study delves into the intricate defense mechanisms of cotton plants against the primary pest Adelphocoris suturalis Jakovlev, which has become a major concern in Bt-cotton fields. Through advanced analytical techniques, this study identified specific changes in the cotton plant’s metabolism and protein expression triggered by A. suturalis feeding. These findings revealed a notable increase in the expression of certain α-linolenic acid metabolism pathway-related proteases and a simultaneous decrease in fructose and mannose biosynthesis-related proteases. These molecular responses shed light on the complex interplay between cotton plants and A. suturalis. Therefore, this study not only deepens our understanding of plant–insect interactions but also provides valuable insights for developing innovative strategies to control this cotton pest, thereby offering potential solutions for sustainable cotton cultivation.

Abstract

The recent dominance of Adelphocoris suturalis Jakovlev as the primary cotton field pest in Bt-cotton-cultivated areas has generated significant interest in cotton pest control research. This study addresses the limited understanding of cotton defense mechanisms triggered by A. suturalis feeding. Utilizing LC-QTOF-MS, we analyzed cotton metabolomic changes induced by A. suturalis, and identified 496 differential positive ions (374 upregulated, 122 downregulated) across 11 categories, such as terpenoids, alkaloids, phenylpropanoids, flavonoids, isoflavones, etc. Subsequent iTRAQ-LC-MS/MS analysis of the cotton proteome revealed 1569 differential proteins enriched in 35 metabolic pathways. Integrated metabolome and proteome analysis highlighted significant upregulation of 17 (89%) proteases in the α-linolenic acid (ALA) metabolism pathway, concomitant with a significant increase in 14 (88%) associated metabolites. Conversely, 19 (73%) proteases in the fructose and mannose biosynthesis pathway were downregulated, with 7 (27%) upregulated proteases corresponding to the downregulation of 8 pathway-associated metabolites. Expression analysis of key regulators in the ALA pathway, including allene oxidase synthase (AOS), phospholipase A (PLA), allene oxidative cyclase (AOC), and 12-oxophytodienoate reductase3 (OPR3), demonstrated significant responses to A. suturalis feeding. Finally, this study pioneers the exploration of molecular mechanisms in the plant–insect relationship, thereby offering insights into potential novel control strategies against this cotton pest.

1. Introduction

Higher plants successfully colonize herbivore-rich environments by developing complex anti-herbivore feeding mechanisms [1]. When plants detect physical and chemical cues from insects, such as oral secretions and oviposition secretions, they can modify their proteins and metabolites [2]. Plants primarily resist herbivorous insects through jasmonate (JA), ethylene, and salicylic acid [3,4]. Specifically, the JA signaling pathway is crucial for this plant resistance [5,6]. Jasmonates (JAs), belonging to the oxylipin family, are key players in plant defense [7]. Initiation with α-linolenic acid (ALA) in the JA biosynthesis pathway produces key anabolic products like JA and its derivatives (e.g., methyl jasmonate, isojasmonic acid, methyl isojasmonate, etc.). In the ALA metabolism pathway, allene oxide synthase (AOS), allene oxide cyclase (AOC), and 12-oxophytodienoate reductase3 (OPR3) are crucial for JA synthesis [8,9,10]. ALA released as a stress signal through lipase activity on chloroplast membranes serves as the substrate for various oxidized lipids, including methyl jasmonate (MeJA) and other JA metabolites [11]. Increasing evidence suggests that the JA signaling pathway is crucial in plant defense, contributing to rice resistance against the leaf folder (Cnaphalocrocis medinalis) [12]. JA also mediates cotton plant resistance induced by Helicoverpa armigera and affects the growth and development of Aphis gossypii [13]. When rice plants are attacked by the brown plant hopper (Nilaparvata lugens), JA activates rice defense responses and gibberellin catabolism [14]. Phospholipase A (PLA), AOS, AOC, and OPR3, the key enzymes in the JA biosynthesis pathway, have been well reported in enhancing plant defense [9,14,15,16,17]. AOS and its family play a vital role in determining the resistance of rice to chewing and phloem-feeding herbivores [14,17].
Cotton (Gossypium hirsutum L.), a significant industrial crop worth 50 billion USD globally [18], has seen widespread adoption of transgenic Bt (Bacillus thuringiensis) cotton, reducing the need for broad-spectrum insecticides in China. However, this shift has resulted in frequent mirid outbreaks [19,20,21]. Adelphocoris suturalis, a mirid bug, has emerged as a top pest in the primary cotton-growing area. Both nymphs and adults of A. suturalis feed on cotton buds, tender shoots, leading to abscission, wilting, abnormal growth, and losses in lint production and quality [22]. These mirid bugs may eventually threaten the future of Bt cotton. Biological control methods, including transgenic plants, entomopathogenic microbes or a combination of both, have shown promising efficacy in pest management [23,24,25,26]. Many studies have focused on A. suturalis ecology and physiology [27,28,29,30]. However, current research on how A. suturalis responds to the mechanism of resistance defense in cotton is still lacking. Therefore, study of key pathways and key regulators will enrich our knowledge about plant defense and can help develop novel strategies to control mirid bugs, and most importantly provide the theoretical basis for transgenic breeding techniques in cultivating cotton against insect pests.
Detection and quantification of metabolites are achieved with techniques like quadrupole time-of-flight mass spectrometry (LC-Q/TOF-MS)-based metabolomics [31,32,33], and for the detection and quantification of proteins, quantitative isobaric tag for relative and absolute quantitation–liquid chromatography–mass spectrometry/mass spectrometry (iTRAQ-LC–MS/MS) proteomics is used. This sensitive method facilitates high-throughput protein identification and quantification [34], with applications in insects and plants like silkworm [35], brown planthopper [36], pine beetle [37], locust [38], and cotton [39]. These techniques provide a comprehensive view of dynamic proteomic and metabolomic variations, thereby enabling the identification of crucial proteins and metabolites essential for plant defenses.
The A. suturalis–cotton interaction forms the scientific basis for this pest control. Our study delves into the molecular-level intricacies of the interaction mechanism between cotton and A. suturalis. The specific objectives of this study were to identify insect-resistant genes and leverage transgenic technology in cotton for effective A. suturalis control. Understanding cotton’s insect resistance mechanism post-A. suturalis feeding holds great significance. This research aims to improve the theoretical framework of integrated pest management for A. suturalis. However, limited reports exist on cotton protein and metabolite changes after A. suturalis feeding [27,28]. Therefore, in this study, we employed proteomic and metabolomic approaches to investigate cotton plant defenses in response to A. suturalis feeding, explored changes, and identified potential molecular mechanisms.

2. Materials and Methods

2.1. Plant Material and Insect Infestation

Potted cotton plants (G. hirsutum CCRI 49) were transplanted into soil pots (15 × 17 cm), then put it into environmentally controlled chambers (50 × 60 × 110 cm) and maintained under either a 14 h (long day) or 10 h (short day) light period at 25 ± 2 °C and 65–70% relative humidity till the plants reached the four-leaf stage for further testing [13]. Subsequently, each cotton plant was placed in a 30 cm squared cage, and nine adult A. suturalis (starved for 12 h before inoculation) were introduced onto the leaves. A preliminary experiment showed that a 24–48 h timeframe is crucial for medium-term resistant defense. After 48 h of insect infestation, all insects were promptly removed. Cotton samples were immediately flash-frozen in liquid nitrogen and stored at −80 °C for further processing.

2.2. Metabolomic Sequencing

The primary aim of metabolomic analysis is to identify statistically significant and biologically relevant metabolites from the extensive pool of detected metabolites, followed by elucidating associated metabolic processes. LC–MS-based targeted metabolomic analysis was carried out as per a previous study [40,41]. Progenesis QI (ver. 2.2) software analyzed the MS data for peak (mz) retention time and ion area. The normalized data were then entered into SIMCA-P V11.0 (Umetrics, Umeå, Sweden) for PCA and OPLS-DA analysis. Univariate analysis was performed using t-tests, and the p-values from the t-tests after false-discovery rate (FDR) correction then produced q-values [41,42]. The results were significant when the q-value was <0.05. OPLS-DA was used to investigate the metabolite pattern changes. The differential metabolites were then selected when VIP values from the OPLS-DA model exceeded 1. FC ≥ 1.2 or ≤0.8 was used to determine differential metabolite variance between groups [41,42,43,44]. The importance of differential metabolites is connected to the FC value and VIP score, and the larger the value, the higher the score and the more important the substance is. Each metabolite’s related pathways were also listed by searching the KEGG (http://www.genome.jp/KEGG/, accessed on 16 June 2017) pathway database, while the metabolite molecular formulae of the matched metabolites were further identified by isotopic distribution measurement. Metabolomic sequencing included ten experimental replicates for each treatment, with group 1 (M48C) indicating uninfested controls grown for 48 h and group 2 (M48E) representing plants infested with insects for 48 h.

2.3. iTRAQ Sequencing

iTRAQ sequencing was conducted at BGI (Shenzhen, China). Using a reported phenol extraction method, cotton plant total proteins were extracted [45,46,47]. Protein concentrations were measured using the Bradford method [48]. The experiment comprised three independent biological replicates, each utilizing 25 μg total protein. Sequencing-grade trypsin (Promega) was used to digest the proteins at a 1:10 (w:w) ratio for 12 h at 37 °C. Subsequently, the samples were labeled using iTRAQ 8-plex kits (AB Sciex Inc., Foster, CA, USA) according to the manufacturer’s instructions. The samples were labeled with P48E (insect feeding, 48 h treatment) and P48C (control without insect feeding, 48 h). LC-MS/MS analysis was carried out as described previously [45,49], on an AB SCIEX Triple TOF 5600 system (Foster, CA, USA) with a Nanospray III source (AB SCIEX, Concord, ON, USA) featuring a pulled quartz tip as the emitter (New Objectives, Woburn, MA, USA). The MS operated with an RP ≥ 30,000 FWHM for TOF MS scans. Mascot version 2.3.02 search engine by Matrix Science, Boston, MA, USA, was used for protein identification and quantitative analysis, as previously reported [50,51,52]. The ExPASy website provided by the Expert Protein Analysis System was utilized for protein sequence prediction and retrieval, as reported previously [53]. Differentially upregulated and downregulated proteins were considered significant if their FC values were ≥1.2 or ≤0.8; q-values obtained through FDR correction of p-values indicated statistical significance at q < 0.05 [39,45,46]. Blast2GO was used to annotate protein functions against the non-redundant NR protein database on NCBI. The KEGG (http://www.genome.jp/KEGG/, accessed on 16 June 2017) and GO (http://geneontology.org/, accessed on 5 July 2017) databases were used for categorizing and grouping identified proteins. For proteomic sequencing, each treatment was performed in triplicate.

2.4. Quantitative Real-Time PCR (qRT-PCR)

For qRT-PCR analysis, six proteins were selected. These proteins included allene oxidase synthase (AOS), phospholipase A (PLA), alkene oxidative cyclase (AOC), alcohol dehydrogenase (ADH), and two 12-oxophytodienoate reductase 3 (OPR3-1 and OPR3-2), and experiments were conducted corresponding to 3 h, 6 h, 12 h, 24 h, 48 h, and 72 h of insect infestation, respectively. Control plants not subjected to infestation, were grown for the same duration. The leaves were promptly frozen in liquid nitrogen and stored at −80 °C until further use.
Total RNA was extracted using the RNAprep Pure Plant Kit (TIANGEN Biotech, BeiJing, China) following the manufacturer’s protocol [28]. First, 1 μg RNA was reverse-transcribed to cDNA via a PrimeScriptTM RT reagent kit (perfect real time) (TaKaRa, Dalian, China) as per the manufacturer’s protocol. The RT-qPCR, conducted on an Eppendorf Mastercycler eprealplex 2.2 (Eppendorf, Hamburg, Germany), utilized GoTaq Qpcr Master Mix (Promega, Madison, WI, USA) with three biological and three technical replicates. The thermal cycle conditions for qRT-PCR were 95 °C for 2 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. Relative gene expression was analyzed through the 2−△△Ct method [54]. The mean Ct value of three technical replications in each sample was used for the gene expression analysis, and the relative expression level of each treatment was the average of the three biological replications, and incorporated two housekeeping genes (GhHis3; GenBank accession AF024716, and Gossypium hirsutum ubiquitin 7 (GhUB7, GenBank accession DQ116441.1) according to the literature [55]. The three-hour control group served as the reference sample for data normalization. All qRT-PCR primer pairs were designed using the online Primer Quest Tool (http://sg.idtdna.com/PrimerQuest/Home/Index, accessed on 12 August 2017) and are detailed in Table S1. Significance in expression level differences was determined by one-way ANOVA, and means were separated using Tukey’s HSD (IBM SPSS Statistics v 22.0).

2.5. Data Analysis

Data analysis and chart preparation were conducted using Microsoft Excel 15.37, Mac Preview (8.1), Adobe Photoshop (2020), and GraphPad Prism 6. Data were analyzed using IBM SPSS Statistics 22.

3. Results

3.1. Metabolite Identification

3.1.1. Identification of Differential Metabolites

A comprehensive metabolomic analysis was conducted on cotton plants infested by A. suturalis for 48 h. The initial mass spectrometry data comprised 7790 positive ions and 4325 negative ions. Subsequent PCA analysis of quality control (QC) samples revealed that PCA1 contributed 32.32% and PCA2 contributed 16.78%, suggesting their suitability for PCA. Importantly, QC samples exhibited clustering, affirming instrument stability and data reliability (Figure S1). In this study, both positive ion (ESI+) and negative ion (ESI-) modes collectively contributed over 50% (ESI+ = 57.08%, ESI− = 52.14%). This indicates that these principal components effectively represent and interpret the significant metabolite changes induced by A. suturalis feeding on cotton. Furthermore, PLS-DA analysis was employed to maximally highlight differences among test groups. The ESI+ model parameters were R2 = 0.86, Q2 = 0.62, while ESI− model parameters were R2 = 0.90, Q2 = 0.69. The minimal difference between R2 and Q2 values (0.2–0.3) suggests the reliability and effectiveness of the current PLS-DA model [56] (Figure S1). Under this standard (VIP ≥ 1, FC ≥ 1.2 or ≤0.8, and q < 0.05), there were 496 differential positive ions (122 downregulated, 374 upregulated) and 53 differential negative ions (7 downregulated, 46 upregulated) (Table S2). A total of 496 cationic metabolites categorized under 11 types were detected in the cotton metabolome (Figure 1A, Tables S3 and S4), including terpenoids, alkaloids, phenylpropanoids, flavonoids or isoflavones, organic acids and their derivatives, phenols and their derivatives, carbohydrates and carbohydrate conjugates, lipids and lipid molecules, amino acid peptides and their analogues, nucleotides and their analogues, etc. These results suggested that metabolites were upregulated, which may induce cotton’s defense against herbivorous insects. All differential cationic metabolites (496) and additional analytical details are listed in the Supplementary Materials (Table S5).
Our results showed that the cotton plant responded to the feeding stress of A. suturalis by increasing its secondary metabolites. These metabolites, including terpenes, alkaloids, phenylpropanoids, and flavonoids, were most significantly upregulated post-A. suturalis feeding (Figure 1B, Table S4). The terpenes germacrene D (q = 0.043) and alpha-barbatene (q = 0.043) showed the highest upregulation, with FC up to 46.8 and VIP up to 3.6. Among alkaloid metabolites, betaine aldehyde (q = 0.043) and acetylcholine (q = 0.036) were upregulated, reaching FCs of 46.8 and 29.4, and VIP scores of 3.6 and 4.0, respectively. The phenylpropanoid metabolites, including coumarin (q = 0.002) and cinnamaldehyde (q = 0.020). showed overall upregulation, reaching FCs of 5.8 and 5.3 and VIP scores of 3.2 and 2.4, respectively. All secondary flavonoids were significantly upregulated, with ononin (q = 0.043) and syringetin (q = 0.0002) exhibiting FCs up to 5.5 and 5.4 and VIP scores up to 2.3 and 3.3, respectively (Figure 1B,C, Table S4).

3.1.2. KEGG Annotation and Functional Analysis of Differential Metabolites

The KEGG database was used to annotate 496 differential metabolites into 65 pathways after cotton consumption compared to the control group. The pathways included the sesquiterpene and triterpene biosynthetic pathway (map00909), arachidonic acid metabolic pathway (map00590), tryptophan metabolic pathway (map00380), etc. (Table 1 and Table S6).

3.1.3. The Two Key Metabolic Pathways

Plant resistance is multifaceted. For example, plants can reduce their own nutrients to affect the growth and development of insects, achieving insect resistance [57]. In addition, α-linolenic acid acts as a stress signal released through lipase activity on chloroplast membranes and serves as the substrate for various oxidized lipids, including methyl jasmonate (MeJA) and other JA metabolites [7]. Therefore, we chose to comprehensively analyze two pathways of interest: the fructose and mannose metabolism pathway and the ALA metabolism pathway. The fructose and mannose metabolism pathway showed six upregulated metabolites (q = 0.028) and eight downregulated metabolites (q = 0.032) (Figure 2A,B), while the ALA metabolism pathway exhibited fourteen upregulated metabolites (q < 0.047) and two downregulated metabolites (q < 0.045) (Figure 2C,D, Table S7).

3.2. iTRAQ Identification of Proteins

Cotton proteomic analysis using iTRAQ involved three replicates, as previously described [39,45,47]. A total of 371,575 and 374,246, spectra were obtained by P48E treatment and P48C control, including 52,750 and 48,103 unique spectra and 17,865 and 16,116 peptide segments, respectively. Furthermore, 12,906 and 11,903 unique peptide proteins were obtained, and 5520 and 5210 unique peptide proteins were identified in P48E and P48C, respectively. Overall, a total of 8302 proteins were identified from the three samples (Table S8).

3.2.1. Functional Categories of Differentially Expressed Proteins

Of all 8302 proteins, there are a total of 1252 differentially expressed proteins (q < 0.05, 775 upregulated and 477 downregulated, Table S9) in the P48E vs. P48C control comparison (Figure 3A). These findings indicate that plants adjust protein levels in response to A. suturalis attacks, activating plant defense responses. GO enrichment analysis of differentially expressed proteins (DEPs) for P48E and its P48C control showed them being classified into three groups: biological process, cellular component, and molecular function (Figure 3B, Table S10).

3.2.2. Pathway Enrichment of Differential Proteins

Furthermore, the metabolic pathway enrichment analysis of cotton post-feeding showed the sequencing results for feeding treatment (P48E) and control (P48C). The analysis identified 1569 differential proteins enriched in 35 metabolic pathways (q < 0.05) (Table S11), with the top 20 pathways depicted in Figure 4A. Subsequent examination indicated a higher number of upregulated proteins in defense-related metabolic pathways compared to downregulated ones (Figure 4B), exemplified by the phenylpropanoid biosynthesis pathway (30 upregulated proteins, 15 downregulated), ALA metabolism pathway (17 upregulated and 2 downregulated), and fructose and mannose metabolism pathways (7 upregulated and 19 downregulated).

3.3. Integrative Metabolomic–Proteomic Analysis

The proteins and corresponding metabolites were considered to be correlated if they were both expressed at the same metabolic pathway. Integrated analysis of the proteome and metabolome revealed significant changes in metabolic pathways in cotton after A. suturalis feeding, including phenylpropanoid biosynthesis, amino sugar and nucleotide sugar metabolism, photosynthesis pathway, fructose and mannose metabolism pathway, glutathione metabolic pathway, pentose phosphate pathway, alpha-linolenic acid metabolism, etc. (Figure 5). Therefore, fructose and mannose metabolism, as well as ALA metabolism, were selected for further analysis.

3.3.1. Integrating Proteomics and Metabolomics for Fructose and Mannose Metabolism Pathway Analysis

In the integrative proteo-metabolomic analysis of fructose and mannose metabolic pathways, we identified 7 upregulated and 19 downregulated proteases (Figure 6, Table S12). The upregulated proteases include phosphofructokinase (q = 0.001), epidermis-specific secreted glycoprotein EP1 (q = 0.038), xylose isomerase (q = 0.001), idonate-5-dehydrogenase (q = 0.002), and fructose-bisphosphate aldolase3 (q = 0.001). The downregulated proteases comprise fructose-bisphosphate aldolase (q = 0.043), triosephosphate isomerase (q = 0.001), fructokinase (q = 0.001), endo-1,4-beta-mannanase (q = 0.017), hexokinase (q = 0.001), triosephosphate isomerase (q = 0.001), and fructokinase-1-6-bisphosphatase (q = 0.001), etc. These differential proteases reached significant levels (q < 0.05), with some at extremely significant levels (q < 0.01).
The metabolomic data indicated six upregulated and eight downregulated metabolites (Figure 6, Table S7). The upregulated metabolites include D-allose (q = 0.028), D-mannose (q = 0.028), L-sorbose (q = 0.028), L-rhamnonate (q = 0.028), L-fuconate (q = 0.028), and D-fructose (q = 0.028). Downregulated metabolites comprised D-allose 6-phosphate (q = 0.032), D-allulose 6-phosphate (q = 0.032), D-mannose 1-phosphate (q = 0.032), D-mannose 6-phosphate (q = 0.032), D-fructose 1-phosphate (q = 0.032), beta-D-fructose 2-phosphate (q = 0.032), sorbose-1-phosphate (q = 0.032), and GDP-4-dehydro-6-deoxy-D-mannose (q = 0.024) (Figure 2 and Figure 6). These differential metabolites reached significant levels (q < 0.05), with some being extremely significant (q < 0.01).
The observed changes in both proteins (seven upregulated and nineteen downregulated) and metabolites (six upregulated and eight downregulated) strongly suggest that plants employ a defense strategy against insect predation by reducing their own nutrient consumption. This subsequently affects the growth and development of A. suturalis, and this is consistent with what was observed during the experiment (e.g., development process is prolonged).

3.3.2. Integrating Proteomics and Metabolomics for ALA Metabolism Pathway Analysis

In the integrative metabolomic–proteomic analysis of the ALA metabolism pathway, 17 upregulated and 2 downregulated proteases were identified (Figure 7, Table S12). The upregulated proteases included salicylate O-methyltransferase (q = 0.032), AOS (q = 0.021–0.038), 3-ketoacyl-CoA thiolase 2 (KAT2; q = 0.040), phospholipase A1-Iigamma (PLA; q = 0.001), 4-coumarate--CoA ligase-like 5 (q = 0.001), acyl-coenzyme A oxidase 2 (ACX2; q = 0.003), acyl-coenzyme A oxidase 4 (ACX4; q = 0.011), OPR3 (q = 0.0110), alkene oxidative cyclase 4 (AOC; q = 0.001), et al. The downregulated proteases included alcohol dehydrogenase class-P (ADH; q = 0.001), and phospholipase A1-IIdelta (PLA; q = 0.011). These differential proteases reached significant levels (q < 0.05), with some at extremely significant levels (q < 0.01).
Metabolomic data for the ALA metabolism pathway showed 14 upregulated and 2 downregulated metabolites (Figure 7, Table S7). The upregulated metabolites included 3-hexenal (q = 0.029); 12,13-EOT (q = 0.022), 9,10-EOT (q = 0.022), colnelenic acid (q = 0.022), etherolenic acid (q = 0.022), 10-OPDA (q = 0.022), 12-OPDA (q = 0.022), JA (q = 0.018), 7-isojasmonic acid (q = 0.018), stearidonic acid (q = 0.047), 3,6-nonadienal (q = 0.041), 9-oxononanoic acid (q = 0.038), MeJA (q = 0.032), and 7-isomethyljasmonate (q = 0.032). The downregulated metabolites were 3-hexenol (q = 0.045) and volicitin (q = 0.031). These differential metabolites reached significant levels (q < 0.05).
Based on the changes seen in both proteins (17 upregulated and 2 downregulated) or metabolites (14 upregulated and 2 downregulated), we can conclude that plants resist insect predation by increasing the production of defense regulators like proteins and metabolites to defend against A. suturalis feeding.

3.4. Analysis of Key ALA Metabolism Pathway Regulators

Upon analysis, significant alterations (upregulated) in the ALA metabolism pathway were observed in both proteomic and metabolomic analyses. The integrated proteomic and metabolomic analysis revealed predominantly downregulated metabolites, except for volicitin and 3-hexenol, whereas most related proteins were upregulated. These proteins included AOS, PLA, AOC, and OPR3, etc. The outcome suggests an increase in upstream regulator proteins, leading to a decrease in downstream metabolites within the ALA metabolism pathway. This protein and metabolite reconfiguration appears to be a plant response to insect-feeding stress, indicating the potential significance of the ALA pathway in cotton defense against herbivores. Key regulatory proteins included AOS, PLA, AOC, and OPR3.
To verify that protease regulators respond to mirid bug feeding, further investigation focused on six selected proteases from the ALA metabolic pathway, comprising four upregulated proteins (Table S12), i.e., AOS (CotAD_35840), PLA (CotAD_52791) and AOC (Cotton_D_gene_10007846), OPR3 (OPR3-1; CotAD_59461, and OPR3-2; Cotton_D_gene_10037325) and one downregulated protein, i.e., alcohol dehydrogenase ADH (CotAD_64340). These proteins underwent qRT-PCR analysis to confirm expression at the transcription level, validating the sequencing results (Figure 8).
The results showed that phospholipase A (PLA) gene expression increased progressively with feeding time, finally peaking at 48 h before decreasing significantly differently from the control group (Figure 8A, p < 0.01). Allene oxidase synthase (AOS) gene expression was highest at 12 h, showing a significant difference from the control group (Figure 8B, p < 0.01). Allene oxide cyclase (AOC) gene expression in cotton leaves gradually increased during A. suturalis feeding, peaking at 24 h, and then declined. AOC expression at 6 h, 12 h, 24 h, 48 h, and 72 h post-insect feeding significantly differed from the control group (Figure 8C, p < 0.01). These results suggest AOC has a relatively faster response in upper plant parts to A. suturalis feeding. OPR3-1 and OPR3-2 displayed similar expression levels, with both reaching significant differences (p < 0.05) at 12 h (p < 0.05) and 24 h (p < 0.01) compared to the control group (Figure 8D,E). ADH was downregulated, reaching significant differences at 3 h, 6 h and 24 h (Figure 8F, p < 0.01).
Through dynamic time variation analysis (Figure 8G), we conclude that PLA expression gradually increases, peaking at 48 h and subsequently decreasing. AOC demonstrated heightened sensitivity and rapid response to A. suturalis feeding, particularly. OPR3-1 and OPR3-2 maintained moderate levels (Figure 8G, pink and blue curves), while AOS experienced notable fluctuations (Figure 8G, green curve). This suggests that plants can adjust protein and gene levels to combat insect feeding. Understanding why plants modify proteins and metabolites post-insect attack is crucial for elucidating the plant–insect relationship, and our results provide the theoretical basis for it.

4. Discussion

Plants in natural environments encounter diverse herbivores, pathogens, and other biotic and abiotic stresses [58,59,60]. To survive these challenges, plants have evolved intricate mechanisms to respond to herbivore attacks, reshaping their transcriptomes, proteomes, and metabolomes upon detecting physical and chemical cues from herbivores [3,60,61]. These responses involve specific alterations in metabolism, gene expression, and plant growth and development patterns [1,62,63,64]. When attacked by herbivore feeding or egg laying, the injury initiates complex reactions that ultimately lead to plant defense [65], while insects can release effectors that disturb host plant defense responses for fitness, and those effectors are crucial components in insects and the host plant [66]. Plants can recognize herbivore-associated molecular patterns (HAMPs) and trigger various types of defense signal transduction, be it against chewing herbivores or piercing–sucking insects [67]. Effectors in oral secretion (OS) have been well reported to trigger or interfere with plant defense [68]. The first reported OS effector was fatty acid conjugates (FACs) of caterpillars [69]. Recently, researchers found an effector, HARP1, from H. armigera oral secretion. HARP1 may affect JA-responsive genes and make the plants more suitable for insect feeding [66]. In recent years, JA has attracted great attention and its functions in plant stress responses against pathogens, herbivorous as well as in development, have been well studied [70]. In our study, we identified five genes (PLA, AOS, AOC, OPR3-1, OPR3-2) enriched in the ALA signaling pathway contributing to increased JA levels (Figure 7). In addition, the ALA metabolism pathway demonstrated a pronounced response to A. suturalis plant feeding at both protein and metabolic levels. We also observed significantly increased expression of metabolites and proteases in the ALA pathway. This is consistent with a previous report observed in rice in response to stress induced by the rice stem borer [40]. Consequently, we hypothesized that the ALA metabolic pathway plays a crucial role in cotton response to mirid bugs, and the upregulated enzymes in the ALA metabolism pathway, including PLA, AOS, AOC and OPR3, positively contribute to JA biosynthesis (Figure 7). This finding aligns with previous studies indicating that AOS, AOC, and their families are vital in determining the resistance of rice to phloem-feeding herbivores [15,70]. However, our results also suggest that downregulated ADH in the ALA metabolism pathway negatively impacts JA biosynthesis (Figure 7).
Over the years, evidence has accumulated that plant primary and secondary metabolism profiles influence insect behavior [66,71]. Some plants may reduce nutrients that are required for insect survival and reproduction [57]. In addition, secondary metabolites may deter insects from feeding or egg laying and can also provide defensive functions and regulate defense signaling pathways against herbivores [65]. In our study, the integrated proteomic–metabolomic data analysis from two key metabolic pathways (fructose and mannose biosynthesis and ALA metabolism) indicated that cotton responds to feeding stress by decreasing nutrient content (e.g., fructose and mannose metabolic pathway) and increasing metabolite and protein content in related pathways (e.g., ALA metabolic pathway). These findings align with previous studies proposing two main hypotheses: active defense and nutrient stress [72,73]. The active defense hypothesis suggests that plants actively responded to insect feeding, resulting in increased metabolite content. In contrast, the nutrient stress hypothesis suggests that insect pests cause a reduction in plant nutrient content, affecting plant-eating insect growth and development [72,73]. Our results demonstrate that the cotton plant respond to A. suturalis feeding stress by increasing its secondary metabolites, including terpenes, alkaloids, phenylpropanoids, and flavonoids, which were most significantly upregulated post-feeding (Figure 1B,C). Correlational studies have identified small molecules with antifeedant or toxic effects on herbivorous insects. Notably, terpenoids, a metabolically diverse class of plant secondary metabolites, contribute to plant defenses [58,74]. Defensive properties of various compounds, such as furanocoumarins, cardenolides, tannins, saponins, glucosinolates, and cyanogenic glycosides, are well established [13,71].
In summary, the defense response of cotton to A. suturalis is a combined manifestation of both the nutrient stress hypothesis and the active defense hypothesis. This study offers a theoretical foundation and identifies candidate target genes for enhancing cotton plant defense resistance through molecular breeding. Therefore, these findings could pave the way for development of novel control strategies against this insect pest. Further studies on those regulators will enrich our knowledge about plant defense and help design novel strategies to control insect pests, most importantly helping us develop novel transgenic breeding techniques in cultivating cotton against insect pests.

5. Conclusions

In this study, we propose that A. suturalis triggers a complex cotton plant resistance process encompassing the upregulation of plant secondary metabolites (terpenes, alkaloid, phenylpropanoid, and flavonoids, Figure 9A), a reduction in nutrient content (metabolites and proteins of fructose and mannose biosynthesis pathway, Figure 9B), and an increase in the key pathway metabolites and proteins (such as those in the ALA metabolic pathway, Figure 9B). Our focus was on elucidating the signaling pathways for A. suturalis-induced plant resistance, particularly the ALA pathway, where most proteases are upregulated, leading to elevated JA levels. We propose that the ALA metabolic pathway significantly contributes to cotton’s defense response, with proteases like AOS, AOC, PLA, and OPR3 identified as key regulators in cotton resistance. Additionally, screening for anti-insect metabolites, proteins, or genes offers valuable insights for biological control of insect-associated pests.
Therefore, it is essential to acknowledge that insect-induced plant resistance is a multifaceted issue. This study specifically addresses cotton resistance under laboratory conditions. Further investigations are warranted, including assessing whether the gene expression increase persists under natural field conditions. Furthermore, long-term effects also merit further exploration.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/insects15040254/s1. Figure S1: PCA score of QC sample (A); PCA score plot (B); PLS-DA score plot (C); ESI+ represents positive ion mode, ESI- represent negative ion mode. Control group: group 1, treatment group: group 2; Table S1: Primer base sequence of cotton for RT-qPCR; Table S2: Overview of differential metabolites; Table S3: Numbers and classification of differentially expressed metabolites of cotton; Table S4: Category and FC analysis of differentially expressed secondary metabolites related to cotton defense; Table S5: Detailed information on 496 metabolites; Table S6: KEGG pathway analysis of different metabolites of cotton induced by A. suturalis; Table S7: Analysis of different metabolites of the two key metabolic pathways; Table S8: Total information of protein identification of each treatment; Table S9: A total of 1252 differentially expressed proteins (DEPs); Table S10: GO enrichment analysis of differentially expressed proteins; Table S11: KEGG enrichment of proteome of cotton plant; Table S12: Analysis of DEPs in different metabolism pathways.

Author Contributions

Conceptualization, J.C. and H.L.; methodology, H.L. and C.M.; software, C.M. and S.Z.; validation, H.L., S.Z. and J.J.; formal analysis, S.Z., X.G. and J.J.; investigation, H.L.; resources, H.H.; data curation, H.L. and X.G.; writing—original draft preparation, H.L.; writing—review and editing, H.L.; visualization, J.C. and J.L.; supervision, H.H.; project administration, J.C. and J.L.; funding acquisition, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Science and Technology Innovation 2030 (2023ZD04062), Key Project of Science and Technology XingMeng—action of durum wheat capacity increasing in Yellow River basin (NMKJXM202201), Hetao College Science and Technology Research Project (HYZY202004) and Henan Province Science and Technology Research project (242102110146).

Data Availability Statement

It may be found in the online version of this article. All data generated and analyzed during the current study are included in this article.

Acknowledgments

We thank the Science and Technology Innovation 2030, Key Project of Science and Technology XingMeng, and Hetao College Science and Technology Research Project for funding and Henan Province Science and Technology Research project and our colleagues Liming Lai and Hongxi Du for editorial assistance and language polishing. We apologize to all the researchers whose relevant work has not been mentioned here.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kessler, A.; Baldwin, I.T. Plant responses to insect herbivory: The emerging molecular analysis. Annu. Rev. Plant Biol. 2002, 53, 299–328. [Google Scholar] [CrossRef] [PubMed]
  2. Wu, J.; Baldwin, I.T. New insights into plant responses to the attack from insect herbivores. Annu. Rev. Genet. 2010, 44, 1–24. [Google Scholar] [CrossRef]
  3. Howe, G.A.; Jander, G. Plant immunity to insect herbivores. Annu. Rev. Plant Biol. 2008, 59, 41–66. [Google Scholar] [CrossRef] [PubMed]
  4. Xu, T.; Zhou, Q.; Chen, W.; Zhang, G.; He, G.F.; Gu, D.X.; Zhang, W.Q. Involvement of Jasmonate-signaling pathway in the herbivore-induced rice plant defense. Chin. Sci. Bull. 2003, 48, 1982–1987. [Google Scholar] [CrossRef]
  5. Cao, J.J.; Li, M.Y.; Chen, J.; Liu, P.; Li, Z. Effects of MeJA on Arabidopsis metabolome under endogenous JA deficiency. Sci. Rep. 2016, 6, 37674. [Google Scholar] [CrossRef]
  6. Townsend, B.J.; Poole, A.; Blake, C.J.; Llewellyn, D.J. Antisense suppression of a (+)-delta-cadinene synthase gene in cotton prevents the induction of this defense response gene during bacterial blight infection but not its constitutive expression. Plant Physiol. 2005, 138, 516–528. [Google Scholar] [CrossRef] [PubMed]
  7. Wasternack, C.; Hause, B. Jasmonates: Biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in Annals of Botany. Ann. Bot. 2013, 111, 1021–1058. [Google Scholar] [CrossRef] [PubMed]
  8. Wasternack, C.; Song, S. Jasmonates: Biosynthesis, metabolism, and signaling by proteins activating and repressing transcription. J. Exp. Bot. 2017, 68, 1303–1321. [Google Scholar] [CrossRef]
  9. Chini, A.; Monte, I.; Zamarreño, A.M.; Hamberg, M.; Lassueur, S.; Reymond, P.; Weiss, S.; Stintzi, A.; Schaller, A.; Porzel, A.; et al. An OPR3-independent pathway uses 4,5-didehydrojasmonate for jasmonate synthesis. Nat. Chem. Biol. 2018, 14, 171–178. [Google Scholar] [CrossRef]
  10. Li, C.; Xu, M.; Cai, X.; Han, Z.; Si, J.; Chen, D. Jasmonate signaling pathway modulates plant defense, growth, and their trade-offs. Int. J. Mol. Sci. 2022, 23, 3945. [Google Scholar] [CrossRef]
  11. Griffiths, G. Jasmonates: Biosynthesis, perception and signal transduction. Essays Biochem. 2020, 64, 501–512. [Google Scholar] [CrossRef] [PubMed]
  12. Zhuang, Y.; Wang, X.; Llorca, L.C.; Lu, J.; Lou, Y.; Li, R. Role of jasmonate signaling in rice resistance to the leaf folder Cnaphalocrocis medinalis. Plant Mol. Biol. 2022, 109, 627–637. [Google Scholar] [CrossRef] [PubMed]
  13. Zheng, S.C.; Liu, W.J.; Luo, J.Y.; Wang, L.S.; Zhu, X.Z.; Gao, X.K.; Hua, H.X.; Cui, J.J. Helicoverpa armigera herbivory negatively impacts Aphis gossypii populations via inducible metabolic changes. Pest Manag. Sci. 2022, 78, 2357–2369. [Google Scholar] [CrossRef] [PubMed]
  14. Jin, G.C.; Qi, J.F.; Zu, H.Y.; Liu, S.T.; Gershenzon, J.; Lou, Y.G.; Baldwin, I.T.; Li, R. Jasmonate-mediated gibberellin catabolism constrains growth during herbivore attack in rice. Plant Cell 2023, 35, 3828–3844. [Google Scholar] [CrossRef] [PubMed]
  15. Zeng, J.; Zhang, T.; Huangfu, J.; Li, R.; Lou, Y. Both allene oxide synthases genes are involved in the biosynthesis of herbivore-induced jasmonic acid and herbivore resistance in rice. Plants 2021, 10, 442. [Google Scholar] [CrossRef]
  16. Zhang, H.; Zhang, Y.X.; Xu, N.; Rui, C.; Fan, Y.P.; Wang, J.; Han, M.; Wang, Q.Q.; Sun, L.Q.; Chen, X.G.; et al. Genome-wide expression analysis of phospholipase A1 (PLA1) gene family suggests phospholipase A1-32 gene responding to abiotic stresses in cotton. Int. J. Biol. Macromol. 2021, 192, 1058–1074. [Google Scholar] [CrossRef] [PubMed]
  17. Sun, T.T.; Chen, Y.; Feng, A.Y.; Zou, W.H.; Wang, D.J.; Lin, P.X.; Chen, Y.L.; You, C.H.; Que, Y.X.; Su, Y.C. The allene oxide synthase gene family in sugarcane and its involvement in disease resistance. Ind. Crops Prod. 2023, 192, 116–136. [Google Scholar] [CrossRef]
  18. Available online: https://www.fao.org/newsroom/detail/world-cotton-day-celebrating-the-role-of-cotton-in-global-development-while-calling-for-developing-the-crop-more-sustainably/en (accessed on 12 January 2024).
  19. Lu, Y.H.; Wu, K.M.; Jiang, Y.Y.; Xia, B.; Li, P.; Feng, H.Q.; Wyckhuys, K.A.; Guo, Y. Mirid bug outbreaks in multiple crops correlated with wide-scale adoption of Bt cotton in China. Science 2010, 328, 1151–1154. [Google Scholar] [CrossRef] [PubMed]
  20. Wu, K.M.; Lu, Y.H.; Feng, H.Q.; Jiang, Y.Y.; Zhao, J.Z. Suppression of cotton bollworm in multiple crops in China in areas with Bt toxin-containing cotton. Science 2008, 321, 1676–1678. [Google Scholar] [CrossRef]
  21. Lu, Y.H.; Wu, K.M. Mirid bugs in China: Pest status and management strategies. Outlooks Pest Manag. 2011, 22, 248–252. [Google Scholar] [CrossRef]
  22. Luo, J.; Liu, X.; Liu, L.; Zhang, P.; Chen, L.Z.; Gao, Q.; Ma, W.H.; Chen, L.Z.; Lei, C. De novo analysis of the Adelphocoris suturalis Jakovlev metathoracic scent glands transcriptome and expression patterns of pheromone biosynthesis-related genes. Gene 2014, 551, 271–278. [Google Scholar] [CrossRef]
  23. Tong, Y.J.; Wu, K.M.; Lu, Y.H.; Gao, X.W. Pathogenicity of Beauveria spp. strains to three species of mirids, Apolygus lucorum, Adelphocoris suturalis and Adelphocoris lineolatus. Agric. Food Sci. 2010, 37, 172176. Available online: https://api.semanticscholar.org/CorpusID:90497852 (accessed on 2 March 2018).
  24. Luo, J.; Liang, S.j.; Li, J.Y.; Xu, Z.P.; Li, L.; Zhu, B.Q.; Li, Z.; Lei, C.L.; Lindsey, K.; Chen, L.Z.; et al. A transgenic strategy for controlling plant bugs (Adelphocoris suturalis) through expression of doublestranded RNA homologous to fatty acyl-coenzyme a reductase in cotton. New Phytol. 2017, 215, 1173–1185. [Google Scholar] [CrossRef]
  25. Ma, M.Q.; Guo, L.B.; Tu, C.J.; Wang, A.L.; Xu, L.T.; Luo, J. Comparative analysis of Adelphocoris suturalis Jakovlev (Hemiptera: Miridae) immune responses to fungal and bacterial pathogens. Front. Physiol. 2021, 18, 646721. [Google Scholar] [CrossRef]
  26. Xue, H.; Zhu, X.Z.; Wang, L.; Zhang, K.X.; Li, D.Y.; Ji, J.C.; Niu, L.; Wu, C.C.; Gao, X.K.; Luo, J.Y.; et al. Gut bacterial diversity in different life cycle stages of Adelphocoris suturalis (Hemiptera: Miridae). Front. Microbiol. 2021, 12, 670383. [Google Scholar] [CrossRef]
  27. Lu, Y.H.; Jiao, Z.; Li, G.; Wyckhuys, K.A.G.; Wu, K.M. Comparative overwintering host range of three Adelphocoris species (Hemiptera: Miridae) in northern China. Crop Protect. 2011, 30, 1455–1460. [Google Scholar] [CrossRef]
  28. Luo, J.; Li, Z.; Ma, C.; Zhang, Z.; Hull, J.J.; Lei, C.; Jin, S.; Chen, L. Knockdown of a metathoracic scent gland desaturase enhances the production of (E)-4-oxo-2-hexenal and suppresses female sexual attractiveness in the plant bug Adelphocoris suturalis. Insect Mol. Biol. 2017, 26, 642–653. [Google Scholar] [CrossRef]
  29. Xue, H.; Huang, X.X.; Chang, G.F.; Ma, W.H.; Hull, J.J.; Chen, L.Z. Reproductive capacity in Adelphocoris suturalis (Hemiptera: Miridae) is regulated by the insulin signaling pathway. Pestic. Biochem. Physiol. 2022, 187, 105195. [Google Scholar] [CrossRef]
  30. Qin, S.; Zhu, B.; Huang, X.; Hull, J.J.; Chen, L.; Luo, J. Functional Role of AsAP in the Reproduction of Adelphocoris suturalis (Hemiptera: Miridae). Insects 2022, 13, 755. [Google Scholar] [CrossRef]
  31. Zhu, Z.J.; Schultz, A.W.; Wang, J.; Johnson, C.H.; Yannone, S.M.; Patti, G.J.; Siuzdak, G. Liquid chromatography quadrupole time-of-flight mass spectrometry characterization of metabolites guided by the METLIN database. Nat. Protoc. 2013, 8, 451–460. [Google Scholar] [CrossRef]
  32. Xia, J.; Sinelnikov, I.V.; Han, B.; Wishart, D.S. MetaboAnalyst 3.0—Making metabolomics more meaningful. Nucleic Acids Res. 2015, 43, W251–W257. [Google Scholar] [CrossRef] [PubMed]
  33. Alseekh, S.; Aharoni, A.; Brotman, Y.; Contrepois, K.; Auria, J.D.; Ewald, J.; Ewald, J.C.; Fraser, P.D.; Giavalisco, P.; Hall, R.D.; et al. Mass spectrometry-based metabolomics: A guide for annotation, quantification and best reporting practices. Nat. Methods 2021, 18, 747–756. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, W.W.; Wang, G.; Baek, S.J.; Shen, R.F. Comparative study of three proteomic quantitative methods, DIGE, cICAT, and iTRAQ, using 2D gel- or LC-MALDI TOF/TOF. J. Proteome Res. 2006, 5, 651–658. [Google Scholar] [CrossRef] [PubMed]
  35. Li, Y.; Wang, X.; Hou, Y.; Zhou, X.Y.; Chen, Q.M.; Guo, C.; Xia, Q.Y.; Zhang, Y.; Zhao, P. Integrative proteomics and metabolomics analysis of insect larva brain: Novel insights into the molecular mechanism of insect wandering behavior. J. Proteome Res. 2016, 15, 193–204. [Google Scholar] [CrossRef] [PubMed]
  36. Huang, H.J.; Liu, C.W.; Huang, X.H.; Zhou, X.; Zhuo, J.C.; Zhang, C.X.; Bao, Y.Y. Screening and functional analyses of Nilaparvata lugens salivary proteome. J. Proteome Res. 2016, 15, 1883–1896. [Google Scholar] [CrossRef]
  37. Keeling, C.I.; Li, M.; Sandhu, H.K.; Henderson, H.; Yuen, M.M.; Bohlmann, J. Quantitative metabolome, proteome and transcriptome analysis of midgut and fat body tissues in the mountain pine beetle, Dendroctonus ponderosae Hopkins, and insights into pheromone biosynthesis. Insect Biochem. Mol. Biol. 2016, 70, 170–183. [Google Scholar] [CrossRef] [PubMed]
  38. Tu, X.; Wang, J.; Hao, K.; Whitman, D.W.; Fan, Y.; Cao, G.; Zhang, Z. Transcriptomic and proteomic analysis of pre-diapause and non-diapause eggs of migratory locust, Locusta migratoria L. (Orthoptera: Acridoidea). Sci. Rep. 2015, 5, 11402. [Google Scholar] [CrossRef]
  39. Liu, J.; Pang, C.Y.; Wei, H.L.; Song, M.Z.; Meng, Y.Y.; Ma, J.H.; Fan, S.L.; Yu, S.X. Data for proteomic profiling of Anthers from a photosensitive male sterile mutant and wild-type cotton (Gossypium hirsutum L.). Data Brief 2015, 4, 500–509. [Google Scholar] [CrossRef]
  40. Liu, Q.S.; Wang, X.Y.; Tzin, V.; Romeis, J.; Peng, Y.F.; Li, Y.H. Combined transcriptome and metabolome analyses to understand the dynamic responses of rice plants to attack by the rice stem borer Chilo suppressalis (Lepidoptera: Crambidae). BMC Plant Biol. 2016, 16, 259. [Google Scholar] [CrossRef]
  41. Zhou, B.; Xiao, J.F.; Tuli, L.; Ressom, H.W. LC-MS-based metabolomics. Mol. Biosyst. 2012, 2, 470–481. [Google Scholar] [CrossRef]
  42. Maiväli, Ü. Chapter 4—Data and Evidence. In Interpreting Biomedical Science; Academic Press: Cambridge, MA, USA, 2015; pp. 159–221. ISBN 9780124186897. [Google Scholar] [CrossRef]
  43. Zhou, C.X.; Zhou, D.H.; Elsheikha, H.M.; Zhao, Y.; Suo, X.; Zhu, X.Q. Metabolomic profiling of mice serum during toxoplasmosis progression using liquid chromatography-mass spectrometry. Sci. Rep. 2016, 6, 19557. [Google Scholar] [CrossRef]
  44. Zhao, S.Y.; Liu, Z.L.; Shu, Y.S.; Wang, M.L.; He, D.; Song, Z.Q.; Zeng, H.L.; Ning, Z.C.; Lu, C.; Lu, A.P.; et al. Chemotaxonomic classification applied to the identification of two closely-related citrus TCMs using UPLC-Q-TOF-MS-based metabolomics. Molecules 2017, 22, 1721. [Google Scholar] [CrossRef]
  45. Chen, T.T.; Zhang, L.; Shang, H.; Liu, S.; Peng, J.; Gong, W.; Shi, Y.; Zhang, S.; Li, J.; Gong, J.; et al. iTRAQ-Based quantitative proteomic analysis of cotton roots and leaves reveals pathways associated with salt stress. PLoS ONE 2016, 11, e0148487. [Google Scholar] [CrossRef]
  46. Saravanan, R.S.; Rose, J.K. A critical evaluation of sample extraction techniques for enhanced proteomic analysis of recalcitrant plant tissues. Proteomics 2004, 4, 2522–2532. [Google Scholar] [CrossRef]
  47. Petricka, J.J.; Schauer, M.A.; Megraw, M.; Breakfield, N.W.; Thompson, J.W.; Georgiev, S.; Soderblom, E.J.; Ohler, U.; Moseley, M.A.; Grossniklaus, U.; et al. The protein expression landscape of the Arabidopsis root. Proc. Natl. Acad. Sci. USA 2012, 109, 6811–6818. [Google Scholar] [CrossRef]
  48. Kruger, N.J. The Bradford method for protein quantitation. Methods Mol. Biol. 1994, 32, 9–15. [Google Scholar] [CrossRef]
  49. Liu, J.F.; Xing, J.Y.; Fang, J.; Ai, P.F.; Cheng, Y.Q. New insight into ovary abortion during ovary development of hazelnut through a combined proteomic and transcriptomic analysis. Sci. Hortic. 2018, 234, 36–48. [Google Scholar] [CrossRef]
  50. Chu, P.; Yan, G.X.; Yang, Q.; Zhai, L.N.; Zhang, C.; Zhang, F.Q.; Guan, R.Z. iTRAQ-based quantitative proteomics analysis of Brassica napus leaves reveals pathways associated with chlorophyll deficiency. J. Proteomics 2015, 113, 244–259. [Google Scholar] [CrossRef]
  51. Zhu, M.M.; Dai, S.J.; McClung, S.; Yan, X.F.; Chen, S.X. Functional differentiation of Brassica napus guard cells and mesophyll cells revealed by comparative proteomics. Mol. Cell. Proteomics 2009, 8, 752–766. [Google Scholar] [CrossRef]
  52. Chen, Y.Z.; Pang, Q.Y.; He, Y.; Zhu, N.; Branstrom, I.; Yan, X.F.; Chen, S. Proteomics and metabolomics of Arabidopsis responses to perturbation of glucosinolate biosynthesis. Mol. Plant 2012, 5, 1138–1150. [Google Scholar] [CrossRef]
  53. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar] [CrossRef]
  54. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using Real-Time Quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  55. Wang, M.; Wang, Q.L.; Zhang, B.H. Evaluation and selection of reliable reference genes for gene expression under abiotic stress in cotton (Gossypium hirsutum L.). Gene 2013, 530, 44–50. [Google Scholar] [CrossRef]
  56. Jonsson, P.; Bruce, S.J.; Moritz, T.; Trygg, J.; Sjöström, M.; Plumb, R.; Granger, J.; Maibaum, E.; Nicholson, J.K.; Holmes, E.; et al. Extraction, interpretation and validation of information for comparing samples in metabolic LC/MS data sets. Analyst 2005, 130, 701–707. [Google Scholar] [CrossRef]
  57. Chiozza, M.V.; O’Neal, M.E.; MacIntosh, G.C. Constitutive and induced differential accumulation of amino acid in leaves of susceptible and resistant soybean plants in response to the soybean aphid (Hemiptera: Aphididae). Environ. Entomol. 2010, 39, 856–864. [Google Scholar] [CrossRef]
  58. Gols, R. Direct and indirect chemical defences against insects in a multitrophic framework. Plant Cell Environ. 2014, 7, 1741–1752. [Google Scholar] [CrossRef]
  59. Erb, M.; Reymond, P. Molecular interactions between plants and insect herbivores. Annu. Rev. Plant Biol. 2019, 29, 527–557. [Google Scholar] [CrossRef]
  60. Fürstenberg-Hägg, J.; Zagrobelny, M.; Bak, S. Plant defense against insect herbivores. Int. J. Mol. Sci. 2013, 14, 10242–10297. [Google Scholar] [CrossRef]
  61. Beran, F.; Petschenka, G. Sequestration of plant defense compounds by insects: From mechanisms to insect-plant coevolution. Annu. Rev. Entomol. 2022, 67, 163–180. [Google Scholar] [CrossRef] [PubMed]
  62. Truong, D.; Heuskin, S.; Delaplace, P.; Francis, F.; Lognay, G. VOC emissions and protein expression mediated by the interactions between herbivorous insects and Arabidopsis plant. Biotechnol. Agron. Soc. Environ. 2014, 18, 455–464. [Google Scholar]
  63. Bede, J.C.; Musser, R.O.; Felton, G.W.; Korth, K.L. Caterpillar herbivory and salivary enzymes decrease transcript levels of Medicago truncatula genes encoding early enzymes in terpenoid biosynthesis. Plant Mol. Biol. 2006, 60, 519–531. [Google Scholar] [CrossRef] [PubMed]
  64. Bonaventure, G. Plants recognize. herbivorous insects by complex signalling networks. Annu. Plant Rev. 2014, 47, 1–35. [Google Scholar] [CrossRef]
  65. Divekar, P.A.; Narayana, S.; Divekar, B.A.; Kumar, R.; Gadratagi, B.G.; Ray, A.; Singh, A.K.; Rani, V.; Singh, V.; Singh, A.K.; et al. Plant secondary metabolites as defense tools against herbivores for sustainable crop protection. Int. J. Mol. Sci. 2022, 23, 2690. [Google Scholar] [CrossRef] [PubMed]
  66. Hogenhout, S.A.; Bos, J.I. Effector proteins that modulate plant--insect interactions. Curr. Opin. Plant. Biol. 2011, 14, 422–428. [Google Scholar] [CrossRef]
  67. Stam, J.M.; Kroes, A.; Li, Y.; Gols, R.; van Loon, J.J.; Poelman, E.H.; Dicke, M. Plant interactions with multiple insect herbivores: From community to genes. Annu. Rev. Plant. Biol. 2014, 65, 689–713. [Google Scholar] [CrossRef] [PubMed]
  68. Alborn, H.T.; Turlings, T.C.J.; Jones, T.H.; Stenhagen, G.; Loughrin, J.H.; Tumlinson, J.H. An Elicitor of plant volatiles from beet armyworm oral secretion. Science 1997, 276, 945–949. [Google Scholar] [CrossRef]
  69. Chen, C.Y.; Liu, Y.Q.; Song, W.M.; Chen, D.Y.; Chen, F.Y.; Chen, X.Y.; Chen, Z.W.; Ge, S.X.; Wang, C.Z.; Zhan, S.; et al. An effector from cotton bollworm oral secretion impairs host plant defense signaling. Proc. Natl. Acad. Sci. USA 2019, 116, 14331–14338. [Google Scholar] [CrossRef]
  70. Wang, Y.; Mostafa, S.; Zeng, W.; Jin, B. Function and mechanism of jasmonic acid in plant responses to abiotic and biotic stresses. Int. J. Mol. Sci. 2021, 22, 8568. [Google Scholar] [CrossRef]
  71. Erb, M.; Kliebenstein, D.J. Plant secondary metabolites as defenses, regulators, and primary metabolites: The blurred functional trichotomy. Plant Physiol. 2020, 184, 39–52. [Google Scholar] [CrossRef]
  72. Haukioja, E.; Neuvonen, S. Induced long-term resistance of birth foliage against defoliators: Defensive or incidental? Ecology 1985, 66, 1303–1308. [Google Scholar] [CrossRef]
  73. Hunter, M.D.; Schultz, J.C. Fetilization mitigates chemical induction and herbivore responses within damaged oak trees. Ecology 1995, 76, 1226–1232. [Google Scholar] [CrossRef]
  74. Ninkuu, V.; Zhang, L.; Yan, J.P.; Fu, Z.C.; Yang, T.F.; Zeng, H.M. Biochemistry of terpenes andrecent advances in plant protection. Int. J. Mol. Sci. 2021, 22, 5710. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Differentially expressed metabolites in cotton. (A) Category of differentially expressed metabolites. (B) FC analysis of differentially expressed secondary metabolites related to cotton defense. (C) VIP analysis of differentially expressed secondary metabolites related to cotton defense.
Figure 1. Differentially expressed metabolites in cotton. (A) Category of differentially expressed metabolites. (B) FC analysis of differentially expressed secondary metabolites related to cotton defense. (C) VIP analysis of differentially expressed secondary metabolites related to cotton defense.
Insects 15 00254 g001
Figure 2. FC and VIP analysis of metabolism pathway. (A) FC and (B) VIP of fructose and mannose metabolism pathway. (C) FC and (D) VIP of ALA metabolism pathway.
Figure 2. FC and VIP analysis of metabolism pathway. (A) FC and (B) VIP of fructose and mannose metabolism pathway. (C) FC and (D) VIP of ALA metabolism pathway.
Insects 15 00254 g002
Figure 3. Protein identification and GO enrichment analysis of DEPs. (A) Up- and downregulated protein in cotton plants after mirid bug feeding. (B) GO annotation and functional classification of differentially expressed proteins in cotton plants, categorized as follows: biological process (A: metabolic process, B: single-organism process, C: cellular process metabolic process, D: response to stimulus, E: localization, F: reproductive process, G: signaling, H: cellular process, I: cellular component organization or biogenesis, J: others); cellular component (K: organelle part, L: organelle, M: cell part, N: membrane part, O: macromolecular complex, P: extracellular region, Q: membrane); and molecular function (R: catalytic activity, S: binding, T: transporter activity, U: structural molecule activity, V: antioxidant activity, W: enzyme regulator activity, X: others).
Figure 3. Protein identification and GO enrichment analysis of DEPs. (A) Up- and downregulated protein in cotton plants after mirid bug feeding. (B) GO annotation and functional classification of differentially expressed proteins in cotton plants, categorized as follows: biological process (A: metabolic process, B: single-organism process, C: cellular process metabolic process, D: response to stimulus, E: localization, F: reproductive process, G: signaling, H: cellular process, I: cellular component organization or biogenesis, J: others); cellular component (K: organelle part, L: organelle, M: cell part, N: membrane part, O: macromolecular complex, P: extracellular region, Q: membrane); and molecular function (R: catalytic activity, S: binding, T: transporter activity, U: structural molecule activity, V: antioxidant activity, W: enzyme regulator activity, X: others).
Insects 15 00254 g003
Figure 4. KEGG pathway enrichment. (A) Top 20 enriched pathways. Arrows represent metabolic pathways of interest. (B) Different metabolite changes of the defense-related metabolic pathways.
Figure 4. KEGG pathway enrichment. (A) Top 20 enriched pathways. Arrows represent metabolic pathways of interest. (B) Different metabolite changes of the defense-related metabolic pathways.
Insects 15 00254 g004
Figure 5. Common metabolic pathways of differential proteins and differential metabolites. Pink represents the number of differentially expressed metabolites (DEMs); blue represents the number of differentially expressed proteins (DEPs).
Figure 5. Common metabolic pathways of differential proteins and differential metabolites. Pink represents the number of differentially expressed metabolites (DEMs); blue represents the number of differentially expressed proteins (DEPs).
Insects 15 00254 g005
Figure 6. Integrated metabolomic and proteomic analysis of the fructose and mannose metabolism pathway. Red represents upregulated protease or metabolites, blue represents downregulated protease or metabolites, and black represents no change happening. The upregulated proteases were xylose isomerase (CotAD_24396, CotAD_63919, CotAD_02520), idonate-5-dehydrogenase (Cotton_D_gene_10030706), epidermis-specific secreted glycoprotein EP1 (CotAD_00556), phosphofructokinase (Cotton_D_gene_10039158), and fructose-bisphosphate aldolase3 (CotAD_58316). The downregulated proteases comprised hexokinase (CotAD_39434), fructokinase (CotAD_45862, CotAD_26511, CotAD_36008, CotAD_69191, Cotton_D_gene_10035332), endo-1,4-beta-mannanase (CotAD_57668), fructokinase-1-6-bisphosphatase(CotAD_39653), fructose-bisphosphate aldolase (CotAD_24238, CotAD_30698, CotAD_30700, CotAD_58944, CotAD_27060, CotAD_24236, CotAD_58035, CotAD_52093) and triosephosphate isomerase (CotAD_15944, CotAD_23341). The numbers in brackets represent FCs.
Figure 6. Integrated metabolomic and proteomic analysis of the fructose and mannose metabolism pathway. Red represents upregulated protease or metabolites, blue represents downregulated protease or metabolites, and black represents no change happening. The upregulated proteases were xylose isomerase (CotAD_24396, CotAD_63919, CotAD_02520), idonate-5-dehydrogenase (Cotton_D_gene_10030706), epidermis-specific secreted glycoprotein EP1 (CotAD_00556), phosphofructokinase (Cotton_D_gene_10039158), and fructose-bisphosphate aldolase3 (CotAD_58316). The downregulated proteases comprised hexokinase (CotAD_39434), fructokinase (CotAD_45862, CotAD_26511, CotAD_36008, CotAD_69191, Cotton_D_gene_10035332), endo-1,4-beta-mannanase (CotAD_57668), fructokinase-1-6-bisphosphatase(CotAD_39653), fructose-bisphosphate aldolase (CotAD_24238, CotAD_30698, CotAD_30700, CotAD_58944, CotAD_27060, CotAD_24236, CotAD_58035, CotAD_52093) and triosephosphate isomerase (CotAD_15944, CotAD_23341). The numbers in brackets represent FCs.
Insects 15 00254 g006
Figure 7. Integrated metabolomic and proteomic analysis of the ALA metabolism pathway. Red represents upregulated protease or metabolites, blue represents downregulated proteases or metabolites, and black represents no change happening. The upregulated proteases were: PLA (CotAD_52791), AOS (CotAD_35840, Cotton_D_gene_10023640, CotAD_58616), AOC (Cotton_D_gene_10007846, Cotton_D_gene_10007844), OPR3 (CotAD_59461, Cotton_D_gene_10037325), 4-coumarate-CoAligase (CotAD_59374), ACX (CotAD_52391, CotAD_12782, Cotton_D_gene_10040584), MFP2 (CotAD_18083), KAT2 (CotAD_51900, CotAD_10744, CotAD_66777), salicylate O-methyltransferase (CotAD_27039). The downregulated proteases comprised PLA (CotAD_64088) and ADH (CotAD_64304). The numbers in brackets represent FCs.
Figure 7. Integrated metabolomic and proteomic analysis of the ALA metabolism pathway. Red represents upregulated protease or metabolites, blue represents downregulated proteases or metabolites, and black represents no change happening. The upregulated proteases were: PLA (CotAD_52791), AOS (CotAD_35840, Cotton_D_gene_10023640, CotAD_58616), AOC (Cotton_D_gene_10007846, Cotton_D_gene_10007844), OPR3 (CotAD_59461, Cotton_D_gene_10037325), 4-coumarate-CoAligase (CotAD_59374), ACX (CotAD_52391, CotAD_12782, Cotton_D_gene_10040584), MFP2 (CotAD_18083), KAT2 (CotAD_51900, CotAD_10744, CotAD_66777), salicylate O-methyltransferase (CotAD_27039). The downregulated proteases comprised PLA (CotAD_64088) and ADH (CotAD_64304). The numbers in brackets represent FCs.
Insects 15 00254 g007
Figure 8. qRT-PCR analysis of the key genes in ALA metabolism pathway in cotton leaves. (A) Key regulator gene PLA, (B) key regulator gene AOS, (C) key regulator gene AOC, (D) key regulator gene OPR3-1, (E) key regulator gene OPR3-2, (F) downregulated gene ADH, (G) dynamic variation in the six ALA metabolism pathway-related genes in cotton plant. “*” represents p < 0.05 and “**” represents p < 0.01.
Figure 8. qRT-PCR analysis of the key genes in ALA metabolism pathway in cotton leaves. (A) Key regulator gene PLA, (B) key regulator gene AOS, (C) key regulator gene AOC, (D) key regulator gene OPR3-1, (E) key regulator gene OPR3-2, (F) downregulated gene ADH, (G) dynamic variation in the six ALA metabolism pathway-related genes in cotton plant. “*” represents p < 0.05 and “**” represents p < 0.01.
Insects 15 00254 g008
Figure 9. Hypothetic model of A. suturalis-induced cotton plant resistance. (A) Changes in secondary metabolites; (B) two key metabolic pathways.
Figure 9. Hypothetic model of A. suturalis-induced cotton plant resistance. (A) Changes in secondary metabolites; (B) two key metabolic pathways.
Insects 15 00254 g009
Table 1. Top 15 KEGG pathway analysis of differential metabolites of cotton induced by A. suturalis.
Table 1. Top 15 KEGG pathway analysis of differential metabolites of cotton induced by A. suturalis.
No.Pathway_IDKEGG Pathway NameTotalDownUp
1map00909Sesquiterpenoid and triterpenoid biosynthesis40040
2map00590Arachidonic acid metabolism28721
3map00380Tryptophan metabolism23518
4map00904Diterpenoid biosynthesis20614
5map00950Isoquinoline alkaloid biosynthesis19514
6map00960 Tropane, piperidine and pyridine alkaloid biosynthesis19217
7map00940Phenylpropanoid biosynthesis18216
8map00350Tyrosine metabolism17215
9map00592alpha-Linolenic acid metabolism16214
10map00902Monoterpenoid biosynthesis15510
11map00051Fructose and mannose metabolism1486
12map00360 Phenylalanine metabolism14212
13map00903Limonene and pinene degradation1439
14map00130 Ubiquinone and other terpenoid-quinone biosynthesis13310
15map00330Arginine and proline metabolism1129
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lu, H.; Zheng, S.; Ma, C.; Gao, X.; Ji, J.; Luo, J.; Hua, H.; Cui, J. Integrated Omics Analysis Reveals Key Pathways in Cotton Defense against Mirid Bug (Adelphocoris suturalis Jakovlev) Feeding. Insects 2024, 15, 254. https://doi.org/10.3390/insects15040254

AMA Style

Lu H, Zheng S, Ma C, Gao X, Ji J, Luo J, Hua H, Cui J. Integrated Omics Analysis Reveals Key Pathways in Cotton Defense against Mirid Bug (Adelphocoris suturalis Jakovlev) Feeding. Insects. 2024; 15(4):254. https://doi.org/10.3390/insects15040254

Chicago/Turabian Style

Lu, Hui, Shuaichao Zheng, Chao Ma, Xueke Gao, Jichao Ji, Junyu Luo, Hongxia Hua, and Jinjie Cui. 2024. "Integrated Omics Analysis Reveals Key Pathways in Cotton Defense against Mirid Bug (Adelphocoris suturalis Jakovlev) Feeding" Insects 15, no. 4: 254. https://doi.org/10.3390/insects15040254

APA Style

Lu, H., Zheng, S., Ma, C., Gao, X., Ji, J., Luo, J., Hua, H., & Cui, J. (2024). Integrated Omics Analysis Reveals Key Pathways in Cotton Defense against Mirid Bug (Adelphocoris suturalis Jakovlev) Feeding. Insects, 15(4), 254. https://doi.org/10.3390/insects15040254

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop