Reference Gene Selection and Evaluation for Gene Expression Studies Using qRT-PCR in the White-Backed Planthopper, Sogatella furcifera (Hemiptera: Delphacidae)

Molecular characterization of Vitellogenin-like1 gene in Sogatella furcifera (Hemiptera: Delphacidae), and its function on reproduction

In this study, a vitellogenin like1 gene (SfVg-like1) in Sogatella furcifera was identified. The open reading frame (ORF) encoded 1,321 amino acid sequence. Structure analysis reveals that the amino acid sequence of SfVg-like1 has 3 conserved LPD_N, DUF1943 and VWFD domains. Phylogenetic analyses showed that SfVg-like1 was clustered in the same branch with the Vg-like1 of Nilaparvata lugens (100% bootstrap value) compared with other Hemiptera insects Vgs associated with vitellogenesis. Temporo-spatial expression analyses showed that SfVg-like1 expressed during all stages, and in both genders. The relative expression levels of SfVg-like1 mRNA were higher in adults than in nymph developmental stages. The knockdown of SfVg-like1 gene resulted in the inhibition of the ovarian development in female adults, whereas the morphology of the testis in male adults was not been affected. The silence of SfVg-like1 could decrease the relative expression levels of target of rapamycin (SfTOR, GenBank MW193765) and vitellogenin (SfVg, GenBank MH271114) genes significantly in female adults. However, the knockdown of SfTOR or SfVg genes in female adults did not affect the transcript level of SfVg-like1. Therefore, it demonstrated that SfVg-like1 might locate on the upstream signaling pathways of SfTOR and SfVg. These results demonstrate that SfVg-like1 is essential for S. furcifera reproduction, and it could be the potential target for the control of this pest.

Selection and validation of optimal reference genes for RT-qPCR analyses in Aphidoletes aphidimyza Rondani (Diptera: Cecidomyiidae)

Aphidoletes aphidimyza is a predator that is an important biological agent used to control agricultural and forestry aphids. Although many studies have investigated its biological and ecological characteristics, few molecular studies have been reported. The current study was performed to identify suitable reference genes to facilitate future gene expression and function analyses via quantitative reverse transcription PCR. Eight reference genes glyceraldehyde-3-phosphate dehydrogenase (GAPDH), RPS13, RPL8, RPS3, α-Tub, β-actin, RPL32, and elongation factor 1 alpha (EF1-α) were selected. Their expression levels were determined under four different experimental conditions (developmental stages, adult tissues, sugar treatment, and starvation treatment) using qRT-PCR technology. The stability was evaluated with five methods (Ct value, geNorm, NormFinder, BestKeeper, and RefFinder). The results showed that GAPDH, RPL32, and EF1-α were ranked as the best reference gene combinations for measuring gene expression levels among different developing stages and in various starvation treatments. RPL8 and RPS3 were recommended to normalize the gene expression levels among different adult tissues. RPL32, β-actin, and EF1-α were recommended sugar-feeding conditions. To validate the utility of the selected reference pair, RPL8, and RPS3, we estimated the tissue-biased expression level of a chemosensory protein gene (AaphCSP1). As expected, AaphCSP1 is highly expressed in the antennae and lowly expressed in the abdomen. These findings will lay the foundation for future research on the molecular physiology and biochemistry of A. aphidimyza.

Evaluation of suitable reference genes for expression analysis using quantitative real-time polymerase chain reaction in the parasitoid Exorista sorbillans (Diptera: Tachinidae)

The parasitoid Exorista sorbillans (Diptera: Tachinidae) is a larval endoparasitoid of the silkworm Bombyx mori, causing severe damage to silkworm cocoon industry. It is also an important natural enemy resource of insect pests in agriculture and forestry. Despite their roles in biocontrol and pest status on sericulture, there has been limited research on the functional studies of dipteran parasitoids. Quantitative real-time polymerase chain reaction (qRT-PCR) is the most commonly used to address gene functions. Using qRT-PCR, stably expressed reference genes under different experimental conditions are required to normalize the expression of target genes. However, no information regarding suitable qRT-PCR reference genes in dipteran parasitoids has been reported. In this study, we evaluate the expression stability of nine commonly used reference genes in insects including eukaryotic translation elongation factor 1δ (eEF1δ), elongation factor 2, 18S ribosomal RNA (18S rRNA), tubulin 3, actin87, ribosomal protein 49 (RP49), ribosomal protein S15, glyceraldehyde-3-phosphate dehydrogenase, and TATA-box binding protein (TBP) in E. sorbillans under different treatments, including tissues, developmental stages, genders, feeding density and pesticide stress, using ∆C_t , BestKeeper, geNorm, Normfinder and RefFinder, respectively. The results showed that the genes RP49, eEF1δ and 18S rRNA were recommended as the most suitable reference genes in E. sorbillans across all experimental conditions. This finding provides the necessary foundation for future functional studies in E. sorbillans and its effective use in both sericulture and pest control.

Mucin Protein of Aedes aegypti Interacts with Dengue Virus 2 and Influences Viral Infection

Dengue, caused by dengue virus (DENV), is the most prevalent vector-borne viral disease, posing a serious health concern to 2.5 billion people worldwide. DENV is primarily transmitted among humans by its mosquito vector Aedes aegypti; hence, the identification of a novel dengue virus receptor in mosquitoes is critical for the development of new anti-mosquito measures. In the current study, we have identified peptides which potentially interact with the surface of the virion particles and facilitate virus infection and movement during their life cycle in the mosquito vector. To identify these candidate proteins, we performed phage-display library screening against domain III of the envelope protein (EDIII), which plays an essential role during host cell receptor binding for viral entry. The mucin protein, which shared sequence similarity with the peptide identified in the screening, was cloned, expressed, and purified for in vitro interaction studies. Using in vitro pulldown and virus overlay protein-binding assay (VOPBA), we confirmed the positive interaction of mucin with purified EDIII and whole virion particles. Finally, blocking of mucin protein with anti-mucin antibodies partially reduced DENV titers in infected mosquitos. Moreover, mucin protein was found to be localized in the midgut of Ae. aegypti. IMPORTANCE Identification of interacting protein partners of DENV in the insect vector Aedes aegypti is crucial for designing vector control-based strategies and for understanding the molecular mechanism DENV uses to modulate the host, gain entry, and survive successfully. Similar proteins can be used in generating transmission-blocking vaccines.

Reproductive Outbreaks of Sogatella furcifera Mediated by Overexpression of the Nuclear Receptor USP under Pressure from Triflumezopyrim

Long-term pesticide-driven selection pressure is one of the main causes of insect outbreaks. In this study, we found that low doses of triflumezopyrim could increase the fecundity of white-backed planthoppers (Sogatella furcifera). By continuously screening 20 generations with a low dose of triflumezopyrim, a triflumezopyrim-resistant strain (Tri-strain, resistance ratio = 20.9-fold) was obtained. The average oviposition quantity and longevity of the Tri-strain (208.77 eggs and 21.31 days, respectively) were significantly higher than those of the susceptible strain (Sus-strain) (164.62 eggs and 17.85 days, respectively). To better understand the mechanism underlying the effects on reproduction, we detected the expression levels of several reproduction-related transcription factors in both the Tri- and Sus-strains. Ultraspiracle (USP) was significantly overexpressed in the Tri-strain. Knockdown of USP by RNAi severely inhibited the moulting process of S. furcifera and disrupted the development of female adult ovaries. Among the potential downstream target genes of USP, Kr-h1 (0.19-fold), Cht8 (0.56-fold) and GPCR A22 (0.31-fold) showed downregulated expression after USP-RNAi. In contrast, the expression of EcR (2.55-fold), which forms heterodimers with USP, was significantly upregulated. Furthermore, RNAi was performed on Kr-h1 in the Tri-strain, and the results show that larval moulting and the development of female adult ovaries were inhibited, consistent with the USP-RNAi results in S. furcifera. These results suggest that the transcription factors USP and Kr-h1 play important roles in the reproductive development of S. furcifera, and overexpression of USP and Kr-h1 in the Tri-resistant strain may result in reproductive outbreaks of pests.

SfDicer2 RNA Interference Inhibits Molting and Wing Expansion in Sogatella furcifera

Simple Summary

