Introduction

Chemical pesticides remain the primary means of controlling many of the world’s most damaging arthropod crop pests (Janssen and van Rijn, 2021; Wu et al., 2018a). However, pesticide applications can result in pest resurgence, increases in pest insect populations that exceed natural, untreated population sizes, following an initial reduction of the pest population (Gandara et al., 2024; Guedes et al., 2016; Hardin et al., 1995; Wu et al., 2020). Two mechanisms have been implicated in pest resurgence - the loss of beneficial insects including natural enemies, and pesticide-enhanced reproduction in the pest insect (Wu et al., 2020). In the case of the latter, several pesticides, such as the insecticides triazophos, deltamethrin and the fungicide jinggangmycin, have been reported to stimulate pest reproduction (Ge et al., 2010; Wu et al., 2018b; Zhang et al., 2014). Pesticide-enhanced pest reproduction has been linked to changes in physiology and biochemistry of pest organisms after exposure to pesticides (Guedes et al., 2016; Wu et al., 2020). However, the molecular mechanisms underlying enhanced reproduction associated with pest resurgence remain poorly resolved.

The brown planthopper (BPH), Nilaparvata lugens (Stål), is a notorious pest of rice crops throughout Asia causing annual losses of ∼300 million dollars across major rice producing countries (Wu et al., 2020; Wu et al., 2018a). BPHs inhibit the growth of rice plants by feeding, and also transmits highly damaging plant viruses including rice grassy stunt virus and rice ragged stunt virus (Sōgawa, 1982). Currently, chemical insecticides play an indispensable role in the control of BPH due to their efficiency, rapid effect, and low cost. However, due to the widespread and intensive use of chemical insecticides, BPH has developed resistance to the majority of compounds used for control (Wu et al., 2018a; Zeng et al., 2023).

Emamectin benzoate (EB) and abamectin are avermectin pesticides, and act as allosteric modulators of insect glutamate gated chloride channels (GluCls), inhibiting muscle contractions that lead to the cessation of insect feeding and subsequent death (Ishaaya et al., 2002). These insecticides exhibit particularly strong activity against Lepidoptera such as the rice leaffolder, Cnaphalocrocis medinalis Guénee, an important foliage-feeding insect which attacks rice during the vegetative stage (Chintalapati et al., 2016). Both BPH and the rice leaffolder are migratory pests with overlapping migratory paths, however, their occurrence period in the field differs by approximately one month, with leaffolders appearing earlier than BPH. Therefore, the use of EB to control rice leaffolder has the potential to impact BPHs arriving later, via exposure to sublethal concentrations of this compound. In this regard, we have observed that when farmers use EB and abamectin to control leaffolders on rice crops in China, BPH outbreaks frequently occur in the same field. While sublethal doses of certain pesticides have been shown to enhance fecundity in BPHs, including the insecticides triazophos and deltamethrin (Ge et al., 2013; Ge et al., 2010; Zhang et al., 2014; Zhang et al., 2022b) and the fungicides carbendazim and jinggangmycin (Wu et al., 2018b), whether avermectins trigger resurgence in BPH via insecticide-enhanced reproduction remains unclear.

Reproduction in insects is influenced by external factors such as light (Wang et al., 2021), temperature (Meiselman et al., 2022), humidity (Roy et al., 2015) and nutrition (Smykal and Raikhel, 2015), and endogenous factors such as the juvenile hormone (JH) (Santos et al., 2019), ecdysone (Hun et al., 2022), insulin (Ling and Raikhel, 2018) and target of rapamycin (TOR) pathways (Ahmed et al., 2020; Du et al., 2022; Lu et al., 2016). Of these, JH has been particularly implicated in insecticide-induced enhanced fecundity, with triazophos and deltamethrin treatments leading to increased circulating JH III titers in BPH females (Wu et al., 2020). JH is synthesized and secreted by the corpora allata in insects (Luo et al., 2021), and can promote reproduction by regulating the synthesis and secretion of vitellogenin (Vg) in the female fat body, and stimulating the absorption of Vg by the developing oocyte (Santos et al., 2019). However, the regulation of JH is complex (Santos et al., 2019) and the key actors involved in JH-mediated pesticide-enhanced reproduction remains an open question. Previous finding reported that the GluCl receptor is involved in JH biosynthesis in cockroach (Chiang et al., 2002a; Chiang et al., 2002b; Liu et al., 2005). However, it is still unknown whether insecticides that target the GluCl can increase reproduction in pest insects and how the GluCl receptor regulates JH production.

In this study, we used a range of diverse approaches to investigate the impact of sublethal doses of avermectins on BPH fecundity, and unravel the molecular mechanisms mediating enhanced reproduction following exposure to this insecticide class. We show that avermectin exposure results in profound changes in the expression of a key suite of genes that in combination regulate JH, resulting in increased JH titer in adult females, which promotes fecundity. Interestingly, we found that GluCl channel modulators can inhibit AstA/AstAR expression and thereby promote JH production.

Results

GluCl allosteric modulators (emamectin benzoate and abamectin) stimulate fecundity of female N. lugens

To investigate whether GluCl modulators affect fecundity in BPH, we first determined the sub-lethal and median lethal doses of emamectin benzoate (EB) on 4^th instar nymphs, as well as on newly emerged males and females of BPH (Table S1). For this we employed two different bioassay methods, the rice seedling dip bioassay method and topical application bioassay method (Wu et al., 2018a; Yao et al., 2017), in order to assess both the systemic and contact toxicity of these insecticides (Table S1). We then systemically treated 4^th instar nymphs of BPH, newly emerged males and females with the estimated LC₁₅ or LC₅₀ concentrations of EB and examined the fecundity of BPH after these individuals mated with treated or untreated individuals. We use the term “t” to represent individuals treated with EB and “ck” to indicate control individuals that were treated with insecticide diluent minus insecticide. First, we tested whether treatment of BPH with EB in nymphal stage impairs the fecundity of BPHs at the adult stage. For 4^th instar nymphs of BPHs, after treatment with the LC₁₅ and LC₅₀ concentrations of EB, the number of eggs laid per female of BPHs in ♀t × ♂t crosses increased by 1.48 and 1.40 times compared with control ♀ck × ♂ck crosses (Figure 1A); the number of eggs laid per female of BPH in ♀t × ♂ck crosses increased by 1.53 and 2.07 times compared with control crosses (Figure 1B); However, the number of eggs laid by per female of BPH in ♀ck × ♂t crosses did not increase significantly compared to control ♀ck × ♂ck crosses (Figure 1C).

Table S1.

Figure 1.

Exposure of 4^th instar nymphs to the LC₁₅ and LC₅₀ concentrations of EB in contact bioassays also significantly stimulated fecundity. After treatment with the LC₁₅ and LC₅₀ concentrations of EB, the number of eggs laid by per female of BPH in ♀t × ♂t crosses increased by 1.18 and 1.26 times compared with the control (♀ck × ♂ck) (Figure 1D); The number of eggs laid per female of BPH in ♀t × ♂ck crosses increased by 1.27 and 1.56 times compared with the control crosses (Figure 1E); However, there was no significant difference in number of eggs laid between ♀ck × ♂t crosses and controls (♀ck × ♂ck) (Figure 1F). These results reveal that EB stimulates the fecundity of females following both the systemic and contact routes of exposure.

