Abstract
The brown planthopper (BPH) is one of the most problematic pests affecting rice (Oryza sativa L.) yields in Asia. Breeding rice varieties containing resistance genes is the most economical and effective means of controlling BPH. In this study, the key factors in resistance to BPH were investigated between the high-resistance rice variety "R26" and the susceptible variety "TN1" using RNA-sequencing. We identified 9527 differentially expressed genes (DEGs) between the rice varieties under BPH-induced stress. Weighted time-course gene co-expression network analysis (WGCNA) indicated that the increased expression of genes is associated with plant hormones, MAPK signaling pathway and biosynthesis of other secondary metabolites, which were involved in disease resistance. A connection network identified a hub gene, OsREM4.1 (BGIOSGA024059), that may affect rice resistance to the BPH. Knocking out OsREM4.1 in rice can lead to a decrease in callose, making it less resistant to BPH. Overall, the expression of differentially expressed genes varies among rice varieties with different resistance in BPH invasion. Inaddition, R26 enhances resistance to BPH by upregulating genes and secondary metabolites related to stress resistance and plant immunity. In summary, our study provides valuable insights into the genome-wide expression profile of DEGs in rice under BPH invasion through high-throughput sequencing, and further suggests that R26 can be used to develop high resistance rice lines in BPH resistant breeding programs.