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.

Understanding the blood-feeding patterns of mosquitoes is essential for evaluating their potential as disease vectors, especially in urban areas where mosquitoes coexist with humans, domestic animals and wildlife. This study aimed to bridge a substantial gap in regional knowledge by identifying the blood meal sources of field-collected mosquitoes in domestic and open green environments from two urbanisations of temperate Argentina, the Área Metropolitana de Buenos Aires (AMBA) and Tandil, using molecular techniques. Female mosquitoes were collected from November 2019 to March 2020 and April–May 2021. A bipartite network analysis was performed for each environment and urbanisation. A total of 103 blood meals from Aedes (2 species) and Culex (7 species) were identified. Among these, five mammal and 18 bird species were recognised as hosts. Aedes mosquitoes exclusively fed on mammals, while Culex mosquitoes exhibited a broader host range including both birds and mammals. In AMBA, the open green environments were composed by more mosquito species than the domestic environments, while both presented similar numbers of vertebrate species. In contrast, in open green environments from Tandil only blood-fed Aedes albifasciatus were collected. For open green environments of AMBA and domestic environments of Tandil, results suggested some degree of host selection. For the three main vectors of diseases in the region, Aedes aegypti, Ae. albifasciatus and Culex pipiens molestus, we present the first molecular evidence of human blood meals in South America. Epidemiological significance of the present findings is discussed.

Reticulitermes chinensis Snyder is an important pest in forestry and construction and is widely distributed in China. We found that Serratia marcescens Bizio strain SM1 has insecticidal activity to R. chinensis, but the pathogenic mechanism of SM1 to R. chinensis is not clear. Therefore, full-length transcriptome sequencing was performed on R. chinensis infected with SM1 and the control group. A total of 230 differentially expressed genes were identified by comparing SM1 infection group and the control group, among which 103 were downregulated and 127 were upregulated. We found downregulated genes in nine metabolic pathway categories, among which carbohydrate metabolism had the most downregulated genes, followed by energy metabolism and amino acid metabolism. We also found that some downregulated genes were related to pattern recognition receptors, cellular immunity, and humoral immunity, indicating that R. chinensis immunity was negatively affected by SM1 infection. In addition, some genes in signal transduction and genetic information processing pathways were downregulated. In this study, high-throughput full-length transcriptome analysis was used to analyse the pathogenic mechanism of SM1 to R. chinensis. The results of this study provide useful information for exploring the relationship between SM1 and R. chinensis, and provide theoretical support for the future application of SM1 and the prevention and treatment of R. chinensis.

Puparia are commonly found in tsetse fly larviposition sites during studies on larval ecology. This chitinous shell is representative of past or ongoing exploitation of these sites by tsetse flies. The morphological characteristics of the puparium are not sufficiently distinctive to allow identification of the species. This study explores the applicability of biomolecular techniques on empty puparia for tsetse fly species identification. Five techniques were compared for DNA extraction from tsetse fly puparia, 1/Chelex® 100 Resin, 2/CTAB, 3/Livak's protocol, 4/DEB + proteinase K and 5/QIAamp® DNA Mini kit, using two homogenisation methods (manual and automated). Using a combination of two primer pairs, Chelex, CTAB, and DEB + K proved the most efficient on fresh puparia with 90, 85, and 70% samples identified, respectively. Shifting from fresh to one- to nine-month-old puparia, the Chelex method gave the best result allowing species identification on puparia up to seven months old. The subsequent testing of the Chelex extraction protocol identified 152 (60%) of 252 field-collected puparia samples at species level. The results show that reliable genetic identification of tsetse flies species can be performed from empty puparia, what can prove of great interest for future ecological studies on larviposition sites. The Chelex technique was the most efficient for DNA extraction, though the age-limit of the samples stood at seven months, beyond which DNA degradation probably compromises the genetic analysis.

