Selection of Reference Genes for Normalization of Gene Expression in Thermobia domestica (Insecta: Zygentoma: Lepismatidae)

Evolution of insect metamorphosis - an update

Metamorphosis endowed the insects with properties that enabled them to conquer the Earth. It is a hormonally controlled morphogenetic process that transforms the larva into the adult. Metamorphosis appeared with the origin of wings and flight. The sesquiterpenoid juvenile hormone (JH) suppresses wing morphogenesis and ensures that metamorphosis takes place at the right ontogenetic time. This review explores the origin of insect metamorphosis and the ancestral function of JH. Fossil record shows that the first Paleozoic winged insects had (hemimetabolous) metamorphosis, and their larvae were likely aquatic. In the primitive wingless silverfish that lacks metamorphosis, JH is essential for late embryogenesis and reproduction. JH production after the embryo dorsal closure promotes hatching and terminal tissue maturation.

E93 is indispensable for reproduction in ametabolous and hemimetabolous insects

Ecdysone-induced protein 93 (E93), known as the 'adult-specifier' transcription factor in insects, triggers metamorphosis in both hemimetabolous and holometabolous insects. Although E93 is conserved in ametabolous insects, its spatiotemporal expression and physiological function remain poorly understood. In this study, we first discover that, in the ametabolous firebrat Thermobia domestica, the previtellogenic ovary exhibits cyclically high E93 expression, and E93 mRNA is broadly distributed in previtellogenic ovarioles. E93 homozygous mutant females of T. domestica exhibit severe fecundity deficiency due to impaired previtellogenic development of the ovarian follicles, likely because E93 induces the expression of genes involved in ECM (extracellular matrix)-receptor interactions during previtellogenesis. Moreover, we reveal that in the hemimetabolous cockroach Blattella germanica, E93 similarly promotes previtellogenic ovarian development. In addition, E93 is also essential for vitellogenesis that is necessary to guarantee ovarian maturation and promotes the vitellogenesis-previtellogenesis switch in the fat body of adult female cockroaches. Our findings deepen the understanding of the roles of E93 in controlling reproduction in insects, and of E93 expression and functional evolution, which are proposed to have made crucial contributions to the origin of insect metamorphosis.

Juvenile hormone signaling is indispensable for late embryogenesis in ametabolous and hemimetabolous insects

BACKGROUND

Juvenile hormone (JH) is an insect-exclusive hormone involved in regulating diverse aspects of insect physiology, and the evolution of its diverse function is widely interesting. Studying embryogenesis in basal wingless insects is important to understand the functional evolution of JH; however, experimental studies in this regard are scarce. In this study, we conducted CRISPR/Cas9-mediated knockout (KO) of genes involved in JH biosynthesis and signaling cascades in the ametabolous firebrat, Thermobia domestica. Additionally, we investigated whether the primitive action of JH is conserved in the hemimetabolous cricket, Gryllus bimaculatus.

RESULTS

We observed that KO of JHAMT, CYP15A1, Met, and Kr-h1 resulted in embryonic lethality in T. domestica. Deprivation of JH or JH signaling arrested the progression of extraembryonic fluid resorption after dorsal closure and hatching, which is consistent with the gene expression pattern showing high Kr-h1 expression in the late embryos of T. domestica. The embryos deficient in JH signaling displayed wrinkled and weak legs. Comparative transcriptome analysis revealed that JH signaling promotes embryonic leg maturation through inducing energy supply and muscle activity, as validated by transmission electron microscopy (TEM). In addition, JH signaling exhibited similar embryonic effects in G. bimaculatus.

CONCLUSIONS

This study reveals the indispensable role of JH signaling in facilitating the maturation of terminal tissues during late embryogenesis, as demonstrated by the regulation of leg development, in ametabolous and hemimetabolous insects. These findings further indicate that the embryonic functions of JH evolved earlier than its anti-metamorphic functions during postembryonic development.

Identification and validation of reference genes for RT-qPCR analysis in Sclerodermus guani (Hymenoptera: Bethylidae)