Next, we examined whether EB treatment of adult BPHs also stimulates reproduction. Indeed, treating newly emerged adults with the LD₁₅ and LD₅₀ concentrations of EB significantly stimulated the number of eggs laid per female (Figure 1-figure supplement 1A). Furthermore, sub-lethal exposure of 4^th instar BPH nymphs to another GluCl allosteric modulator, abamectin (LC₁₅ and LC₅₀ concentrations) was also found to significantly enhance reproduction (Figure 1-figure supplement 1B).

Figure 1-figure supplement 1.

To examine if EB also stimulates egg-laying in other insect species we conducted bioassays on the small brown planthopper, Laodelphax striatellus, the white backed planthopper, Sogatella furcifera and fruit flies Drosophila melanogaster. In contrast to our findings in BPHs, we found that sub-lethal doses (LC₁₅ and LC₅₀) of EB inhibits fecundity of female L. striatellus, (Figure1-figure supplement 2A-C) and has no impact on the fecundity of S. furcifera, (Figure1-figure supplement 2D-F). In addition, we found that sublethal doses (LC₁₅, LC₃₀ or LC₅₀) of EB also inhibit fecundity in D. melanogaster (Figure1-figure supplement 2G and H). These results indicate that the stimulation of reproduction by EB in BPH is species-specific and does not even extend to related insect species.

Figure 1-figure supplement 2.

The impact of EB treatment on BPH reproductive fitness

To better understand the effects of EB on the reproductive fitness of BPH, the duration of the preoviposition period, emergence rate, female ratio, female longevity and female weight were evaluated following exposure using the systemic bioassay. The preoviposition period of females treated with the LC₅₀ of EB decreased significantly compared with the control (Figure 1-figure supplement 3A). In contrast no significant effects of EB on emergence rate and female ratio were observed (Figure 1-figure supplement 3B and C). Female survival analysis showed that exposure of 4^th instar nymphs to the LC₅₀ of EB has no impact on female longevity (Figure 1-figure supplement 3D). Interestingly, occurrence of brachypterism (short-wing phenotype) in females and female weight were significantly increased after EB exposure (Figure 1-figure supplement 3E and F). It was reported that short-winged females usually display higher fecundity than long-winged females in BPHs and two different insulin receptors determine the alternative wing morphs (Xu et al., 2015). However, whether insulin signaling is involved in the fecundity increase of BPHs in our experiments needs further investigation.

Figure 1-figure supplement 3.

EB promotes ovarian maturation in BPH

To investigate the cause of increased egg-laying following EB exposure we examined if EB influences ovarian maturation in BPH. We dissected and compared the ovaries of females treated with the LC₅₀ of EB at 1, 3, 5 and 7 days after eclosion (DAE) with control females. At 3, 5 and 7 DAE, the number of eggs in the ovary of BPH in the EB treated group were significantly higher than that of controls (Figure 2A and B). We also observed that EB treatment of female adults also increases the number of mature eggs in the ovary (Figure 2-figure supplement 1). Hence, we found that EB does not affect normal egg developmental stages. However, our results suggest that EB treatment can promote oogenesis and as a result the insects both produce more eggs in the ovary and more eggs were laid.

Figure 2.

Figure 2-figure supplement 1.

Unlike the model insect Drosophila melanogaster, studies on oogenesis in the BPH are relatively limited. Therefore, we used immunohistochemistry to stain eggs at different stages of BPH development. We captured images with a laser confocal microscope, to document each stage of oogenesis in this species. We found that oogenesis in the BPH differs from that in D. melanogaster and belongs to the telotrophic meroistic type, where the oocyte is connected to nutrient-providing cells in the anterior trophic region of the ovariole through a nutritive cord. Oogenesis in the BPHs occurs in three main compartments (Figure 2-figure supplement 2):

1. Terminal Filament (TF): This region, composed of terminal filament cells, is similar to the terminal filament structure found in Drosophila.

2. Germarium: This region is where germline stem cells (GSCs) differentiate into oocytes and follicle cells. It typically consists of cap cells adjacent to the terminal filament, germline stem cells, escort cells, and cystoblasts.

3. Egg Chamber: The egg chamber in the BPH primarily contains eggs at six developmental stages. The distinct layering of follicle cells, visible through cytoskeletal staining, is the main feature distinguishing the six different developmental stages listed in the figure legend (see also Figure 2-figure supplement 2).

We next explored whether EB treatment could enhance or impair oogenesis in BPH. However, dissection of various developmental stages revealed that emamectin benzoate treatment has no significant negative effects on developmental stages and morphology of oogenesis in BPH (Figure 2C and D). Considering the number of eggs laid by EB treated females was larger than in control females (Figure 1 and Figure 1-figure supplement 1), our data indicated that EB treatment of BPHs can both promote oogenesis and fecundity (egg-laying) in BPHs.

Figure 2-figure supplement 2.

Stage 1: Connected to the germarium, the egg is trapezoidal. Stage 2: The egg is larger than in Stage 1 and is square to rectangular. Stage 3: The egg begins to elongate toward both poles, taking on an oval shape. Stage 4: The egg continues to enlarge and elongate, becoming spindle-shaped. Due to rapid growth, follicle cells on the egg surface start dividing, indicating nuclear division. Stage 5: The egg further enlarges, and most follicle cells display binucleation, the most notable feature of this stage. Small amounts of lipids start entering the egg. Stage 5 is the final stage before maturation, after which the egg enters the pedicel, a structure primarily formed by muscle fibers. The expanded mature egg stretches the muscle fibers in the pedicel, increasing the spacing between them. Stage 6: This is the mature egg stage, which occurs after the egg passes through the pedicel. During this stage, the egg undergoes rapid growth, continuously absorbing vitellogenin and lipids (Nile red staining confirmed that Stage 6 has the highest lipid content). This absorption causes the egg’s surface color to transition from the transparency seen in Stage 5 to mature eggs’ dark, opaque appearance. Once the egg enters the pedicel and becomes completely opaque, it reaches full maturity. The mature egg then moves into the lateral oviduct for oviposition, completing the egg-laying process.

EB exposure increases circulating sugars and the storage of macromolecules in BPHs

Nutritional status is an important indicator of reproductive fitness. Thus, to investigate whether EB affects intermediary metabolism and energy storage of BPHs, glycogen, triacylglyceride (TAG), total protein content, cholesterol and four circulating carbohydrates were quantified in 4^th instar BPH nymphs following exposure to the LC₅₀ of EB.

