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

Selection of a suitable reference gene for gene-expression studies in Trichomonas gallinae under various biotic and abiotic stress conditions

Trichomonas gallinae, a globally distributed protozoan parasite, significantly affects the pigeon-breeding industry. T. gallinae infection mainly causes yellow ulcerative nodules on the upper respiratory tract and crop mucosa of pigeons, impeding normal breathing and feeding and ultimately causing death. Real-time quantitative PCR (qPCR) is a crucial technique for gene-expression analysis in molecular biology. Reference-gene selection for normalization is critical for ensuring this technique's accuracy. However, no systematic screening or validation of T. gallinae reference genes has been reported. This study quantified the transcript levels of ten candidate reference genes in T. gallinae isolates with different genotypes and culture conditions using qPCR. Using the geNorm, NormFinder, and BestKeeper algorithms, we assessed these reference genes' stabilities and ranked them using RankAggreg analysis. The most stable reference gene was tubulin beta chain (TUBB), while the widely used reference genes TUBG and GAPDH demonstrated poor stability. Additionally, we evaluated these candidate reference genes' stabilities using the T. gallinae TgaAtg8 gene. On using TUBB as a reference gene, TgaAtg8's expression profiles in T. gallinae isolates with different genotypes remained relatively consistent under various culture conditions. Conversely, using ACTB as a reference gene distorted the data. These findings provide valuable reference-gene-selection guidance for functional gene research and gene-expression analysis in T. gallinae.

Evaluation and Validation of Reference Genes for Gene Expression Analysis Using qRT-PCR in the Sugarcane Stem Borer Chilo sacchariphagus (Lepidoptera: Pyralidae)

Chilo sacchariphagus (Lepidoptera: Pyralidae) is an economically important sugarcane pest. Although numerous studies were conducted on the physiological responses in C. sacchariphagus, little is known regarding the genes regulating these physiological processes. Gene expression analysis by qRT-PCR can offer a significant indication for functional gene studies. To our knowledge, the reference genes of C. sacchariphagus have not been screened or evaluated, which hinders the functional gene study. In the present study, the stability of seven reference genes (β-ACT, GAPDH, BTF3, 28S, RPL7, EF1α, and SDHA) was evaluated in C. sacchariphagus under different experimental conditions, including tissues (antenna, head, thorax, abdomen, leg, and wing), temperatures (4 °C, 25 °C, and 37 °C) and sexes (male and female), through RefFinder, which integrates four algorithms (Normfinder, BestKeeper, ΔCt method, and geNorm). The findings suggested that the combination of β-ACT and RPL7 is ideal to analyze gene expressions in different tissues and at distinct temperatures, and EF1α and SDHA were suitable reference genes for comparing gene expressions between sexes. Finally, the expression profiles of CsacPBP1 gene were evaluated, and the outcomes further confirm the importance of selecting fitting reference genes for normalization of qRT-PCR data. This study represents the first kind in screening out suitable reference genes for gene expression analysis in C. sacchariphagus. Information from this study is poised to galvanize future inquiry into the gene expression of C. sacchariphagus, an economically important pest of sugarcane.

Identification of specific reference gene for normalization of RT-qPCR data in rhythmic gene expression studies of the effect of developmental hormone antagonist in postembryonic development in Bombyx mori

Bombyx mori is a lepidopteran holometabolous insect with distinct developmental stages: egg, larvae, pupae, and adult. The lepidopteran insect undergoes major modifications in the central nervous system (CNS) so as to adapt to the lifestyle of these distinct stages with specific habitats and functions from voraciously feeding larval stages to flying reproductive adults via dormant pupal stages. Such transitions are linked to transcriptional, epigenetic, and translational complexities. Therefore, studying rhythmic gene expression in CNS of various developmental stages and the effects of antagonists on developmental hormones requires a very stable reference gene (RG). To facilitate rhythmic gene expression studies using reverse transcription quantitative polymerase chain reaction (RT-qPCR) in B. mori and the effect of developmental hormone juvenile hormone (JH) and 20-hydroxy ecdysone hormone (20 HE), antagonists Precocene 1 and testosterone, respectively, were used. Eight candidate RGs, namely, Translational initiation factor 3 subunit 4 (TI3S4), Translational initiation factor 3 subunit 5 (TI3S5), Ribosomal protein subunit 7 (RPs7), TATA-binding protein association factor (TAF13), Translational initiation factor 4 A (TI4A), Ribosomal protein (RPL32), Elongation factor 1 (EF1), and Arginine kinase (AK), were assessed in the CNS of B. mori. The postembryonic developmental (PED) stages used were the fifth late larval instar, early pupa, mid pupa, late pupa, and adult. The assessments were done at four different time points, Zeitgeber time (ZT) 0, 6, 12, and 18, to find stability towards 24-h rhythmic expression. RefFinder, geNorm, and Ct value analysis were performed. RefFinder and geNORM studies suggested stability order as TI3S4 > TI3S5 > RPs7, but Ct value evaluation showed stability order as TI3S5 > TI3S4 > RPs7. We therefore demonstrated that TI3S4, TI3S5, and RPs7 can be used as RG in various PED stages in CNS of B. mori (Strain: CB-hybrid, PM×CSR2) towards studies with effects of JH and 20 HE antagonists.