Endoribonuclease 2 (Dicer2) plays various physiological roles in the RNA interference (RNAi) pathway by fragmenting double-stranded RNA to generate small interfering RNA, which then mediates gene silencing. In this study, the role of Dicer2 in the regulation of molting and wing expansion in Sogatella furcifera (white-backed planthopper) was investigated. In particular, SfDicer2-mediated RNAi resulted in wing deformities and lethal modifications in S. furcifera, which are attributable to the significant inhibition of chitin synthesis and degradation and wing expansion genes. This study provides insights into the biological functions of Dicer2 in insects, which can aid in RNAi-mediated pest control.

Abstract

Endoribonuclease 2 (Dicer2) is a key nicking endonuclease involved in the small interfering RNA biosynthesis, and it plays important roles in gene regulation and antiviral immunity. The Dicer2 sequence was obtained using the transcriptomic and genomic information of Sogatella furcifera (Horváth), and the spatiotemporal characteristics and functions of molting and wing expansion regulation were studied using real-time quantitative polymerase chain reaction and RNA interference (RNAi) technology. The expression of SfDicer2 fluctuated during the nymphal stage of S. furcifera. Its expression decreased significantly over the course of molting. SfDicer2 exhibited the highest transcript level in the nymphal stage and adult fat body. After SfDicer2 was silenced, the total mortality rate was 42.69%; 18.32% of the insects died because of their inability to molt. Compared with the effects of dsGFP or water, 44.38% of the insects subjected to the silencing of SfDicer2 exhibited wing deformities after successful eclosion. After SfDicer2 RNAi, the expression of chitinase, chitin deacetylase, trehalase, chitin synthase 1, and wing expansion-related genes was significantly inhibited. These findings indicate that SfDicer2 controls molting by affecting genes associated with chitin synthesis and degradation and regulates wing expansion by altering the expression of wing expansion-related genes in S. furcifera.

Silencing the Autophagy-Related Genes ATG3 and ATG9 Promotes SRBSDV Propagation and Transmission in Sogatella furcifera

Autophagy plays diverse roles in the interaction among pathogen, vector, and host. In the plant virus and insect vector system, autophagy can be an antiviral/pro-viral factor to suppress/promote virus propagation and transmission. Here, we report the antiviral role of autophagy-related genes ATG3 and ATG9 in the white-backed planthopper (Sogatella furcifera) during the process of transmitting the southern rice black-streaked dwarf virus (SRBSDV). In this study, we annotated two autophagy-related genes, SfATG3 and SfATG9, from the female S. furcifera transcriptome. The cDNA of SfATG3 and SfATG9 comprised an open reading frame (ORF) of 999 bp and 2295 bp that encodes a protein of 332 and 764 amino acid residues, respectively. SfATG3 has two conserved domains and SfATG9 has one conserved domain. In S. furcifera females exposed to SRBSDV, expression of autophagy-related genes was significantly activated and shared similar temporal patterns to those of SRBSDV S9-1 and S10, all peaking at 4 d post viral exposure. Silencing the expression of SfATG3 and SfATG9 promoted SRBSDV propagation and transmission. This study provides evidence for the first time that S. furcifera autophagy-related genes ATG3 and ATG9 play an antiviral role to suppress SRBSDV propagation and transmission.

Endosymbionts Reduce Microbiome Diversity and Modify Host Metabolism and Fecundity in the Planthopper Sogatella furcifera

Endosymbionts can strongly affect bacterial microbiota in pests. The white-backed planthopper Sogatella furcifera, a notorious pest in rice, is usually co-infected with Cardinium and Wolbachia, but the effects of these endosymbionts together or individually on the host microbiome and fecundity are unclear. Here, we established three S. furcifera lines (Cardinium and Wolbachia double-infected, Cardinium single-infected, and both-uninfected lines) backcrossed to a common nuclear background and found that single and double infections reduced bacterial diversity and changed bacterial community structure across nymph and adult stages and across adult tissues. The endosymbionts differed in densities between adults and nymphs as well as across adult tissues, with the distribution of Cardinium affected by Wolbachia. Both the single infection and particularly the double infection reduced host fecundity. Lines also differed in levels of metabolites, some of which may influence fecundity (e.g., arginine biosynthesis and nicotinamide metabolism). Cardinium in the single-infected line upregulated metabolic levels, while Wolbachia in the double-infected line appeared to mainly downregulate them. Association analysis pointed to possible connections between various bacteria and differential metabolites. These results reveal that Cardinium by itself and in combination with Wolbachia affect bacterial microbiota and levels of metabolites, with likely effects on host fecundity. Many of the effects of these metabolically limited endosymbionts that are dependent on the hosts may be exerted through manipulation of the microbiome.

IMPORTANCE Endosymbionts can profoundly affect the nutrition, immunity, development, and reproduction of insect hosts, but the effects of multiple endosymbiont infections on microbiota and the interaction of these effects with insect host fitness are not well known. By establishing S. furcifera lines with different endosymbiont infection status, we found that Cardinium and the combined Cardinium + Wolbachia infections differentially reduced bacterial diversity as well as changing bacterial community structure and affecting metabolism, which may connect to negative fitness effects of the endosymbionts on their host. These results established the connections between reduced bacterial diversity, decreased fecundity and metabolic responses in S. furcifera.

Identifying optimal reference genes for gene expression studies in Eurasian spruce bark beetle, Ips typographus (Coleoptera: Curculionidae: Scolytinae)

Eurasian spruce bark beetle (Ips typographus [L.]) causes substantial damage to spruce forests worldwide. Undoubtedly, more aggressive measures are necessary to restrict the enduring loss. Finishing genome sequencing is a landmark achievement for deploying molecular techniques (i.e., RNA interference) to manage this pest. Gene expression studies assist in understanding insect physiology and deployment of molecular approaches for pest management. RT-qPCR is a valuable technique for such studies. However, accuracy and reliability depend on suitable reference genes. With the genome sequence available and the growing requirement of molecular tools for aggressive forest pest management, it is crucial to find suitable reference genes in Ips typographus under different experimental conditions. Hence, we evaluated the stability of twelve candidate reference genes under diverse experimental conditions such as biotic (developmental, sex and tissues) and abiotic factors (i.e., temperature and juvenile hormone treatment) to identify the reference genes. Our results revealed that ribosomal protein 3a (RPS3-a) was the best reference gene across all the experimental conditions, with minor exceptions. However, the stability of the reference gene can differ based on experiments. Nevertheless, present study provides a comprehensive list of reference genes under different experimental conditions for Ips typographus and contributes to "future genomic and functional genomic research".

Selection and Validation of Reference Genes for Quantitative Real-Time PCR Normalization in Athetis dissimilis (Lepidoptera: Noctuidae) Under Different Conditions