The fall armyworm (FAW), Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae) is a highly destructive polyphagous pest that primarily damages maize. Maize is considered a most versatile crop for growing intercrops due to the wide row it needs. Maize–pea intercropping is preferred by small and marginal farmers worldwide due to various advantages including higher yield and improved economic benefits. However, the success of this intercropping system may be hampered if pea could sustain the FAW population. Thus, to clarify the fitness and potential effect of S. frugiperda on pea, we analysed the survival and development of S. frugiperda fed on pea leaves in the laboratory and constructed age-stage and two-sex life tables. Results showed that FAW successfully completed its life cycle when fed on pea and produced fertile offspring. The pre-adult duration was significantly higher on pea than maize. The net reproductive rate, intrinsic and finite rate of population increase on pea (135.06 offspring per individual, 0.12 offspring per individual per day and 1.13 times per day) were all significantly different from those on maize (417.64 offspring per individual, 0.19 offspring per individual per day and 1.21 times per day). The probability of survival of S. frugiperda at each stage was lower when fed on pea leaves than that of maize-fed larvae. Due to the overlapping growth periods of the maize and pea, S. frugiperda can easily proliferate throughout the year by shifting between adjacent crops. Thus, this study revealed the adaptability of S. frugiperda for pea and provides the foundation for further assessment of FAW risk to other inter-crops.

The eggplant fruit and shoot borer (EFSB) is a devastating pest of eggplants (Solanum aethiopicum L. and Solanum melongena L.) in Ghana, causing significant economic losses. Although initially thought to be the Leucinodes orbonalis Guenee species found in Asia, recent European and Mediterranean Plant Protection Organization reports suggest its absence in Africa. However, eight Leucinodes species have been recently described in Africa, including two new species, Leucinodes africensis sp. n. and Leucinodes laisalis Walker, which were intercepted in eggplant fruits exported from Ghana to the United Kingdom. Despite the reported absence of L. orbonalis in Africa, it remains on the pest list of Ghana as a species known to attack eggplants. To accurately determine the identity of the EFSB complex occurring on eggplant in Southern Ghana, molecular and morphological taxonomic tools were employed, and adult male populations were monitored in on-farm conditions. Our results revealed the presence of two EFSB species, L. africensis and L. laisalis, in the shoot and fruits of eggplants, with L. africensis being the dominant species and widely distributed in Southern Ghana. Notably, L. africensis males were attracted to the pheromone lure of L. orbonalis despite the two species being biologically distinct. This study provides crucial information on correctly identifying the EFSB species attacking eggplants in Southern Ghana and has significant implications for developing management interventions against these pests and their effects on international eggplant trade.

Hemiptera insects exhibit a close relationship to plants and demonstrate a diverse range of dietary preferences, encompassing phytophagy as the predominant feeding habit while a minority engages in carnivorous or haematophagous behaviour. To counteract the challenges posed by phytophagous insects, plants have developed an array of toxic compounds, causing significant evolutionary selection pressure on these insects. In this study, we employed a comparative genomics approach to analyse the expansion and contraction of gene families specific to phytophagous insect lineages, along with their adaptive evolutionary traits, utilising representative species from the Hemiptera order. Our investigation revealed substantial expansions of gene families within the phytophagous lineages, especially in the Pentatomomorpha branch represented by Oncopeltus fasciatus and Riptortus pedestris. Notably, these expansions of gene families encoding enzymes are potentially involved in hemipteran-plant interactions. Moreover, the adaptive evolutionary analysis of these lineages revealed a higher prevalence of adaptively evolved genes in the Pentatomomorpha branch. The observed branch-specific gene expansions and adaptive evolution likely contribute significantly to the diversification of species within Hemiptera. These results help enhance our understanding of the genomic characteristics of the evolution of different feeding habits in hemipteran insects.

Drosophila suzukii is a pest of global concern due to its great impact on several crops. The entomopathogenic nematode Steinernema carpocapsae was highly virulent to the larvae of the fly although some immune mechanisms were triggered along the infection course. Thus, to understand the gene activation profile we performed a comparative transcriptome of D. suzukii larvae infected with S. carpocapsae and Xenorhabdus nematophila to map the differentially expressed genes involved in the defence response. The analysis exposed the induction of genes involved in the humoral response such as the antimicrobial peptides and pattern-recognition receptors while there was a suppression of the cellular defence. Besides, genes involved in melanisation, and clot formation were downregulated hindering the encapsulation response and wound healing. After the infection, larvae were in a stress condition with an enrichment of metabolic and transport functionalities. Concerning the stress response, we observed variations of the heat-shock proteins, detoxification, and peroxidase enzymes. These findings set a genetical comprehensive knowledge of the host-pathogen relation of D. suzukii challenged with S. carpocapsae which could support further comparative studies with entomopathogenic nematodes.

