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

Zygentoma occupies a key evolutionary position for understanding the evolution of insect metamorphosis but has received little attention in terms of genetic analysis. To develop functional genomic studies in this insect, we evaluated five candidate internal reference genes for quantitative RT-PCR (qPCR) studies from Thermobia domestica, a representative species of Zygentoma, including Actin 5C (Actin5C), Elongation factor-1 alpha (EF1A), Ribosome protein S26 (RPS26), Ribosome protein L32 (RPL32), and Superoxide dismutase 2 (SOD2), at different developmental stages, in various body parts, and under dsRNA microinjection and starvation stresses, using four algorithms (delta Ct, geNorm, NormFinder and BestKeeper) and a comparative algorithm (RefFinder). Specific suitable reference genes were recommended across specific experimental conditions, and the combination of RPS26 and RPL32 was appropriate for all tested samples. Employing our selected reference gene combination, we investigated the gene expression pattern of Myoglianin (Myo), a crucial gene-regulating insect metamorphosis, in ametabolous T. domestica, and demonstrated the efficiency of RNA interference (RNAi) in firebrat nymphs. This study provides a basis for reliable quantitative studies of genes and greatly benefits evolutionary and functional genomics studies in Zygentoma.