Gene expression studies in organisms are often conducted using reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), and the accuracy of RT-qPCR results relies on the stability of reference genes. We examined ten candidate reference genes in Sclerodermus guani, a parasitoid wasp that is a natural enemy of long-horned beetle pests in forestry, including ACT, EF1α, Hsc70, Hsp70, SRSF7, α-tubulin, RPL7A, 18S, 28S, and SOD1, regarding variable biotic and abiotic factors such as body part, life stage, hormone, diet, and temperature. Data were analysed using four dedicated algorithms (ΔCt, BestKeeper, NormFinder, and geNorm) and one comparative tool (RefFinder). Our results showed that the most stable reference genes were RPL7A and EF1α regarding the body part, SRSF7 and Hsc70 regarding the diet, RPL7A and α-tubulin regarding the hormone, SRSF7 and RPL7A regarding the life stage, and SRSF7 and α-tubulin regarding temperature. To ascertain the applicability of specific reference genes, the expression level of the target gene (ACPase) was estimated regarding the body part using the most stable reference genes, RPL7A and EF1α, and the least stable one, SOD1. The highest expression level of ACPase was observed in the abdomen, and the validity of RPL7A and EF1α was confirmed. This study provides, for the first time, an extensive list of reliable reference genes for molecular biology studies in S. guani.

Depletion of lipid storage droplet-1 delays endoreplication progression and induces cell death in Drosophila salivary gland

Perilipins are evolutionarily conserved from insects to mammals. Drosophila lipid storage droplet-1 (LSD-1) is a lipid storage droplet membrane surface-binding protein family member and a counterpart to mammalian perilipin 1 and is known to play a role in lipolysis. However, the function of LSD-1 during specific tissue development remains under investigation. This study demonstrated the role of LSD-1 in salivary gland development. Knockdown of Lsd-1 in the salivary gland was established using the GAL4/UAS system. The third-instar larvae of knockdown flies had small salivary glands containing cells with smaller nuclei. The null mutant Drosophila also showed the same phenotype. The depletion of LSD-1 expression induced a delay of endoreplication due to decreasing CycE expression and increasing DNA damage. Lsd-1 genetically interacted with Myc in the third-instar larvae. These results demonstrate that LSD-1 is involved in cell cycle and cell death programs in the salivary gland, providing novel insight into the effects of LSD-1 in regulating salivary gland development and the interaction between LSD-1 and Myc.

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.

The embryonic role of juvenile hormone in the firebrat, Thermobia domestica, reveals its function before its involvement in metamorphosis

To gain insights into how juvenile hormone (JH) came to regulate insect metamorphosis, we studied its function in the ametabolous firebrat, Thermobia domestica. Highest levels of JH occur during late embryogenesis, with only low levels thereafter. Loss-of-function and gain-of-function experiments show that JH acts on embryonic tissues to suppress morphogenesis and cell determination and to promote their terminal differentiation. Similar embryonic actions of JH on hemimetabolous insects with short germ band embryos indicate that JH's embryonic role preceded its derived function as the postembryonic regulator of metamorphosis. The postembryonic expansion of JH function likely followed the evolution of flight. Archaic flying insects were considered to lack metamorphosis because tiny, movable wings were evident on the thoraces of young juveniles and their positive allometric growth eventually allowed them to support flight in late juveniles. Like in Thermobia, we assume that these juveniles lacked JH. However, a postembryonic reappearance of JH during wing morphogenesis in the young juvenile likely redirected wing development to make a wing pad rather than a wing. Maintenance of JH then allowed wing pad growth and its disappearance in the mature juvenile then allowed wing differentiation. Subsequent modification of JH action for hemi- and holometabolous lifestyles are discussed.

Evaluation and normalization of a set of reliable reference genes for quantitative sgk-1 gene expression analysis in Caenorhabditis elegans-focused cancer research