Paranosema locustae is an environmentally friendly parasitic predator with promising applications in locust control. In this study, transcriptome sequencing was conducted on gonadal tissues of Locusta migratoria males and females infected and uninfected with P. locustae at different developmental stages. A total of 18,635 differentially expressed genes (DEGs) were identified in female ovary tissue transcriptomes, with the highest number of DEGs observed at 1 day post-eclosion (7141). In male testis tissue transcriptomes, a total of 32,954 DEGs were identified, with the highest number observed at 9 days post-eclosion (11,245). Venn analysis revealed 25 common DEGs among female groups and 205 common DEGs among male groups. Gene ontology and Kyoto Encyclopaedia of Genes and Genome analyses indicated that DEGs were mainly enriched in basic metabolism such as amino acid metabolism, carbohydrate metabolism, lipid metabolism, and immune response processes. Protein–protein interaction analysis results indicated that L. migratoria regulates the expression of immune- and reproductive-related genes to meet the body's demands in different developmental stages after P. locustae infection. Immune- and reproductive-related genes in L. migratoria gonadal tissue were screened based on database annotation information and relevant literature. Genes such as Tsf, Hex1, Apolp-III, Serpin, Defense, Hsp70, Hsp90, JHBP, JHE, JHEH1, JHAMT, and VgR play important roles in the balance between immune response and reproduction in gonadal tissues. For transcriptome validation, Tsf, Hex1, and ApoLp-III were selected and verified by quantitative real-time polymerase chain reaction (qRT-PCR). Correlation analysis revealed that the qRT-PCR expression patterns were consistent with the RNA-Seq results. These findings contribute to further understanding the interaction mechanisms between locusts and P. locustae.

Aphis spiraecola Patch is one of the most economically important tree fruit pests worldwide. The pyrethroid insecticide lambda-cyhalothrin is commonly used to control A. spiraecola. In this 2-year study, we quantified the resistance level of A. spiraecola to lambda-cyhalothrin in different regions of the Shaanxi province, China. The results showed that A. spiraecola had reached extremely high resistance levels with a 174-fold resistance ratio (RR) found in the Xunyi region. In addition, we compared the enzymatic activity and expression level of P450 genes among eight A. spiraecola populations. The P450 activity of A. spiraecola was significantly increased in five regions (Xunyi, Liquan, Fengxiang, Luochuan, and Xinping) compared to susceptible strain (SS). The expression levels of CYP6CY7, CYP6CY14, CYP6CY22, P4504C1-like, P4506a13, CYP4CZ1, CYP380C47, and CYP4CJ2 genes were significantly increased under lambda-cyhalothrin treatment and in the resistant field populations. A L1014F mutation in the sodium channel gene was found and the mutation rate was positively correlated with the LC50 of lambda-cyhalothrin. In conclusion, the levels of lambda-cyhalothrin resistance of A. spiraecola field populations were associated with P450s and L1014F mutations. Our combined findings provide evidence on the resistance mechanism of A. spiraecola to lambda-cyhalothrin and give a theoretical basis for rational and effective control of this pest species.

