RNA-Seq and Electrical Penetration Graph Revealed the Role of Grh1-Mediated Activation of Defense Mechanisms towards Green Rice Leafhopper (Nephotettix cincticeps Uhler) Resistance in Rice (Oryza sativa L.)

The green rice leafhopper (GRH, Nephotettix cincticeps Uhler) is one of the most important insect pests causing serious damage to rice production and yield loss in East Asia. Prior to performing RNA-Seq analysis, we conducted an electrical penetration graph (EPG) test to investigate the feeding behavior of GRH on Ilpum (recurrent parent, GRH-susceptible cultivar), a near-isogenic line (NIL carrying Grh1) compared to the Grh1 donor parent (Shingwang). Then, we conducted a transcriptome-wide analysis of GRH-responsive genes in Ilpum and NIL, which was followed by the validation of RNA-Seq data by qPCR. On the one hand, EPG results showed differential feeding behaviors of GRH between Ilpum and NIL. The phloem-like feeding pattern was detected in Ilpum, whereas the EPG test indicated a xylem-like feeding habit of GRH on NIL. In addition, we observed a high death rate of GRH on NIL (92%) compared to Ilpum (28%) 72 h post infestation, attributed to GRH failure to suck the phloem sap of NIL. On the other hand, RNA-Seq data revealed that Ilpum and NIL GRH-treated plants generated 1,766,347 and 3,676,765 counts per million mapped (CPM) reads, respectively. The alignment of reads indicated that more than 75% of reads were mapped to the reference genome, and 8859 genes and 15,815,400 transcripts were obtained. Of this number, 3424 differentially expressed genes (DEGs, 1605 upregulated in Ilpum and downregulated in NIL; 1819 genes upregulated in NIL and downregulated in Ilpum) were identified. According to the quantile normalization of the fragments per kilobase of transcript per million mapped reads (FPKM) values, followed by the Student’s t-test (p < 0.05), we identified 3283 DEGs in Ilpum (1935 upregulated and 1348 downregulated) and 2599 DEGs in NIL (1621 upregulated and 978 downregulated) with at least a log(2) (logarithm base 2) twofold change (Log(2)FC ≥2) in the expression level upon GRH infestation. Upregulated genes in NIL exceeded by 13.3% those recorded in Ilpum. The majority of genes associated with the metabolism of carbohydrates, amino acids, lipids, nucleotides, the activity of coenzymes, the action of phytohormones, protein modification, homeostasis, the transport of solutes, and the uptake of nutrients, among others, were abundantly upregulated in NIL (carrying Grh1). However, a high number of upregulated genes involved in photosynthesis, cellular respiration, secondary metabolism, redox homeostasis, protein biosynthesis, protein translocation, and external stimuli response related genes were found in Ilpum. Therefore, all data suggest that Grh1-mediated resistance against GRH in rice would involve a transcriptome-wide reprogramming, resulting in the activation of bZIP, MYB, NAC, bHLH, WRKY, and GRAS transcription factors, coupled with the induction of the pathogen-pattern triggered immunity (PTI), systemic acquired resistance (SAR), symbiotic signaling pathway, and the activation of genes associated with the response mechanisms against viruses. This comprehensive transcriptome profile of GRH-responsive genes gives new insights into the molecular response mechanisms underlying GRH (insect pest)–rice (plant) interaction.