Reference genes are the key to study gene expression patterns using quantitative real-time PCR (qRT-PCR). No studies on the reference genes of Athetis dissimilis, an important agricultural pest, have been reported. In order to determine the reference genes for qRT-PCR normalization in A. dissimilis under different conditions, 10 candidate genes [18S ribosomal protein (18S), 28S ribosomal protein (28S), arginine kinase (AK), elongation factor 1 alpha (EF1-α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L32 (RPL32), ribosomal protein L40 (RPL40), alpha-tubulin (α-TUB), beta-actin (β-ACT), and beta-tubulin (β-TUB)] of A. dissimilis were selected to evaluate their stability as reference genes under different biotic and abiotic conditions by using five tools, geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder. Furthermore, CSP1 and superoxide dismutase (SOD) were used as target genes to validate the candidate reference genes. The results showed that different reference genes were needed under different experimental conditions, among which, EF-1α, RPL40, and 18S are most suitable reference genes for studying genes related development stages of A. dissimilis, RPL40 and α-TUB for larval tissues, α-TUB and 28S for adult tissues, EF-1α and β-ACT for insecticidal treatments, β-ACT and 28S for temperature treatments, EF-1α and β-ACT for starvation treatments, RPL40 and 18S for dietary treatments, and 18S, 28S, and α-TUB for all the samples. These results provide suitable reference genes for studying gene expression in A. dissimilis under different experimental conditions, and also lay the foundation for further research into the function of related genes in A. dissimilis.

Comparison of Transcriptome Responses between Sogatella furcifera Females That Acquired Southern Rice Black-Streaked Dwarf Virus and Not

Simple Summary

The southern rice black-streaked dwarf virus (SRBSDV) is transmitted horizontally by the planthopper, Sogatella furcifera. During feeding on virus-infected plants, S. furcifera may acquire or fail to acquire SRBSDV. In this study, the responses were compared among the S. furcifera successfully acquiring the virus, those failing to acquire the virus, and those not exposed to SRBSDV (the control). A total of 1043 and 2932 differentially expressed genes (DEGs) were obtained in S. furcifera females that acquired SRBSDV and that failed to, in comparison with the control, respectively. Functionally, these DEGs are primarily involved in diverse signaling pathways related to primary metabolism and innate immunity, such as apoptosis. Additional bioassays confirmed the activation of apoptosis in S. furcifera by SRBSDV exposure. Interestingly, we also found that six female-specific genes were also upregulated in S. furcifera females exposed to SRBSDV. Our results further the understanding of the interactions between the vector S. furcifera females and SRBSDV at the molecular level.

Abstract

The southern rice black-streaked dwarf virus (SRBSDV) is transmitted horizontally by Sogatella furcifera in a persistent, propagative manner. Exposure of S. furcifera females to SRBSDV-infected rice plants may trigger transcriptomic changes in the insects, the transcriptomes of females that acquired SRBSDV and those that failed to, as well as females fed on healthy rice plants as control, were sequenced and compared. Nine transcriptomic libraries were constructed, from which a total of 53,084 genes were assembled. Among the genes, 1043 and 2932 were differentially expressed genes (DEGs) in S. furcifera females that acquired SRBSDV and that failed to, in comparison with the control, respectively. Functional enrichment analysis showed that DEGs identified in S. furcifera females exposed to SRBSDV are primarily involved in diverse signaling pathways related to primary metabolism and innate immunity. The DEGs in the S. furcifera females that failed to acquire the virus significantly outnumbered that in the insects that acquired the virus, and the virus exposure activated the humoral and cellular immune responses of the vectors, especially the apoptosis. The key gene in apoptosis encoding caspase 1 was upregulated by SRBSDV exposure, especially in S. furcifera females that failed to acquire the virus. Analysis of caspase 1 activity validated that SRBSDV exposure induced caspase 1 accumulation. Surprisingly, the expression of six female-specific genes was also upregulated by SRBSDV exposure, which was confirmed by RT-qPCR analysis. This study provides evidence to explain the differential virus acquisition at the transcriptome level.

Selection and Validation of Reference Genes for RT-qPCR Normalization in Bradysia odoriphaga (Diptera: Sciaridae) Under Insecticides Stress

Bradysia odoriphaga (Diptera: Sciaridae) is the most serious root maggot pest which causes substantial damage to the Chinese chive. Organophosphate (OP) and neonicotinoid insecticides are widely used chemical pesticides and play important roles in controlling B. odoriphaga. However, a strong selection pressure following repeated pesticide applications has led to the development of resistant populations of this insect. To understand the insecticide resistance mechanism in B. odoriphaga, gene expression analysis might be required. Appropriate reference gene selection is a critical prerequisite for gene expression studies, as the expression stability of reference genes can be affected by experimental conditions, resulting in biased or erroneous results. The present study shows the expression profile of nine commonly used reference genes [elongation factor 1α (EF-1α), actin2 (ACT), elongation factor 2α (EF-2α), glucose-6-phosphate dehydrogenase (G6PDH), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L10 (RPL10), ribosomal protein S3 (RPS3), ubiquitin-conjugating enzyme (UBC), and α-tubulin (TUB)] was systematically analyzed under insecticide stress. Moreover, we also evaluated their expression stability in other experimental conditions, including developmental stages, sexes, and tissues. Five programs (NormFinder, geNorm, BestKeeper, RefFinder, and ΔCt) were used to validate the suitability of candidate reference genes. The results revealed that the most appropriate sets of reference genes were RPL10 and ACT across phoxim; ACT and TUB across chlorpyrifos and chlorfluazuron; EF1α and TUB across imidacloprid; EF1α and EF2α across developmental stages; RPL10 and TUB across larvae; EF1α and ACT across tissues, and ACT and G6PDH across sex. These results will facilitate the standardization of RT-qPCR and contribute to further research on B. odoriphaga gene function under insecticides stress.

Expression Analyses of Vitellogenin and Target of Rapamycin of Sogatella furcifera (Hemiptera: Delphacidae), and Their Effects on Reproduction

Vitellogenin (Vg) and the target of rapamycin (TOR) are important genes involved in insect reproduction regulation. In this work, the full lengths of the Vg (SfVg) and TOR (SfTOR) genes of the white-backed planthopper Sogatella furcifera were cloned. The expression pattern in females showed that SfVg was highly expressed in fat bodies, and SfTOR was highly expressed in Malpighian tubules. After silencing SfVg or SfTOR, female adults did not deposit eggs. Their ovarian development was delayed, and yolk protein deposition in the oocytes was reduced. However, wild-type females mated with SfTOR-silenced males could lay eggs and produce offsprings normally. The dissections of testes and accessory glands of males with SfTOR knockdown showed that their development was not affected. Therefore, the silencing of the SfVg or SfTOR genes can effectively inhibit female reproduction, but SfTOR knockdown has no significant effect on male reproductive capacity. Furthermore, silencing SfTOR can cause SfVg expression to decrease significantly. All of the above results revealed that SfVg and SfTOR are essential for white-backed planthopper reproduction and may provide a potential target for pest control.

Reference Gene Selection for Normalizing Gene Expression in Ips Sexdentatus (Coleoptera: Curculionidae: Scolytinae) Under Different Experimental Conditions

Ips sexdentatus (Coleoptera: Curculionidae: Scolytinae) is one of the most destructive and economically important forest pests. A better understanding of molecular mechanisms underlying its adaptation to toxic host compounds may unleash the potential for future management of this pest. Gene expression studies could be considered as one of the key experimental approaches for such purposes. A suitable reference gene selection is fundamental for quantitative gene expression analysis and functional genomics studies in I. sexdentatus. Twelve commonly used reference genes in Coleopterans were screened under different experimental conditions to obtain accurate and reliable normalization of gene expression data. The majority of the 12 reference genes showed a relatively stable expression pattern among developmental stages, tissue-specific, and sex-specific stages; however, some variabilities were observed during varied temperature incubation. Under developmental conditions, the Tubulin beta-1 chain (β-Tubulin) was the most stable reference gene, followed by translation elongation factor (eEF2) and ribosomal protein S3 (RPS3). In sex-specific conditions, RPS3, β-Tubulin, and eEF2 were the most stable reference genes. In contrast, different sets of genes were shown higher stability in terms of expression under tissue-specific conditions, i.e., RPS3 and eEF2 in head tissue, V-ATPase-A and eEF2 in the fat body, V-ATPase-A and eEF2 in the gut. Under varied temperatures, β-Tubulin and V-ATPase-A were most stable, whereas ubiquitin (UbiQ) and V-ATPase-A displayed the highest expression stability after Juvenile Hormone III treatment. The findings were validated further using real-time quantitative reverse transcription PCR (RT-qPCR)-based target gene expression analysis. Nevertheless, the present study delivers a catalog of reference genes under varied experimental conditions for the coleopteran forest pest I. sexdentatus and paves the way for future gene expression and functional genomic studies on this species.

Cloning and Functional Verification of CYP408A3 and CYP6CS3 Related to Chlorpyrifos Resistance in the Sogatella furcifera (Horváth) (Hemiptera: Delphacidae)

The white-back planthopper (WBPH), Sogatella furcifera, mainly harms rice and occurs in most rice regions in China and Asia. With the use of chemical pesticides, S. furcifera has developed varying degrees of resistance to a variety of pesticides. In our study, a chlorpyrifos-resistant population (44.25-fold) was built through six generations of screening with a sublethal dose of chlorpyrifos (LD50) from a field population. The expression levels of ten selected resistance-related P450 genes were analyzed by RT-qPCR and found that CYP408A3 and CYP6CS3 were significantly more expressed in the third instar nymphs of the XY17-G5 and XY17-G6 populations, about 25-fold more than the Sus-Lab strain, respectively (p < 0.01). To elucidate their molecular function in the development of resistance towards chlorpyrifos, we cloned two P450 full lengths and predicted their tertiary protein structures. CYP408A3 and CYP6CS3 were also downregulated after injecting dsCYP408A3, dsCYP6CS3, or their mixture compared to the control group. Moreover, the mortality rates of the dsCYP6CS3 (91.7%) and the mixture injection treatment (93.3%) treated by the LC50 concentration of chlorpyrifos were significantly higher than the blank control group (51.7%) and dsCYP408A3 injection treatment (69.3%) at 72 h (p < 0.01). Meanwhile, the P450 enzyme activities in the dsRNA treatments were lower than that in the control (XY17-G6) (p < 0.01). Therefore, the P450 gene CYP6CS3 may be one of the main genes in the development of chlorpyrifos resistance in S. furcifera.

Selection and Validation of Reference Genes for Gene Expression Analysis in Tuta absoluta Meyrick (Lepidoptera: Gelechiidae)

Simple Summary

Reference genes are critical for standardizing expression data of RT-qPCR across samples of organisms under different experimental conditions. However, most commonly used reference genes may not be stably expressed leading to a risk of misinterpretation of results. In our study, nine reference genes were evaluated in Tuta absoluta (a destructive pest of tomato) at different developmental stages, tissues, 20-hydroxyecdysone (20E) and insecticide treatments. Finally, the expression profile of indicator gene EcR after 20E treatment was evaluated to verify the accuracy of the results. This study is essential for improving accuracy and reliability to normalize gene expression data in T. absoluta and provides a useful strategy for other insects.

Abstract

The tomato leaf miner, Tuta absoluta is a destructive pest of tomato. The leaf-mining activities of its larvae can cause significant yield losses. Real-time quantitative polymerase chain reaction (RT-qPCR) is commonly used to measure gene expression, and the selection of stable reference genes for calibration and standardization is critical for accurate use of RT-qPCR. We studied the stable expression of nine common housekeeping genes in T. absoluta. These were examined at different developmental stages, in larval tissues, as well as those induced by exposure to 20E and insecticides. Four dedicated algorithms (geNorm, BestKeeper, NormFinder, and ΔCt method) and online tool (RefFinder) were used to analyze and rank the tested reference genes. Based on the standardized gene expression data of target gene ecdysone receptor (EcR), the applicability of specific reference genes was verified. The results clarify that the optimal internal reference genes vary greatly under different experimental conditions. GAPDH and RPS11 were the best reference genes for developmental stages; RPL28 and RPL10 for different tissues; EF1α and RPL28 for 20E treatment; EF1α and RPL7A for insecticide treatments. The most suitable reference genes in all experimental conditions are EF1α and RPL28.

Selection and Validation of Reference Genes For qRT-PCR Analysis of Rhopalosiphum padi (Hemiptera: Aphididae)

Rhopalosiphum padi (L.) (Hemiptera: Aphididae) is an important cosmopolitan pest in cereal crops. Reference genes can significantly affect qRT-PCR results. Therefore, selecting appropriate reference genes is a key prerequisite for qRT-PCR analyses. This study was conducted to identify suitable qRT-PCR reference genes in R. padi. We systematically analyzed the expression profiles of 11 commonly used reference genes. The ΔCt method, the BestKeeper, NormFinder, geNorm algorithms, and the RefFinder online tool were used to evaluate the suitability of these genes under diverse experimental conditions. The data indicated that the most appropriate sets of reference genes were β-actin and GAPDH (for developmental stages), AK and TATA (for populations), RPS18 and RPL13 (for tissues), TATA and GAPDH (for wing dimorphism), EF-1α and RPS6 (for antibiotic treatments), GAPDH and β-actin (for insecticide treatments), GAPDH, TATA, RPS18 (for starvation-induced stress), TATA, RPS6, and AK (for temperatures), and TATA and GAPDH (for all conditions). Our study findings, which revealed the reference genes suitable for various experimental conditions, will facilitate the standardization of qRT-PCR programs, while also improving the accuracy of qRT-PCR analyses, with implications for future research on R. padi gene functions.

Status of insecticide resistance and biochemical characterization of chlorpyrifos resistance in Sogatella furcifera (Hemiptera:Delphacidae) in Sichuan Province, China

The white-backed planthopper, Sogatella furcifera (Horváth) (Hemiptera, Delphacidae), is an energetic rice insect pest in rice production or rice-growing areas. Due to excessive use of the chemical insecticide, S. furcifera has produced the high resistance to some frequently used insecticides. In this paper, the resistance levels of S. furcifera from the eight different areas of Sichuan Province against the five chemicals were monitored by using the rice seedling dipping during 2017-2018 to understand the resistance levels. The results showed that most of all populations have developed low or moderate level of resistance for chlorpyrifos (3.4 to 44.3-fold) and thiamethoxam (3.9- to 15.5-fold), the populations in the LS (1.7 to 5.4- fold)and WS (1.6 to 5.0- fold) regions were still sensitive or low resistance levels compared with other local populations. Almost all populations displayed the susceptible to imidacloprid (0.9- to 5.0-fold), buprofezin (0.9- to 4.3-fold) or low levels of resistance to pymetrozine (1.5- to 6.8-fold). The synergism experiment indicated that P450 enzymes may be important contributed to the metabolic detoxification of chlorpyrifos. The cross-resistance bioassay showed that there was no cross-resistance between chlorpyrifos and triflumezopyrim, but for sulfoxaflor, in the XY17 population. The relative expression level of twelve insecticide resistant-related P450 genes were analyzed by using qRT-PCR and found that CYP4C77, CYP418A1, CYP418A2, CYP408A3 and CYP6ER4 were significantly more expressed in the 3rd-instar nymph of the XY17 and XY18 field populations. To determine the main resistant-related P450 gene for chlorpyrifos, the relative expression level of five P450 genes were detected by using qRT-PCR from the G2 and G4 generation of XY17 under the pressure with LC₅₀ of chlorpyrifos. The results showed that CYP6ER4 was significantly up-regulated expression in XY17 G2 and G4 generations population over 700-fold (P < 0.01). The full length and proteins tertiary structure were also cloned and predicted. Meanwhile, the function of CYP6ER4 was analyzed by RNA interference and the results indicated that the relative expression of CYP6ER4 in the XY17 (G4) population after injected dsRNA was lower than that in the dsGFP injected group. Moreover, the mortality rates of the S. furcifera treated with the LC₅₀ concentration of chlorpyrifos after dsRNA microinjection was significantly higher than that of the dsGFP injected group 72 h after treatment (P < 0.01). Therefore, the overexpression of CYP6ER4 may be one of the primary factors in the development of chlorpyrifos resistance in S. furcifera.

Selection and validation of experimental condition-specific reference genes for qRT-PCR in Metopolophium dirhodum (Walker) (Hemiptera: Aphididae)

Metopolophium dirhodum (Walker) (Hemiptera: Aphididae) is one of the most common aphid pests of winter cereals. To facilitate accurate gene expression analyses with qRT-PCR assays, the expression stability of candidate reference genes under specific experimental conditions must be verified before they can be used to normalize target gene expression levels. In this study, 10 candidate reference genes in M. dirhodum were analyzed by qRT-PCR under various experimental conditions. Their expression stability was evaluated with delta Ct, BestKeeper, geNorm, and NormFinder methods, and the final stability ranking was determined with RefFinder. The results indicate that the most appropriate sets of internal controls were SDHB and RPL8 across geographic population; RPL8, Actin, and GAPDH across developmental stage; SDHB and NADH across body part; RPL8 and Actin across wing dimorphism and temperature; RPL4 and EF1A across starvation stress; AK and RPL4 across insecticide treatments; RPL8 and NADH across antibiotic treatments; RPL8, RPL4, Actin, and NADH across all samples. The results of this study provide useful insights for establishing a standardized qRT-PCR procedure for M. dirhodum and may be relevant for identifying appropriate reference genes for molecular analyses of related insects.

Identification of Ginger (Zingiber officinale Roscoe) Reference Genes for Gene Expression Analysis

Quantitative real-time PCR (qRT-PCR) is widely used in the detection of gene expression level. However, there is no suitable ginger reference gene for qPCR analysis. Therefore, it is the primary task to select and validate the appropriate ginger reference gene to normalize the expression of target genes. In this study, 14 candidate reference genes were selected and analyzed in different tissues (leaf, and rhizome), different development stages, different varieties, and abiotic stress (ABA and salt stress). Expression stability was calculated using geNorm and NormFinder, Bestkeeper, and RefFinder. For abiotic stress and total conditions, 28S and COX were identified as the most stable genes. In addition, RPII was the most stable in the different development stages and different varieties. TEF2 and RPL2 were the least stably expressed in the tissue and all the conditions. In order to verify the feasibility of these genes as reference genes, we used the most stable and least stable reference genes to normalize the expression levels of ZoSPS genes under different conditions. This work can provide theoretical support for future research on ginger gene expression.

A novel negative thermotaxis behavior in rice planthoppers is regulated by TRPA1 channel

BACKGROUND

In the field, we observed that many white-backed planthoppers (Sogatella furcifera, WBPHs) stayed in the top region of rice plants exposed to direct sunshine. It was known that WBPHs frequently took flight when the ground temperature was about 25 °C, then climbed to and flew in a dense layer corresponding to an air temperature of about 16 °C in the sky. Its migration height was higher than the top of the surface temperature inversion. It is still unclear whether WBPHs prefer warm or cold regions, and therefore we studied the thermal responses of WBPHs and other insects using a simulated system.

RESULTS

We found that WBPHs preferred a cold region to a warm one, unexpectedly below their comfort temperature zone. After comparative analysis with other insect species, such as small brown planthoppers, brown planthoppers, Trialeurodes vaporariorum (stinkbugs, a predator of planthoppers) and Bemisia tabaci (whitefly), only three planthoppers showed cold preference behavior. RNA interference experiments revealed that this behavior of WBPHs can be regulated by the transient receptor potential (TRP) channels TRPA1 channel. Furthermore, podocarpic acid, an agonist of TRPA1, weakened the cold preference, whereas A-967079, an antagonist of TRPA1, had the opposite effect.

CONCLUSION

We reported a novel cold preference (negative thermotaxis) in rice planthoppers that was regulated in WBPHs by the TRPA1 channel. Cold preference of rice planthoppers is probably related to its choice behavior of the special migratory temperature layer. Our results expanded a new perspective to develop novel strategies for behavioral manipulation and management of rice planthoppers. © 2020 Society of Chemical Industry.

Full-length transcriptome analysis of Phytolacca americana and its congener P. icosandra and gene expression normalization in three Phytolaccaceae species

BACKGROUND

Phytolaccaceae species in China are not only ornamental plants but also perennial herbs that are closely related to human health. However, both large-scale full-length cDNA sequencing and reference gene validation of Phytolaccaceae members are still lacking. Therefore, single-molecule real-time sequencing technology was employed to generate full-length transcriptome in invasive Phytolacca americana and non-invasive exotic P. icosandra. Based on the transcriptome data, RT-qPCR was employed to evaluate the gene expression stability in the two plant species and another indigenous congener P. acinosa.

RESULTS

Total of 19.96 Gb and 19.75 Gb clean reads of P. americana and P. icosandra were generated, including 200,857 and 208,865 full length non-chimeric (FLNC) reads, respectively. Transcript clustering analysis of FLNC reads identified 89,082 and 98,448 consensus isoforms, including 86,989 and 96,764 high-quality ones. After removing redundant reads, 46,369 and 50,220 transcripts were obtained. Based on structure analysis, total 1675 and 1908 alternative splicing variants, 25,641 and 31,800 simple sequence repeats (SSR) as well as 34,971 and 36,841 complete coding sequences were detected separately. Furthermore, 3574 and 3833 lncRNA were predicted and 41,676 and 45,050 transcripts were annotated respectively. Subsequently, seven reference genes in the two plant species and a native species P. acinosa were selected and evaluated by RT-qPCR for gene expression analysis. When tested in different tissues (leaves, stems, roots and flowers), 18S rRNA showed the highest stability in P. americana, whether infested by Spodoptera litura or not. EF2 had the most stable expression in P. icosandra, while EF1-α was the most appropriate one when attacked by S. litura. EF1-α showed the highest stability in P.acinosa, whereas GAPDH was recommended when infested by S. litura. Moreover, EF1-α was the most stable one among the three plant species whenever germinating seeds or flowers only were considered.

CONCLUSION

Full-length transcriptome of P. americana and P. icosandra were produced individually. Based on the transcriptome data, the expression stability of seven candidate reference genes under different experimental conditions was evaluated. These results would facilitate further exploration of functional and comparative genomic studies in Phytolaccaceae and provide insights into invasion success of P. americana.

Suitable reference genes for RT-qPCR analysis in Dichelops melacanthus (Hemiptera: Pentatomidae)

The relative quantification of gene expression is mainly realized through reverse transcription-quantitative PCR (RT-qPCR). However, the accuracy of this technique is deeply influenced by the expression stability of the reference genes used for data normalization. Therefore, the selection of suitable reference genes for a given experimental condition is a prerequisite in gene expression studies. Dichelops melacanthus (Hemiptera: Pentatomidae) is an important phloem sap-sucking insect pest of soybean, wheat, and maize in Brazil. Most of the genetic and molecular biology studies require gene expression analysis. Nevertheless, there are no reports about reference genes for RT-qPCR data normalization in D. melacanthus. In this study, we evaluated the expression stability of nine candidate reference genes (nadh, sdhb, gapdh, fau, ef1a, rpl9, ube4a, gus and rps23) in different developmental stages, body parts, sex, starvation-induced stress and dsRNA exposure by RefFinder software that integrates the statistical algorithms geNorm, NormFinder, BestKeeper, and ΔCt method. Our results showed that ef1a and nadh are the most stable reference genes for developmental stages, fau and rps23 for sex, ube4a and rps23 for body parts, rpl9 and fau for starvation stress, and nadh and sdhb for dsRNA exposure treatment. The reference genes selected in this work will be useful for further RT-qPCR analyses on D. melacanthus, facilitating future gene expression studies that can provide a better understanding of the developmental, physiological, and molecular processes of this important insect pest. Moreover, the knowledge gained from these studies can be helpful to design effective and sustainable pest management strategies.

Selection of Reference Genes for RT-qPCR Analysis Under Extrinsic Conditions in the Hawthorn Spider Mite, Amphitetranychus viennensis

Hawthorn spider mite, Amphitetranychus viennensis Zacher, is an economically important arthropod pest for fruit trees and woody ornamental plants. Extensive and repetitive use of synthetic acaricides has led to the development of resistance in A. viennensis. To understand the molecular basis of pesticide resistance, and to develop genetic-based control alternatives (e.g., RNAi-based biopesticides), a standardized protocol for real-time quantitative reverse transcription PCR (RT-qPCR) is needed. In the proceeding phase of this research, we screened for the internal references for RT-qPCR analysis from a pool of A. viennensis housekeeping genes under the intrinsic conditions, including developmental stage, sex, and diapause. Here, we continued our efforts to search for the reference genes under an array of extrinsic conditions, including temperature, humidity, photoperiod, host plant, and dietary RNAi. The stability of these candidate reference genes was investigated using geNorm, NormFinder, BestKeeper, and ΔCt method, respectively. Finally, RefFinder, a statistical platform integrating all four algorisms, provided a comprehensive list of genes for each extrinsic condition: (1) EF1A, α-tubulin and Actin3 were the best candidates for temperature, (2) GAPDH, 18S, and Actin3 were the most stable genes for humidity, (3) V-ATPase B, Actin3, and 18S were the top reference genes for photoperiod, (4) GAPDH, V-ATPase B, and α-tubulin were recommended for host plants, and (5) GAPDH, V-ATPase B, and RPS9 were the top choices for dietary RNAi. Overall, V-ATPase B, GAPDH, and Actin3 were the most commonly selected reference genes in A. viennensis regardless of the experimental conditions, including both intrinsic and extrinsic. Information present here lays the foundation for the genomic and functional genomic research in A. viennensis.

Screening and evaluation of the stability of expression of reference genes in Lymantria dispar (Lepidoptera: Erebidae) using qRT-PCR

The quantitative real-time polymerase chain reaction (qRT-PCR) has rapidly become the most sensitive and accurate method for the quantitative analysis of gene expression. Normalization of gene expression to that of relatively stably expressed housekeeping genes is required to facilitate the study of gene expression and to obtain more accurate RT-PCR data. However, no studies of the stability of expression of housekeeping genes in Lymantria dispar have been reported. In the present study, BestKeeper, GeNorm and NormFinder statistical software was used to evaluate the expression of thirteen candidate reference genes in L. dispar under different conditions. The expression levels of candidate reference genes were determined for two biological factors (developmental stages and tissues) and four abiotic treatments (temperature, insecticide, CO₂ and starvation). The results showed that the best candidate reference genes in L. dispar were TUB, AK, RPS15 for developmental stages, RPL32 and GAPDH for tissues, ACTB and EF1-α for CO₂ stress, GAPDH and RPL32 for temperature stress, RPS3 and GAPDH for insecticide stress, and GAPDH and RPS3 for starvation stress. In summary, EF1-α and TUB are preferential housekeeping genes in L. dispar under various conditions. These results provide a basis for the further study of functional genes of L. dispar.

Molecular Characterization of Vitellogenin and Its Receptor in Sogatella furcifera, and Their Function in Oocyte Maturation

The yolk protein precursor, vitellogenin (Vg), provides nutrition for embryonic development whereas the vitellogenin receptor (VgR) is responsible for the uptake of yolk protein by maturing oocytes. These two proteins are key reproduction-related proteins in insects. We cloned and characterized Vg and VgR genes in Sogatella furcifera, and investigated their function in oocyte maturation. Cloned SfVg and SfVgR have open reading frames of 6,114 and 5,796 bp, encoding 2,037 and 1,931 amino acid residues, respectively. Structural analysis indicates that SfVg has the three conserved LPD_N, DUF1943, and VWFD domains, SfVgR contains all conservative motifs of the LDLR superfamily. Both genes were highly expressed in adult females; SfVg was most highly expressed in the fat body whereas SfVgR was mainly expressed in the ovary. Knockdown of either gene reduced yolk protein deposition in oocytes and arrested oocyte maturation. However, silencing one of these two genes did not affect the transcript level of the other. These results demonstrate the role of SfVgR in transporting SfVg into oocytes. Both SfVg and SfVgR are essential for oocyte maturation in S. furcifera and both genes could potentially be targeted as means of controlling this pest.

Selection of Reference Genes for RT-qPCR Analysis Under Intrinsic Conditions in the Hawthorn Spider Mite, Amphitetranychus viennensis (Acarina: Tetranychidae)

Hawthorn spider mite, Amphitetranychus viennensis Zacher, is one of the most devastating pests of deciduous fruit trees. The overall goal of this research is to develop a standardized protocol for real-time quantitative reverse transcription PCR (RT-qPCR) analysis in A. viennensis following the MIQE (minimum information for publication of Quantitative real time PCR experiments) guidelines. Based on the previous knowledge, we hypothesized that internal references for RT-qPCR analysis reside in housekeeping genes (HKGs). To test this hypothesis, we examined the stability of nine HKGs from A. viennensis, including 18S ribosomal RNA (18S), 28S ribosomal RNA (28S), Elongation factor 1-alpha (EF1A), Actin3, V-ATP vacuolar-type H⁺-ATPase (V-ATPase), α-tubulin (α-tubulin), Ribosomal protein L13 (RPL13), 40S ribosomal protein S9 (RPS9), and Glyceraldehyde 3-phosphate dehydrogenase (GAPDH). The expression profile of these candidates under intrinsic conditions was evaluated by a panel of computational programs, including geNorm, Normfinder, BestKeeper, and ΔCt method. Based on RefFinder, a comprehensive software integrating all four above-mentioned algorithms, V-ATPase, Actin3, and GAPDH are the top three reference genes, which are stably expressed across all the intrinsic conditions, including developmental stage, sex, and diapause. In addition, we compared reference genes recommended for different developmental stages among the nine cell-content feeding arthropods, including four spider mites, A. viennensis, Tetranychus urticae, Tetranychus cinnabarinus, and Panonychus citri, and five hemipterans, Myzus persicae, Aphis gossypii, Toxoptera citricida, Lipaphis erysimi, and Sogatella furcifera. Not surprisingly, rRNAs and ribosomal proteins, the most abundant RNA species, is the top choice, and follows by EF1A, Actin, GAPDH, and tubulin. Information present here lays the foundation for the genomic and functional genomic research in cell-content feeding arthropods in general and A. viennensis in particular.

Determination of Suitable RT-qPCR Reference Genes for Studies of Gene Functions in Laodelphax striatellus (Fallén)

The reverse transcription quantitative polymerase chain reaction (RT-qPCR) has been widely used to determine gene functions in Laodelphax striatellus (Fallén) (small brown planthopper). Selection of suitable reference gene(s) for normalizations of RT-qPCR data is critical for reliable results. To date, reports on identification of suitable L. striatellus reference genes are still very limited. L. striatellus is a destructive rice pest and it can transmit multiple viruses, including Rice black-streaked dwarf virus (RBSDV), Rice stripe virus (RSV), and Maize rough dwarf virus (MRDV), to many important cereal crops worldwide. In this study, we examined the stablity of seven selected candidate reference genes in L. striatellus at different developmental stages, in different tissues, in RBSDV- or RSV-infected L. striatellus or in RBSDV-infected and Lssynaptojanin 1 (LsSYNJ1)-silenced L. striatellus. The RT-qPCR data representing individual candidate genes were analyzed using five different methods: the delta Ct method, geNorm, NormFinder, BestKeeper, and the RefFinder algorithm, respectively. The most stable reference gene for the specific condition was selected according to a comprehensive analysis using the RefFinder method. Ribosomal protein L5 (LsRPL5) and LsRPL8 are the most stably expressed genes in L. striatellus at different developmental stages. Alpha-1-tubulin (Lsα-TUB) is the most stably expressed reference gene in different tissues of RBSDV viruliferous (RBSDV-V) or non-viruliferous (RBSDV-NV) L. striatellus. LsRPL8 is the most stably expressed reference gene in RBSDV-V or RSV viruliferous (RSV-V) L. striatellus, while beta-tubulin (Lsβ-TUB) is the most stably expressed reference gene in RBSDV-V and LsSYNJ1-silenced L. striatellus. The selected reference genes were further investigated during analyses of RBSDV P5-1 and P10 gene expression in different tissues from RBSDV-V or RBSDV-NV L. striatellus. The stably expressed reference genes identified in this study will benefit future gene function studies using L. striatellus.

Transcriptome Analysis of Sogatella furcifera (Homoptera: Delphacidae) in Response to Sulfoxaflor and Functional Verification of Resistance-Related P450 Genes

The white-back planthopper (WBPH), Sogatella furcifera, is a major rice pest in China and in some other rice-growing countries of Asia. The extensive use of pesticides has resulted in severe resistance of S. furcifera to variety of chemical insecticides. Sulfoxaflor is a new diamide insecticide that acts on nicotinic acetylcholine receptors (nAChRs) in insects. The aim of this study was to explore the key genes related to the development of resistance to sulfoxaflor in S. furcifera and to verify their functions. Transcriptomes were compared between white-back planthoppers from a susceptible laboratory strain (Sus-Lab) and Sus-Lab screened with the sublethal LC₂₅ dose of sulfoxaflor for six generations (SF-Sel). Two P450 genes (CYP6FD1 and CYP4FD2) and three transcription factors (NlE78sf, C2H2ZF1 and C2H2ZF3) with upregulated expression verified by qRT-PCR were detected in the Sus-Lab and SF-Sel strains. The functions of CYP6FD1 and CYP4FD2 were analyzed by RNA interference, and the relative normalized expressions of CYP6FD1 and CYP4FD2 in the SF-Sel population were lower than under dsGFP treatment after dsRNA injection. Moreover, the mortality rates of SF-Sel population treated with the LC₅₀ concentration of sulfoxaflor after the injecting of dsRNA targeting CYP6FD1 and CYP4FD2 were significantly higher than in the dsGFP group from 72 h to 96 h (p < 0.05), and mortality in the CYP6FD1 knockdown group was clearly higher than that of the CYP4FD2 knockdown group. The interaction between the tertiary structures of CYP6FD1 and CYP4FD2 and sulfoxaflor was also predicted, and CYP6FD1 showed a stronger metabolic ability to process sulfoxaflor. Therefore, overexpression of CYP6FD1 and CYP4FD2 may be one of the primary factors in the development of sulfoxaflor resistance in S. furcifera.

Transcriptome profiling and dimorphic expression of sex-related genes in fifth-instar nymphs of Sogatella furcifera, an important rice pest

Sogatella furcifera is an important rice pest. In order to understand the molecular basis of the sex determination in this pest, we performed de novo transcriptome sequencing of six cDNA libraries (three biological replicates) of female and male fifth-instar nymphs. Total 65,199 unigenes were obtained, with an average length of 971.5 bp and N50 length of 1708 bp. 20,287 open reading frames (ORFs) were predicted and annotated. Total 1019 differentially expressed genes with 873 upregulated and 146 downregulated were found in male compared to female. Total 164 sex-determining genes were identified, including the key sex-determining genes in fruit flies, such as Sxl, tra, dsx, etc. It implied that the sex determination mechanisms of S. furcifera may be the same as that of fruit flies. This study provided transcriptome resource as a fundamental support for future functional studies to elucidate the sex determination regulatory networks governing sexual dimorphism of S. furcifera.

Three odorant-binding proteins are involved in the behavioral response of Sogatella furcifera to rice plant volatiles

Plant volatiles play an important role in regulating insect behavior. Odorant binding proteins (OBPs) are involved in the first step of the olfactory signal transduction pathway and plant volatiles recognition. Sogatella furcifera is one of the most destructive pests of rice crops. Understanding the functions of S. furcifera OBPs (SfurOBPs) in the host plant location and the behavioral responses of S. furcifera to rice plant volatiles could lead to improved, more environmentally-friendly, methods for controlling this pest. We found that SfurOBP1 displayed only weak binding with all the tested volatiles. SfurOBP2, SfurOBP3 and SfurOBP11 had different binding affinities to β-ionone. SfurOBP2 and SfurOBP11 had strong binding affinities to β-caryophyllene (K_i = 2.23 µM) and plant alcohol (K_i = 2.98 µM), respectively. The results of Y-olfactometer experiments indicate that S. furcifera was significantly repelled by octanal and n-octane but strongly attracted by (+)-limonene, acetophenone, 2-heptanone, n-hendecane, α-farnesene and β-ionone. Furthermore, the dsRNA-mediated gene silencing of SfurOBP2, SfurOBP3 and SfurOBP11 shifted the olfactory behavior of S. furcifera for β-ionone, α-farnesene and plant alcohol, respectively. These results suggest that the SfurOBPs are involved in the recognition of rice plant volatiles, and several potential repellants and lures for controlling this pest.

Selection of Reference Genes for Optimal Normalization of Quantitative Real-Time Polymerase Chain Reaction Results for Diaphorina citri Adults

The Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Liviidae), can cause direct damage to citrus trees and is the main vector for the devastating disease, citrus greening disease or huanglongbing. Most molecular studies on this important insect pest use real-time reverse-transcription quantitative polymerase chain reaction (RT-qPCR) to quantify gene expression, including analyzing molecular basis for insecticide resistance in field populations. One critical factor to cause inaccuracy in RT-qPCR results is the lack of appropriate internal reference genes for optimal data normalization. In this study, the expression levels of 10 selected reference genes were evaluated in different tissue samples of psyllid adults and in the insects treated with different temperatures and insecticides. Data were analyzed using different computational algorithms, including Delta Ct, BestKeeper, NormFinder, geNorm, and RefFinder. According to our results, at least two reference genes should be used for the normalization of RT-qPCR data in this insect. The best choices of reference genes for different samples are as follows: ACT1 and Ferritin for different tissue samples, RPS20 and Ferritin for samples treated with different temperatures, TBP and EF1α for samples treated with imidacloprid, and Ferritin and TBP for samples treated with beta-cypermethrin. The reference genes identified in this study should be useful for future studies to analyze the expression patterns of target genes, especially for genes linked with temperature adaptability and insecticide resistance in this insect species in the future.