Deformed wing virus (DWV), which is an important honey bee virus transmitted by Varroa destructor (V. destructor), causes colony losses in honey bee colonies. This study aimed to investigate the prevalence and genetic diversity of DWV in honey bees in Türkiye and to determine the role of V. destructor in the transmission of the genetic variants of DWV. Honey bee samples were collected from 62 apiaries, by simple random sampling, during March 2022 and April 2023. The presence of V. destructor in collected bee samples was examined using a stereo microscope. Real-time RT-PCR was used for the detection of DWV-A and DWV-B (Varroa destructor virus-1 (VDV-1)) viruses. Genetic characterisation of the positive samples was conducted by sequencing polyprotein genomic region. Considering the V. destructor infestation rate of 3% as relevant, out of the 62 apiaries examined, 17 (27.4%) were positive. However, DWV-A and VDV-1 specific RNA was not detected in V. destructor samples. VDV-1 specific RNA was detected in 6.5% (4/62) of the apiaries, whereas DWV-A was not detected in the sampled apiaries. Phylogenetic analysis showed that isolates detected in this study were located in a separate cluster from previously characterised DWV-A and VDV-1 isolates. According to RDP4 and GARD analyses, DWV-VDV-1 recombination breakpoints were detected in field isolates. To the best our knowledge, this is the first report of the presence of VDV-1-DWV recombinants in Türkiye. Further studies are needed to determine the impact of VDV-1-DWV recombinants and their virological and antigenic properties.

Lateralisation is a well-established phenomenon observed in an increasing number of insect species. This study aims to obtain basic details on lateralisation in courtship and mating behaviour in Ostrinia furnacalis, the Asian corn borer. We conducted laboratory investigations to observe lateralisation in courtship and mating behaviours in adult O. furnacalis. Our goal was also to detect lateralised mating behaviour variations during sexual interactions and to elucidate how these variances might influence the mating success of males. Our findings reveal two distinct lateralised traits: male approaches from the right or left side of the female and the direction of male turning displays. Specifically, males approaching females from their right side predominantly exhibited left-biased 180° turning displays, while males approaching females from the left-side primarily displayed right-biased 180° turning displays. Notably, left-biased males, executing a 180° turn for end-to-end genital contact, initiated copulation with fewer attempts and began copulation earlier than their right-biased approaches with left-biased 180° turning displays. Furthermore, mating success was higher when males subsequently approached the right side of females during sexual encounters. Left-biased 180° turning males exhibited a higher number of successful mating interactions. These observations provide the first report on lateralisation in the reproductive behaviour of O. furnacalis under controlled laboratory conditions and hold promise for establishing reliable benchmarks for assessing and monitoring the quality of mass-produced individuals in pest control efforts.

Cadra cautella is a serious insect pest of stored figs and dates. The irradiation sensitivity of different development stages of C. cautella and large-scale testing of the proposed irradiation quarantine doses (50–500 Gy), were investigated. The impact of a PI dose of 400 Gy on the physiochemical and microbiological quality of dry dates (Bartamoda cv.) stored at room temperature was also investigated. An irradiation dose of 100 Gy prevented egg hatching in the F1 generation when 1–3 days old eggs were irradiated. Irradiation doses of 200 and 300 Gy prevented adult emergence when 2nd and 4th instar larvae were irradiated. When the pupae stage was irradiated, an irradiation dose of 400 Gy prevented the hatchability of F1 generation, indicating that this stage was the most radio-tolerant. The results of large-scale testing of the proposed phytosanitary irradiation dose (400 Gy) applied to 18, 0000 pupae resulted in no reproduction (zero hatching of F1 generation). There were no significant differences in the physiochemical properties of stored dates during the storage period at room temperature. Stable ESR signal intensity was recorded for 6 months in all parts of the irradiated fruits, and the intensity was highest in the kernel. The PI dose of 400 Gy also slightly reduced all microorganisms' counts. In conclusion, the dose level of 400 Gy stopped the reproduction potential of C. cautella. and they maintained the quality characteristics of dry date Bartamoda fruits during storage at room temperature for 6 months in tightly closed packages.