We found that EB exposure has no impact on glycogen levels (Figure 2-figure supplement 2A). The amount of TAG in EB-treated BPHs was 27% higher (p < 0.05) than those in controls, but only in BPH of the late fifth instar (5L) stage, with no significant differences observed in subsequent developmental stages (Figure 2-figure supplement 2B). The amount of total protein content in EB-treated BPH was higher than the control groups in the case of all developmental stages from 5L nymph to 7DAE (Figure 2-figure supplement 2C). EB exposure also increased cholesterol levels at 4 and 5 DAE (Figure 2-figure supplement 2D). Compared with the solvent control, EB treatment caused significant increases (p < 0.05) in the levels of sucrose, glucose, fructose, and trehalose (Figure 2-figure supplement 2E-H). Thus, collectively, these data provide evidence that EB exposure leads to energy mobilization and the metabolism of carbohydrates and lipids in BPHs.

Figure 2-figure supplement 2.

EB stimulates egg-laying that is mediated by the JH signaling pathway

Given the important role of juvenile hormone (JH) in vitellogenesis and egg development in insects (Jing et al., 2021; Ling and Raikhel, 2021; Luo et al., 2021; Santos et al., 2019; Zheng et al., 2022), we asked whether EB-treatment could influence the titer of JH in BPH. As measured by ELISA, the juvenile hormone titer of BPH nymphs treated with the LC₅₀ concentration of EB was significantly lower than that of controls in systemic bioassays during the middle and late stages of the 4^th instar (Figure 3A). However, at 2, 3 and 4 DAE, the JH titer in the EB treated group was significantly higher than that of the control (Figure 3A). Interestingly, the titer of another important insect hormone, the steroid ecdysone, was not significantly different between EB-treated BPH and solvent-treated controls (Figure 3-figure supplement 1). To independently validate the results of ELISA assays, we employed HPLC-MS/MS to measure JH titer in BPH following EB exposure (Bownes and Rembold, 1987; Guo et al., 2020; Jing et al., 2021). The results showed that the JH III titer significantly decreased after EB-treatment at the late 4^th instar nymph stage (Figure 3A and B), but significantly increased at the third day after eclosion (3 DAE) (Figure 3A and C). We also investigated the effects of EB treatment on the JH titer of female adults. The data indicates that the JH titer was also significantly increased in the EB-treated female adults compared with controls (Figure 3-figure supplement 3A). However, the steroid ecdysone, was also not significantly different between EB-treated BPH and controls (Figure 3-figure supplement 3B).

Figure 3.

Figure 3-figure supplement 1.

Figure 3-figure supplement 2.

Figure 3-figure supplement 3.

To further investigate the role of JH in EB-enhanced fecundity in BPHs, we treated BPHs with methoprene and pyriproxyfen, JH analogues or biologically active JH mimic respectively, to determine whether they can stimulate fecundity in BPH. Both compounds significantly increased egg-laying in BPHs (Figure 3D). Taken together these results reveal that EB stimulates an increase in JH titer and that JH signaling enhances fecundity in BPHs.

Since we found that EB could induce JH synthesis in the BPH, we asked whether EB could influence the expression of genes that are involved in JH biosynthesis or degradation. For this we treated 4^th instar nymphs or newly emerged female adults of BPHs with the LC₅₀ concentration of EB using systemic bioassay, and then collected early (5E), middle (5M) and late (5L) stage of 5^th instar nymph and 1-7 DAE female adults for analysis. Quantitative PCR was then used to examine the expression of key genes previously implicated in the regulation of JH (Gospocic et al., 2021; Nouzova et al., 2018).

Farnesoic acid O-methyltransferase (FAMeT) (Liu et al., 2010), is known as an important enzyme in the JH biosynthetic pathway, catalyzing methylation of farnesoic acid (FA) to methyl farnesoate (MF) (Nouzova et al., 2018). We found that this gene was significantly upregulated in 5M instar nymph to 4 DAE after EB-treatment (1.0-fold to 3.0-fold) (Figure 3E). Juvenile hormone acid methyltransferase (JHAMT) (Nouzova et al., 2018; Shinoda and Itoyama, 2003), which is involved in the biosynthesis of JH, was also upregulated in EB-treated BPH at 1 to 5 DAE compared to controls (1.5-fold to 3.0-fold) (Figure 3F). We furthermore found that EB treatment in female adults stimulates JHAMT expression (Figure 3-figure supplement 3C).

Methoprene-tolerant (Met), belongs to the basic helix–loop–helix Per/Arnt/Sim (bHLH-PAS) family of transcription factors and is the intracellular (nuclear) receptor of JH (Riddiford, 2020; Zhu et al., 2019). The levels of met mRNA increased in EB-treated BPH at the 5M and 5L instar nymph and 1 to 5 DAE adult stages compared to controls (1.7-fold to 2.9-fold) (Figure 3G).

Krüppel homolog 1 (Kr-h1), a transcriptional target of JH signaling, is reported to be sensitive to levels of JH and its expression levels are directly correlated with JH titers (Gospocic et al., 2021; Jindra et al., 2013; Zhang et al., 2022a). We found that Kr-h1 was significantly upregulated in the adults of EB-treated BPH at the 5M, 5L nymph and 4 to 5 DAE stages (4.7-fold to 27.2-fold) when 4^th instar nymph and female adults were treated with EB (Figure 3H and Figure 3-figure supplement 3D).

Similarly, the expression levels of vitellogenin (Vg), a key downstream component of JH signaling triggering ovarian development in insects including BPH (Shen et al., 2019), was markedly increased in females at 2–5 DAE by EB (1.7-fold to 5.5-fold) (Figure 3I). During the maturation process, eggs require uptake of vg, and an increase in vg content can accelerate egg maturation, producing more mature eggs. Our molecular data suggest that EB treatment leads to an upregulation of vg expression. Based on these findings, we conclude that the increase in egg-laying caused by EB treatment is due to the upregulation of vg (Figure 3I), which raises Vg content, promoting the uptake of Vg by maturing eggs and resulting in the production of more mature eggs.

Juvenile hormone esterase (JHE) is the primary JH-specific degradation enzyme that plays a key role in regulating JH titers (Kamita et al., 2003). Interestingly, we observed a significant upregulation of JHE mRNA levels in the early and middle 5^th instar nymph stage followed by downregulation in 5L instar nymphs of EB-treated BPH (1.3-fold to 2.6-fold) (Figure 3J).

In combination these results reveal that EB has a profound impact on the expression of key genes involved in the biosynthesis of JH or in downstream signaling pathway genes that might promote egg development and female fecundity.

To further understand whether these JH pathway-related genes were involved in egg-laying behavior in BPH, we performed RNAi experiments to downregulate the expression of Kr-h1 (Figure 3-figure supplement 2A). We found that silencing of Kr-h1 diminished vg gene expression (Figure 3-figure supplement 2B). Importantly, diminishment of JHAMT, Met and Kr-h1 genes in female BPHs was also found to suppress egg-laying (Figure 3K-M). However, this phenotype was rescued by EB treatment, which restored egg-laying to normal levels in BPH injected with JHAMT, Met and Kr-h1 dsRNA (Figure 3K-M). Note that this rescue is possible since knockdown of the genes is incomplete when using dsRNA injection (and residual gene expression allows for EB action). Together these results provide a mechanistic understanding of how EB enhances fecundity in BPH by modulating the expression of key genes involved in JH synthesis (Figure 3N).