Transcriptomic and Expression Analysis of the Salivary Glands in Brown Planthoppers, Nilaparvata lugens (Hemiptera: Delphacidae)

The brown planthopper (BPH), Nilaparvata lugens Stål (Hemiptera: Delphacidae), is a serious rice pest because of its destructive feeding. The salivary glands of the BPH play a key role in successful feeding. In this study, we explored the salivary gland transcriptome (sialotranscriptome) of adult BPHs using Illumina sequencing and a total of 55,913 transcripts and 45,421 unigenes were obtained. We identified one reference gene RPL9 (Ribosomal protein L9) and 19 salivary protein genes from the BPH sialotranscripome, which were categorized as those involved in sugar metabolism, extra-oral digestion of cell wall components, detoxification, and suppression of plant defenses. Tissue expression profiles of 19 salivary protein genes analysis revealed that the expression level of alpha-glucosidase family 31 had no difference in five tissues, suggesting that it may have functions in the whole-body parts. Glucose dehydrogenase (flavine adenine dinucleotide, quinone)-like was expressed highly in the salivary gland, which might play putative role in insect feeding. Glucose dehydrogenase (acceptor) was expressed the highest level in head without salivary gland. Other salivary protein genes were highly expressed in gut or malpighian tubule, suggesting that they may play roles in digestive and detoxification mechanism. Moreover, we detected RPL9 is one of the genes that is most consistently present for acquisition of gene expression in different tissues. Thus, RPL9 can be a new reference for expression studies of BPH. The obtained BPH sialotranscripome provides a list of genes that have potential roles in feeding and interaction between BPHs and rice plants.