Multiple signaling pathways have been discovered to play a role in aging and longevity, including the insulin/IGF-1 signaling system, AMPK pathway, TOR signaling, JNK pathway, and germline signaling. Mammalian serum and glucocorticoid-inducible kinase 1 (sgk-1), which has been associated with various disorders including hypertension, obesity, and tumor growth, limits survival in C. elegans by reducing DAF-16/FoxO activity while suppressing FoxO3 activity in human cell culture. C. elegans provides significant protection for a number of genes associated with human cancer. The best known of these are the lin-35/pRb (mammalian ortholog pRb) and CEP-1 (mammalian ortholog p53) genes. Therefore, in this study, we aimed to investigate the expression analyzes of sgk-1, which is overexpressed in many types of mammalian cancer, in mutant lin-35 and to demonstrate the validation of reference genes in wild-type N2 and mutant lin-35 for C. elegans-focused cancer research. To develop functional genomic studies in C. elegans, we evaluated the expression stability of five candidate reference genes (act-1, ama-1, cdc-42, pmp-3, iscu-1) by quantitative real-time PCR using five algorithms (geNorm, NormFinder, Delta Ct method, BestKeeper, RefFinder) in N2 and lin-35 worms. According to our findings, act-1 and cdc-42 were effective in accurately normalizing the levels of gene expression in N2 and lin-35. act-1 and cdc-42 also displayed the most consistent expression patterns, therefore they were utilized to standardize expression level of sgk-1. Furthermore, our results clearly showed that sgk-1 was upregulated in lin-35 worms compared to N2 worms. Our results highlight the importance of definitive validation using mostly expressed reference genes.

Selection of reference genes for quantitative real-time PCR normalization in the coffee white stem borer, Xylotrechus quadripes Chevrolat (Coleoptera: Cerambycidae)

The coffee white stem borer, Xylotrechus quadripes Chevrolat (Coleoptera: Cerambycidae), is a major destructive pest of Coffea arabica L. (Gentianales: Rubiaceae), widely planted in many Asian countries, including China. Quantitative real-time polymerase chain reaction (qRT-PCR) is a common method for quantitative analysis of gene transcription levels. To obtain accurate and reliable qRT-PCR results, it is necessary to select suitable reference genes to different experimental conditions for normalizing the target gene expression. However, the stability of the expression of reference genes in X. quadripes has rarely been studied. In this study, the expression stability of nine candidate reference genes were investigated under biotic and abiotic conditions for use in qRT-PCR's normalization. By integrating the results of four algorithms of NormFinder, BestKeeper, geNorm, and RefFinder, the optimal reference gene combinations in different experimental conditions were performed as follows: RPL10a and EIF3D were the optimal reference genes for developmental stage samples, EIF4E, RPL10a, and RPS27a for tissue samples, V-ATP and EF1α for the sex samples, EIF3D and V-ATP for temperature treatment, RPS27a and RPL10a for insecticide stress, and RPL10a, RPS27a, and EF1α for all the samples. This study will help to obtain the stable internal reference genes under biotic and abiotic conditions and lay the foundation for in-depth functional research of target genes or genomics on olfactory molecular mechanisms, temperature adaptability, and insecticide resistance in X. quadripes.

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.

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.

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.

Screening Potential Reference Genes in Tuta absoluta with Real-Time Quantitative PCR Analysis under Different Experimental Conditions

Tuta absoluta is one of the most significant invasive pests affecting tomato plants worldwide. RT-qPCR has emerged as one of the most sensitive and accurate methods for detecting gene expression data. The screening of stable internal reference genes is the most critical step for studying the molecular mechanisms of environmental adaptability. The stable reference genes expressed in T. absoluta under specific experimental conditions have not yet been clarified. In this study, seven candidate reference genes (RPL27, RPS13, RPS15, EF1-α, TUB, TBP, and β-actin) and their optimal numbers were evaluated under biotic (developmental stages and adult tissues) and abiotic (insecticide, temperature, and plant VOC) conditions using four software programs. Our results identified the following reference genes and numbers as optimal: three genes (EF1-α, RPS13, and RPL27) for different developmental stages (egg, larva, pupa, unmated adult), two genes (RPS13 and TBP) for adult tissues (antenna, head, thorax, abdomen, leg), two genes (TBP and RPS13) for insecticides (Bacillus thuringiensis, chlorpyrifos, abamectin-aminomethyl, and chlorantraniliprole), two genes (RPL27 and TUB) for temperature-induced stresses (0, 25, and 40 °C), and two genes (RPS13 and TUB) for VOC-induced stresses (nonanal, α-phellandrene, and tomato leaves). Our results provide a reference for selecting appropriate reference genes for further study of the functional genes of T. absoluta under different experimental conditions.