EB induces JH biosynthesis through the peptidergic AstA/AstAR signaling pathway

The timing and level of JH biosynthesis can be precisely regulated by neuropeptides, stimulatory allatotropins (ATs) and inhibitory allatostatins (Asts), in many insects (Bellés et al., 1999; Kataoka et al., 1989; Kramer et al., 1991; Lorenz et al., 1995; Stay and Tobe, 2007; Verlinden et al., 2015; Wegener and Chen, 2022; Woodhead et al., 1989). Insects can, in a species-specific manner, produce one type of AT and three types of Asts: FGL-amide Ast (AstA) (Woodhead et al., 1989; Yin et al., 2006); myoinhibitory peptide (MIP or AstB) (Stay and Tobe, 2007) and PISCF Ast (AstC) (Wang et al., 2012). In some species, there also exists two paralogue genes of AstCs which are named AstCC and AstCCC (Veenstra, 2009, 2016). Interestingly, the allatostatic activity of these three types of Ast peptides varies between insect species so that in each species only one type of Ast (for example AstA) controls JH production (Nässel, 2002; Stay and Tobe, 2007; Wang et al., 2012; Wegener and Chen, 2022).

Analysis of the BPH neural transcriptome sequence data has revealed the presence of one AT, four types of Asts and four corresponding receptors, allatotropin receptor (A16, ATR), AstA receptor (A2, AstAR), AstB (MIP) receptor (A10, AstBR or MIPR) and AstC receptor (A1, AstCR) (Tanaka et al., 2014). We cloned the five neuropeptide genes (AT, AstA, AstB/MIP, AstCC and AstCCC) and confirmed the sequence obtained from transcriptome data (Figure 4-figure supplement 1) (Tanaka et al., 2014). Interestingly, we found that AstC is missing in the genome of BPH and only AstCC and AstCCC are present (Figure 4-figure supplement 1). Next, we also cloned their corresponding receptors (Veenstra, 2009) including ATR (A16), AstAR (A2), AstBR (A10) and AstCR (A1) which might be activated by AstCC and/or AstCCC (Audsley et al., 2013; Veenstra, 2009; Zhang et al., 2022a). Sequence alignments and phylogenetic analysis are shown in Figure 4-figure supplement 2.

Quantitative PCR was then used to examine if EB treatment could influence the expression of the genes encoding these neuropeptides and their receptors. Treating 4^th instar or female adults BPH with the LC₅₀ concentration of EB significantly increased the expression of AT, ATR and AstCCC while resulting in the downregulation of AstA, AstB/mip, AstCC, AstAR and AstBR/mipr at the adult stage (Figure 4A and B, Figure 4-figure supplement 3A-G). Among these, AstA and AstAR were the most downregulated genes after EB treatment (Figure 4A and B, Figure 4-figure supplement 3H) and thus the AstA/AstAR signaling system was selected for subsequent functional analysis. Silencing of AstAR in female BPH using RNAi (Figure 4C), resulted in an increased number of eggs laid per female compared with dsgfp-injected controls (Figure 4D). Interestingly, silencing AstAR also resulted in the upregulation of JHAMT, Met and Kr-h1 which are involved in the JH biosynthesis/signaling (Figure 4E-H). However, JHE was not influenced by AstAR silencing (Figure 4I). We next investigated whether silencing the AstAR gene would influence the JH titer in BPHs. A significantly increased JH titer was observed in AstAR silenced BPH compared with controls (Figure 4J). Thus, our data strongly suggest that AstA is a key inhibitor of JH production in BPH.

Figure 4.

Figure 4-figure supplement 1.

Figure 4-figure supplement 2.

Figure 4-figure supplement 3.

Finally, we investigated whether injection of mature Ast and AT peptides could influence the number of eggs laid and JH titer in BPH. We synthesized one AT, six AstAs (AstA1 to AstA6), one AstCC and one AstCCC peptide according to their determined sequences (Figure 4-figure supplement 1). Indeed, we found that injection of AstA1 and AstA6 in female adults reduced the number of eggs laid per female over 48 h (Figure 4K). We also observed that AstA1 injection decreased the JH titer as measured 16h and 48h after injection (Figure 4L), and AT injection increased the JH titer after 2h, but levels returned to normal 4h after injection (Figure 4M). Collectively, these results provide compelling evidence that EB induces reproduction in BPH through the AstA/AstAR and JH signaling pathways (Figure 4N) and further supports the role of AstA and AT in regulation of JH titer in this species.

EB-enhanced fecundity in BPHs is dependent on its molecular target protein, the GluCl channel

EB and abamectin, are allosteric modulators, which target glutamate-gated chloride channels (GluCls) (Sparks et al., 2020; Sparks et al., 2021; Wu et al., 2017). Hence, we examined whether EB-stimulated fecundity in BPH is mediated by its molecular target GluCl. The full length GluCl coding sequence from BPH was cloned and sequenced (Figure 5-figure supplement 1) and the impact of EB on GluCl gene expression examined using quantitative PCR. Treatment of BPH with the LC₅₀ concentration of EB significantly downregulated GluCl gene expression at the 5E and 5M nymph stages, while upregulating GluCl gene expression at 2 DAE and 5 DAE in the adult stage (Figure 5A). We also found that EB treatment in female adults also increases the GluCl gene expression (Figure 5-figure supplement 3). To examine the role of GluCl gene in BPH fecundity, RNAi was used to knockdown expression of this gene in female adult BPH (Figure 5B). A significant decrease in the number of eggs laid by per female was observed in dsGluCl-injected insects compared with dsgfp-injected insects (Figure 5C). However, treatment with EB was found to rescue the decreased egg-laying phenotype induced by dsGluCl injection (Figure 5C). Note that this rescue is possible since knockdown of the genes is incomplete when using dsRNA injection (and residual gene expression allows for EB action). To investigate the mechanism by which GluCl expression modulates fecundity we examined if silencing GluCl influences JH titer and JH-related gene expression. Indeed, we observed that RNAi knockdown of GluCl leads to a decrease in JH titer (Figure 5D) and down-regulation of genes including JHAMT which is responsible for JH synthesis, and the JH signaling downstream genes Met and Kr-h1 (Figure 5E-G). In contrast, expression of JHE was not changed in the GluCl knockdown insects (Figure 5H and I). We also examined whether silencing GluCl impacts the AstA/AstAR signaling pathway in female adults. Silencing GluCl in female adults was found to have no impact on the expression of AT, AstA, AstB, AstCC, AstAR, and AstBR. However, the expression of AstCCC and AstCR was significantly upregulated in dsGluCl-injected insects (Figure 5-figure supplement 2A-H). These results suggest that EB activates GluCl which induces JH biosynthesis and release, which in turn stimulates reproduction in BPH (Figure 5J).