Selection of Reference Genes for the Normalization of RT-qPCR Data in Gene Expression Studies in Insects: A Systematic Review

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for quantifying expression levels of targeted genes during various biological processes in numerous areas of clinical and biological research. Selection of appropriate reference genes for RT-qPCR normalization is an elementary prerequisite for reliable measurements of gene expression levels. Here, by analyzing datasets published between 2008 and 2017, we summarized the current trends in reference gene selection for insect gene expression studies that employed the most widely used SYBR Green method for RT-qPCR normalization. We curated 90 representative papers, mainly published in 2013–2017, in which a total of 78 insect species were investigated in 100 experiments. Furthermore, top five journals, top 10 frequently used reference genes, and top 10 experimental factors have been determined. The relationships between the numbers of the reference genes, experimental factors, analysis tools on the one hand and publication date (year) on the other hand was investigated by linear regression. We found that the more recently the paper was published, the more experimental factors it tended to explore, and more analysis tools it used. However, linear regression analysis did not reveal a significant correlation between the number of reference genes and the study publication date. Taken together, this meta-analysis will be of great help to researchers that plan gene expression studies in insects, especially the non-model ones, as it provides a summary of appropriate reference genes for expression studies, considers the optimal number of reference genes, and reviews the average number of experimental factors and analysis tools per study.

