Next Article in Journal
Next-Generation Sequencing Workflow for NSCLC Critical Samples Using a Targeted Sequencing Approach by Ion Torrent PGM™ Platform
Next Article in Special Issue
Phenolic Phytoalexins in Rice: Biological Functions and Biosynthesis
Previous Article in Journal
Current Status of Herbal Medicines in Chronic Liver Disease Therapy: The Biological Effects, Molecular Targets and Future Prospects
Previous Article in Special Issue
Carbon Fluxes between Primary Metabolism and Phenolic Pathway in Plant Tissues under Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Omics-Based Comparative Transcriptional Profiling of Two Contrasting Rice Genotypes during Early Infestation by Small Brown Planthopper

1
College of Chemistry and Life Sciences, Zhejiang Normal University, Jinhua 321004, China
2
Analysis Center of Agrobiology and Environmental Sciences, Zhejiang University, Hangzhou 310058, China
3
State Key Laboratory Breeding Base for Zhejiang Sustainable Pest and Disease Control, Ministry of China Key Laboratory of Biotechnology in Plant Protection, Institute of Virology and Biotechnology, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2015, 16(12), 28746-28764; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms161226128
Submission received: 20 September 2015 / Revised: 2 November 2015 / Accepted: 4 November 2015 / Published: 3 December 2015
(This article belongs to the Special Issue Molecular Research in Plant Secondary Metabolism 2015)

Abstract

:
The small brown planthopper (SBPH) is one of the destructive pests of rice. Although different biochemical pathways that are involved in rice responding to planthopper infestation have been documented, it is unclear which individual metabolic pathways are responsive to planthopper infestation. In this study, an omics-based comparative transcriptional profiling of two contrasting rice genotypes, an SBPH-resistant and an SBPH-susceptible rice line, was assessed for rice individual metabolic pathways responsive to SBPH infestation. When exposed to SBPH, 166 metabolic pathways were differentially regulated; of these, more than one-third of metabolic pathways displayed similar change patterns between these two contrasting rice genotypes; the difference of change pattern between these two contrasting rice genotypes mostly lies in biosynthetic pathways and the obvious difference of change pattern lies in energy metabolism pathways. Combining the Pathway Tools Omics Viewer with the web tool Venn, 21 and 6 metabolic pathways which potentially associated with SBPH resistance and susceptibility, respectively were identified. This study presents an omics-based comparative transcriptional profiling of SBPH-resistant and SBPH-susceptible rice plants during early infestation by SBPH, which will be very informative in studying rice-insect interaction. The results will provide insight into how rice plants respond to early infestation by SBPH from the biochemical pathways perspective.

1. Introduction