Figure 5.

Figure 5-figure supplement 1.

Figure 5-figure supplement 2.

Figure 5-figure supplement 3.

Discussion

Pesticide-induced resurgence of pest insects is a serious problem in rice and several other crops (Wu et al., 2020). However, the mechanisms underpinning pesticide-enhanced reproduction in insects remain poorly understood. Here we reveal that a suite of molecular actors underlie this trait that, in combination, mediate profound physiological changes in the reproductive fitness of female BPHs. Our data provide fundamental insights into the molecular mechanisms by which xenobiotics modify insect reproduction and have applied implications for the control of a highly damaging crop pest. We discuss these topics below.

Sublethal doses of the GluCl modulators, EB and abamectin, stimulates fecundity in BPH

We show that in both contact and systemic assays the avermectins EB and abamectin stimulate reproduction in BPHs. Thus, insecticide-enhanced reproduction is likely a key factor in the BPH resurgence observed when farmers use EB and abamectin to control leaffolders in rice crops in China (Yao et al., 2017). Although this is the first report of sublethal doses of avermectins enhancing insect fecundity, our findings are consistent with previous studies, which have shown that certain insecticides, herbicides and fungicides stimulate BPH reproduction (Cheng et al., 2014; Ge et al., 2011; Ge et al., 2013; Wan et al., 2013; Wang et al., 2010; Wu et al., 2020; Yang et al., 2019; Zhang et al., 2014; Zhao et al., 2011). Intriguingly, we show that EB only induces fecundity in female adults and is specific to BPH, with EB exposure failing to enhance reproduction in two related species, the small brown planthopper, L. striatellus and the white backed planthopper, S. furcifera, or the model insect D. melanogaster. Thus, the mechanisms underpinning this trait appear to act exclusively on female BPHs and may be specific to this species. Nevertheless, our findings are of general interest since they link the molecular EB target, the GluCl channel and juvenile hormone signaling to fecundity.

Pesticides may stimulate insect reproduction through a variety of physiological and molecular mechanisms. Our data reveal that exposure to sub-lethal concentrations of EB results in profound changes to female BPH fitness, leading to increases in female weight, contents of total protein, cholesterol, and sugar, as well as egg production and decreases in duration of the preoviposition period. Some of these findings exhibit parallels with previous studies, which demonstrated that treating third-instar BPH nymphs with either deltamethrin, triazophos, or imidacloprid led to increased soluble sugar levels in the corresponding adults (Yin et al., 2014). Such metabolites provide the energy that drives BPH reproduction and hence resurgence (Wu et al., 2020). Thus, together with prior work, our results suggest that pesticides associated with resurgence stimulate nutritional intake in BPH to fuel enhanced energy-intensive reproduction.

The JH signaling pathway plays an essential role in EB-induced fecundity in BPH

JH is a pleiotropic hormone which plays important roles in development and reproduction in insects (Riddiford, 2008; Santos et al., 2019). Circulating JH titers are regulated by factors that control JH production in the corpora allata including biosynthetic enzymes and catabolic enzymes that regulate JH levels. Our results show that EB increases circulating JH III titers in BPH females over 2–4 days after eclosion (DAE) and promotes ovarian development. Previous studies have reported that triazophos and deltamethrin treatments also lead to increased circulating JH III titers in BPH females over 1–3 days post emergence. Similarly, jinggangmycin treatments were found to lead to increased JH titers (by approximately 45–50%) in BPH females over two days post emergence (Xu et al., 2016). Thus, our findings, in combination with these previous studies, demonstrate that insecticide treatments can have dramatic effects on the regulation of key insect hormones involved in pest reproduction, which can in turn drive pesticide resurgence.

Although increased JH titers following pesticide exposure have been correlated with reduced levels of active JH esterase during the first three days PE (Ge et al., 2010), the type and number of mechanisms mediating the observed increase in hormone titer has remained an open question. Our data reveal that elevated JH titer in EB-exposed BPH is associated with the upregulation of genes that encode biosynthetic enzymes for JH (JHAMT) and downstream signaling genes that can induce vg gene expression (kr-h1). Using RNAi we provide functional evidence of the role of these genes in the regulation of JH III and fecundity of female BPHs, and demonstrate that EB can restore the reduction in egg production resulting from the knockdown of JHAMT, met and kr-h1 expression.

JHMAT is an enzyme that catalyzes the conversion of inactive precursors of JH to active JH in the final stages of JH biosynthesis (Nouzova et al., 2018; Shinoda and Itoyama, 2003). Interestingly, while it has not been previously implicated in pesticide resurgence, treatment of the stored product pest Sitotroga cerealella with diallyl trisulfide, an insecticidal compound in garlic essential oil, was found to increase JHAMT mRNA levels (Shah et al., 2022). Because JHMAT is the key rate-limiting enzyme in regulation of JH titer our results suggest that its enhanced expression is a key molecular mechanism of pesticide resurgence in BPH.

Met is a ligand-activated member of the basic helix–loop–helix Per/Arnt/Sim (bHLH-PAS) transcription factors and is the intracellular receptor for JH (Riddiford, 2020; Zhu et al., 2019). Kr-h1 is a zinc finger protein that acts downstream of Met and is expressed in response to JH signaling. Although the genes encoding these proteins have not been previously linked to pesticide resurgence, our finding that they are upregulated following EB exposure, and demonstration of their role of in promoting fecundity, is consistent with previous studies (Santos et al., 2019). Specifically, treatment of BPH with JH III or the insecticidal analogs methoprene or pyriproxifen was found to induce the expression of Kr-h1 (Jin et al., 2014). Furthermore, knockdown of Met and Kr-h1 in BPH brachypterous females was found to result in delayed ovariole maturation and this was significantly more pronounced than the response observed in BPH treated separately with dsNlMet or dsNlKr-h1 (Lin et al., 2015). This finding provides evidence of a possible interaction between Met and Kr-h1 and, in combination with our data, suggests that Met and Kr-h1 may act in concert to mediate EB-enhanced fecundity. However, our results showed that EB treatment can mildly increase (about 2-fold) expression of the Met gene in brown planthoppers (Figure 3G). Our data furthermore indicated that Met and FAMeT expression levels were not drastically influenced by EB compared with kr-h1 and vg (Figure 3H and I). It should be mentioned that JH action does not directly result in the increase of Met. However, we cannot rule out the possibility of other factors (indirect effects), induced by EB treatment, increases the mRNA expression level of Met. One recent paper reported that downregulation of transcription factor CncC will increase met expression in beetles (see Figure 6A in the reference) (Jiang et al., 2023). And many studies have reported that insecticide treatment activates the CncC gene signaling pathway, which regulates detoxification gene expression (Amezian et al., 2023; Fu et al., 2024; Hu et al., 2021). Hence, it is possible that EB might influence the CncC gene pathway, which in turn induces met expression. This EB effect on met upregulation may be similar to the upregulation of GluCl and some other secondary effects (see below).