Reference genes selection for quantitative gene expression studies in tea green leafhoppers, Empoasca onukii Matsuda

Empoasca onukii Matsuda is one of the most devastating pests of the tea plant (Camellia sinensis). Still, the presumed expression stability of its reference genes (RGs) has not been analyzed. RGs are essential for accurate and reliable gene expression analysis, so this absence has hampered the study of the insect’s molecular biology. To find candidate RGs for normalizing gene expression data, we cloned ten common housekeeping genes from E. onukii. Using the ΔCt method, geNorm, NormFinder and BestKeeper, we screened the RGs that were appropriate for quantifying the mRNA transcription of cellular responses under five experimental conditions. We identified the combinations of α-TUB and G6PDH, α-TUB and UBC, two RGs (α-TUB and β-TUB1) or three RGs (α-TUB, RPL13 and GAPDH), AK and UBC, or RPL13 and α-TUB as the best for analyzing gene expression in E. onukii adults of both sexes in different tissues, nymphs at different developmental stages, nymphs exposed to different temperatures or nymphs exposed to photoperiod stress. Finally, the E. onukii cysteine proteinase (Eocyp) was chosen as the target gene to validate the rationality of the proposed RGs. In conclusion, our study suggests a series of RGs with which to study the gene expression profiles of E. onukii that have been manipulated (biotically or abiotically) using reverse transcription quantitative polymerase chain reaction. The results offer a solid foundation for further studies of the molecular biology of E. onukii.