Changes in the distribution of species due to global climate change have a critically significant impact on the increase in the spread of invasive species. An in-depth study of the distribution patterns of invasive species and the factors influencing them can help to better predict and combat invasive alien species. Rhynchophorus ferrugineus Olivier is an invasive species that primarily harms plants of Trachycarpus H. Wendl. The pest invades trees in three main ways: by laying eggs and incubating them in the crown of the plant, on roots at the surface and at the base of the trunk or petiole. Most of the plants in the genus Trachycarpus are taller, and the damage is concentrated in the middle and upper parts of the plant, making control more difficult. In this paper, we combine 19 bioclimatic variables based on the MaxEnt model to project the current and future distributions of R. ferrugineus under three typical emission scenarios (2.6 W m−2 (SSP1-2.6), 4.5 W m−2 (SSP2-4.5) and 8.5 W m−2 (SSP5-8.5)) in the 2050s and 2090s. Among the 19 bioclimatic variables, five variables were screened out by contribution rates, namely annual mean temperature (BIO 1), precipitation of driest quarter (BIO 17), minimum temperature of coldest month (BIO 6), mean diurnal range (BIO 2) and precipitation of wettest quarter (BIO 16). These five variables are key environmental variables that influence habitat suitability for R. ferrugineus and are representative in reflecting its potential habitat. The results showed that R. ferrugineus is now widely distributed in the southeastern coastal area of China (high suitability zone), concentrating in the provinces of Hainan, Guangdong, Fujian, Guangxi and Taiwan. In the future, the area of high and low suitability zones will increase and the area of medium suitability zones will decrease. The area of low suitability zone will still be in the largest proportion. This study aims to provide a theoretical reference for the future control of R. ferrugineus from the perspective of geographic distribution.

The Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Psyllidae), is the key vector insect transmitting the Candidatus Liberibacter asiaticus (CLas) bacterium that causes the devastating citrus greening disease (Huanglongbing, HLB) worldwide. The D. citri salivary glands (SG) exhibit an important barrier against the transmission of HLB pathogen. However, knowledge on the molecular mechanism of SG defence against CLas infection is still limited. In the present study, we compared the SG transcriptomic response of CLas-free and CLas-infected D. citri using an illumine paired-end RNA sequencing. In total of 861 differentially expressed genes (DEGs) in the SG upon CLas infection, including 202 upregulated DEGs and 659 downregulated DEGs were identified. Functional annotation analysis showed that most of the DEGs were associated with cellular processes, metabolic processes, and the immune response. Gene ontology and Kyoto Encyclopaedia of Genes and Genomes enrichment analyses revealed that these DEGs were enriched in pathways involving carbohydrate metabolism, amino acid metabolism, the immune system, the digestive system, the lysosome, and endocytosis. A total of 16 DEGs were randomly selected to further validate the accuracy of RNA-Seq dataset by reverse-transcription quantitative polymerase chain reaction. This study provides substantial transcriptomic information regarding the SG of D. citri in response to CLas infection, which may shed light on the molecular interaction between D. citri and CLas, and provides new ideas for the prevention and control of citrus psyllid.

Spodoptera littoralis (Boisd) (Lepidoptera: Noctuidae) is a highly polyphagous insect that significantly reduces agricultural production of several food staples. We evaluated performance of S. littoralis on several meridic diets based on various maize hybrids, including Oteel, Simon, Valbum, SC703, and SC704. Growth, feeding behaviours, and activity of digestive enzymes of S. littoralis were examined under laboratory conditions. In addition, selected biochemical characteristics of maize hybrid seeds were evaluated, including starch, protein, anthocyanin, as well as phenolic and flavonoid contents, to examine relationships between plant properties and digestive performance of S. littoralis. Performance of S. littoralis on maize hybrids, as measured by nutritional indices, was related to both proteolytic and amylolytic activities quantified using gut extracts. Larval S. littoralis reared on SC703 exhibited the highest efficiency of conversion of digested food, while the lowest was recorded in those fed on the Oteel hybrid. S. littoralis reared on SC703 and Oteel also exhibited the highest and lowest relative growth rates, respectively. The highest levels of proteolytic activity in S. littoralis were measured from larvae reared on the SC703 hybrid, while the lowest levels occurred on the Oteel and Valbum hybrids. Amylolytic activity was lowest in larvae reared on SC703 and Valbum hybrids and highest in larvae reared on the Oteel hybrid. Our results suggest that the SC703 hybrid was the most suitable host for S. littoralis, while the Oteel hybrid demonstrated the greatest level of tolerance against S. littoralis of those evaluated. We discuss the potential utility of maize hybrids exhibiting tolerance traits against this cosmopolitan pest with reference to cultivation of tolerant varieties and identification of specific tolerance traits.