Figure 6.

It must be mentioned that although we exposed the BPH to EB at 4th instar, we make the insect feed on the EB-treated rice plants for four days. After that, the insect will develop into 5th instar, the final nymphal stage of brown planthopper. Since BPH do not have a pupal stage, this might cause the EB presented to the insects last a longer time even in the adult stage. Besides this, we found that EB treatment will increase the weight of adult females, which indicates that EB might increase food intake in BPHs that might produce more insulin peptide. Hence, insulin might increase the JH synthesis at the adult stage.

EB-induced fecundity in BPHs is dependent on the allatostatin A signaling pathway

In addition to regulatory proteins that promote JH production, insects utilize peptides that inhibit JH biosynthesis. These include the allatostatins: FGLamides (FGLa; AstA), the W(X)6Wamides (AstB), and the PISCFs (AstC) (Nässel, 2002; Stay et al., 1996; Stay and Tobe, 2007; Tanaka et al., 2014; Wegener and Chen, 2022; Zhang et al., 2021). Interestingly, our results showed that EB exposure results in a marked downregulation of the expression of the genes encoding allatostatin AstA and its receptor AstAR. We provide evidence for the functional impact of this on JH synthesis and BPH fecundity by: i) demonstrating that RNAi knockdown of AstAR expression results in increased JH titer and enhanced egg production, and, ii) showing that injection of female BPH with synthetic AstA peptide reduces JH titer and decreases egg production. Thus, our data provide unequivocal evidence that AstA is a key inhibitor of JH production in BPH. This finding is consistent with previous work which has shown that FGLa/ASTs (AstA) inhibit JH biosynthesis in cockroaches, and termites (Woodhead et al., 1989; Yagi et al., 2005). To our knowledge, our study is the first report of insecticides down-regulating the expression of the neuropeptide receptor, AstAR, and linking this to increases in JH titer and enhanced reproduction in insects.

Interestingly knockdown of AstAR resulted in significant increases in the expression of genes involved in JH synthesis/signaling including JHAMT, Met and Kr-h1. Related to this finding, previous work has shown that knockdown of the AstA receptor gene, Dar-2, in D. melanogaster results in changes in the expression of genes encoding Drosophila insulin-like peptides (DILPs) and adipokinetic hormone (AKH) (Hentze et al., 2015). Together with our findings, this demonstrates that AstA receptors may modulate the expression of numerous downstream genes involves in metabolism, energy store and reproduction. In the case of pesticide resurgence our results imply significant cross-talk in the expression of genes that inhibit JH production and those that promote it.

The GluCl plays an essential role in EB-induced fecundity in BPH

EB and abamectin are allosteric modulators of GluCls (Sparks et al., 2020; Sparks et al., 2021; Wu et al., 2017). Our data revealed that EB exposure increases the expression of the GluCl in adult stage of BPH, and knockdown of GluCl expression resulted in a reduction in both JH levels and egg production and egg laying. There are several studies showing that insects treated with insecticides display increases in the expression of target genes. For example, the relative expression level of the ryanodine receptor gene of the rice stem borer, Chilo suppressalis was increased 10-fold after treatment with chlorantraniliprole, an insecticide which targets the ryanodine receptor (Peng et al., 2017). Besides this, in Drosophila, starvation (and low insulin) elevates the transcription level of the short neuropeptide F and tachykinin receptors (Ko et al., 2015; Root et al., 2011). In BPHs, reduction in mRNA and protein expression of a nicotinic acetylcholine receptor α8 subunit is associated with resistance to imidacloprid (Zhang et al., 2015). RNA interference knockdown of the α8 gene decreased the sensitivity of N. lugens to imidacloprid (Zhang et al., 2015). Hence, expression of receptor genes may be regulated by diverse factors, including insecticide treatment. In our case, we found that EB upregulates its target gene GluCl. However, we do not suggest that EB functions as transcriptional activator of GluCl and we have no data to suggest a mechanism for how EB treatment changes the expression of GluCl in the BPH. Considering that EB exposure in our experiments is lasting several days, it might be an indirect (or secondary) effect caused by other factors downstream GluCl that affect channel expression Possibly the allosteric interaction with GluCl by EB renders the channel dysfunctional and the cellular response is to upregulate the channel/receptor to compensate.

Interestingly, the GluCl has been reported to inhibit the biosynthesis of JH in the cockroach, Diploptera punctata (Liu et al., 2005). Recent work has also reported that modulation of glutamatergic signals may contribute to the photoperiodic control of reproduction in bean bug, Riptortus pedestris (Hasebe and Shiga, 2022). Furthermore, work on D. punctata has revealed that application of the GluCl channel agonist ivermectin, which like EB belongs to the avermectin family, caused a decline in JH synthesis in the corpus allatum (Liu et al., 2005). While the inhibitory effect of ivermectin observed in this previous study differs from the activating effect of EB we observed in our study, it is consistent with our finding of a role for the GluCl channel in the regulation of JH regulation. Interestingly, we found that knockdown of the GluCl gene expression results in the down-regulation of JHAMT, Met and Kr-h1, further revealing significant convergent relationships between genes underpinning pesticide resurgence.

Conclusion

Our study has revealed that exposure of BPHs to EB affects a diverse suite of genes that act in combination to enhance JH titer and thus increase female fecundity. A schematic of how these factors promote ovarian maturation in the adult stage of N. lugens through the AstA and JH signaling pathway is provided in Figure 6. Our findings provide the foundation for further work to understand how these genes interact and the mechanisms by which their expression is activated or repressed by EB. Furthermore, our findings provide fundamental insights into the molecular response in insects to xenobiotic stress and illustrate that pesticides can have unexpected and negative impacts on pest populations. In this regard our findings also have applied implications for the control of a highly damaging crop pest. Previous studies have reported that avermectins such as abamectin are toxic to the wolf spider Pardosa pseudoannulata, which is the main predator of BPH in rice crops (Chen et al., 2017; Sogawa, 2015). Thus, these insecticides both stimulate reproduction in BPH while killing their natural enemies providing a ‘perfect storm’ for inducing damaging BPH outbreaks. Based on these findings, to avoid BPH resurgence, we suggest that the agents, EB and abamectin, should not be (or rarely be) applied to rice plants at growth stages when BPHs are present. On a more optimistic note, our findings have identified numerous genes that play key roles in BPH reproduction and thus represent promising targets for the development of novel controls against this important pest.

Materials and methods

Insects

BPH was initially collected from Wanan, JiangXi Province in 2020, reared on ‘Taichung Native 1’ (TN1) rice seedlings in the laboratory without exposure to any insecticides. The strain was maintained in a climatic chamber at 27 ± 1°C, with relative humidity of 70 ± 10% and a light: dark = 16 h: 8 h photoperiod.

Chemicals