Selection and evaluation of reference genes for expression analysis using quantitative real-time PCR in the Asian Ladybird Harmonia axyridis (Coleoptera: Coccinellidae)

Harmonia axyridis (Coleoptera: Coccinellidae) is a polyphagous insect that is an important biological agent used to control agricultural and forestry pests. The role of functional genes in H. axyridis based on quantitative real-time PCR (qRT-PCR) is increasingly well understood to investigate biology, physiology, feeding behavior and the role of important genes in physiological processes. Quantitative real-time PCR (qRT-PCR) is a powerful and reliable technique to quantify gene expression. Using qRT-PCR, expression levels of target genes are determined based on the levels of internal reference genes; therefore, reference genes need to be stably expressed under specific experimental conditions. However, there have been no studies on the stability of reference genes used in H. axyridis. In this study, we systematically investigated expression profiles of nine candidate reference genes from H. axyridis, including β-actin (ACTIN); elongation factor 1 α (EF1A); ribosomal proteins L10, L18, L28, S13, and S15 (RPL10, RPL18, RPL28, RPS13 and RPS15); glyceraldehyde-3-phosphate dehydrogenase (GAPDH); and superoxide dismutase (SOD). Four analytical methods (geNorm, BestKeeper, NormFinder, and the ΔCt method) were used to evaluate the suitability of these genes as internal reference genes for three biotic factors (developmental stage, tissue, and sex) and two abiotic treatments (temperature and photoperiod). RefFinder, a comprehensive evaluation platform integrating the four analytical methods, was used to rank the overall stability of these reference genes. Among the nine candidate genes, different reference genes were identified as having the most stable expression across biotic and abiotic factors. Genes encoding ribosomal proteins typically had the most stable expression, though EF1A was the most stable across developmental stages and photoperiods. To validate the suitability of these reference genes, heat shock protein 90 (HSP90) was chosen as a target. Significant up-regulation in HSP90 expression level in response to both low and high temperature was observed when using the most suitable reference genes but not when using an arbitrarily selected reference gene. The reference genes identified in this study will provide the basis for future functional genomics research in H. axyridis and will also facilitate the establishment of a standardized qRT-PCR program for other related insects.