Biochemistry and transcriptomic analyses of Phthorimaea absoluta (Lepidoptera: Gelechiidae) response to insecticides

Phthorimaea absoluta is an invasive solanaceous plant pest with highly devastating effects on tomato plant. Heavy reliance on insecticide use to tackle the pest has been linked to insecticide resistance selection in P. absoluta populations. To underline insights on P. absoluta insecticide resistance mechanisms to diamides and avermectins, we evaluated the transcriptomic profile of parental (field-collected) and F8 (lab-reared) populations. Furthermore, to screen for the presence of organophosphate and pyrethroid resistance, we assessed the gene expression levels of acetylcholinesterase (ace1) and para-type voltage-gated sodium channel (VGSG) genes in the F1 to F8 lab-reared progeny of diamide and avermectin exposed P. absoluta field-collected populations. The VGSG gene showed up-regulation in 12.5% and down-regulation in 87.5% of the screened populations, while ace1 gene showed up-regulation in 37.5% and down-regulation in 62.5% of the screened populations. Gene ontology of the differentially expressed genes from both parental and eighth generations of diamide-sprayed P. absoluta populations revealed three genes involved in the metabolic detoxification of diamides in P. absoluta. Therefore, our study showed that the detoxification enzymes found could be responsible for P. absoluta diamide-based resistance, while behavioural resistance, which is stimulus-dependent, could be attributed to P. absoluta avermectin resistance.

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.

Nanoparticle-Shielded dsRNA Delivery for Enhancing RNAi Efficiency in Cotton Spotted Bollworm Earias vittella (Lepidoptera: Nolidae)

The spotted bollworm Earias vittella (Lepidoptera: Nolidae) is a polyphagous pest with enormous economic significance, primarily affecting cotton and okra. However, the lack of gene sequence information on this pest has a significant constraint on molecular investigations and the formulation of superior pest management strategies. An RNA-seq-based transcriptome study was conducted to alleviate such limitations, and de novo assembly was performed to obtain transcript sequences of this pest. Reference gene identification across E. vittella developmental stages and RNAi treatments were conducted using its sequence information, which resulted in identifying transcription elongation factor (TEF), V-type proton ATPase (V-ATPase), and Glyceraldehyde -3-phosphate dehydrogenase (GAPDH) as the most suitable reference genes for normalization in RT-qPCR-based gene expression studies. The present study also identified important developmental, RNAi pathway, and RNAi target genes and performed life-stage developmental expression analysis using RT-qPCR to select the optimal targets for RNAi. We found that naked dsRNA degradation in the E. vittella hemolymph is the primary reason for poor RNAi. A total of six genes including Juvenile hormone methyl transferase (JHAMT), Chitin synthase (CHS), Aminopeptidase (AMN), Cadherin (CAD), Alpha-amylase (AMY), and V-type proton ATPase (V-ATPase) were selected and knocked down significantly with three different nanoparticles encapsulated dsRNA conjugates, i.e., Chitosan-dsRNA, carbon quantum dots-dsRNA (CQD-dsRNA), and Lipofectamine-dsRNA conjugate. These results demonstrate that feeding nanoparticle-shielded dsRNA silences target genes and suggests that nanoparticle-based RNAi can efficiently manage this pest.

The uncommon function and mechanism of the common enzyme glyceraldehyde-3-phosphate dehydrogenase in the metamorphosis of Helicoverpa armigera

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), a key enzyme in glycolysis, is commonly used as an internal reference gene in humans, mice, and insects. However, the function of GAPDH in insect development, especially in metamorphosis, has not been reported. In the present study, Helicoverpa armigera and Spodoptera frugiperda ovarian cell lines (Sf9 cells) were used as materials to study the function and molecular mechanism of GAPDH in larval metamorphosis. The results showed that HaGAPDH was more closely related to GAPDH of S. frugiperda and Spodoptera litura. The transcript peaks of HaGAPDH in sixth instar larvae were 6L-3 (epidermal and midgut) and 6L-1 (fat body) days, and 20E and methoprene significantly upregulated the transcripts of HaGAPDH of larvae in qRT-PCR. HaGAPDH–GFP–His was specifically localized in mitochondria in Sf9 cells. Knockdown of HaGAPDH by RNA interference (RNAi) in sixth instar larvae resulted in weight loss, increased mortality, and decreases in the pupation rate and emergence rates. HaGAPDH is directly bound to soluble trehalase (HaTreh1) physically and under 20E treatment in yeast two-hybrid, coimmunoprecipitation, and colocalization experiments. In addition, knockdown of HaGAPDH increased the Treh1 activity, which in turn decreased the trehalose content but increased the glucose content in larvae. Therefore, these data demonstrated that GAPDH controlled the glucose content within the normal range to ensure glucose metabolism and metamorphosis by directly binding with HaTreh1.

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 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.

