Comparative transcriptome and histological analyses of wheat in response to phytotoxic aphid Schizaphis graminum and non-phytotoxic aphid Sitobion avenae feeding
BMC Plant Biology volume 19, Article number: 547 (2019)
Infestation of the phytotoxic aphid Schizaphis graminum can rapidly induce leaf chlorosis in susceptible plants, but this effect is not observed with the nonphytotoxic aphid Sitobion avenae. However, few studies have attempted to identify the different defence responses induced in wheat by S. graminum and S. avenae feeding and the mechanisms underlying the activation of chlorosis by S. graminum feeding.
S. graminum feeding significantly reduced the chlorophyll content of wheat leaves, and these effects were not observed with S. avenae. A transcriptomic analysis showed that the expression levels of genes involved in the salicylic acid, jasmonic acid and ethylene signalling defence pathways were significantly upregulated by both S. avenae and S. graminum feeding; however, more plant defence genes were activated by S. graminum feeding than S. avenae feeding. The transcript levels of genes encoding cell wall-modifying proteins were significantly increased after S. graminum feeding, but only a few of these genes were induced by S. avenae. Furthermore, various reactive oxygen species-scavenging genes, such as 66 peroxidase (POD) and 8 ascorbate peroxidase (APx) genes, were significantly upregulated after S. graminum feeding, whereas only 15 POD and one APx genes were induced by S. avenae feeding. The activity of four antioxidant enzymes was also significantly upregulated by S. graminum feeding. Cytological examination showed that S. graminum feeding induced substantial hydrogen peroxide (H2O2) accumulation in wheat leaves. The chlorosis symptoms and the loss of chlorophyll observed in wheat leaves after S. graminum feeding were reduced and inhibited by the scavenging of H2O2 by dimethylthiourea, which indicated that H2O2 plays important role in the induction of chlorosis by S. graminum feeding.
S. graminum and S. avenae feeding induces the JA, SA and ET signalling pathways, but S. graminum activated stronger plant defence responses than S. avenae. S. graminum feeding triggers strong ROS-scavenging activity and massive H2O2 production in wheat leaves, and the accumulation of H2O2 induced by S. graminum feeding is involved in the activation of chlorosis in wheat leaves. These results enhance our understanding of mechanisms underlying aphid-wheat interactions and provide clues for the development of aphid-resistant wheat varieties.
Through interactions with insects over a hundred million years, plants have evolved complex and accurate defence mechanisms against herbivores. In response to herbivory, plants can perceive damage-associated molecular patterns (DAMPs) or herbivory-associated molecular patterns (HAMPs) in insect oral secretions and subsequently induce direct and indirect plant defence responses [1, 2]. Direct defences involve the production and accumulation of plant defensive chemicals, such as plant secondary metabolites (PSMs), proteinase inhibitors, polyphenol oxidases and other defensive proteins, which are induced by herbivory and reduce herbivore performance [3, 4]. Indirect defences include the synthesis and release of complex blends of volatiles that attract parasitoids and predators of the herbivores [5, 6].
Insects from different feeding guilds tend to elicit distinct plant defensive strategies of plants [7, 8]. For example, leaf-eating beetles (Coleoptera) or caterpillars (Lepidoptera) cause extensive tissue damage during herbivory, which usually activates the jasmonic acid (JA)-mediated defence pathway in plants [9, 10]. Different from leaf chewing insects, hemipterans, such as aphids and whiteflies, have highly modified piecing-sucking mouthparts (stylets) that can penetrate the extracellular pathway and feed on the nutrients from phloem sap provided by sieve elements (SEs). Although the stylets puncture through most parts of plant cells during the probing track, these insects cause less damage to cells than leaf-chewing herbivores . Many studies have demonstrated that hemipteran feeding results in the induction of the salicylic acid (SA)-dependent defence pathway in plants [12, 13].
The grain aphid Sitobion avenae and the greenbug Schizaphis graminum (Hemiptera: Aphididae) are considered two important pests of wheat and other cereals worldwide, as they suck phloem sap and serve as vectors of the barley yellow dwarf virus (BYDV), resulting in significant yield losses [14, 15]. Aphids are classified as nonphytotoxic or phytotoxic according to the symptoms and damage caused by their feeding . Similar to most aphid species, S. avenae are nonphytotoxic, and no typical symptoms of plant damage are rapidly induced by their feeding processes. However, infestation by S. graminum, a phytotoxic species, can immediately induce obvious leaf chlorosis in susceptible plants, resulting in the deterioration of plant quality and even in plant death [17, 18]. Several studies have characterized the plant defence responses induced by infestation of S. avenae and S. graminum. For example, Zhao et al. found that S. avenae feeding induces the expression of several genes involved in both the SA- and JA-mediated defence pathways in wheat , and Zhu-Salzman et al. demonstrated that S. graminum feeding strongly induces the expression of genes involved in the SA-dependent pathway in sorghum (Sorghum bicolor) but weakly increases the expression levels of JA-related defence genes, such as lipoxidase (LOX) and proteinase inhibitors (PIs) . Using a cDNA microarray, several transcripts with significantly different expression patterns were identified in sorghum seedlings after S. graminum feeding, and these herbivore-responsive genes were mainly associated with photosynthesis, the biosynthesis of defence molecules, cell wall fortification, oxidative bursts and stress .
Previous studies have also shown that S. graminum feeding significantly increases the concentrations of amino acids, particularly essential amino acids, in the phloem of wheat and barley, and the enhancement of the nutritional quality of plants induced by S. graminum feeding improves the aphid performance . It has been proposed that leaf-senescence-like changes triggered by S. graminum feeding are associated with the nutritional enhancement of host plants [18, 23]. Feeding damage caused by phytotoxic aphids is usually observed in susceptible hosts; therefore, it has been hypothesized that the nutritional enhancement of plants derived from senescence-like feeding damage is a strategy used by phytotoxic aphids to potentially counteract the negative effects of induced plant resistance, which would eventually improve aphid fitness on the host .
However, few studies have attempted to identify the molecular basis of the necrosis symptoms in wheat triggered by the phytotoxic aphid S. graminum and the different defence responses induced by S. graminum and S. avenae feeding. In this study, we integrated gene expression profiling through high-throughput RNA sequencing (RNA-Seq) with cytological examination to reveal the responses of wheat leaves to S. graminum and S. avenae feeding, to compare the differences in the metabolic pathways affected by the two cereal aphids and to uncover the mechanism underlying the induction of damage symptoms by S. graminum.
Damage symptoms and changes in the chlorophyll content of wheat leaves after S. avenae and S. graminum feeding
As shown in Fig. 1a-c, no obvious damage symptoms were detected in leaves 48 h after S. avenae infestation compared with leaves without aphid infestation, whereas S. graminum feeding caused severe chlorosis in wheat leaves. Delayed fluorescence, which was used as a direct indicator of the chlorophyll content, was also measured in the wheat leaves 48 h after aphid infestation. As demonstrated in Fig. 1d-f, the untreated leaves and the S. avenae-infested leaves exhibited strong signals of delayed fluorescence, whereas low signals of delayed fluorescence were detected in the wheat leaves infested with S. graminum, which suggested the occurrence of chlorophyll degradation.
The chlorophyll content of the leaves after aphid feeding was further investigated. As shown in Fig. 2, no significant differences in the chlorophyll content were found between the S. avenae-infested and control plants. However, the chlorophyll content was significantly decreased to 0.46 ± 0.068 (F2,6 = 10.494, P = 0.011) after S. graminum feeding.
Transcriptome data from aphid-infested wheat leaves
The transcriptomes of wheat leaves infested with the two cereal aphid species were compared in this study. A total of 62.98 Gb of clean data were obtained from the nine leaf samples, and each of these samples contained ≥6.91 Gb with Q30 quality scores of ≥94.82% (Additional file 1: Table S1). Subsequently, 83.5 to 94.3% of the clean reads from each sample were aligned onto the wheat reference genome and matched to either unique or multiple genomic locations (Table 1).
Identification and functional annotation of DEGs
The gene expression levels of each replicate were assessed through principal component analysis (PCA). Replicates from the same group were clustered closely together, which suggested that the repeatability of each treatment group was satisfactory, and the samples from the S. avenae- and S. graminum-infested groups (Sa48h and Sg48h, respectively) clustered far from the control samples, which indicated that aphid feeding induced significant changes in gene expression (Fig. 3). A total of 12,8195 transcripts were detected across all the samples (Additional file 2: Data S1). Gene expression levels with an adjusted P value < 0.00001 and |Log2Fold Change| ≥ 1 were selected as DEGs for further analysis. Forty-eight hours of S. avenae feeding significantly upregulated 1718 genes and significantly downregulated 172 genes in wheat leaves (Fig. 4a). In addition, 7893 and 5098 genes were significantly upregulated and downregulated, respectively, after 48 h of S. graminum feeding (Fig. 4b).
To investigate the differences in plant responses to infestation by S. graminum and S. avenae, the DEGs in wheat leaves induced by these aphids were also compared in our study. The results showed that the expression levels of 857 genes in wheat leaves were significantly upregulated by both S. graminum and S. avenae feeding. Additionally, 11,046 and 1914 transcripts were specifically and significantly upregulated after S. graminum and S. avenae feeding for 48 h, respectively (Fig. 5a). In contrast, a total of 128 transcripts were significantly downregulated after both S. graminum and S. avenae infestation, and 7036 and 861 genes were only significantly downregulated after S. graminum and S. avenae feeding, respectively (Fig. 5b). This finding suggested that the global response of wheat to S. graminum feeding is distinct from that of wheat to S. avenae feeding.
All the DEGs were subjected to GO term and KEGG pathway analyses to identify the major DEG-associated metabolic pathways. The top 30 enriched GO terms and 20 most enriched KEGG pathways are shown in Fig. 6 and Figure S1 (Additional file 3). As shown in Fig. 6, within the biological process category, the DEGs induced by S. graminum were mainly enriched in the electron transport, small molecule metabolic process and carbohydrate metabolic process terms, and the DEGs induced by S. avenae were mainly enriched in protein phosphorylation, protein modification process and phosphorus metabolic process. Within the molecular function category, the largest proportion of DEGs induced by S. graminum was enriched in the catalytic activity and oxidoreductase activity terms, and the majority of the DEGs activated by S. avenae were enriched in the catalytic activity, protein kinase activity and phosphotransferase activity.
Transcripts related to photosynthesis, sucrose and starch metabolism and nitrogen metabolism
S. graminum feeding negatively affected the photosynthesis process of wheat, and many genes associated with light-harvesting and photosystem-associated genes, such as chlorophyll a-b binding proteins, ferrochelatase, and photosystem I and II proteins, were significantly downregulated (Table 2). The expression levels of ribulose bisphosphate carboxylase oxygenase (RuBisCO) and carbonic anhydrase genes with roles in the Calvin cycle were also significantly reduced after S. graminum feeding. However, few genes involved in photosynthesis were significantly regulated in the S. avenae-infested plants. The transcriptional profiles of some genes involved in sucrose and starch metabolism were also investigated. The sucrose synthase 3 gene, one trehalose-6-phosphate synthase gene and six beta-glucosidase genes were significantly downregulated in S. graminum- and S. avenae-infested leaves, and the transcript levels of sucrose-phosphatase genes were significantly upregulated in wheat leaves infested with S. graminum and S. avenae. S. graminum feeding but not S. avenae feeding also significantly affected nitrogen metabolism. The transcript levels of nitrate reductase in the leaves were strongly downregulated by S. graminum, and glutamate dehydrogenase was significantly upregulated in wheat leaves infested with S. graminum. However, few genes involved in nitrogen metabolism were modulated by S. avenae feeding.
Effects of S. graminum and S. avenae feeding on transcripts related to the SA, JA, and ET signalling pathways involved in plant defence
To characterize how plant defence responses are modulated in response to S. graminum and S. avenae feeding, genes known to be involved in the SA, JA, and ET-defence pathways were examined . The transcriptome data in Table 3 showed that 21 PAL genes involved in SA biosynthesis were significantly upregulated in response to S. graminum and S. avenae feeding. Furthermore, 13 PR genes responding to SA were significantly upregulated by S. graminum and S. avenae feeding. Additionally, under S. graminum and S. avenae feeding, one AOC (3.7-fold), three AOS (2.7- to 7.9-fold) and five LOX (2.8- to 8.8-fold) genes involved in JA biosynthesis were significantly upregulated, and the expression levels of three PI genes, which are JA-responsive defence genes, significantly increased. Two ACS (5.1 to 7.8-fold) and three ACO (4.0 to 5.0-fold) genes, which are involved in the ET signalling pathway, and 11 genes encoding ethylene-responsive transcription factors involved in ET biosynthesis were significantly upregulated in response to S. graminum feeding, but only one ACO (4.4-fold) gene and two genes encoding ethylene-responsive transcription factors (1.8-fold and induced) were upregulated after S. avenae feeding.
Although the transcript levels of some defence genes were significantly upregulated in response to both S. graminum and S. avenae feeding, a higher number of DEGs involved in SA-, JA-, and ET-mediated defence pathways were induced by S. graminum feeding than by S. avenae feeding (Additional file 4: Data S2). For example, 37 PR genes (downstream of SA) and 10 PI genes (downstream of JA) were significantly upregulated in response to S. graminum feeding, but only 17 PR genes and four PI genes were significantly upregulated in response to S. avenae feeding. Additionally, S. graminum feeding induced greater fold changes in these two types of genes than S. avenae feeding, which indicated that the former triggered a stronger defence response than the latter.
Effects of S. graminum feeding on transcripts associated with plant cell wall modification proteins (PCMDPs) in wheat leaves
S. graminum feeding induced the expression of many genes encoding enzymatic or non-enzymatic proteins related to plant cell wall dynamics (Table 4). For example, the transcript levels of callose synthases were significantly upregulated after S. graminum and S. avenae feeding. The transcript levels of 19 PGs (1.92 to 6.45-fold; induced) and four PEMs (3.37 to 3.97-fold; induced) were significantly increased after 48 h of S. graminum feeding, but no PGs or PEMs were significantly induced in response to S. avenae feeding. Similarly, the transcript levels of six genes encoding beta-expansin, which is a non-enzymatic protein that plays important roles in cell wall loosening, were significantly upregulated (5.30 to 8.70-fold; induced) in wheat leaves infested with S. graminum, but no expansin genes were significantly regulated after 48 h of S. avenae feeding.
Effects of S. graminum and S. avenae feeding on the transcript levels and activities of antioxidant enzymes involved in ROS scavenging in wheat leaves
In plants, herbivore attacks usually trigger oxidative responses . Plants possess a battery of ROS scavengers, such as POD, SOD, and CAT enzymes, and these enzymes can protect cells from oxidative damage . As shown in Fig. 7, 74 PODs were significantly up- or down- regulated in response to S. graminum feeding, and 66 of these PODs were significantly upregulated. However, only 15 PODs were significantly upregulated by S. avenae feeding. Similarly, the expression levels of 12 APx genes were significantly modulated by S. graminum feeding, but the expression levels of only two APx genes were significantly affected by S. avenae feeding. Additionally, five CAT, eight SOD and seven glutathione peroxidase (GPx) genes were significantly regulated by S. graminum feeding, but not by S. avenae feeding (Additional file 5: Data S3). The increased number of ROS scavengers induced by S. graminum feeding suggested that S. graminum feeding induces stronger oxidative stress in wheat leaves than S. avenae.
Compared with the control, the activity of POD was significantly increased after 24 h (t4 = − 4.387, P = 0.012) of S. graminum feeding and reached a peak at 48 h (t4 = − 9.981, P = 0.001) (Fig. 8). The activity of POD in leaves was also significantly increased 72 h after S. avenae feeding (t4 = − 3.353, P = 0.028). Furthermore, the activities of SOD, CAT and APx were significantly increased after 24 h of S. graminum feeding (t4 = − 12.295, P < 0.001; t4 = − 2.789, P = 0.049; t4 = − 7.761, P = 0.001), whereas S. avenae feeding had no significant effects on the activities of SOD, CAT and APx (t4 = − 1.560, P = 0.194; t4 = − 0.600, P = 0.581; t4 = 0.048, P = 0.964).
Cytological examination of callose deposition and ROS accumulation in aphid-infested leaves
To detect whether callose was deposited at the feeding sites, aphid-infested leaves were stained with aniline blue. As shown in Fig. 9, no callose deposits were observed in vascular tissues without aphid infestation (Fig. 9a). However, in S. avenae- and S. graminum-infested tissues, callose deposits were clearly detected as bright blue fluorescence directly at the feeding sites (Fig. 9b and c).
H2O2 accumulation has been shown to be induced by wounding and by pathogen and herbivore attacks in plants and is involved in plant defence responses as a signal molecule . To record the accumulation of H2O2 after aphid infestation, S. avenae- and S. graminum-infested leaves were examined after cytological staining with DAB, which was used to detect the production of H2O2. As shown in Fig. 9d-e, no obvious DAB staining was observed in the non-infested leaves, and a small brown-stained area was detected in the wheat leaves after S. avenae feeing. However, H2O2 was clearly detected in the areas of S. graminum feeding, which indicated that S. graminum feeding induced a massive accumulation of H2O2 (Fig. 9f). The H2O2 contents in the feeding sites of wheat leaves after S. avenae and S. graminum infestations were also examined. As indicated in Fig. 10, the concentration of H2O2 in wheat leaves infested with S. graminum (143.19 ± 31.15 μmol g− 1 FW; F2,6 = 7.345; P = 0.024) was significantly higher than that in S. avenae-infested and control leaves. In contrast, S. avenae feeding had no significant effects on the H2O2 content compared with the control (Fig. 10). The changes in the H2O2 content in response to aphid feeding were consistent with the DAB staining results.
Scavenging of H2O2 using DMTU reduces S. graminum feeding-induced damage on wheat leaves
To further investigate the role of H2O2 accumulation on the damage induced by S. graminum feeding, wheat seedlings infested with aphids were treated with 5 mM DMTU (an H2O2 scavenger). The DAB staining results shown in Fig. 11 demonstrated that DMTU treatment inhibited the S. graminum feeding-induced production of H2O2 in wheat leaves and the symptoms of damage in wheat leaves caused by S. graminum feeding. The delayed fluorescence and chlorophyll content were also assessed, and the results showed that the DMTU-treated infested leaves showed decreased chlorophyll degradation and that the chlorophyll content in the DMTU-treated leaves was significantly higher than that in the non-DMTU-treated leaves infested with S. graminum (F2, 6 = 13.93, P = 0.0056).
A previous study showed that S. graminum feeding led to obvious feeding damage and the loss of chlorophyll in aphid-susceptible winter wheat accession Beijing 837 . Similarly, serious chlorosis symptoms were observed on another winter wheat accession, Zhongmai 175, after 48 h of S. graminum feeding in this study, and this effect was also accompanied by a significant reduction in the total chlorophyll content of the wheat leaves, further demonstrating the phytotoxic effects of S. graminum on susceptible wheat plants. To further compare the similarities and differences between the responses to S. graminum and S. avenae feeding at the molecular level, a comparative transcriptome analysis of wheat leaves after aphid feeding was performed. We found that more than 20,000 genes were significantly regulated in wheat infested with S. graminum, but only 1700 genes were significantly modulated after 48 h of S. avenae feeding, which indicated that the physiological changes induced by S. graminum are notably different from those induced by S. avenae and that various metabolic pathways are involved in the development of damage caused by S. graminum feeding. Moreover, many genes involved in plant photosynthesis were strongly downregulated after S. graminum feeding, and this finding provides molecular evidence showing that chlorosis is induced by S. graminum.
S. graminum feeding induces stronger plant defence responses than S. avenae
Piercing-sucking hemipteran insects, such as aphids and whiteflies, mainly induce SA-mediated defence signal pathways [28, 29]. However, some studies have also demonstrated that genes involved in both the JA and SA defence response pathways, such as LOX, PIs, PAL, and PR1, are significantly upregulated in response to aphid feeding [20, 30, 31]. Similarly, we found that both S. avenae and S. graminum feeding significantly increased the expression levels of genes related to the SA, JA and ET signalling pathways. Plant defence responses activated by aphids are closely associated with the plant species, aphid density and infestation time [32, 33]. In the future, we will further identify the wheat defence responses induced by S. graminum and S. avenae feeding under various aphid densities and feeding periods. Although both of these cereal aphids induced both SA-, JA- and ET-dependent defence pathways, S. graminum feeding induced the expression of more genes involved in plant defence pathways in this study. For example, 30 PAL genes were upregulated in response to S. graminum feeding, but only four PAL genes were upregulated in response to S. avenae feeding. The transcript levels of five ET-responsive genes were upregulated in response to S. graminum, but only one was modulated in response to S. avenae feeding. Zhang et al. demonstrated that the fold changes in the expression levels of PR genes and the SA contents in wheat leaves were significantly greater after S. graminum feeding than after S. avenae feeding . Argandoña et al. also suggested that S. graminum induced more ethylene production than the non-phytotoxic aphid Rhopalosiphum padi . The stronger defence responses activated by S. graminum feeding than by S. avenae feeding might be responsible for the induction of chlorosis in wheat.
Genes encoding plant cell-modifying proteins are significantly upregulated in response to S. graminum feeding
Multiple modifications can be triggered in cell walls in response to microbial and insect attack . Callose deposition in crop plants is observed in response to biotrophic fungal infection at papillae sites and in sieve elements in response to aphids [36, 37]. It has been proposed that callose deposition impedes fungal attacks at the sites of attempted penetration in epidermal cells and thereby supports pathogen resistance . The transcript levels of callose synthases were significantly upregulated in response to S. graminum and S. avenae feeding. Obvious callose accumulation was also observed in wheat leaves after S. avenae and S. graminum feeding. However, the role of callose in aphid-plant interactions remains unknown. It has been hypothesized that callose deposition induced by aphids is involved in the sealing of the sieve pores as a phloem defence mechanism that impedes mass flow and prevents the flow of nutrients to piercing-sucking insects .
Many studies have demonstrated that the damage symptoms induced by pathogens and herbivores in plants are caused by the secretion of plant cell wall enzymes during the process of pathogen infection and mirid Lygus hesperus feeding [40, 41]. Additionally, the activities of PEMs and PGs have also been detected in S. graminum watery saliva, and the injection of these commercial enzymes in plant leaves causes damage symptoms similar to those induced by S. graminum feeding [42, 43]. However, pectinase activity has also been detected in the saliva of non-phytotoxic aphids, such as S. avenae and A. pisum [43, 44]. Transmission electron microscopy has shown that stylets predominantly penetrate between the layers of cellulose fibres and not via the middle lamella pectin layer . The role of the saliva pectinases of S. graminum in the induction of chlorosis remains unclear. Interestingly, in our study, many plant-derived enzymes and proteins involved in plant cell wall modifications, such as PGs, PEMs and expansins, were induced by S. graminum feeding, but none were induced by S. avenae feeding. The upregulation of PG and PEM activity might result in the degradation of the cell wall around aphid feeding sites. Plant cells exploit complicated mechanisms for sensing the loss of cell wall integrity (CWI) during biotic stress and activate a variety of defence responses . For instance, the production of oligogalacturonic acid (OGA) fragments derived from the degradation of plant cell walls has been shown to trigger oxidative bursts, hypersensitive responses (HR), and other downstream defence responses in many plant species as a host-derived DAMP [47, 48], and these effects might further promote the induction of damage symptoms. In addition, ethylene production has been shown to be involved in the induction of plant cell wall-modifying proteins and the death of plant tissues. The role of PCMDPs and the ethylene pathway in the induction of chlorosis symptoms by S. graminum needs to be further investigated.
S. graminum feeding induces strong ROS-scavenging activity in wheat leaves
Aphid feeding usually leads to oxidative stress in host plants . Oxidative stress is controlled by cellular antioxidant mechanisms in which multiple enzymatic scavengers, such as POD, APx, and CAT, are utilized by the cell to limit damage from reactive oxygen species [50, 51]. The transcriptomic and enzymatic results showed that both S. graminum and S. avenae feeding increased the transcript levels and enzyme activities of ROS scavengers in wheat leaves. However, the expression levels and activity of antioxidants, particularly POD, induced by S. graminum feeding were notably higher than those induced by S. avenae, which suggested that S. graminum infestation results in strong oxidative stress and substantial H2O2 accumulation. Although ROS scavengers were significantly upregulated in response to S. graminum infestation, ROS production can exceed the cellular antioxidant capacity, resulting in oxidative damage to cellular components and cell death in leaves.
Induction of high H2O2 accumulation by S. graminum feeding is involved in leaf chlorosis
H2O2 is involved in the activation of HR, which is characterized by the rapid death of cells in the region surrounding the site of pathogen infection site [52, 53]. To further investigate the roles of H2O2 in the induction of feeding damage caused by S. graminum feeding, the accumulation of H2O2 in wheat leaves was detected. We found that S. graminum feeding induced the obvious accumulation of H2O2 at feeding sites, but S. avenae feeding had no significant effects on H2O2 production. In addition, seedlings treated with the H2O2 scavenger DMTU showed reductions in the chlorosis symptoms and chlorophyll loss triggered by S. graminum feeding. These results demonstrate that H2O2 accumulation plays important roles in the induction of chlorosis in wheat leaves in response to S. graminum feeding.
Aphid saliva is known to be involved in the induction of plant defence responses, and the eliciting activity of watery saliva of other aphid species such as M. persicae and S. avenae has been investigated [54, 55]. The transient overexpression of Mp10, a salivary protein of M. persicae, induces plant defence responses and obvious chlorosis in N. benthamiana . Specific elicitors or pathogen-like toxins in S. graminum saliva are likely involved in the induction of chlorosis. Comparative analyses of the salivary proteomes of four differentially virulent S. graminum biotypes revealed six salivary proteins with significant proteomic variation, and these proteins might thus be involved in the induction of feeding damage in plants . Further research is required to identify the virulence factors in the salivary proteins of S. graminum and the mechanism underlying the induction of chlorosis.
In summary, the transcriptomic profiling of wheat performed in this study revealed similarities and differences among the responses of wheat to feeding by the phytotoxic aphid S. graminum and the non-phytotoxic aphid S. avenae. Both aphids induced the JA, SA and ET signalling pathways, but S. graminum triggered stronger plant defence responses and greater ROS-scavenging activity than S. avenae. A cytological analysis showed that aphid feeding induced callose deposition in wheat leaves and that substantial H2O2 accumulation was induced by S. graminum feeding. Our results also demonstrated that H2O2 plays vital roles in the induction of chlorosis in wheat leaves in responses to S. graminum feeding. Our future studies will focus on the mechanisms of H2O2 accumulation induced by S. graminum feeding and the roles of salivary proteins of S. graminum in the induction of chlorosis symptoms in wheat.
Plants and aphids
Seeds of Triticum aestivum var. Zhongmai 175 were germinated in distilled water for 3–4 days at a temperature of 25 ± 1 °C in a Petri dish. Healthy seedlings of similar sizes were planted in 7.2 × 7.2 cm plastic plots filled with organic soil and grown under controlled environmental conditions in climate chambers with a temperature of 20 ± 1 °C, a 40–60% relative humidity and a 14-h-light/10-h-dark photoperiod. Clones of S. graminum and S. avenae were maintained on the wheat plants (Zhongmai 175) as described previously .
At the two-leaf stage (12-day old plants), 20 apterous adult S. graminum or S. avenae were confined on the first leaf of wheat seedlings using a clip cage as described previously . New-born nymphs produced by aphid adults were carefully removed every 12 h using a brush. After 48 h of feeding, all the aphids were removed, and leaf tissues of approximately 2.5 × 2.5 cm from the aphid feeding sites of each plant were harvested flash frozen with liquid nitrogen and stored at − 80 °C until further processing for RNA extraction. Detection of delayed fluorescence and histological staining were conducted immediately after sample collection. Three leaf sections covering the aphid feeding sites were collected from three independent plants and pooled to form one biological replicate. Three biological replicates were performed for each treatment.
Changes in chlorophyll levels in wheat leaves after aphid infestation
Delayed fluorescence is associated with extremely weak light emitted by chlorophyll molecules in plants and can reflect the chlorophyll content, providing a powerful tool for studying stress reactions in plants. The chlorophyll content in wheat leaves were each infested by 20 aphids as described above was first detected using the NightShade LB 985 In vivo Plant Imaging System (Berthold Technologies, Bad Wildbad, Germany). After 48 h of aphid feeding, the leaves were cut and immediately illuminated for 30 s with an LED panel. After the light was switched off, the delayed fluorescence was measured immediately using the NightShade system. The exposure time was set to 30 s using 4-by-4 pixel binning. The total chlorophyll content in wheat leaves after aphid infestation was also examined using a Chlorophyll Assay Kit (Solarbio, Beijing, China) according to the manufacturer’s instructions. In brief, 0.1 g of fresh leaf tissues was ground to a fine powder and extracted with 2 mL of 80% acetone (v/v) at 4 °C for overnight. The homogenate was centrifuged at 4000 g for 10 min at 4 °C, and the supernatant was used for the chlorophyll assay. The amounts of chlorophyll were detected spectrophotometrically, by reading the absorbance at 645 and 663 nm (DU800, Beckman, USA), and then calculated as described previously .
RNA preparation and sequencing library construction
Wheat leaves were first infested with aphids for 48 h as described above. The total RNA from the wheat leaves was extracted with the TRIzol reagent (Invitrogen) according to the manufacturer’s recommended protocol. The RNA concentration was measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Inc., USA), and the RNA integrity was evaluated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Samples with an RNA integrity number (RIN) ≥ 7.0 were used in the subsequent analysis. Libraries were constructed using the TruSeq Stranded mRNA LT Sample Prep Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions and were sequenced on the Illumina sequencing platform (Illumina HiSeq 4000), which generated 150-bp paired-end reads were generated.
RNA sequencing and data analysis
The raw data (raw reads) were filtered to obtain high-quality reads by removing the reads containing adaptor sequences, more than 3% ambiguous bases (noted as N) or more than 50% low-quality bases (Phred quality score Q < 30). The resulting high-quality clean reads were then mapped to the reference genome (https://plants.ensembl.org/Triticum_aestivum/Info/Index) using TopHat2 with the default values . The fragments per kilobase of exon model per million mapped reads (FPKM) data were used to estimate the transcript expression levels in all the samples, and the genes with more than 1 FPKM in at least one sample of wheat leaves were used for further analysis . The differentially expressed genes (DEGs) between the control and treated samples were screened using DESeq based on the following criteria: adjusted p value (q value) threshold < 0.00001 and |log2 FC| ≥ 1. A PCA of nine different leaf samples was performed based on pairwise comparisons using the DESeq package of R . GO enrichment and KEGG pathway enrichment analyses of the DEGs were performed using R based on the hypergeometric distribution . To annotate the functions of the transcripts, the unigenes were blasted against the Nr database using the BLAST programme with an E-value ≤1e-5.
Assays of antioxidant enzymes in wheat leaves after aphid infestation
The leaves were infested with 20 apterous adult S. graminum or S. avenae as described previously. The activities of peroxidase (POD), superoxide dismutase (SOD), catalase (CAT) and ascorbate peroxidase (APx) in the wheat leaves at different time points after the aphid infestation (12 h, 24 h, 48 h and 72 h) were examined using corresponding kits (Jiancheng Bioengineering Institute, Nanjing, China) according to the manufacturer’s instructions. Briefly, 0.2 g of fresh leaf tissues was ground with 1.5 mL of ice-cooled 50 mM Na-phosphate buffer (pH = 7.8) containing 0.1 mM EDTA and 1.0% (w/v) polyvinylpyrrolidone. The homogenate was centrifuged at 15,000 g for 30 min at 4 °C, and the supernatant was immediately collected for further enzyme assays. The activities of POD, SOD, CAT, and APx were measured by following the changes in absorbance at 470 nm, 560 nm, 240 nm and 290 nm, respectively, according to a previous study .
Detection of H2O2 and callose accumulation in wheat leaves after aphid infestation
The detection of H2O2 in wheat leaves by 3′-diaminobenzidine (DAB) staining was performed according to the histochemical methods described by Wang et al.  with some modifications. In brief, leaf segments previously infested with 20 apterous adults of S. avenae or S. graminum were immersed in 1 mg mL− 1 DAB solution (10 mmol L− 1 Na2HPO4, pH 3.8), and incubated in the dark overnight at room temperature. Then, the leaves were decolorized in boiling 95% ethanol for 10 min and hyalinized in saturated chloral hydrate. The stained leaves were imaged using an Olympus BX-63 microscope (Olympus Corporation, Japan). The endogenous H2O2 content in the wheat leaves after aphid feeding was determined using the protocols reported by Ferguson et al. . For the visualization of callose, the leaves were first fixated, destained overnight in 1:3 acetic acid/ethanol (v/v) solution and washed in 150 mM K2HPO4 for 30 min. The leaves were subsequently incubated for 6 h with 150 mM K2HPO4 and 0.01% aniline blue for staining, and the callose depositions were observed and photographed with an Olympus SZX-16 fluorescence microscope (Olympus Corporation, Japan) using a DAPI filter.
Wheat seedlings treated with DMTU solution
The leaves of wheat seedlings were treated with 5 mM dimethylthiourea (DMTU, a scavenger of H2O2) solution or deionized water (control) for 24 h and then infested with 20 apterous adults of S. graminum for 48 h. Assessments of DAB staining and delayed fluorescence analyses and an assessment of the chlorophyll content of the wheat leaves were then performed as described previously.
All the data were analysed using SPSS Statistics 20.0 software (SPSS Inc., Chicago, IL., USA), and the differences between or among groups were examined through an independent-samples t-test or one-way analysis of variance (Duncan). P values less than 0.05 were considered statistically significant.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
1-aminocyclopropane-1-carboxylate oxidase 1
Allene oxide cyclase
Allene oxide synthase
Barley yellow dwarf virus
Cell wall integrity
Damage associated molecular patterns
Differentially expressed genes
Adjusted p value
Fragments per kilobase of exon model per million mapped reads
- H2O2 :
Herbivory associated molecular patterns
- PR protein:
High-throughput RNA sequencing
Reactive oxygen species
Ribulose bisphosphate carboxylase oxygenase
Howe GA, Jander G. Plant immunity to insect herbivores. Annu Rev Plant Biol. 2008;59:41–66.
Bonaventure G. Perception of insect feeding by plants. Plant Biol. 2012;14(6):872–80.
Steppuhn A, Baldwin IT. Resistance management in a native plant: nicotine prevents herbivores from compensating for plant protease inhibitors. Ecol Lett. 2007;10(6):499–511.
Hartl M, Giri AP, Kaur H, Baldwin IT. Serine protease inhibitors specifically defend Solanum nigrum against generalist herbivores but do not influence plant growth and development. Plant Cell. 2010;22(12):4158–75.
Paré PW, Tumlinson JH. Plant volatiles as a defense against insect herbivores. Plant Physiol. 1999;121(2):325–32.
Kessler A, Baldwin IT. Defensive function of herbivore-induced plant volatile emissions in nature. Science. 2001;291(5511):2141–4.
Leitner M, Boland W, Mithöfer A. Direct and indirect defences induced by piercing-sucking and chewing herbivores in Medicago truncatula. New Phytol. 2005;167(2):597–606.
Heidel AJ, Baldwin IT. Microarray analysis of salicylic acid-and jasmonic acid-signalling in responses of Nicotiana attenuata to attack by insects from multiple feeding guilds. Plant Cell Environ. 2004;27(11):1362–73.
Mithöfer A, Wanner G, Boland W. Effects of feeding Spodoptera littoralis on lima bean leaves. II. Continuous mechanical wounding resembling insect feeding is sufficient to elicit herbivory-related volatile emission. Plant Physiol. 2005;137(3):1160–8.
Vos IA, Verhage A, Schuurink RC, Watt LG, Pieterse CMJ, van Wees SCM. Onset of herbivore-induced resistance in systemic tissue primed for jasmonate-dependent defenses is activated by abscisic acid. Front Plant Sci. 2013;4:539.
Kaloshian I, Walling LL. Hemipterans as plant pathogens. Annu Rev Phytopathol. 2005;43:491–521.
Coppola V, et al. Transcriptomic and proteomic analysis of a compatible tomato-aphid interaction reveals a predominant salicylic acid-dependent plant response. BMC Genomics. 2013;14(1):515.
Donovan MP, Nabity PD, DeLucia EH. Salicylic acid-mediated reductions in yield in Nicotiana attenuata challenged by aphid herbivory. Arthropod Plant Interact. 2013;7(1):45–52.
Hu XS, Liu XF, Thieme T, Zhang GS, Liu TX, Zhao HY. Testing the fecundity advantage hypothesis with Sitobion avenae, Rhopalosiphum padi, and Schizaphis graminum (Hemiptera: Aphididae) feeding on ten wheat accessions. Sci Rep. 2015;5:18549.
Blackman RL, Eastop VF. Aphids on the world’s crops: an identification and information guide. New York: John Wiley & Sons Ltd.; 2000.
Miles PW. Aphid saliva. Biol Rev. 1999;74:41–85.
Al-Mousawi AH, Richardson PE, Burton RL. Ultrastructural studies of greenbug (Hemiptera: Aphididae) feeding damage to susceptible and resistant wheat cultivars. Ann Entomol Soc Am. 1983;76(6):964–71.
Zhang Y, Fan J, Fu Y, Francis F, Chen JL. Plant-mediated interactions between two cereal aphid species: promotion of aphid performance and attraction of more parasitoids by infestation of wheat with phytotoxic aphid Schizaphis graminum. J Agr Food Chem. 2019;67(10):2763–73.
Zhao LY, Chen JL, Cheng DF, Sun JR, Liu Y, Tian Z. Biochemical and molecular characterizations of Sitobion avenae-induced wheat defense responses. Crop Prot. 2009;28(5):435–42.
Zhu-Salzman K, Salzman RA, Ahn JE, Koiwa H. Transcriptional regulation of sorghum defense determinants against a phloem-feeding aphid. Plant Physiol. 2004;134(1):420–31.
Park SJ, Huang Y, Ayoubi P. Identification of expression profiles of sorghum genes in response to greenbug phloem-feeding using cDNA subtraction and microarray analysis. Planta. 2006;223(5):932–47.
Sandström J, Telang A, Moran N. Nutritional enhancement of host plants by aphids-a comparison of three aphid species on grasses. J Insect Physiol. 2000;46(1):33–40.
Machado-Assefh CR, Lucatti AF, Alvarez AE. Induced senescence promotes the feeding activities and nymph development of Myzus persicae (Hemiptera: Aphididae) on potato plants. J Insect Sci. 2014;14:155.
Morkunas I, Mai VC, Gabryś B. Phytohormonal signaling in plant responses to aphid feeding. Acta Physiol Plant. 2011;33(6):2057–73.
Bi JL, Felton GW. Foliar oxidative stress and insect herbivory: primary compounds, secondary metabolites, and reactive oxygen species as components of induced resistance. J Chem Ecol. 1995;21(10):1511–30.
Blokhina O, Virolainen E, Fagerstedt KV. Antioxidants, oxidative damage and oxygen deprivation stress: a review. Ann Bot. 2002;91(2):174–94.
Choudhury S, Panda P, Sahoo L, Panda SK. Reactive oxygen species signaling in plants under abiotic stress. Plant Signal Behav. 2013;8(4):e23681.
Züst T, Agrawal AA. Mechanisms and evolution of plant resistance to aphids. Nat Plants. 2016;2(1):15206.
Wang XW, Li P, Liu SS. Whitefly interactions with plants. Curr Opin Insect Sci. 2017;19:70–5.
Reddy SK, Weng Y, Rudd JC, Akhunova A, Liu SY. Transcriptomics of induced defense responses to greenbug aphid feeding in near isogenic wheat lines. Plant Sci. 2013;212:26–36.
Moran PJ, Thompson GA. Molecular responses to aphid feeding in Arabidopsis in relation to plant defense pathways. Plant Physiol. 2001;125(2):1074–85.
Stewart SA, Hodge S, Bennett M, Mansfield JW, Powell G. Aphid induction of phytohormones in Medicago truncatula is dependent upon time post-infestation, aphid density and the genotypes of both plant and insect. Arthropod Plant Interact. 2016;10(1):41–53.
Hodge S, Bennett M, Mansfield JW, Powell G. Aphid-induction of defence-related metabolites in Arabidopsis thaliana is dependent upon density, aphid species and duration of infestation. Arthropod Plant Interact. 2019;13(3):387–99.
Argandoña VH, Chaman M, Cardemil L, Muñoz O, Zúñiga GE, Corcuera LJ. Ethylene production and peroxidase activity in aphid-infested barley. J Chem Ecol. 2001;27(1):53–68.
Malinovsky FG, Fangel JU, Willats WGT. The role of the cell wall in plant immunity. Front Plant Sci. 2014;5:178.
Voigt CA. Callose-mediated resistance to pathogenic intruders in plant defense-related papillae. Front Plant Sci. 2014;5:168.
Saheed SA, Cierlik I, Larsson KAE, Delp G, Bradley G, Jonsson LM, Botha CE. Stronger induction of callose deposition in barley by Russian wheat aphid than bird cherry-oat aphid is not associated with differences in callose synthase or β-1, 3-glucanase transcript abundance. Physiol Plantarum. 2009;135(2):150–61.
Luna E, Pastor V, Robert J, Flors V, Mauch-Mani B, Ton J. Callose deposition: a multifaceted plant defense response. Mol Plant-Microbe Interact. 2011;24(2):183–93.
Will T, van Bel AJE. Physical and chemical interactions between aphids and plants. J Exp Bot. 2006;57(4):729–37.
Kubicek CP, Starr TL, Glass NL. Plant cell wall–degrading enzymes and their secretion in plant-pathogenic fungi. Annu Rev Phytopathol. 2014;52:427–51.
de la Paz C-MM, et al. Polygalacturonase causes lygus-like damage on plants: cloning and identification of western tarnished plant bug (Lygus hesperus) polygalacturonases secreted during feeding. Arthropod Plant Interact. 2008;2(4):215.
Ma RZ, Reese JC, Black WC IV, Bramel-Cox P. Chlorophyll loss in a greenbug-susceptible sorghum due to pectinases and pectin fragments. J Kansas Entomol Soc. 1998;71(1):51–60.
Cherqui A, Tjallingii WF. Salivary proteins of aphids, a pilot study on identification, separation and immunolocalisation. J Insect Physiol. 2000;46(8):1177–86.
Liu Y, Wang WL, Guo GX, Ji XL. Volatile emission in wheat and parasitism by Aphidius avenae after exogenous application of salivary enzymes of Sitobion avenae. Entomol Exp Appl. 2009;130(3):215–21.
Tjallingii WF, Esch TH. Fine structure of aphid stylet routes in plant tissues in correlation with EPG signals. Physiol Entomol. 1993;18(3):317–28.
Hématy K, Cherk C, Somerville S. Host–pathogen warfare at the plant cell wall. Curr Opin Plant Biol. 2009;12(4):406–13.
Ferrari S, Savatin DV, Sicilia F, Gramegna G, Cervone F, De Lorenzo G. Oligogalacturonides: plant damage-associated molecular patterns and regulators of growth and development. Front Plant Sci. 2013;4:49.
Lee S, et al. Oligogalacturonic acid and chitosan reduce stomatal aperture by inducing the evolution of reactive oxygen species from guard cells of tomato and Commelina communis. Plant Physiol. 1999;121(1):147–52.
Bednarski W, Borowiak-Sobkowiak B, Wilkaniec B, Samardakiewicz S, Morkunas I. Oxidative stress in pea seedling leaves in response to Acyrthosiphon pisum infestation. Phytochemistry. 2013;93:49–62.
Mittler R. Oxidative stress, antioxidants and stress tolerance. Trends Plant Sci. 2002;7(9):405–10.
Demidchik V. Mechanisms of oxidative stress in plants: from classical chemistry to cell biology. Environ Exp Bot. 2015;109:212–28.
Thordal-Christensen H, Zhang Z, Wei Y, Collinge DB. Subcellular localization of H2O2 in plants. H2O2 accumulation in papillae and hypersensitive response during the barley-powdery mildew interaction. Plant J. 1997;11(6):1187–94.
Govrin EM, Levine A. The hypersensitive response facilitates plant infection by the necrotrophic pathogen Botrytis cinerea. Curr Biol. 2000;10(13):751–7.
De Vos M, Jander G. Myzus persicae (green peach aphid) salivary components induce defence responses in Arabidopsis thaliana. Plant Cell Environ. 2009;32(11):1548–60.
Zhang Y, Fan J, Francis F, Chen JL. Watery saliva secreted by the grain aphid Sitobion avenae stimulates aphid resistance in wheat. J Agr Food Chem. 2017;65(40):8798–805.
Bos JIB, Prince D, Pitino M, Maffei ME, Win J, Hogenhout SA. A functional genomics approach identifies candidate effectors from the aphid species Myzus persicae (green peach aphid). PLoS Genet. 2010;6(11):e1001216.
Nicholson SJ, Puterka GJ. Variation in the salivary proteomes of differentially virulent greenbug (Schizaphis graminum Rondani) biotypes. J Proteome. 2014;105:186–203.
Arnon DI. Copper enzymes in isolated chloroplasts: polyphenoloxidase in Beta vulgaris. Plant Physiol. 1949;24(1):1–15.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletion and gene fusions. Genome Biol. 2013;14(4):36.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.
Anders S, Huber W. Differential expression of RNA-Seq data at the gene level–the DESeq package. Heidelberg: European Molecular Biology Laboratory (EMBL); 2012.
Kanehisa M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:480–4.
Jiang L, Yang H. Prometryne-induced oxidative stress and impact on antioxidant enzymes in wheat. Ecotoxicol Environ Saf. 2009;72(6):1687–93.
Wang CF, Huang LL, Buchenauer H, Han QM, Zhang HC, Kang ZS. Histochemical studies on the accumulation of reactive oxygen species (O2 − and H2O2) in the incompatible and compatible interaction of wheat - Puccinia striiformis f. sp. tritici. Physiol Mol Plant Pathol. 2007;71(4–6):230–9.
Ferguson IB, Watkins CB, Harman JE. Inhibition by calcium of senescence of detached cucumber cotyledons effect on ethylene and hydroperoxide production. Plant Physiol. 1983;71(1):182–6.
We would like to thank technician Ms. Yanxia Liu for aphid rearing in our laboratory.
This study was funded by the National Natural Science Foundation of China (31871979, 31901881), National Key R&D Plan in China (2017YFD0201701, 2016YFD0300701), and the Cooperation Project between Belgium and China from MOST (2014DF32270). The funding organizations played no role in the design of study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The seeds of wheat, Zhongmai 175 were bred by institute of Crop Sciences, Chinese Academy of Agricultural and Sciences and bought from Henan Shengyuan Seed Industry Technology Co., Ltd., Xuchang, China. Cereal aphids Sitobion avenae and Schizaphis graminum used in our study were collected in Langfang, Hebei Province, China, and has been rear in our greenhouse for more than 9 years. The authors declare that all the experiments performed in this study comply with the institutional, national, or international guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of transcriptome data.
List of all DEGs after infestation by aphids in comparison with controls.
DEGs involved in SA, JA and ET- dependent defence pathways in response to aphid feeding.
Transcript levels of genes involved in ROS scavenge after aphid feeding.
About this article
Cite this article
Zhang, Y., Fu, Y., Fan, J. et al. Comparative transcriptome and histological analyses of wheat in response to phytotoxic aphid Schizaphis graminum and non-phytotoxic aphid Sitobion avenae feeding. BMC Plant Biol 19, 547 (2019). https://doi.org/10.1186/s12870-019-2148-5