Comparative analysis of diet-associated responses in two rice planthopper species
BMC Genomics volume 21, Article number: 565 (2020)
Host adaptation is the primary determinant of insect diversification. However, knowledge of different host ranges in closely related species remains scarce. The brown planthopper (Nilaparvata lugens, BPH) and the small brown planthopper (Laodelphax striatellus, SBPH) are the most destructive insect pests within the family Delphacidae. These two species differ in their host range (SBPH can well colonize rice and wheat plants, whereas BPH survives on only rice plants), but the underlying mechanism of this difference remains unknown. High-throughput sequencing provides a powerful approach for analyzing the association between changes in gene expression and the physiological responses of insects. Therefore, gut transcriptomes were performed to elucidate the genes associated with host adaptation in planthoppers. The comparative analysis of planthopper responses to different diets will improve our knowledge of host adaptation regarding herbivorous insects.
In the present study, we analyzed the change in gene expression of SBPHs that were transferred from rice plants to wheat plants over the short term (rSBPH vs tSBPH) or were colonized on wheat plants over the long term (rSBPH vs wSBPH). The results showed that the majority of differentially expressed genes in SBPH showed similar changes in expression for short-term transfer and long-term colonization. Based on a comparative analysis of BPH and SBPH after transfer, the genes associated with sugar transporters and heat-shock proteins showed similar variation. However, most of the genes were differentially regulated between the two species. The detoxification-related genes were upregulated in SBPH after transfer from the rice plants to the wheat plants, but these genes were downregulated in BPH under the same conditions. In contrast, ribosomal-related genes were downregulated in SBPH after transfer, but these genes were upregulated in BPH under the same conditions.
The results of this study provide evidence that host plants played a dominant role in shaping gene expression and that the low fitness of BPH on wheat plants might be determined within 24 h after transfer. This study deepens our understanding of different host ranges for the two planthopper species, which may provide a potential strategy for pest management.
Most herbivorous insects are restricted to a narrow range of hosts, whereas other insects are adapted to a wider host range . Generally, successful host adaptation involves several essential traits of herbivorous insects, including the abilities to detect the correct plants, acquire nutrients while avoiding intoxication, and overcome plant defenses . Absence of these essential traits have been reported to result in the failure of an insect to feed on specific plants [3,4,5,6]. To date, studies of diet-associated responses have mainly focused on a specific insect fed on different diets . The research findings have shown that the genes associated with detoxification, digestion, and transport are significantly influenced by host transfer . However, knowledge of adaptation to disparate hosts in different insect species is scarce, and most studies have primarily focused on differences in detoxification across species in response to plant secondary metabolites [8,9,10]. Comparative transcriptomic analysis of closely related taxa provides an ideal approach for revealing the distinctions in host adaptation across herbivorous species.
The brown planthopper (Nilaparvata lugens, BPH) and the small brown planthopper (Laodelphax striatellus, SBPH) are closely related insect pests in the family Delphacidae. Although they both use rice plants as their primary food source, their host ranges differ. BPH is a monophagous insect pest restricted to rice plants, whereas SBPH is an oligophagous insect that can feed on rice, wheat, and other gramineous plants . In a wheat-rice rotation system, SBPH, which is able to overwinter in temperate zones, shifts between rice and wheat plants each year . In contrast, the northern border of the overwintering areas of BPH is at approximately 21–25°N. The migratory BPH is seldom in contact with wheat plants in the locations where BPH overwinters . The driving factors that are potentially implicated in the different host ranges of planthoppers have been investigated for decades. Denno & Roderick suggested that a chemical barrier rather than a physical barrier might prevent BPHs from feeding on non-adaptive plants, as the insects could still settle and insert their stylets into the plant tissues . SBPHs perform better on rice than on wheat plants [15, 16]. Liu et al. found that feeding on different plants significantly influenced the activity of detoxification enzymes in SBPH yet had little impact on BPH . Additionally, the presence of endosymbionts, which provide essential amino acids, was also critical for planthopper colonization . The sap compositions of rice and wheat plants differ, particularly in their amino acid concentrations [19, 20]. These differences in nutrient content influence planthopper performance. Recently, the genomic information for both planthoppers and their host plants (Oryza sativa and Triticum aestivum) has become available [18, 21, 22], providing an opportunity to study the mechanism of adaptation to different hosts in these insects.
The gut is the main organ where digestion and detoxification occur, and it exerts an important influence on insect feeding on specific plants. Through omics technologies, changes in gut physiology in response to different host plants have been investigated in several herbivorous species, especially the lepidopterans . The expression patterns of gut-associated transcripts, including digestive and detoxifying enzymes, transporters, and peritrophic genes, have been reported to be substantially altered when larvae are exposed to novel diets [9, 23, 24]. For generalist larvae, there is profound transcriptional variation in the gut to overcome the detrimental effects of plant secondary metabolites [24, 25]. For specialist larvae, a specific detoxification system avoids the activation of general stress responses and minimizes the metabolic costs [9, 25], which might restrain a specialist from adapting to a novel phytotoxin . In our previous work, we determined that CYP4DE1 was critical for SBPH to feed on wheat hosts, with increased mortality observed in dsCYP4DE1-treated SBPHs on wheat plants but not on rice plants . CYP4DE1 expression was upregulated 12 h after transfer to wheat plants. Interestingly, the induction of CYP4DE1 was specifically detected in the gut but not other tissues . This finding reveals the potential role of the gut in planthoppers feeding on different hosts.
Gene regulation upon an initial exposure (short-term transfer) greatly influences the subsequent adaptation of herbivores to a novel environment . To understand the successful colonization of SBPH on rice and wheat plants, we profiled the gut gene expression patterns of SBPHs colonizing rice plants (rSBPH), colonizing wheat plants (wSBPH), and transferred from rice plants to wheat plants (tSBPH). By comparing the differentially expressed genes (DEGs) of SBPH between short-term transfer (rSBPH vs tSBPH) and long-term colonization (rSBPH vs wSBPH), the rapid response of SBPHs to the change in host plant and the leading role of the host plant in shaping gene expression were documented. To elucidate the differences in host adaptation between BPH and SBPH, the gut transcriptomes of BPHs colonizing rice plants (rBPH) and transferred to wheat plants (tBPH) were also profiled. Based on a comparative analysis of short-term transfer in BPHs (rBPH vs tBPH) and SBPHs (rSBPH vs tSBPH), the potential mechanism of host adaptation in the two planthoppers was examined. This study increases our knowledge of diet-associated responses in herbivorous insects, which might suggest a potential strategy for pest management.
Performance of planthoppers on rice and wheat plants
Rice plants, but not wheat plants, can be colonized by BPH with only rice-colonizing BPH strains (rBPH) produced under laboratory conditions. In contrast, SBPH can successfully colonize both rice and wheat plants, and two SBPH strains (rSBPH and wSBPH) were maintained in our laboratory for over 30 generations. According to the survival analysis, more than 90% of rBPH survived on rice plants for 12 days, which was significantly higher than that on wheat plants (Fig. 1a). Additionally, rBPH survived longer on wheat plants (LT50 = 6.1 days) than those provided with water only (LT50 = 3.3 days). These results indicate that BPH could ingest wheat sap and survive on wheat plants for a short time but not for an extended period. For SBPH, both rSBPH and wSBPH successfully survived on rice and wheat plants (Fig. 1b), which is consistent with previous reports .
Overview of RNA sequencing data
To explore the mechanism underlying the different performances of the two planthoppers, rice-colonized planthoppers were transferred to wheat plants for 24 h (short-term transfer) or reared on wheat plants for over 30 generations (long-term colonization). The guts of planthoppers that colonized rice plants (rSBPH and rBPH) or wheat plants (wSBPH) or were transferred from rice plants to wheat plants (tSBPH and tBPH) were then isolated and underwent high-throughput sequencing. A total of 15 libraries (5 treatments and 3 biological replicates for each treatment) were generated, with clean reads exceeding 45 million in each library. The clean reads were mapped to their reference genomes [18, 21]. For BPH, 75–83% of clean reads were mapped to the reference genome. For SBPH, 60–66% of clean reads were mapped to the reference genome. According to the saturation analysis, the number of detected genes decreased as the number of reads increased, and the library capacity reached saturation when the number of sequence reads approached 20.0 million (Fig. S1, Supporting information). Furthermore, principal component analysis (PCA) demonstrated that the expression patterns of wSBPH and tSBPH were closely related, indicating that host plants exerted a nonnegligible influence on gene expression (Fig. S2, Supporting information).
Analysis of differentially expressed genes (DEGs)
Gene expression changes were analyzed by comparing rice colony planthoppers to transfer planthoppers (tSBPH vs rSBPH and tBPH vs rBPH) and rice colony planthoppers to wheat-colony planthoppers (wSBPH vs rSBPH) using a threshold change of > 2-fold and an FDR-adjusted p-value of < 0.05. For rice colony planthoppers transferred to wheat plants, a total of 2877 and 2638 genes were differentially expressed in SBPH and BPH, respectively (Fig. 2). There were 2372 genes upregulated and 505 genes downregulated when rSBPH was transferred to wheat hosts (Table S1). Among these DEGs, genes participating in signal transduction were particularly upregulated. CYP4DE1, which mediates wheat adaptation and ethiprole tolerance  in SBPH, was also significantly induced after transfer (Fig. S3, Supporting information). In contrast, 71 genes related to ribosomal proteins and 48 genes related to oxidative phosphorylation were significantly downregulated, indicating decreased protein production and energy metabolism (Fig. S4, Supporting information). In BPH, the number of genes downregulated (2171 genes) exceeded the number of genes upregulated (467 genes) (Table S2). The majority of genes that were significantly downregulated were associated with intestinal mucins, serine proteinases, and sugar transporters. In addition, reduced expression was also found for detoxification-related genes (Fig. S3, Supporting information), which included 9 ABC transporters, 8 P450s, 5 UGTs, and 1 GST. In contrast to SBPH, the majority of ribosomal proteins were upregulated in BPH (Fig. S4, Supporting information). Cuticular proteins, which form the insect cuticle and are involved in insect molting, were dramatically upregulated after BPH was transferred to wheat (Table S2).
In the comparison of rSBPH and wSBPH, a total of 2516 DEGs were identified (Fig. 2). Strikingly, 90.9% of DEGs (2288 genes) showed higher expression in wSBPH than in rSBPH (Table S3), with genes related to peroxisomal biogenesis factor, nucleotide exchange factor, peptide transporter, and CYP6FK1 exhibiting the most dramatic changes. Similar to the patterns of rSBPH transferred to wheat hosts, 37 genes participating in signal transduction were significantly enriched. Among the 228 downregulated genes, the most dramatic changes in zinc metalloproteinase, UGT, and alpha-glucosidase were observed. Other downregulated genes participating in chitin metabolism, carbohydrate derivative metabolism, starch and sucrose metabolism, and oxidative phosphorylation were significantly enriched (Table S3).
Classification of SBPH genes associated with diet changes
To elucidate the successful colonization of SBPH on rice and wheat plants, the DEGs of SBPH for short-term transfer (rSBPH vs tSBPH) and long-term colonization (rSBPH vs wSBPH) were analyzed. Based on the gene expression changes in response to different diets, the DEGs of SBPH were classified into four types (Fig. 3): I) genes changed in the same direction for short-term transfer and long-term colonization, II) genes changed in opposite direction for short-term transfer and long-term colonization, III) genes changed in response to short-term transfer but not long-term colonization, and IV) genes changed in response to long-term colonization but not short-term transfer.
There were 1558 genes associated with a type I response (Fig. 3; Table S4). The enrichment analysis showed that genes participating in signal transduction and immune system were significantly overrepresented. Only 22 genes were associated with a type II response (Fig. 3; Table S4). Four genes were downregulated after transfer but were dramatically upregulated during colonization, whereas the other 18 genes showed the reciprocal expression pattern. There were 1297 genes associated with a type III response (Fig. 3; Table S4). The enrichment analysis showed that the ribosome pathway, oxidative phosphorylation pathway, and retrograde endocannabinoid signaling pathway were significantly overrepresented. It was noteworthy that the expression level of ribosome proteins was initially suppressed but recovered when SBPH colonized wheat plants over the long term. A total of 936 genes were associated with a type IV response (Fig. 3; Table S4). However, we failed to find GO terms or KEGG terms that were significantly enriched. Genes such as integrin alpha-PS4-like, integumentary mucin, and proliferation-associated protein showed a higher expression level in wSBPH.
Comparative genomics in response to host transfer
To comprehend the different diet-associated responses across species, 6139 gene families with only one ortholog in BPH and SBPH were selected and compared. A total of 1995 gene families were differentially expressed in at least one planthopper species after transfer (Table S5), among which 370 genes were responsive to host transfer in both planthoppers. Interestingly, only 22 genes changed (14 genes upregulated and 8 genes downregulated) in the same direction in both species, including heat-shock protein, prophenoloxidase activating factor, MAP kinase-interacting serine/threonine-protein kinase, and small nuclear ribonucleoprotein. Nonetheless, other 348 genes showed different expression patterns between BPH and SBPH. Among these, 293 genes were upregulated in SBPH after transfer, but downregulated in BPH; 55 genes, including 25 ribosomal proteins, were downregulated in SBPH after transfer but were upregulated in BPH.
To confirm the validity of the transcriptomic data, 15 SBPH genes and 18 BPH genes were selected for qPCR analysis. Thirteen SBPH genes (Fig. 4) and 17 BPH genes (Fig. 5) showed a concordant direction of change for the qPCR and transcriptomic results, indicating acceptable accuracy of the DEG transcriptomic results. The heat-shock proteins were significantly upregulated after BPH and SBPH transfer to wheat. The expression levels of the ABC transporters and cytochrome P450 were significantly increased after SBPH transfer to wheat but were significantly decreased in BPH under the same conditions. The ribosomal proteins were significantly downregulated in SBPH after transfer to wheat but were significantly upregulated in BPH after transfer. It is worth noting that two trehalose transporters in BPH showed significantly different changes in expression. The trehalose transporter NLU013658.1 was dramatically downregulated after BPH was transferred to wheat, but NLU003716.1 was dramatically upregulated under the same conditions. In SBPH, the nucleotide exchange factor, peroxisomal biogenesis factor, and peptide transporter were significantly upregulated after transfer, but the venom serine carboxypeptidase-like and maltase were significantly downregulated. In BPH, genes such as cryptosporidial mucin, serine proteinase stubble, and peptide methionine sulfoxide reductase were significantly downregulated after transfer, but the cuticle protein 16.5-like, chemosensory protein, and lipid storage droplets surface-binding protein were significantly upregulated.
The wSBPH and rSBPH used in this study were reared under laboratory conditions for more than 30 generations, whereas the tSBPH was collected by transferring rSBPH to wheat plants within 24 h. Surprisingly, the expression pattern of tSBPH was more similar to wSBPH than rSBPH (Fig. S2, Supporting information), indicating a rapid response by SBPH to the change in host plant and a leading role of the host plant in shaping gene expression. Additionally, the transfer of rSBPH to wheat plants gave rise to broad transcriptional readjustments, and the majority of DEGs were changed in the same direction in response to short-term transfer and long-term colonization. SBPH shifts between rice and wheat plants each year, whereas BPH seldom comes in contact with wheat plants [12, 13]. The frequent host transfer might help SBPH rapidly respond to new hosts, resulting in successful adaptation to wheat. A similar phenomenon was reported in Tetranychus urticae, which rapidly adapted to pre-exposed hosts . In contrast, different expression patterns were revealed in 94.1% of homologous genes in BPH and SBPH when planthoppers were transferred to wheat plants (Table S5). Specifically, homologous genes associated with ribosomal proteins and detoxification were differentially regulated in the two species. With the majority of genes changed in the “wrong” direction, the low fitness of BPH on wheat plants might be determined within 24 h after transfer.
The genes involved in detoxification were significantly influenced when planthoppers were transferred to wheat plants (Fig. S3, Supporting information). In previous work, CYP4DE1 was found to be responsible for wheat adaptation, and a knockdown of CYP4DE1 significantly lowered the intestinal cell viability of SBPH reared on wheat plants . In this study, we found that CYP4DE1 and other P450 genes were significantly upregulated after rSBPH was transferred to wheat plants, indicating that the induction of P450s might be critical. Similarly, 9 ABC transporters, which transport a large diversity of substrates across lipid membranes and out of cells, were also significantly upregulated in response to wheat plants. The involvement of ABC transporters in xenobiotic transport and insecticide resistance have been well documented . These proteins have also been reported to participate in host adaptation in Manduca sexta , Aphis nerii , and Chrysomela populi . The increased expression of ABC transporters signified that SBPH might accelerate the secretion of metabolites when feeding on wheat plants. Comparative analysis showed that the majority of differentially expressed detoxifying genes were upregulated when SBPH was transferred to wheat plants. However, these genes showed reduced expression in BPH (Fig. S3, Supporting information). We presumed that this might be associated with a resource-based metabolic trade-off . When coping with such an adverse environment, it might become more critical for BPH to increase basic metabolism than to invest in the detoxification process. Further studies are needed to determine the variation in detoxification between BPH and SBPH on wheat plants.
Large numbers of ribosomal proteins were differentially expressed when the planthoppers were transferred to wheat plants. However, the majority of these genes showed different patterns of change between SBPH and BPH (Fig. S4, Supporting information). Most genes associated with standard roles in protein translation and ribosomal proteins are generally deemed to be stably expressed and are used as housekeeping genes . However, recent studies have demonstrated that ribosomal proteins were differentially expressed in response to host transfers in Bemisia tabaci , Helicoverpa armigera , Polygonia c-album , and Cryptolaemus montrouzieri . In this study, the expression level of ribosomal proteins was significantly downregulated when rSBPH was transferred to wheat but upregulated in rBPH under the same conditions. In an insect-plant model, regulation of ribosomal proteins is believed to counteract the ribosome inactivating proteins (RIPs), which are produced by host plants and inhibit protein synthesis in insects . Different expression patterns of ribosomal proteins in the two planthoppers might influence the ability of the insect to counteract plant RIPs; this deserves further investigation. Additionally, ribosomal proteins play fundamental roles in the cellular process of translation . The large-scale change of ribosomal proteins inevitably influences the regulation of other genes, which might partially explain the distinct variation in gene expression between the two planthopper species.
Genes associated with the amide/peptide biosynthetic process, amino sugar metabolic process, and aminoglycan metabolic process were significantly enriched in both SBPH and BPH after host transfer. Planthoppers are piercing-sucking herbivores that depend on phloem sap containing high concentrations of sugars, amino acids and inorganic ions as their food sources [15, 37]. Previous work has demonstrated that the sucrose concentration of rice phloem sap and wheat phloem sap is very similar, but the latter has a higher concentration of total amino acids compared with the former [19, 20]. In addition, the predominance of amides between the two phloem saps was also very different [19, 20]. In previous work, we found that the BPH genome lacks the ability to carry out de novo synthesis of some essential amino acids, requiring an additional supply from endosymbionts or plant phloem sap . When planthoppers were sustained on different hosts, their endosymbiotic communities were changed . Our presumption is that the difference in phloem sap composition or endosymbionts might influence nutritional status, which might result in a change in basic metabolism (particularly amino acid metabolic processes). This hypothesis has been partially verified by planthoppers reared on resistant rice, in which the physiological adaptation of insects to a novel rice variety was closely linked to a change in amino acid metabolism .
Sugar transporters were significantly differentially expressed when planthoppers were transferred to or colonized on wheat plants (Fig. S5, Supporting information). Sugar transporters mediate the movement of sugars into and out of cells in a diverse array of organisms and are vital for the utilization of ingested sugars as a nutritional resource and for the maintenance of osmotic balance in insects . We found that the majority of differentially expressed sugar transporters belonged to facilitated trehalose transporters. Trehalose is the major hemolymph sugar in most insects . Its concentration in hemolymph is determined by the balance between discharge from the fat body and uptake by other tissues via facilitated trehalose transporters . In Anopheles gambiae, suppression of the trehalose transporter significantly diminished the hemolymph trehalose concentration, making the insect more sensitive to stress . Similarly, the trehalose transporter from Polypedilum vanderplanki enhanced the desiccation tolerance of cell lines . Involvement of the trehalose transporter in host adaptation has only been reported in a few insects . However, increasing evidence demonstrates that trehalose, the regulatory product of the trehalose transporter, plays a key role in insects transitioning to different diets. For example, the trehalose concentration in hemolymph affected the food preference and carbohydrate intake in Heliothis zea and Manduca sexta [47, 48], and trehalose homeostasis was critical for Drosophila adapting to different dietary conditions . In this study, the majority of sugar transporters were suppressed in SBPH and BPH reared on wheat plants. Accordingly, we inferred that the two planthoppers might use a similar strategy to maintain sugar homeostasis. However, it was also worth noting that two of the trehalose transporters in BPH were dramatically upregulated. BPH harbors many sugar transporter genes, and they differ in their function . Further studies are needed to illustrate the function of sugar transporters in planthoppers coping with different diets.
Heat-shock proteins were prominently upregulated when the two planthopper species were transferred to or colonized wheat plants (Fig. S5, Supporting information). Heat-shock proteins are well-known stress proteins that respond to an array of stresses including thermal hardening, oxidative hardening, chemical pesticides, and desiccation [51, 52]. The induced expression of these proteins as a significant molecular chaperone prevents the irreversible denaturation of proteins and enhances the coping capacity of insects in the face of stress . Previously, the influence of diet quality on heat-shock proteins expression was described in Drosophila melanogaster, and diet-induced heat-shock proteins increased the heat and desiccation tolerance of flies . Feeding on wheat was a biotic stressor for BPH and SBPH , which might result in increased heat-shock protein expression.
Overall, SBPH successfully coped with wheat hosts, and the majority of DEGs responded similarly to short-term transfer and long-term colonization. Compared with SBPH, BPH showed distinct gene expression changes after transfer to the wheat plants. The different pattern of changes in gene expression between BPH and SBPH might account for their adaptive differences on wheat plants. Specifically, genes associated with sugar transporters and heat-shock proteins showed similar expression trends for BPH and SBPH. Other genes associated with detoxification, ribosomal proteins, and amino acid metabolism were regulated differently between the two planthoppers. Our work increases our knowledge of planthoppers adapting to rice and wheat hosts, and might be useful in pest management.
The SBPH and BPH populations used in this study were originally collected from a rice field at Huajiachi Campus, Zhejiang University, Hangzhou, China. This rice field (30.271 ° N, 120.199 ° E) was specifically used for field experiments. The rSBPH and wSBPH were generated by rearing SBPHs on fresh rice (Xiushui 134) and wheat (Luyuan 502) plants for over 30 generations. The rBPH have been reared on rice (Xiushui 134) plants for more than 40 generations. All insects were reared in large cages, and more than 300 individuals were used for propagation in each generation.
Survival analysis of rSBPH on different hosts has been done previously . However, the performance of wSBPH, as well as rBPH, on wheat plants remained unknown. For BPH, 3rd instar nymphs (rBPH) were reared on rice plants (4–5-leaf stage) or were transferred to wheat plants (4–5-leaf stage). The rBPH that provided with water only were used as a control. For SBPH, rSBPH were reared on rice plants (4–5-leaf stage) or were transferred to wheat plants (4–5-leaf stage), and wSBPH were reared on wheat plants (4–5-leaf stage). We recorded the survival rates of each treatment every day. Three biological replicates were conducted for each treatment. Each replicate contains a group of 40–50 insects (approximately 4–5 insects per plant). The log-rank test in SPSS 19.0 (SPSS, Chicago, IL, USA) was used to determine the statistical significance of survival distributions for each treatment group.
For rSBPH, wSBPH, and rBPH populations, newly emerged 5th instar nymphs were collected and continually reared on their original hosts for 24 h. For tSBPH and tBPH populations, newly emerged 5th instar nymphs (maintained on rice plants) were collected and transferred to wheat plants for 24 h. The nymphs were anesthetized on ice and their guts were carefully dissected under a stereomicroscope (COIC, Chongqing, China). Total RNA extraction was performed using TRIzol Total RNA Isolation Kit (Takara, Dalian, China) according to the manufacturer’s instructions. Three replicates (all replicates were collected from the same insect strain) were performed for each treatment, and each replicate contained a pool of 100 guts.
The RNA samples were sent to Beijing Genomics Institute (BGI, Beijing, China) for transcriptomic sequencing as we previously described . A BGISEQ-500 platform (BGI, Beijing, China) was used for library sequencing. After filtering the low quality reads using the internal software, the clean reads from each cDNA library were aligned to reference genome sequences of Nilaparvata lugens (BioProjects: PRJNA398259; Link: https://www.ncbi.nlm.nih.gov/genome/2941?genome_assembly_id=209423)  and Laodelphax striatellus (BioProjects: PRJNA393384; Link: https://www.ncbi.nlm.nih.gov/genome/12198?genome_assembly_id=449717)  using Hierarchical Indexing for Spliced Alignment of Transcripts (HISAT). The clean data of rSBPH, rBPH, wSBPH, tSBPH, and tBPH have been submitted to the database of the NCBI Sequence Read Archive under the accession number PRJNA564687.
The DEGseq software [55, 56] was used to analyze the DEGs, and genes with log2-ratio > 1 and adjusted p value < 0.05 were identified. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and Gene Ontology (GO) were based on the KEGG pathway database (http://www.genome.jp/kegg/) and GO Database (http://www.geneontology.org/).
Gene family identification
Previous work demonstrated that genes associated with detoxification, digestion, and transport played critical roles in host transfer . To analyze these genes in BPH and SBPH, the target genes were firstly identified according to their gene annotations. Then, the identified genes were used as queries to BLAST against the genomic databases of Nilaparvata lugens  and Laodelphax striatellus  using local BLAST (ftp://ftp.ncbi.nlm.nih.gov/blast). Finally, we used the non-redundant (nr) protein database in National Center for Biotechnology Information (NCBI) to validate the gene annotation.
Comparative genomic analysis
The protein families of BPH and SBPH were analyzed using the ORTHOMCL software 2.0.3 (http://orthomcl.org/common/downloads/software/v2.0/) according to the method we previously described . Briefly, the deduced amino acid sequences of BPH and SBPH were firstly downloaded from NCBI (BioProjects: PRJNA398259 and PRJNA393384). Then, all-vs-all BLASTp algorithm was used to find the homologous pairs of sequences with an E-value of <1e-5. The BLASTp result was subsequently converted into a normalized similarity matrix. The matrix was then analyzed by a Markov chain clustering (MCL) for clustering . We use an inflation factor of 1.5 to control the cluster tightness.
qPCR analysis was used to confirm the validity of transcriptomic data. Sample collection and RNA extraction was performed as described above. The concentration of each RNA sample was determined by NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). Reverse transcription was performed by HiScript II Q RT SuperMix (Vazyme, Nanjing, China), and 1 μg RNA was used in a 20-μl reaction system. RNA with no-reverse-transcriptase was used as a negative control. The primers used for qPCR were designed by Primer Premier 6.0 and listed in Table S6. We use two housekeeping genes (β-actin and GAPDH) as internal controls. qPCR was performed using the ABI 7500 Real-Time PCR System (Applied Biosystems, Carlsbad, CA) and the SYBR Green Supermix Kit (Yeasen, Shanghai, China). The program was run under the following conditions: denaturation at 95 °C for 5 min, followed by 40 cycles at 95 °C for 10 s and 60 °C for 30 s. qPCR result was calculated according to a relative quantitative method (2-∆∆Ct) . Three independent biological replicates were performed.
Availability of data and materials
All sequencing data generated in this study were submitted to the NCBI Sequence Read Archive under accession number PRJNA564687. Other related data are available within the manuscript and its additional files. The genome data of Nilaparvata lugens  and Laodelphax striatellus  were downloaded from NCBI under the BioProject: PRJNA398259 and PRJNA393384.
Differentially expressed gene
Markov chain clustering
Principal component analysis
Ribosome inactivating protein
Brown planthopper that colony on rice
Small brown planthopper that colony on rice
Small brown planthopper
Brown planthopper that transferred from rice to wheat
Small brown planthopper that transferred from rice to wheat
Small brown planthopper that colonized on wheat
Schoonhoven LM, Van Loon B, van Loon JJ, Dicke M. Insect-plant biology: Oxford University press on demand; 2005.
Powell G, Tosh CR, Hardie J. Host plant selection by aphids: behavioral, evolutionary, and applied perspectives. Annu Rev Entomol. 2006;51:309–30.
Michael L, Sophie M, Clark AG, Géraldine G, Catherine B, Matzkin LM, Emilie G, Takuji YY, Hiroshi K, Ryusuke N. Mutations in the neverland gene turned Drosophila pachea into an obligate specialist species. Science. 2012;337(6102):1658–61.
Ying Z, Aardema ML, Medina EM, Molly S, Peter A. Parallel molecular evolution in an herbivore community. Science. 2012;337(6102):1634–7.
Chris B, Zimmer CT, Riveron JM, Wilding CS, Wondji CS, Martin K, Field LM, Williamson MS, Ralf N. Gene amplification and microsatellite polymorphism underlie a recent insect host shift. Proc Natl Acad Sci U S A. 2013;110(48):19460–5.
Lavista-Llanos S, Svatoš A, Kai M, Riemensperger T, Birman S, Stensmyr MC, Hansson BS. Dopamine drives Drosophila sechellia adaptation to its toxic host. eLife. 2014;3:e03785.
Seo B, Kwon Y, Jung J, Kim G. Electrical penetration graphic waveforms in relation to the actual positions of the stylet tips of Nilaparvata lugens in rice tissue. J Asia Pac Entomol. 2009;12(2):89–95.
Schweizer F, Vogel H, Reymond P. Arabidopsis glucosinolates trigger a contrasting transcriptomic response in a generalist and a specialist herbivore. Insect Biochem Mol Biol. 2017;85:21–31.
Roy AIW, Vogel H, Chattington S, Larsson MC, Anderson P, Heckel DG, Schlyter F. Diet dependent metabolic responses in three generalist insect herbivores Spodoptera spp. Insect Biochem Mol Biol. 2016;71:91–105.
Ragland GJ, Almskaar K, Vertacnik KL, Gough HM, Feder JL, Hahn DA, Schwarz D. Differences in performance and transcriptome-wide gene expression associated with Rhagoletis (Diptera: Tephritidae) larvae feeding in alternate host fruit environments. Mol Ecol. 2015;24(11):2759–76.
Fu Q, Huang SW. Diagnosis and control of rice pests and diseases (in Chinese). Beijing: Beijing Jindun Press; 2005. p. 93–7.
Otuka A, Matsumura M, Sanada-Morimura S, Takeuchi H, Watanabe T, Ohtsu R, Inoue H. The 2008 overseas mass migration of the small brown planthopper, Laodelphax striatellus, and subsequent outbreak of rice stripe disease in western Japan. Appl Entomol Zool. 2010;45(2):259–66.
Cheng SN, Chen JC, Xi X, Yang LM, Zhu ZL. Studies on the migrations of brown planthopper Nilaparvata lugens Stål. Acta Entomol Sin. 1979;22:1–21.
Denno RF, Roderick GK. Population biology of planthoppers. Annu Rev Entomol. 1990;35(1):489–520.
Huang HJ, Cui JR, Guo Y, Sun JT, Hong XY. Roles of LsCYP4DE1 in wheat adaptation and ethiprole tolerance in Laodelphax striatellus. Insect Biochem Mol Biol. 2018;101:14–23.
Li W, Guo HF, Wang RF, Liu BS, Zhong WF, Fang JC. Comparison of the life tables of Laodelphax striatellus (Homoptera:Delphacidae) on different host plants. Acta Entomol Sin. 2009;52:531–6.
Liu YK, Wang WX, Qiang FU, Lai FX, Luo J. Effects of host plants on activities of detoxification and protective enzymes in three rice planthoppers. Chin J Rice Sci. 2011;25(6):659–66.
Xue J, Zhou X, Zhang CX, Yu LL, Fan HW, Wang Z, Xu HJ, Xi Y, Zhu ZR, Zhou WW. Genomes of the rice pest brown planthopper and its endosymbionts reveal complex complementary contributions for host adaptation. Genome Biol. 2014;15(12):521.
Hayashi H, Chino M. Collection of pure phloem sap from wheat and its chemical composition. Plant Cell Physiol. 1986;27(7):1387–93.
Hayashi H, Chino M. Chemical composition of phloem sap from the uppermost internode of the rice plant. Plant Cell Physiol. 1990;31(2):247–51.
Zhu J, Jiang F, Wang X, Yang P, Bao Y, Zhao W, Wang W, Lu H, Wang Q, Cui N. Genome sequence of the small brown planthopper Laodelphax striatellus. Gigascience. 2017;6(12):1–12.
Huang HJ, Zhang CX, Hong XY. How does saliva function in planthopper–host interactions? Arch Insect Biochem Physiol. 2019;100:e21537.
Jing X, Vogel H, Grebenok RJ, Zhu-Salzman K, Behmer ST. Dietary sterols/steroids and the generalist caterpillar Helicoverpa zea: physiology, biochemistry and midgut gene expression. Insect Biochem Mol Biol. 2012;42(11):835–45.
Herde M, Howe GA. Host plant-specific remodeling of midgut physiology in the generalist insect herbivore Trichoplusia ni. Insect Biochem Mol Biol. 2014;50(1):58–67.
Geetha G, Omprakash M, Thasso G, Silke A, Sebastian BC, Ian Thomas B. Unbiased transcriptional comparisons of generalist and specialist herbivores feeding on progressively defenseless Nicotiana attenuata plants. PLoS One. 2010;5(1):e8735.
Vogel H, Musser RO. Transcriptome responses in herbivorous insects towards host plant and toxin feeding. Annu Plant Rev. 2014;47:197–233.
Wybouw N, Zhurov V, Martel C, Bruinsma KA, Hendrickx F, Grbić V, Van Leeuwen T. Adaptation of a polyphagous herbivore to a novel host plant extensively shapes the transcriptome of herbivore and host. Mol Ecol. 2015;24(18):4647–63.
Qiao H, Luo J, Lai FX, Fu Q, Wang HD, Dai DJ. Fitness of small brown planthopper (Laodelpax striatellus) on different plants. Chin J Rice Sci. 2009;23(1):71–8.
Dermauw W, Leeuwen TV. The ABC gene family in arthropods: comparative genomics and role ininsecticide transport and resistance. Insect Biochem Mol Biol. 2014;45(1):89–110.
Koenig C, Bretschneider A, Heckel DG, Grosse-Wilde E, Hansson BS, Vogel H. The plastic response of Manduca sexta to host and non-host plants. Insect Biochem Mol Biol. 2015;63:72–85.
Ssl B, Rinker DC, Gerardo NM, Abbot P. Transcriptional profile and differential fitness in a specialist milkweed insect across host plants varying in toxicity. Mol Ecol. 2017;26(1):6742–61.
Strauss AS, Wang D, Stock M, Gretscher RR, Groth M, Boland W, Burse A. Tissue-specific transcript profiling for ABC transporters in the sequestering larvae of the phytophagous leaf beetle Chrysomela populi. PLoS One. 2014;9(6):e98637.
Després L, Stalinski R, Faucon F, Navratil V, Viari A, Paris M, Tetreau G, Poupardin R, Riaz MA, Bonin A. Chemical and biological insecticides select distinct gene expression patterns in Aedes aegypti mosquito. Biol Lett. 2014;10(12):20140716.
Ramakrishnan V. Ribosome structure and the mechanism of translation. Cell. 2002;108(4):557–72.
Alon M, Elbaz M, Ben-Zvi MM, Feldmesser E, Vainstein A, Morin S. Insights into the transcriptomics of polyphagy: Bemisia tabaci adaptability to phenylpropanoids involves coordinated expression of defense and metabolic genes. Insect Biochem Mol Biol. 2012;42(4):251–63.
Celorio-Mancera MP, Heckel DG, Vogel H. Transcriptional analysis of physiological pathways in a generalist herbivore: responses to different host plants and plant structures by the cotton bollworm, Helicoverpa armigera. Entomol Exp Appl. 2012;144(1):123–33.
Maria CM, Wheat CW, Heiko V, Lina SD, Niklas J, Sören N. Mechanisms of macroevolution: polyphagous plasticity in butterfly larvae revealed by RNA-Seq. Mol Ecol. 2013;22(19):4884–95.
Li HS, Pan C, Clercq PD, Ślipiński A, Pang H. Variation in life history traits and transcriptome associated with adaptation to diet shifts in the ladybird Cryptolaemus montrouzieri. BMC Genomics. 2016;17(1):281.
Puri M, Kaur I, Perugini MA, Gupta RC. Ribosome-inactivating proteins: current status and biomedical applications. Drug Discov Today. 2012;17(13–14):774–83.
Xu HX, Zheng XS, Yang YJ, Tian JG, Fu Q, Ye GY, Lu ZX. Changes in endosymbiotic bacteria of brown planthoppers during the process of adaptation to different resistant rice varieties. Environ Entomol. 2015;44(3):582–7.
Chen YH, Bernal CC, Tan J, Horgan FG, Fitzgerald MA. Planthopper “adaptation” to resistant rice varieties: changes in amino acid composition over time. J Insect Physiol. 2011;57(10):1375–84.
Price DRG, Tibbles K, Shigenobu S, Smertenko A, Russell CW, Douglas AE, Fitches E, Gatehouse AMR, Gatehouse JA. Sugar transporters of the major facilitator superfamily in aphids; from gene prediction to functional characterization. (special issue: the aphid genome). Insect Mol Biol. 2010;19(s2):97–112.
Thompson S. Trehalose-the insect blood sugar. Adv Insect Physiol. 2003;31:203–85.
Takahiro K, Ayako S, Yasushi K, Yuichi N, Ken-Ichi I, Daisuke T, Masahiko W, Takashi O. Trehalose transporter 1, a facilitated and high-capacity trehalose transporter, allows exogenous trehalose uptake into cells. Proc Natl Acad Sci U S A. 2007;104(28):11585–90.
Kun L, Yuemei D, Yuzheng H, Rasgon JL, Peter A. Impact of trehalose transporter knockdown on Anopheles gambiae stress adaptation and susceptibility to Plasmodium falciparum infection. Proc Natl Acad Sci U S A. 2013;110(43):17504–9.
Nilay C, Menze MA, Heidi E, Halong V, Yarmush ML, Hand SC, Mehmet T. Trehalose transporter from African chironomid larvae improves desiccation tolerance of Chinese hamster ovary cells. Cryobiology. 2012;64(2):91–6.
Friedman S, Waldbauer GP, Eertmoed JE, Naeem M, Ghent AW. Blood trehalose levels have a role in the control of dietary self-selection by Heliothis zea larvae. J Insect Physiol. 1991;37(12):919–28.
Thompson SN, Redak RA, Wang LW. Altered dietary nutrient intake maintain metabolic homeostasis in parasitized larvae of the insect Manduca sexta L. J Exp Biol. 2001;204(23):4065–80.
Yasugi T, Yamada T, Nishimura T. Adaptation to dietary conditions by trehalose metabolism in Drosophila. Sci Rep. 2017;7(1):1619.
Kikuta S, Kikawada T, Hagiwara-Komoda Y, Nakashima N, Noda H. Sugar transporter genes of the brown planthopper, Nilaparvata lugens: a facilitated glucose/fructose transporter. Insect Biochem Mol Biol. 2010;40(11):805–13.
King AM, MacRae TH. Insect heat shock proteins during stress and diapause. Annu Rev Entomol. 2015;60(1):59–75.
Zhao L, Jones W. Expression of heat shock protein genes in insect stress response. Invertebr Surviv J. 2012;9:93–101.
Andersen LH, Kristensen TN, Loeschcke V, Toft SR, Mayntz D. Protein and carbohydrate composition of larval food affects tolerance to thermal stress and desiccation in adult Drosophila melanogaster. J Insect Physiol. 2010;56(4):336–40.
Kang X, Lu YJ, Yang K, Huo SM, Hong XY. Co-infection of Wolbachia and Spiroplasma in spider mite Tetranychus truncatus increases male fitness. Insect Sci. 2019. https://doi.org/10.1111/1744-7917.12696.
de la Paz C-MM, Ahn S-J, Vogel H, Heckel DG. Transcriptional responses underlying the hormetic and detrimental effects of the plant secondary metabolite gossypol on the generalist herbivore Helicoverpa armigera. BMC Genomics. 2011;12(1):575.
Likun W, Zhixing F, Xi W, Xiaowo W, Xuegong Z. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.
Huang HJ, Xue J, Zhuo JC, Cheng RL, Xu HJ, Zhang CX. Comparative analysis of the transcriptional responses to low and high temperatures in three rice planthopper species. Mol Ecol. 2017;26:2726–37.
Li L, Stoeckert CJ Jr, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13(9):2178–89.
TS KJL. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta CT) method. Methods. 2001;25:402–8.
The authors thank Jie Chen for assistance with rice cultivation and management.
This work has received funding from the National Natural Science Foundation of China under grant agreement 31801734 (H.J. Huang) “Study on the mechanism of saliva-mediated interaction between small brown planthopper and host plants” and the Natural Science Foundation of Jiangsu Province under grant agreement BK20180520 (H.J. Huang) “Screen the critical genes that facilitate small brown planthopper adapted to different hosts”. The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
The SBPH and BPH populations were originally collected from rice field at an experimental station of Zhejiang University. No permissions were necessary to collect the specimens in this field. The authors declare that the collections of specimens comply with institutional, national, or international guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Saturation analysis of 15 sequencing libraries.
Principal component analysis (PCA) of gene expression patterns in SBPH (A) and BPH (B). The first two principal components (PC1 and PC2) based on transcriptomic results are shown with each plot representing one sample.
The expression pattern of detoxification-related genes. The differentially expressed genes associated with ABC transporters, cytochrome P450s, UDP-glucuronosyltransferases (UGTs), and esterases in SBPH (A) and BPH (B) are illustrated in the heat map.
The expression pattern of ribosomal proteins. The differentially expressed genes associated with ribosomal proteins in SBPH (A) and BPH (B) are illustrated in the heat map.
The expression pattern of sugar transporters and heat-shock proteins. The differentially expressed genes associated with trehalose transporters, sugar transporters, and heat-shock proteins in SBPH (A, C) and BPH (B, D) are illustrated in the heat map.
Differentially expressed genes when rice-colonized SBPH (rSBPH) were transferred to wheat plants.
Differentially expressed genes when rice-colonized BPH (rBPH) were transferred to wheat plants.
Differentially expressed genes in rice-colonized SBPH (rSBPH) and wheat-colonized SBPH (wSBPH).
Classification of the identified differentially expressed genes in SBPH.
Comparative genomic analysis of genes in response to host transfer in two planthopper species.
Primers used for qPCR.
About this article
Cite this article
Huang, HJ., Cui, JR. & Hong, XY. Comparative analysis of diet-associated responses in two rice planthopper species. BMC Genomics 21, 565 (2020). https://doi.org/10.1186/s12864-020-06976-2
- Nilaparvata lugens
- Laodelphax striatellus
- Comparative transcriptome
- Host adaptation
- Diet-associated responses