Rice (Oryza sativa L.) is an important food crop and a source of calories for billions of people [1]. Many insects attack rice. Of these insects, the brown planthopper (BPH; Nilaparvata lugens Stål), the small brown planthopper (SBPH; Laodelphax striatellus Fallén) and the white-backed planthopper (WBPH; Sogatella furcifera Horváth) are three major planthoppers attacking rice, causing significant loss of rice yield [2]. Understanding the defense responses of rice to planthopper infestation is important for developing strategies for such insect control. Studies of the molecular responses of rice to these species of sucking insects have shown that rice responses to these species of sucking insects are regulated through a complex network of salicylic acid (SA)- and jasmonic acid (JA)/ethylene (ET)-dependent signaling pathways [3,4,5,6,7,8]. Gene expression profiles and analyses of physiological responses suggested that the genes involved in several different pathways including cell defense, cellular transport, metabolism, signal transduction, macromolecular degradation and biogenesis, cellular communication and plant defenses were found to be significantly differentially regulated by BPH/SBPH feeding [3,7,9,10,11,12]. Transcription factors (TFs) play a key role in regulation of transcriptional expression in biological processes [13]. Microarray analysis suggested that rice TFs were involved in defense response to BPH attack [8,14] and over-expression of rice TFs in transgenic rice plants increased resistance to WBPH [15]. Although gene expression profiles resulting from planthopper infestation have been studied by several groups by using two contrasting rice genotypes, the defense mechanisms of rice against these species of sucking insects remain to be largely explored.
In attempt to gain insight into the molecular mechanisms of rice resistance to SBPH, in this study, two contrasting rice genotypes, an SBPH-resistant introgression line (Pf9279-4; R) and its SBPH-susceptible parent (02428; S), were assessed for rice plant responses to SBPH attack at the molecular level. Previously, two-dimensional fluorescence difference gel electrophoresis analysis of these two contrasting rice genotypes that sampled at different time-points after SBPH infestation showed that the differential protein expression between these two contrasting rice genotypes were obvious within 6 h of infestation by SBPH. We therefore examined the comparative transcriptional profiling of these two contrasting rice genotypes which infested with SBPH for 0 h and 6 h by using Affymetrix GeneChip analysis. Although different biochemical pathways including cell defense, cellular transport, metabolism, signal transduction, macromolecular degradation and biogenesis, cellular communication and plant defenses have been documented [3,7,9,10,11,12], it is unclear which individual pathways are involved in the response to planthopper infestation in rice. The increasing number of sequenced plant genomes and advances in bioinformatics strategies offer immense opportunities to study biological processes related to biotic and/or abiotic stress at the cellular and whole-organism level by using a novel systems-level approach [16]. Among the bioinformatics strategies, because of the unique ability to display the differentially regulated individual biochemical pathways with omics data painted onto the pathway, the Pathway Tools Omics Viewer (http://pathway.iplantcollaborative.org/) was ultimately employed to visualize the expression patterns and further identify the individual pathways that are responsive to SBPH feeding in rice.
The Pathway Tools Omics Viewer, an integrated tool in the Gramene “Cyc” Pathways databases (http://www.gramene.org), provides visual analysis of whole-organism datasets and maps gene expression data sets to the metabolic pathways [17]. The ability to visualize the expression patterns facilitates the ability to pinpoint and interpret expression patterns of interest [16]. Of the Gramene “Cyc” Pathways databases, the RiceCyc pathway database now features 306 metabolic pathways; 2103 enzymatic reactions, 87 transport reactions, 1543 compounds and metabolites and 6040 protein coding genes are mapped to these reactions and pathways (http://pathway.iplantcollaborative.org/). In the present study, to understand the global gene expression profiles, the transcriptomic data was uploaded to visualize expression patterns by using the Pathway Tools Omics Viewer. To identify the SBPH resistance and susceptibility-related metabolic pathways, by using Pathway Tools Omics Viewer, the individual metabolic pathways that were differentially regulated were displayed with omics data painted onto the pathway; by combining Pathway Tools Omics Viewer with the web tool Venn (http://bioinfogp.cnb.csic.es/), 21 and 6 metabolic pathways potentially related to SBPH resistance and susceptibility in rice were identified, respectively. This study presents an omics-based comparative transcriptional profiling of SBPH-resistant and SBPH-susceptible rice plants during early infestation by SBPH and aims to identify individual pathways that are responsive to early infestation by SBPH in rice. The identified biochemical pathways will provide insight into how rice plants respond to early infestation by SBPH and the results will help us understand the general defense responses of rice plants to insect infestation from the perspective of biochemical pathways.

2. Results

2.1. The Global Transcriptional Profiling of Rice In Response to Small Brown Planthopper (SBPH) Infestation

According to the results of our previous comparative two-dimensional fluorescence difference gel electrophoresis analysis, significant differences in protein expression level between these two contrasting rice genotypes sampled at different time-points after SBPH infestation were obvious within 6 h of infestation. To gain better understanding on how rice plants respond to SBPH attack at the molecular level, the transcriptional profiling of SBPH-resistant introgression line (Pf9279-4) and its SBPH-susceptible parent (02428) that had been fed upon by SBPH for 0 and 6 h was thusly examined via Affymetrix GeneChip analysis. To analyze the molecular mechanism underlying rice plant SBPH resistance, the expression levels of all rice genes in the resistant and susceptible rice plants that had been fed upon by SBPH for 0 and 6 h were examined and the comparison between before and after infestation was made and named as S6_S0 and R6_R0 in the resistant and susceptible rice plants, respectively. To gain an understanding of the basis of defense, comparisons of gene expression were also made between resistant and susceptible rice plants and named as R0_S0 and R6_S6 before and after infestation, respectively. The microarray experiments identified a total of 57,381 probe sets expressed in the four comparisons. In the comparison of rice plants that had been fed upon by SBPH for 0 and 6 h, a total of 2730 differentially expressed genes were identified in the resistant rice plants (R6_R0), while a total of 3024 differentially expressed genes were identified in the susceptible rice plants (S6_S0). In the resistant and susceptible comparison, a total of 2299 differentially expressed genes were identified when the rice plants had been fed upon by SBPH for 0 h (R0_S0), while a total of 2891 differentially expressed genes were identified when the rice plants had been fed upon by SBPH for 6 h (R6_S6).

2.2. The Global Metabolic Map in Response to SBPH Infestation

As shown in Table S1, in the resistant and susceptible comparison, the number of the genes that appear in the overview of metabolic map is 407 and 552 in R0_S0 and R6_S6, respectively; the ratio of the genes that appear in the overview of the metabolic map to all the differentially expressed genes is about 17.7% and 19.1% in R0_S0 and R6_S6, respectively. In the comparison of rice plants that had been fed upon by SBPH for 6 and 0 h, the number of the genes that appear in the overview of metabolic map is 671 and 707 in R6_R0 and S6_S0, respectively; the ratio of the genes that appear in the overview of the metabolic map to all the differentially expressed genes is about 24.6% and 23.4% in R6_R0 and S6_S0, respectively. The cellular overview diagram was organized as follows. The functionally related pathways are grouped together, with shaded background boxes to delineate these groupings; biosynthetic pathways were positioned on the left, energy metabolism pathways in the middle and degradation pathways on the right [17]. In the overview diagram in this study, red for data values that exceed 2; purple for data values less than the inverse of 2, and black for data values in between; the 2 is the cutoff value of fold change that was employed to discriminate differential expression of genes. As shown in Table S1 (all the categories of the metabolic pathways in the cellular overview diagram are also shown in bold as follows in this paragraph), in response to SBPH infestation, other biosynthesis differentially altered only in the R6_R0; respiration differentially altered in the R6_S6 and S6_S0 as well as R0_S0; inorganic nutrients metabolism differentially altered in both R6_R0 and R0_S0; C1 compounds utilization and assimilation (utilization and assimilation of compounds containing one carbon) differentially altered in the R6_S6 and R6_R0 as well as R0_S0. In the R0_S0, in addition to other biosynthesis, energy metabolism pathways including pentose phosphate pathways and respiration, degradation pathways including hormones degradation, secondary metabolites degradation, carbohydrates degradation and degradation/utilization/assimilation-other were not differentially altered in all four comparisons. The remaining metabolic pathways as follows were differentially altered in all four comparisons. Biosynthetic pathways include nucleosides and nucleotides biosynthesis; cofactors, prosthetic groups, electron carriers biosynthesis; carbohydrates biosynthesis; secondary metabolites biosynthesis; fatty acids and lipids biosynthesis; amino acids biosynthesis; amines and polyamines biosynthesis; hormones biosynthesis and cell structures biosynthesis. Energy metabolism pathways include fermentation; glycolysis; photosynthesis and TCA cycle (tricarboxylicacidcycle acid cycle). Degradation pathways include amino acids degradation; alcohols degradation; detoxification; carbohydrates degradation; fatty acid and lipids degradation as well as nucleosides and nucleotides degradation.

2.3. An Overview of the Differentially Regulated Metabolic Pathways

The Pathway Tools Omics Viewer provides the capability to display the output as an animation for time-series experiments but not for single experiment time step. As shown in Table S1, for one comparative transcriptome analysis of interest, e.g., R6_R0, the metabolic map generated by Pathway Tools Omics Viewer is static but not animated. To analyze the molecular mechanism underlying rice plant SBPH resistance, the metabolic pathways that were differentially regulated between R6_R0 and S6_S0 were examined. To more readily visualize the difference in the differential regulation of metabolic pathways responsive to SBPH infestation between the resistant and susceptible rice plants, the cellular overview diagram of R6_R0 and S6_S0 was then taken as a combination for an animation. As viewed in the Pathway Tools Omics Viewer, where the difference exists between R6_R0 and S6_S0, the line between nodes in the cellular overview diagram will flash in red/purple (Movie S1). With the help of the animated pictures, we accurately and rapidly traced and marked the difference between R6_R0 and S6_S0 and thereafter shown that in one picture. As shown in Table S2, in response to SBPH infestation, the difference between R6_R0 and S6_S0 mostly lies in biosynthetic pathways; while the obvious difference between R6_R0 and S6_S0 lies in energy metabolism pathways; after attack by SBPH, respiration was down-regulated and glycolysis was up-regulated in the susceptible rice plants, whereas, respiration was not differentially altered and glycolysis was partially down-regulated in the resistant rice plants. As for degradation pathways, the major difference between R6_R0 and S6_S0 lies in inorganic nutrients metabolism, which was down-regulated in the resistant rice plants but was not differentially altered in the susceptible rice plants. Another difference between R6_R0 and S6_S0 lies in degradation/utilization/assimilation-other, which almost up-regulated in the resistant rice plants but was not differentially altered in the susceptible rice plants.
To gain an understanding of the difference in the basis of defense in response to SBPH infestation, the R6_S6 and R0_S0 were also taken as the combination for an animation and marked the difference between R6_S6 and R0_S0 and then shown in one picture. As shown in Table S2, if the value of the ratio of the differentially regulated metabolic pathways between R6_S6 and R0_S0 to the all metabolic pathways contained in the cellular overview diagram is taken to discriminate whether the difference between R6_S6 and R0_S0 is significant or not, degradation pathways were the most significant, and energy metabolism pathways were the second most significant. As for biosynthetic pathways, the difference between R6_S6 and R0_S0 lies in carbohydrates biosynthesis, amino acids biosynthesis, amines and polyamines biosynthesis, fatty acids and lipids biosynthesis, hormones biosynthesis and cell structures biosynthesis. As for energy metabolism pathways, fermentation and glycolysis showed a difference between R6_S6 and R0_S0. As for degradation pathways, detoxification, C1 compounds utilization and assimilation, inorganic nutrients metabolism, hormones degradation, carbohydrates degradation, amino acids degradation, degradation/utilization/assimilation-other, fatty acid and lipids degradation displayed difference between R6_S6 and R0_S0.

2.4. The Global Batch Identification of Individual Metabolic Pathways by Pathway Tools Omics Viewer

To further explore the defense responses of rice plants to SBPH attack in more detail, the individual metabolic pathways were afterward screened by Pathway Tools Omics Viewer. Besides the capabilities for overlaying transcriptomic data onto the overview as a whole, the Pathway Tools Omics Viewer can also generate a table of magnified views of all pathways with omics data painted onto the pathway [17]. Therefore, we then chose the alternative display to generate a table containing all individual pathways with omics data painted onto the pathway by Pathway Tools Omics Viewer (Table S3). As shown in Supplement Table S4, in the resistant and susceptible comparison, a total of 82 and 108 metabolic pathways that associated with the differentially expressed genes were identified in R0_S0 and R6_S6, respectively. In the comparison of rice plants that had been fed upon by SBPH for 6 and 0 h, a total of 149 and 128 metabolic pathways that associated with the differentially expressed genes were determined in R6_R0 and S6_S0, respectively. Among the up-regulated metabolic pathways in response to the SBPH infestation in these two contrasting rice genotypes, 26 were shared, 13 and 11 were unique in R6_R0 and S6_S0, respectively; while the number of metabolic pathways that were down-regulated in both rice genotypes was 38, and unique in R6_R0 and S6_S0 were 23 and 5, respectively. In summary, a total of 190 differentially regulated metabolic pathways were determined in the four comparisons by using Pathway Tools Omics Viewer, accounting for nearly two third of the metabolic pathways that now features in the RiceCyc database. As shown in Table 1, in response to SBPH infestation, biosynthesis was the largest class of the metabolic pathways that were differentially regulated in both resistant and susceptible rice plants (R6_R0 and S6_S0) while degradation/utilization/assimilation was the second. Within biosynthesis grouping, secondary metabolites biosynthesis, amino acids biosynthesis and cofactors, prosthetic groups, electron carriers biosynthesis as well as carbohydrates biosynthesis were the largest subclass of the metabolic pathways that were differentially regulated. Within degradation/utilization/assimilation grouping, amino acids degradation and carbohydrates degradation were the largest subclass of the metabolic pathways that were differentially regulated. As for the R6_S6 and R0_S0, the largest class and its subclass of the metabolic pathways that were differentially regulated were the same as that of R6_R0 and S6_S0.
Table 1. Number of the differentially regulated metabolic pathways in response to SBPH attack.
Table 1. Number of the differentially regulated metabolic pathways in response to SBPH attack.
CategoryR6_R0S6_S0R6_S6R0_S0
Activation/Inactivation/Interconversion
Inactivation2 2
activation1
Biosynthesis
Secondary metabolites biosynthesis30201816
Amino acids biosynthesis20191010
Cofactors, prosthetic groups, electron carriers biosynthesis1416159
Carbohydrates biosynthesis1313106
Nucleosides and nucleotides biosynthesis9955
Cell Structures biosynthesis6654
Fatty acids and lipids biosynthesis6877
Hormones biosynthesis5645
Amines and polyamines biosynthesis4222
Aromatic compounds biosynthesis2122
Aminoacyl-tRNA charging 1111
Degradation/Utilization/Assimilation
Amino acids degradation1416158
Carbohydrates degradation109106
Inorganic nutrients metabolism51 1
Fatty acid and lipids degradation4244
Detoxification3312
C1 compounds utilization and assimilation3221
Secondary metabolites degradation3221
Alcohols degradation23
Nucleosides and nucleotides degradation2112
Hormones degradation111
Polymeric compounds degradation111
Amines and polyamines degradation11
Aldehyde degradation 11
Degradation/Utilization/Assimilation-Other 1
Aromatic compounds degradation1
Nicotine degradation1
Generation of Precursor Metabolites and Energy
Fermentation4333
Glycolysis221
Pentose phosphate pathways221
Photosynthesis221
Respiration112
TCA cycle 1122
Acetyl-CoA biosynthesis 1
Methanogenesis 1
C1 compounds utilization and assimilation, utilization and assimilation of compounds containing one carbon; TCA cycle, tricarboxylicacidcycle acid cycle.

2.5. The Metabolic Pathways Potentially Related to SBPH Resistance in Rice

With the help of the Pathway Tools Omics Viewer, the interpretation of large-scale biological data sets can be analyzed and derived from transcriptome profiling experiments. For example, hypotheses about how the phenotypic changes which resulted from different kinds of factors such as stresses based on metabolic explanations can be generated [16,18]. Meanwhile, another web tool Venn (http://bioinfogp.cnb.csic.es/), can calculate the intersection(s) of lists of elements and generate an output displaying which elements are in each intersection or are unique to a certain list [19]. Therefore, based on the analysis by the Pathway Tools Omics Viewer, the web tool Venn was next employed to identify the metabolic pathways potentially related to SBPH resistance. Generally, plant defenses are classified as induced if gene expression changes after herbivore infestation [20,21,22]. Similarly, the metabolic pathways potentially related to SBPH resistance in rice were then examined in the light of the following classification:
Class I: in response to SBPH infestation, the metabolic pathway was differentially regulated in the resistant rice plants (R6_R0) but not in the susceptible rice plants (S6_S0); meanwhile, the basis of defense showed difference between resistant and susceptible rice plants after attack by SBPH (R6_S6) but not before attack by SBPH (R0_S0).
Class II: in response to SBPH infestation, the metabolic pathway was differentially regulated in the resistant rice plants (R6_R0) but not in the susceptible rice plants (S6_S0); meanwhile, the basis of defense showed difference between resistant and susceptible rice plants before attack by SBPH (R0_S0) and after attack by SBPH (R6_S6).
Class III: in response to SBPH infestation, the metabolic pathway was differentially regulated in both resistant rice plants (R6_R0) and susceptible rice plants (S6_S0) and showed opposing differential changes in S6_S0 and R6_R0; meanwhile, the basis of defense showed difference between resistant and susceptible rice plants after attack by SBPH (R6_S6) and or not before attack by SBPH (R0_S0).
As shown in Figure 1 and Table S3, there were 190 metabolic pathways in the 15 categories of the Venn diagram. The number of each category was shown in the Venn diagram. For example, the number of metabolic pathway that was unique in S6_S0, R6_R0, R0_S0 and R6_S6 was 10, 16, 2 and 9, respectively. The individual metabolic pathway of each category in the Venn diagram was shown in Table S3. According to the classification of the metabolic pathways potentially related to SBPH resistance in rice, as shown in Table 2, a total of 8, 11 and 2 metabolic pathways were identified in Classes I to III, respectively (Click the individual metabolic pathway in Supplementary Tables where the individual metabolic pathway mentioned can navigate the individual metabolic pathway webpage).
Figure 1. The Venn diagram showing the differentially regulated metabolic pathways identified by the Omics Viewer. The numbers in the Venn diagram are the number of the differentially regulated metabolic pathway that belongs to each category in the Venn diagram.
Figure 1. The Venn diagram showing the differentially regulated metabolic pathways identified by the Omics Viewer. The numbers in the Venn diagram are the number of the differentially regulated metabolic pathway that belongs to each category in the Venn diagram.
Ijms 16 26128 g001
Table 2. The identified SBPH resistance-related metabolic pathways and its differential regulation in the four comparisons.
Table 2. The identified SBPH resistance-related metabolic pathways and its differential regulation in the four comparisons.
PathwayS6_S0R6_R0R0_S0R6_S6
Class I
Betanidin degradation down lower
Cytokinins degradation down lower
Glutamate degradation III down lower
IAA conjugate biosynthesis I up higher
IAA conjugate biosynthesis II up higher
Spermine biosynthesis up higher
Very long chain fatty acid biosynthesis up higher
Momilactone biosynthesis up higher
Class II
Flavonoid biosynthesis downhigherhigher
Mixed acid fermentation downhigherhigher
Pinobanksin biosynthesis downhigherhigher
Aminopropanol biosynthesis downlowerlower
Salicylate biosynthesis downlowerlower
Serine biosynthesis downlowerlower
Threonine degradation II downlowerlower
Threonine degradation III (to methylglyoxal) downlowerlower
Reductive TCA cycle I uphigherhigher
13-LOX and 13-HPL pathway uplowerlower
Divinyl ether biosynthesis II (13-LOX) uplowerlower
Class III
Ureide biosynthesisupdownlowerlower
phenylalanine degradation IIIupdown lower
IAA, indole-3-acetic acid; 13-LOX, 13-lipoxygenas; 13-HPL, 13-hydroperoxide lyase; TCA cycle, tricarboxylicacidcycle acid cycle.

2.6. Construction the Network of the Metabolic Pathways Potentially Related to SBPH Resistance in the SBPH-Resistant Rice Plant

Based on the classification hierarchy for individual metabolic pathways from the rice pathway database (http://pathway.iplantcollaborative.org/), the network of the metabolic pathways potentially related to SBPH resistance in the SBPH-resistant rice plant was subsequently constructed. As shown in Table 2, Table S5 and Figure 2, the primary metabolic pathways related to fatty acid (very long chain fatty acid) and amino acid (glutamate, serine, threonine and phenylalanine) and the secondary metabolic pathways related to phytoalexin (momilactone A); amine and polyamine (ureide and spermine); alkaloid (betanidin); hormone (jasmonic acid, salicylate, cytokinin and auxin); phenylpropanoid derivatives (flavonoid and pinobanksin) had been exploited potentially related to SBPH resistance in rice. The C1 compounds utilization and assimilation (reductive TCA cycle I) and fermentation (mixed acid fermentation) had also been exploited to be involved in resistance to SBPH.
Figure 2. The network of the metabolic pathways potentially related to SBPH resistance in the SBPH-resistant rice plant. The up-regulated metabolic pathway was marked with red box and the red arrows denote the steps in the up-regulated metabolic pathways; the down-regulated metabolic pathway was marked with green box and the green arrows denote the steps in the down-regulated metabolic pathways. Abreviations: IAA, indole-3-acetic acid; SAM, S-adenosyl-methionine; Dc-SAM, decarboxylated S-adenosylmethionine; VLCFA, the very long chain fatty acid; TCA cycle, tricarboxylicacidcycle acid cycle.
Figure 2. The network of the metabolic pathways potentially related to SBPH resistance in the SBPH-resistant rice plant. The up-regulated metabolic pathway was marked with red box and the red arrows denote the steps in the up-regulated metabolic pathways; the down-regulated metabolic pathway was marked with green box and the green arrows denote the steps in the down-regulated metabolic pathways. Abreviations: IAA, indole-3-acetic acid; SAM, S-adenosyl-methionine; Dc-SAM, decarboxylated S-adenosylmethionine; VLCFA, the very long chain fatty acid; TCA cycle, tricarboxylicacidcycle acid cycle.
Ijms 16 26128 g002

2.7. The Metabolic Pathways Potentially Related to SBPH Susceptibility in Rice

As in the identification of the metabolic pathways potentially related to SBPH resistance in rice, the metabolic pathways potentially related to SBPH susceptibility in rice were examined in light of the following classification:
Class I: in response to SBPH infestation, the metabolic pathway was differentially regulated in the susceptible rice plants (S6_S0) but not in the resistant rice plants (R6_R0); meanwhile, the basis of defense showed difference between resistant and susceptible rice plants after attack by SBPH (R6_S6) but not before attack by SBPH (R0_S0).
Class II: in response to SBPH infestation, the metabolic pathway was differentially regulated in the susceptible rice plants (S6_S0) but not in the resistant rice plants (R6_R0); meanwhile, the basis of defense showed a difference between resistant and susceptible rice plants before attack by SBPH (R0_S0) and after attack by SBPH (R6_S6).
Class III: in response to SBPH infestation, the metabolic pathway was differentially regulated in both susceptible rice plants (S6_S0) and resistant rice plants (R6_R0) and showed opposing differential changes in S6_S0 and R6_R0; meanwhile, the basis of defense showed a difference between susceptible and resistant rice plants after attack by SBPH (R6_S6) and or not before attack by SBPH (R0_S0). According to the classification of the metabolic pathways potentially related to SBPH susceptibility in rice, as shown in Table 3, a total of four and two metabolic pathways were identified in Classes II to III, respectively. The phenylalanine degradation III and ureide biosynthesis that were the metabolic pathways potentially related to SBPH susceptibility in rice were the same as that of the metabolic pathways potentially related to SBPH resistance in rice. As shown in Table 3 and Table S5, the primary metabolic pathways related to amino acid (lysine degradation I, phenylalanine degradation III) and fatty acids and lipids biosynthesis (phospholipid biosynthesis II) and the secondary metabolic pathways related to cofactors, prosthetic groups, electron carriers biosynthesis (methylerythritol phosphate pathway, thiamine biosynthesis ) as well as amine and polyamine (ureide biosynthesis) had been utilized potentially related to SBPH susceptibility in rice.
Table 3. The identified SBPH susceptibility-related metabolic pathways and its differential regulation in the four comparisons.
Table 3. The identified SBPH susceptibility-related metabolic pathways and its differential regulation in the four comparisons.
PathwayS6_S0R6_R0R0_S0R6_S6
Class II
methylerythritol phosphate pathwaydown lowerlower
thiamine biosynthesisdown lowerlower
lysine degradation Iup higherhigher
phospholipid biosynthesis IIup lowerlower
Class III
phenylalanine degradation IIIupdown lower
ureide biosynthesisupdownlowerlower

2.8. The Changes of Free Amino Acid Levels in Response to SBPH Infestation

The result of the identification of metabolic pathways potentially related to SBPH resistance and susceptibility in rice showed that several kinds of amino acids were found to be correlated with SBPH resistance or susceptibility in rice, therefore, the free amino acid levels were measured. As shown in the Table 4, in response to SBPH infestation, total free amino acid levels decreased 47% in the susceptible rice plants but there was no change in the resistant rice plants; glutamate was the most abundant amino acids and its level increased in the resistant rice plants but its level decreased in the susceptible rice plants; isoleucine level decreased greatly in the susceptible rice plants while cystathionine level decreased greatly in the resistant rice plants. The γ-amino-butyric acid showed the largest change in concentration in the resistant rice plants in response to SBPH infestation (an increase of about 900%) and the tyrosine was the second (an increase of about 300%). The tyrosine showed the largest change in concentration in the susceptible rice plants in response to SBPH infestation (an increase of about 140%). The level of arginine, phenylalanine, leucine, alanine, glycine, serine, alanine, threonine and leucine decreased in the susceptible rice plants but increased in the resistant rice plants; the level of a-amino-butyric acid decreased in the resistant plant but increased in the susceptible rice plants; the level of lysine, aspartic acid, valine, citrulline, methionine, cystine, ethanolamine, ß-aminoisobutyric acid, cystathionine and ornithine decreased in both resistant and susceptible rice plants; in addition to the level of tyrosine and γ-amino-butyric acid, the level of a-aminoadipic acid and proline also increased in both resistant and susceptible rice plants.
Table 4. The concentrations of free amino acids in the resistant (R) and susceptible (S) rice plants before (0 h) and after (6 h) attack by SBPH.
Table 4. The concentrations of free amino acids in the resistant (R) and susceptible (S) rice plants before (0 h) and after (6 h) attack by SBPH.
Free Amino AcidConcentration of Free Amino Acids (µg per g FW) xPercent Change of Concentration
S0R0S6R6S6_S0R6_R0
Isoleucine27.7 ± 2.743.5 ± 4.21.9 ± 0.412.5 ± 1.4−93.1−71.3
Citrulline12 ± 1.313 ± 1.42.4 ± 0.37 ± 0.8−80.0−46.2
Lysine27 ± 2.831.1 ± 3.46.3 ± 0.69.8 ± 0.9−76.7−68.5
Methionine11.80 ± 0.593.60 ± 0.28 11.03 ± 0.670.42 ± 0.04−69.5−96.2
Cystine11.44 ± 0.995.97 ± 0.4524.47 ± 1.7910.14 ± 0.91−47.8−58.6
Valine23.4 ± 1.923.1 ± 2.513.6 ± 1.516.8 ± 1.3−41.9−27.3
Aspartic acid26.95 ± 3.5915.68 ± 1.2323.70 ± 2.2521.65 ± 1.94−41.8−8.6
Cystathionine2.73 ± 0.251.74 ± 0.134.58 ± 0.680.39 ± 0.03−36.3−91.5
Ornithine2.1 ± 0.13.6 ± 0.11.4 ± 0.12.9 ± 0.2−33.3−19.4
β-Aminoisobutyric Acid3.6 ± 0.35 ± 0.42.7 ± 0.21.8 ± 0.2−25.0−64.0
Ethanolamine11 ± 1.39.1 ± 18.7 ± 0.97.5 ± 0.8−20.9−17.6
Phenylalanine18.5 ± 1.87.5 ± 0.33.2 ± 0.37.9 ± 0.5−82.75.3
Arginine90.4 ± 10.930 ± 3.418.6 ± 1.846.7 ± 4.8−79.455.7
β-Alanine6.4 ± 0.35.5 ± 0.43.9 ± 0.48.1 ± 0.6−39.147.3
Serine7.19 ± 0.714.43 ± 0.312.87 ± 0.265.46 ± 0.45−38.490.2
Glycine7.7 ± 19.3 ± 1.65.2 ± 0.411.1 ± 1.2−32.519.4
Leucine2.2 ± 0.16.3 ± 0.21.5 ± 0.111.3 ± 1.3−31.879.4
Glutamic Acid114.3 ± 8.8143.8 ± 14.281.9 ± 8.9160.4 ± 15.4−28.311.5
Threonine2.9 ± 0.31.9 ± 0.22.1 ± 0.23.6 ± 0.3−27.689.5
Alanine10 ± 1.110.9 ± 1.28.9 ± 0.420 ± 0.8−11.083.5
α-Amino-n-butyric Acid2.5 ± 0.23.6 ± 0.24 ± 0.22.9 ± 0.260.0−19.4
Tyrosine4.7 ± 0.33.7 ± 0.311.3 ± 1.214.7 ± 1.2140.4297.3
α-Aminoadipic Acid7.9 ± 0.611.5 ± 1.414.3 ± 0.828.4 ± 1.281.0147.0
γ-Amino-n-butyric Acid4.4 ± 0.42.3 ± 0.26.2 ± 0.522.9 ± 240.9895.7
Proline5 ± 0.32.6 ± 0.25.7 ± 0.35.3 ± 0.414.0103.8
Total443.7 ± 40.7434.1 ± 40.9235.3 ± 21.1439.6 ± 37.9−47.01.3
x Results are mean values of three replicates ± standard error. FW means fresh weight.

3. Discussion

In this study, we visualized the rice transcriptomic data of microarray expression profiling by using Pathway Tools Omics Viewer. By using the capabilities of Pathway Tools Omics Viewer for overlaying transcriptomic data onto the overview as a whole and animation of the comparative transcriptome analysis, our data suggested that the difference of change pattern between these two contrasting rice genotypes mostly lies in biosynthetic pathways when exposed to SBPH; while the obvious difference of change pattern between these two contrasting rice genotypes lies in energy metabolism pathways (Table S2). By using the capabilities of Pathway Tools Omics Viewer for generation a table containing all individual pathways with omics data painted onto the pathway, 166 metabolic pathways that were differentially regulated when exposed to SBPH were identified; of these, 64 metabolic pathways displayed similar change pattern in both two contrasting rice genotypes when exposed to SBPH (Table S4); in response to SBPH infestation, biosynthesis was the largest class of the differentially regulated metabolic pathways in both resistant and susceptible rice plants (R6_R0 and S6_S0; Table 1). By combination Pathway Tools Omics Viewer and the web tool Venn, 21 and 6 metabolic pathways were considered to be potentially associated with SBPH resistance and susceptibility in rice; these 21 SBPH resistance-related metabolic pathways were related to cuticle, phytoalexin, amine and polyamine, alkaloid, hormone, amino acid, fermentation and C1 compounds utilization and assimilation; in response to SBPH infestation, the metabolic pathways derived from phenylalanine including betanidin degradation, flavonoid biosynthesis, pinobanksin biosynthesis, salicylate biosynthesis and phenylalanine degradation III were down-regulated in the resistant rice plants (Figure 2); these six SBPH susceptibility-related metabolic pathways were related to amino acid, lipid, cofactors, prosthetic groups, electron carriers biosynthesis as well as amine and polyamine. Of the metabolic pathways that were considered to be potentially associated with SBPH resistance and susceptibility in rice, the change pattern of the phenylalanine degradation III and ureide biosynthesis in the resistant rice plants displayed opposing differential changes in the susceptible rice plants (Table 2 and Table 3), indicating these two metabolic pathways are highly potentially responsive to SBPH infestation.
The plant epidermis is considered to have an important role in defense against insect attack [23]. Observations implicated that rice plant epidermis plays a role in selection of food plant by BPH and that a high ratio of long to short carbon-chain compounds in BPH resistant rice varieties largely determined planthopper feeding responses [24]. These observations provide support for the viewpoint that increased activity in response to chemical characteristics of epidermis wax may play an important role in field resistance [24]. In this study, the epidermis wax related metabolic pathway, the very long chain fatty acid biosynthesis, was found to be up-regulated in the SBPH-resistant rice plants. The very long chain fatty acids are fatty acids of 20 to 36 carbons and are required for plant cuticle biosynthesis [25,26]. Recently, the very long chain fatty acid biosynthesis pathway has been found to be associated with plant defense; well-organized cuticle layers constitute the outermost layer of epidermal cells and thereafter act as the first natural barrier when encountered by pathogens [27]. Similarly, the very long chain fatty acid biosynthesis in the SBPH resistant rice plants was up-regulated when exposed to SBPH infestation.
It is reported that some phloem feeding insects induce the SA pathway while suppressing the plant’s JA response [28]. In response to planthopper attack, these antagonistic phenomena of SA and JA have also been observed. However, these antagonistic phenomena associated with rice resistance to planthopper seem to be in a contrasting manner. In response to BPH infestation, JA activity was suppressed while levels of SA were much higher accumulated [28]. In contrast, in this study, in the SBPH resistant rice plants, the JA biosynthesis was up-regulated while the SA biosynthesis and its branch pinobanksin biosynthesis were down-regulated. Growing evidence shows that plant defense against their different kinds of attackers is also in a hormones-dependent manner [29]. Besides JA and SA, other plant hormones such as auxin and cytokinin (CK) also serve as modulators of the hormone signaling backbone [29]. Auxin signaling is critical in regulation plant growth and development [30]. SBPH infestation seems to significantly modify auxin metabolism to decrease plant growth capacity. In the SBPH resistant rice plants, the inactivation of IAA (IAA, indole-3-acetic acid; IAA conjugate biosynthesis I and II), showed the strongest response among the 21 metabolic pathways potentially related to SBPH resistance in rice and was significantly up-regulated after attack by SBPH (Table S5). In rice-WBPH interactions, genes involved in auxin signaling such as auxin-responsive protein, auxin efflux carrier component 4, and auxin-responsive protein IAA14 were also strongly suppressed [31]. It has been reported that CK mediates insect resistance and also has an important role in the activation of JA biosynthesis [32]. For example, after wounding, the level of CK increased in transgenic tobacco (Nicotiana tabacum cv. Xanthinc) plants over-expressing a small GTP-binding protein, which resulted in increased rates of JA production [33]. In this study, we also observed that CK degradation was down-regulated and JA biosynthesis was up-regulated in the SBPH-resistant plant.
Betanidin is synthesized from tyrosine and is water-soluble nitrogen-containing pigment [34]. Pigment accumulation has been reported to be induced in response to insect feeding and the corresponding biosynthetic genes maybe directed against herbivores [35]. In the SBPH-resistant rice plants, betanidin degradation was down-regulated. It is therefore significant that tyrosine, a precursor of betanidin, showed the second most significant change in the SBPH-resistant rice plants, increasing about 300%. Phytoalexins are antimicrobial specialised compounds which are produced by plants as a response to pathogen attack [36]. Momilactones from rice are recognized as the major diterpenoid phytoalexins and growth inhibitors [37,38]. It has been reported that WBPH attack resulted in increased levels of momilactone A within 24 h [39]. Of the 21 metabolic pathways potentially related to SBPH resistance in rice, momilactone biosynthesis was the second most significantly altered pathway after attack by SBPH and was up-regulated in the SBPH-resistant rice plants (Table S5). In plants, flavonoids are thought to have many functions including attracting pollinating insects [40] and have been reported to stimulate the probing of rice WBPH and SBPH [41,42]. In accordance with the resistance, in the SBPH-resistant rice plants, flavonoid biosynthesis was down-regulated.
Amino acids and amino acid derivatives have been considered to have an important role in defense response [43]. Glutamate was the most abundant amino acid in the resistant rice plants before and after attack by SBPH (Table 4). Arginine is synthesized from glutamate and is the precursor of spermine and γ-amino-butyric acid [44,45]. The level of glutamate and arginine increased in the resistant rice plants but decreased in the susceptible rice plants after attack by SBPH. The γ-amino-butyric acid was the amino acid showing the largest change in concentration in the resistant rice plants in response to SBPH infestation (an increase of about 900%). Several lines of evidence indicate that increased γ-amino-butyric acid level is related to plant defense [46]. Spermine has also been implicated in plant defenses and its intermediate putrescine has been reported to directly affect insect growth and development and therefore could also act as a defense mechanism against some herbivorous insects [47]. Thus, the coordinated and linked changes pattern of glutamate and arginine level maybe enhance spermine and γ-amino-butyric acid biosynthesis and act as part of the defense against SBPH infestation. Meanwhile, before attack by SBPH, the level of isoleucine was the second highest, whereas, after attack by SBPH, the second highest level of amino acid was replaced by arginine. The two most abundant amino acid accounts for nearly one half of the total amino acid content when exposed to SBPH. Furthermore, although the functional significance of the phosphorylated pathway of serine biosynthesis in plants is not yet known, it is clear that in the phosphorylated pathway of serine biosynthesisthe glutamate functions as the substrate for the transfer of an amino group from glutamate to 3-phosphoserine [48]; in the resistant rice plants, the phosphorylated pathway of serine biosynthesiswas down-regulated in the resistant rice plants, meaning less glutamate used in the phosphorylated pathway of serine biosynthesis. In addition, the glutamate degradation III was down-regulated in the resistant rice plants. In the last step of synthesis of isoleucine, an amino group from glutamate was transferred to produce the isoleucine (http://pathway.iplantcollaborative.org); in the resistant rice plants, isoleucine level which was the second most abundant before attack showed the third significant decreased after attack, meaning less glutamate used in synthesis of isoleucine. Collectively, the increased N-rich glutamate and arginine in the resistant rice plants means that the N-rich soluble molecules accumulated. The accumulated N-rich soluble molecules have been considered as a means of nitrogen (N) and the carbon (C) store after wounding, which helps to prevent N losses during the wounding response and could work as a source for new growth after the wound recovery phase [49,50]. Additionally, C store was also modulated in the resistant rice plants during the SBPH response by down-regulation of the mixed acid fermentation in which the products are mixed acids, particularly equal amounts of carbon dioxide (CO2) [51,52] and by up-regulation of the reductive TCA cycle I in which two molecules of CO2 can be fixed [53].
Threonine has been conferred as the amino acid that resistant to Hyaloperonospora. arabidopsidis, presumably by altering the pathogen’s ability to grow under that condition; mutations in threonine biosynthesis gene displayed increased resistance to H. arabidopsidis and increased accumulation of threonine [54]. Similarly, the threonine level increased and its expression (aminopropanol biosynthesis, threonine degradation II and threonine degradation III (to methylglyoxal)) was down-regulated in the resistant rice plants when exposed to SBPH. Hopperburn partially results from the amides accumulation in infested rice leaf blade; the amides accumulate in the leaf-blade tissues because translocation systems of amides are not functioning [55]. There are two major forms of organic nitrogen compounds that can be transported amides (such as asparagine and glutamine) and ureides (such as allantoin and allantoate) [56]. In response to SBPH infestation, the ureide biosynthesis was down-regulated in the resistant rice plants but up-regulated in the susceptible rice plants, presumably indicating the amide transportation systems are not functioning in the susceptible rice plants.
Overall, more than one third of the 166 differentially regulated metabolic pathways displayed similar change pattern in both two contrasting rice genotypes when exposed to SBPH. Twenty-one differentially regulated metabolic pathways were associated with SBPH resistance in rice. Briefly, SBPH infestation modulated phytohormones such as SA, JA, auxin and cytokinin; particularly, SA and JA showed antagonistic phenomena; inactivation of auxin enhanced significantly and showed the strongest response, meaning a shift from growth and development to defense. Following SBPH infestation, the major diterpenoid phytoalexin momilactone A metabolism enhanced significantly and showed the second strongest response; the N-rich glutamate and arginine levels increased and were the most abundant in resistant rice plants, together with the down-regulation of mixed acid fermentation and up-regulation of the C1 compounds utilization and assimilation, suggesting a N and C store to prevent N and C losses and work as a source for new growth after the phase of wound recovery. In summary, SBPH infestation caused changes in the primary and secondary metabolism; rice plants protect themselves against SBPH at several levels. Of these, one is the well-organized cuticle layers as first barrier encountered by SBPH. The second is the biosynthesis of inducible chemical defences such as the major diterpenoid phytoalexin momilactone A. The third is a shift from growth and development to defense, enhancing the capacity of rice plants to main high fitness in the face of enemies. The fourth is the modulation of N and C store and work as a source for new growth after the phase of wound recovery. In infested resistant rice plants, metabolic processes presumably not only protect them against SBPH but also enhance the capacity of rice plants to repair the tissue that injured by SBPH attack.

4. Materials and Methods

4.1. Plant and Insect Materials

Two rice genotypes, Pf9279-4 and O. sativa (cv. 02428), the SBPH-resistant and SBPH-susceptible rice line, respectively, were selected for this study. Pf9279-4 is an introgression line derived from the asymmetric somatic hybridization between Oryza officinalis and O. sativa (cv. 02428). Pf9279-4 bears a major QTL Qsbph3d and two minor QTLs Qsbph7a and Qsbph12b [57]. The SBPH insects used for infestation were collected from rice fields in Jiangsu Province, China, that was confirmed to be non-viruliferous or viruliferous by dotimmunobinding assay and PCR detection [58]. The non-viruliferous population and viruliferous SBPH population were independently maintained on 02428 plants in the insect-proof in greenhouse in Zhejiang Academy of Agricultural Sciences. For simulation the natural reaction of rice attacked by SBPH in open field in this study, second to fourth instar nymphae of the mixture of the non-viruliferous and viruliferous SBPH population were transferred from 02428 and employed for infestation experiments.

4.2. SBPH Infestation Experiments

The resistant and susceptible rice plants were planted in pots (50 cm in length and 30 cm in width). When the seedlings were at the fourth leaf stage, weak plants were removed from the pots. In the insect-proof in greenhouse, second to fourth instar nymphs of the mixture of non-viruliferous and viruliferous SBPH population were introduced at a density of ten insects per seedling [57]. At each selected time-point, the two outmost layers of leaf sheaths were separated from the rice plants which had been fed and not fed by SBPH infestation, frozen immediately in liquid nitrogen and stored at −80 °C. The rice plants that had not been fed upon by SBPH were performed as a control experiment and its two outmost layers of leaf sheaths were separated and collected as controls.

4.3. Affymetrix GeneChip Analysis

According to the results of our previously comparative two-dimensional fluorescence difference gel electrophoresis analysis, significant differences in protein expression level between these two contrasting rice genotypes that sampled at different times after SBPH infestation were obvious within 6 h of infestation. Therefore, the two outmost layers of leaf sheaths of these two contrasting rice genotypes that had been fed upon by SBPH for 6 h and controls were used for microarray experiments. Beijing Biochip Co., Ltd. carried out the microarray experiments and data analysis. RNA purification, cDNA probe preparation, target hybridization, washing, staining, scanning, image analysis and the initial data processing were performed according to Affymetrix GeneChip® Expression Analysis Technical Manual. Three biological replications were performed. Expression levels were compared with those in Affymetrix® Microarray Suite version 5.0 (MAS 5.0, Sacramento, CA, USA). Significantly differentially expressed genes were screened in the Significant Analysis of Microarray software (SAM, version 3.02, Palo Alto, CA, USA) [59]. Significantly differentially expressed genes were determined with a fold change of  ≥2.0 and a q-value of <0.05 in the SAM output result. Checks of annotations from Affymetrix for significant differentially expressed genes were performed by using Blastx program in NCBI public database [60] and performed in TIGR Rice Genome Annotation version 4.0 (http://rice.tigr.org) [61]; corrected if necessary.

4.4. Validating of Rice Transcriptome Data for Pathway Tools Omics Viewer

The Pathway Tools Omics Viewer maps the expression values is in Gene ID-dependent manner [17]. Therefore, batch translation all probe set ids of the all differentially expressed genes to corresponding TIGR locus ids by using the Gene List Suite tool at the Plant Expression Database (http://www.plexdb.org) or by using the tool Omics Validator at the Gramene “Cyc” Pathways databases (http://pathway.iplantcollaborative.org/). For single experiment time step, listing all differentially expressed gene in a microsoft excel worksheet. The table in the microsoft excel worksheet contains 2 columns, the first row contains column header following TIGR locus id and the relative or absolute expression level of each differentially expressed gene, one gene per row. In the Pathway Tools Omics Viewer, for column numbering purposes, the first column, which contains the TIGR locus id, is column number 0 and the first potential data column is column number 1 [17]. Therefore, for time series, listing all differentially expressed gene in a microsoft excel worksheet, the first row contains column header following TIGR locus id and column numbers (with each column number corresponding to the relative or absolute expression level of each differentially expressed gene for a single timepoint), one gene per row. The microarray expression data is then saved as a tab-delimited data file and stored on the local computer [17].

4.5. Visualization of Rice Transcriptome Data Using Pathway Tools Omics Viewer

This section was performed as follows according to the guidelines at Pathway Tools Omics Viewer webpage (http://pathway.iplantcollaborative.org). It is worth noting that the rice pathway database should be first selected when the rice transcriptomic data was visualized [17]. Open the Pathway Tools Omics Viewer webpage (http://pathway.iplantcollaborative.org), click on the “change” link. In the dialog that pops up, select the Oryza sativa Japonica Group, then click OK to exit the dialog. The Pathway Tools Omics Viewer is next selected. Upload the transcriptomic data that validated as noted above. In the column number fields, enter 1 for single experiment time step and enter a list of column numbers (with each column number corresponding to a single timepoint) for animated time series. In the Color Scheme, “three color display with specified threshold” was selected and the specified threshold was entered in the specified threshold field. The other fields should be default selection. Then click the “Submit ”at the bottom of the Pathway Tools Omics Viewer webpage. When the Oryza sativa Japonica group metabolic map was generated, a metabolite can be identified by moving the mouse over the metabolite icon; the metabolite page or a related pathway page can be navigated by clicking on the metabolite icon [17].

4.6. Animation the Comparative Transcriptome Analysis Using Pathway Tools Omics Viewer

For single experiment time step, the metabolic map generated by Pathway Tools Omics Viewer is static but not animated. It is more ready and direct to visualize the expression patterns for single experiment time step if the metabolic map generated by Pathway Tools Omics Viewer is animated. However, in the Pathway Tools Omics Viewer, one comparative transcriptome analysis of interest, e.g., the differentially expressed genes between before and after infestation in the resistant rice plants, was considered as single experiment time step. Meanwhile, the single experiment time step can’t mimic animated time series by enter two 1 in the column number fields or copy the relative or absolute expression level of each differentially expressed gene in the second column for a third column in the validating of rice transcriptome data for Pathway Tools Omics Viewer process and enter 1 and 2 in the column number fields in the visualization of rice transcriptome data using Pathway Tools Omics Viewer process. To get animated pictures for comparative transcriptome analysis and more readily visualize the expression patterns, in this study, the individual metabolic pathways that were differentially regulated between before and after infestation in the resistant (R6_R0) and susceptible rice plants (S6_S0) respectively was taken as one combination for animated time series; while the individual metabolic pathways that were differentially regulated between resistant and susceptible rice plants before (R0_S0) and after infestation (R6_S6) respectively was taken as another combination for animated time series. Then mark the differentially regulated individual metabolic pathways between R6_S6 and R0_S0 or between R6_R0 and S6_S0 with green box.

4.7. Identification of the Differentially Regulated Individual Metabolic Pathways by the Pathway Tools Omics Viewer

Although a pathway in the metabolic map generated by pathway tools can be navigated, in most cases, e.g., comparative transcriptome analysis of interest, it is more ideal and direct to get clear understanding of the differentially regulated metabolic pathways in more details if the differentially expressed genes and its differential expression values can be mapped to the individual metabolic pathway. Notably, the Pathway Tools Omics Viewer provides an alternative display to generate a table containing all individual metabolic pathways [17]. Therefore, to get the differentially regulated individual metabolic pathways with omics data painted onto the metabolic pathway, in the “Display Type”, the third option, “generate a table of individual pathways exceeding threshold”, was selected. Then do the others as do for visualization of rice transcriptomic data using Pathway Tools Omics Viewer as described above. According to the fold changes of the differentially expressed genes that belong to the corresponding individual metabolic pathways derived from the Pathway Tools Omics Viewer, the individual metabolic pathway was then classified into three groups, if all the differentially expressed genes were up-regulated or higher, the individual metabolic pathway was up-regulated or higher; if all the differentially expressed genes were down-regulated or lower, the individual metabolic pathway was down-regulated or lower; if the differentially expressed genes were both up-regulated and down-regulated or higher and lower, the individual metabolic pathway was both up-regulated and down-regulated or both higher and lower and marked up/down or higher/lower in this report.

4.8. Identification of Metabolic Pathways Potentially Related to SBPH Resistance by Pathway Tools Omics Viewer

To identify the most relevant metabolic pathways associated with SBPH resistance and susceptibility in rice, the web tool Venn (http://bioinfogp.cnb.csic.es/) was next employed. The differentially regulated individual metabolic pathways in the four comparisons that identified by the Pathway Tools Omics Viewer were represented in a Venn diagram. In the Venn diagram webpage, click the numbers in the Venn diagram to see the results and copy each result of 15 categories to a microsoft excel worksheet. The table in the microsoft excel worksheet contains five columns, the first row contains column header following metabolic pathway and four comparisons, one metabolic pathway per row. According to the result of the individual metabolic pathway that was classified into three groups mentioned as above, mark the individual metabolic pathway with up or down or up/down as well as higher or lower or higher/lower. By searching the metabolic pathway in the Gramene database or click the link of the metabolic pathway in the result that displayed as a table containing individual metabolic pathway, download the classification hierarchy for individual metabolic pathways from the rice pathway database (http://pathway.iplantcollaborative.org/).

4.9. Measurement of Free Amino Acids

The result of the identification of the metabolic pathway associated with resistance and susceptibility to SBPH showed that several kinds of the amino acids were correlated with resistance and susceptibility to SBPH, therefore, the concentrations of the free amino acids were then determined as described by Lei [62]. Briefly, leaf sheath power (0.2 g) was sampled, sulfosalicylic acid (3%) was then added to a volume of 4 mL for protein precipitation, and extracted in 80 °C water for 30 min. After cooling to room temperature, the extract was centrifuged at 13,500× g at 4 °C for 10 min. The extracted supernatant was then filtered through a 0.45-um membrane filter. Finally, the supernatant (100 µL) was determined by using the automatic amino acid analyzer L-8900 (Hitachi Construction Machinery Co., Ltd., Tokyo, Japan).

Supplementary Materials

Acknowledgments

We sincerely thank Professor Qiusheng Yan and Xueqin Zhang (China National Rice Research Institute) for kindly providing the SBPH-resistant and susceptible rice lines. This research was funded by Research on Public Welfare Technology Application Projects of Zhejiang Province (Grant No. 2014C32015), Chinese National Basic Research Program (973 Program, Grant No. 2014CB138400 and 2014CB160309), the National Natural Science Foundation of China (Grant No. 31301289), the Zhejiang Provincial Natural Science Foundation of China under Grant No.LY12C13001 and the Ningbo Natural Science Foundation of China under Grant No. 2014A610200.

Author Contributions

Weilin Zhang, Chengqi Yan and Jianping Chen conducted the initial literature review and contributed to the study design and drafted the article; Li Mei measured the free amino acid levels; Ling Yang, Mei Li and Bojun Ma analyzed the experiment results, prepared figures and tables.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Choi, M.S.; Kim, Y.H.; Park, H.M.; Seo, B.Y.; Jung, J.K.; Kim, S.T.; Kim, M.C.; Shin, D.B.; Yun, H.T.; Choi, I.S.; et al. Expression of BrD1, a plant defensin from Brassica rapa, confers resistance against brown planthopper (Nilaparvata lugens) in transgenic rices. Mol. Cells 2009, 28, 131–137. [Google Scholar] [CrossRef] [PubMed]
  2. Brar, D.S.; Virk, P.S.; Jena, K.K.; Khush, G.S. Breeding for resistance to planthoppers in rice. In Planthoppers: New Threats to the Sustainability of Intensive Rice Production Systems in Asia; Heong, K., Hardy, B., Eds.; International Rice Research Institute: Los Baños, Philippines, 2009; pp. 401–428. [Google Scholar]
  3. Yuan, H.; Chen, X.; Zhu, L.; He, G. Identification of genes responsive to brown planthopper Nilaparvata lugens Stål (Homoptera: Delphacidae) feeding in rice. Planta 2005, 221, 105–112. [Google Scholar] [CrossRef] [PubMed]
  4. Du, B.; Zhang, W.; Liu, B.; Hu, J.; Wei, Z.; Shi, Z.; He, R.; Zhu, L.; Chen, R.; Han, B.; et al. Identification and characterization of Bph14, a gene conferring resistance to brown planthopper in rice. Proc. Natl. Acad. Sci. USA 2009, 106, 22163–22168. [Google Scholar] [CrossRef] [PubMed]
  5. Hu, J.; Zhou, J.; Peng, X.; Xu, H.; Liu, C.; Du, B.; Yuan, H.; Zhu, L.; He, G. The Bphi008a gene interacts with the ethylene pathway and transcriptionally regulates MAPK genes in the response of rice to brown planthopper feeding. Plant Physiol. 2011, 156, 856–872. [Google Scholar] [CrossRef] [PubMed]
  6. Lu, J.; Ju, H.; Zhou, G.; Zhu, C.; Erb, M.; Wang, X.; Wang, P.; Lou, Y. An EAR-motif-containing ERF transcription factor affects herbivore-induced signaling, defense and resistance in rice. Plant J. 2011, 68, 583–596. [Google Scholar] [CrossRef] [PubMed]
  7. Zheng, W.; Ma, L.; Zhao, J.; Li, Z.; Sun, F.; Lu, X. Comparative transcriptome analysis of two rice varieties in response to rice stripe virus and small brown planthoppers during early interaction. PLoS ONE 2013, 8, e82126. [Google Scholar] [CrossRef] [PubMed]
  8. Lv, W.; Du, B.; Shangguan, X.; Zhao, Y.; Pan, Y.; Zhu, L.; He, Y.; He, G. BAC and RNA sequencing reveal the brown planthopper resistance gene BPH15 in a recombination cold spot that mediates a unique defense mechanism. BMC Genom. 2014, 15, 674. [Google Scholar] [CrossRef] [PubMed]
  9. Zhang, F.; Zhu, L.; He, G. Differential gene expression in response to brown planthopper feeding in rice. J. Plant Physiol. 2004, 161, 53–62. [Google Scholar] [CrossRef] [PubMed]
  10. Hao, P.; Liu, C.; Wang, Y.; Chen, R.; Tang, M.; Du, B.; Zhu, L.; He, G. Herbivore-induced callose deposition on the sieve plates of rice: An important mechanism for host resistance. Plant Physiol. 2008, 146, 1810–1820. [Google Scholar] [CrossRef] [PubMed]
  11. Wang, Y.; Wang, X.; Yuan, H.; Chen, R.; Zhu, L.; He, R.; He, G. Responses of two contrasting genotypes of rice to brown planthopper. Mol. Plant Microbe Interact. 2008, 21, 122–132. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, Y.; Li, H.; Si, Y.; Zhang, H.; Guo, H.; Miao, X. Microarray analysis of broad-spectrum resistance derived from an indica cultivar Rathu Heenati. Planta 2012, 235, 829–840. [Google Scholar] [CrossRef] [PubMed]
  13. Riechmann, J.L.; Ratcliffe, O.J. A genomic perspective on plant transcription factors. Curr. Opin. Plant Biol. 2000, 3, 423–434. [Google Scholar] [CrossRef]
  14. Wang, Y.; Guo, H.; Li, H.; Zhang, H.; Miao, X. Identification of transcription factors potential related to brown planthopper resistance in rice via microarray expression profiling. BMC Genom. 2012, 13, 687. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, H.; Hao, J.; Chen, X.; Hao, Z.; Wang, X.; Lou, Y.; Peng, Y.; Guo, Z. Overexpression of rice WRKY89 enhances ultraviolet B tolerance and disease resistance in rice plants. Plant Mol. Biol. 2007, 65, 799–815. [Google Scholar] [CrossRef] [PubMed]
  16. Dharmawardhana, P.; Ren, L.; Amarasinghe, V.; Monaco, M.; Thomason, J.; Ravenscroft, D.; McCouch, S.; Ware, D.; Jaiswal, P. A genome scale metabolic network for rice and accompanying analysis of tryptophan, auxin and serotonin biosynthesis regulation under biotic stress. Rice 2013, 6, 15. [Google Scholar] [CrossRef] [PubMed]
  17. Paley, S.M.; Karp, P.D. The Pathway Tools cellular overview diagram and Omics Viewer. Nucleic Acids Res. 2006, 34, 3771–3778. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, P.; Dreher, K.; Karthikeyan, A.; Chi, A.; Pujar, A.; Caspi, R.; Karp, P.; Kirkup, V.; Latendresse, M.; Lee, C.; et al. Creation of a genome-wide metabolic pathway database for Populus trichocarpa using a new approach for reconstruction and curation of metabolic pathways for plants. Plant Physiol. 2010, 153, 1479–1491. [Google Scholar] [CrossRef] [PubMed]
  19. Oliveros, J. VENNY. An interactive tool for comparing lists with Venn Diagrams. Available online: http://bioinfogp.cnb.csic.es/tools/venny_old/index.html (accessed on 13 November 2014).
  20. Arimura, G.; Tashiro, K.; Kuhara, S.; Nishioka, T.; Ozawa, R.; Takabayashi, J. Gene responses in bean leaves induced by herbivory and by herbivore-induced volatiles. Biochem. Biophys. Res. Commun. 2000, 277, 305–310. [Google Scholar] [CrossRef] [PubMed]
  21. Agrawal, A.A. Induced responses to herbivory and increased plant performance. Science 1998, 279, 1201–1202. [Google Scholar] [CrossRef] [PubMed]
  22. Dicke, M. Induced responses to herbivory by R. Karban and I.T. Baldwin. Trends Ecol. Evol. 1998, 13, 83. [Google Scholar] [CrossRef]
  23. Chapman, R. The role of the leaf surface in food selection by acridids and other insects. Colloq. Int. Centre Natl. Rech. Sci. 1977, 265, 133–149. [Google Scholar]
  24. Woodhead, S.; Padgham, D.E. The effect of plant surface characteristics on resistance of rice to the brown planthopper, Nilaparvata lugens. Entomol. Exp. Appl. 1988, 47, 15–22. [Google Scholar] [CrossRef]
  25. Samuels, L.; Kunst, L.; Jetter, R. Sealing plant surfaces: Cuticular wax formation by epidermal cells. Annu. Rev. Plant Biol. 2008, 59, 683–707. [Google Scholar] [CrossRef] [PubMed]
  26. Shepherd, T.; Wynne, G.D. The effects of stress on plant cuticular waxes. New Phytol. 2006, 171, 469–499. [Google Scholar] [CrossRef] [PubMed]
  27. Raffaele, S.; Leger, A.; Roby, D. Very long chain fatty acid and lipid signaling in the response of plants to pathogens. Plant Signal. Behav. 2009, 4, 94–99. [Google Scholar] [CrossRef] [PubMed]
  28. Zhou, G.; Qi, J.; Ren, N.; Cheng, J.; Erb, M.; Mao, B.; Lou, Y. Silencing OsHI-LOX makes rice more susceptible to chewing herbivores, but enhances resistance to a phloem feeder. Plant J. 2009, 60, 638–648. [Google Scholar] [CrossRef] [PubMed]
  29. Robert-Seilaniantz, A.; Grant, M.; Jones, J.D. Hormone crosstalk in plant disease and defense: More than just jasmonate-salicylate antagonism. Annu. Rev. Phytopathol. 2011, 49, 317–343. [Google Scholar] [CrossRef] [PubMed]
  30. Pangesti, N.; Pineda, A.; Pieterse, C.M.; Dicke, M.; van Loon, J.J. Two-way plant mediated interactions between root-associated microbes and insects: From ecology to mechanisms. Front. Plant Sci. 2013, 4, 414. [Google Scholar] [CrossRef] [PubMed]
  31. Yang, Y.; Xu, J.; Leng, Y.; Xiong, G.; Hu, J.; Zhang, G.; Huang, L.; Wang, L.; Guo, L.; Li, J.; et al. Quantitative trait loci identification, fine mapping and gene expression profiling for ovicidal response to whitebacked planthopper (Sogatella furcifera Horváth) in rice (Oryza sativa L.). BMC Plant Biol. 2014, 14, 145. [Google Scholar] [CrossRef] [PubMed]
  32. Erb, M.; Meldau, S.; Howe, G.A. Role of phytohormones in insect-specific plant reactions. Trends Plant Sci. 2012, 17, 250–259. [Google Scholar] [CrossRef] [PubMed]
  33. Sano, H.; Seo, S.; Koizumi, N.; Niki, T.; Iwamura, H.; Ohashi, Y. Regulation by cytokinins of endogenous levels of jasmonic and salicylic acids in mechanically wounded tobacco plants. Plant Cell Physiol. 1996, 37, 762–769. [Google Scholar] [CrossRef]
  34. Strack, D.; Vogt, T.; Schliemann, W. Recent advances in betalain research. Phytochemistry 2003, 62, 247–269. [Google Scholar] [CrossRef]
  35. Casique-Arroyo, G.; Martinez-Gallardo, N.; de La Gonzalez, V.L.; Delano-Frier, J.P. Betacyanin biosynthetic genes and enzymes are differentially induced by (a)biotic stress in Amaranthus hypochondriacus. PLoS ONE 2014, 9, e99012. [Google Scholar] [CrossRef] [PubMed]
  36. Ahuja, I.; Kissen, R.; Bones, A.M. Phytoalexins in defense against pathogens. Trends Plant Sci. 2012, 17, 73–90. [Google Scholar] [CrossRef] [PubMed]
  37. Kato, T.; Kabuto, C.; Sasaki, N.; Tsunagawa, M.; Aizawa, H.; Fujita, K.; Kato, Y.; Kitahara, Y. Momilactones, growth inhibitors from rice, oryza sativa L. Tetrahedron Lett. 1973, 14, 3861–3864. [Google Scholar] [CrossRef]
  38. Yamane, H. Biosynthesis of phytoalexins and regulatory mechanisms of it in rice. Biosci. Biotechnol. Biochem. 2013, 77, 1141–1148. [Google Scholar] [CrossRef] [PubMed]
  39. Kanno, H.; Hasegawa, M.; Kodama, O. Accumulation of salicylic acid, jasmonic acid and phytoalexins in rice, Oryza sativa, infested by the white-backed planthopper, Sogatella furcifera (Hemiptera: Delphacidae). Appl. Entomol. Zool. 2012, 47, 27–34. [Google Scholar] [CrossRef]
  40. Verhoeyen, M.E.; Bovy, A.; Collins, G.; Muir, S.; Robinson, S.; de Vos, C.H.; Colliver, S. Increasing antioxidant levels in tomatoes through modification of the flavonoid biosynthetic pathway. J. Exp. Bot. 2002, 53, 2099–2106. [Google Scholar] [CrossRef] [PubMed]
  41. Adjei-Afriyie, F.; Kim, C.S.; Takemura, M.; Ishikawa, M.; Horiike, M. Isolation and identification of the probing stimulants in the rice plant for the white-back planthopper, Sogatella furcifera (Homoptera: Delphacidae). Biosci. Biotechnol. Biochem. 2000, 64, 443–446. [Google Scholar] [CrossRef] [PubMed]
  42. Adjei-Afriyie, F.; Kim, C.S.; Takemura, M.; Ishikawa, M.; Tebayashi, S.; Horiike, M. Probing stimulants from the rice plant towards the smaller brown planthopper, Laodelphax striatellus (Fallén) (Homoptera: Delphacidae). Z. Naturforschung C 2000, 55, 1038–1043. [Google Scholar] [CrossRef]
  43. Steinbrenner, A.D.; Gomez, S.; Osorio, S.; Fernie, A.R.; Orians, C.M. Herbivore-induced changes in tomato (Solanum lycopersicum) primary metabolism: A whole plant perspective. J. Chem. Ecol. 2011, 37, 1294–1303. [Google Scholar] [CrossRef] [PubMed]
  44. Gill, S.S.; Tuteja, N. Polyamines and abiotic stress tolerance in plants. Plant Signal. Behav. 2010, 5, 26–33. [Google Scholar] [CrossRef] [PubMed]
  45. Kusano, T.; Yamaguchi, K.; Berberich, T.; Takahashi, Y. Advances in polyamine research in 2007. J. Plant Res. 2007, 120, 345–350. [Google Scholar] [CrossRef] [PubMed]
  46. Huang, T.; Jander, G.; de Vos, M. Non-protein amino acids in plant defense against insect herbivores: Representative cases and opportunities for further functional analysis. Phytochemistry 2011, 72, 1531–1537. [Google Scholar] [CrossRef] [PubMed]
  47. Malik, N.S.; Perez, J.L.; Kunta, M.; Patt, J.M.; Mangan, R.L. Changes in free amino acids and polyamine levels in Satsuma leaves in response to Asian citrus psyllid infestation and water stress. Insect Sci. 2014, 21, 707–716. [Google Scholar] [CrossRef] [PubMed]
  48. Cascales-Minana, B.; Munoz-Bertomeu, J.; Flores-Tornero, M.; Anoman, A.D.; Pertusa, J.; Alaiz, M.; Osorio, S.; Fernie, A.R.; Segura, J.; Ros, R. The phosphorylated pathway of serine biosynthesis is essential both for male gametophyte and embryo development and for root growth in Arabidopsis. Plant Cell 2013, 25, 2084–2101. [Google Scholar] [CrossRef] [PubMed]
  49. Zhou, L.; Thornburg, R. Wound-inducible genes in plants. In Inducible Gene Expression; Reynolds, P.H.S., Ed.; Book News: Portland, OR, USA, 1999; pp. 127–158. [Google Scholar]
  50. Sardans, J.; Gargallo-Garriga, A.; Perez-Trujillo, M.; Parella, T.J.; Seco, R.; Filella, I.; Penuelas, J. Metabolic responses of Quercus ilex seedlings to wounding analysed with nuclear magnetic resonance profiling. Plant Biol. 2014, 16, 395–403. [Google Scholar] [CrossRef] [PubMed]
  51. Huang, J.J.; Chen, K.W.; Ruaan, M.K. Mixed acid fermentation of glucose as a mechanism of emphysematous urinary tract infection. J. Urol. 1991, 146, 148–151. [Google Scholar] [PubMed]
  52. Trchounian, A.; Gary, S.R. Novel insights into the bioenergetics of mixed-acid fermentation: Can hydrogen and proton cycles combine to help maintain a proton motive force? IUBMB Life 2014, 66, 1–7. [Google Scholar] [CrossRef] [PubMed]
  53. Hugler, M.; Wirsen, C.O.; Fuchs, G.; Taylor, C.D.; Sievert, S.M. Evidence for autotrophic CO2 fixation via the reductive tricarboxylic acid cycle by members of the epsilon subdivision of proteobacteria. J. Bacteriol. 2005, 187, 3020–3027. [Google Scholar] [CrossRef] [PubMed]
  54. Stuttmann, J.; Hubberten, H.M.; Rietz, S.; Kaur, J.; Muskett, P.; Guerois, R.; Bednarek, P.; Hoefgen, R.; Parker, J.E. Perturbation of Arabidopsis amino acid metabolism causes incompatibility with the adapted biotrophic pathogen Hyaloperonospora arabidopsidis. Plant Cell 2011, 23, 2788–2803. [Google Scholar] [CrossRef] [PubMed]
  55. Yoshida, R.; Oritani, T.; Nishi, A. Studies on nitrogen metabolism in crop plants. VIII. Occurrence of kinetic-like factor in root exudates of rice plant. Proc. Crop Sci. Soc. Jpn. 1970, 39, 363–369. [Google Scholar] [CrossRef]
  56. Schubert, K.R. Enzymes of purine biosynthesis and catabolism in glycine max: I. Comparison of activities with N2 fixation and composition of xylem exudates during nodule development. Plant Physiol. 1981, 68, 1115–1122. [Google Scholar] [CrossRef] [PubMed]
  57. Zhang, W.; Dong, Y.; Yang, L.; Ma, B.; Ma, R.; Huang, F.; Wang, C.; Hu, H.; Li, C.; Yan, C.; et al. Small brown planthopper resistance loci in wild rice (Oryza officinalis). Mol. Genet. Genom. 2014, 289, 373–382. [Google Scholar] [CrossRef] [PubMed]
  58. Qin, S.W.; Gao, D.M.; Chen, S.X. Studies on techniques of rapid detecting rice stripe virus in Laodelphax striatellus. Zhejiang J. Agric. Sci. 1994, 6, 226–229. [Google Scholar]
  59. Tusher, V.G.; Tibshirani, R.; Chu, G. Significance analysis of microarrays applied to the ionizing radiation response. Proc. Natl. Acad. Sci. USA 2001, 98, 5116–5121. [Google Scholar] [CrossRef] [PubMed]
  60. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [PubMed]
  61. Yuan, Q.; Ouyang, S.; Wang, A.; Zhu, W.; Maiti, R.; Lin, H.; Hamilton, J.; Haas, B.; Sultana, R.; Cheung, F.; et al. The institute for genomic research Osa1 rice genome annotation database. Plant Physiol. 2005, 138, 18–26. [Google Scholar] [CrossRef] [PubMed]
  62. Lei, Y.; Korpelainen, H.; Li, C. Physiological and biochemical responses to high Mn concentrations in two contrasting Populus cathayana populations. Chemosphere 2007, 68, 686–694. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Zhang, W.; Yang, L.; Li, M.; Ma, B.; Yan, C.; Chen, J. Omics-Based Comparative Transcriptional Profiling of Two Contrasting Rice Genotypes during Early Infestation by Small Brown Planthopper. Int. J. Mol. Sci. 2015, 16, 28746-28764. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms161226128

AMA Style

Zhang W, Yang L, Li M, Ma B, Yan C, Chen J. Omics-Based Comparative Transcriptional Profiling of Two Contrasting Rice Genotypes during Early Infestation by Small Brown Planthopper. International Journal of Molecular Sciences. 2015; 16(12):28746-28764. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms161226128

Chicago/Turabian Style

Zhang, Weilin, Ling Yang, Mei Li, Bojun Ma, Chengqi Yan, and Jianping Chen. 2015. "Omics-Based Comparative Transcriptional Profiling of Two Contrasting Rice Genotypes during Early Infestation by Small Brown Planthopper" International Journal of Molecular Sciences 16, no. 12: 28746-28764. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms161226128

Article Metrics

Back to TopTop