- Research article
- Open Access
Effects of sugarcane aphid herbivory on transcriptional responses of resistant and susceptible sorghum
BMC Genomics volume 19, Article number: 774 (2018)
Sugarcane aphid (Melanaphis sacchari) outbreaks in sorghum that were first reported in 2013 are now the most significant threat to this crop in all major sorghum production areas in the U.S. The outcomes of interactions between sugarcane aphid and sorghum and thus the severity of the outbreaks depend on sorghum genotype and potentially also on the phenology of sorghum. Mechanisms underlying these interactions are not known, however. Thus, the goal of this research was to characterize transcriptional changes in a commercially available resistant and a susceptible genotype of sorghum at 2- and 6-wk post-emergence exposed to M. sacchari herbivory. The effects of sorghum age and genotype on the daily change in aphid densities were also evaluated in separate greenhouse experiments.
A higher number of diffentially expressed genes (DEGs) was recovered from the 2-wk plants exposed to aphid herbivory compared to the 6-wk plants across genotypes. Further, gene ontology and pathway analysis indicated a suite of transcriptional changes in the resistant genotype that were weak or absent in the susceptible sorghum. Specifically, the aphid-resistant genotype exposed to M. sacchari up-regulated several genes involved in defense, which was particularly evident in the 2-wk plants that showed the most robust transcriptional responses. These transcriptional changes in the younger resistant sorghum were characterized by induction of hormone-signaling pathways, pathways coding for secondary metabolites, glutathion metabolism, and plant-pathogen interaction. Furthermore, the 2-wk resistant plants appeared to compensate for the effects of oxidative stress induced by sugarcane aphid herbivory with elevated expression of genes involved in detoxification. These transcriptional responses were reflected in the aphid population growth, which was significantly faster in the susceptible and older sorghum than in the resistant and younger plants.
This experiment provided the first insights into molecular mechanisms underlying lower population growth of M. sacchari on the resistant sorghum genotype. Further, it appears that the younger resistant sorghum was able to mount a robust defense response following aphid herbivory, which was much weaker in the older sorghum. Several pathways and specific genes provide specific clues into the mechanisms underlying host plant resistance to this invasive insect.
Plants have evolved diverse strategies to defend, tolerate, or avoid insect herbivores [1, 2] that provide valuable resources for deploying in pest management [3, 4]. When plants detect the injury inflicted by herbivores, they fine-tune their responses through a flexible and finely balanced network that in turn activates specific defense responses [4, 5] with a large intraspecific variation in signaling cascades and production of metabolites . The genetic diversity of herbivore-resistant traits is also a valuable resource for studying the mechanism of resistance and the ultimate deployment of them in crop breeding programs [4, 6]. In fact, resistant cultivars and hybrids have been employed as an effective method for controlling economically damaging arthropods in many different crops, including sorghum [Sorghum bicolor (L.) Moench] [7,8,9,10,11], which is vulnerable to significant injury caused by diverse insect pests . The employment of host-plant resistance in controlling greenbug, Schizaphis graminum Rodandi (Hemiptera: Aphididae), in sorghum is an example of an effective integrated pest management in this crop [13,14,15,16]. As a result of the selected resistance sources that were bred into commercial sorghum hybrids, the frequency of chemical control to suppress the pest was significantly reduced [17, 18].
Invasive pests pose unique challenges to crop protection, which was recently exemplified by an invasion of sugarcane aphid, Melanaphis sacchari Zehtner (Hemiptera: Aphididae), a new invasive pest threatening sorghum production in the U.S. [19, 20]. While the species has been present in the U.S. in sugarcane for over a century, an invasion of a new haplotype  absent in the country prior to the outbreaks is thought to be responsible for the eruptive population dynamics of M. sacchari in sorghum, which were first reported in 2013 [22, 23]. Similar to many other invasive insects [24,25,26], the genetic diversity of M. sacchari is generally low [27, 28], and sugarcane aphids colonizing sorghum in the U.S. appear to be genetically identical .
Current tactics used to suppress sugarcane aphids rely primarily on insecticides, but the rapid aphid population growth and high dispersion can hinder the effectiveness of the insecticides . Therefore, host plant resistance has been successfully deployed to alleviate the impact of sugarcane aphids on sorghum and has proven to be an effective suppression tactic in the field [30,31,32]. In fact, several sorghum varieties have high levels of resistance to sugarcane aphids [22, 33]. For example, Szczepaniec  reported that host plant resistance in the resistant genotype was the most effective factor in population dynamics of M. sacchari when compared to insecticide seed treatment or planting date. Specifically, densities of aphids were 2–2.5 times greater on the susceptible sorghum compared to the resistant hybrid. Further, it has been demonstrated in the field that populations of these aphids can reach exponential rate of increase within two weeks of colonizing post-bloom sorghum . While suppressive effects of host plant resistance on aphid population growth were not surprising, the differences in aphid population growth between younger sorghum and sorghum close to reproductive stage have important consequences to ecology and management of the aphid. Sugarcane aphids colonize sorghum in vegetative stages in the southern U.S., where they can overwinter on non-crop vegetation , while wind-aided migration is responsible for aphid colonization of sorghum in reproductive stages in the northern states (e.g., High Plains) . Mechanisms underlying sorghum resistance to sugarcane aphids across sorghum phenology are not well understood, however.
The molecular mechanisms underlying plant-aphid interactions, in general, have been studied extensively using molecular mapping and QTLs screening [34,35,36]. Over the last few years, however, as it became apparent that transcriptional re-programming is the backbone of many plant defense responses, the research has shifted toward transcriptomic analysis of herbivore-plant interactions. Aphid-resistant cultivars respond to aphid feeding by up-regulation of diverse genes that contribute to aphid resistance [37,38,39]. These genes have been found to be involved in phytohormone signaling pathways, particularly salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) [40,41,42]. Further, the production of plant secondary metabolites which are compounds that play a key role in defense and are often toxic to insects is another common induced response of plants exposed to insect herbivory [43, 44].
A comprehensive understanding of the molecular interactions between the susceptible and resistant sorghum and sugarcane aphid can provide insights into the mechanism of resistance and advance the potential to apply them to crop protection. Therefore, the goal of this research was to characterize transcriptional responses of a resistant and susceptible sorghum to herbivory by M. sacchari. Further, owing to the fact that these aphids can colonize sorghum at vegetative or at reproductive stages of sorghum development depending on the geographic location, we also explored if the interaction between sorghum resistance and sorghum age had an effect on gene expression and on population growth of the aphids. We hypothesized that expression of genes involved in defense pathways will be greater in the resistant and in younger sorghum than in the susceptible and in older sorghum and that the aphid population growth will reflect these differences.
Results and discussion
Mapping results and DEGs in resistant and susceptible sorghum exposed to aphid herbivory
Sequencing of 24 libraries generated 51.6 to 82.8 million reads from individual samples, and 43.9 to 68.8 million reads were uniquely mapped to sorghum reference genome sbicolor_454 v3.1.1. On average, 91% of the reads mapped to exonic regions (Additional file 1). A higher number of DEGs was recovered from the 2-wk plants exposed to aphid herbivory compared to the 6-wk plants across genotypes. The number of DEGs was also affected by genotype such that the resistant genotype had a higher number of DEGs in response to aphid herbivory at both 2- and 6-wk post-emergence than the susceptible genotype (Fig. 1). Of the 5955 DEGs in the 2-wk resistant genotype exposed to aphid herbivory, 54.9% were up-regulated and 45% were down-regulated. On the other hand, 62.6% DEGs were up-regulated and 37.3% were down-regulated in the 2-wk susceptible sorghum in response to the aphid (Additional files 2 and 3). Fewer responses were noted in the 6-wk sorghum. Specifically, the resistant genotype up-regulated 2376 genes while the susceptible genotype up-regulated only 985 genes when aphids were feeding on the plants (Additional files 4 and 5). The overlaps between different sets of DEGs representing changes in the 2- and 6-wk susceptible and resistant genotypes in response to aphid herbivory were further illustrated in Venn diagrams (Fig. 2). In the 2-wk sorghum, 2430 DEGs were shared between both genotypes in response to aphid herbivory, while 3525 and 986 DEGs were only expressed in the resistant and susceptible genotypes, respectively. This result indicates that the resistant 2-wk old sorghum responds to the aphid feeding by eliciting a higher number of unique responses than the susceptible genotype. On the other hand, 869 DEGs were shared between the two genotypes exposed to the aphids at 6-wk post-emergence, and 3820 and 116 DEGs were uniquely expressed in resistant and susceptible genotype, respectively. The higher numbers of DEGs in the younger resistant genotype in response to M. sacchari indicate that these plants were undergoing significant metabolic changes. This outcome was not surprising and has been reported previously. For example, Lee et al.  found higher number of DEGs in resistant soybean in response to soybean aphid, Aphis glycines Matsumura (Hemiptera: Aphididae) in comparison to an aphid-susceptible soybean line.
Gene ontology (GO) analysis
Aphid herbivory affected a higher number of significant GO terms in the 2-wk plants than in the 6-wk plants regardless of genotype (Fig. 3). In the resistant 2-wk genotype, up-regulated genes associated with metabolic process, biosynthetic process, macromolecule process, and oxidation reduction were included in the enriched biological process GO terms. “Oxidoreductase activity” was also the prominent enriched GO term in the molecular function category in the 2-wk resistant genotype (Fig. 3a). Similarly, we noted enrichment in polysaccharide biosynthesis, carbohydrate biosynthesis, metabolic process, and oxidation reduction GO terms in up-regulated genes in the 2-wk susceptible genotype (Fig. 3b). The number of enriched GO terms in both up- and down-regulated genes was higher in the 6-wk resistant genotype (Fig. 3c) than in the 6-wk susceptible genotype (Fig. 3d). The majority of the down-regulated genes were enriched in several nucleotide binding GO terms, and it was noteworthy that “oxidoreductase activity” was the only GO term that was down-regulated in the 6-wk susceptible genotype (Fig. 3). On the other hand, genes involved in detoxification, which is one of the defenses linked to plant response to aphid herbivory [37, 45, 46] was up-regulated in both 2- and 6-wk resistant genotype and in the 2-wk susceptible genotype (Fig. 3a, b and c). The up-regulation of genes involved in oxidative-reduction process has been reported in several plants in response to herbivory, such as in barly in response to Diuraphis noxia Kurdjumov (Hemiptera: Aphididae) , in switchgrass exposed to greenbug , in soybean at both compatible and incompatible interactions with soybean aphid [37, 49], and in cotton in response to Bemisia tabaci Gennadius (Hemiptera: Aleyrodidae) . Our outcomes corroborate these studies and suggest that oxidative responses that are elicited in the resistant genotype are may be the key components of sorghum resistance to sugarcane aphids.
KEGG pathway analysis overview
Pathway analysis revealed significant effects of the genotype on a number of enriched pathways (Fig. 4). Specifically, 3270 genes that mapped to 24 pathways were up-regulated in the 2-wk resistant genotype in response to aphid herbivory, while 1275 up-regulated genes from 10 pathways were elicited in the 2-wk susceptible genotype (Fig. 4a, b). Fewer pathways were enriched in the 6-wk old sorghum.
In general, we found that a significantly higher number of pathways were up regulated rather than down-regulated in the 2-wk old plants regardless of the genotype (Fig. 4a, b). It is noteworthy that in our study plant hormone signal transduction and glutathione metabolism pathways, both shown previously to be highly relevant to plant defense [4, 50,51,52], were only enriched in the 2-wk resistant genotype along with a plant-pathogen interaction pathway (Fig. 4a). Induction of specific genes within these pathways suggests that their higher expression may be driving resistance to the aphids in this genotype, which appears the strongest in younger plants. In the susceptible 2-wk plants, biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, and biosynthesis of flavonoids associated with plant defense and stress responses were up-regulated in response to aphid herbivory (Fig. 4b).
Aphids induced expression of genes in defense pathways in the 6-wk resistant genotype as well. These pathways included biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, plant-pathogen interaction, and linoleic acid metabolism. However, the magnitude of these inductions was lower than in the younger sorghum (Fig. 4c). Biosynthesis of secondary metabolites and phenylpropanoid pathway were also enriched in the 6-wk susceptible sorghum (Fig. 4d), but the number of up-regulated genes was lower in these pathways compared with the 2-wk plants regardless of genotype. Tu et al.  reported that synthesis of linoleic acid in response to thrips, Odontothrips loti Haliday (Thysanoptera: Thripidae) increased only in the resistant alfalfa cultivar and not in the susceptible cultivar. Moreover, Ponzio et al.  also reported that aphid infestation in Brassica nigra caused the increasing of the linoleic acid. These secondary metabolites have been previously shown to be induced by herbivory in general [43, 55]. For example, Dubey et al.  reported the induction of transcripts related to secondary metabolite processes including phenylpropanoid and flavonoid biosynthesis in cotton exposed to A. gossypii. Further, Liang et al.  reported induction of genes involved in signal transduction, flavenoid metabolism, and plant-pathogen interactions that were likely invovled in resistance of cucumber to A. gossypii.
Specific genes and pathways relevant to sorghum resistance to M. sacchari
Based on the GO enrichment analysis and KEGG pathway enrichment analysis, we identified multiple defense-related genes that were differentially expressed in response to aphid herbivory. These aphid-responsive genes were grouped according to their function and discussed further. Specifically, we focused on signaling related genes, which act as mediators of herbivory-induced changes that occur during the early phases of the attack [4, 46, 57, 58] and phytohormone signaling pathways. We also elaborate on pathways and genes involved in secondary metabolism and stress response, and transcription factors (TFs), which are essential for the regulation of gene expression .
Signaling-related genes: Calcium sensors, ROS, and MAPKs
A number of genes from the Ca2+ sensors group were up-regulated in response to aphid herbivory regardless of genotype or sorghum age. These genes encode calmodulin-like (CMLs) and calmodulin (CAM) proteins that are important in herbivore defense pathways in plants [57, 60]. We noted several genes from this family to be affected by exposure to the aphids. For example, gene Sobic.001G393500 from the CAM gene family showed higher expression in response to aphid herbivory regardless of genotype and age, while another gene from this family, Sobic.008G159100 was up-regulated in the 6-wk sorghum (Additional files 2, 3, 4 and 5). Further, the majority of mitogen-activated protein kinases (MAPKs) genes including LPR receptor-like, and serine/threonine-protein kinase were up-regulated in response to aphid herbivory in the resistant genotype regardless of age. The activation of MAPKs, in turn, can phosphorylate transcriptional regulators and ultimately affect gene transcription . The serine/threonine-protein kinase OXI1 gene (Sobic.006G128200), involved in signaling pathway by producing reactive oxygen species (ROS), was up-regulated in response to aphid across sorghum genotypes and age. The OXI1 gene in Arabidopsis has been reported to be required for resistance to downy mildew infection . Moreover, genes from the protein phosphatase 2C family (PP2C) involved in phosphorylation were also up-regulated in response to the aphid herbivory. Plant PP2Cs have been reported to be relevant to plant stress signaling [63, 64].
Signaling-related genes: Plant hormone signal transduction
In hormone signal transduction pathways, the majority of DEGs (25 genes) were up-regulated in the 2-wk resistant genotype in response to aphid herbivory (Fig. 5). These genes belong to the SA, JA, abscisic acid (ABA), auxin, cytokinins (CKs), and ET signal transduction pathways. Among these phytohormones, salicylic acid plays an important role in mediating induced plant responses to pathogens in particular [40, 65, 66], but has been implicated in herbivore responses as well [67, 68]. In this pathway, Sobic.001G143000 gene [NONEXPRESSOR OF PATHOGENESIS-RELATED GENE 1 (NPR1)], had a 10.6- and 2.4-fold increase in expression in response to aphid herbivory in both 2- and 6-wk resistant genotype, respectively (Additional files 2 and 4). In the susceptible genotype, three NPR1 genes were up-regulated only in the 2-wk plants exposed to the aphids (Additional file 3). The NPR1 is a SA receptor protein and the transcriptional regulation of NPR1 is a key factor in regulation of defense signaling responses downstream of SA [69, 70].
Several DEGs in the JA signaling pathway that regulates direct and indirect plant responses against herbivores [71,72,73] were also induced in our experiment. These DEGs included genes from the jasmonate ZIM-domain (JAZ) family, which function as repressors of JA signal output and high levels of JAZ proteins control the JA responses by repressing their transcription factors [73, 74]. These genes had more than 2-fold change in response to aphid herbivory in the 2- and 6-wk resistant and susceptible genotype (Fig. 5). They are located up-stream of MYC2 transcription factor and act as the negative regulator of MYC2 function. The MYC2 transcription factor, on the other hand, increases JA-mediated resistance to insect pests . These interactions are examples of negative crosstalk between JA and SA [76, 77], and the outcomes of our experiments indicate an antagonistic interaction between these hormones in sorghum exposed to aphid herbivory. We noted the greatest expression of NPR1 genes along with up-regulation of JAZ genes that act as negative regulator of the JA-signaling pathway in both 2- and 6-wk resistant and in the 2-wk susceptible genotype. On the other hand, JAZ genes were upregulated while NPR1 gene expression was unaffected in the 6-wk susceptible genotype, which had the weakest overall transcriptional responses to aphid herbivory. As has been reported previously [40, 78], interactions between JA and SA are complex and can be host-  and herbivore-specific . Moreover, it has been suggested that insects can manipulate plant defenses by modulating the JA/SA crosstalk [81,82,83]. Therefore, we speculate that up-regulation of JAZ genes may be mediated by sugarcane aphids in order to suppress host defense strategies. This was previously demonstrated by Moran, Thompson  who reported that green peach aphid (Myzus persicae Sulzer, Hemiptera: Aphididae) and greenbug aphid feeding on Arabidopsis and sorghum respectively, up-regulated SA-dependent pathway and reduced activation of JA-dependent defenses in host plants.
Abscisic acid is a plant hormone involved primarily in the perception of abiotic stresses [85,86,87], which can be elicited by aphid herbivory [88, 89]. Four genes in the ABA signaling cascade were up-regulated in the resistant 2-wk sorghum including Sobic.006G279100, Sobic.001G350700 genes encoding SNF1-related protein kinase 2 (SnRK2, a positive regulator of ABA), and Sobic.003G198200, Sobic.001G424400 genes coding for type 2C protein phosphates (PP2C, a negative regulator of ABA) (Additional file 2). Five genes in this pathway were also up-regulated in the resistant 6-wk sorhgum, including bZIP23 transcription factor (Sobic.004G309600), which showed 3.2-fold induction in response to aphid herbivory (Additional file 4). Zong et al.  reported that OsbZIP23 acts as a central regulator in ABA signaling, biosynthesis, and drought resistance in rice. Therefore, the effects of aphid feeding on inducing of ABA signaling are likely related to induced water stress [91, 92].
Genes involved in auxin signaling pathway were also differentially expressed in both resistant and susceptible genotypes (Fig. 5). Auxin is involved in regulation of plant growth and development, and activation of auxin signaling is usually linked to increased disease susceptibility [93, 94]. The majority of these genes were up-regulated in 2-wk sorghum regardless of genotype (Fig. 5). However, there were key differences between the resistant and susceptible 2-wk old plants. Specifically, two members of AUX/IAA gene family i.e. Sobic.009G069700 and Sobic.002G055900 that act as a transcriptional repressor were up-regulated in the 2-wk resistant genotype in response to aphid herbivory, but only Sobic.002G055900 gene was up-regulated in the susceptible 2-wk genotype (Additional files 2 and 3). While these trends in expression are potentially important, it is not clear if these subtle transcriptional responses can explain why sugarcane aphid performance differs among genotypes.
The SA-regulated auxin homeostasis, in general, has interesting implications for the trade-off between growth and defense [95, 96] such that higher SA levels reduce the pool of active indole-3-acetic acid  and subsequently defense is prioritized over growth. On the other hand, in absence of pathogens, auxin-mediated suppression of SA responses can have positive consequences for plant growth . In our experiment, Sobic.002G361500 encoding GH3–8, an auxin-responsive gene involved in processes that halt accumulation of free IAA, showed a 6.6- and 3.9-fold increases in expression in response to aphid herbivory in 2-wk resistant and susceptible genotype, respectively (Additional files 2 and 3). This suggests that up-regulation of these genes may have triggered a change in allocation of resources from growth to defense and conferred greater resistance to the aphids in the younger plants. It has been also reported that auxin-regulated genes such as AUX/IAA and GH, are involved in an inhibitory feedback loop that maintains cellular auxin homeostasis that is required for adjusting auxin concentration for the optimal plant growth and development [99,100,101]. Induction of these genes in the 2-wk sorghum regardless of genotype but not in the 6-wk old sorghum suggests that the younger plants may be diverting their resources from primary to secondary metabolism thus increasing their resistance to the herbivore. However, more research is needed to unequivocally link these transcriptional responses to decreased aphid performance in the younger and resistant sorghum.
Cytokinins are plant hormones that act as key regulators of plant growth-defense trade-off and are also involved in plant-biotic interactions [98, 102]. We found a number of significantly induced genes in cytokinine signaling pathway, several of which were induced in the absence of aphids and suppressed when aphids were feeding on the plants. Specifically, in the 2-wk aphid-free resistant sorghum, three positive cytokinine regulator genes (Sobic.009G202900, Sobic.003G292600, and Sobic.004G330900) were up-regulated, indicating that plants were regulating their growth processes actively (Additional file 2). On the other hand, Sobic.006G090300 gene, a two-component response regulator ORR1 that acts as a negative regulator of cytokinine  was up-regulated in 2-wk resistant genotype exposed to aphid herbivory, suggesting that the younger resistant plants were suppressing their growth-related processes (Additional file 2). This outcome provides another example of the younger resistant sorghum switching from growth to defense in response to the aphids and highlights the lack of parallel responses in the susceptible sorghum. These results indicate a potential trade-off between growth and defense, which are both metabolically costly processes that usually cannot be carried out concurrently at high levels .
Ethylene is one of induced defense related phytohormones that is produced in response to multiple stresses [103, 104]. Ethylene expression decreases with plant age, with a marked decline at the onset of reproduction associated with a significant reduction in the ability of plants to induce defenses . This was consistent with our results – we found ethylene gene expression only in the younger plants. In ethylene signaling pathway, two genes i.e. Sobic.007G210700 and Sobic.009G050400, involved in the ethylene pathway were both up-regulated in the resistant 2-wk genotype in response to aphid and only Sobic.007G210700 was up-regulated in the susceptible 2-wk genotype (Additional files 2 and 3). The release of ethylene in response to aphid herbivory has been reported in other host plant species [2, 106, 107].
Secondary metabolites related to oxidative stress response genes
A number of DEGs affected by aphid herbivory in our study are involved in metabolite synthesis through phenylpropanoid and flavonoid biosynthesis pathways. In addition, several other genes that were induced regulate plant responses to oxidative stress caused by herbivory. A number of genes in phenylpropanoid biosynthesis pathway showed differential expression in both susceptible and resistant genotypes exposed to the aphids (Fig. 6). Specifically, two genes Sobic.004G220400 and Sobic.004220700, both encoding phenylalanine ammonia-lyase (PAL) showed more than 2-fold increase in response to aphid herbivory in both the 2- and 6-wk resistant and susceptible genotypes (Additional files 2, 3, 4, 5). Phenylalanine ammonia-lyase catalyzes the first step of the phenylpropanoid pathway, which is a key reaction in the control of lignin, flavonoid and salicylic acid biosynthesis . Lv et al.  reported that expression of PAL was significantly increased in cotton and corn seedlings damaged by mechanical wounding or by cotton aphid and corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae) herbivory. Further, Chaman et al.  reported a positive relationship between SA concentration, PAL activity, and resistance of barley to aphid infestations. Therefore, induction of these genes across sorghum genotypes and age in our experiment indicates that they play a key role in sorghum response to aphid herbivory regardless of sorghum age or genotype.
Genes that were involved in response to oxidative stress belong to glutathione S-transferase, L-ascorbate peroxidase, peroxidase, and cytochrome P450 gene family. The up-regulation of a number of genes involved in glutathione metabolism pathway was mostly found in the resistant 2-wk sorghum in response to aphid herbivory (Additional file 2). It has been proposed that the activation of glutathione metabolism is a natural cell response to stress that can be induced by aphids in general [56, 110], and the increase in expression of enzymes involved in glutathione synthesis is correlated with resistance to various biotic challenges  and in the detoxification of ROS . Further, among the DEGs involved in phenylpropanoid pathway in our experiment, a higher number of genes encoding peroxidase showed differential expression in response to aphid herbivory in younger plants in comparison with older ones regardless of genotype. Peroxidases not only code for ROS-detoxifying enzymes but are also involved in oxidative signal transduction, regulating the redox and Ca2+ homeostasis as well as the expression of defense genes [113,114,115]. Higher expression of peroxidases has been reported previously in barley exposed to Diuraphis noxia (Hemiptera: Aphididae) , and in tomato fed upon by Macrosyphum euphorbiae (Hemiptera: Aphididae) . Likewise, the cytochrome P450 superfamily genes were another abundant gene family that was induced in response to oxidative stress imposed by the aphids in our study. This gene family plays a crucial role in plant defense via multiple biosynthetic and detoxification pathways . Prochaska et al. , for example, reported that genes from cytochrome P450 family were up-regulated in the resistant genotype of soybean in response to aphid herbivory. Notably, we found these genes to be induced the strongest in the 2-wk resistant sorghum exposed to sugarcane aphids (Additional file 2), which may explain the superior host plant resistance observed in this genotype shortly after emergence.
DEGs encoding transcription factors
Our results indicated differential expression of 71 genes encoding transcription factors (TFs) including WRKY, MYB, bHLH, NAC, and AP2/ERF. The WRKY family of TFs is unique to plants and is involved in regulation of defense responses to pathogens and insects [117, 118]. The involvement of selected WRKY family in response to biotic stresses has been reported in Arabidopsis [119, 120], and we also tracked the expression of the orthologs of these genes in sorghum. Specifically, we found the highest number of differentially expressed TFs (52) in the resistant genotype and in the younger plants, and 18 of these TFs belonged to the WRKY family. In the 6-wk resistant genotype, a total of 42 TFs were up regulated in response to aphid herbivory, and 13 of these were WRKY TFs. A comparable number of WRKY genes were up-regulated in response to aphid herbivory in the 2- and 6-wk susceptible plants (Fig. 7). It is noteworthy that four WRKY transcription factors were found to be expressed exclusively in the 2-wk resistant genotype exposed to the aphids, including WRKY1 (Sobic.001G095500), WRKY19 (Sobic.009G238200), WRKY28 (Sobic.003G199400) and WRKY72 (Sobic.005G117400) (Additional file 2). The involvement of the homologs of these genes in hormone signaling pathways and response to abiotic stresses has been reported in Arabidopsis [42, 78, 119, 121]. For example, the homolog of WRKY19 in Arabidopsis, i.e., AT4G11070 has been shown to be a key regulator of cross-talk between SA and JA signaling pathways . Our outcomes corroborate these reports – the sorghum orthologs of AtWRKY70 i.e. Sobic.003G337900, and Sobic.008G060300 were only expressed in response to aphid herbivory in 2-wk plants regardless of genotype, with higher expression in the resistant genotype. Similarly, the AtWRKY40 ortholog in sorghum, i.e. Sobic.004G065900, was up-regulated in response to aphid herbivory regardless of genotype and age. The up-regulation of AtWRKY40 has been reported in Arabidopsis in response to cabbage aphid, Brevicoryne brassica (Hemiptera: Aphididae) . Three TFs with AP2 domain i.e., Sobic.002G415600, Sobic.003G380900, and Sobic.004G310600, prominent TF regulators of host defense [123, 124], were up-regulated in response to aphid herbivory across genotypes and sorghum age. Other types of TFs identified among the DEGs were bZIP, bHLH, MYB, NAC, zinc finger, several of which have been reported to be important for plant adaptive response to biotic stress [125, 126].
Real-time quantitative RT-PCR (RT-qPCR) of selected DEGs
All genes selected for validation of RNA-seq data showed a similar expression trend (up-regulation or down-regulation) in the RT-qPCR analysis. High correlations (all above R2 > 0.70) were found between RNA-seq and RT-qPCR results (Additional file 6), indicating that the measured changes in gene expression detected by RNA-seq reflected the actual transcriptome differences between the samples.
Aphid population growth on resistant and susceptible sorghum 2- and 6-wk post-emergence
Sorghum variety and growth stage had a significant effect on the mean daily change in M. sacchari densities (Fig. 8). The average change in aphid densities in susceptible sorghum was significantly greater than in resistant sorghum (F1,32 = 92.87; p-value < 0.001), and aphid populations increased faster in the 6-wk old sorghum than in the 2-wk old sorghum (F1,32 = 15.01; p-value < 0.001). However, there was no interactive effect of these factors on the average change in aphid densities (F1,32 = 1.36; p-value = 0.253). Aphid populations grew the slowest in the young resistant sorghum, and their numbers increased the most rapidly in the older susceptible sorghum. Rapid population growth of M. sacchari in susceptible sorghum in reproductive stages has been reported previously , and our outcomes corroborate those findings. Further, these results indicate that the robust transcriptional response to aphid herbivory in the younger resistant genotype is phenotypically reflected in the significantly diminished aphid performance on these plants.
The host plant resistance reported previously for the resistant genotype and manifested in significantly lower aphid densities in the greenhouse and field experiments  was reflected at the transcriptome level in our experiment. Specifically, transcriptome analysis of the susceptible and resistant sorghum genotypes revealed a suite of genes and pathways elicited in response to M. sacchari herbivory, particularly in the resistant and younger sorghum. These differentially expressed genes were involved in diverse processes including signaling, metabolism, and oxidative stress. We detected a higher number of DEGs and expression of a higher number of defense related genes in the resistant genotype in response to aphid herbivory than in the susceptible genotype. Further, the 2-wk old resistant genotype, in particular, mounted the most robust response to M. sachari feeding, and several genes and pathways provided insight into the mechanisms of aphid resistance.
The enriched GO terms and pathways highlighted the general differences between the two genotypes in response to aphid herbivory. The oxidation reduction GO term was highly enriched in the resistant genotype compared to the susceptible genotype. We also found that pathways coding for secondary metabolites, plant hormone signal transduction, plant-pathogen interaction, and glutathione metabolism were over-expressed in the resistant genotype. The resistant genotype also activated more genes involved in SA, JA, ET, and ABA signaling pathways. Our results suggested a potential cross-talk between SA and JA in the resistant genotype such that up-regulation of SA pathway genes negatively affected JA pathways genes through higher expression of JAZ genes that act as a negative regulator of JA biosynthesis pathway genes. Aphid herbivory also affected the balance between growth and defense related genes to the benefit of defense. Similar differences were noted when gene expression was compared between the 2-wk and 6-wk old sorghum – the number of DEGs, enriched GO terms, and significantly affected pathways were higher in the 2-wk old plants than in the 6-wk old sorghum. On the other hand, the lower resistance of older plants may be related to the developmental stage of these plants. Sorghum at 6-wk after planting are usually at panicle initiation stage, and are especially sensitive to biotic and abiotic stress during this phase .
In conclusion, sugarcane aphid herbivory elicited the most robust transcriptional changes in the younger resistant genotype, and several pathways and specific genes provide insights into the mechanisms underlying host plant resistance to this invasive insect. Despite the apparent cross talk between SA and JA elicited by aphid herbivory, the younger resistant plants regulated their responses to abiotic stress, diverted resources to defense over growth, and induced strong responses to oxidative stress at the highest level compared to all other treatments. Lack of parallel responses in the older sorghum plants has relevant implications for aphid-sorghum interactions, and suggests lower host plant resistance in sorghum colonized at the onset of the reproductive stage. Further, unique induction of four WRKY genes only in the 2-wk resistant genotype highlights the likely importance of these TFs in the superior resistance of these plants to aphids, and should be considered in future studies. Further research should explore functional role of genes we found to be induced by M. sacchari to further clarify their impact on aphid-plant interactions.
Sorghum growth and M. sacchari colony
All experiments were conducted at the Texas A&M AgriLife Research greenhouse complex located at the Plant Stress Laboratory in Bushland, TX. Two sorghum genotypes, aphid-susceptible (DKS 44–20, Dekalb ®, Dekalb, IL), and aphid-resistant (DKS 37–07, Dekalb ®, Dekalb, IL) at 2- and 6-wk post-emergence grown under a long-day photoperiod regime were used in this study. The two genotypes were grown under controlled greenhouse conditions with a temperature range of 21 ± 3 °C  to 30 ± 3 °C . Plants were grown in pots containing Redi-earth plug and seedling soil mix (Sun Grow Horticulture Inc., Agawam, MA) in 60 × 60 × 60 cm insect-proof rearing tents (MegaView Science Education Services Co., Taipei, Taiwan). Plants assigned to the 6-wk treatment were grown in 1 × 1 × 1 m cages constructed from PVC pipes and insect-proof mesh.
Sugarcane aphid colony originated from aphids collected from sorghum growing on a research farm affiliated with Texas A&M AgriLife Research in Bushland, TX, and no permits were required to enter the land or collect the aphids. The colony was maintained in the greenhouse for one year prior to the experiment. In order to confirm that the colony was genetically homogeneous, universal DNA primers were used to amplify a 710 bp fragment of the mitochondrial cytochrome c oxidase subunit I gene (COI) following Vrijenhoek  methods. Two samples from two cages harboring the colony were collected for this analysis. The PCR products were purified using QIAGEN MinElute PCR purification kit (QIAGEN, Valencia, CA, USA) and sequenced on ABI 3130xl genetic analyzer by the Institute for Plant Genomics and Biotechnology at Texas A&M University. The sequencing data from four samples were compared with each other using EMBOSS Needle [pairwise sequence alignment (Nucleotide)] (http://www.ebi.ac.uk/Tools/psa/emboss_needle/nucleotide.html). There was no variation within and between a population that confirmed the integrity of the aphid colony. Further, the COI sequencing confirmed that the aphids were of different haplotype than previously described in sugarcane in the U.S. .
Sample collection, RNA isolation, and sequencing
The experiment was a 2 × 2 × 2 factorial design with two levels of the sorghum genotype factor (susceptible, resistant), two levels of the sorghum age (2-wk, 6-wk), and two levels of aphid herbivory (present, absent). An individual pot was considered a biological replicate and three biological replications and three technical replications (i.e. individual plant in each pot) were included in this experiment. Half of the susceptible and half of the resistant plants were randomly assigned to one of the two sorghum age treatments and the herbivory treatment, which consisted of exposing two leaves in each plant to 15 aphids for 24 h. The entire above-ground tissue of each 2-wk plant was collected for RNA-seq. Six-wk old plants were exposed to aphids in the same manner, with the exception of tissue collection for RNA extraction – only leaves exposed to the aphids and not entire plants were collected. Owing to the extensive size of the 6-wk old sorghum that surpassed 50 cm in height and had more than 10 fully unfurled leaves only the tissue that was directly exposed to aphid herbivory was collected and used in RNA extractions. Aphids on the 2-wk old plants, on the other hand, were not confined to the single leaf to which they were introduced, and hence the entire plant (10–15 cm in height with three leaves) was collected. Samples were immediately frozen in liquid nitrogen (LN2) and stored at − 80 °C until further sample processing.
Frozen leaf samples were ground under LN2 into a fine powder using a mortar and pestle. RNA was extracted from 100 mg of frozen ground tissue using the miRNeasy Mini Kit (Qiagen, Valencia, CA, USA) and subsequently treated with DNase Max Kit (Qiagen, Valencia, CA, USA). RNA quality and quantity were assessed using the NanoVue Plus spectrophotometer (GE, Healthcare, Piscataway, NJ, USA). Three equimolar RNA samples extracted from three technical replicates for each biological replicate were pooled after extraction for RNA-seq template preparation. The quality of each pooled RNA sample was assessed with an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) prior to RNA-seq template preparation by Texas AgriLife Research Genomic and Bioinformatics Services. Paired-end (PE) sequencing was performed on an Illumina HiSeq4000 using a 75 bp paired-end strategy on three biological replicates. RNA-seq template preparation and sequencing was completed at Texas AgriLife Research Genomics and Bioinformatics Center. The quality of RNA-seq reads was assessed using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/).
Gene expression analysis
Sequence reads were imported into the CLC Genomics Workbench version 11 (Qiagen, Valencia, CA, USA) and mapped to the S. bicolor reference genome [Sbicolor_454 v3.1.1, www.phytozome.jgi.doe.gov, . Based on the total read counts for each annotated gene, differential gene expression analyses were conducted using the Empirical Analysis of DGE tool, which implements the ‘Exact Test’ for two-group comparisons . The transcriptional response of sorghum plants at 2- and 6-wk post-emergence and exposed to aphid herbivory was compared to aphid-free plants (the control) to identify DEGs. Differentially expressed genes were defined as having a fold change ≥1.9 or ≤ − 1.9 with a false discovery rate (FDR) corrected p-value < 0.05. For functional annotation, GO analysis was performed using AgriGO gene ontology analysis tools  to determine overrepresented GO categories in the up- and down-regulated DEGs, and significantly enriched GO terms were identified. Pathway analysis was performed using the KEGG database  to underline the pathways to which the up- and down-regulated DEGs contribute, and pathway enrichment analysis completed using the KOBAS server (version v.3) . The overlaps between different sets of DEGs were generated with VENNY.2.1 (http://bioinfogp.cnb.csic.es/tools/venny/).
RT-qPCR analysis of selected DEGs
Five DEGs involved in plant-hormone signal transduction pathway were selected for validation of RNA-seq data using RT-qPCR. The same RNA samples from three biological replicates that were used for the sequencing were used for the validation of RNA-seq data. The total RNA was reverse-transcribed using Omniscript RT Kit (QIAGEN, Valencia, CA, USA). Primer pairs were designed using NCBI Primer-BLAST (http://www.ncbi.nlm. nih.gov/tools/primer-blast/), which are listed in Additional file 7. The RT-qPCR reaction was performed in a total volume of 20 μl, containing 1 μl of diluted cDNA, 0.5 μl of reverse and forward primers, 8 μl of ddH2O, 10 μl of SYBR Green Master Mix (Applied Biosystems, Foster City, CA, USA). Samples were run on a Applied Biosystems ViiA 7 Real-Time PCR system (Applied Biosystems, Foster City, CA, USA) according to the standard protocol as follows: 2 min at 50 °C, 10 min at 95 °C, followed by 40 cycles of 15 s at 95 °C for and 1 min at 60 °C. A melt curve was generated by the end of each PCR reaction to verify the formation of single peak and to exclude the possibility of primer dimer and non-specific product formation. All reactions were performed in duplicate, including the non-template control reactions. The relative quantitation of gene expression of RT-qPCR was measured via the 2−ΔΔCt method , with CYP (Cyclophilin/Peptidylprolyl Isomerase) gene as the endogenous reference gene that showed the highest stability among three tested reference genes. The correlation between the RNA-seq data and the RT-qPCR results was determined by Pearson’s correlation coefficient.
Aphid population growth evaluation
In another set of experiments, differences in aphid densities on 2- and 6-wk susceptible and resistant sorghum genotypes was evaluated in early spring 2017. Experiments were conducted in the greenhouse following the protocol described above (Sorghum growth and M. sacchari colony). The experimental design was a randomized complete design with two levels of sorghum growth stage (2- and 6-wk old sorghum), and two levels of the sorghum variety factor (Susceptible – DKS 44–20, and Resistant – DKS 37–07). The entire experiment was conducted twice (time block) and each of the treatment combinations was replicated a total of nine times across the two time blocks (n = 36). When sorghum reached the assigned growth stage, 15 sugarcane aphids were moved using a paintbrush and clippings of infested sorghum onto two leaves on each of the plants. The aphids were counted and replenished as needed 1 h following the introductions, and aphid densities following a 24 h acclimation period were used as the initial densities. Aphids were counted on entire plants every three days thereafter for 14 d. The mean daily change in aphid densities per plant was the response variable recorded and analyzed.
Statistical analyses for greenhouse evaluation
The effects of sorghum age and variety on the average daily change in aphid densities was analyzed using two-way ANOVA with sorghum age and variety combination as a fixed factor and time block as a random factor. Owing to lack of a significant interaction between treatments and time block data from the two greenhouse experiments were combined. Effects of sorghum variety, age, and the interaction between the two factors were included in the analyses, and means separation tests were performed were appropriate with a Bonferroni adjustment. All analyses were performed in R .
Walling LL. The myriad plant responses to herbivores. J Plant Growth Regul. 2000;19(2):195–216.
Howe GA, Jander G. Plant immunity to insect herbivores. Annu Rev Plant Biol. 2008;59:41–66.
Koch KG, Chapman K, Louis J, Heng-Moss T, Sarath G. Plant tolerance: a unique approach to control hemipteran pests. Front Plant Sci. 2016. doi:org/https://doi.org/10.3389/fpls.2016.01363.
Wu J, Baldwin IT. New insights into plant responses to the attack from insect herbivores. Annu Rev Genet. 2010;44:1–24.
Foyer CH, Rasool B, Davey JW, Hancock RD. Cross-tolerance to biotic and abiotic stresses in plants: a focus on resistance to aphid infestation. J Exp Bot. 2016;67(7):2025–37.
Goggin FL, Zhu-Salzman K. Editorial overview: pests and resistance: social networking—studying the web of plant–insect interactions to improve host plant resistance. Curr Opin Insect Sci. 2015;9:v–viii.
Hawley CJ, Peairs FB, Randolph TL. Categories of resistance at different growth stages in halt, a winter wheat resistant to the Russian wheat aphid (Homoptera: Aphididae). J Econ Entomol. 2003;96(1):214–9.
Smith CM, Boyko EV. The molecular bases of plant resistance and defense responses to aphid feeding: current status. Entomol Exp Appl. 2007;122(1):1–16.
Hill CB, Li Y, Hartman GL. A single dominant gene for resistance to the soybean aphid in the soybean cultivar Dowling. Crop Sci. 2006;46(4):1601–5.
Omoigui L, Ekeuro G, Kamara A, Bello L, Timko M, Ogunwolu G. New sources of aphids [Aphis craccivora (Koch)] resistance in cowpea germplasm using phenotypic and molecular marker approaches. Euphytica. 2017;213(8):178–92.
Sharma HC, Bhagwat VR, Daware DG, Pawar DB, Munghate RS, Sharma SP, et al. Identification of sorghum genotypes with resistance to the sugarcane aphid Melanaphis sacchari under natural and artificial infestation. Plant Breed. 2014;133(1):36–44.
Reddy PS, Bhagwat VR, Prasad GS, Tonapi VA. Breeding for insect resistance in sorghum and millets. In: Arora R, Sandhu S, editors. Breeding insect resistant crops for sustainable agriculture. Singapore: Springer; 2017. p. 231–64.
Bayoumy MH, Perumal R, Michaud J. Comparative life histories of greenbugs and sugarcane aphids (Hemiptera: Aphididae) coinfesting susceptible and resistant sorghums. J Econ Entomol. 2015;109(1):385–91.
Reese JC, Schwenke JR, Lamont PS, Zehr DD. Importance and quantification of plant tolerance in crop pest management programs for aphids: greenbug resistance in sorghum. J Agric Entomol. 1994;11(3):255–70.
Johnson J, Rosenow D, Teetes G, Phillips J. Registration of 19 greenbug resistant sorghum germplasm lines1 (Reg. Nos. GP 79 to GP 97). Crop Sci. 1982;22(6):1272.
Kofoid K, Harvey T. Registration of greenbug resistant sorghum germplasm lines KS 116 a/B through KS 120 a/B. Crop Sci. 2005;45(2):802–4.
Punnuri S, Huang Y. Identification and confirmation of greenbug resistance loci in an advanced mapping population of sorghum. J Agric Sci. 2017;155(10):1610–22.
Huang Y. Improvement of crop protection against greenbug using the worldwide sorghum germplasm collection and genomics-based approaches. Plant Genet Resour. 2011;9(2):317–20.
Zapata SD, Villanueva R, Sekula D, Esparza-Diaz G, Duke K, Mutaleb M, editors. The economic impact of the sugarcane aphid on sorghum production. Southern Agricultural Economics Association Annual Meeting. San Antonio: Southern Agricultural Economics Association; 2016.
Villanueva R, Brewer M, Way M, Biles S, Sekula D, Bynum E, et al. Sugarcane aphid: a new pest of sorghum: Texas a&M AgriLife extension: Ento-035; 2014.http://denton.agrilife.org/files/2013/08/ENTO-035-The-Sugarcane-Aphid-2014.pdf
Nibouche S, Costet L, Holt JR, Jacobson A, Pekarcik A, Sadeyen J, et al. Invasion of sorghum in the Americas by a new sugarcane aphid (Melanaphis sacchari) superclone. PLoS One. 2018;13(4):e0196124.
Bowling RD, Brewer MJ, Kerns DL, Gordy J, Seiter N, Elliott NE, et al. Sugarcane aphid (Hemiptera: Aphididae): a new pest on sorghum in North America. J Integr Pest Manag. 2016;7(1):1–13.
Brewer MJ, Gordy JW, Kerns DL, Woolley JB, Rooney WL, Bowling RD. Sugarcane aphid population growth, plant injury, and natural enemies on selected grain sorghum hybrids in Texas and Louisiana. J Econ Entomol. 2017;110(5):2109–18.
Tsutsui ND, Suarez AV, Holway DA, Case TJ. Reduced genetic variation and the success of an invasive species. Proc Natl Acad Sci. 2000;97(11):5948–53.
Bossdorf O, Auge H, Lafuma L, Rogers WE, Siemann E, Prati D. Phenotypic and genetic differentiation between native and introduced plant populations. Oecologia. 2005;144(1):1–11.
Grapputo A, Boman S, Lindstroem L, Lyytinen A, Mappes J. The voyage of an invasive species across continents: genetic diversity of north American and European Colorado potato beetle populations. Mol Ecol. 2005;14(14):4207–19.
Nibouche S, Fartek B, Mississipi S, Delatte H, Reynaud B, Costet L. Low genetic diversity in Melanaphis sacchari aphid populations at the worldwide scale. PLoS One. 2014;9(8):e106067.
Medina RF, Armstrong SJ, Harrison K. Genetic population structure of sugarcane aphid, Melanaphis sacchari, in sorghum, sugarcane, and Johnsongrass in the continental USA. Entomol Exp Appl. 2017;162(3):358–65.
Harris-Shultz K, Ni X, Wadl PA, Wang X, Wang H, Huang F, et al. Microsatellite markers reveal a predominant sugarcane aphid (Homoptera: Aphididae) clone is found on sorghum in seven states and one territory of the USA. Crop Sci. 2017;57(4):2064–72.
Smith CM. Plant resistance in arthropod pest management systems. In: Smith CM, editor. Plant resistance to arthropods: molecular and conventional approaches. The Netherlands: Springer Science & Business Media; 2005. p. 381–412.
Smith CM, Chuang WP. Plant resistance to aphid feeding: behavioral, physiological, genetic and molecular cues regulate aphid host selection and feeding. Pest Manag Sci. 2014;70(4):528–40.
Szczepaniec A. Interactive effects of crop variety, insecticide seed treatment, and planting date on population dynamics of sugarcane aphid (Melanaphis sacchari) and their predators in late-colonized sorghum. Crop Prot. 2018;109:72–9.
Mbulwe L, Peterson GC, Scott-Armstrong J, Rooney WL. Registration of sorghum germplasm Tx3408 and Tx3409 with tolerance to sugarcane aphid [Melanaphis sacchari (Zehntner)]. J Plant Regist. 2016;10(1):51–6.
Züst T, Agrawal AA. Mechanisms and evolution of plant resistance to aphids. Nat Plants. 2016. https://doi.org/10.1038/NPLANTS.2015.206.
Dogimont C, Bendahmane A, Chovelon V, Boissot N. Host plant resistance to aphids in cultivated crops: genetic and molecular bases, and interactions with aphid populations. C R Biol. 2010;333(6–7):566–73.
Boissot N, Thomas S, Sauvion N, Marchal C, Pavis C, Dogimont C. Mapping and validation of QTLs for resistance to aphids and whiteflies in melon. Theor Appl Genet. 2010;121(1):9–20.
Lee S, Cassone BJ, Wijeratne A, Jun T-H, Michel AP, Mian MR. Transcriptomic dynamics in soybean near-isogenic lines differing in alleles for an aphid resistance gene, following infestation by soybean aphid biotype 2. BMC Genomics. 2017;18(1):472–83.
Liang D, Liu M, Hu Q, He M, Qi X, Xu Q, et al. Identification of differentially expressed genes related to aphid resistance in cucumber (Cucumis sativus L.). Sci Rep. 2015. https://doi.org/10.1038/srep09645.
Prochaska TJ, Donze-Reiner T, Marchi-Werle L, Palmer N, Hunt TE, Sarath G, et al. Transcriptional responses of tolerant and susceptible soybeans to soybean aphid (Aphis glycines Matsumura) herbivory. Arthropod Plant Interact. 2015;9(4):347–59.
Pieterse CM, Van der Does D, Zamioudis C, Leon-Reyes A, Van Wees SC. Hormonal modulation of plant immunity. Annu Rev Cell Dev Biol. 2012;28(1):489–521.
Erb M, Meldau S, Howe GA. Role of phytohormones in insect-specific plant reactions. Trends Plant Sci. 2012;17(5):250–9.
Bari R, Jones JD. Role of plant hormones in plant defence responses. Plant Mol Biol. 2009;69(4):473–88.
Bennett RN, Wallsgrove RM. Secondary metabolites in plant defence mechanisms. New Phytol. 1994;127(4):617–33.
Kessler A, Baldwin IT. Plant responses to insect herbivory: the emerging molecular analysis. Annu Rev Plant Biol. 2002;53(1):299–328.
Toan TN, Thanh Huyen TTV, Chung MV. Antioxidative defense response to aphid-induced oxidative stress in Glycine max (L.) Merr. cv.“Nam Dan”. Vietnam J Sci Tech. 2016;54(6):719–28.
Li J, Zhu L, Hull JJ, Liang S, Daniell H, Jin S, et al. Transcriptome analysis reveals a comprehensive insect resistance response mechanism in cotton to infestation by the phloem feeding insect Bemisia tabaci (whitefly). Plant Biotechnol J. 2016;14(10):1956–75.
Gutsche A, Heng-Moss T, Sarath G, Twigg P, Xia Y, Lu G, et al. Gene expression profiling of tolerant barley in response to Diuraphis noxia (Hemiptera: Aphididae) feeding. Bull Entomol Res. 2009;99(2):163–73.
Donze-Reiner T, Palmer NA, Scully ED, Prochaska TJ, Koch KG, Heng-Moss T, et al. Transcriptional analysis of defense mechanisms in upland tetraploid switchgrass to greenbugs. BMC Plant Biol. 2017;17(1):46–66.
Bansal R, Mian M, Mittapalli O, Michel AP. RNA-Seq reveals a xenobiotic stress response in the soybean aphid, Aphis glycines, when fed aphid-resistant soybean. BMC Genomics. 2014;15(1):972–85.
Bidart-Bouzat MG, Imeh-Nathaniel A. Global change effects on plant chemical defenses against insect herbivores. J Integr Plant Biol. 2008;50(11):1339–54.
Patra B, Schluttenhofer C, Wu Y, Pattanaik S, Yuan L. Transcriptional regulation of secondary metabolite biosynthesis in plants. Biochim Biophys Acta. 2013;1829(11):1236–47.
Vannozzi A, Dry IB, Fasoli M, Zenoni S, Lucchin M. Genome-wide analysis of the grapevine stilbene synthase multigenic family: genomic organization and expression profiles upon biotic and abiotic stresses. BMC Plant Biol. 2012;12(1):130–51.
Tu X, Liu Z, Zhang Z. Comparative transcriptomic analysis of resistant and susceptible alfalfa cultivars (Medicago sativa L.) after thrips infestation. BMC Genomics. 2018. https://doi.org/10.1186/s12864-018-4495-2.
Ponzio C, Papazian S, Albrectsen BR, Dicke M, Gols R. Dual herbivore attack and herbivore density affect metabolic profiles of Brassica nigra leaves. Plant Cell Environ. 2017;40(8):1356–67.
Wink M. Plant breeding: importance of plant secondary metabolites for protection against pathogens and herbivores. Theor Appl Genet. 1988;75(2):225–33.
Dubey NK, Goel R, Ranjan A, Idris A, Singh SK, Bag SK, et al. Comparative transcriptome analysis of Gossypium hirsutum L. in response to sap sucking insects: aphid and whitefly. BMC Genomics. 2013;14(1):241–60.
Scholz S, Heyer M, Vadassery J, Mithöfer A. A role for calmodulin-like proteins in herbivore defense path-ways in plants. Endocytobiosis Cell Res. 2016;27(1):1–12.
Coppola V, Coppola M, Rocco M, Digilio MC, D’Ambrosio C, Renzone G, 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–32.
Yu H, Luscombe NM, Qian J, Gerstein M. Genomic analysis of gene expression relationships in transcriptional regulatory networks. Trends Genet. 2003;19(8):422–7.
Ranty B, Aldon D, Cotelle V, Galaud J-P, Thuleau P, Mazars C. Calcium sensors as key hubs in plant responses to biotic and abiotic stresses. Front Plant Sci. 2016. doi:org/https://doi.org/10.3389/fpls.2016.00327.
Tena G, Boudsocq M, Sheen J. Protein kinase signaling networks in plant innate immunity. Curr Opin Plant Biol. 2011;14(5):519–29.
Rentel MC, Lecourieux D, Ouaked F, Usher SL, Petersen L, Okamoto H, et al. OXI1 kinase is necessary for oxidative burst-mediated signalling in Arabidopsis. Nature. 2004;427(6977):858–561.
Fuchs S, Grill E, Meskiene I, Schweighofer A. Type 2C protein phosphatases in plants. FEBS J. 2013;280(2):681–93.
Moorhead GB, Trinkle-Mulcahy L, Ulke-Lemée A. Emerging roles of nuclear protein phosphatases. Nat Rev Mol Cell Biol. 2007;8(3):234–44.
Schweiger R, Heise AM, Persicke M, Müller C. Interactions between the jasmonic and salicylic acid pathway modulate the plant metabolome and affect herbivores of different feeding types. Plant Cell Environ. 2014;37(7):1574–85.
Vos I, Pieterse C, Wees S. Costs and benefits of hormone-regulated plant defences. Plant Pathol. 2013;62(S1):43–55.
Li Q, Xie Q-G, Smith-Becker J, Navarre DA, Kaloshian I. Mi-1-mediated aphid resistance involves salicylic acid and mitogen-activated protein kinase signaling cascades. Mol Plant-Microbe Interact. 2006;19(6):655–64.
Maffei ME, Mithöfer A, Boland W. Before gene expression: early events in plant–insect interaction. Trends Plant Sci. 2007;12(7):310–6.
Wu Y, Zhang D, Chu JY, Boyle P, Wang Y, Brindle ID, et al. The Arabidopsis NPR1 protein is a receptor for the plant defense hormone salicylic acid. Cell Rep. 2012;1(6):639–47.
Medina-Puche L, Castelló MJ, Canet JV, Lamilla J, Colombo ML, Tornero P. β-Carbonic anhydrases play a role in salicylic acid perception in Arabidopsis. PLoS One. 2017;12(7):e0181820.
Dar TA, Uddin M, Khan MMA, Hakeem K, Jaleel H. Jasmonates counter plant stress: a review. Environ Exp Bot. 2015;115:49–57.
Koo AJ. Metabolism of the plant hormone jasmonate: a sentinel for tissue damage and master regulator of stress response. Phytochem Rev. 2017;17(1):1–30.
Lortzing T, Steppuhn A. Jasmonate signalling in plants shapes plant–insect interaction ecology. Curr Opin Insect Sci. 2016;14:32–9.
Kazan K, Manners JM. MYC2: the master in action. Mol Plant. 2013;6(3):686–703.
Dombrecht B, Xue GP, Sprague SJ, Kirkegaard JA, Ross JJ, Reid JB, et al. MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis. Plant Cell. 2007;19(7):2225–45.
Wei J, van Loon JJ, Gols R, Menzel TR, Li N, Kang L, et al. Reciprocal crosstalk between jasmonate and salicylate defence-signalling pathways modulates plant volatile emission and herbivore host-selection behaviour. J Exp Bot. 2014;65(12):3289–98.
Thaler JS, Humphrey PT, Whiteman NK. Evolution of jasmonate and salicylate signal crosstalk. Trends Plant Sci. 2012;17(5):260–70.
Caarls L, Pieterse CM, Van Wees S. How salicylic acid takes transcriptional control over jasmonic acid signaling. Front Plant Sci. 2015. doi:org/https://doi.org/10.3389/fpls.2015.00170.
Szczepaniec A, Raupp MJ, Parker RD, Kerns D, Eubanks MD. Neonicotinoid insecticides alter induced defenses and increase susceptibility to spider mites in distantly related crop plants. PLoS One. 2013;8(5):e62620.
De Vos M, Van Oosten VR, Van Poecke RM, Van Pelt JA, Pozo MJ, Mueller MJ, et al. Signal signature and transcriptome changes of Arabidopsis during pathogen and insect attack. Mol Plant-Microbe Interact. 2005;18(9):923–37.
Walling LL. Avoiding effective defenses: strategies employed by phloem-feeding insects. Plant Physiol. 2008;146(3):859–66.
Zhang P, Zhu X, Huang F, Liu Y, Zhang J, Lu Y, et al. Suppression of jasmonic acid-dependent defense in cotton plant by the mealybug Phenacoccus solenopsis. PLoS One. 2011;6(7):e22378.
Zarate SI, Kempema LA, Walling LL. Silverleaf whitefly induces salicylic acid defenses and suppresses effectual jasmonic acid defenses. Plant Physiol. 2007;143(2):866–75.
Moran PJ, Thompson GA. Molecular responses to aphid feeding in Arabidopsis in relation to plant defense pathways. Plant Physiol. 2001;125(2):1074–85.
Tuteja N. Abscisic acid and abiotic stress signaling. Plant Signal Behav. 2007;2(3):135–8.
Danquah A, de Zelicourt A, Colcombet J, Hirt H. The role of ABA and MAPK signaling pathways in plant abiotic stress responses. Biotechnol Adv. 2014;32(1):40–52.
Dinh ST, Baldwin IT, Galis I. The HERBIVORE ELICITOR-REGULATED1 (HER1) gene enhances abscisic acid levels and defenses against herbivores in Nicotiana attenuata plants. Plant Physiol. 2013. doi:org/https://doi.org/10.1104/pp.113.221150.
Vos IA, Verhage A, Schuurink RC, Watt LG, Pieterse CM, Van Wees S. Onset of herbivore-induced resistance in systemic tissue primed for jasmonate-dependent defenses is activated by abscisic acid. Front Plant Sci. 2013. doi:org/https://doi.org/10.3389/fpls.2013.00539.
Bodenhausen N, Reymond P. Signaling pathways controlling induced resistance to insect herbivores in Arabidopsis. Mol Plant-Microbe Interact. 2007;20(11):1406–20.
Zong W, Tang N, Yang J, Peng L, Ma S, Xu Y, et al. Feedback regulation of ABA signaling and biosynthesis by a bZIP transcription factor targets drought resistance related genes. Plant Physiol. 2016. https://doi.org/10.1104/pp.16.00469.
Hillwig MS, Chiozza M, Casteel CL, Lau ST, Hohenstein J, Hernández E, et al. Abscisic acid deficiency increases defence responses against Myzus persicae in Arabidopsis. Mol Plant Pathol. 2016;17(2):225–35.
Zhu-Salzman K, Salzman RA, Ahn J-E, Koiwa H. Transcriptional regulation of sorghum defense determinants against a phloem-feeding aphid. Plant Physiol. 2004;134(1):420–31.
Yue J, Hu X, Huang J. Origin of plant auxin biosynthesis. Trends Plant Sci. 2014;19(12):764–70.
Berens ML, Berry HM, Mine A, Argueso CT, Tsuda K. Evolution of hormone signaling networks in plant defense. Annu Rev Phytopathol. 2017;55:401–25.
Attaran E, Major IT, Cruz JA, Rosa BA, Koo AJ, Chen J, et al. Temporal dynamics of growth and photosynthesis suppression in response to jasmonate signaling. Plant Physiol. 2014;165(3):1302–14.
Huot B, Yao J, Montgomery BL, He SY. Growth–defense tradeoffs in plants: a balancing act to optimize fitness. Mol Plant. 2014;7(8):1267–87.
Ding X, Cao Y, Huang L, Zhao J, Xu C, Li X, et al. Activation of the indole-3-acetic acid–amido synthetase GH3-8 suppresses expansin expression and promotes salicylate-and jasmonate-independent basal immunity in rice. Plant Cell. 2008;20(1):228–40.
Naseem M, Kaltdorf M, Dandekar T. The nexus between growth and defence signalling: auxin and cytokinin modulate plant immune response pathways. J Exp Bot. 2015;66(16):4885–96.
Grones P, Friml J. Auxin transporters and binding proteins at a glance. J Cell Sci. 2015;128:1–7.
Park J-E, Park J-Y, Kim Y-S, Staswick PE, Jeon J, Yun J, et al. GH3-mediated auxin homeostasis links growth regulation with stress adaptation response in Arabidopsis. J Biol Chem. 2007;282(13):10036–46.
Kazan K, Manners JM. Linking development to defense: auxin in plant–pathogen interactions. Trends Plant Sci. 2009;14(7):373–82.
Giron D, Frago E, Glevarec G, Pieterse CM, Dicke M. Cytokinins as key regulators in plant–microbe–insect interactions: connecting plant growth and defence. Funct Ecol. 2013;27(3):599–609.
Broekgaarden C, Voorrips RE, Dicke M, Vosman B. Transcriptional responses of Brassica nigra to feeding by specialist insects of different feeding guilds. Insect Sci. 2011;18(3):259–72.
Dubois M, Van den Broeck L, Inzé D. The pivotal role of ethylene in plant growth. Trends Plant Sci. 2018;23(4):311–23.
Diezel C, Allmann S, Baldwin IT. Mechanisms of optimal defense patterns in Nicotiana attenuata: flowering attenuates herbivory-elicited ethylene and jasmonate signaling. J Integr Plant Biol. 2011;53(12):971–83.
Lu J, Li J, Ju H, Liu X, Erb M, Wang X, et al. Contrasting effects of ethylene biosynthesis on induced plant resistance against a chewing and a piercing-sucking herbivore in rice. Mol Plant. 2014;7(11):1670–82.
Louis J, Basu S, Varsani S, Castano-Duque L, Jiang V, Williams WP, et al. Ethylene contributes to maize insect resistance1-mediated maize defense against the phloem sap-sucking corn leaf aphid. Plant Physiol. 2015;169(1):313–24.
Lv M, Kong H, Liu H, Lu Y, Zhang C, Liu J, et al. Induction of phenylalanine ammonia-lyase (PAL) in insect damaged and neighboring undamaged cotton and maize seedlings. Int J Pest Manag. 2017;63(2):166–71.
Chaman ME, Copaja SV, Argandoña VH. Relationships between salicylic acid content, phenylalanine ammonia-lyase (PAL) activity, and resistance of barley to aphid infestation. J Agric Food Chem. 2003;51(8):2227–31.
Sytykiewicz H, Chrzanowski G, Czerniewicz P, Sprawka I, Łukasik I, Goławska S, et al. Expression profiling of selected glutathione transferase genes in Zea mays (L.) seedlings infested with cereal aphids. PLoS One. 2014;9(11):e111863.
Hernández J, Barba-Espín G, Diaz-Vivancos P. Glutathione-Mediated Biotic Stress Tolerance in Plants. Glutathione in Plant Growth, Development, and Stress Tolerance. Switzerland: Springer; 2017. p. 309–29.
Zechmann B. Compartment-specific importance of glutathione during abiotic and biotic stress. Front Plant Sci. 2014. doi:org/https://doi.org/10.3389/fpls.2014.00566.
Kawano T. Roles of the reactive oxygen species-generating peroxidase reactions in plant defense and growth induction. Plant Cell Rep. 2003;21(9):829–37.
Gulsen O, Eickhoff T, Heng-Moss T, Shearman R, Baxendale F, Sarath G, et al. Characterization of peroxidase changes in resistant and susceptible warm-season turfgrasses challenged by Blissus occiduus. Arthropod Plant Interact. 2010;4(1):45–55.
Mai VC, 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.
Jun X, X-y WANG, W-z GUO. The cytochrome P450 superfamily: key players in plant development and defense. J Integr Agric. 2015;14(9):1673–86.
Rushton PJ, Somssich IE, Ringler P, Shen QJ. WRKY transcription factors. Trends Plant Sci. 2010;15(5):247–58.
Satapathy L, Kumar D, Mukhopadhyay KWRKY. Transcription factors: involvement in plant–pathogen interactions. In: Shukla P, editor. Recent advances in applied microbiology: Springer; 2017. p. 229–46.
Banerjee A, Roychoudhury A. WRKY proteins: signaling and regulation of expression during abiotic stress responses. Singapore: Sci World J. 2015. doi:org/https://doi.org/10.1155/2015/807560.
Kuśnierczyk A, Winge P, Jørstad TS, Troczyńska J, Rossiter JT, Bones AM. Towards global understanding of plant defence against aphids–timing and dynamics of early Arabidopsis defence responses to cabbage aphid (Brevicoryne brassicae) attack. Plant Cell Environ. 2008;31(8):1097–115.
Li J, Brader G, Kariola T, Tapio Palva E. WRKY70 modulates the selection of signaling pathways in plant defense. Plant J. 2006;46(3):477–91.
Higashi K, Ishiga Y, Inagaki Y, Toyoda K, Shiraishi T, Ichinose Y. Modulation of defense signal transduction by flagellin-induced WRKY41 transcription factor in Arabidopsis thaliana. Mol Genet Genomics. 2008;279(3):303–12.
Huang P-Y, Catinot J, Zimmerli L. Ethylene response factors in Arabidopsis immunity. J Exp Bot. 2016;67(5):1231–41.
Tsuda K, Somssich IE. Transcriptional networks in plant immunity. New Phytol. 2015;206(3):932–47.
Alves MS, Dadalto SP, Gonçalves AB, De Souza GB, Barros VA, Fietto LG. Plant bZIP transcription factors responsive to pathogens: a review. Int J Mol Sci. 2013;14(4):7815–28.
Nuruzzaman M, Sharoni AM, Kikuchi S. Roles of NAC transcription factors in the regulation of biotic and abiotic stress responses in plants. Front Microbiol. 2013. doi:org/https://doi.org/10.3389/fmicb.2013.00248.
Kelly J. Chapter 1-growth and development. In: Espinoza L, Kelly J, editors. Grain Sorghum production handbook. Cooperative Extention service: University of Arkansas. Fayetteville, AR. 2004. https://www.uaex.edu/publications/pdf/mp297/MP297.pdf.
Chisholm ST, Coaker G, Day B, Staskawicz BJ. Host-microbe interactions: shaping the evolution of the plant immune response. Cell. 2006;124(4):803–14.
Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294–9.
Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, et al. The Sorghum bicolor genome and the diversification of grasses. Nature. 2009;457(7229):551–6.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38:W64–70.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39(2):W316–W22.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8.
Team RC. R: a language and environment for statistical computing. Vienna: R Foundation for statistical computing. 2012. ISBN 3-900051-07-02012
We would like to thank Eva Huff for assistance in the greenhouse, Li Paetzold for help in the laboratory, and Kristyne A. Varela for assistance in the sample preparation. We would also like to acknowledge three anonymous reviewers whose comments and suggestions greatly improved this manuscript.
This research was supported by funding from USDA NIFA (TEX09638), USDA-ARS, Areawide Pest Management of the Invasive Sugarcane Aphid in Grain Sorghum (3072–22000-016-00D), and Texas A&M AgriLife Genomics and Bioinformatics Center. The funders had no role in the design of the study, collection, analysis, interpretation of the data, or in writing the manuscript.
Availability of data and materials
All raw sequencing reads have been submitted to the NCBI Sequence Read Archive and are available under BioProject ID: PRJNA451032 (http://www.ncbi.nlm.nih.gov/bioproject/451032).
Ethics approval and consent to participate
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 RNA-seq reads from susceptible and resistant sorghum genotypes mapped to the sorghum genome (Sbicolor_454 version v.3.1.1). Unique RNA-seq reads mapping to exons, introns, and intergenic regions are shown as the percentage of total reads distributed to these annotated regions of the sorghum genome. (DOCX 23 kb)
Differentially expressed genes (DEGs) in 2-wk resistant sorghum genotype in response to aphid herbivory. DEGs were defined as having a fold change ≥1.9 or ≤ − 1.9 with a false discovery rate (FDR) adjusted p-value < 0.05. (XLSX 431 kb)
Differentially expressed genes (DEGs) in 2-wk susceptible sorghum genotype in response to aphid herbivory. DEGs were defined as having a fold change ≥1.9 or ≤ − 1.9 with a false discovery rate (FDR) adjusted p-value < 0.05. (XLSX 258 kb)
Differentially expressed genes (DEGs) in 6-wk resistant sorghum genotype in response to aphid herbivory. DEGs were defined as having a fold change ≥1.9 or ≤ − 1.9 with a false discovery rate (FDR) adjusted p-value < 0.05. (XLSX 339 kb)
Differentially expressed genes (DEGs) in 6-wk susceptible sorghum genotype in response to aphid herbivory. DEGs were defined as having a fold change ≥1.9 or ≤ − 1.9 with a false discovery rate (FDR) adjusted p-value < 0.05. (XLSX 106 kb)
Real-time quantitative RT-PCR (RT-qPCR) validation of gene expression changes detected by RNA-seq. Data are shown as a fold change of the mean expression levels from three biological replicates. The correlation between the RNA-seq data and the RT-qPCR results is shown by Pearson’s correlation coefficient (R2). (PDF 140 kb)
Primers sequences used for RT-qPCR amplification of the five differentially expressed genes selected for validation. (DOCX 16 kb)
About this article
Cite this article
Kiani, M., Szczepaniec, A. Effects of sugarcane aphid herbivory on transcriptional responses of resistant and susceptible sorghum. BMC Genomics 19, 774 (2018). https://doi.org/10.1186/s12864-018-5095-x
- Melanaphis sacchari
- Sorghum bicolor
- Plant-insect interaction
- Plant defenses