Selection of reference genes for tissue/organ samples on day 3 fifth-instar larvae in silkworm, Bombyx mori

The silkworm, Bombyx mori, is one of the world's most economically important insect. Surveying variations in gene expression among multiple tissue/organ samples will provide clues for gene function assignments and will be helpful for identifying genes related to economic traits or specific cellular processes. To ensure their accuracy, commonly used gene expression quantification methods require a set of stable reference genes for data normalization. In this study, 24 candidate reference genes were assessed in 10 tissue/organ samples of day 3 fifth-instar B. mori larvae using geNorm and NormFinder. The results revealed that, using the combination of the expression of BGIBMGA003186 and BGIBMGA008209 was the optimum choice for normalizing the expression data of the B. mori tissue/organ samples. The most stable gene, BGIBMGA003186, is recommended if just one reference gene is used. Moreover, the commonly used reference gene encoding cytoplasmic actin was the least appropriate reference gene of the samples investigated. The reliability of the selected reference genes was further confirmed by evaluating the expression profiles of two cathepsin genes. Our results may be useful for future studies involving the quantification of relative gene expression levels of different tissue/organ samples in B. mori.

Proteomic analysis of watery saliva secreted by white-backed planthopper, Sogatella furcifera

The white-backed planthopper, Sogatella furcifera, is a phloem sap feeder that secretes watery and gelling saliva during feeding. In this study, we identified the major proteins in watery saliva of S. furcifera by shotgun LC-MS/MS analysis combined with transcriptomic analysis. A total of 161 proteins were identified, which were divided into 8 function categories, including enzymes, transporter, calcium ion binding protein, salivary sheath protein, cytoskeleton protein, DNA-, RNA-, and protein-binding or regulating proteins, other non-enzyme proteins and unknown proteins. Gene expression pattern of 11 secretory proteins were analyzed by real time quantitative-PCR. We detected the mucin-like protein, which had a unique expression level in salivary gland, most likely to be a candidate effector involved in regulation of plant defense. This study identified the watery saliva component of S. furcifera and it provided a list of proteins which may play a role in interaction between S. furcifera and rice.

Reference gene selection to determine differences in mitochondrial gene expressions in phosphine-susceptible and phosphine-resistant strains of Cryptolestes ferrugineus, using qRT-PCR

Cryptolestes ferrugineus is a serious pest of stored grain and has developed high levels of resistance to phosphine fumigants in many countries. Measuring differences in expression levels of certain 'resistant' genes by quantitative real-time PCR (qRT-PCR) may provide insights into molecular mechanisms underlying resistance to phosphine in C. ferrugineus, but reliable qRT-PCR results depend on suitable reference genes (RGs). We evaluated the stability of nine candidate RGs across different developmental stages and phosphine strains of C. ferrugineus, using four softwares. The results showed that RPS13 and EF1α were the most stable RGs, whereas α-TUB was the least under developmental stages. Across the different strains, RPS13 and γ-TUB were the most stable RGs, whereas CycA and GAPDH were the least. We confirmed the reliability of the selected RGs by qRT-PCR analyses of the mitochondrial cox1 gene. Expression of cox1 was not significantly different in the phosphine-resistant strain compared with the phosphine-susceptible strain, but three mitochondrial genes (nad3, atp6 and cob) were significantly down-regulated. These results suggest that alterations in the expressions of these three genes may be associated with phosphine resistance in C. ferrugineus. The findings will facilitate future functional genomics studies on the development and phosphine resistance in C. ferrugineus.

Selection and Evaluation of Reference Genes for Expression Analysis Using qRT-PCR in Chilo suppressalis (Lepidoptera: Pyralidae)

Quantitative real-time polymerase chain reaction (qRT-PCR) is a valuable tool for estimating gene expression; however, the validity is largely dependent on the selection of stable reference genes. The suitability of various reference genes for qRT-PCR analysis was evaluated in, Chilo suppressalis (Walker). The ΔCt method, geNorm, NormFinder, and BestKeeper were used to evaluate the suitability of nine candidate reference genes for normalizing gene expression in larval tissues and organs and during high and low temperature stress. The ΔCt method, geNorm, and NormFinder produced similar stability rankings; H3, UBI, and EF1 were the most stable reference genes for monitoring gene expression in larval tissue and organs, and EF1, TUB, and AK were the optimal genes for thermal stress. However, for thermal stress, RPS11 was the most stable gene based on BestKeeper. To validate these recommendations, the expression profile of the gene encoding heat shock protein 60 (Hsp60) was investigated. Hsp60 transcript levels showed significant differences when normalized to the most versus least stable reference genes. These results further confirm the importance of testing reference genes using the selected experimental parameters. The reference genes identified in the present study will improve the quality of gene expression data obtained for C. suppressalis and will facilitate future studies aimed at understanding the biology of this important insect pest.

Transcriptomic and Expression Analysis of the Salivary Glands in White-Backed Planthoppers, Sogatella furcifera

The white-backed planthopper (WBPH), Sogatella furcifera (Horváth), is one of the serious rice pests because of its destructive feeding. The salivary glands of the WBPH play an important role in the feeding behaviour. Currently, however, very little is known about the salivary glands at the molecular level. We sequenced the salivary gland transcriptome (sialotranscripome) of adult WBPHs using the Illumina sequencing. A total of 65,595 transcripts and 51,842 unigenes were obtained from salivary glands. According to annotations against the Nr database, many of the unigenes identified were associated with the most studied enzymes in hemipteran saliva. In the present study, we identified 32 salivary protein genes from the WBPH sialotranscripome, which were categorized as those involved in sugar metabolism, detoxification, suppression of plant defense responses, immunity-related responses, general digestion, and other phytophagy processes. Tissue expression profiles analysis revealed that four of 32 salivary protein genes (multicopper oxidase 4, multicopper oxidase 6, carboxylesterase and uridine phosphorylase 1 isform X2) were primarily expressed in the salivary gland, suggesting that they played putative role in insect-rice interactions. 13 of 32 salivary protein genes were primarily expressed in gut, which might play putative role in digestive and detoxify mechanism. Development expression profiles analysis revealed that the expression level of 26 of 32 salivary protein genes had no significant difference, suggesting that they may play roles in every developmental stages of salivary gland of WBPH. The other six genes have a high expression level in the salivary gland of adult. 31 of 32 genes (except putative acetylcholinesterase 1) have no significant difference in male and female adult, suggesting that their expression level have no difference between sexes. This report analysis of the sialotranscripome for the WBPH, and the transcriptome provides a foundational list of the genes involved in feeding. Our data will be useful to investigate the mechanisms of interaction between the WBPH and the host plant.

Identification and evaluation of reference genes in the Chinese white wax scale insect Ericerus pela

BACKGROUND

The Chinese white wax scale insect, Ericerus pela, is a well-known resource insect. The females and males are dramatically distinct at each developmental stage. We sought to identify suitable reference genes to use as internal controls in molecular research on E. plea.

RESULTS

geNorm, RefFinder and Normfinder analyses showed that ßTub-2 was the best reference gene throughout different developmental stages; SdhA-1 was the most stable reference gene in different tissues, and ßTub-1 was the most reliable reference gene under treatment with different temperatures. The results also showed that the optimal number of reference genes for analyzing target gene expression levels in the three experimental conditions was two.

CONCLUSIONS

The identified reference genes are suitable reference genes for normalization in RT-qPCR of E. pela samples.