Bin Yu1
Mengjia Geng1
Yu Xue1Qingqing Yu1Bojie Lu2Miao Liu1Yuhan Shao1Chenxi Li1
Jingang Xu1Jintao Li1
Wei Hu3
Hengmin Tang1
Peng Li1
Qingsong Liu1,4*
Shengli Jing1*- 1College of Life Sciences, Xinyang Normal University, Xinyang, China
- 2Hubei Provincial Key Laboratory for Protection and Application of Special Plant Germplasm in Wuling Area of China, College of Life Sciences, South-Central Minzu University, Wuhan, China
- 3Guangdong Provincial Key Laboratory of New Technology in Rice Breeding, Rice Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou, China
- 4State Key Laboratory of Cotton Bio-breeding and Integrated Utilization, State Key Laboratory of Crop Stress Adaptation and Improvement, Key Laboratory of Plant Stress Biology, School of Life Sciences, Henan University, Kaifeng, China
Introduction: The brown planthopper (BPH) poses a significant threat to rice production in Asia. The use of resistant rice varieties has been effective in managing this pest. However, the adaptability of BPH to resistant rice varieties has led to the emergence of virulent populations, such as biotype Y BPH. YHY15 rice, which carries the BPH resistance gene Bph15, exhibits notable resistance to biotype 1 BPH but is susceptible to biotype Y BPH. Limited information exists regarding how resistant rice plants defend against BPH populations with varying levels of virulence.
Methods: In this study, we integrated miRNA and mRNA expression profiling analyses to study the differential responses of YHY15 rice to both avirulent (biotype 1) and virulent (biotype Y) BPH.
Results: YHY15 rice demonstrated a rapid response to biotype Y BPH infestation, with significant transcriptional changes occurring within 6 hours. The biotype Y-responsive genes were notably enriched in photosynthetic processes. Accordingly, biotype Y BPH infestation induced more intense transcriptional responses, affecting miRNA expression, defenserelated metabolic pathways, phytohormone signaling, and multiple transcription factors. Additionally, callose deposition was enhanced in biotype Y BPH-infested rice seedlings.
Discussion: These findings provide comprehensive insights into the defense mechanisms of resistant rice plants against virulent BPH, and may potentially guide the development of insect-resistant rice varieties.
1 Introduction
Rice (Oryza sativa L.) was domesticated approximately 10,000 years ago in the lower Yangtze Valley in China. From there, it spreads across Asia, Africa, Europe, and the Americas, and now serves as a staple crop for more than half of the world’s population (Cheng et al., 2013; Gutaker et al., 2020). However, rice production suffers from numerous pests and pathogens. Among them, the brown planthopper (BPH; Nilaparvata lugens Stål) is considered extremely destructive (Bottrell and Schoenly, 2012; Cheng et al., 2013; Jing et al., 2017). BPH feeds on phloem sap and causes dwarfing, wilting, browning, drying, and ultimately death in severe cases (Wang et al., 2008). In addition, BPH serves as a vector for viral diseases, resulting in significant yield shortfalls and economic losses (Cheng et al., 2013).
Over the course of an evolutionary arms race between these two species, rice has developed sophisticated defensive mechanisms against BPH (Cheng et al., 2013), including both basic defenses and resistance (R) gene-mediated defenses. To date, approximately 40 major BPH resistance genes have been identified in cultivated and wild rice species (Wang et al., 2023). Molecular cloning and functional characterization of BPH resistance genes have clarified the molecular mechanisms of rice resistance to BPH (Zheng et al., 2021). One mechanism involves the occlusion of sieve tubes with callose, which is a common plant defense against sap-sucking insects. This mechanism is the most effective in rice varieties carrying BPH resistance genes such as B5 (carrying Bph14 and Bph15), RI35 (carrying Bph14), and YHY15 (carrying Bph15) (Hao et al., 2008). In susceptible varieties such as Taichung Native1 (TN1), β-1,3-glucanases (which are only weakly induced in resistant plants) decompose the deposited callose and thereby facilitate continuous feeding by BPH (Hao et al., 2008). Bph14 encodes a typical CC-NB-LRR protein which interacts with the transcription factors (TFs) OsWRKY46 and OsWRKY72 to increase the expression of the receptor-like cytoplasmic kinase gene RLCK281 and the callose synthase gene by binding to their promoters (Du et al., 2009; Hu et al., 2017). Bph15 is located on the short arm of chromosome 4 and is composed of a gene cluster encoding three plant lectin receptor-like kinase proteins (LecRLKs) (Yang et al., 2004; Lv et al., 2014; Xiao et al., 2016). However, the precise mechanism by which Bph15 confers resistance to BPH is unknown.
Phytohormones are thought to play pivotal roles in the interaction of rice plants and BPH herbivory, including salicylic acid (SA), ethylene (ET), jasmonic acid (JA), cytokinin (CK), brassinosteroid (BR), and abscisic acid (ABA). The SA pathway contributes to the immune response against piercing–sucking insects and is involved in Bph6-, Bph9-, Bph14-, and Bph29-mediated resistance in rice (Du et al., 2009; Wang et al., 2015; Zhao et al., 2016; Hu et al., 2017; Guo et al., 2018). Antagonistically, ET negatively regulates BPH resistance in rice (Lu et al., 2011, 2014; Ma et al., 2020). However, the role of JA in rice resistance to BPH remains controversial as the silencing of different genes related to JA biosynthesis and signaling results in diverse impacts on BPH resistance. For example, silencing 9-lipoxygenase (OsLOX9/OsHI-LOX, a JA biosynthesis-related gene) enhances BPH resistance, suggesting that JA negatively regulates BPH resistance (Zhou et al., 2009). On the other hand, silencing coronatine insensitive 1 (OsCOI1, a JA receptor gene) has no effect on BPH resistance (Ye et al., 2012). According to studies in allene oxide cyclase (AOC, a JA biosynthesis-related enzyme)- and Myelocytomatosis protein 2 (MYC2, a bHLH TF in the JA pathway)-knockout mutant rice, JA appears to be a positive regulator of BPH resistance (Xu et al., 2021). The exogenous application of JA similarly suggests this (Guo et al., 2018). Recent research also suggests that CK may positively regulate BPH resistance in a JA-dependent manner (Zhang et al., 2022). Conversely, BR promotes BPH susceptibility by modulating SA and JA signaling (Pan et al., 2018). ABA enhances BPH resistance by promoting callose formation (Dinh et al., 2013; Liu et al., 2017) and synergizes with JA to stimulate the expression of TFs in BPH-infested rice (Li et al., 2022). Furthermore, a coordinated CK, SA, and JA signaling network has been found to be activated in Bph6-near isogenic lines (NILs) (Guo et al., 2018). Taken together, these findings suggest that the various phytohormones play diverse roles in the BPH defense response, and that there is complex crosstalk between them.
BPH biotypes with increased virulence have emerged in response to pressures imposed by these defense mechanisms, which are capable of overcoming resistance conferred by major resistance genes. Biotype 1 BPH, which exhibits low virulence on resistant rice varieties, is widely distributed across southeast Asia and primarily parasitizes susceptible varieties such as TN1 (Alam & Cohen, 1998). Rearing biotype 1 BPH on resistant rice variety YHY15 (carrying Bph15) for several years resulted in the development of highly-virulent biotype Y BPH, which are able to overcome resistance conferred by Bph15 (Jing et al., 2011). Similarly, after force-feeding 40 generations of local BPH using resistant IR56 rice (carrying Bph3), the resulting BPH population (IR56-BPH) was able to overcome Bph3-conferred resistance (Zheng et al., 2016). In addition, the Mudgo BPH population has been reported to cause substantial damage to Mudgo rice plants (carrying Bph1) (Ji et al., 2013; Wan et al., 2019). In response to selective pressures imposed by these resistance genes, BPH populations accumulate adaptations over generations which eventually allow them to overcome such resistance mechanisms. In effect, this process constitutes a loss of resistance in formerly-resistant rice varieties against them. However, the effectiveness of resistance strategies used by different rice varieties against BPH populations with varying levels of virulence requires further investigation. Such combinations of avirulent/virulent BPH and resistant rice provide ideal models for studying resistance adaptation mechanisms.
Just as plants have evolved intricate defense mechanisms to protect themselves against herbivorous insects, insects have in turn developed strategies to overcome plant defenses (Liu Q. et al., 2021). Interactions between plants and insects involve an array of molecular, biochemical, and physiological processes occurring at multiple levels. Multi-omics analyses integrate genomics, transcriptomics, proteomics, and metabolomics, as well as other “omics” approaches, in order to clarify the intricate signaling pathways, molecular responses, and biochemical processes involved in the dynamic interplay between plants and insects (Wang et al., 2020; Shi et al., 2023). Transcriptional profiling has aided our understanding of the defense mechanisms utilized by rice against BPH. For example, research suggests that BPH infestation results in the upregulation of genes involved in signaling, oxidative stress, pathogen-related response, and macromolecule degradation, as well as the downregulation of genes associated with flavonoid biosynthesis, photosynthesis, and cell growth (Zhang et al., 2004; Yuan et al., 2005; Wang et al., 2008). In addition, microarray analyses of BPH-infested Rathu Heenati and TN1 rice underscore the importance of TFs and phytohormones in the defense response (Wang et al., 2012; Li et al., 2017).
MicroRNAs (miRNAs) are small (approximately 21 nt in length) regulatory RNAs produced through endonucleolytic processing of hairpin precursors (Axtell and Meyers, 2018). miRNAs regulate gene expression by binding to complementary sequences in mRNA molecules, resulting in degradation and/or translational inhibition (Bartel, 2009). In plants, miRNAs are involved in various processes such as phytohormone signaling; abiotic and biotic stress response (Zhang et al., 2013, 2016; Li et al., 2017; Natarajan et al., 2018; Salvador-Guirao et al., 2018; Zhang et al., 2018; Liu X. et al., 2021); and leaf, flower, shoot, root, and vascular tissue development (Marin et al., 2010; Peng et al., 2014; Mangrauthia et al., 2017). However, only a few miRNAs have been found to play roles in the insect-plant interaction (Jing et al., 2023). For example, in an investigation of resistant and susceptible rice varieties, the BPH-responsive miRNAs miR156 and miR396 were found to negatively regulate BPH resistance by regulating JA and flavonoid biosynthesis, respectively (Wu et al., 2017; Ge et al., 2018; Dai et al., 2019). Studies of BPH-responsive mRNA and miRNA transcriptomes have uncovered certain universal responses of rice to BPH infestation (Li et al., 2017; Wu et al., 2017; Nanda et al., 2020). Additionally, integrated expression profiling of miRNAs and target genes associated with the BPH-rice interaction has been conducted (Tan et al., 2020). Such combined miRNA and mRNA analyses will be key to unraveling the transcriptional responses of rice to BPH infestation.
In this study, we employed high-throughput sequencing to analyze the mRNA and miRNA expression profiles of YHY15 rice seedlings infested with either biotype 1 BPH (avirulent population) or biotype Y BPH (virulent population). Furthermore, we combined sequence analysis and physiological assays to reveal the underlying resistance mechanisms of rice against BPH. The findings presented in this work will provide a valuable resource for further genome-wide investigations of BPH-responsive genes, as well as studies of Bph15-mediated resistance. Moreover, these findings improve our understanding of the intricate interactions between rice and BPH, and may be used in the development of effective BPH management strategies.
2 Materials and methods
2.1 Plant and insect materials
In this study, we used two rice varieties (Taichung Native1 [TN1] and YHY15) and two BPH populations (biotype 1 and biotype Y). TN1 is a susceptible rice cultivar, while YHY15 is a recombinant inbred line (RIL) derived from the RI93 × TN1 F2 population carrying resistance gene Bph15 (Yang et al., 2004). Biotype 1 BPH (avirulent) originated from Wuhan University (China) and are reared on TN1. Biotype Y BPH (virulent) were developed by rearing biotype 1 BPH on YHY15 plants beginning in January 2007 (Jing et al., 2011). All rice plants were grown from seeds sown in sponges (6 cm diameter, 2 cm height), with 8 rice plants per cup. Rice plants were reared in a controlled-environment incubator maintained at 30 ± 2 °C during daytime hours (16 h, 06:00–22:00) and 28 ± 2 °C during nighttime hours (8 h, 22:00–06:00). Rice plants were grown for approximately 2-5 weeks following sowing, depending on experimental needs. The BPH populations were reared at Xinyang Normal University (China) under the following conditions: 26 ± 1 °C, 16 h light/8 h dark cycle. Third instar BPH nymphs were used for the infestation experiments.
2.2 Evaluation of BPH resistance in rice
At the two-leaf (2-week-old) stage, YHY15 rice plants were infested with third instar biotype 1 or biotype Y BPH nymphs at a rate of 15 nymphs per seedling. The growth status of each plant was photographically recorded daily until all biotype Y-infested seedlings died. Each experiment consisted of three biological replicates.
2.3 Measurement of BPH weight gain and honeydew excretion
BPH weight gain and honeydew excretion were measured as described previously (Shi et al., 2021). Briefly, newly-emerged adult female BPH were weighed using an electronic balance (Mettler Toledo, MS105DU, Switzerland) and subsequently separated into pre-weighed parafilm sachets (2 × 2.5 cm) fixed to the leaf sheaths of 4-week-old rice plants. After 48 h, the insects were carefully removed from the sachets, and both the insect and the honeydew in each sachet were separately weighed. Weight gain was calculated by comparing each insect’s weight before and after feeding, and the weight gain ratio was calculated by dividing the weight gain by the initial weight. Both the weight gain and honeydew excretion assays were conducted using at least 37 replicates.
2.4 BPH infestation and sample collection
Three-week-old YHY15 rice seedlings were infested with third-instar biotype 1 (RT) and biotype Y (RY) BPH at a rate of 15 nymphs per seedling. Each treatment consisted of three biological replicates, with six seedlings per replicate. For RNA-seq and miRNA-seq analysis, leaf sheaths were collected from non-infested controls (0 h), and during early (6 h) and late (48 h) infestation. According to their time of collection, samples of non-infested rice plants were labeled ‘R0’, samples of biotype 1-infested rice plants were labeled as ‘RT6’ or ‘RT48’, and samples of biotype Y-infested rice plants were labeled as ‘RY6’ or ‘RY48’. Each sampled leaf sheath blade was excised, frozen in liquid nitrogen, and stored at -80 °C for further use.
2.5 mRNA transcriptome sequencing and analysis
2.5.1 RNA extraction, quantification, and qualification
Total RNA was isolated from rice samples using RNA Trizol reagent (Life Technologies, NY, USA), following the manufacturer’s instructions. RNA degradation and contamination were evaluated using 1% agarose gels. RNA purity was quantified using a NanoPhotometer spectrophotometer (Implen, CA, USA). RNA integrity was assessed using an RNA Nano 6000 Assay Kit for the Bioanalyzer 2100 system (Agilent Technologies, CA, USA).
2.5.2 Library construction, quality control, and sequencing
RNA libraries were generated using 1 µg of RNA per sample and the NEBNext Ultra RNA Library Prep Kit for Illumina. Index codes were incorporated to distinguish between samples. Briefly, mRNA was purified using poly-T oligo-attached magnetic beads, followed by fragmentation using divalent cations. First-strand cDNA synthesis was carried out using random hexamer primers and M-MuLV Reverse Transcriptase. The second strand was synthesized using DNA Polymerase I and RNase H. Overhangs were blunted using exonucleases or polymerases and NEBNext adaptors were ligated following adenylation. The library fragments (250~300 bp) were purified using an AMPure XP system. The USER Enzyme was applied to size-selected, adaptor-ligated cDNA prior to PCR. PCR was carried out using Phusion High-Fidelity DNA polymerase, Universal PCR primers, and Index (X) Primer. The PCR products were purified using an AMPure XP system, and library quality was evaluated using an Agilent Bioanalyzer 2100 system. Index-coded samples were clustered with a cBot Cluster Generation System using a TruSeq PE Cluster Kit v3-cBot-HS (Illumina), according to the manufacturer’s instructions. After clustering, the libraries were sequenced on an Illumina NovaSeq platform, generating 150 bp paired-end reads.
2.5.3 Data analysis
Raw data (fastq format) were first processed using in-house perl scripts. In this step, low-quality reads, adapter sequences, and poly-Ns were removed. Subsequently, the Q20, Q30, and GC content of the clean reads were calculated. All downstream analyses utilized only clean, high-quality data. Next, the reference genome and gene model annotation files were downloaded. An index of the reference genome was constructed, and clean paired-end reads were aligned to the reference genome, using Hisat2 (v2.0.5). The mapped reads from each sample were assembled with StringTie (v1.3.3b) using a reference-based approach. featureCounts (v1.5.0-p3) was used to count the number of reads mapped to each gene. The Transcripts Per Kilobase Million (TPM) of each gene was quantified based on the gene length and the number of reads mapped to the gene. Differentially expressed genes (DEGs) were identified using the R (v1.16.1) package DESeq2, with three biological replicates per treatment, according to the following criteria: P-value < 0.05, false discovery rate (FDR) < 5, and absolute value of log2 fold change (FC) ≥ 1. The DEGs underwent additional screening through soft clustering using the Mfuzz package, employing a fuzzy c-means algorithm, as previously reported (Kumar and Futschik, 2007). DEGs exhibiting similar expression patterns were categorized into 20 clusters, and the genes within these clusters were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. GO annotations were downloaded from NCBI (http://www.ncbi.nlm.nih.gov/) and GO (http://www.geneontology.org/). The KEGG database was used to identify BPH-responsive pathways. Fisher’s exact tests were applied to identify significant GO and KEGG categories according to the absolute values of P < 0.05 and FDR < 0.05.
2.6 miRNA transcriptome sequencing and analysis
2.6.1 RNA extraction, quantification, and qualification
RNA extraction, quantification, and qualification were conducted as described in section 2.5.1.
2.6.2 Library construction, quality control, and sequencing
Briefly, 3’ and 5’ adaptors were ligated to the 3’ and 5’ ends of small RNAs, respectively. Next, first strand cDNA was synthesized after hybridization with the reverse transcription primer. The double-stranded cDNA library was generated through PCR enrichment. After purification and size selection, libraries with 18~40 bp insertions were selected for Illumina sequencing with SE50. The library was quantified with Qubit and real-time PCR and the library size distribution was evaluated with a Bioanalyzer. Quantified libraries were pooled and sequenced on an Illumina platform, according to the effective library concentration and amount of data required.
2.6.3 Data analysis
Raw data (fastq format) were first processed using in-house perl and python scripts. In this step, low-quality reads; reads containing poly-Ns, 5’ adapter sequences, or poly-As/Ts/Gs/Cs; and reads missing 3’ adapter sequences or insert tags were removed. Subsequently, the Q20, Q30, and GC content of the clean reads were calculated. All downstream analyses utilized only clean, high-quality data. Small RNA tags were mapped to the reference sequence using Bowtie (Langmead et al., 2009), either without mismatches or with only one mismatch, to analyze their expression and distribution. Mapped small RNA tags were used to identify known miRNAs, with miRBase (v22.0) used as a reference. miRNAs were identified using a modified version of mirdeep2 (Friedlander et al., 2012) and srna-tools-cli was used to draw the secondary structures. Custom scripts were used to obtain miRNA counts as well as to determine base bias at the first position of identified miRNAs of a certain length and at each position of all identified miRNAs. To remove tags originating from protein-coding genes, repeat sequences, rRNA, tRNA, snRNA, and snoRNA, small RNA tags were mapped to RepeatMasker, the Rfam database, or species-specific data. DEGs were identified using DESeq2. P-values were adjusted using the Benjamini & Hochberg method. A corrected P-value of 0.05 was selected as the threshold for determining significantly differential expression.
2.7 Analysis of transcriptional signatures of phytohormone responses
To identify the transcriptional signatures of BPH-responsive phytohormone responses, Hormonometer was used to compare gene expression in rice with gene expression in phytohormone-treated Arabidopsis thaliana (Volodarsky et al., 2009; Liu et al., 2016). Specifically, only orthologous genes detected in the RNA-seq analysis which were related to the Arabidopsis thaliana probe set identifiers were selected for further analyses.
2.8 qRT–PCR validation of DEGs
First-strand cDNA was synthesized using a PrimeScript RT Reagent Kit with gDNA Eraser (Takara, Japan), according to the manufacturer’s instructions. qRT-PCR assays of candidate genes were conducted using a PrimeScript™ RT reagent Kit with gDNA Eraser (Perfect Real Time) (Takara, RR047A, China). qRT-PCR was carried out on a CFX96 Real-Time System (Bio-Rad, CA, USA) according to the following protocol: 95 °C for 5 min, followed by 40 cycles at 95 °C for 5 s, 60 °C for 30 s, and 72 °C for 30 s. Relative gene expression was calculated with the 2−11Ct method, using PP2A as the reference gene. All primer sequences are listed in Supplementary Table 1.
2.9 Callose staining and evaluation
Callose staining and evaluation were performed as described in a previous study (Hu et al., 2017). Briefly, fresh sheaths collected from two-leaf stage rice seedlings infested with either biotype 1 or Y BPH for 48 h were fixed in an ethanol:acetic acid (3:1 v/v) solution for 5 h. The fixative was changed frequently to ensure thorough fixing and clearing. The samples were then rehydrated successively in 70% ethanol for 2 h, in 50% ethanol for 2 h, and in water overnight. After rinsing three times with water, the samples were treated with 10% NaOH for 1 h to make the tissues transparent. After rinsing four times with water, the samples were incubated in 150 mM K2HPO4 (pH 9.5) containing 0.01% aniline blue for 4 h. Finally, the samples were mounted on a slide and callose deposits were observed with a positive fluorescence microscope (Nikon, Eclipse 80i, Japan) under the UV channel.
2.10 Statistical analysis
Statistical analyses were conducted using R (v4.0.4) and SPSS (v22.0) (IBM SPSS, Somers, NY, USA). Two-sided Student’s t-tests were used to determine statistically significant differences between groups. All bioinformatics analyses were conducted using R packages.
3 Results
3.1 YHY15 exhibits differential resistance to biotype 1 and biotype Y BPH
YHY15 rice, which contains the BPH resistance gene Bph15, exhibits robust resistance against avirulent BPH biotypes (i.e., biotype 1) (Yang et al., 2004). In this study, YHY15 seedlings were subjected to infestation by either biotype 1 or biotype Y BPH. The results suggest that YHY15 exhibited distinct responses under the two infestation scenarios. In response to infestation with biotype Y BPH, YHY15 seedlings exhibited signs of withering at 4 days and eventually wilted completely by 7 days. In contrast, seedlings infested with biotype 1 BPH remained healthy and continued to grow vigorously (Figure 1A). To evaluate the ability of YHY15 rice to affect the biology of BPH insects, we measured the weight gain and amount of honeydew produced by the two BPH populations. As expected, biotype Y exhibited a significantly higher weight gain ratio (mean = 41.4%) than biotype 1 (mean = −2.2%) (Figure 1B). In addition, biotype Y produced significantly more honeydew (mean = 16.49 mg) than biotype 1 (mean = 0.57 mg) (Figure 1C). These findings indicate that YHY15 rice plants are highly resistant to biotype 1 BPH but are susceptible to biotype Y BPH, as previously reported (Yang et al., 2004; Jing et al., 2011; Guan et al., 2022). Overall, these results provide strong evidence for the contrasting responses of YHY15 rice to BPH with different levels of virulence, and highlight the efficacy of the Bph15 gene in conferring resistance against specific BPH strains.
Figure 1 Resistance of YHY15 rice plants against biotype 1 and biotype Y BPH. (A) YHY15 seedlings infested with biotype 1 or biotype Y BPH. DAI, days after infestation. Scale bar = 5 cm. (B) BPH weight gain ratio on YHY15 seedlings. (C) BPH honeydew excretion on YHY15 seedlings. Data are shown as means ± SD of 38 biological replicates in (B) and (C). Biotype 1, biotype 1 BPH-infested YHY15 seedlings; Biotype Y, biotype Y BPH-infested YHY15 seedlings. Asterisks indicate statistically significant differences between YY and TY (two-tailed Student’s t-test, **P < 0.01).
3.2 Overview of the miRNA- and RNA-seq results
To better understand the mechanisms underlying the differential resistance displayed by YHY15 rice to virulent and avirulent BPH, we performed RNA and miRNA sequencing analyses using leaf sheaths from YHY15 seedlings infested with either biotype 1 (RT) or Y (RY) BPH for either 6 (RT6, RY6) or 48 h (RT48, RY48). Un-infested rice plants were used as controls and named as R0. After constructing and sequencing the mRNA/miRNA libraries, high-quality raw sequence reads were normalized and subjected to further analysis.
Out of the 738 identified miRNAs, miR396, miR167, miR166, miR162, miR159, miR156, miR820, miR408, miR1425, miR1862, miR444, and miR827 exhibited the highest relative abundance (Supplementary Table 2). Principal component analysis (PCA) of the miRNA-seq data revealed significant variation between the R0 group and the RY6 or RY48 group, indicating that biotype Y BPH infestation led to distinct fluctuations in the miRNA profiles (Figure 2A). Among the 50 differentially expressed miRNAs (DEMs) detected in R0 vs RY6, 13 were upregulated and 37 were downregulated (Figure 2B; Supplementary Table 3). As infection progressed, more DEMs (64) were identified in R0 vs RY48, 16 of which were upregulated and 48 were downregulated (Figure 2B; Supplementary Table 3). In total, 38 DEMs were shared between R0 vs RY6 and R0 vs RY48, with 12 unique DEMs detected at the earlier time point and 26 at the later time point (Figure 2C). These results suggest that many more miRNAs are involved in the interaction between YHY15 and biotype Y BPH than in the interaction between YHY15 and biotype 1 BPH. Among these, miR156, miR5076, miR1856, miR398, miR5072, miR5079, miR408, miR2873, and miR169 exhibited significant variation in amplitude (Table 1). Many of these miRNAs have been reported to play specific roles in plant developmental processes and biotic and abiotic stress responses (Sharma et al., 2015; Zhao et al., 2017; Lin et al., 2018; Liu et al., 2019; Liebsch and Palatnik, 2020; Li et al., 2020; Liu X. et al., 2021; Gao et al., 2022; Pachamuthu and Hari Sundar, 2022; Zhao et al., 2022). Their high abundance and variable expression patterns suggest that they may contribute to BPH resistance in YHY15 seedlings, and therefore subsequent analyses were conducted on these miRNAs.
Figure 2 Overview of miRNA-seq and RNA-seq results. (A) Principal component analysis (PCA) of miRNA-seq data from five comparisons. (B) Number of miRNAs up- or downregulated in all comparisons (∣log2 fold change∣ > 1, P < 0.05). (C) Venn diagrams of differentially expressed miRNAs in all comparisons. (D) PCA of RNA-seq data from five comparisons. (E) Number of mRNAs up- or downregulated in all comparisons (|log2 fold change| > 1, P < 0.05). (F) Venn diagrams of differentially expressed mRNAs in all comparisons. There are four comparisons: RT6/R0, RT48/R0, RY6/R0, RY48/R0. R0, uninfected controls; RT6, YHY15 seedlings infested with biotype 1 BPH for 6 h; RT48, YHY15 seedlings infested with biotype 1 BPH for 48 h; RY6, YHY15 seedlings infested with biotype Y BPH for 6 h; RY48, YHY15 seedlings infested with biotype Y BPH for 48 h.
Table 1 Candidate BPH resistance-related DEMs.
In comparison, PCA of the RNA-seq data revealed distinct separation among all groups (Figure 2D; Supplementary Table 4). Through comparison of RNA expression levels across different groups, we identified 1027 and 1154 differential expression analysis (DEGs) in R0 vs RT6 and R0 vs RT48, respectively. Notably, a larger number of DEGs were observed in R0 vs RY6 and R0 vs RY48 (6817 and 6273, respectively), suggesting that infestation with biotype Y BPH results in more severe effects than infestation with biotype 1 BPH (Figure 2E). A Venn diagram was used to compare the expression patterns of DEGs in avirulent/virulent BPH-infested rice at the early (6 h) and late (48 h) stages of infestation, and the identified DEGs were subjected to GO and KEGG enrichment analyses (Figure 2F; Supplementary Figure 1; Supplementary Table 5). We identified 730 overlapping DEGs in R0 vs RT6 and R0 vs RY6 at the early infestation stage, which were enriched in cytochrome P450 (KEGG), protein modification-related biological processes (BP, GO), extracellular region (CC, GO), and monooxygenase activity (MF, GO) (Figure 2F; Supplementary Figures 1A, B). Specifically, genes in R0 vs RY6 were enriched in ribosome, photosynthesis proteins, carbon fixation in photosynthetic organisms, and DNA replication proteins (Figure 2F; Supplementary Figures 2A, B). However, genes in R0 vs RT6 did not exhibit these enrichments. We identified 827 overlapping DEGs in R0 vs RT48 and R0 vs RY48 at the late infestation stage, which were enriched in photosynthesis proteins, porphyrin metabolism, zeatin biosynthesis, and cytochrome P450, suggesting disturbance to the photosynthetic system (KEGG) (Figure 2F; Supplementary Figures 3A, B). Additionally, genes in R0 vs RY48 were enriched in ribosome-related pathways, DNA replication proteins, photosynthesis proteins, and carbon fixation in photosynthetic organisms (Figure 2F; Supplementary Figures 4A, B). These results indicate that infestation with biotype Y BPH results in damage to the photosynthetic system earlier than infestation with biotype 1 BPH, which is consistent with the phenotypic observations (Figure 1A).
3.3 Transcriptional responses in YHY15 rice to BPH infestation
To elucidate the genome-wide responses of YHY15 rice to BPH infestation, the transcriptome was analyzed over the entire infection time course. According to Mfuzz analysis, distinct temporal patterns were clustered into 20 groups (Figure 3; Supplementary Table 6). Notably, clusters 1 to 4 exhibited significantly lower levels of transcription in RY6/48 compared to R0, whereas this trend was less pronounced in RT6/48 (Figure 3; Supplementary Table 6). The genes in these clusters were subjected to GO and KEGG enrichment analyses to explore their functions. Genes in cluster 1 were enriched in photosynthesis proteins, carbon fixation in photosynthetic organisms, terpenoid backbone biosynthesis, and porphyrin metabolism (Figure 4A; Supplementary Table 7). Accordingly, GO analysis confirmed that many of these genes were related to the photosynthesis biological process, as well as photosynthetic cellular components such as thylakoid, photosynthetic membrane, and photosystem (Supplementary Figure 5A; Supplementary Table 7). These results suggest that infestation with biotype Y BPH results in significantly greater physical impairment than does infestation with biotype 1 BPH, which was consistent with the growth observations (Figure 1A).
Figure 3 Clustering and time-course expression of genes after BPH infestation. The 20 distinct temporal gene expression patterns were computed with Mfuzz. The x axis represents the five treatment groups, and the y axis represents log2-transformed normalized intensity ratios for each group. R0, uninfected controls; RT6, YHY15 seedlings infested with biotype 1 BPH for 6 h; RT48, YHY15 seedlings infested with biotype 1 BPH for 48 h; RY6, YHY15 seedlings infested with biotype Y BPH for 6 h; RY48, YHY15 seedlings infested with biotype Y BPH for 48 h.
Figure 4 KEGG pathway enrichment analyses of representative clusters. (A–F) KEGG pathway enrichment analyses of cluster 1, 6, 7, 10, 11, and 17, respectively.
In clusters 5 to 10, gene expression levels increased significantly in the RY groups, but only slightly in the RT groups (Figure 3; Supplementary Table 6). Genes in cluster 6 (Figure 4B; Supplementary Figure 5B; Supplementary Table 7) and cluster 9 (data not shown) were enriched in membrane trafficking (KEGG), cell membrane systems (GO), and protein transporter (GO), suggesting that defense against BPH relies on membrane transport systems and cell secretion. In cluster 7, genes were enriched in glycosyltransferases, GTP-binding proteins, phenylpropanoid biosynthesis, and two-component system (Figure 4C; Supplementary Figure 5C; Supplementary Table 7). The phenylpropanoid pathway is involved in the production of various metabolites, including flavonoids, lignin, lignans, and cinnamic acid amide, among others (Dong and Lin, 2021). Genes in cluster 10 were enriched in ribosomes, DNA replication proteins, translation factors, and other transcriptional activity-related factors (Figure 4D; Supplementary Figure 5D; Supplementary Table 7). Overall, the transcriptional system, transmembrane transport system, exosome, and phenylpropanoid pathway are notably stimulated by BPH infestation, particularly infestation with biotype Y BPH, suggesting that they may be crucial for rice resistance to insects.
The expression patterns of several clusters appeared to be desynchronized, indicative of different responses to infestation by either biotype 1 or biotype Y BPH. In clusters 11 to 14, gene expression was increased at the early (6 h) and late (48 h) stages during infestation with biotype 1 BPH, but remained relatively stable during infestation with biotype Y BPH (Figure 3; Supplementary Table 6). These clusters were mainly enriched in zeatin biosynthesis, nitrogen metabolism, cytochrome P450, and plant hormone signal transduction (Figure 4E; Supplementary Figure 5E; Supplementary Table 7). In plants, cytochrome P450 catalyzes various primary and secondary metabolic reactions and is involved in the synthesis and metabolism of terpenoids, alkaloids, sterols, fatty acids, plant hormones, signaling molecules, pigments, flavonoids, and isoflavones, among others (Hansen et al., 2021). Such multifunctionality makes cytochrome P450 important in the plant defense against pests, diseases, and abiotic stressors. In addition, zeatin is the primary active component of the phytohormone CK. KEGG analysis revealed that the cytochrome P450-phytohormone signal transduction pathway was activated in YHY15 rice under biotype 1 BPH infestation. In contrast, gene expression in clusters 15 to 17 was increased at the early (6 h) and late (48 h) stages during infestation with biotype Y BPH (Figure 3; Supplementary Table 6). These genes were enriched in proteasome, membrane trafficking, steroid biosynthesis, protein processing in endoplasmic reticulum, exosome, and GTP-binding proteins (Figure 4F; Supplementary Figure 5F; Supplementary Table 7). These results suggest that cell secretion and extracellular materials are crucial for defense against insects. Finally, clusters 18 to 20 did not exhibit any notable tendencies (Figure 3).
3.4 Phytohormonal responses in YHY15 rice to BPH infestation
Both the clustering analysis and KEGG/GO enrichment analysis indicated that phytohormone signaling was differentially affected by infestation with either biotype 1 or Y BPH (Figures 3, 4E). Consequently, we compared the expression of BPH-responsive genes in rice with those induced by exogenous phytohormone application in A. thaliana. According to Hormonometer analysis, a total of 15,603 A. thaliana orthologs of rice genes were selected for comparison (Supplementary Table 8). The results revealed that genes associated with JA-dependent signaling were significantly altered by infestation with biotype 1 BPH, but only moderately changed by infestation with biotype Y BPH (Figure 5). Genes associated with ABA exhibited a similar pattern, except for being negatively correlated with genes whose expression in A. thaliana was elicited 0.5 h after ABA treatment (Figure 5). In contrast, the expression of genes involved in the ET and BR pathways exhibited negative correlations with those responding to phytohormone application in A. thaliana, although with more moderate suppression in biotype Y BPH-infested rice (Figure 5). These results show that genes associated with the JA, ABA, ET, and BR pathways are the most responsive to BPH infestation, particularly in the case of biotype 1 BPH, indicating their crucial roles in the rice defense against insects.
Figure 5 Comparison of transcriptomic phytohormone signatures between BPH-infested rice and phytohormone-treated Arabidopsis thaliana. Red shading indicates positive correlations and blue shading indicates negative correlations between BPH-infested rice and phytohormone-treated A. thaliana. MJ, methyl jasmonate; ACC, 1-aminocyclopropane-1-caroxylic acid (precursor of ethylene); ABA, abscisic acid; IAA, indole-3-acetic acid; GA3, gibberellic acid 3; ZT, zeatin; BR, brassinosteroid; SA, salicylic acid. R0, uninfected controls; RT6, YHY15 seedlings infested with biotype 1 BPH for 6 h; RT48, YHY15 seedlings infested with biotype 1 BPH for 48 h; RY6, YHY15 seedlings infested with biotype Y BPH for 6 h; RY48, YHY15 seedlings infested with biotype Y BPH for 48 h.
Previous research suggests that JA is involved in BPH resistance (Xu et al., 2021). We found that the JA-dependent pathway was differentially induced by infestation with either biotype 1 or Y BPH, implying that JA plays a crucial role in rice defense (Figure 5). We further analyzed the expression of genes involved in the JA pathway, and found that most of the genes involved in the biosynthesis of JA and jasmonoyl-isoleucine (JA-Ile) (e.g., LOX and JAR) were significantly suppressed by BPH infestation, especially infestation with biotype Y BPH (Figures 6A, B). Accordingly, most JAZs, which inhibit JA signaling, were induced in RT6/48 and RY6/48 (Figures 6A, B). In addition, MYC2 was downregulated in YHY15 rice infested with biotype Y BPH, suggesting that the JA-dependent pathway was inactivated upon infestation with the virulent BPH (Ge et al., 2018). Finally, the expression of many TFs was altered in BPH-infested rice, especially rice infested with biotype Y BPH, including those belonging to the JA-responsive bHLH, ERF, WRKY, and MYB families (Supplementary Figure 6). These distinct response patterns imply that these TFs have functions in BPH defense, although further research is required to verify these functions. The expression profiles of the identified TFs can be found in Supplementary Table 9.
Figure 6 BPH-induced responses in the jasmonic acid (JA) pathway. (A) Overview of the JA pathway. (B) Heat map of the expression of genes associated with the JA pathway. Asterisks indicate statistically significant differences in gene expression at different time points in either biotype 1 or biotype Y BPH-infested rice relative to control (0 h) (*padj < 0.05, **padj < 0.01, via the Benjamini and Hochberg adjustment method). PLA, phospholipase A1; LOX, lipoxygenase; AOS, allene oxide synthase; AOC, allene oxide cyclase; OPR, 12-oxophytodienoate reductase; OPCL, OPC8-CoA ligase; ACX, acyl-CoA oxidase; MFP, multifunctional protein; KAT, 3-ketoacyl-CoA thiolase; JAR, jasmonate resistant; JAZ, jasmonate-ZIM domain. Transcription factor families: myelocytomatosis protein 2 (MYC2), basic helix-loop-helix (bHLH), ethylene responsive factor (ERF), WRKYGOK (WRKY), MYB. R0, uninfected controls; RT6, YHY15 seedlings infested with biotype 1 BPH for 6 h; RT48, YHY15 seedlings infested with biotype 1 BPH for 48 h; RY6, YHY15 seedlings infested with biotype Y BPH for 6 h; RY48, YHY15 seedlings infested with biotype Y BPH for 48 h.
3.5 Integrated miRNA and mRNA transcriptional analyses
The functions of the identified miRNAs were initially evaluated by scanning the literature. It has been reported that miR156 negatively regulates BPH resistance in rice by altering the expression of genes related to JA biosynthesis and signaling (Ge et al., 2018). We observed that miR156 expression was significantly decreased in biotype Y BPH-infested rice. Consequently, WRKY53, several MPKs (JA biosynthesis repressors), and several JAZs (negative regulators of JA signaling) were upregulated (Table 2). These transcriptional alterations resulted in low JA expression and improved BPH resistance. According to other research, miR396 suppresses BPH resistance through the miR396b–growth regulating factor 8 (GRF8)–flavanone 3-dioxygenase (F3H)–flavonoid pathway (Dai et al., 2019). Therefore, we examined the expression of miR396, GRF8, and F3H. Mature OsmiR396a and OsmiR396b share the same sequence (Dai et al., 2019), and we observed that miR396a was downregulated by infestation with biotype Y BPH. In contrast, GRF8 and F3H-1 expression increased in RY6 and RY48. Notably, GRF1 and two flavanone 3-dioxygenase 2-like genes exhibited the same trend (Table 2). These results suggest that the miR396b–GRF8–F3H–flavonoid pathway was activated in biotype Y BPH- infested rice. Genetic experiments indicate that miR398b negatively regulates pathogen-associated molecular pattern (PAMP)-induced callose deposition (Li et al., 2010). As miR398 expression was significantly downregulated in RY6/48, we further analyzed the transcription of callose deposition-related genes. Surprisingly, a callose synthase gene was activated following biotype Y BPH infestation, while a series of callose degradation genes were suppressed (Table 2). These results suggest that infestation with biotype Y BPH may lead to callose deposition in rice. Taken together, it appears that the expression of these miRNAs and their regulation of target genes may contribute the ability of YHY15 rice to resist infestation with biotype Y BPH.
Table 2 The miRNA-mRNA interactions related to plant resistance.
3.6 Callose deposition is activated by BPH infestation
The combined miRNA-seq and RNA-seq analysis indicated that callose deposition might play a vital role in the rice response to BPH infestation. Callose deposition prevents plant hoppers from ingesting phloem sap (Hao et al., 2008; Du et al., 2009). Specifically, callose deposition ensures that the phloem sieve tubes remain occluded. To determine the role of callose in the rice response to BPH infestation, the expression of genes related to callose synthesis and degradation was examined. Consistent with the RNA-seq results, the transcription of callose synthase 3-like (CAL1) increased, while several 1,3-beta-glucanase (BG) genes decreased, in RY48 (Figures 7A, B). To confirm this result, callose deposition was measured during BPH infestation. The outermost sheaths of biotype Y BPH-infested seedlings exhibited larger and more numerous callose spots than those from biotype 1 BPH-infested and uninfested seedlings (Figures 7C, D). These results suggest that BPH infestation, particularly infestation with biotype Y BPH, results in the activation of callose synthesis and the suppression of callose degradation.
Figure 7 Induction of callose deposition in YHY15 rice plants following BPH infestation. (A) Heat map of the expression of genes associated with callose deposition and degradation. Asterisks indicate statistically significant differences in gene expression at different time points in either biotype 1 or biotype Y BPH-infested rice relative to control (0 h) (*padj < 0.05, **padj < 0.01, via the Benjamini and Hochberg adjustment method). (B) Expression patterns of callose deposition-related genes in BPH-infested rice plants. UBQ was used as a control. (C) Callose deposition (CD) in BPH-infested YHY15 plants. The images were taken 48 h after BPH infestation. Samples were collected from the outermost sheath. Callose deposition values are the means of 20 biological replicates. Asterisks indicate statistically significant differences between BPH-infested and uninfested rice plants (one-way ANOVA, *P < 0.05). (D) Callose deposition in YHY15 plants infested with either biotype 1 or biotype Y BPH. Scale bar = 20 μm. The images were taken 48 h after BPH infestation. Uninfested YHY15 plants in the same condition served as control. R0, uninfected controls; RT6, YHY15 seedlings infested with biotype 1 BPH for 6 h; RT48, YHY15 seedlings infested with biotype 1 BPH for 48 h; RY6, YHY15 seedlings infested with biotype Y BPH for 6 h; RY48, YHY15 seedlings infested with biotype Y BPH for 48 h.
4 Discussion
The management of rice plant resistance to BPH infestation will be crucial for effective pest management. Understanding the defense strategies employed by resistant rice varieties against both avirulent and virulent BPH populations can provide valuable insights for developing effective pest control strategies. Transcriptomics, which can quantify changes in gene expression and associated regulatory mechanisms, can aid in unraveling these complexities. In plants, miRNAs are involved in various developmental processes and play significant roles in abiotic and biotic stress responses (Zhang et al., 2013; Kumar, 2014; Zhang et al., 2016; Li et al., 2017; Mangrauthia et al., 2017; Natarajan et al., 2018). Specifically, miRNAs regulate targeted gene expression by binding complementary sequences in mRNA molecules, resulting in degradation and/or inhibited translation (Bartel, 2009; Marin et al., 2010; Peng et al., 2014). Integrated mRNA and miRNA transcriptomics analyses have been used to identify miRNA-mRNA networks associated with developmental processes and insect defense (Tan et al., 2020; Li et al., 2021; Mei et al., 2021). However, little attention has been paid to defense strategies employed by resistant rice against differentially virulent BPH populations. In this study, we employed an integrated mRNA and miRNA transcriptomics approach to characterize the defense responses of resistant YHY15 rice (contains Bph15) to infestation with avirulent (biotype 1) and virulent (biotype Y) BPH. Our results revealed that YHY15 rice seedlings exhibited distinct responses under the two infestation scenarios (Figure 1). YHY15 rice is highly resistant to biotype 1 BPH and susceptible to biotype Y BPH, which is in accordance with previous reports (Yang et al., 2004; Jing et al., 2011; Guan et al., 2022).
The BPH-resistance gene Bph15 has been widely applied in rice breeding programs, although the molecular mechanisms underlying Bph15-mediated resistance remain unclear (Lv et al., 2014; Zheng et al., 2021). Previous RNA-seq studies of Bph15 introgression lines and recipient lines before and after BPH infestation have identified key defense mechanisms associated with phytohormone signaling, mitogen-activated protein kinase (MAPK) cascades, receptor kinases, protein post-translational modifications, TFs, Ca2+ signaling, and pathogenesis-related proteins (Lv et al., 2014). In addition, 20 upregulated and 3 downregulated miRNAs were identified in resistant rice variety P15 (containing Bph15) compared to susceptible rice variety PC (recipient line) (Wu et al., 2017). Combined with the mRNA transcriptome data, the 67 potential targets of these miRNAs were related to resistance responses to avirulent BPH, including abiotic and biotic stimuli, regulation of plant hormones (GA, SA, ET, and CK), cellulose biosynthesis, amino acid biosynthesis, and protein folding (Cheng et al., 2013; Lv et al., 2014; Wu et al., 2017). Here, comprehensive analysis of miRNA-seq and mRNA-seq data shed light on the underlying mechanisms responsible for the contrasting responses of YHY15 rice to BPH biotypes with different virulence.
Plant defense against BPH is a dynamic and sophisticated process which involves many levels of organizational and functional complexity (Barah and Bones, 2015; Erb and Reymond, 2019). During feeding, the BPH stylet transiently punctures the epidermis and then penetrates the plant cell wall. Subsequently, the insect salivates into the cells and ingests the phloem sap (Hao et al., 2008). According to electronic penetration graph (EPG) waveform recordings, BPH begin feeding on phloem sap within 1-3 hours of settling on rice plants (Hao et al., 2008). A greater number of miRNAs were found to be upregulated in biotype 1 BPH-infested resistant P15 rice (a Bph15 introgression line) than in susceptible PC rice (recipient line) during the early infestation stage (6 h) when the plants had not yet been severely damaged (Wu et al., 2017). In another report, the inducible BPH defense responses (indicated by upregulated DEGs and DEMs) were more robust during the early feeding stages (e.g., 6, 12, and 24 h) in resistant BPH6G rice (BPH6-transgenic rice) than in susceptible Nipponbare rice (wild type, WT) (Tan et al., 2020). Moreover, a miRNA profiling was conducted on resistant IR56 rice (carrying Bph3) under separate infestations of a virulent IR56-BPH and an avirulent TN1-BPH (Nanda et al., 2020). This study revealed that BPH feeding caused significant alterations in miRNA expression profiles of IR56 rice, with a greater number of miRNAs showing downregulation when IR56 rice was infested with TN1-BPH. However, the distinct mechanisms underlying rice plant responses to BPH of varying levels of virulence remains unclear. Here, resistant YHY15 rice plants were exposed to avirulent (biotype 1) and virulent (biotype Y) BPH. Notably, DEMs were only identified in rice plants infested with biotype Y BPH (Figure 2A). In addition, many more DEGs were identified in rice plants infested with biotype Y BPH, regardless of the infestation time (Figure 2E). Together, these results suggest that biotype Y BPH elicit more intense defense responses in YHY15 rice. DaEGs related with cytochrome P450 (KEGG) were more common at the early infestation stage (6 h) in rice infested with both types of BPH (Figure 2F; Supplementary Figures 1A, B). At the later infestation stage (48 h), overlapping DEGs in R0 vs RT48 and R0 vs RY48 were enriched in photosynthesis-related pathways (Figure 2F; Supplementary Figures 3A, B). Specifically, DEGs in R0 vs RY6 were mainly enriched in photosynthetic organisms and DNA replication proteins (Figure 2F; Supplementary Figures 2A, B), indicating that damage to the photosynthetic system occurs earlier in biotype Y BPH-infested rice. Taken together, compared with biotype 1 BPH, infestation with biotype Y BPH results in more serious damage and induces more intense transcriptional responses in YHY15 rice.
In plants, including rice, phytohormone signaling is widely known to be involved in insect defense and resistance (Liu et al., 2016; Jing et al., 2023; Liu Q. et al., 2021). Among phytohormones, JA is critical to the regulation of plant defenses against insect herbivores (Thaler et al., 2012; Howe et al., 2018). Other phytohormones, such as SA, ET, BR, ABA, and CK are also involved in plant responses to herbivory through cross-talk with JA signaling (Bruessow et al., 2010; Schäfer et al., 2015; Pan et al., 2018; Ma et al., 2020; Li et al., 2022; Zhang et al., 2022). In this study, we found that infestation with differentially virulent BPH populations resulted in noticeable effects on the expression of genes associated with phytohormone signaling in YHY15 rice plants. Specifically, infestation with biotype 1 BPH induced JA- and ABA-related signaling pathways, but suppressed ET- and BR-related signaling (Figure 5). While, infestation with biotype Y BPH induced more moderate expression of JA- and ABA-related pathways. These results suggest that although infestation with virulent biotype Y BPH results in extensive damage to YHY15 rice plants, they exhibit weaker JA- and ABA-mediated defense responses. These findings are in line with previous research suggesting that IR56 BPH can overcome Bph3-mediated resistance by suppressing the transcription of defense-responsive MAPK pathways, phytohormone biosynthesis, and secondary metabolite production (Nanda et al., 2018). Previous studies also suggest that oral secretions, digestive and detoxifying enzymes, and endosymbionts can help insect herbivores adapt to host plants (Simon et al., 2015; Yates and Michel, 2018). We therefore speculate that specific BPH effectors may interact with Bph15 and affect Bph15-mediated immunity in rice. However, further studies will be required to test this hypothesis.
The DEMs and their associated mRNAs identified in this study may potentially play roles in the response of rice to BPH infestation. Among the identified DEMs, miR156, miR5076, miR1856, miR398, miR5072, miR5079, miR408, miR2873, and miR169 exhibited significant variation in amplitude (Supplementary Table 3). In addition, many of them have been implicated in various plant developmental processes and stress responses (Li et al., 2017b) (Sharma et al., 2015; Zhao et al., 2017; Lin et al., 2018; Liu et al., 2019; Liebsch and Palatnik, 2020; Li et al., 2020; Liu X. et al., 2021; Gao et al., 2022; Pachamuthu and Hari Sundar, 2022; Zhao et al., 2022). The altered expression of these miRNAs may enhance BPH resistance in YHY15 rice. Combined with mRNA transcriptome data (Table 2), we found that the miR156-JA, miR396b–GRF8–F3H–flavonoid, and miR398b-callose deposition pathways may contribute to the resistance of YHY15 rice to biotype Y BPH. Physiological tests verified the increased deposition of callose in biotype Y BPH-infested rice (Figure 7). Notably, these miRNA-mediated responses only occurred in rice infested with biotype Y BPH, which suggests that they may contribute to the differential resistance of rice against BPH populations with varying levels of virulence.
Finally, research suggests that numerous TFs are involved in the rice response to BPH infestation. For example, OsMYB30 (an R2R3 MYB TF) induces the expression of phenylalanine ammonia-lyase (PAL) enzymes, thereby improving BPH resistance in rice (He et al., 2020). OsMYB22 promotes rice resistance by affecting flavonoid biosynthesis (Sun et al., 2023). The bHLH protein MYC2 is involved in JA-mediated insect resistance (Schweizer et al., 2013; Xu et al., 2021). OsWRKY45, OsWRKY53, OsWRKY70, and OsWRKY89 also mediate herbivore resistance (Chujo et al., 2007; Wang et al., 2007; Hu et al., 2015; Li et al., 2015; Huangfu et al., 2016). A recent study has demonstrated the pivotal involvement of OsWRKY71 in Bph15-mediated resistance (Li et al., 2023). Here, we identified several differentially expressed TFs associated with BPH infestation. Notably, several predominant TF families, including bHLH, MYB, ERF, WRKY, bZIP, NAC, C2H2, TALE, G2-like, HD-ZIP, MYB-related, HSF, and NF-Y, were differentially responsive to BPH infestation. The expression of most TFs was altered in biotype Y BPH-infested rice, indicating their specific roles in defense against the virulent biotype. Notably, the expression of certain other TFs was disturbed specifically in response to biotype 1 BPH infestation. These findings highlight the important functions of TFs in BPH defense and warrant further research to uncover their specific roles.
5 Conclusion
In conclusion, our study provides valuable insights into the differential defense strategies employed by resistant YHY15 rice (carrying BPH resistance gene Bph15) against avirulent (biotype 1) and virulent (biotype Y) BPH. The BPH defense response was found to involve the modulation of miRNAs, TFs, phytohormone signaling pathways, and the induction of callose deposition. These responses were most noticeable in biotype Y BPH-infested rice plants. These findings contribute to the elucidation of the molecular intricacies underlying rice-BPH interactions and pave the way for further research into the specific genes, pathways, and regulatory elements involved in plant defense against diverse BPH populations. Studying these defense mechanisms will aid our understanding of the intricate interactions between rice and BPH and allow the development of targeted pest control strategies for improved rice cultivation. It is worth noting that the infestation-induced defense responses of YHY15 rice do not appear to alter the survivability of biotype Y BPH, implying the existence of corresponding adaptive responses in the virulent biotype. This result will be the subject of in-depth exploration in future studies.
Data availability statement
All raw RNA sequencing data generated in this study have been deposited under the NCBI SRA database under BioProject PRJNA997052 and PRJNA994598.
Ethics statement
The manuscript presents research on animals that do not require ethical approval for their study.
Author contributions
BY: Writing – original draft, Writing – review & editing. MG: Investigation, Writing – original draft, Writing – review & editing. YX: Investigation, Writing – review & editing. QY: Writing – review & editing. BL: Formal analysis, Writing – review & editing. ML: Writing – review & editing. YS: Writing – review & editing. CL: Writing – review & editing. JX: Writing – review & editing. JL: Writing – review & editing. WH: Writing – review & editing. HT: Writing – review & editing. PL: Writing – review & editing. QL: Investigation, Supervision, Formal analysis, Writing – review & editing. SJ: Supervision, Project administration, Resources, Funding acquisition, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by grants from the National Natural Science Foundation of China (U1704111, 32372548 and 31401732), ZHONGYUAN YINGCAI JIHUA (ZYYCYU202012165), the Open Project Funding of the State Key Laboratory of Crop Stress Adaptation and Improvement (2023KF10), Young Elite Scientists Sponsorship Program by CAST (2023QNRC001), and the National Key Research and Development Program of China (2022YFD1401600).
Acknowledgments
We thank Prof. Guangcun He (Wuhan University) for kindly providing the insects and plants.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1366515/full#supplementary-material
References
Alam, S. N., Cohen, M. B. (1998). Durability of brown planthopper, Nilaparvata lugens, resistance in rice variety IR64 in greenhouse selection studies. Entomol. Exp. Appl. 89, 71–78. doi: 10.1046/j.1570-7458.1998.00383.x
Axtell, M. J., Meyers, B. C. (2018). Revisiting criteria for plant microRNA annotation in the Era of big data. Plant Cell 30, 272–284. doi: 10.1105/tpc.17.00851
Barah, P., Bones, A. M. (2015). Multidimensional approaches for studying plant defence against insects: from ecology to omics and synthetic biology. J. Exp. Bot. 66, 479–493. doi: 10.1093/jxb/eru489
Bartel, D. P. (2009). MicroRNAs: target recognition and regulatory functions. Cell 136, 215–233. doi: 10.1016/j.cell.2009.01.002
Bottrell, D. G., Schoenly, K. G. (2012). Resurrecting the ghost of green revolutions past: The brown planthopper as a recurring threat to high-yielding rice production in tropical Asia. J. Asia Pac. Entomol. 15, 122–140. doi: 10.1016/j.aspen.2011.09.004
Bruessow, F., Gouhier-Darimont, C., Buchala, A., Metraux, J. P., Reymond, P. (2010). Insect eggs suppress plant defence against chewing herbivores. Plant J. 62, 876–885. doi: 10.1111/tpj.2010.62.issue-5
Cheng, X. Y., Zhu, L. L., He, G. C. (2013). Towards understanding of molecular interactions between rice and the brown planthopper. Mol. Plant 6, 621–634. doi: 10.1093/mp/sst030
Chujo, T., Takai, R., Akimoto-Tomiyama, C., Ando, S., Minami, E., Nagamura, Y., et al. (2007). Involvement of the elicitor-induced gene OsWRKY53 in the expression of defense-related genes in rice. Biochim. Biophys. Acta 1769, 497–505. doi: 10.1016/j.bbaexp.2007.04.006
Dai, Z. Y., Tan, J., Zhou, C., Yang, X. F., Yang, F., Zhang, S. J., et al. (2019). The OsmiR396-OsGRF8-OsF3H-flavonoid pathway mediates resistance to the brown planthopper in rice (Oryza sativa). Plant Biotechnol. J. 17, 1657–1669. doi: 10.1111/pbi.13091
Dinh, S. T., Baldwin, I. T., Galis, I. (2013). The HERBIVORE ELICITOR-REGULATED1 gene enhances abscisic acid levels and defenses against herbivores in Nicotiana attenuata plants. Plant Physiol. 162, 2106–2124. doi: 10.1104/pp.113.221150
Dong, N. Q., Lin, H. X. (2021). Contribution of phenylpropanoid metabolism to plant development and plant-environment interactions. J. Integr. Plant Biol. 63, 180–209. doi: 10.1111/jipb.13054
Du, B., Zhang, W. L., Liu, B. F., Hu, J., Wei, Z., Shi, Z. Y., et al. (2009). Identification and characterization of Bph14, a gene conferring resistance to brown planthopper in rice. Proc. Natl. Acad. Sci. 106, 22163–22168. doi: 10.1073/pnas.0912139106
Erb, M., Reymond, P. (2019). Molecular interactions between plants and insect herbivores. Annu. Rev. Plant Biol. 70, 527–557. doi: 10.1146/annurev-arplant-050718-095910
Friedlander, M. R., Mackowiak, S. D., Li, N., Chen, W., Rajewsky, N. (2012). miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 40, 37–52. doi: 10.1093/nar/gkr688
Gao, Y., Feng, B. H., Gao, C. X., Zhang, H. Q., Wen, F. T., Tao, L. X., et al. (2022). The evolution and functional roles of miR408 and its targets in plants. Int. J. Mol. Sci. 23, 530–530. doi: 10.3390/ijms23010530
Ge, Y. F., Han, J. Y., Zhou, G. X., Xu, Y. M., Ding, Y., Shi, M., et al. (2018). Silencing of miR156 confers enhanced resistance to brown planthopper in rice. Planta 248, 813–826. doi: 10.1007/s00425-018-2942-6
Guan, W., Shan, J. H., Gao, M. Y., Guo, J. P., Wu, D., Zhang, Q., et al. (2022). Bulked segregant RNA sequencing revealed difference between virulent and avirulent brown planthoppers. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.843227
Guo, J. P., Xu, C. X., Wu, D., Zhao, Y., Qiu, Y. F., Wang, X. X., et al. (2018). Bph6 encodes an exocyst-localized protein and confers broad resistance to planthoppers in rice. Nat. Genet. 50, 297–306. doi: 10.1038/s41588-018-0039-6
Gutaker, R. M., Groen, S. C., Bellis, E. S., Choi, J. Y., Pires, I. S., Bocinsky, R. K., et al. (2020). Genomic history and ecology of the geographic spread of rice. Nat. Plants 6, 492–502. doi: 10.1038/s41477-020-0659-6
Hansen, C. C., Nelson, D. R., Møller, B. L., Werck-Reichhart, D. (2021). Plant cytochrome P450 plasticity and evolution. Mol. Plant 14, 1244–1265. doi: 10.1016/j.molp.2021.06.028
Hao, P. Y., Liu, C. X., Wang, Y. Y., Chen, R. Z., Tang, M., Du, B., et al. (2008). Herbivore-induced callose deposition on the sieve plates of rice: an important mechanism for host resistance. Plant Physiol. 146, 1810–1820. doi: 10.1104/pp.107.111484
He, J., Liu, Y. Q., Yuan, D. Y., Duan, M. J., Liu, Y. L., Shen, Z. J., et al. (2020). 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. 117, 271–277. doi: 10.1073/pnas.1902771116
Howe, G. A., Major, I. T., Koo, A. J. (2018). Modularity in jasmonate signaling for multistress resilience. Annu. Rev. Plant Biol. 69, 387–415. doi: 10.1146/annurev-arplant-042817-040047
Hu, L., Wu, Y., Wu, D., Rao, W. W., Guo, J. P., Ma, Y. H., et al. (2017). The coiled-coil and nucleotide binding domains of BROWN PLANTHOPPER RESISTANCE14 function in signaling and resistance against planthopper in rice. Plant Cell. 29, 3157–3185. doi: 10.1105/tpc.17.00263
Hu, L. F., Ye, M., Li, R., Zhang, T. F., Zhang, G. X., Wang, Q., et al. (2015). The rice transcription factor WRKY53 suppresses herbivore-induced defenses by acting as a negative feedback modulator of mitogen-activated protein kinase activity. Plant Physiol. 169, 2907–2921. doi: 10.1104/pp.15.01090
Huangfu, J. Y., Li, J. C., Li, R., Ye, M., Kuai, P., Zhang, T. F., et al. (2016). The transcription factor OsWRKY45 negatively modulates the resistance of rice to the brown planthopper Nilaparvata lugens. Int. J. Mol. Sci. 17, 697. doi: 10.3390/ijms17060697
Ji, R., Yu, H. X., Fu, Q., Chen, H. D., Ye, W. F., Li, S. H., et al. (2013). Comparative transcriptome analysis of salivary glands of two populations of rice brown planthopper, Nilaparvata lugens, that differ in virulence. PloS One 8, e79612. doi: 10.1371/journal.pone.0079612
Jing, S. L., Liu, B., Peng, L., Peng, X., Zhu, L., Fu, Q., et al. (2011). Development and use of EST-SSR markers for assessing genetic diversity in the brown planthopper (Nilaparvata lugens Stål). Bull. Entomol. Res. 102, 113–122. doi: 10.1017/S0007485311000435
Jing, S. L., Xu, J. G., Tang, H. M., Li, P., Yu, B., Liu, Q. S. (2023). The roles of small RNAs in rice-brown planthopper interactions. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1326726
Jing, S. L., Zhao, Y., Du, B., Chen, R. Z., Zhu, L. L., He, G. C. (2017). Genomics of interaction between the brown planthopper and rice. Curr. Opin. Insect Sci. 19, 82–87. doi: 10.1016/j.cois.2017.03.005
Kumar, R. (2014). Role of microRNAs in biotic and abiotic stress responses in crop plants. Appl. Biochem. Biotechnol. 174, 93–115. doi: 10.1007/s12010-014-0914-2
Kumar, L., Futschik, M. E. (2007). Mfuzz: a software package for soft clustering of microarray data. Bioinformation 2, 5–7. doi: 10.6026/bioinformation
Langmead, B., Trapnell, C., Pop, M., Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25. doi: 10.1186/gb-2009-10-3-r25
Li, J. T., Chen, L., Ding, X., Fan, W. Y., Liu, J. L. (2022). Transcriptome Analysis reveals crosstalk between the abscisic acid and jasmonic acid signaling pathways in rice-mediated defense against. Nilaparvata lugens. Int. J. Mol. Sci. 23, 6319–6319. doi: 10.3390/ijms23116319
Li, C. Y., Luo, C., Zhou, Z. H., Wang, R., Ling, F., Xiao, L. T., et al. (2017). Gene expression and plant hormone levels in two contrasting rice genotypes responding to brown planthopper infestation. BMC Plant Biol. 17, 57. doi: 10.1186/s12870-017-1005-7
Li, X. P., Ma, X. C., Wang, H., Zhu, Y., Liu, X. X., Li, T. T., et al. (2020). Osa-miR162a fine-tunes rice resistance to Magnaporthe oryzae and Yield. Rice (New York N.Y.). 13, 38. doi: 10.1186/s12284-020-00396-2
Li, H. Y., Meng, H. L., Sun, X. Q., Deng, J., Shi, T. X., Zhu, L. W., et al. (2021). Integrated microRNA and transcriptome profiling reveal key miRNA-mRNA interaction pairs associated with seed development in Tartary buckwheat (Fagopyrum tataricum). BMC Plant Biol. 21, 132. doi: 10.1186/s12870-021-02914-w
Li, R., Zhang, J., Li, J. C., Zhou, G. X., Wang, Q., Bian, W. B., et al. (2015). Prioritizing plant defence over growth through WRKY regulation facilitates infestation by non-target herbivores. eLife 4, e04805. doi: 10.7554/eLife.04805.021
Li, X. Z., Zhang, J., Shangguan, X. X., Yin, J. J., Zhu, L. L., Hu, J., et al. (2023). Knockout of OsWRKY71 impairs Bph15-mediated resistance against brown planthopper in rice. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1260526
Li, Y., Zhang, Q. Q., Zhang, J. G., Wu, L., Qi, Y. J., Zhou, J. M. (2010). Identification of MicroRNAs involved in pathogen-associated molecular pattern-triggered plant innate immunity. Plant Physiol. 152, 2222–2231. doi: 10.1104/pp.109.151803
Li, Y., Zhao, S. L., Li, J. L., Hu, X. H., Wang, H., Cao, X. L., et al. (2017b). Osa-miR169 negatively regulates rice immunity against the blast fungus Magnaporthe oryzae. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.00002
Liebsch, D., Palatnik, J. F. (2020). MicroRNA miR396, GRF transcription factors and GIF co-regulators: a conserved plant growth regulatory module with potential for breeding and biotechnology. Curr. Opin. Plant Biol. 53, 31–42. doi: 10.1016/j.pbi.2019.09.008
Lin, W. Y., Lin, Y. Y., Chiang, S. F., Syu, C., Hsieh, L. C., Chiou, T. J. (2018). Evolution of microRNA827 targeting in the plant kingdom. New Phytol. 217, 1712–1725. doi: 10.1111/nph.14938
Liu, J. L., Du, H. T., Ding, X., Zhou, Y. D., Xie, P. F., Wu, J. C. (2017). Mechanisms of callose deposition in rice regulated by exogenous abscisic acid and its involvement in rice resistance to Nilaparvata lugens Stål (Hemiptera: Delphacidae). Pest Manage. Sci. 73, 2559–2568. doi: 10.1002/ps.4655
Liu, Q. S., Hu, X. Y., Su, S. L., Ning, Y. S., Peng, Y. F., Ye, G. Y., et al. (2021). Cooperative herbivory between two important pests of rice. Nat. Commun. 12, 6772. doi: 10.1038/s41467-021-27021-0
Liu, X., Huang, S., Xie, H. T. (2021). Advances in the regulation of plant development and stress response by miR167. Front. Biosci. 26, 655–665. doi: 10.52586/4974
Liu, M. M., Shi, Z. Y., Zhang, X. H., Wang, M. X., Zhang, K. Z., Liu, J. Y., et al. (2019). Inducible overexpression of Ideal Plant Architecture1 improves both yield and disease resistance in rice. Nat. plants. 5, 389–400. doi: 10.1038/s41477-019-0383-2
Liu, Q. S., Wang, X. Y., Tzin, V., Romeis, J., Peng, Y. F., Li, Y. H. (2016). Combined transcriptome and metabolome analyses to understand the dynamic responses of rice plants to attack by the rice stem borer Chilo suppressalis (Lepidoptera: Crambidae). BMC Plant Biol. 16, 259. doi: 10.1186/s12870-016-0946-6
Lu, J., Ju, H. P., Zhou, G. X., Zhu, C. S., Erb, M., Wang, X. P., et al. (2011). An EAR-motif-containing ERF transcription factor affects herbivore-induced signaling, defense and resistance in rice. Plant J. 68, 583–596. doi: 10.1111/j.1365-313X.2011.04709.x
Lu, J., Li, J. C., Ju, H. P., Liu, X. L., Erb, M., Wang, X., et al. (2014). Contrasting effects of ethylene biosynthesis on induced plant resistance against a chewing and a piercing-sucking herbivore in rice. Mol. Plant 7, 1670–1682. doi: 10.1093/mp/ssu085
Lv, W. T., Du, B., Shangguan, X. X., Zhao, Y., Pan, Y. F., Zhu, L. L., et al. (2014). BAC and RNA sequencing reveal the brown planthopper resistance gene BPH15 in a recombination cold spot that mediates a unique defense mechanism. BMC Genom. 15, 674. doi: 10.1186/1471-2164-15-674
Ma, F. L., Yang, X. F., Shi, Z. Y., Miao, X. X. (2020). Novel crosstalk between ethylene- and jasmonic acid-pathway responses to a piercing-sucking insect in rice. New Phytol. 225, 474–487. doi: 10.1111/nph.16111
Mangrauthia, S. K., Bhogireddy, S., Agarwal, S., Prasanth, V. V., Voleti, S. R., Neelamraju, S., et al. (2017). Genome-wide changes in microRNA expression during short and prolonged heat stress and recovery in contrasting rice cultivars. J. Exp. Bot. 68, 2399–2412. doi: 10.1093/jxb/erx111
Marin, E., Jouannet, V., Herz, A., Lokerse, A. S., Weijers, D., Vaucheret, H., et al. (2010). miR390, Arabidopsis TAS3 tasiRNAs, and their AUXIN RESPONSE FACTOR targets define an autoregulatory network quantitatively regulating lateral root growth. Plant Cell. 22, 1104–1117. doi: 10.1105/tpc.109.072553
Mei, M., Wei, J., Ai, W. F., Zhang, L. J., Lu, X. J. (2021). Integrated RNA and miRNA sequencing analysis reveals a complex regulatory network of Magnolia sieboldii seed germination. Sci. Rep. 11, 10842. doi: 10.1038/s41598-021-90270-y
Nanda, S., Wan, P. J., Yuan, S. Y., Lai, F. X., Wang, W. X., Fu, Q. (2018). Differential responses of OsMPKs in IR56 rice to two BPH populations of different virulence levels. Int. J. Mol. Sci. 19, 4030. doi: 10.3390/ijms19124030
Nanda, S., Yuan, S. Y., Lai, F. X., Wang, W. X., Fu, Q., Wan, P. J. (2020). Identification and analysis of miRNAs in IR56 rice in response to BPH infestations of different virulence levels. Sci. Rep. 10, 19093. doi: 10.1038/s41598-020-76198-9
Natarajan, B., Kalsi, H. S., Godbole, P., Malankar, N., Thiagarayaselvam, A., Siddappa, S., et al. (2018). MiRNA160 is associated with local defense and systemic acquired resistance against Phytophthora infestans infection in potato. J. Exp. Bot. 69, 2023–2036. doi: 10.1093/jxb/ery025
Pachamuthu, K., Hari Sundar, V. (2022). Nitrate-dependent regulation of miR444-OsMADS27 signalling cascade controls root development in rice. J. Exp. Bot. 73, 3511–3530. doi: 10.1093/jxb/erac083
Pan, G., Liu, Y. Q., Ji, L. S., Zhang, X., He, J., Huang, J., et al. (2018). Brassinosteroids mediate susceptibility to brown planthopper by integrating with the salicylic acid and jasmonic acid pathways in rice. J. Exp. Bot. 69, 4433–4442. doi: 10.1093/jxb/ery223
Peng, T., Sun, H. Z., Qiao, M. M., Zhao, Y. F., Du, Y. X., Zhang, J., et al. (2014). Differentially expressed microRNA cohorts in seed development may contribute to poor grain filling of inferior spikelets in rice. BMC Plant Biol. 14, 196. doi: 10.1186/s12870-014-0196-4
Salvador-Guirao, R., Hsing, Y., San Segundo, B. (2018). The polycistronic miR166k-166h positively regulates rice immunity via post-transcriptional control of EIN2. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00337
Schäfer, M., Meza-Canales, I. D., Brütting, C., Baldwin, I. T., Meldau, S. (2015). Cytokinin concentrations and CHASE-DOMAIN CONTAINING HIS KINASE 2 (NaCHK2)- and NaCHK3-mediated perception modulate herbivory-induced defense signaling and defenses in Nicotiana attenuata. New Phytol. 207, 645–658. doi: 10.1111/nph.13404
Schweizer, F., Fernández-Calvo, P., Zander, M., Diez-Diaz, M., Fonseca, S., Glauser, G., et al. (2013). Arabidopsis basic helix-loop-helix transcription factors MYC2, MYC3, and MYC4 regulate glucosinolate biosynthesis, insect performance, and feeding behavior. Plant Cell. 25, 3117–3132. doi: 10.1105/tpc.113.115139
Sharma, N., Tripathi, A., Sanan-Mishra, N. (2015). Profiling the expression domains of a rice-specific microRNA under stress. Front. Plant Sci. 6. doi: 10.3389/fpls.2015.00333
Shi, S. J., Wang, H. Y., Nie, L. Y., Tan, D., Zhou, C., Zhang, Q., et al. (2021). Bph30 confers resistance to brown planthopper by fortifying sclerenchyma in rice leaf sheaths. Mol. Plant 14, 1714–1732. doi: 10.1016/j.molp.2021.07.004
Shi, S. J., Zha, W. J., Yu, X.Y., Wu, Y., Li, S. H., Xu, H. S., et al. (2023). Integrated transcriptomics and metabolomics analysis provide insight into the resistance response of rice against brown planthopper. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1213257
Simon, J. C., d’Alencon, E., Guy, E., Jacquin-Joly, E., Jaquiéry, J., Nouhaud, P., et al. (2015). Genomics of adaptation to host-plants in herbivorous insects. Brief. Funct. Genomics 14, 413–423. doi: 10.1093/bfgp/elv015
Sun, B., Shen, Y. J., Chen, S., Shi, Z. Y., Li, H. C., Miao, X. (2023). A novel transcriptional repressor complex MYB22-TOPLESS-HDAC1 promotes rice resistance to brown planthopper by repressing F3′H expression. New Phytol. 239, 720–738. doi: 10.1111/nph.18958
Tan, J. Y., Wu, Y., Guo, J. P., Li, H. M., Zhu, L. L., Chen, R. Z., et al. (2020). A combined microRNA and transcriptome analyses illuminates the resistance response of rice against brown planthopper. BMC Genom. 21, 144. doi: 10.1186/s12864-020-6556-6
Thaler, J., Humphrey, P., Whiteman, N. (2012). Evolution of jasmonate and salicylate signal crosstalk. Trends Plant Sci. 17, 260–270. doi: 10.1016/j.tplants.2012.02.010
Volodarsky, D., Leviatan, N., Otcheretianski, A., Fluhr, R. (2009). HORMONOMETER: a tool for discerning transcript signatures of hormone action in the Arabidopsis transcriptome. Plant Physiol. 150, 1796–1805. doi: 10.1104/pp.109.138289
Wan, P. J., Zhou, R. N., Nanda, S., He, J. C., Yuan, S. Y., Wang, W. X., et al. (2019). Phenotypic and transcriptomic responses of two Nilaparvata lugens populations to the Mudgo rice containing. Bph1. Sci. Rep. 9, 14049. doi: 10.1038/s41598-019-50632-z
Wang, Y., Cao, L. M., Zhang, Y. X., Cao, C. X., Liu, F., Huang, F. K., et al. (2015). Map-based cloning and characterization of BPH29, a B3 domain-containing recessive gene conferring brown planthopper resistance in rice. J. Exp. Bot. 66, 6035–6045. doi: 10.1093/jxb/erv318
Wang, Y. B., Guo, H. M., Li, H. C., Zhang, H., Miao, X. X. (2012). Identification of transcription factors potential related to brown planthopper resistance in rice via microarray expression profiling. BMC Genom. 13, 687. doi: 10.1186/1471-2164-13-687
Wang, H. H., Hao, J. J., Chen, X. J., Hao, Z. N., Wang, X., Lou, Y. G., et al. (2007). Overexpression of rice WRKY89 enhances ultraviolet B tolerance and disease resistance in rice plants. Plant Mol. Biol. 65, 799–815. doi: 10.1007/s11103-007-9244-x
Wang, Y. Q., Liu, Q. S., Du, L. X., Hallerman, E. M., Li, Y. H. (2020). Transcriptomic and metabolomic responses of rice plants to Cnaphalocrocis medinalis caterpillar infestation. Insects 11, 705–705. doi: 10.3390/insects11100705
Wang, H. Y., Shi, S. J., Hua, W. (2023). Advances of herbivore-secreted elicitors and effectors in plant-insect interactions. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1176048
Wang, Y. Y., Wang, X. L., Yuan, H. Y., Chen, R. Z., Zhu, L. L., He, R. F., et al. (2008). Responses of two contrasting genotypes of rice to brown planthopper. Mol. Plant Microbe Interact. 21, 122–132. doi: 10.1094/MPMI-21-1-0122
Wu, Y., Lv, W. T., Hu, L., Rao, W. W., Zeng, Y., Zhu, L. L., et al. (2017). Identification and analysis of brown planthopper-responsive microRNAs in resistant and susceptible rice plants. Sci. Rep. 7, 8712. doi: 10.1038/s41598-017-09143-y
Xiao, C., Hu, J., Ao, Y. T., Chen, M. X., Gao, G. J., Zhang, Q. L., et al. (2016). Development and evaluation of near-isogenic lines for brown planthopper resistance in rice cv. 9311. Sci. Rep. 6, 38159. doi: 10.1038/srep38159
Xu, J., Wang, X. J., Zu, H. Y., Zeng, X., Baldwin, I. T., Lou, Y. G., et al. (2021). Molecular dissection of rice phytohormone signaling involved in resistance to a piercing-sucking herbivore. New Phytol. 230, 1639–1652. doi: 10.1111/nph.17251
Yang, H. Y., You, A. Q., Yang, Z. F., Zhang, F. T., He, R. F., Zhu, L. L., et al. (2004). High-resolution genetic mapping at the BPH15 locus for brown planthopper resistance in rice (Oryza sativa L.). Theor. Appl. Genet. 110, 182–191. doi: 10.1007/s00122-004-1844-0
Yates, A. D., Michel, A. (2018). Mechanisms of aphid adaptation to host plant resistance. Curr. Opin. Insect Sci. 26, 41–49. doi: 10.1016/j.cois.2018.01.003
Ye, M., Luo, S. M., Xie, J. F., Li, Y. F., Xu, T., Liu, Y., et al. (2012). Silencing COI1 in rice increases susceptibility to chewing insects and impairs inducible defense. PloS One 7, e36214. doi: 10.1371/journal.pone.0036214
Yuan, H. Y., Chen, X. P., Zhu, L. L., He, G. C. (2005). Identification of genes responsive to brown planthopper Nilaparvata lugens Stål (Homoptera: Delphacidae) feeding in rice. Planta 221, 105–112. doi: 10.1007/s00425-004-1422-3
Zhang, X., Bao, Y. L., Shan, D. Q., Wang, Z. H., Song, X. N., Wang, Z. Y., et al. (2018). Magnaporthe oryzae induces the expression of a microRNA to suppress the immune response in rice. Plant Physiol. 177, 352–368. doi: 10.1104/pp.17.01665
Zhang, C., Ding, Z. M., Wu, K. C., Yang, L., Li, Y., Yang, Z., et al. (2016). Suppression of jasmonic acid-mediated defense by viral-inducible microRNA319 facilitates virus infection in rice. Mol. Plant 9, 1302–1314. doi: 10.1016/j.molp.2016.06.014
Zhang, X., Liu, D. M., Gao, D., Zhao, W. N., Du, H. Y., Qiu, Z. Y., et al. (2022). Cytokinin confers brown planthopper resistance by elevating jasmonic acid pathway in rice. Int. J. Mol. Sci. 23, 5946–5946. doi: 10.3390/ijms23115946
Zhang, S. H., Yue, Y., Sheng, L., Wu, Y. Z., Fan, G. H., Li, A., et al. (2013). PASmiR: a literature-curated database for miRNA molecular regulation in plant response to abiotic stress. BMC Plant Biol. 13, 33. doi: 10.1186/1471-2229-13-33
Zhang, F. T., Zhu, L. L., He, G. C. (2004). Differential gene expression in response to brown planthopper feeding in rice. J. Plant Physiol. 161, 53–62. doi: 10.1078/0176-1617-01179
Zhao, Y., Huang, J., Wang, Z. Z., Jing, S. L., Wang, Y., Ouyang, Y. D., et al. (2016). Allelic diversity in an NLR gene BPH9 enables rice to combat planthopper variation. Proc. Natl. Acad. Sci. U. S. A. 113, 12850–12855. doi: 10.1073/pnas.1614862113
Zhao, C., Ma, J. J., Zhang, Y. H., Yang, S. X., Feng, X. C., Yan, J. (2022). The miR166 mediated regulatory module controls plant height by regulating gibberellic acid biosynthesis and catabolism in soybean. J. Integr. Plant Biol. 64, 995–1006. doi: 10.1111/jipb.13253
Zhao, Y. F., Wen, H. L., Teotia, S., Du, Y. X., Zhang, J., Li, J. Z., et al. (2017). Suppression of microRNA159 impacts multiple agronomic traits in rice (Oryza sativa L.). BMC Plant Biol. 17, 215. doi: 10.1186/s12870-017-1171-7
Zheng, Y., He, J. C., Wan, P. J., Lai, F. X., Sun, Y. Q., Lin, J. J., et al. (2016). Virulence characteristics of Nilaparvata lugens (Stål) reared on resistant rice variety IR56. Chin. J. Rice Sci. 30, 552–558. doi: 10.16819/j.1001-7216.2016.6016
Zheng, X. H., Zhu, L. L., He, G. C. (2021). Genetic and molecular understanding of host rice resistance and Nilaparvata lugens adaptation. Curr. Opin. Insect Sci. 45, 14–20. doi: 10.1016/j.cois.2020.11.005
Keywords: rice, Bph15, brown planthopper, virulent populations, resistance mechanism
Citation: Yu B, Geng M, Xue Y, Yu Q, Lu B, Liu M, Shao Y, Li C, Xu J, Li J, Hu W, Tang H, Li P, Liu Q and Jing S (2024) Combined miRNA and mRNA sequencing reveals the defensive strategies of resistant YHY15 rice against differentially virulent brown planthoppers. Front. Plant Sci. 15:1366515. doi: 10.3389/fpls.2024.1366515
Received: 06 January 2024; Accepted: 26 February 2024;
Published: 18 March 2024.
Edited by:
Lida Zhang, Shanghai Jiao Tong University, ChinaReviewed by:
Zhenying Shi, Chinese Academy of Sciences (CAS), ChinaZhaohai Wang, Jiangxi Agricultural University, China
Copyright © 2024 Yu, Geng, Xue, Yu, Lu, Liu, Shao, Li, Xu, Li, Hu, Tang, Li, Liu and Jing. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qingsong Liu, qingsongliu@henu.edu.cn; Shengli Jing, shljing@xynu.edu.cn