- Research
- Open access
- Published:
Transcriptomics identifies key defense mechanisms in rice resistant to both leaf-feeding and phloem feeding herbivores
BMC Plant Biology volume 21, Article number: 306 (2021)
Abstract
Background
Outbreaks of insect pests in paddy fields cause heavy losses in global rice yield annually, a threat projected to be aggravated by ongoing climate warming. Although significant progress has been made in the screening and cloning of insect resistance genes in rice germplasm and their introgression into modern cultivars, improved rice resistance is only effective against either chewing or phloem-feeding insects.
Results
In this study, the results from standard and modified seedbox screening, settlement preference and honeydew excretion tests consistently showed that Qingliu, a previously known leaffolder-resistant rice variety, is also moderately resistant to brown planthopper (BPH). High-throughput RNA sequencing showed a higher number of differentially expressed genes (DEGs) at the infestation site, with 2720 DEGs in leaves vs 181 DEGs in sheaths for leaffolder herbivory and 450 DEGs in sheaths vs 212 DEGs in leaves for BPH infestation. The leaf-specific transcriptome revealed that Qingliu responds to leaffolder feeding by activating jasmonic acid biosynthesis genes and genes regulating the shikimate and phenylpropanoid pathways that are essential for the biosynthesis of salicylic acid, melatonin, flavonoids and lignin defensive compounds. The sheath-specific transcriptome revealed that Qingliu responds to BPH infestation by inducing salicylic acid-responsive genes and those controlling cellular signaling cascades. Taken together these genes could play a role in triggering defense mechanisms such as cell wall modifications and cuticular wax formation.
Conclusions
This study highlighted the key defensive responses of a rarely observed rice variety Qingliu that has resistance to attacks by two different feeding guilds of herbivores. The leaffolders are leaf-feeder while the BPHs are phloem feeders, consequently Qingliu is considered to have dual resistance. Although the defense responses of Qingliu to both insect pest types appear largely dissimilar, the phenylpropanoid pathway (or more specifically phenylalanine ammonia-lyase genes) could be a convergent upstream pathway. However, this possibility requires further studies. This information is valuable for breeding programs aiming to generate broad spectrum insect resistance in rice cultivars.
Background
Rice (Oryza sativa L.) is cultivated on 160 million hectares of land in over a hundred countries, and approximately 491 million tons of milled rice are produced annually [1]. Rice is the most important staple food for half of the world’s population, especially in Asia and Africa [2]. However, insect pests pose a long-standing threat to rice farming worldwide [3]. There are over 100 species of rice-feeding insect pests, of which approximately 20 species are considered major pests that severely affect rice yield [4]. Recently, frequent and widespread outbreaks of major insect pests of rice, including the brown planthopper (BPH; Nilaparvata lugens) and rice leaffolder (Cnaphalocrocis medinalis), have caused losses of hundreds of millions of dollars annually and threatened food security [3, 5]. Insect pests cause severe damage to rice plants, as these pests thrive in the warm and humid environment of fertile paddy fields [6]. Moreover, insect pest infestations in paddy fields are predicted to become more rampant and serious in a warming climate that increases the metabolic rates and population growth of insects [7].
For over 50 years, spraying chemical insecticides has been the main strategy for managing crop pests [8]. However, the misuse of insecticides compromises the food safety and disrupts ecological and environmental integrity, which could lead to secondary pest outbreaks [9]. Therefore, policy makers in different countries are beginning to adopt environmentally friendly pest management strategies to minimize insecticide use [10]. For example, integrated pest management (IPM) is a long-term, environmentally friendly strategy that aims to control pests through a combination of approaches such as cultural practices, biological control and the use of resistant varieties [11]. Regarding the latter, a concerted effort has been made over the past decades to screen for insect resistance in rice germplasm, identify resistance genes, characterize the molecular mechanisms of host resistance, and introgress the resistance genes into modern cultivars [12]. To date, 34 resistance loci against the BPH have been identified in rice germplasm, of which 14 loci have been cloned and 5 loci have been functionally characterized (Bph3, Bph14, Bph18, Bph26 and Bph29), while fewer than 10 genes have been introgressed into modern rice cultivars [12, 13].
Rice-feeding insect pests have different feeding guilds, but they can generally be divided into two groups: chewing insects and piercing-sucking insects [14]. Chewing insects, such as leaffolders and stem borers, have mouthparts that tear off and chew the foliage, causing extensive mechanical wounding to host plants [15]. Mechanical wounding prompts jasmonic acid (JA) production and signaling pathways, activating defense mechanisms such as the oxidative burst, induction of JA-responsive genes and generation of trypsin protease inhibitors that retard the feeding and growth of chewing insects [12]. On the other hand, piercing-sucking insect pests, such as planthoppers and leafhoppers, have specialized slender, tapering stylet mouthparts that penetrate the apoplastic space of epidermal cells and suck phloem sap [16]. Hoppers are also vectors of plant viruses, causing secondary damage to plants [17]. Previous studies indicated that salicylic acid (SA) biosynthesis and signaling are induced in rice plants in response to feeding by the BPH [18, 19].
Several studies of mutant rice lines defective in JA biosynthesis and signaling have confirmed the involvement of different phytohormone pathways in their defenses against chewing and phloem-feeding insects [20]. For example, knockout of a JA-biosynthetic lipoxygenase gene (LOX) increased the susceptibility of rice plants to striped stem borers (chewing insects), but improved their resistance to BPHs (phloem-feeding insects), presumably due to decreased JA and increased SA and hydrogen peroxide levels [21]. In addition, knockout of the coronatine-insensitive 1 gene (COI1; encodes JA receptor) increased the susceptibility of rice plants to chewing insects but did not alter rice resistance/susceptibility to BPHs, implying that JA signaling is not implicated in the rice response to the BPH [22]. Taken together, the antagonistic crosstalk between JA and SA pathways is a factor compromising the identification of plants with dual resistance to chewing and piercing-sucking insects [23].
Previous studies consistently showed that the Qingliu (QL) rice variety is resistant to leaffolder infestation [24,25,26]. Qingliu is a local native indica rice variety in Taiwan where its cultivation record can be dated back to the Qing dynasty era around 1742 [27]. In the past, this native variety was grown on a large scale because of its short life cycle, high yield and seed quality [27]. However, given that Qingliu is excessively tall and tends to lodge, its cultivation was eventually phased out in favor of modern semidwarf rice varieties developed during the Green Revolution [27]. Recently, a comparative proteomic study showed that the enzymes involved in flavonoid and JA biosynthesis were expressed at higher levels in Qingliu than in the susceptible check Taichung Native 1 (TN1) before and during leaffolder herbivory, potentially contributing to Qingliu resistance to chewing insects [28]. However, the RNA sequencing (RNA-Seq) profiles of Qingliu in response to leaffolder herbivory are not yet available. The objectives of this study are to (i) assess whether Qingliu has resistance to phloem-feeding BPHs and (ii) highlight the transcriptomic changes in key defense-related processes in Qingliu subjected to leaffolder and BPH infestation.
Results
Qingliu is moderately resistant to brown planthopper infestation
The standard seedbox screening test (SSST) showed that Qingliu was relatively less severely infested by BPH (average damage score = 5.0) than the two susceptible cultivars TN1 (8.3) and TNG67 (7.0) (two-way ANOVA and LSD post hoc test, p-value < 0.001) (Table 1). In addition, Qingliu also had the same average damage score as two moderately resistant checks: Mudgo and H105. As expected, Baiqiaowan, a leaffolder-resistant check, was highly susceptible to BPH infestation (9.0) (Table 1). In addition, the modified seedbox screening test (MSST) revealed a consistent result with the SSST (two-way ANOVA and LSD post hoc test, p-value < 0.001) (Table 1). These results indicated that Qingliu is moderately resistant to BPH infestation. In the choice test between Qingliu and TN1, a lower number of BPHs settled on Qingliu than on TN1 after 3, 6, 24 and 48 h of insect treatment (z-test, p-value = 0.0002, 0.0003, 0.0124, 0.0213 at the four time points) (Fig. 1a).
In the honeydew excretion test, individual BPHs feeding on Qingliu excreted a lower volume of phloem-derived excretion than those feeding on TN1 (two-way ANOVA, p-value = 0.0100) (Fig. 1b). This suggests that BPHs feeding on Qingliu either had less access to phloem sap or lower phloem sap intake rate than those feeding on TN1. Moreover, individual BPHs feeding on either Qingliu or TN1 excreted a similar volume of xylem-derived excretions (two-way ANOVA, p-value = 0.1443) (Fig. 1b). This implies that the BPHs feeding on Qingliu and TN1 either had equal access to xylem sap or similar xylem sap intake rates.
A greater transcriptional change occurs at the infestation site of Qingliu than in other tissues in response to leaffolder and brown planthopper infestation
To discover the molecular responses of Qingliu potentially conferring the observed dual resistance to chewing and piercing-sucking insects, RNA-Seq profiles of the leaf and sheath tissues were generated under uninfested control, leaffolder-infested (24 h) and BPH-infested (24 h) conditions. A total of thirty RNA-Seq samples were generated with five biological replicates in each experimental group (2 tissues × 3 conditions × 5 biological replicates). RNA-Seq produced a large number of raw reads per sample, falling within the range of 39,700,000-52,700,000 (Table S1). A total of 36,000,000-49,200,000 clean reads per sample (90–93% of raw reads) passed the read trimming and filtering steps (Table S1). At least 92.71% of the clean reads per sample were mapped to the rice reference genome with at least 89.16% of clean reads mapped to a single genomic locus, confirming the quality of the RNA-Seq samples (i.e., free of microbial contamination) (Table S1).
Principal component analysis (PCA) was carried out to assess the homogeneity of the biological replicates (Fig. S1). At principal component 1 (13.25%), which showed the largest variation among the RNA-Seq samples, all fifteen leaf samples were segregated from the fifteen sheath samples. The second largest variation among the RNA-Seq samples (7.66%) was contributed by the relatively high heterogeneity among the sheath samples compared with that of the leaf samples. In each tissue, the samples of three different experimental conditions did not segregate considerably, implying the mild effect of both insect infestations on transcriptomic changes in Qingliu. Nevertheless, the biological replicates of each experimental group formed their own clusters (Fig. S1).
At 24 h of leaffolder infestation, a higher number of DEGs (2720 vs 181 DEGs) were identified in the leaves than in the sheaths (Fig. 2; Table S2a, b). This finding was expected because the leaf is the infestation site of leaffolders. Among the DEGs in leaves, 1767 genes (64.96%) were upregulated and 953 genes (35.04%) were downregulated (Fig. 2; Table S2a). In contrast, more DEGs were detected in the sheaths than in the leaves at 24 h of BPH infestation: 450 vs 212 DEGs (Fig. 2; Table S2c, d). This was also expected because the sheath is the infestation site of the BPH. Among the DEGs in sheaths, 292 genes (64.89%) were upregulated and 158 genes (35.11%) were downregulated (Fig. 2; Table S2d). In addition, DEGs commonly regulated by both types of insect pests were also identified: 24 upregulated and 75 downregulated genes in leaves and 27 upregulated and 5 downregulated genes in sheaths (Fig. 2; Table S2e, f). For example, several defense-associated genes encoding Bowman-Birk trypsin inhibitors (Os01g0124200), peroxidases (Os03g0235000, Os10g0109600), BURP and universal stress protein A (UspA) domain-containing proteins (Os05g0217800, Os05g0217700, Os03g0305400) were commonly induced in leaves under separate leaffolder and BPH infestations (Table S2e).
Nine defense-related or sample-specific DEGs were shortlisted for qRT-PCR validation (Fig. 3). The qRT-PCR results of eight genes showed a strong correlation with the RNA-Seq results (except AK064067), attesting to the robustness of the generated transcriptomes.
RNA-Seq reveals the involvement of JA biosynthesis and SA signaling in the Qingliu response to leaffolder herbivory and SA signaling against brown planthopper infestation
To determine the defense phytohormones controlling Qingliu defense responses against either leaffolder or BPH infestation, the DEGs from plants infested with each insect were separately mapped to 798, 615 and 183 previously reported SA-, JA- and ethylene (ACC)-responsive genes, respectively [29]. A total of 80 SA-, 95 JA- and 25 ACC-responsive genes were mapped to the DEGs (Fig. 4; Table S3).
Notably, 86.3, 89.5 and 84.0% of the total mapped SA-, JA- and ACC-responsive genes, respectively, were induced by leaffolder feeding on the leaves, indicating the involvement of these three defense phytohormones in the modulation of the Qingliu defense response against this chewing insect (Fig. 4; Table S3). Induced expression of marker genes associated with the SA response were identified. These included abnormal inflorescence meristem 1 OsAIM1 (Os02g0274100) that is involved in SA biosynthesis and two genes involved in SA signaling, nonexpressor of pathogenesis-related genes OsNPR1/3 (Os03g667100 and Os01g0194300 (Table S2). For JA, although the expression of three JA biosynthesis genes was upregulated: lipoxygenases OsLOX2.1/2.2/2.3 (Os08g0508800, Os12g0559200, Os08g0509100), seven genes encoding the repressors of JA signaling: jasmonate ZIM-domain proteins OsJAZs (Os03g0180900, Os10g0392400, Os07g0615200, Os03g0402800, Os09g0439200, Os03g0180800, Os04g0395800) were also upregulated (Table S2a). This result suggested the involvement of SA biosynthesis and signaling as well as JA biosynthesis in the defense response of Qingliu against leaffolder infestation.
Next, 11.3, 8.4 and 16% of the total mapped SA-, JA- and ACC-responsive genes, respectively, were induced in the BPH-infested sheaths (Fig. 4; Table S3). However, two genes encoding repressors of JA signaling, OsJAZs (Os03g0181100, Os03g0180800), were induced (Table S2d). This finding implicated the involvement of SA signaling in the defense response of Qingliu against BPH infestation. Indeed, among the nine SA-responsive genes upregulated in the BPH-infested sheaths, four genes were associated with transcription factor (Os08g0474000), cell wall modification (Os04g0604300, Os12g0443500) and lipid metabolism (Os03g0718800) (Fig. 4a; Table S3). In the coming sections, the relevance of these components in Qingliu defense response against the BPH will be explained based on Gene Ontology (GO) enrichment and MapMan annotation results.
A defense response involving JA biosynthesis and the shikimate pathway is activated in the leaves of Qingliu plants upon leaffolder herbivory
To understand the defense response of Qingliu against leaffolder herbivory, GO enrichment analysis was performed separately on the up- and downregulated genes in each tissue. A higher number of enriched GO terms was found in the leaves (infestation site) than in the sheaths: upregulated genes in leaves were annotated to 68 GO terms, downregulated genes in leaves to 21 GO terms, upregulated genes in sheaths to 15 GO terms and downregulated genes in sheaths to 0 GO terms (Table S4). After removing the interrelated GO terms by manually curating the GO hierarchical tree graphs, upregulated genes in leaves were annotated to 22 GO terms, downregulated genes in leaves to 4 GO terms (mainly related to chlorophyll biosynthesis) and upregulated genes in sheaths to 1 GO term (Fig. S2).
The 22 enriched GO terms of upregulated genes in leaffolder-infested leaves consisted of 12 biological process terms, 2 cellular component terms and 8 molecular function terms (Fig. 5a; Table S4a). Based on the biological process annotation, it can be deduced that phosphorylation cascades (GO:0006468: protein amino acid phosphorylation) were activated in Qingliu as an early response to leaffolder infestation. The signaling cascades triggered JA biosynthesis and defense mechanisms in Qingliu (GO:0031407: oxylipin metabolic process, GO:0006952: defense response, GO:0009607: response to biotic stimulus, GO:0009611: response to wounding). Notably, the defense mechanisms involved an induction of the genes involved in cell wall catabolism, biosynthesis of aromatic amino acids (tryptophan, phenylalanine and tyrosine) via the shikimate pathway, and production of phenylpropanoid and terpenoid secondary metabolites (GO:0016998: cell wall macromolecule catabolic process, GO:0016052: carbohydrate catabolic process, GO:0019438: aromatic compound biosynthetic process, GO:0006568: tryptophan metabolic process, GO:0046417: chorismate metabolic process, GO:0006721: terpenoid metabolic process).
To demonstrate the role of phosphorylation cascades in activating the defense mechanisms of Qingliu against leaffolder infestation, coexpression analysis was performed using the 101 protein kinase genes annotated to GO:0006468: protein amino acid phosphorylation (Fig. 5a; Table S4a). At hierarchy = 0 and mutual rank (MR) = 14, four groups of coexpressed protein kinase genes were identified: (1) Os07g0542400, Os07g0537500, Os07g0537600, Os08g0203100; (2) Os02g0165100, Os05g0127300; (3) Os03g0773300, Os06g0602500; and (4) Os01g0664200, Os04g0514800. Hence, these four groups of coexpressed genes were separately used as guide genes to generate coexpression networks. At hierarchy = 2 and MR = 7, two biologically meaningful coexpression networks were generated from the first two groups of guide genes (Fig. 5b, c). The first coexpression network highlighted upregulated genes encoding (i) a lipoxygenase (Os08g0508800) involved in JA biosynthesis, and (ii) two peroxidases (Os03g0235000, Os08g0113000) annotated to three Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of osa01110: biosynthesis of secondary metabolites, osa00360: phenylalanine metabolism and osa00940: phenylpropanoid biosynthesis (Fig. 5b). The second coexpression network revealed induced expression of genes encoding another lipoxygenase (Os03g0700700) and a 4-coumarate:coenzyme A ligase (Os08g0448000) involved in lignin biosynthesis for defense against wounding (Fig. 5c). Notably, a phenylalanine ammonia-lyase gene (Os02g0626600) also was found in the coexpression network (Fig. 5c).
A defense response involving cell wall modification and lipid metabolism is activated in the sheaths of Qingliu plants upon brown planthopper infestation
To decipher the defense response of Qingliu against BPH infestation, GO enrichment analysis was conducted separately on the up- and downregulated genes in each tissue. More enriched GO terms were discovered in the sheaths (infestation site) than in the leaves: upregulated genes in sheaths were annotated to 44 GO terms, downregulated genes in sheaths to 0 GO terms, upregulated genes in leaves to 2 GO terms and downregulated genes in leaves to 4 GO terms (Table S5). After removing the interrelated GO terms, upregulated genes in sheaths were annotated to 12 GO terms, upregulated genes in leaves to 1 GO term and downregulated genes in leaves to 1 GO term (Fig. S3).
The 12 enriched GO terms of upregulated genes in BPH-infested sheaths consisted of 5 biological process terms, 4 cellular component terms and 3 molecular function terms (Fig. 6a). The biological process annotation revealed that Qingliu responded to BPH infestation by undergoing transcriptional reprogramming (GO:0006355: regulation of transcription, DNA-dependent). This transcriptional modification activated defense mechanisms that regulate cell wall modification and lipid metabolism/transport (GO:0006950: response to stress, GO:0006073: cellular glucan metabolic process, GO:0006629: lipid metabolic process, GO:0006869: lipid transport). As expected, the molecular function and cellular component annotation supported the aforementioned defense response (GO:0003700: transcription factor activity, GO:0005634: nucleus, GO:0004553: hydrolase activity, GO:0016762: xyloglucan:xyloglucosyl transferase activity, GO:0005618: cell wall) (Fig. 6a; Table S5a).
To validate the defense response of Qingliu against BPH infestation shown by the GO enrichment, the DEGs in the sheaths were mapped to MapMan pathways (Fig. 6b-d; Table S6). As a result, eighteen genes encoding receptor kinases were responsive to BPH infestation with the following five genes being upregulated: Os01g0741200, Os01g0821900, Os02g0190500, Os03g0717000 and Os08g0374701 (Fig. 6b; Table S6a). The activated receptor kinases potentially trigger calcium, phophoinositide and G-protein signaling pathways that transmit signals to the nucleus via a series of protein phosphorylation events (Fig. 6c; Table S6b). In agreement with the GO enrichment results, this signal activates transcription factors (Fig. 6c; Table S6b) that control the expression of genes required for defense mechanisms, particularly cell wall modification and lipid metabolism (Fig. 6d; Table S6c).
Discussion
Dual resistance of Qingliu to leaffolders and brown planthoppers
Screening rice germplasm for dual resistance to chewing and phloem-feeding insect pests is an important goal in breeding programs aiming to generate broad-spectrum insect-resistant rice cultivars [12]. However, this effort could be hampered by the antagonistic roles of JA and SA in the mediation of rice defense mechanisms against chewing and phloem-feeding insects [23]. For example, previous studies reported that Baiqiaowan is resistant to leaffolders because it sustained less leaf roll damage and suppressed larval growth more than susceptible rice lines [24, 26]. Nevertheless, the SSST and MSST results in this study showed that Baiqiaowan is highly susceptible to BPH infestation and cannot be considered to have dual resistance (Table 1).
Previous studies consistently showed that Qingliu is resistant to leaffolder infestation [24,25,26]. This work began with natural insect infestation screening tests of 270 rice genotypes conducted at an experimental farm in Taiwan [24]. Among the plants tested, Qingliu had the lowest percentage of leaf roll damage [24]. A previous study showed that in both intact plants and detached leaves, third-instar leaffolder caterpillars that fed on Qingliu recorded a higher mortality rate and lower relative growth rate than those fed on the leaffolder-susceptible TN1 [26]. A number of unique attributes were previously identified in Qingliu, such as a high leaf trichome density, that could impede leaffolder infestation [25].
Although Qingliu is known to be a resistant to chewing insect pests, its resistance to phloem-feeding insect pests such as the BPH has not been reported. In this study, Qingliu was shown to be moderately resistant to the BPH (Fig. 1; Table 1). For example, based on the damage scores derived from the SSST and MSST experiments, Qingliu appeared to be moderately resistant to BPH infestation (Table 1). In the choice test, the piercing-sucking insect showed a lower settling preference for Qingliu than for TN1 (Fig. 1a). Furthermore, honeydew excretion was found to be lower for the individual BPHs feeding on Qingliu than TN1 (Fig. 1b).
This study showed that Qingliu displays a rarely observed dual resistance to leaffolder and BPH. However, according to the growth-defense tradeoff hypothesis, plant defense activation generally comes at the expense of plant growth due to limited metabolic resources that can be allocated either for growth or defense [30]. For example, growth inhibition by SA can be observed in constitutive defense mutants, which have a dwarf phenotype due to elevated SA accumulation or signaling [31]. Over the years it has been demonstrated that Qingliu has promising growth and yield performance under uninfested conditions [27], but it is unknown whether these indicators will drop drastically under infested conditions due to defense activation. Therefore, further phenotypic and genetic characterization of Qingliu is required prior to generating genetic markers for use in breeding programs aiming to generate broad-spectrum insect-resistant rice cultivars.
The importance of shikimate and phenylpropanoid pathways in Qingliu-leaffolder interactions
The elicitors in the salivary secretions of leaf-chewing insects trigger host defenses via an intricate wound-signaling pathway in which JA biosynthesis plays an important role [18, 32]. This explains the the induction of lipoxygenases (JA biosynthesis genes; Os08g0508800, Os12g0559200, Os08g0509100) in the leaffolder-infested leaves (Table S2a). Although JA biosynthesis in rice plants increases under leaffolder herbivory, a previous study showed that Qingliu produces a higher level of JA than the leaffolder-susceptible TN1 [26]. JA biosynthesis was reported to increase rice resistance to leaffolders by activating the phenylpropanoid pathway [33]. Here, increased transcript levels were observed for biosynthesis genes in the shikimate pathway that produce three essential aromatic amino acids: tryptophan, phenylalanine and tyrosine (Fig. 5a; Table S4a). A recent comparative proteomic study showed that enzymes (Q10NY1, Q5VNW0) catalyzing the synthesis or chorismate, the precursor of aromatic amino acids were at higher levels in Qingliu than in TN1 [28]. In addition to being used for protein synthesis, these essential aromatic amino acids also serve as precursors of various defensive compounds, such as SA, polyphenolic compounds and lignin, via the phenylpropanoid pathway in plants [32, 34].
Tryptophan-derived secondary metabolites in plants generally include indole alkaloids, indole glucosinolates, camalexin, auxin, serotonin and melatonin [35]. In response to leaffolder herbivory, the induction of genes encoding biosynthetic enzymes for serotonin (Os08g0140300, Os08g0140500) and melatonin (Os08g0157500) was observed (Table S2a). A recent study showed that suppression of serotonin biosynthesis enhanced the resistance of rice to BPH and striped stem borers, presumably by elevating phenylalanine-derived SA levels [36]. Based on this finding, the observed increased biosynthesis of serotonin in response to leaffolder feeding might serve as a precursor for melatonin biosynthesis. Melatonin in plants delays senescence by enhancing the activity of antioxidant enzymes and increasing flavonoid biosynthesis [37].
Phenylalanine and tyrosine are channeled into the phenylpropanoid pathway to produce defensive secondary metabolites – flavonoids [38, 39]. Phenylalanine ammonia lyase (PAL) catalyzes the first and committed step in the phenylpropanoid pathway, synthesizing trans-cinnamate from phenylalanine [40]. The trans-cinnamate is then hydroxylated to form p-coumarate [41]. Recent studies indicate that the monocot grass family Poaceae contains a bifunctional phenylalanine/tyrosine ammonia lyase (PTAL), which can synthesize p-coumarate from tyrosine, bypassing the trans-cinnamate intermediate [35, 39]. In the RNA-Seq expression profile of Qingliu in response to leaffolder herbivory, the upregulation of genes encoding OsPAL6 (Os04g0518400) was found together with the following biosynthesis enzymes for flavonoids: chalcone synthase (CHS; Os07g0526400), chalcone isomerase (CHI; Os06g0203600), chalcone reductase (Os04g0167800), flavanone 3-hydroxylase (Os03g0122300) and leucoanthocyanidin dioxygenase (Os01g0832600, Os03g0289800, Os05g0127500) (Table S2a). These results are supported by a recent comparative proteomic study, which showed that PAL, CHS and CHI were expressed at higher levels in Qingliu than in TN1 before and/or during infestation [28].
Lignin is also produced by the phenylalanine/tyrosine phenylpropanoid pathway in monocotyledonous plants [42, 43]. Given protective role of the plant cell wall, lignin biosynthesis and deposition can be triggered in response to wounding caused by chewing insect pests [44]. In the RNA-Seq results, the induction of genes encoding the biosynthesis enzymes for lignin monomers: 4-coumarate:CoA ligase (Os08g0448000, Os02g0177600), caffeic acid O-methyltransferase (Os08g0157500), caffeoyl CoA O-methyltransferase (Os08g0498400, Os08g0498100), coniferaldehyde 5-hydroxylase (Os10g0512400), cinnamyl alcohol dehydrogenase (Os02g0187800) and cinnamoyl-CoA reductase (Os08g0441500) was observed (Table S2a). The lignin monomers are eventually transported into the apoplast where they are polymerized by peroxidase and laccase [45]. In agreement with this finding, a study reported that the insect-specific toxin scorpion peptide LqhIT2 improved rice resistance to leaffolders by promoting JA biosynthesis, which in turn increased lignin deposition in the cell walls of leaf tissues [33].
The PAL pathway is the main route of SA biosynthesis in rice [46] and in this study the induction of marker genes associated with SA biosynthesis was found: OsAIM1 (Os02g0274100) and SA signaling: OsNPR1/3 (Os03g0667100, Os01g0194300) (Table S2a). A previous study showed that Qingliu has higher levels of SA at 6 h of leaffolder feeding, and the SA level in Qingliu is constitutively higher than that in leaffolder-susceptible TN1 prior to insect feeding [26]. The elevated SA biosynthesis in Qingliu under leaffolder herbivory might account for the observed suppressive effect on JA signaling due to the well-known antagonistic interaction between SA and JA [19]. This explains the paradoxical observation of the upregulation of genes involved in JA biosynthesis and negative regulators of JA signaling in Qingliu in response to leaffolder infestation (Table S2a). While the functional importance of SA biosynthesis and signaling is well-documented against phloem-feeding insect pests, the roles of SA against chewing insect pests still require further investigation [18, 47, 48]. Unless there is evidence suggesting currently unknown functions of SA signaling, the increased formation of SA through the PAL pathway might merely serve to repress JA signaling, whereas accumulation of phenylpropanoid-derived secondary metabolites mainly account for Qingliu resistance to leaffolders.
The importance of cell wall modification and lipid metabolism in Qingliu-brown planthopper interactions
Compared with leaffolder larvae that cause extensive leaf damage, BPHs that insert their needle-like stylets through the apoplastic space between epidermal cells to reach sap-containing phloem cells cause little mechanical damage to host tissues [12]. By distinguishing the degree of tissue damage at the feeding site, rice plants are able to employ distinct defense mechanisms against these two types of insect pests with different feeding guilds [49]. To date, five BPH resistance loci have been functionally characterized by map-based gene cloning: Bph3, Bph14, Bph18, Bph26 and Bph29 [18, 47, 50,51,52]. Except for Bph29, which encodes a resistance gene containing a B3 DNA-binding domain (a potential transcription factor), the other four loci encode either coiled-coil, nucleotide-binding site and leucine-rich repeat proteins (CC-NBS-LRR) or lectin receptor kinases (LecRK), which activate the SA signaling pathway and induce callose deposition in phloem cells [18, 47, 50,51,52].
In this study, induced expression of five receptor kinase genes was found in BPH-infested sheaths at 24 h: Os01g0821900, Os08g0374701, Os03g0717000, Os02g0190500 and Os01g0741200 (Table S6a). Three of the receptor kinase genes have been reported as upstream regulators in signal transduction against abiotic and biotic stresses [53,54,55]. For example, rice receptor-like cytoplasmic kinase 253 (OsRLCK253; Os08g0374701) interacts with stress-associated proteins 1/11 (OsSAP1/11) in the plasma membrane, nuclear membrane and nucleus, and their overexpression in Arabidopsis improves water deficit and salinity tolerance [54]. An ortholog of rice lysin motif-type receptor-like kinase 4 (OsLysM-RLK4; Os01g0741200) in Arabidopsis triggers plant innate immunity against pathogens by functioning as a chitin recognition receptor to turn on the chitin signaling pathway [55].
Genes associated with calcium signaling (Os09g0482800, Os04g0492800, Os06g0683400, Os01g0949500, Os01g0955100), phosphoinositide lipid signaling (Os02g0285300) and G-protein signaling (Os02g0719000, Os05g0513800, Os05g0454200) were induced in BPH-infested sheaths at 24 h (Table S6b). These signaling pathways could lead to the activation of transcription factors in the nucleus. As expected, transcript levels of 28 genes mainly encoding transcription factors of the AP2/ERF (8 genes), MYB (5 genes), bZIP (3 genes), bHLH (3 genes) and WRKY (2 genes) families (Table S7a) were induced. Among the eight AP2/ERF transcription factor genes, four genes belonged to the OsDREB (dehydration response element binding) subfamily: OsDREB1A/C/E/H. These OsDREB transcription factors and OsbHLH148 are well known to control the rice response to drought, cold and salt stresses [56, 57]. Interestingly, the induction of a major OsMYB4 transcription factor gene was observed, where a previous study comprehensively showed that the constitutive expression of OsMYB4 in Arabidopsis improves tolerance to drought, salt, ultraviolet radiation, ozone, viruses, bacteria and fungi [58]. The induction of OsWRKY76 was also observed – its overexpression in rice plants was reported to increase susceptibility to Magnaporthe oryzae but enhance tolerance to cold stress [59].
Transcription factors reprogram the expression of genes controlling various downstream defense mechanisms in BPH-infested rice sheaths. The cell wall is the main “warzone” in the rice-BPH interaction; therefore, it is unsurprising that cell wall-related GO terms were enriched in the infested sheaths (Fig. 6) [60]. The salivary secretions of the BPH may contain cell wall modifying enzymes (CWMEs) that facilitate the penetration of stylets through the extracellular matrix of epidermal cells [61]. To counter this, rice plants also employ another set of CWMEs to fortify the cell wall, thwarting the feeding process of phloem sap-sucking pests [60]. Consistently, the expression of genes encoding enzymes essential for the biosynthesis, degradation and modification of major cell wall components, including pectin, cellulose, hemicellulose, and callose was induced in the BPH-infested rice sheaths at 24 h (Table S7b). For example, pectin lyase (Os12g0554800) catalyzes the depolymerization of pectin to oligogalacturonides, which function as damage-associated molecular patterns (DAMPs), activating the defense response through the mitogen-activated protein kinase (MAPK) signaling pathway [62, 63]. Xyloglucan, a hemicellulose, is a direct target degraded by pathogens for successful colonization of host plants [64]. To prevent the degradation of xyloglucan, the expression of four cell wall-localized xyloglucan endotransglucosylases/hydrolases (Os04g0604300, Os02g0696500, Os06g0696400, Os03g0239000), which serve to modify the xyloglucan chains, was induced in the BPH-infested sheaths of Qingliu (Table S7b) [64]. In addition, depolymerized xyloglucan oligomers can also function as DAMPs, triggering the MAPK signaling cascade [65]. Instead of callose synthase, induced expression of eight genes encoding β-1,3-glucanase was found, indicating that callose deposition in the phloem might not be the mechanism underlying the resistance of Qingliu to BPH (Table S7b) [66].
Lipid metabolism including lipid-mediated signaling and wax biosynthesis was another process regulated in BPH-infested rice sheaths (Table S7c). For lipid signaling, increased expression of genes encoding phospholipase C (Os02g0588500) and two phospholipase D enzymes (Os08g0401800, Os02g0120200) was found (Table S7c). Under environmental stress, phospholipase enzymes catalyze the hydrolysis of membrane phospholipids to generate signal messengers such as phosphatidic acid, diacylglycerol, phosphoinositides and inositol polyphosphates [67]. Increased expression of a fatty acid 2-hydroxylase gene (OsFAH1; Os12g0628400), which functions in the biosynthesis of sphingolipids, was also identified (Table S7c) [68]. Proper production of sphingolipids protects the integrity of plasma membrane microdomains required for reactive oxygen species generation through the Rac1-RbohB pathway, which confers rice immunity against Magnaporthe oryzae [68]. Intriguingly, two fatty acid elongase genes (Os02g0205500, Os05g0568000) were induced in the infested sheaths (Table S7c). In the endoplasmic reticulum, the fatty acid elongase complex catalyzes the elongation of C16 and C18 fatty acids to very-long-chain fatty acids (VLCFAs) with C26 to C34 chains [69, 70]. VLCFAs are converted to cuticular wax that forms the outermost hydrophobic layer of aerial organs for the protection of plants against abiotic and biotic stresses [71]. The cuticular wax produced in the endoplasmic reticulum needs to be transported to the cuticle [71]. Thus, increased expression of five genes encoding lipid transfer proteins (Os03g0718800, Os12g0114500, Os06g0643500, Os07g0175600, Os10g0505500) was identified (Table S7c).
Conclusions
Dual resistance of rice to chewing and piercing-sucking insect pests is a rarely observed phenotype due to the involvement of different defense mechanisms in each type of resistance and the antagonistic crosstalk between these mechanisms. Here, Qingliu was reported to display dual resistance to leaffolder (chewing insect) and BPH (piercing-sucking insect). It is vital to discover the genetic loci in Qingliu that contribute to the dual insect resistance phenotype as this information will be useful for breeding programs aiming to generate broad spectrum insect-resistant rice cultivars. To achieve this objective, tissue-specific RNA-Seq profiling was carried out to decipher the key defense mechanisms of Qingliu under leaffolder and BPH infestation respectively. As a result, transcriptional activation of JA biosynthesis and the interconnected shikimate and phenylpropanoid pathways was found in the leaffolder-infested leaves of Qingliu, which may increase the production of various important defense-related secondary metabolites, including SA, melatonin, flavonoids and lignin. On the other hand, transcriptional activation of SA signaling and cellular signaling transduction was found in the BPH-infested sheaths of Qingliu, which may lead to cell wall modification and formation of cuticular wax that perturb BPH feeding. Although the defense responses of Qingliu to both types of insect pests seem largely unrelated, the phenylpropanoid pathway (or more specifically PAL) could be a possible upstream convergent route conferring dual resistance, which warrant further studies.
Methods
Plant materials
Seeds of the Qingliu rice variety were obtained from the Taichung District Agricultural Research and Extension Station, COA, Changhua, Taiwan; seeds of TN1 were obtained from the Department of Agronomy, National Taiwan University, and seeds of Mudgo, H105, TNG67 and Baiqiaowan were obtained from the Chiayi Agricultural Experiment Station, Taiwan Agricultural Research Institute. Qingliu, an indica rice variety, is resistant to C. medinalis [24]. The indica rice cultivar TN1 is highly susceptible to C. medinalis and N. lugens and was used as a susceptibility check [26, 72]. Rice seeds of Qingliu and TN1 were sterilized for 20 mins, rinsed in sterile water and sown on water-moistened filter paper in Petri dishes. After 2 days of incubation at 37 °C, germinated seeds were transferred to 600 ml beakers containing 1X Kimura B nutrient solution [73]. Seedlings were grown in a growth chamber set at 30/25 °C (day/night) with a 12/12 h (day/night) photoperiod. The 1X Kimura B nutrient solution was changed every 3 days. After 8–10 days (at approximately the two-leaf stage), seedlings of uniform size were transplanted individually into a plastic pot (64 mm diameter at the base, 95 mm diameter at the aperture, and height of 165 mm) containing rice paddy soil from Taoyuan, Taiwan, and grown in environmentally controlled growth chambers. Potted seedlings were treated with soluble fertilizer (33, 31, and 27 kg ha− 1 NPK, respectively) and watered every 2 days. Plants were used for experiments 30 days after germination.
Insect colonies
Cnaphalocrocis medinalis
The C. medinalis colony was originally collected from the Taichung District Agricultural Research and Extension Station, COA, Changhua, Taiwan, and reared on maize seedlings (White Pearl, Known-You Seed Co., Taiwan) in mesh cages (BugDorm-4, MegaView, Taiwan) as described by Guo et al. 2019 [26]. The 3rd instar larvae were used in the experiments.
Nilaparvata lugens
The N. lugens biotype 1 colony was obtained from the Chiayi Agricultural Experiment Station, Taiwan Agricultural Research Institute, COA, Chiayi, Taiwan. N. lugens were reared on TN1 seedlings in mesh cages (BugDorm-4, MegaView, Taiwan) in a growth chamber set at 30/25 °C (day/night) with a 12/12 h (day/night) photoperiod.
Standard seedbox screening test (SSST)
A SSST was used to evaluate the resistance of Qingliu to N. lugens [74]. Four varieties were used as BPH-resistant and BPH-susceptible checks in the SSST and the following modified seedbox screening test (MSST): Mudgo and H105 were BPH-moderately resistant checks [75]; TN1 and TNG67 were BPH-susceptible checks [72, 76], in addition to Baiqiaowan (a leaffolder-resistant check) [24]. A total of twenty-four seeds of each rice variety were sown in lines. Only twenty seedlings from each variety were selected to perform the SSST. At 14 days after sowing, seedlings were infested with 2nd- to 3rd-instar N. lugens nymphs at a density of 10 nymphs per seedling. The damage rating used was based on the standard evaluation system used for N. lugens infestation in rice, where 0 = no injury, 1 = slightly damaged, 3 = 1st and 2nd leaves of plants show yellowing symptoms, 5 = 10 to 25% of plants show pronounced yellowing and stunting or wilting symptoms, 7 = more than 50% of plants wilted, and 9 = all plants wilted or dead [77]. The resistance of rice lines was also determined based on the 0–9 scale where 0–3 was classified as resistant, 4–6 as moderately resistant and 7–9 as susceptible [74, 77]. This SSST experiment was repeated three times.
Modified seedbox screening test (MSST)
The MSST experimental setup was modified from the standard evaluation system of the IRRI [77]. Briefly, six rice varieties were planted in a square plastic box (1 m × 1 m). In each box, each variety had four hills with 3–5 seedlings. At the tillering stage (ca. 30–35 days after sowing), rice plants were infested with gravid female N. lugens adults at a density of 0.5–1 adults per plant. When the susceptible checks (TN1) were moderately wilted due to N. lugens infestation, the condition of each seedling was scored using the same standard evaluation system as that used in the SSST every 2 to 3 days. This MSST experiment was repeated three times.
Settlement preference test
Five germinated seeds of Qingliu and five germinated seeds of TN1 were transplanted to each plastic pot containing 1X Kimura B nutrient solution. Seeds of both varieties were transplanted on a metal platform in an alternating arrangement along the circumference of a 55 mm filter paper (ADVANTEC) with a 17 mm linear distance apart from one another. At 14 days after transplanting, 100 third instar N. lugens nymphs were released onto the filter paper of each hydroponic pot (a total of four pots). The pots were enclosed with plastic covers containing mesh cloth windows to prevent insects from escaping. The number of N. lugens nymphs settling on each rice variety was recorded at 3, 6, 24 and 48 h of infestation (N = 400).
Honeydew excretion experiment
The honeydew excretion experiment using treated bromo-cresol green filter paper was adapted from Jena et al. (2017) [78]. Briefly, the 55 mm filter paper (ADVANTEC) was treated with 0.1% bromocresol green (Alfa Aesar) for 5 min. The tillers of rice plants except the major stem were removed before the experiment. A cardboard with a hole in the middle to accommodate the main stem was placed on each 30-day-old potted plant. Then, the dried bromo-cresol green filter paper was placed on the cardboard and covered by a plastic cup (22 mm diameter at the base, 45 mm diameter at the aperture, and height of 45 mm) (upside down) with a cotton plug to prevent N. lugens escape. One gravid female N. lugens starved for 1 h was released on each plant and allowed to feed. After feeding for 24 h, the filter paper was collected, scanned and analyzed using ImageJ [79]. Due to the chemical features of bromo-cresol green, blue-rimmed spots represent phloem-based honeydew, whereas white spots represent xylem-based honeydew. This experiment was repeated three times with a sample size of 8–10.
Plant treatment, RNA isolation and sequencing library preparation
Each 30-day-old Qingliu rice plant was infested with fifteen 3rd instar BPH nymphs or one 3rd instar leaffolder larva for 24 h. The uninfested plants were used as controls. Each treatment (BPH infestation, leaffolder infestation and uninfested control) had 15 individual plants. Leaf and sheath tissues of each plant were collected separately and immediately frozen in liquid nitrogen and stored at − 80 °C. Total RNA was extracted with TRIzol reagent (Ambion) and purified using Quick-RNA Miniprep (Zymo Research) for each sample according to the manufacturer’s instructions. RNA integrity and quality were checked with a NanoDrop ND 1000 Spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA). The total RNA of three individual plant tissues from the same treatment was pooled as one replicate. Each treatment had five biological replicates. Sequencing libraries were generated using the KAPA mRNA HyperPrep Kit (KAPA Biosystems, Roche, Basel, Switzerland) following the manufacturer’s recommendations. Library preparation and high-throughput RNA sequencing were performed by the Biotools Company (New Taipei City, Taiwan).
Bioinformatics analysis
The original data obtained by high-throughput sequencing (Illumina NovaSeq 6000 platform) were transformed into raw sequenced reads by CASAVA base calling and stored in FASTQ format. FastQC and MultiQC were used to check the quality of fastq files [80]. The obtained raw paired-end reads were filtered by Trimmomatic (v0.38) to discard low-quality reads and trim adaptor sequences and eliminate poor-quality bases with the following parameters: LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:30 [81]. The obtained high-quality data (clean reads) were used for subsequent analysis. Read pairs from each sample were aligned to the rice reference genome (Os-Nipponbare-Reference-IRGSP-1.0) by HISAT2 software (v2.1.0) [82, 83]. FeatureCounts (v1.6.0) was used to count the read numbers mapped to individual genes [84]. For gene expression, “relative log expression” normalization (RLE) was performed using DESeq2 (v1.22.1) [85, 86]. DEG identification analysis between two conditions was performed in R using DESeq2. The resulting p-values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate (FDR) [87]. Genes that had an absolute log2-fold change higher than one (LFC1) and an adjusted p-value less than 0.05 were considered DEGs. GO enrichment and pathway analysis of DEGs were conducted with agriGO (v2.0) and MapMan (v3.6.0) with default settings [88, 89]. Gene coexpression analysis was performed using the RiceFREND server [90]. Phytohormone-responsive genes were extracted from Garg et al. (2012) [29], and heat maps to illustrate the expression change patterns were generated using the pheatmap package in R [91].
Real-time RT-PCR
Total RNA samples from another independent experimental batch were used for qRT-PCR. For each sample, 1 μg of total RNA was reverse transcribed using the iScript cDNA Synthesis Kit (Bio-Rad, USA). Each qPCR consisted of 3 μL PCR-grade water, 5 μL iQ SYBR Green Supermix (Bio-Rad), 0.5 μL forward primer, 0.5 μL reverse primer and 1 μL cDNA template. qPCR was conducted on a CFX Connect Real-Time PCR Detection System (Bio-Rad) following the qPCR kit protocol. Melting curves were generated to confirm primer specificity. Relative expression of a gene of interest [ΔCt; ΔCt = Ct (gene of interest)-Ct (reference gene)] in each sample was measured in four biological replicates and three technical repeats. Primer details for nine genes of interest and a reference gene ubiquitin (Os03g0234200) are shown in Table S8. The reference gene ubiquitin (Os03g0234200) has been used as control for qPCR analysis in previous studies [92,93,94].
Data analyses
Statistical analysis was carried out using the z-test or analysis of variance (ANOVA) as appropriate, and least significant difference (LSD) was used to test for differences between samples with alpha = 0.05. Data were analyzed using the free statistical software R (Version 3.5.1) [95].
Availability of data and materials
The RNA-Seq datasets generated and/or analyzed during the current study are available in the NCBI Sequence Read Archive (SRA) under accession number PRJNA689251 [https://www.ncbi.nlm.nih.gov/bioproject/PRJNA689251]. The reviewer link to the RNA-Seq datasets is provided [https://dataview.ncbi.nlm.nih.gov/object/PRJNA689251?reviewer=2b89mnmtv4988j74d7umd9hpq8] and the datasets will be released once the reviewers and Editor recommend acceptance, i.e. before the manuscript is accepted for publication. All data generated or analyzed during this study are included in this published article and its supplementary information files.
Abbreviations
- SSST:
-
Standard seedbox screening test
- MSST:
-
Modified seedbox screening test
- BPH:
-
Brown planthopper
- DEG:
-
Differentially expressed gene
- RNA-Seq:
-
RNA sequencing
- IPM:
-
Integrated pest management
- JA:
-
Jasmonic acid
- SA:
-
Salicylic acid
- ACC:
-
Ethylene
- QL:
-
Qingliu
- TN1:
-
Taichung native 1
- PCA:
-
Principal component analysis
- qRT-PCR:
-
Real-time reverse transcription polymerase chain reaction
- GO:
-
Gene ontology
- KEGG:
-
Kyoto encyclopedia of genes and genomes
- LOX:
-
Lipoxygenase
- COI:
-
Coronatine-insensitive
- AIM:
-
Abnormal inflorescence meristem
- NPR:
-
Non expressor of pathogenesis-related
- JAZ:
-
Jasmonate ZIM-domain
- PAL:
-
Phenylalanine ammonia-lyase
- CHS:
-
Chalcone synthase
- CHI:
-
Chalcone isomerase
- CC-NBS-LRR:
-
Coiled-coil, nucleotide-binding site and leucine-rich repeat
- CWME:
-
Cell-wall modifying enzyme
- DAMP:
-
Damage-associated molecular pattern
- MAPK:
-
Mitogen-activated protein kinase
- VLCFA:
-
Very-long-chain fatty acid
- FAH:
-
Fatty acid 2-hydroxylase
References
Datta A, Ullah H, Ferdous Z. Water management in rice. In: Chauhan BS, Jabran K, Mahajan G, editors. Rice production worldwide. Switzerland: Springer; 2017. p. 255–78. https://doi.org/10.1007/978-3-319-47516-5_11.
Hedge S, Hedge V. Assessment of global rice production and export opportunity for economic development in Ethiopia. Int J Sci Res. 2013;2(6):257–60.
Ali MP, Bari MN, Haque SS, Kabir MMM, Afrin S, Nowrin F, et al. Establishing next-generation pest control services in rice fields: eco-agriculture. Sci Rep. 2019;9(1):10180. https://doi.org/10.1038/s41598-019-46688-6.
Pathak MD, Khan ZR. Insect pests of rice. Philippines: International Rice Research Institute; 1994.
Heong KL, Wong L, Reyes JHD. Addressing planthopper threats to Asian rice farming and food security: fixing insecticide misuse. Philippines: Asian Development Bank; 2013.
Child RE. Insect damage as a function of climate. In: Padfield T, Borchersen K, editors. Museum microclimates. Denmark: National Museum of Denmark; 2007. p. 57–60.
Deutsch CA, Tewksbury JJ, Tigchelaar M, Battisti DS, Merrill SC, Huey RB, et al. Increase in crop losses to insect pests in a warming climate. Science. 2018;361(6405):916–9. https://doi.org/10.1126/science.aat3466.
Ali MP, Bari MN, Ahmed N, Kabir MMM, Afrin S, Zaman MAU, et al. Rice production without insecticide in smallholder farmer's field. Front Environ Sci. 2017;5:16. https://doi.org/10.3389/fenvs.2017.00016.
Xu H, Yang Y, Lu Y, Zheng X, Tian J, Lai F, et al. Sustainable management of rice insect pests by non-chemical-insecticide technologies in China. Rice Sci. 2017;24(2):61–72.
Ranga Rao G, Kumari B, Sahrawat K, Wani S. Integrated pest management (IPM) for reducing pesticide residues in crops and natural resources. In: Chakravarthy A, editor. New horizons in insect science: towards sustainable pest management. New Delhi: Springer; 2015. p. 397–412.
Van Nguyen N, Ferrero A. Meeting the challenges of global rice production. Paddy Water Environ. 2006;4(1):1–9. https://doi.org/10.1007/s10333-005-0031-5.
Du B, Chen R, Guo J, He G. Current understanding of the genomic, genetic, and molecular control of insect resistance in rice. Mol Breeding. 2020;40(2):24. https://doi.org/10.1007/s11032-020-1103-3.
Fujita D, Kohli A, Horgan FG. Rice resistance to planthoppers and leafhoppers. Crit Rev Plant Sci. 2013;32(3):162–91. https://doi.org/10.1080/07352689.2012.735986.
Al Atalah B, Smagghe G, Van Damme EJ. Orysata, a jacalin-related lectin from rice, could protect plants against biting-chewing and piercing-sucking insects. Plant Sci. 2014;221–222:21–8.
Walling LL. The myriad plant responses to herbivores. J Plant Growth Regul. 2000;19(2):195–216. https://doi.org/10.1007/s003440000026.
Zhou Z, Zang Y, Li Y, Zhang Y, Wang P, Luo X. Rice plant-hopper infestation detection and classification algorithms based on fractal dimension values and fuzzy C-means. Math Comput Model. 2013;58(3–4):701–9. https://doi.org/10.1016/j.mcm.2011.10.028.
Nault LR, Ammar ED. Leafhopper and planthopper transmission of plant viruses. Annu Rev Entomol. 1989;34(1):503–29. https://doi.org/10.1146/annurev.en.34.010189.002443.
Du B, Zhang W, Liu B, Hu J, Wei Z, Shi Z, et al. Identification and characterization of Bph14, a gene conferring resistance to brown planthopper in rice. Proc Natl Acad Sci U S A. 2009;106(52):22163–8. https://doi.org/10.1073/pnas.0912139106.
Li C, Luo C, Zhou Z, Wang R, Ling F, Xiao L, et al. Gene expression and plant hormone levels in two contrasting rice genotypes responding to brown planthopper infestation. BMC Plant Biol. 2017;17(1):57. https://doi.org/10.1186/s12870-017-1005-7.
Guo HM, Li HC, Zhou SR, Xue HW, Miao XX. Cis-12-oxo-phytodienoic acid stimulates rice defense response to a piercing-sucking insect. Mol Plant. 2014;7(11):1683–92. https://doi.org/10.1093/mp/ssu098.
Zhou G, Qi J, Ren N, Cheng J, Erb M, Mao B, et al. Silencing OsHI-LOX makes rice more susceptible to chewing herbivores, but enhances resistance to a phloem feeder. Plant J. 2009;60(4):638–48. https://doi.org/10.1111/j.1365-313X.2009.03988.x.
Ye M, Luo SM, Xie JF, Li YF, Xu T, Liu Y, et al. Silencing COI1 in rice increases susceptibility to chewing insects and impairs inducible defense. PLoS One. 2012;7(4):e36214. https://doi.org/10.1371/journal.pone.0036214.
Spoel SH, Koornneef A, Claessens SMC, Korzelius JP, Van Pelt JA, Mueller MJ, et al. NPR1 modulates cross-talk between salicylate- and jasmonate-dependent defense pathways through a novel function in the cytosol. Plant Cell. 2003;15(3):760–70. https://doi.org/10.1105/tpc.009159.
Liao CT, Chen CL. Oviposition preference and larval performance of Cnaphalocrocis medinalis (Lepidoptera: Pyralidae) on rice genotypes. J Econ Entomol. 2017;110(3):1291–7. https://doi.org/10.1093/jee/tox071.
Tu KY, Tsai SF, Guo TW, Lin HH, Yang ZW, Liao CT, et al. The role of plant abiotic factors on the interactions between Cnaphalocrocis medinalis (Lepidoptera: Crambidae) and its host plant. Environ Entomol. 2018;47(4):857–66. https://doi.org/10.1093/ee/nvy066.
Guo TW, Liao CT, Chuang WP. Defensive responses of rice cultivars resistant to Cnaphalocrocis medinalis (Lepidoptera: Crambidae). Arthropod Plant Interact. 2019;13(4):611–20. https://doi.org/10.1007/s11829-019-09673-4.
Tsai CH. Writing and displaying agricultural knowledge: the case study of rice breeds transition of northern Taiwan in 19th century. Bull Taiwan Hist Res. 2013;6:3–51.
Cheah BH, Lin HH, Chien HJ, Liao CT, Liu LYD, Lai CC, et al. SWATH-MS-based quantitative proteomics reveals a uniquely intricate defense response in Cnaphalocrocis medinalis-resistant rice. Sci Rep. 2020;10(1):6597. https://doi.org/10.1038/s41598-020-63470-1.
Garg R, Tyagi AK, Jain M. Microarray analysis reveals overlapping and specific transcriptional responses to different plant hormones in rice. Plant Signal Behav. 2012;7(8):951–6. https://doi.org/10.4161/psb.20910.
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. https://doi.org/10.1093/mp/ssu049.
Zhang Y, Goritschnig S, Dong X, Li X. A gain-of-function mutation in a plant disease resistance gene leads to constitutive activation of downstream signal transduction pathways in suppressor of npr1-1, constitutive 1. Plant Cell. 2003;15(11):2636–46. https://doi.org/10.1105/tpc.015842.
Yadav V, Wang Z, Wei C, Amo A, Ahmed B, Yang X, et al. Phenylpropanoid pathway engineering: an emerging approach towards plant defense. Pathogens. 2020;9(4):312. https://doi.org/10.3390/pathogens9040312.
Tianpei X, Li D, Qiu P, Luo J, Zhu Y, Li S. Scorpion peptide LqhIT2 activates phenylpropanoid pathways via jasmonate to increase rice resistance to rice leafrollers. Plant Sci. 2015;230:1–11. https://doi.org/10.1016/j.plantsci.2014.10.005.
Tohge T, Watanabe M, Hoefgen R, Fernie AR. Shikimate and phenylalanine biosynthesis in the green lineage. Front Plant Sci. 2013;4:62.
Barros J, Dixon RA. Plant phenylalanine/tyrosine ammonia-lyases. Trends Plant Sci. 2020;25(1):66–79. https://doi.org/10.1016/j.tplants.2019.09.011.
Lu H, Luo T, Fu H, Wang L, Tan Y, Huang J, et al. Resistance of rice to insect pests mediated by suppression of serotonin biosynthesis. Nat Plants. 2018;4(6):338–44. https://doi.org/10.1038/s41477-018-0152-7.
Liang D, Shen Y, Ni Z, Wang Q, Lei Z, Xu N, et al. Exogenous melatonin application delays senescence of kiwifruit leaves by regulating the antioxidant capacity and biosynthesis of flavonoids. Front Plant Sci. 2018;9:426. https://doi.org/10.3389/fpls.2018.00426.
War AR, Paulraj MG, Ahmad T, Buhroo AA, Hussain B, Ignacimuthu S, et al. Mechanisms of plant defense against insect herbivores. Plant Signal Behav. 2012;7(10):1306–20. https://doi.org/10.4161/psb.21663.
Barros J, Serrani-Yarce JC, Chen F, Baxter D, Venables BJ, Dixon RA. Role of bifunctional ammonia-lyase in grass cell wall biosynthesis. Nat Plants. 2016;2(6):16050. https://doi.org/10.1038/nplants.2016.50.
Camm EL, Towers GHN. Phenylalanine ammonia lyase. Phytochemistry. 1973;12(5):961–73. https://doi.org/10.1016/0031-9422(73)85001-0.
Hasegawa S, Maier VP. Cinnamate hydroxylation and the enzymes leading from phenylpyruvate to p-coumarate synthesis in grapefruit tissues. Phytochemistry. 1972;11(4):1365–70. https://doi.org/10.1016/S0031-9422(00)90089-X.
Barros J, Serk H, Granlund I, Pesquet E. The cell biology of lignification in higher plants. Ann Bot. 2015;115(7):1053–74. https://doi.org/10.1093/aob/mcv046.
Liu Q, Luo L, Zheng L. Lignins: biosynthesis and biological functions in plants. Int J Mol Sci. 2018;19(2):335. https://doi.org/10.3390/ijms19020335.
Santiago R, Barros-Rios J, Malvar RA. Impact of cell wall composition on maize resistance to pests and diseases. Int J Mol Sci. 2013;14(4):6960–80. https://doi.org/10.3390/ijms14046960.
Alejandro S, Lee Y, Tohge T, Sudre D, Osorio S, Park J, et al. AtABCG29 is a monolignol transporter involved in lignin biosynthesis. Curr Biol. 2012;22(13):1207–12. https://doi.org/10.1016/j.cub.2012.04.064.
Lefevere H, Bauters L, Gheysen G. Salicylic acid biosynthesis in plants. Front Plant Sci. 2020;11:338. https://doi.org/10.3389/fpls.2020.00338.
Wang Y, Cao L, Zhang Y, Cao C, Liu F, Huang F, et al. Map-based cloning and characterization of BPH29, a B3 domain-containing recessive gene conferring brown planthopper resistance in rice. J Exp Bot. 2015;66(19):6035–45. https://doi.org/10.1093/jxb/erv318.
He J, Liu Y, Yuan D, Duan M, Liu Y, Shen Z, et al. An R2R3 MYB transcription factor confers brown planthopper resistance by regulating the phenylalanine ammonia-lyase pathway in rice. Proc Natl Acad Sci U S A. 2020;117(1):271–7. https://doi.org/10.1073/pnas.1902771116.
Howe GA, Jander G. Plant immunity to insect herbivores. Annu Rev Plant Biol. 2008;59(1):41–66. https://doi.org/10.1146/annurev.arplant.59.032607.092825.
Tamura Y, Hattori M, Yoshioka H, Yoshioka M, Takahashi A, Wu J, et al. Map-based cloning and characterization of a brown planthopper resistance gene BPH26 from Oryza sativa L. ssp. indica cultivar ADR52. Sci Rep. 2014;4(1):5872.
Liu Y, Wu H, Chen H, Liu Y, He J, Kang H, et al. A gene cluster encoding lectin receptor kinases confers broad-spectrum and durable insect resistance in rice. Nat Biotechnol. 2015;33(3):301–5. https://doi.org/10.1038/nbt.3069.
Ji H, Kim S, Kim Y, Suh J, Park H, Sreenivasulu N, et al. Map-based cloning and characterization of the BPH18 gene from wild rice conferring resistance to brown planthopper (BPH) insect pest. Sci Rep. 2016;6(1):34376. https://doi.org/10.1038/srep34376.
van der Knaap E, Song W, Ruan D, Sauter M, Ronald PC, Kende H. Expression of a gibberellin-induced leucine-rich repeat receptor-like protein kinase in Deepwater rice and its interaction with kinase-associated protein phosphatase. Plant Physiol. 1999;120(2):559–69. https://doi.org/10.1104/pp.120.2.559.
Giri J, Vij S, Dansana PK, Tyagi AK. Rice A20/AN1 zinc-finger containing stress-associated proteins (SAP1/11) and a receptor-like cytoplasmic kinase (OsRLCK253) interact via A20 zinc-finger and confer abiotic stress tolerance in transgenic Arabidopsis plants. New Phytol. 2011;191(3):721–32. https://doi.org/10.1111/j.1469-8137.2011.03740.x.
Wan J, Tanaka K, Zhang X-C, Son GH, Brechenmacher L, Nguyen THN, et al. LYK4, a lysin motif receptor-like kinase, is important for chitin signaling and plant innate immunity in Arabidopsis. Plant Physiol. 2012;160(1):396–406. https://doi.org/10.1104/pp.112.201699.
Dubouzet JG, Sakuma Y, Ito Y, Kasuga M, Dubouzet EG, Miura S, et al. OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt- and cold-responsive gene expression. Plant J. 2003;33(4):751–63. https://doi.org/10.1046/j.1365-313X.2003.01661.x.
Seo JS, Joo J, Kim MJ, Kim YK, Nahm BH, Song SI, et al. OsbHLH148, a basic helix-loop-helix protein, interacts with OsJAZ proteins in a jasmonate signaling pathway leading to drought tolerance in rice. Plant J. 2011;65(6):907–21. https://doi.org/10.1111/j.1365-313X.2010.04477.x.
Vannini C, Iriti M, Bracale M, Locatelli F, Faoro F, Croce P, et al. The ectopic expression of the rice Osmyb4 gene in Arabidopsis increases tolerance to abiotic, environmental and biotic stresses. Physiol Mol Plant Pathol. 2006;69(1-3):26–42. https://doi.org/10.1016/j.pmpp.2006.12.005.
Yokotani N, Sato Y, Tanabe S, Chujo T, Shimizu T, Okada K, et al. WRKY76 is a rice transcriptional repressor playing opposite roles in blast disease resistance and cold stress tolerance. J Exp Bot. 2013;64(16):5085–97. https://doi.org/10.1093/jxb/ert298.
Silva-Sanzana C, Estevez JM, Blanco-Herrera F. Influence of cell wall polymers and their modifying enzymes during plant–aphid interactions. J Exp Bot. 2020;71(13):3854–64. https://doi.org/10.1093/jxb/erz550.
van Bel AJE, Will T. Functional evaluation of proteins in watery and gel saliva of aphids. Front Plant Sci. 2016;7:1840.
Kubicek CP, Starr TL, Glass NL. Plant cell wall–degrading enzymes and their secretion in plant-pathogenic fungi. Annu Rev Phytopathol. 2014;52(1):427–51. https://doi.org/10.1146/annurev-phyto-102313-045831.
Bacete L, Mélida H, Miedes E, Molina A. Plant cell wall-mediated immunity: cell wall changes trigger disease resistance responses. Plant J. 2018;93(4):614–36. https://doi.org/10.1111/tpj.13807.
Miedes E, Lorences EP. The implication of xyloglucan endotransglucosylase/hydrolase (XTHs) in tomato fruit infection by Penicillium expansum link. A J Agric Food Chem. 2007;55(22):9021–6. https://doi.org/10.1021/jf0718244.
Claverie J, Balacey S, Lemaître-Guillier C, Brulé D, Chiltz A, Granet L, et al. The cell wall-derived xyloglucan is a new DAMP triggering plant immunity in Vitis vinifera and Arabidopsis thaliana. Front Plant Sci. 2018;9:1725. https://doi.org/10.3389/fpls.2018.01725.
Hao P, Liu C, Wang Y, Chen R, Tang M, Du B, et al. Herbivore-induced callose deposition on the sieve plates of rice: an important mechanism for host resistance. Plant Physiol. 2008;146(4):1810–20. https://doi.org/10.1104/pp.107.111484.
Singh A, Pandey A, Baranwal V, Kapoor S, Pandey GK. Comprehensive expression analysis of rice phospholipase D gene family during abiotic stresses and development. Plant Signal Behav. 2012;7(7):847–55. https://doi.org/10.4161/psb.20385.
Nagano M, Ishikawa T, Fujiwara M, Fukao Y, Kawano Y, Kawai-Yamada M, et al. Plasma membrane microdomains are essential for Rac1-RbohB/H-mediated immunity in rice. Plant Cell. 2016;28(8):1966–83. https://doi.org/10.1105/tpc.16.00201.
Todd J, Post-Beittenmiller D, Jaworski JG. KCS1 encodes a fatty acid elongase 3-ketoacyl-CoA synthase affecting wax biosynthesis in Arabidopsis thaliana. Plant J. 1999;17(2):119–30. https://doi.org/10.1046/j.1365-313X.1999.00352.x.
Wang X, Guan Y, Zhang D, Dong X, Tian L, Qu LQ. A β-ketoacyl-CoA synthase is involved in rice leaf cuticular wax synthesis and requires a CER2-LIKE protein as a cofactor. Plant Physiol. 2017;173(2):944–55. https://doi.org/10.1104/pp.16.01527.
Xue D, Zhang X, Lu X, Chen G, Chen Z. Molecular and evolutionary mechanisms of cuticular wax for plant drought tolerance. Front Plant Sci. 2017;8:621. https://doi.org/10.3389/fpls.2017.00621.
Wan P, Zhou R, Nanda S, He J, Yuan S, Wang W, et al. Phenotypic and transcriptomic responses of two Nilaparvata lugens populations to the Mudgo rice containing Bph1. Sci Rep. 2019;9(1):14049. https://doi.org/10.1038/s41598-019-50632-z.
Yoshida S. Routine procedure for growing rice plants in culture solution. In: Yoshida S, Forno DA, Cook JH, Gomez KA, editors. Laboratory manual for physiological studies of rice. Philippines: International Rice Research Institute; 1976. p. 61–6.
Velusamy R, Heinrichs EA, Medrano FG. Greenhouse techniques to identify field resistance to the brown planthopper, Nilaparvata lugens (Stål) (Homoptera: Delphacidae), in rice cultivars. Crop Prot. 1986;5(5):328–33. https://doi.org/10.1016/0261-2194(86)90112-2.
Huang SH, Cheng CH, Chen CN, Wu WJ. Bioassay of the biotypes of brown plant hopper in Taiwan. J Taiwan Agric Res. 2016;65:84–91.
Lin HY, Wu YP, Hour AL, Ho SW, Wei FJ, Hsing Y, et al. Genetic diversity of rice germplasm used in Taiwan breeding programs. Bot Stud. 2012;53:363–76.
IRRI. Standard evaluation system for rice. Philippines: International Rice Research Institute; 2002.
Jena KK, Hechanova SL, Verdeprado H, Prahalada GD, Kim SR. Development of 25 near-isogenic lines (NILs) with ten BPH resistance genes in rice (Oryza sativa L.): production, resistance spectrum, and molecular analysis. Theor Appl Genet. 2017;130(11):2345–60. https://doi.org/10.1007/s00122-017-2963-8.
Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5. https://doi.org/10.1038/nmeth.2089.
Ewels P, Magnusson M, Lundin S, Kaller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8. https://doi.org/10.1093/bioinformatics/btw354.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Kawahara Y, de la Bastide M, Hamilton JP, Kanamori H, McCombie WR, Ouyang S, et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013;6(1):4. https://doi.org/10.1186/1939-8433-6-4.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. https://doi.org/10.1038/nmeth.3317.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30. https://doi.org/10.1093/bioinformatics/btt656.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106. https://doi.org/10.1186/gb-2010-11-10-r106.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.
Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J Royal Statist Soc Series B. 1995;57:289–300.
Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37(6):914–39. https://doi.org/10.1111/j.1365-313X.2004.02016.x.
Tian T, Liu Y, Yan H, You Q, Yi X, Du Z, et al. agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017;45(W1):W122–W9. https://doi.org/10.1093/nar/gkx382.
Sato Y, Namiki N, Takehisa H, Kamatsuki K, Minami H, Ikawa H, et al. RiceFREND: a platform for retrieving coexpressed gene networks in rice. Nucleic Acids Res. 2013;41(Database issue):D1214–21. https://doi.org/10.1093/nar/gks1122.
Kolde R. pheatmap: Pretty heatmap [Software]. 2015.
Dong W, Thomas N, Ronald PC, Goyer A. Overexpression of thiamin biosynthesis genes in rice increases leaf and unpolished grain thiamin content but not resistance to Xanthomonas oryzae pv. oryzae. Front Plant Sci. 2016;7:616.
Wu T, Zhang M, Zhang H, Huang K, Chen M, Chen C, et al. Identification and characterization of EDT1 conferring drought tolerance in rice. J Plant Biol. 2019;62(1):39–47. https://doi.org/10.1007/s12374-018-0203-7.
Wang W, Liu B, Xu M, Jamil M, Wang G. ABA-induced CCCH tandem zinc finger protein OsC3H47 decreases ABA sensitivity and promotes drought tolerance in Oryza sativa. Biochem Biophys Res Commun. 2015;464(1):33–7. https://doi.org/10.1016/j.bbrc.2015.05.087.
Development R. Core team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for statistical. Computing. 2010.
Acknowledgments
We thank Dr. Shu-Jen Wang (National Taiwan University) for generously providing the TN1 seeds. We thank Dr. Men-Chi Chang (National Taiwan University) and Dr. Liang-Jwu Chen (National Chung Hsing University) for advice on data analysis and discussion.
Funding
This work was supported by grants from the Ministry of Science and Technology, Taiwan, Grant Numbers 106–2311-B-002-025 & 107–2311-B-002-018-MY3 (WPC.)
Author information
Authors and Affiliations
Contributions
WPC conceived and designed the experiment; YL, YHK and SCL conducted the experiments; YL, BHC and YFF analyzed the data; BHC, CTL, SHH, YFL and WPC supervised the research and contributed to the discussion of the results; and BHC, YFL and WPC wrote the paper. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no conflicts of interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Fig. S1.
Principal component analysis (PCA) plot of thirty RNA-Seq samples. The samples consisted of three experimental conditions (control, BPH or leaffolder infestation) and two tissues (leaves, sheaths) with five biological replicates for each combination. Fig. S2. Shortlisting nonredundant GO terms of leaffolder-responsive DEGs by manually curating the GO hierarchical tree graphs. GO enrichment of a upregulated genes in the leaffolder vs control_leaf comparison, b downregulated genes in the leaffolder vs control_leaf comparison and c upregulated genes in the leaffolder vs control_sheath comparison. Fig. S3. Shortlisting nonredundant GO terms of BPH-responsive DEGs by manually curating the GO hierarchical tree graphs. GO enrichment of a upregulated genes in the BPH vs control_leaf comparison, b downregulated genes in the BPH vs control_leaf comparison and c upregulated genes in the BPH vs control_sheath comparison.
Additional file 2: Table S1.
Statistics of raw reads and clean reads per RNA-Seq sample mapped to the rice reference genome.
Additional file 3: Table S2.
List of differentially expressed genes. Up- and downregulated genes are shown for a leaffolder vs control_leaf comparison, b leaffolder vs control_sheath comparison, c BPH vs control_leaf comparison, d BPH vs control_sheath comparison, common DEGs in e leaves and f sheaths under separate leaffolder and BPH infestations.
Additional file 4: Table S3.
The expression changes of defensive phytohormone-responsive genes in Qingliu under leaffolder and BPH infestations.
Additional file 5: Table S4.
Details of GO enrichment results for the DEGs in Qingliu at 24 h of leaffolder herbivory. a Upregulated genes in the leaffolder vs control_leaf comparison, b downregulated genes in the leaffolder vs control_leaf comparison, c upregulated genes in the leaffolder vs control_sheath comparison and d downregulated genes in the leaffolder vs control_sheath comparison.
Additional file 6: Table S5.
Details of GO enrichment results for the DEGs in Qingliu at 24 h of BPH infestation. a Upregulated genes in the BPH vs control_sheath comparison, b downregulated genes in the BPH vs control_sheath comparison, c upregulated genes in the BPH vs control_leaf comparison and d downregulated genes in the BPH vs control_leaf comparison.
Additional file 7: Table S6.
MapMan pathway annotation results for the DEGs in the BPH-infested sheaths of Qingliu at 24 h. The mapping results of a receptor kinases, b biotic stress and c metabolism overview are shown.
Additional file 8: Table S7.
Merged GO and MapMan results for the DEGs in the BPH-infested sheaths of Qingliu at 24 h. a transcription factor, b cell wall modification and c lipid metabolism.
Additional file 9: Table S8.
Primers used for qRT-PCR amplifications.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Li, Y., Cheah, B.H., Fang, YF. et al. Transcriptomics identifies key defense mechanisms in rice resistant to both leaf-feeding and phloem feeding herbivores. BMC Plant Biol 21, 306 (2021). https://doi.org/10.1186/s12870-021-03068-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-021-03068-5