Reference Genes for Expression Analyses by qRT-PCR in Phthorimaea operculella (Lepidoptera: Gelechiidae)

Due to a lack of effective internal references, studies on functional genes in Phthorimaea operculella, a serious Lepidopteran pest attacking potatoes worldwide, have been greatly limited. To select suitable endogenous controls, ten housekeeping genes of actin (ACT), α-tubulin (α-TUB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), elongation factor 1α (EF1α), 18S and 28S ribosomal RNA (18S, 28S), ribosomal protein genes RPL4, RPL13 and RPL27 and superoxide dismutase (SOD) were tested. Their expression levels were determined under three different experimental conditions (developmental stages, tissues/organs and temperatures) using qRT-PCR technology. The stability was evaluated with five methods (Ct value, geNorm, NormFinder, BestKeeper and RefFinder). The results clarified that RPL13, EF1α and RPL27 are ranked as the best reference gene combination for measuring gene expression levels among different developing stages and under various temperatures; EF1α and RPL13 are recommended to normalize the gene expression levels among diverse tissues. EF1α and RPL13 are the best reference genes in all the experimental conditions. To validate the utility of the selected reference pair, EF1α and RPL13, we estimated the tissue-biased expression level of chitin synthase A gene (PoChSA). As expected, PoChSA was abundantly expressed in ectodermally derived epidermal cells, and lowly transcribed in the midgut. These findings will lay the foundation for future research on the molecular physiology and biochemistry of P. operculella.

Evaluation of Optimal Reference Genes for qRT-PCR Analysis in Hyphantria cunea (Drury)

The relative quantification of gene expression is mainly achieved through reverse transcription-quantitative PCR (qRT-PCR); however, its reliability and precision rely on proper data normalization using one or more optimal reference genes. Hyphantria cunea (Drury) has been an invasive pest of forest trees, ornamental plants, and fruit trees in China for many years. Currently, the molecular physiological role of reference genes in H. cunea is unclear, which hinders functional gene study. Therefore, eight common reference genes, RPS26, RPL13, UBI, AK, RPS15, EIF4A, β-actin, α-tub, were selected to evaluate levels of gene expression stability when subjected to varied experimental conditions, including developmental stage and gender, different tissues, larvae reared on different hosts and different larval density. The geNorm, BestKeeper, ΔCt method, and NormFinder statistical algorithms were used to normalize gene transcription data. Furthermore, the stability/suitability of these candidates was ranked overall by RefFinder. This study provides a comprehensive evaluation of reference genes in H. cunea and could help select reference genes for other Lepidoptera species.

Evaluation of Reference Genes in Glenea cantor (Fabricius) by Using qRT-PCR

Kapok is the main host of Glenea cantor (Fabricius), which causes serious damage and is difficult to control. In severe cases, it often causes the kapok trees to die continuously, which seriously affects the results of urban landscaping. To provide reference for the functional research on related genes in G. cantor, we screened the stable expression of candidate reference genes at different developmental stages (i.e., eggs, larvae, pupae, and adults), in various adult tissues (i.e., head, thorax, abdomen, feet, antennae, and wings), and sexes (i.e., male pupae, female pupae, male adults, and female adults). In this study, 12 candidate reference genes (i.e., ACTINLIKE, ACTININ, TUB, RPL36, RPL32, RPS20, TBP, GAPDH, 18S rRNA, EF1A1, EF1A2, and UBQ) were evaluated using different adult tissues, developmental stages, and sexes. RefFinder, geNorm, NormFinder, and BestKeeper were used to evaluate and comprehensively analyze the stability of the expression of the candidate reference genes. The results show that RPL32 and EF1A1 were the most suitable reference genes in the different adult tissues, and RPL36 and EF1A1 were best at the different developmental stages. RPL36 and EF1A2 were the best fit for the qRT-PCR reference genes in the different sexes, while RPL36 and EF1A1 were the most appropriate qRT-PCR reference genes in all samples. Results from geNorm showed that the optimal number of reference genes was two. We also surveyed the expression of cellulase at the different developmental stages and in the different adult tissues. Results further verified the reliability of the reference genes, and confirmed the best reference genes under the different experimental conditions. This study provides a useful tool for molecular biological studies on G. cantor.