Many parasitoids alter their reproductive behaviour in response to the quality of encountered hosts. They make adaptive decisions concerning whether to parasitise a potential host, the number of eggs laid on an accepted host, and the allocation of sex to their offspring. Here we present evidence that Goniozus jacintae Farrugia (Hymenoptera: Bethylidae), a gregarious ectoparasitoid of larval tortricids, adjusts its reproductive response to the size and developmental stage of larvae of the light brown apple moth (LBAM), Epiphyas postvittana (Walker) (Lepidoptera: Tortricidae). Goniozus jacintae parasitises instars 3–6 of LBAM, but most readily parasitises the later, larger, instars. Brood sizes were bigger on larger hosts and brood sex ratios were female biased (proportion of males = 0.23) with extremely low variance (never >1 male in a brood at emergence), perhaps the most precise of all studied bethylids. Host size did not influence brood development time, which averaged 19.64 days, or the body size of male offspring. However, the size of females was positively correlated with host size and negatively correlated with brood size. The sizes of individual males and females were positively related to the average amount of host resource available to individuals within each brood, suggesting that adult body size is affected by scramble competition among feeding larvae. Average brood sizes were: 3rd instar host, 1.3 (SE ± 0.075); 4th instar, 2.8 (SE ± 0.18); 5th instar, 4.7 (SE ± 0.23); 6th instar, 5.4 (SE ± 0.28). The largest brood size observed was 8 individuals (7 females, 1 male) on the 6th instar of LBAM. These results suggest that later instars would give the highest yield to optimise mass-rearing of G. jacintae if used for augmentative biological pest control.

Since its introduction in Spain in 2004, Aedes albopictus has rapidly spread across the country. Its aggressive biting behaviour causes nuisance, limiting outdoor activities. Also, its role as a vector of several arboviruses implies a major public health risk, with several cases of autochthonous dengue having been reported nationwide over the past few years. Control strategies usually focus on interventions in breeding sites. As such, accurate knowledge of the main larval habitats becomes a major priority in infested areas. A detailed identification of breeding sites of Ae. albopictus was carried out in the outdoors of 60 residential properties during July–August 2022 in El Vedat de Torrent (Valencia, Eastern Spain), an area recently colonised by this species. A total of 1444 real and potential breeding sites were examined. The most abundant potential larval habitat were plant pot plates (6.48 units/house), although a low infestation level was found, both for larvae (2.06% positivity, x̄ = 30.5 larvae/container), and pupae (0.51%, x̄ = 2.5 pupae/container). A total of 7715 larvae and 205 pupae were found in a disused flooded water pool depuration system. Animal drinkers, buckets and irrigation water containers were found to be the most common positive containers. No statistical difference was observed among the different container materials. A general statistical increase of 1 larva per 11.7 ml of water in breeding sites was detected. Breeding sites of other species such as Culex pipiens (n = 2) and Culex modestus (n = 1) were also rarely found in this residential area. To our knowledge, this is the first aedic index study carried out in Europe, and it provides valuable information about the main domestic breeding habitats of Ae. albopictus, which can greatly improve control programmes.

The coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Curculionidae), is a major destructive insect pest of coffee, which impacts the coffee crops negatively. As a draft genome has been completed for this insect, most molecular studies on gene transcriptional levels under different experimental conditions will be conducted using real-time reverse-transcription quantitative polymerase chain reactions (RT-qPCR). However, the lack of suitable internal reference genes will affect the accuracy of RT-qPCR results. In this study, the expression stability of nine candidate reference genes was evaluated under different developmental stages, temperature stress, and Beauveria bassiana infection. Data analyses were completed by four commonly used programs, BestKeeper, NormFinder, geNorm, and RefFinder. The result showed that RPL3 and EF1α combination were recommended as the most stable reference genes for developmental stages. EF1α and RPS3a combination were the top two stable reference genes for B. bassiana infection. RPS3a and RPL3 combination performed as the optimal reference genes both in temperature stress and all samples. Our results should provide a good foundation for the expression profile analyses of target genes in the future, especially for molecular studies on insect genetic development, temperature adaptability, and immune mechanism to entomogenous fungi in H. hampei.