Emamectin benzoate (95.2%) was obtained from Hebei Weiyuan Co., Ltd. (Hebei, China). Abamectin (96.8%) was obtained from Hebei Weiyuan Co., Ltd. (Hebei, China). Pyriproxyfen (98%) was obtained from Shanghai Shengnong Co., Ltd. (Shanghai, China). Methoprene (S)-(+) (mx7457-100mg) was purchased from Shanghai MaoKang Biotechnology Co., Ltd., (Shanghai, China). Juvenile hormone standard sample (J912305-10mg) was purchased from Shanghai Macklin Biotechnology Co., Ltd., (Shanghai, China).

Bioassay

Different life stages of insects were used to perform bioassay to investigate the effects of insecticide on nymphs and adults. To test whether treatment of the nymph stage of insects would promote reproduction in female, we used 4^th instar nymphs of BPH or 3^rd instar nymph of Laodelphax striatellus and Sogatella furcifera to perform bioassays. To test whether treatment of the adult stage of insects would promote reproduction in female, we used newly emerged male and female BPH.

Systemic route

The rice-seeding dipping bioassay method was used to evaluate the susceptibility of BPH, L. striatellus and S. furcifera to EB. Technical-grade insecticides were dissolved in acetone as stock solution then diluted in a series of six concentrations with water containing 0.1% Triton. Selected rice seedlings at the 6-8 cm growth stage were dipped in insecticide solutions for 30 s and then air-dried at room temperature. The roots of the rice seedlings were wrapped with cotton strips and placed seedlings placed in a plastic cup 5 cm in diameter. Fifteen insects were introduced into each plastic cup for each replicate. The top of the cup was then selaed with gauze to prevent escape. All experiments comprised at least three biological replicates. Control rice seedings were treated with 0.1% Triton X-100 water solution only. All treatments were maintained under standard conditions of 27 ± 1 °C and 70–80% relative humidity with a 16 h light/8 h dark photoperiod. Mortality was assessed after 4 d for N. lugens or 2 d for L. striatellus and S. furcifera after treatment with insecticides. The insects were considered dead if they were unable to move after a gentle prodding with a fine brush.

For Drosophila larvae bioassay, we adopted a method described previously in our lab with minor modifications (Huang et al., 2020). Briefly, twenty third instar larvae were placed in fly vials containing fly food (based on corn powder, brown sugar, yeast and agar) supplemented with EB of different concentrations. Four concentrations (LC₁₀, LC₃₀ and LC₅₀) were tested together with a negative (no insecticide) control. For Drosophila adult bioassays, we selected virgin females three days after eclosion. Several concentrations were overlaid onto fly food in standard Drosophila vials and allowed to dry overnight at room temperature. 15 adult flies (three days after eclosion) were then added to each vial and mortality assessed after 2 d. Four replicates were carried out for each concentration. Control mortality was assessed using vials containing food with solvent minus insecticide.

Contact route

For topical bioassays working insecticide solutions were prepared in acetone. 4^th instar nymphs or newly emerged males/females were anesthetized with carbon dioxide for 5 s, and then 0.04 μl/insect test solution applied topically to the dorsal plates with a capillary micro-applicator (Yao et al., 2017). Insects were then placed in an artificial climate incubator with a temperature of 27±1°C, a photoperiod of 16:8 h (L:D), and a humidity of 70%±10%. Mortality was determined 2 d after treatment. Data with over 20% mortality in the control treatment were discarded, and the assay was repeated at least three times.

Fecundity assays

Fourth instar nymphs of BPH and 3^rd instar nymph of L. striatellus and S. furcifera were treated with the LC₁₅ and LC₅₀ of EB, and then transferred to fresh rice seedlings. After eclosion, the adults were used in the following experiment. Newly emerged treated adults and untreated adults were paired to produce four groups: untreated males and untreated females (♂ck×♀ck; ck indicates untreated); untreated males and treated females (♂ck×♀t; t indicates insecticide treated); treated males and untreated females (♂t×♀ck); treated males and treated females (♂t×♀t). Each group comprised at least 10 mating pairs. All pairs were transferred to glass cylinders (diameter 5 cm and height 20 cm) containing rice plants (25 days old almost 20 cm high) as a food source for eleven days. The number of eggs and nymphs in plants were counted by a microscope.

For Drosophila egg-laying assay, we adopted our previous method (Wu et al., 2019). Briefly, insecticide-treated virgin females were paired with untreated males for three days and then the mated females transferred into the Drosophila ovipositing apparatus. The number of eggs were counted after 16 hours.

Fitness analysis

The fitness of EB-treated BPH were analyzed using methods reported previously (Zeng et al., 2023). We selected two groups, (♂ck×♀ck) and (♂ck×♀t), to study the effects of the LC₅₀ concentration of EB on BPH fitness. In the case of systemic exposure, 4^th instar nymphs of BPH were treated with the LC₅₀ of EB for 4 days and then transferred to tubes containing untreated rice plants for individual rearing. The rice plants were replaced every three days with untreated plants. The emergence ratio, female ratio, preoviposition period, female longevity, brachypterism female ratio and female weight were calculated.

Examination of ovarian development

Adult females from ♂ck×♀ck control and ♂ck×♀t group on 1, 3, 5, 7 DAE were dissected for observe the ovarian development. The mature eggs in ovary were photographed and recorded. Each group has at least fifteen replicates.

To examination whether EB treated impairs egg maturation, we dissected untreated or EB-treated ovaries and fixed them in 4% paraformaldehyde (PFA) in phosphate buffered saline (PBS) for 30 min at room temperature. After four washes of 10 min (4 × 10 min) in PAT3 (PBS with 0.5% Triton X-100 and 0.5% bovine serum albumin), the ovaries were then incubated with DAPI (4’,6-diamidino-2-phenylindole, 100 nM) and Actin-stain 670 Fluorescent Phalloidin (200 nM). Imaging was performed using Zeiss LSM980 confocal laser microscope.

Measurements of glycogen, triglyceride, total protein content, cholesterol and four sugars

The content of glycogen, triglyceride, cholesterol and total protein was determined by spectrophotometry at 620 nm, 510 nm, 510 nm and 562 nm respectively using the glycogen assay kit (A043-1-1), triglyceride reagent kit (A110-1-1), cholesterol assay kit (A111-1-1) and total protein assay kit (A045-2-2) obtained from Nanjing Jiancheng Bioengineering Institute following to the manufacturer’s instructions. The determined results were normalized to the protein content in the sample, which was determined using BCA Protein Assay Reagent Kit (Thermo Scientific, Waltham, USA). Each sample contained tissue extracts from five adult female BPH, with three biological replicates per sample.

To assess the content of four sugars (sucrose, glucose, fructose and trehalose) in the extract of BPH tissue, the same extraction method was used as above. Sugar content was quantified using the colorimetric method by scanning at 290 nm, 505 nm, 285 nm and 620 nm respectively using the sucrose reagent kit (A099-1-1), glucose reagent kit (F006-1-1), fructose reagent kit (A085-1-1) and trehalose reagent kit (A150-1-1) obtained from Nanjing Jiancheng Bioengineering Institute based on the manufacturer’s instructions. Each sample contained tissue extracts from five adult female N. lugens, with three biological replicates per sample.

Determination of Juvenile hormone III and ecdysone titers of BPH by ELISA

The titer of Juvenile hormone III in BPH was measured using the Juvenile hormone ELISA Kit (Lot Number: HLE92086, Shanghai HaLing Biotechnology Co., Ltd., Shanghai, China) which employs the dual-antibody sandwich ELISA method. The titer of ecdysone in BPH were measured using the ecdysone ELISA Kit (Lot Number: ZK14705, Shanghai ZhenKe Biotechnology Co., Ltd., Shanghai, China). At least three biological replicates were employed for each treatment.

Determination of Juvenile hormone III titer in BPH using HPLC-MS/MS

The whole bodies of 5 individuals BPH were mixed with 1 ml of n-hexane, followed by centrifugation at 10,000×g for 10 min, the upper hexane layer was then dried with nitrogen, dissolved in methanol and sonicated for 10 min, after centrifugation at 10000×g for 10 min, the supernatant was collected through the organic filter membrane of 0.22 μm into 2 ml vials for JH LJ determination. JH LJ standard sample (J912305-10 mg) purchased from (Shanghai McLean Biochemical Technology Co. Ltd), dissolved in methyl alcohol as stock solution 10,000 mg/L was diluted in a series of six concentration gradients to serve as a reference. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) was then carried out using UPLC Xevo TQ-S Micro (Waters technology), quantitative method according to the external standard, the chromatographic column was EC-C18 (4.6 mm×150 mm, 2.7 μm), column temperature was 30°C, injection volume was 20 μl, elution flow rate was 0.3 ml/min, and the mobile phase was acetonitrile:formic acid water (90:10), detection wavelength was 218 nm, the peak height was used for quantification.

Cloning, sequence and phylogenetic analysis

The NCBI database and BLAST program were used to carry out sequence alignment and analysis. Open Reading Frames (ORFs) were predicted with EditSeq. Primers were designed using the primer design tool in NCBI. Total RNA Extraction was extracted from 30 adults BPH using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. cDNA was synthesized using the Biotech M-MLV reverse transcription kit. Full-length gene sequences were amplified by PCR using cDNA as template and Rapid Taq Master Mix (Vazyme Biotech, Cat# P222-02). The PCR product was purified on a 1% agarose gel, cloned into pClone007 Simple Vector Kit (Tsingke Biotech, Cat# TSV-007S), and then sequenced using the 3730 XL DNA analyzer (Applied Biosystems, Carlsbad, CA, USA). Table S2 contains a list of the primers used in this study.

Table S2

The exon and intron architectures of AT, AstA, AstB, AstCC and AstACCC were predicted based on the alignments of putative ORFs against their corresponding genomic sequences. Sequence similarity/annotations and orthologous gene searches were performed using BLAST programs available in NCBI. Multiple alignments of the complete amino acid sequences were performed with Clustal Omega (http://www.ebi.ac.uk/Tools/msa/clustalo). Phylogeny was conducted using the maximum likelihood technique to create phylogenetic trees and these were bootstrapped with 1000 replications were used using the MEGA 6 software (Tamura et al., 2013).

RNA interference

Double-stranded RNA (dsRNA) of gfp (green fluorescent protein), JHAMT (juvenile hormone acid O-methyltransferase,), Met (methoprene-tolerant), Kr-h1 (krüppel homolog 1), AstAR (allatostatin-A receptor) and GluCl (glutamate-gated chloride channel) was prepared using Ambion’s MEGAscript T7 kit instructions following the manufacturer’s instructions. The primer sequences for double-stranded RNA synthesis are listed in Table S2. Newly emerged females were injected with 40 nl (5,000 ng/μl) of double-stranded RNA of gfp (dsgfp) or double-stranded RNA of the target genes in the conjunctive part between prothorax and mesothorax of insects. In the RNAi experiments, BPH were then treated with the LC₅₀ of EB 24h after dsRNA injection and the whole body sampled for qRT-PCR analysis.

Quantitative RT-PCR

Fourth instar nymphs of BPH were treated with EB after which total RNA was extracted from 5^th instar nymphs and 1-7 day post-eclosion females of N. lugens using the methods detailed above. The HiScript® II RT SuperMix for qPCR (+gDNA wiper) kit from Vazyme, Nanjing, China, was used to create first-strand cDNA. Primer3 (http://bioinfo.ut.ee/primer3/) was used to design real-time quantitative PCR (qPCR) primers listed in Table S2. mRNA levels of candidate genes were detected by qPCR using the UltraSYBR Mixture (with ROX) Kit (CWBIO, Beijing, China). Each reaction contained 2 μL of cDNA template (500 ng), 1 μL each forward and reverse qPCR primer (10 μM), 10 μL of UltraSYBR mixture buffer, and 6 μL of RNase-free water. Q-PCR was run on an ABI 7500 Real-Time PCR System (Applied Biosystems, Foster City, CA) under the following conditions: 5 min at 95°C, followed by 40 cycles of 95°C for 10 s and 60°C for 30 s. Three independent biological replicates and four technical replicates were used in each qPCR experiment. The housekeeping genes 18S ribosomal RNA of BPH were selected to normalize the expression of candidate genes. The 2^−ΔΔCt method (Ct represents the cycle threshold) was used to measure relative expression levels (Livak and Schmittgen, 2001). Three biological replicates were used for statistical comparison between samples. Table S2 contains a list of the primers used in this study.

Statistics

PoloPlus v2.0 (LeOra Sofware 2008) was used to calculate the lethal concentration (LC₅₀) and 95% fiducial limits (95% F.L.). GraphPad Prism 8.0 software (GraphPad Software Inc., San Diego, USA) was used to generate graphs and perform statistical analysis of data. Data presented in this study were first verified for normal distribution using the D’Agostino–Pearson normality test. One-way analysis of variance (ANOVA) with Duncan’s multiple range test was used to test differences among multiple groups of normally distributed groups of data. Student’s t test was used to test the differences between two groups. If not normally distributed, Mann–Whitney test was used for pairwise comparisons, and Kruskal–Wallis test was used for comparisons among multiple groups, followed by Dunn’s multiple comparisons. All data are presented as mean ± s.e.m. The sample sizes and statistical tests used for each experiment are stated in the figures or figure legends.

Acknowledgements

This research was funded by the National Key R&D Program of China to SFW (2022YFD1700200) (https://service.most.gov.cn/) and National Natural Science Foundation of China to SFW (32022011) (https://www.nsfc.gov.cn/).

Funding: This work was supported by the National Natural Science Foundation of China to SFW (No.32022011 & 31772205) (https://isisn.nsfc.gov.cn/egrantindex/funcindex/prjsearch-list?locale=zh_CN).

The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.