Reference Gene Selection for Insect Expression Studies Using Quantitative Real-Time PCR: The Head of the Honeybee, Apis mellifera, After a Bacterial Challenge

Uncovering hidden dangers: The combined toxicity of abamectin and lambda-cyhalothrin on honey bees

Studying the toxic effects of pesticides on bees has consistently been a prominent area of interest for researchers. Nonetheless, existing research has predominantly concentrated on individual toxicity assessments, leaving a gap in our understanding of mixed toxicity. This study delves into the individual and combined toxic effects of abamectin (ABA) and lambda-cyhalothrin (LCY) on honey bees (Apis mellifera) in laboratory settings. We discovered that ABA (96 h-LC50 value of 0.079 mg/L) exhibited greater acute toxicity to honey bees compared to LCY (96 h-LC50 value of 9.177 mg/L). Moreover, the mixture of ABA and LCY presented an acute antagonistic effect on honey bees. Additionally, our results indicated that exposure to LCY, at medium concentration, led to a reduction in the abundance of gut core bacterium Snodgrassella. However, an increase in the abundance of Bifidobacterium was noted when exposed to a medium concentration of LCY and its mixture with ABA. Transcriptomic analysis revealed significant regulation of certain genes in the medium concentration of all three treatments compared to the control group, primarily enriching in metabolism and immune-related pathways. Following chronic exposure to field-relevant concentrations of ABA, LCY, and their mixture, there were significant alterations in the activities of immunity-related enzyme polyphenol oxidase (PPO) and detoxification enzymes glutathione S-transferase (GST) and carboxylesterase (CarE). Additionally, the expression of four genes (abaecin, cyp9e2, cyp302a1, and GstD1) associated with immune and detoxification metabolism was significantly altered. These findings suggest a potential health risk posed by the insecticides ABA and LCY to honey bees. Despite exhibiting acute antagonistic effect, mixed exposure still induced damage to bees at all levels. This study advances our knowledge of the potential adverse effects of individual or combined exposure to these two pesticides on non-target pollinators and offers crucial guidance for the use of insecticides in agricultural production.

Determination of internal controls for quantitative gene expression of Spodoptera litura under microbial pesticide stress

Quantitative real-time polymerase chain reaction (qRT-PCR) has become a commonly used method for the quantification of gene expression. However, accurate qRT-PCR analysis requires a valid internal reference for data normalization. To determine the valid reference characterized with low expression variability among Spodoptera litura samples after microbial pesticide treatments, nine housekeeping genes, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), arginine kinase, ubiquitin C, actin-5C (ACT5C), actin, ribosomal protein S13 (RPS13), tubulin, acidic ribosomal protein P0 (RPLP0) and ubiquinol-cytochrome c reductase, were evaluated for their suitability using geNorm, Normfinder, BestKeeper, RefFinder and the comparative delta CT methods in this study. S. litura larvae after direct treatment (larvae were immersed in biopesticides), indirect treatment (larvae were fed with biopesticide immersed artificial diets) and comprehensive treatment (larvae were treated with the first two treatments in sequence), respectively with Metarhizium anisopliae, Empedobacter brevis and Bacillus thuringiensis, were investigated. The results indicated that the best sets of internal references were as follows: RPLP0 and ACT5C for direct treatment conditions; RPLP0 and RPS13 for indirect treatment conditions; RPS13 and GAPDH for comprehensive treatment conditions; RPS13 and RPLP0 for all the samples. These results provide valuable bases for further genetic researches in S. litura.

A longitudinal field study of commercial honey bees shows that non-native probiotics do not rescue antibiotic treatment, and are generally not beneficial

Probiotics are widely used in agriculture including commercial beekeeping, but there is little evidence supporting their effectiveness. Antibiotic treatments can greatly distort the gut microbiome, reducing its protective abilities and facilitating the growth of antibiotic resistant pathogens. Commercial beekeepers regularly apply antibiotics to combat bacterial infections, often followed by an application of non-native probiotics advertised to ease the impact of antibiotic-induced gut dysbiosis. We tested whether probiotics affect the gut microbiome or disease prevalence, or rescue the negative effects of antibiotic induced gut dysbiosis. We found no difference in the gut microbiome or disease markers by probiotic application or antibiotic recovery associated with probiotic treatment. A colony-level application of the antibiotics oxytetracycline and tylosin produced an immediate decrease in gut microbiome size, and over the longer-term, very different and persistent dysbiotic effects on the composition and membership of the hindgut microbiome. Our results demonstrate the lack of probiotic effect or antibiotic rescue, detail the duration and character of dysbiotic states resulting from different antibiotics, and highlight the importance of the gut microbiome for honeybee health.

Efficient RNA interference method during caste differentiation with hormone treatment in the termite Reticulitermes speratus (Isoptera: Rhinotermitidae)

Unveiling the proximate mechanism of caste differentiation is crucial for understanding insect social evolution, and gene function analysis is an important tool in this endeavor. The RNA interference (RNAi) technique is useful in termites, but its knockdown effects may differ among species. One of the most important model species in the field of termite sociogenomics is Reticulitermes speratus Kolbe (Isoptera: Rhinotermitidae). Presoldier and worker differentiation of this species can be artificially induced by juvenile hormone and 20-hydroxyecdysone application, respectively. However, appropriate RNAi technique of genes expressed during caste differentiation has never been considered. To clarify this issue, first, we injected nine different volumes of nuclease-free water (NFW, 0-404.8 nL) into workers and found that survival and caste differentiation rates were strongly reduced by the application of the top three largest volumes. Second, we injected double-stranded (ds) RNA of ecdysone receptor homolog (RsEcR) (2.0 µg/151.8 nL NFW) into workers with hormone treatments. The expression levels of RsEcR were significantly reduced at 9 days after dsRNA injection. RsEcR RNAi strongly affected both molting events during presoldier and worker differentiation induced by hormone treatments. The present results highlight the need for caution regarding injection volumes for RNAi experiments using hormone treatments. We suggest that the injection of dsRNA solution (2 µg; approximately 100-200 nL) is suitable for RNAi experiments during caste differentiation induced by hormone application in R. speratus.

Expression of Molecular Markers of Resilience against Varroa destructor and Bee Viruses in Ethiopian Honey Bees (Apis mellifera simensis) Focussing on Olfactory Sensing and the RNA Interference Machinery

Varroa destructor mites and the viruses it vectors are two major factors leading to high losses of honey bees (Apis mellifera) colonies worldwide. However, honey bees in some African countries show resilience to varroa infestation and/or virus infections, although little is known about the mechanisms underlying this resilience. In this study, we investigated the expression profiles of some key molecular markers involved in olfactory sensing and RNA interference, as these processes may contribute to the bees' resilience to varroa infestation and virus infection, respectively. We found significantly higher gene expression of the odorant binding protein, OBP14, in the antennae of Ethiopian bees compared to Belgian bees. This result suggests the potential of OBP14 as a molecular marker of resilience to mite infestation. Scanning electron microscopy showed no significant differences in the antennal sensilla occurrence and distribution, suggesting that resilience arises from molecular processes rather than morphological adaptations. In addition, seven RNAi genes were upregulated in the Ethiopian honey bees and three of them-Dicer-Drosha, Argonaute 2, and TRBP2-were positively correlated with the viral load. We can conclude that the antiviral immune response was triggered when bees were experiencing severe viral infection and that this might contribute to the bees' resilience to viruses.

RNAi-mediated silencing of AccCYP6k1 revealed its role in the metabolic detoxification of Apis cerana cerana

Insect cytochrome P450 monooxygenases (P450s or CYPs) perform important functions in the metabolic detoxification of both endogenous and exogenous substrates. However, the mechanism of action of the P450 genes in bees is unclear. In this study, we investigated the effects of AccCYP6k1 on the metabolism and detoxification of Apis cerana cerana. Spatiotemporal expression profiling revealed that the expression of AccCYP6k1 was the highest in foragers (A15) and was mainly expressed in the leg, midgut and head. RT-qPCR results showed that AccCYP6k1 exhibited different expression patterns following exposure to xenobiotics. In addition, silencing AccCYP6k1 increased the pesticides sensitivity and affected the detoxification system and antioxidant process of A. cerana cerana. In brief, the induced expression of AccCYP6k1 is related to the resistance of A. cerana cerana, while knockdown AccCYP6k1 affect the pesticides resistance and metabolic detoxification system of A. cerana cerana. These findings not only support the theoretical basis of metabolic detoxification in bees but also provide a better understanding of P450-mediated resistance to pesticides in insects.

Non-optimal ambient temperatures aggravate insecticide toxicity and affect honey bees Apis mellifera L. gene regulation

In this study, we conducted a transcriptional analysis of five honey bee genes to examine their functional involvement vis-à-vis ambient temperatures and exposure to imidacloprid. In a 15-day cage experiment, three cohorts of one-day-old sister bees emerged in incubators, were distributed into cages, and maintained at three different temperatures (26 °C, 32 °C, 38 °C). Each cohort was fed a protein patty and three concentrations of imidacloprid-tainted sugar (0 ppb, 5 ppb and 20 ppb) ad libitum. Honey bee mortality, syrup and patty consumption were monitored daily over 15 days. Bees were sampled every three days for a total of five time points. RT-qPCR was used to longitudinally assess gene regulation of Vg, mrjp1, Rsod, AChE-2 and Trx-1 using RNA extracted from whole bee bodies. Kaplan-Meier models show that bees kept at both non-optimal temperatures (26 °C and 38 °C) were more susceptible to imidacloprid, with significantly higher mortality (P < 0.001 and P < 0.01, respectively) compared to the control. At 32 °C, no differences in mortality (P = 0.3) were recorded among treatments. In both imidacloprid treatment groups and the control, the expression of Vg and mrjp1 was significantly downregulated at 26 °C and 38 °C compared to the optimal temperature of 32 °C, indicating major influence of ambient temperature on the regulation of these genes. Within the ambient temperature groups, both imidacloprid treatments exclusively downregulated Vg and mrjp1 at 26 °C. AChE-2 and the poorly characterized Rsod gene were both consistently upregulated at the highest temperature (38 °C) compared to the ideal temperature (32 °C) in all treatment groups. Trx-1 showed no effect to both temperature and imidacloprid treatments and was regulated in an age-related manner. Overall, our results indicate that ambient temperatures amplify imidacloprid toxicity and affect honey bee gene regulation.

Characterisation of the heat shock protein Tid and its involvement in stress response regulation in Apis cerana

Objective: The impact of various environmental stresses on native Apis cerana cerana fitness has attracted intense attention in China. However, the defence responses of A. cerana cerana to different stressors are poorly understood. Here, we aimed to elucidate the regulatory mechanism mediated by the tumorous imaginal discs (Tid) protein of A. cerana cerana (AccTid) in response to stressors. Methods: We used some bioinformatics softwares to analyse the characterisation of Tid. Then, qRT-PCR, RNA interference and heat resistance detection assays were used to explore the function of Tid in stress response in A. cerana cerana. Results: AccTid is a homologous gene of human Tid1 and Drosophila Tid56, contains a conserved J domain and belongs to the heat shock protein DnaJA subfamily. The level of AccTid induced expression was increased under temperature increases from 40°C to 43°C and 46°C, and AccTid knockdown decreased the heat resistance of A. cerana cerana, indicating that the upregulation of AccTid plays an important role when A. cerana cerana is exposed to heat stress. Interestingly, contrary to the results of heat stress treatment, the transcriptional level of AccTid was inhibited by cold, H₂O₂ and some agrochemical stresses and showed no significant change under ultraviolet ray and sodium arsenite stress. These results suggested that the requirement of A. cerana cerana for Tid differs markedly under different stress conditions. In addition, knockdown of AccTid increased the mRNA levels of some Hsps and antioxidant genes. The upregulation of these Hsps and antioxidant genes may be a functional complement of AccTid knockdown. Conclusion: AccTid plays a crucial role in A. cerana cerana stress responses and may mediate oxidative damage caused by various stresses. Our findings will offer fundamental knowledge for further investigations of the defence mechanism of A. cerana cerana against environmental stresses.

Evaluation of Reference Genes for Real-Time Quantitative PCR Analysis in Tissues from Bumble Bees (Bombus Terrestris) of Different Lines

Bumble bees are important alternative pollinators and model insects due to their highly developed sociality and colony management. In order to better understand their molecular mechanisms, studies focusing on the genetic and molecular aspects of their development and behavior are needed. Although quantitative real-time polymerase chain reaction (qRT-PCR) can be used to quantify the relative expression of target genes, internal reference genes (which are stably expressed across different lines and tissues) must first be identified to ensure the accurate normalization of target genes. In order to contribute to molecular studies on bumble bees, we used Bombus terrestris to determine the expression stability of eight reference genes (β-actin (ACT), Arginine Kinase (AK), Phospholipase A2 (PLA2), Elongation factor 1 alpha (EF-1), Ribosomal proteins (S5, S18, S28) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH)) in five different lines and several tissues (ovary, thorax, fat body, and head) using RT-qPCR procedures and four analysis programs (RefFinder, NormFinder, BestKeeper, and geNorm). In general, the S28, S5, and S18 ribosomal protein genes and the PLA2 and EF-1 genes showed the highest stability and were therefore identified as suitable reference genes for the bumble bee species and their defined lines and tissues. Our results also emphasized the need to evaluate the stability of candidate reference genes for any differently designed lines and tissue conditions in bumble bee species.

Validation of quantitative real-time PCR reference genes and spatial expression profiles of detoxication-related genes under pesticide induction in honey bee, Apis mellifera

Recently, pesticides have been suggested to be one of the factors responsible for the large-scale decline in honey bee populations, including colony collapse disorder. The identification of the genes that respond to pesticide exposure based on their expression is essential for understanding the xenobiotic detoxification metabolism in honey bees. For the accurate determination of target gene expression by quantitative real-time PCR, the expression stability of reference genes should be validated in honey bees exposed to various pesticides. Therefore, in this study, to select the optimal reference genes, we analyzed the amplification efficiencies of five candidate reference genes (RPS5, RPS18, GAPDH, ARF1, and RAD1a) and their expression stability values using four programs (geNorm, NormFinder, BestKeeper, and RefFinder) across samples of five body parts (head, thorax, gut, fat body, and carcass) from honey bees exposed to seven pesticides (acetamiprid, imidacloprid, flupyradifurone, fenitrothion, carbaryl, amitraz, and bifenthrin). Among these five candidate genes, a combination of RAD1a and RPS18 was suggested for target gene normalization. Subsequently, expression levels of six genes (AChE1, CYP9Q1, CYP9Q2, CYP9Q3, CAT, and SOD1) were normalized with a combination of RAD1a and RPS18 in the different body parts from honey bees exposed to pesticides. Among the six genes in the five body parts, the expression of SOD1 in the head, fat body, and carcass was significantly induced by six pesticides. In addition, among seven pesticides, flupyradifurone statistically induced expression levels of five genes in the fat body.

Identification of the cuticle protein AccCPR2 gene in Apis cerana cerana and its response to environmental stress

Cuticular proteins (CPs) are known to play important roles in insect development and defence responses. The loss of CP genes can lead to changes in insect morphology and sensitivity to the external environment. In this study, we identified the AccCPR2 gene, which belongs to the CPR family (including the R&R consensus motif) of CPs, and explored its function in the response of Apis cerana cerana to adverse external stresses. Our results demonstrated that AccCPR2 was highly expressed in the late pupal stage and epidermis, and the expression of AccCPR2 may be induced or inhibited under different stressors. RNA interference experiments showed that knockdown of AccCPR2 reduced the activity of antioxidant enzymes, led to the accumulation of oxidative damage and suppressed the expression of several antioxidant genes. In addition, knockdown of AccCPR2 also reduced the pesticide resistance of A. cerana cerana. The overexpression of AccCPR2 in a prokaryotic system further confirmed its role in resistance to various stresses. In summary, AccCPR2 may play pivotal roles in the normal development and environmental stress response of A. cerana cerana. This study also enriched the theoretical knowledge of the resistance biology of bees.

Activating transcription factor 2 (AccATF2) regulates tolerance to oxidative stress in Apis cerana cerana

Activating transcription factor 2 (ATF2), a basic leucine zipper (bZIP) transcription factor, plays a crucial role in immune and DNA damage response in mammals. However, the function of ATF2 in insects remains unknown. Here, we isolated the ATF2 gene from Apis cerana cerana (AccATF2) and found that AccATF2 was a main regulator of the honeybee response to oxidative stress. Our results showed that AccATF2 was highly expressed in the head, thorax and integument. AccATF2 was expressed throughout the development period of honeybees, and the highest AccATF2 transcript level was noted in brown-eyed pupae, indicating its indispensable roles in honeybee survival. Antioxidant function analysis showed that AccATF2 expression was markedly induced in response to oxidative stress caused by various environmental stresses. AccATF2 overexpression substantially enhanced the tolerance to oxidative stress of Escherichia coli cells compared with control cells. AccATF2 knockdown significantly increased the production of malondialdehyde (MDA), the transcription of antioxidant genes and the activity of antioxidant enzymes in honeybees, suggesting that AccATF2 knockdown resulted in oxidative damage to honeybees. Moreover, AccATF2 knockdown decreased honeybee resistance to oxidative stress caused by high temperature. Overall, AccATF2 plays an important role in maintaining redox homeostasis and protecting honeybees from oxidative stress caused by various environmental stimuli. Our discoveries add to a growing understanding of how honeybees cope with various adverse environmental conditions to ensure their survival.

Worker bees (Apis mellifera) deprived of pollen in the first week of adulthood exhibit signs of premature aging

Pollinator populations, including bees, are in rapid decline in many parts of the world, raising concerns over the future of ecosystems and food production. Among the factors involved in these declines, poor nutrition deserves attention. The diet consumed by adult worker honeybees (Apis mellifera) is crucial for their behavioral maturation, i.e., the progressive division of labor they perform, such as nurse bees initially and later in life as foragers. Poor pollen nutrition is known to reduce the workers' lifespan, but the underlying physiological and genetic mechanisms are not fully understood. Here we investigate how the lack of pollen in the diet of workers during their first week of adult life can affect age-related phenotypes. During the first seven days of adult life, newly emerged workers were fed either a pollen-deprived (PD) diet mimicking that of an older bee, or a control pollen-rich (PR) diet, as typically consumed by young bees. The PD-fed bees showed alterations in their fat body transcriptome, such as a switch from a protein-lipid based metabolism to a carbohydrate-based metabolism, and a reduced expression of genes involved with immune response. The absence of pollen in the diet also led to an accumulation of oxidative stress markers in fat body tissue and alterations in the cuticular hydrocarbon profiles, which became similar to those of chronologically older bees. Together, our data indicate that the absence of pollen during first week of adulthood triggers the premature onset of an aging-related worker phenotype.

Antibiotics Alter the Expression of Genes Related to Behavioral Development in Honey Bees (Hymenoptera: Apidae)

Honey bees, as many species of social insects, display a division of labor among colony members based on behavioral specializations related to age. Adult worker honey bees perform a series of tasks in the hive when they are young (such as brood care or nursing) and at ca. 2-3 wk of age, shift to foraging for nectar and pollen outside the hive. The transition to foraging involves changes in metabolism and neuroendocrine activities. These changes are associated with a suite of developmental genes. It was recently demonstrated that antibiotics influence behavioral development by accelerating or delaying the onset of foraging depending on timing of antibiotic exposure. To understand the mechanisms of these changes, we conducted a study on the effects of antibiotics on expression of candidate genes known to regulate behavioral development. We demonstrate a delay in the typical changes in gene expression over the lifetime of the individuals that were exposed to antibiotics during immature stage and adulthood. Additionally, we show an acceleration in the typical changes in gene expression on individuals that were expose to antibiotics only during immature stage. These results show that timing of antibiotic exposure alter the typical regulation of behavioral development by metabolic and neuroendocrine processes.

Distinct molecular impact patterns of abamectin on Apis mellifera ligustica and Apis cerana cerana

The effects of insecticides on bee health are a topic of intensive research. Although abamectin is toxic to bees, the molecular impact of abamectin needs to be clarified. Here, we found that Apis cerana cerana exhibited a higher mortality rate when exposed to abamectin than Apis mellifera ligustica. In addition, A. cerana cerana had markedly higher numbers of differentially expressed genes (DEGs), differentially expressed proteins (DEPs) and differentially expressed metabolites (DEMs) than A. mellifera ligustica during exposure to abamectin. These results indicate that abamectin exposure exerts stronger effects on A. cerana cerana than on A. mellifera ligustica. In addition, six DEGs, two DEPs and two DEMs overlapped between the two bee species under abamectin exposure; however, some genes or proteins from the zinc finger protein, superoxide dismutase and peroxiredoxin families and the energy metabolism pathway were only unregulated in A. cerana cerana, which indicates a significant difference in the impact of abamectin on the two bee species. Despite these differences, several of the same gene families, such as heat shock proteins, cytochrome P450, odorant-binding proteins and cuticle proteins, and pathways, including the carbohydrate metabolism, immune system, lipid metabolism, amino acid metabolism, sensory system, locomotion and development pathways, were influenced by abamectin exposure in both A. cerana cerana and A. mellifera ligustica. Together, our results indicate that abamectin causes adverse effects on bees and thus poses a risk to bee populations and that abamectin exposure affects A. cerana cerana more strongly than A. mellifera ligustica. These findings improve our understanding of the behavioural and physiological effects of abamectin on bees.

Roles of the protein disulphide isomerases AccPDIA1 and AccPDIA3 in response to oxidant stress in Apis cerana cerana

Protein disulphide isomerase (PDI) plays an important role in a variety of physiological processes through its oxidoreductase activity and molecular chaperone activity. In this study, we cloned two PDI family members, AccPDIA1 and AccPDIA3, from Apis cerana cerana. AccPDIA1 and AccPDIA3 had typical sequence features of PDI family members and were constitutively expressed in A. cerana cerana. The expression levels of AccPDIA1 and AccPDIA3 were generally upregulated after treatment with a variety of environmental stress factors. Inhibition assays showed that E. coli expressing recombinant AccPDIA1 and AccPDIA3 proteins was more resistant to oxidative stress than control E. coli. In addition, silencing AccPDIA1 or AccPDIA3 in A. cerana cerana resulted in significant changes in the expression levels of several antioxidant-related genes as well as the enzymatic activities of peroxidase (POD), superoxide dismutase (SOD) and catalase (CAT) and reduced the survival rate of A. cerana cerana under oxidative stress caused by high temperature. In conclusion, our results suggest that AccPDIA1 and AccPDIA3 may play important roles in the antioxidant activities of A. cerana cerana.

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.

Gut microbiome drives individual memory variation in bumblebees

The potential of the gut microbiome as a driver of individual cognitive differences in natural populations of animals remains unexplored. Here, using metagenomic sequencing of individual bumblebee hindguts, we find a positive correlation between the abundance of Lactobacillus Firm-5 cluster and memory retention on a visual discrimination task. Supplementation with the Firm-5 species Lactobacillus apis, but not other non-Firm-5 bacterial species, enhances bees' memory. Untargeted metabolomics after L. apis supplementation show increased LPA (14:0) glycerophospholipid in the haemolymph. Oral administration of the LPA increases long-term memory significantly. Based on our findings and metagenomic/metabolomic analyses, we propose a molecular pathway for this gut-brain interaction. Our results provide insights into proximate and ultimate causes of cognitive differences in natural bumblebee populations.

Identification and characterization of an Apis cerana cerana nucleoside diphosphate kinase (AccNDPK) associated with oxidative stress

Nucleoside diphosphate kinases (NDPKs) are widespread nucleotide-metabolizing enzymes that are involved in a variety of biological processes, including responses to oxidative stress. Although studies have been conducted on NDPKs in mammals and some plants, there is scant research on insect NDPKs, especially in honey bees. In the present study, we isolated AccNDPK from Apis cerana cerana. Sequence analysis showed that AccNDPK has high homology with many NDPKs and contains a highly conserved NDPK active site motif. Based on phylogenetic analysis, AccNDPK has a relatively recent evolutionary relationship with NDPKs in other hymenopteran insects. AccNDPK was found to be highly expressed in newly emerged honey bees and muscle tissues, and RT-qPCR analysis and bacteriostatic assays showed that the expression level of AccNDPK is affected by abnormal temperature, UV light, H₂O₂, heavy metals, and various pesticides. After AccNDPK knockdown, antioxidant-related genes, including AccCAT, AccCYP4G11, AccGSTS4, AccTpx1 and AccMsrA, were upregulated, whereas AccGSTD, AccGST1, AccHSP22.6 and AccTrx1 were downregulated. Furthermore, catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) activities were significantly increased, and the tolerance of bees to oxidative stress caused by cyhalothrin was reduced by silencing of AccNDPK. Given these findings, we speculate that AccNDPK plays an important role in the oxidative stress response of A. cerana cerana.

A potential probiotic Leuconostoc mesenteroides TBE-8 for honey bee

An isolated bacterium TBE-8, was identified as Leuconostoc mesenteroides according to the sequences of 16S rDNA and the 16S-23S rDNA intergenic spacer region. The probiotic properties of the L. mesenteroides TBE-8 strain were characterized and revealed that TBE-8 could utilize various carbohydrates, exhibited high tolerance to sucrose's osmotic pressure and acidic conditions, and could mitigate the impact of the bee pathogen Paenibacillus larvae. In addition, we found that the TBE-8 broth increased the expression of the nutrition-related genes major royal jelly protein 1 and vitellogenin in bees by approximately 1400- and 20-fold, respectively. The expression of genes encoding two antibacterial peptides, hymenoptaecin and apidaecin, in the bee abdomen was significantly increased by 17- and 7-fold in bees fed with the TBE-8 fermented broth. Furthermore, we fed four-frame bee colonies with 50% sucrose syrup containing TBE-8 and can detect the presence of approximately 2 × 106 16S rDNA copies of TBE-8 in the guts of all bees in 24 h, and the retention of TBE-8 in the bee gut for at least 5 days. These findings indicate that the L. mesenteroides TBE-8 has high potential as a bee probiotic and could enhance the health of bee colonies.

Vairimorpha (Nosema) ceranae Infection Alters Honey Bee Microbiota Composition and Sustains the Survival of Adult Honey Bees

Simple Summary

The gut microbiota, in addition to the hosts and the pathogens, has become the third factor involved in gut disease developments, including honey bees. Interestingly, various studies reported positive associations between the gut bacteria and the most commonly found microsporidian pathogen instead of negative associations. To investigate the positive associations, a prebiotic that also exists in honey was added in the trials. Bees fed the prebiotics have slightly higher pathogen counts but lower mortalities. Microbiota analyses suggested that bees with the infection have a microbiota composition similar to that of bees with a longer lifespan, and the prebiotic seemed to enhance the similarities. Since microsporidia typically cause chronic infections, the positive associations may serve to sustain the host lifespans which is the optimal outcome for the pathogen that the survived bees can withstand pathogen proliferation and transmit the pathogens. Although the mechanisms underlying the associations were not revealed, this study indicated that nosema disease management in bees through changes in microbiota may shorten the lifespans or enhance both the infection and the bee population. Such results have appeared in recent field studies. More studies will be needed for the disease management using bee gut microbiota.

Abstract

Vairimorpha (Nosema) ceranae is the most common eukaryotic gut pathogen in honey bees. Infection is typically chronic but may result in mortality. Gut microbiota is a factor that was recently noted for gut infectious disease development. Interestingly, studies identified positive, instead of negative, associations between core bacteria of honey bee microbiota and V. ceranae infection. To investigate the effects of the positive associations, we added isomaltooligosaccharide (IMO), a prebiotic sugar also found in honey, to enhance the positive associations, and we then investigated the infection and the gut microbiota alterations using qPCR and 16S rRNA gene sequencing. We found that infected bees fed IMO had significantly higher V. ceranae spore counts but lower mortalities. In microbiota comparisons, V. ceranae infections alone significantly enhanced the overall microbiota population in the honey bee hindgut and feces; all monitored core bacteria significantly increased in the quantities but not all in the population ratios. The microbiota alterations caused by the infection were enhanced with IMO, and these alterations were similar to the differences found in bees that naturally have longer lifespans. Although our results did not clarify the causations of the positive associations between the infections and microbiota, the associations seemed to sustain the host survival and benefit the pathogen. Enhancing indigenous gut microbe to control nosema disease may result in an increment of bee populations but not the control of the pathogen. This interaction between the pathogen and microbiota potentially enhances disease transmission and avoids the social immune responses that diseased bees die prematurely to curb the disease from spreading within colonies.

Expression stability of candidate RT-qPCR housekeeping genes in Spodoptera frugiperda (Lepidoptera: Noctuidae)

Reverse-transcription quantitative polymerase chain reaction (RT-qPCR) is commonly used to quantify gene expression. For normalization, the expression of each gene is compared with a reference "housekeeping" gene that is stably expressed under relevant stress. Unfortunately, there have been no reports on the stability of such reference genes under various treatments of the Spodoptera frugiperda. In this study, we used five tools (RefFinder, GeNorm, NormFinder, BestKeeper, and ΔCt methods) to evaluate the stability of 12 candidate reference genes (RPS18, β-tubulin, GAPDH, RPS7, RPS15, RPL7, RPL32, Actin-5C, EF1-α, EF1-γ, RPL27, and ACE) in different instars, tissues, and treatments (high and low temperature, UV-A, and emamectin benzoate). Several ribosomal proteins (RPS7, RPS15, RPL32, RPS18, and RPL7), GAPDH, Actin-5C, and β-tubulin, were relatively stable, suggesting that they are ideal housekeeping genes for various treatments. ACE was extremely unstable under various experimental treatments, rendering it unsuitable as an internal reference. This study identified the reference housekeeping genes stably expressed by S. frugiperda under different treatments, thus setting a foundation for further exploration of the physiological and biochemical mechanisms.

Pollen nutrition fosters honeybee tolerance to pesticides

A reduction in floral resource abundance and diversity is generally observed in agro-ecosystems, along with widespread exposure to pesticides. Therefore, a better understanding on how the availability and quality of pollen diets can modulate honeybee sensitivity to pesticides is required. For that purpose, we evaluated the toxicity of acute exposure and chronic exposures to field realistic and higher concentrations of azoxystrobin (fungicide) and sulfoxaflor (insecticide) in honeybees provided with pollen diets of differing qualities (named S and BQ pollens). We found that pollen intake reduced the toxicity of the acute doses of pesticides. Contrary to azoxystrobin, chronic exposures to sulfoxaflor increased by 1.5- to 12-fold bee mortality, which was reduced by pollen intake. Most importantly, the risk of death upon exposure to a high concentration of sulfoxaflor was significantly lower for the S pollen diet when compared with the BQ pollen diet. This reduced pesticide toxicity was associated with a higher gene expression of vitellogenin, a glycoprotein that promotes bee longevity, a faster sulfoxaflor metabolization and a lower concentration of the phytochemical p-coumaric acid, known to upregulate detoxification enzymes. Thus, our study revealed that pollen quality can influence the ability of bees to metabolize pesticides and withstand their detrimental effects, providing another strong argument for the restoration of suitable foraging habitat.

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.

Transcriptome-level assessment of the impact of deformed wing virus on honey bee larvae

Deformed wing virus (DWV) prevalence is high in honey bee (Apis mellifera) populations. The virus infects honey bees through vertical and horizontal transmission, leading to behavioural changes, wing deformity, and early mortality. To better understand the impacts of viral infection in the larval stage of honey bees, artificially reared honey bee larvae were infected with DWV (1.55 × 10¹⁰ copies/per larva). No significant mortality occurred in infected honey bee larvae, while the survival rates decreased significantly at the pupal stage. Examination of DWV replication revealed that viral replication began at 2 days post inoculation (d.p.i.), increased dramatically to 4 d.p.i., and then continuously increased in the pupal stage. To better understand the impact of DWV on the larval stage, DWV-infected and control groups were subjected to transcriptomic analysis at 4 d.p.i. Two hundred fifty-five differentially expressed genes (DEGs) (fold change ≥ 2 or ≤ -2) were identified. Of these DEGs, 168 genes were downregulated, and 87 genes were upregulated. Gene Ontology (GO) analysis showed that 141 DEGs (55.3%) were categorized into molecular functions, cellular components and biological processes. One hundred eleven genes (38 upregulated and 73 downregulated) were annotated by KO (KEGG Orthology) pathway mapping and involved metabolic pathways, biosynthesis of secondary metabolites and glycine, serine and threonine metabolism pathways. Validation of DEGs was performed, and the related gene expression levels showed a similar tendency to the DEG predictions at 4 d.p.i.; cell wall integrity and stress response component 1 (wsc1), cuticular protein and myo-inositol 2-dehydrogenase (iolG) were significantly upregulated, and small conductance calcium-activated potassium channel protein (SK) was significantly downregulated at 4 d.p.i. Related gene expression levels at different d.p.i. revealed that these DEGs were significantly regulated from the larval stage to the pupal stage, indicating the potential impacts of gene expression levels from the larval to the pupal stages. Taken together, DWV infection in the honey bee larval stage potentially influences the gene expression levels from larvae to pupae and reduces the survival rate of the pupal stage. This information emphasizes the consequences of DWV prevalence in honey bee larvae for apiculture.

Identification of chemosensory genes by antennal transcriptome analysis and expression profiles of odorant-binding proteins in parasitoid wasp Aulacocentrum confusum

The endoparasitoid wasp, Aulacocentrum confusum (Hymenoptera: Braconidae), is a preponderant natural enemy of the larvae of Glyphodes pyloalis Walker (Lepidoptera: Pyralidae), which is a destructive pest of mulberry trees. We first constructed the antennal transcriptome database of A. confusum. In total, we obtained 48,262,304 clean reads from the dataset and assembled 24,324 unigenes. A total of 12,690 (52.17%) unigenes indicated significant similarity (E-value < 10^-5) compared to known protein sequences of other species from the NCBI non-redundant protein database. Gene ontology (GO) and cluster of orthologous groups (COG) analyses were used to determine the functional categories of these genes. A total of 84 putative chemosensory genes were identified from the antennal transcriptome of A. confusum, including 11 putative odorant-binding protein (OBP) genes, six chemosensory protein (CSP) genes, 44 olfactory receptor (OR) genes (including one olfactory co-receptor, Orco), 19 ionotropic receptor (IR) genes, and four sensory neuron membrane protein (SNMP) genes. Results of qPCR assays indicated that among of 11 AconOBPs, nine AconOBP genes were significantly expressed in the antennae of A. confusum adults. AconOBP8 was significantly expressed in the abdomen and AconOBP10 was highly expressed in the thorax. These findings can build a basis for further study on the processes of chemosensory perception in A. confusum at the molecular level.

Golden orb-weaving spider (Trichonephila clavipes) silk genes with sex-biased expression and atypical architectures

Spider silks are renowned for their high-performance mechanical properties. Contributing to these properties are proteins encoded by the spidroin (spider fibroin) gene family. Spidroins have been discovered mostly through cDNA studies of females based on the presence of conserved terminal regions and a repetitive central region. Recently, genome sequencing of the golden orb-web weaver, Trichonephila clavipes, provided a complete picture of spidroin diversity. Here, we refine the annotation of T. clavipes spidroin genes including the reclassification of some as non-spidroins. We rename these non-spidroins as spidroin-like (SpL) genes because they have repetitive sequences and amino acid compositions like spidroins, but entirely lack the archetypal terminal domains of spidroins. Insight into the function of these spidroin and SpL genes was then examined through tissue- and sex-specific gene expression studies. Using qPCR, we show that some silk genes are upregulated in male silk glands compared to females, despite males producing less silk in general. We also find that an enigmatic spidroin that lacks a spidroin C-terminal domain is highly expressed in silk glands, suggesting that spidroins could assemble into fibers without a canonical terminal region. Further, we show that two SpL genes are expressed in silk glands, with one gene highly evolutionarily conserved across species, providing evidence that particular SpL genes are important to silk production. Together, these findings challenge long-standing paradigms regarding the evolutionary and functional significance of the proteins and conserved motifs essential for producing spider silks.

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.

A Potential Fungal Probiotic Aureobasidium melanogenum CK-CsC for the Western Honey Bee, Apis mellifera

Aureobasidium melanogenum has been used as an animal feed additive for improving thehealth of pets, however, it has not yet been applied in honey bees. Here, a fungal strain CK-CsC isolated from bee bread pollen, was identified as A. melanogenum. Following characterizing CK-CsC fermentation broth, the 4-days fermentation broth (SYM medium or bee pollen) of the CK-CsC was used to feed newly emerged adult honey bees in cages under laboratory-controlled conditions for analysis of survival, gene expression of nutrient and antibacterial peptide, and gut microbiota of honey bees. It was found that the CK-CsC fermentation broth (SYM medium or bee pollen) is nontoxic to honey bees, and can regularly increase nutrient gene expression of honey bees. However, significant mortality of bees was observed after bees were fed on the supernatant liquid of the fermentation broth. Notably, this mortality can be lowered by the simultaneous consumption of bee pollen. The honey bees that were fed bee pollen exhibited more γ-Proteobacteria, Bacteriodetes, and Actinobacteria in their gut flora than did the honey bees fed only crude supernatant liquid extract. These findings indicate that A. melanogenum CK-CsC has high potential as a bee probiotic when it was fermented with bee pollen.

Genes regulating development and behavior exhibited altered expression in Drosophila melanogaster exposed to bisphenol A: use of real-time quantitative PCR (qRT-PCR) and droplet digital PCR (ddPCR) in genotoxicity study

Toxicity of bisphenol A on morphological and life-history traits of model insect Drosophila melanogaster was reported in our previous work. In the present study, we have analyzed the adversity of bisphenol A on the reproductive behavior of adult and on the expression of selected genes in the larva and adult stage of fruit fly exposed to bisphenol A (0.007 g/2 ml. or 3.5 mg/ml), in addition to determination of LC50 value of bisphenol A in larva and pupal stage. We employed both the quantitative reverse transcriptase PCR and droplet digital PCR for analyzing the expression profile of seven genes namely, decapentaplegic, vestigial, wingless, foraging, insulin-like receptor, doublesex, and fruitless. We found bisphenol A has more adverse effects on male sexual behavior than females. Moreover, we observed significant downregulation of all the selected genes in treated larvae except, fruitless in male where it showed significant upregulation. On contrary among the treated adult flies, significant downregulation of all target genes in both sexes is evident, except, doublesex and fruitless in males which showed significant upregulation. We did not observe any deviation of male: female sex ratio from 1:1 under bisphenol A exposure. All these results suggest bisphenol A adversely affects the optimum functioning of genes which are involved in the regulation of metabolic pathways, behavioral pattern, stress response, endocrine homeostasis, neural functioning, and the development of the specific organ in Drosophila melanogaster. Our result not only provides a foundation to study further the bisphenol A toxicity on different pivotal genes in Drosophila but also suggests the use of the droplet digital PCR technology in toxicity measurement at the molecular level in eukaryotic model systems.

Thiamethoxam exposure deregulates short ORF gene expression in the honey bee and compromises immune response to bacteria

Maximizing crop yields relies on the use of agrochemicals to control insect pests. One of the most widely used classes of insecticides are neonicotinoids that interfere with signalling of the neurotransmitter acetylcholine, but these can also disrupt crop-pollination services provided by bees. Here, we analysed whether chronic low dose long-term exposure to the neonicotinoid thiamethoxam alters gene expression and alternative splicing in brains of Africanized honey bees, Apis mellifera, as adaptation to altered neuronal signalling. We find differentially regulated genes that show concentration-dependent responses to thiamethoxam, but no changes in alternative splicing. Most differentially expressed genes have no annotated function but encode short Open Reading Frames, a characteristic feature of anti-microbial peptides. As this suggested that immune responses may be compromised by thiamethoxam exposure, we tested the impact of thiamethoxam on bee immunity by injecting bacteria. We show that intrinsically sub-lethal thiamethoxam exposure makes bees more vulnerable to normally non-pathogenic bacteria. Our findings imply a synergistic mechanism for the observed bee population declines that concern agriculturists, conservation ecologists and the public.

Variation in Expression of Reference Genes across Life Stages of a Bee, Megachile rotundata

Simple Summary

Reference genes are key to normalizing expression data across samples of organisms collected after different treatments are applied, but often, reference genes are not properly validated for this purpose. In this report, we screened several genes for a solitary bee, Megachile rotundata, and identified two (RPS18, and RPL8) with very stable expression levels across all life stages of the bee, and under a variety of environmental conditions, including during and after diapause. These genes should make good reference genes. We also identified other genes with stable expression, even if used only for a limited number of developmental stages. This information is important for future gene expression studies on these bees, but it also demonstrates the importance of validating reference genes in general.

Abstract

The alfalfa leafcutting bee, Megachile rotundata is widely used in the western United States as a pollinator for alfalfa seed production. Unfortunately, immatures experience high mortality in agriculturally managed populations. Quantified gene expression could be used to identify how this bee responds during different life stages to pathogens, environmental toxins, and other stresses, but stably expressed reference genes are needed to normalize transcription data. We evaluated twelve candidate genes for their transcription stability across different life stages, including during and after diapause. RPS18 and RPL8 were the two most stably expressed genes, followed by RPS5 and RPL27A. These genes were also very stable even during and after diapause, while the most variable genes being APN, PMIIM, NPC2, and Cr-PII had increased expression levels during larval growth and were also variable during and after diapause. The four reference genes we identified in M. rotundata may prove useful for transcriptomic studies in other bees as well, such as honey bees.

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.

Transcriptomic and metabolomic landscape of the molecular effects of glyphosate commercial formulation on Apis mellifera ligustica and Apis cerana cerana

Understanding the causes of the decline in bee population has attracted intensive attention worldwide. The indiscriminate use of agrochemicals is a persistent problem due to their physiological and behavioural damage to bees. Glyphosate and its commercial formulation stand out due to their wide use in agricultural areas and non-crop areas, such as parks, railroads, roadsides, industrial sites, and recreational and residential areas, but the mode of action of glyphosate on bees at the molecular level remains largely unelucidated. Here, we found that the numbers of differentially expressed genes and metabolites under glyphosate commercial formulation (GCF) stress were significantly higher in Apis cerana cerana than in Apis mellifera ligustica. Despite these differences, the number of differentially expressed transcripts increased following an increase in the GCF treatment time in both A. cerana cerana and A. mellifera ligustica. GCF exerted adverse impacts on the immune system, digestive system, nervous system, amino acid metabolism, carbohydrate metabolism, growth and development of both bee species by influencing their key genes and metabolites to some extent. The expression of many genes involved in immunity, agrochemical detoxification and resistance, such as antimicrobial peptides, cuticle proteins and cytochrome P450 families, was upregulated by GCF in both bee species. Collectively, our results indicate that both A. cerana cerana and A. mellifera ligustica strive to mitigate the pernicious effects caused by GCF by regulating detoxification and immune systems. Moreover, A. cerana cerana might be better able to withstand the toxic effects of GCF with lower fitness costs than A. mellifera ligustica. Our work will contribute to elucidating the deleterious physiological and behavioural impacts of GCF on bees.

Black queen cell virus and Nosema ceranae coinfection in Africanized honey bees from southern Brazil

ABSTRACT: Bees are fundamental in several aspects, especially in relation to plant biodiversity and pollination. Recently, immense losses are being faced in the number of Brazilian colonies, mainly in southern states of the country, which has a strong beekeeping activity. There are indications that, among the reasons for the losses, pathogens that affect the health of bees may be involved. Among them, the microsporidium Nosema and the black queen cell virus (BQCV) stand out for their prevalence. In this study, 92 colonies of 17 apiaries from southern Brazil were evaluated for infection by Nosema ceranae, Nosema apis and BQCV. Nucleic acid extractions and cDNA synthesis were performed from adult bee samples, followed by Reverse Transcription Polymerase Chain Reaction (RT-PCR) and multiplex PCR. Eight BQCV positive samples were subjected to sequencing. The results showed that N. ceranae and BQCV are circulating in the Southern region of the country, which may be the reason for the loss of colonies. N. apis was not found. N. ceranae was found in 57.6% (53/92) of the colonies and BQCV in 32.6% (30/92). Co-infection was found in 25% (23/92) of the colonies studied, a factor that is suggested to be reducing the hosts’ longevity due to the synergistic action of the pathogens. The samples submitted to sequencing indicated similarity of 96.8 to 100% between them, in addition to strong similarity with sequences from Asia, United States, Germany and Peru. This study reports the circulation of N. ceranae and BQCV in apiaries in southern Brazil, in addition to being the first phylogenetic analysis of the Brazilian BQCV sequence.

Development of Sarcophaga dux (diptera: Sarcophagidae) at constant temperatures and differential gene expression for age estimation of the pupae

Sarcophaga dux (Diptera: Sarcophagidae) is a necrophagous flesh fly species with potential forensic value for estimating minimum postmortem interval (PMI_min). The basic developmental data and precise age estimates of the pupae are significant for PMI_min estimation in forensic investigations. In the present study, we investigated the development data of that species at seven constant temperatures varying from 16 °C to 34 °C, including body length changes of the larve, developmental duration and accumulated degree hours of the preadults. Several reference genes for relative quantification of the differentially expressed genes (DEGs) were firstly selected and evaluated in the pupae of different ages under different temperatures. The DEGs of the insects during the pupal period at different constant temperatures (34, 25 and 16 °C) were further analyzed for more precise age estimation. The results showed that the developmental durations of the preadults at 16, 19, 22, 25, 28, 31 and 34 °C were 1478.6 ± 18.3 h, 726.1 ± 15.8 h, 538.5 ± 0.9 h, 394.1 ± 9.5 h, 375.6 ± 10.8 h, 284.1 ± 7.3 h, and 252.5 ± 6.1 h, respectively. The developmental threshold temperature the flies was 12.27 ± 0.35 °C, and the thermal summation constant was 5341.71 ± 249.29° hours. The most reliable reference genes during the pupal period at different temperatures were found: GST1 and 18S rRNA for the 34 °C group, GST1 and RPL49 for 25 °C, and 18S rRNA and 28S rRNA for 16 °C. The four differential expression genes (Hsp60, A-alpha, ARP, and RPL8) have the potential to be used for more precise age estimation of pupal S. dux. This work provides important basic developmental data and a more precise age estimation method for pupal S. dux, and improves the value of this species for PMI_min estimation in forensic investigations.

Evaluation of reference genes for quantitative real-time PCR normalization in the scarab beetle Holotrichia oblita

Quantitative real-time polymerase chain reaction (qPT-PCR) is commonly used to analyze gene expression, however, the accuracy of the normalized results is affected by the expression stability of reference genes. Holotrichia oblita (Coleoptera: Scarabaeidae) causes serious damage to crops. Reliable reference genes in H. oblita are needed for qRT-PCR analysis. Therefore, we evaluated 13 reference genes under biotic and abiotic conditions. RefFinder provided a comprehensive stability ranking, and geNorm suggested the optimal number of reference genes for normalization. RPL13a and RPL18 were the most suitable reference genes for developmental stages, tissues, and temperature treatments; RPL13a and RPS3 were the most suitable for pesticide and photoperiod treatments; RPS18 and RPL18 were the most suitable for the two sexes. We validated the normalized results using odorant-binding protein genes as target genes in different tissues. Compared with the selected suitable reference genes, the expression of OBP1 in antennae, abdomen, and wings, and OBP2 in antennae and wings were overestimated due to the instability of ACTb. These results identified several reliable reference genes in H. oblita for normalization, and are valuable for future molecular studies.

Selection of Reference Genes for Quantitative Real-Time PCR in Aquatica leii (Coleoptera: Lampyridae) Under Five Different Experimental Conditions

Aquatic fireflies are important indicators of the quality of freshwater environments and key models for research on insect adaptation to freshwater environments. For these investigations, gene expression analyses using quantitative real-time PCR are heavily dependent on reliable reference genes. In this study, based on a transcriptome assembly and annotation for the aquatic firefly Aquatica leii at the adult and larval stages, 10 candidate reference genes (α-tubulin, β-tubulin, β-actin, EF1A, SDHA, UBQ, GST, GAPDH, RPS31, and RPL13A) were identified for analyses of expression stability. Quantitative real-time PCR analyses for each candidate reference genes in A. leii was conducted for four developmental stages, four adult tissue types, two adult sexes, and two ecological stressors [adults exposed to five temperatures and larvae exposed to four concentrations of benzo(a)pyrene]. Results were evaluated by three independent algorithms (geNorm, NormFinder, and BestKeeper) and one comparative algorithm (RefFinder). The expression stability of candidate reference genes in A. leii differed under various conditions. Reference genes with the most stable expressions levels in different tissues, temperatures, sexes, developmental stages, and concentrations of benzo(a)pyrene were α-tubulin, GST, β-actin, β-tubulin, and α-tubulin, respectively. Furthermore, the optimal normalization factors (NFs) for the quantification of the expression levels of target genes by quantitative real-time PCR analyses of A. leii were identified for each experimental group. In particular, NF = 2 for different tissues (α-tubulin + β-tubulin), different sexes (β-actin + EF1A), and larvae exposed to different concentrations of benzo(a)pyrene (α-tubulin + EF1A); NF = 3 for developmental stages (GST + GAPDH + SDHA) and adults exposed to different temperatures (β-tubulin + EFA + GST). In addition, we surveyed the expression profiles of two target genes (CYP3A and CSP8) in larvae exposed to different concentrations of benzo(a)pyrene and in different adult tissues. The results further validated the reliability of the reference genes. The optimal reference genes for various experimental conditions identified in these analyses provide a useful tool for ecological studies of aquatic fireflies.

Tyramine and its receptor TYR1 linked behavior QTL to reproductive physiology in honey bee workers (Apis mellifera)

Honey bees (Apis mellifera) provide an excellent model for studying how complex social behavior evolves and is regulated. Social behavioral traits such as the division of labor have been mapped to specific genomic regions in quantitative trait locus (QTL) studies. However, relating genomic mapping to gene function and regulatory mechanism remains a big challenge for geneticists. In honey bee workers, division of labor is known to be regulated by reproductive physiology, but the genetic basis of this regulation remains unknown. In this case, QTL studies have identified tyramine receptor 1 (TYR1) as a candidate gene in region pln2, which is associated with multiple worker social traits and reproductive anatomy. Tyramine (TA), a neurotransmitter, regulates physiology and behavior in diverse insect species including honey bees. Here, we examine directly the effects of TYR1 and TA on worker reproductive physiology, including ovariole number, ovary function and the production of vitellogenin (VG, an egg yolk precursor). First, we used a pharmacology approach to demonstrate that TA affects ovariole number during worker larval development and increases ovary maturation during the adult stage. Second, we used a gene knockdown approach to show that TYR1 regulates vg transcription in adult workers. Finally, we estimated correlations in gene expression and propose that TYR1 may regulate vg transcription by coordinating hormonal and nutritional signals. Taken together, our results suggest TYR1 and TA play important roles in regulating worker reproductive physiology, which in turn regulates social behavior. Our study exemplifies a successful forward-genetic strategy going from QTL mapping to gene function.

Identification and temporal expression profiles of cuticular proteins in the endoparasitoid wasp, Microplitis mediator

Recently, parasitoid wasp species Microplitis mediator has evoked increasing research attention due to its possible use in the control of Lepidoptera insects. Because insect development involves changes in cuticle composition, identification and expression analysis of M. mediator cuticular proteins may clarify the mechanisms involved in parasite development processes. We found 70 cuticular proteins from the M. mediator transcriptome and divided them into seven distinct families. Expression profiling indicated that most of these cuticular protein genes have expression peaks specific for one particular developmental stage of M. mediator. Eggs and pupae have the highest number of transcriptionally active cuticular protein genes (47 and 52 respectively). Only 12 of these genes maintained high expression activity during late larval development. Functional analysis of two larval proteins, MmCPR3 and MmCPR14, suggested their important role in the proper organization of the cuticle layers of larvae. During M. mediator larval development, normal cuticle formation can be supported by a limited number of cuticular proteins.

Evaluation of the expression stability of reference genes in Apis mellifera under pyrethroid treatment

Honeybees (Apis mellifera L.), which unquestionably play an economically important role in pollination and agricultural production, are at risk of decline. To study changes in gene expression in insects upon exposure to pesticides or other external stimuli, appropriate reference genes are required for data normalization. Since there is no such gene that is absolutely invariable under all experimental conditions, the aim of this study was to identify the most stable targets suitable for subsequent normalization in quantitative experiments based on real-time polymerase chain reaction in honeybee research. Here, we evaluated the expression of fifteen candidate housekeeping genes from three breeding lines of honeybees treated with pyrethroids to identify the most stable genes. The tested insects were exposed to deltamethrin or lambda-cyhalothrin, and then, changes in the accumulation of selected transcripts were assessed, followed by statistical analyses. We concluded that AmRPL32, AmACT and AmRPL13a were the commonly recorded most stable genes in honeybees treated with the selected pyrethroids.

Molecular Assessment of Genes Linked to Immune Response Traits of Honey Bees in Conventional and Organically Managed Apiaries

Simple Summary

Honey bees play a critical role in agriculture as they provide pollination services to many agricultural crops. However, honey bee populations continue to decline due to exposure to pesticides, habitat destruction, pests, diseases and beekeeping practices. In this study, we assessed selected biological parameters associated with honey bee health in two beekeeping practices (organic and conventional). We compared total protein content in young worker bees from organically and conventionally managed apiaries. We also assessed differential gene expression at two levels of Varroa mite infestations (0% and 5%) in selected genes involved in nutrition and cellular defense (vitellogenin (Vg), malvolio (Mvl), prophenoloxidase (PPO)-, genes involved in lifespan (superoxide dismutase (Sod 1), superoxide dismutase 2 (Sod2)) and immune function genes (immune deficiency (Imd), spaetzle (Spz). Total soluble protein in young adult worker bees was similar in both beekeeping practices. The genes PPO, Vg and Mvl were upregulated in young adult bees with no mite infestations from organically managed apiaries and could mount an immune response through Spz and Sod 1 when challenged by 5% Varroa mite infestation levels. Overall, these findings provide useful insights into the genetic response of honey bees under two beekeeping practices and could help improve honey bee health.

Abstract

Honey bees are of great economic importance, not only for honey production but also for crop pollination. However, honey bee populations continue to decline mainly due to exposure to pesticides, pathogens and beekeeping practices. In this study, total soluble protein was measured, total RNA was extracted and first-strand cDNAs were generated. Quantitative PCR was used to assess the relative expression (transcript abundances) of immune function-related genes in honey bees collected from organically and conventionally managed hives. Honey bees collected from conventionally managed hives with 0% Varroa mite infestation levels displayed an upregulated expression of the prophenoloxidase gene (cellular defense). Similarly, honey bees collected from organically managed hives had increased levels of the vitellogenin gene (immune function and longevity). The gene expression for malvolio (sucrose responsiveness) was highest in organically managed hives with 0% Varroa mite infestations. Young adult bees collected from organically managed hives with 5% Varroa mite infestation levels had upregulated expressions of the gene spaetzle, whereas bees from similarly infested, conventionally managed hives did not, suggesting that honey bees from organically managed hives could mount an immune response. In young adult bees collected from organically managed hives only, the expression of the immune deficiency gene (antimicrobial defense) was upregulated. The relative gene expression for superoxide dismutase 1 increased in young adult bees collected from hives with 5% Varroa mite infestation levels as expected. However, for superoxide dismutase 2, there was a high level of gene expression in adult bees from both conventionally managed hives with 0% Varroa mite infestation levels and organically managed hives with 5% Varroa mite infestations. The gene CYP9Q3 (pesticide detoxification) that metabolizes coumaphos and fluvalinate was upregulated in adult bees collected from organically managed bees. Overall, these findings provide useful insights into the genetic response of honey bees to some environmental stressors and could be an important component of best beekeeping practices that intend to enhance honey bee health.

Reference gene selection for qRT-PCR analysis of season- and tissue-specific gene expression profiles in the honey bee Apis mellifera

Honey bees are both important pollinators and model insects due to their highly developed sociality and colony management. To better understand the molecular mechanisms underlying honey bee colony management, it is important to investigate the expression of genes putatively involved in colony physiology. Although quantitative real-time PCR (qRT-PCR) can be used to quantify the relative expression of target genes, internal reference genes (which are stably expressed across different conditions) must first be identified to ensure accurate normalisation of target genes. To identify reliable reference genes in honey bee (Apis mellifera) colonies, therefore, we evaluated seven candidate genes (ACT, EIF, EF1, RPN2, RPS5, RPS18 and GAPDH) in samples collected from three honey bee tissue types (head, thorax and abdomen) across all four seasons using three analysis programmes (NormFinder, BestKeeper and geNorm). Subsequently, we validated various normalisation methods using each of the seven reference genes and a combination of multiple genes by calculating the expression of catalase (CAT). Although the genes ranked as the most stable gene were slightly different on conditions and analysis methods, our results suggest that RPS5, RPS18 and GAPDH represent optimal honey bee reference genes for target gene normalisation in qRT-PCR analysis of various honey bee tissue samples collected across seasons.

Screening and Validation of Reference Genes for RT-qPCR Under Different Honey Bee Viral Infections and dsRNA Treatment

Honey bee viruses are one of the most important pathogens that have contributed to the decrease in honey bee colony health. To analyze the infection dynamics of honey bee viruses, quantification of viral gene expression by RT-qPCR is necessary. However, suitable reference genes have not been reported from viral and RNAi studies of honey bee. Here, we evaluated the expression of 11 common reference genes (ache2, rps18, β-actin, tbp, tif, rpl32, gadph, ubc, α-tubulin, rpl14, and rpsa) from Apis mellifera (Am) and Apis cerana (Ac) under Israeli acute paralysis virus (IAPV), chronic bee paralysis virus (CBPV), and Chinese sacbrood virus (CSBV) infection as well as dsRNA-PGRP-SA treatment, and we confirmed their validation by evaluating the levels of the defensin 1 and prophenoloxidase (ppo) genes during viral infection. Our results showed that the expression of selected genes varied under different viral infections. ache2, rps18, β-actin, tbp, and tif can be used to normalize expression levels in Apis mellifera under IAPV infection, while the combination of actin and tif is suitable for CBPV-infected experiments. The combination of rpl14, tif, rpsa, ubc, and ache2 as well as more reference genes is suitable for CSBV treatment in Apis cerana. Rpl14, tif, rps18, ubc, and α-tubulin were the most stable reference genes under dsRNA treatment in Apis mellifera. Furthermore, the geNorm and NormFinder algorithms showed that tif was the best suitable reference gene for these four treatments. This study screened and validated suitable reference genes for the quantification of viral levels in honey bee, as well as for RNAi experiments.

Role of Apis cerana cerana N-terminal asparagine amidohydrolase (AccNtan1) in oxidative stress

N-Terminal asparagine amidohydrolase is a component of the ubiquitin-dependent N-end rule pathway of protein degradation that has been implicated in a variety of physiological functions, including the sensing of heme, oxygen, nitric oxide, selective elimination of misfolded proteins and the repair of DNA. We identified the Apis cerana cerana N-terminal asparagine amidohydrolase (AccNtan1) gene from A. cerana cerana and investigated its role in oxidation resistance. Multiple sequence alignments and phylogenetic analysis revealed that N-terminal asparagine amidohydrolase is highly conserved in insect species. Quantitative real-time polymerase chain reaction analysis indicated that the expression levels of AccNtan1 were significantly lower in the wing, honey sac and abdomen than in other tissues and were significantly higher in early stages of development, including the larva, prepupa and pink-eyed pupa stages, than in later stages. We further observed that AccNtan1 expression was induced by several environmental stressors, including aberrant temperature, H2O2, UV, heavy metals and pesticides. Moreover, a bacteriostatic assay suggested that overexpression of AccNtan1 enhances the resistance of bacteria to oxidative stress. In addition, knockdown of AccNtan1 using RNA interference significantly affected the expression levels of most antioxidant genes and the activity levels of several antioxidant enzymes. Thus, we hypothesize that AccNtan1 plays important roles in environmental stress responses and antioxidative processes.

Overexpression of Glutathione S-Transferase Genes in Field λ-Cyhalothrin-Resistant Population of Cydia pomonella: Reference Gene Selection and Expression Analysis

Analysis of the glutathione S-transferase (GST) gene expression in an insecticide-resistant strain of Cydia pomonella using real-time quantitative polymerase chain reaction is a key step toward more mechanism studies that require suitable reference genes with stable expression. Here, nine commonly used reference genes were selected, and their expression stabilities were analyzed. Results showed that EF-1α was the most stable reference gene in all of the experimental sets. The combinations of EF-1α and 18S, EF-1α and RPL12, and EF-1α and GAPDH were sufficient for normalization of gene expression analysis accurately in developmental stages, tissues, and larvae exposed to sublethal dose of λ-cyhalothrin, respectively. Additionally, the suitability of particular reference genes was verified by analyzing the spatiotemporal and insecticide-induced expression profiles of CpGSTe3, CpGSTd3, and CpGSTd4, which were overexpressed in a λ-cyhalothrin-resistant population from northeast China. These genes were used to confer the practicability of reference genes chosen in this study.

Analyses of the function of DnaJ family proteins reveal an underlying regulatory mechanism of heat tolerance in honeybee

There is clear evidence of severe honeybee declines in recent years, and parallel declines of plant community and crop productivity that rely on them. Different stresses, including heat stress, are among the primary drivers of this decline. However, the mechanisms by which honeybees respond to heat stress are elusive. Though heat shock proteins (Hsps) play important roles in heat stress response, the function of DnaJs (a subfamily of Hsps) is unclear. Here, we aimed to determine the underlying regulatory mechanism of honeybees to heat stress mediated by DnaJs. We found that several DnaJ genes, including DnaJA1, DnaJB12 and DnaJC8, are key for honeybee heat tolerance. DnaJA1 and DnaJB12 are cytoplasmic proteins, and DnaJC8 is a nuclear protein. The expression of DnaJA1, DnaJB12 and DnaJC8 was induced at different levels under short-term and long-term heat stress. Phenotypic analysis indicated that DnaJA1, DnaJB12 and DnaJC8 knockdown attenuated honeybee heat resistance. In addition, DnaJA1 participated in the heat stress response by upregulating many heat-inducible genes at the transcriptome-wide level, especially LOC108002668 and LOC107995148. Importantly, the upregulation of LOC108002668 and LOC107995148 was significantly repressed under heat stress when DnaJA1 was knocked down. We also found that knockdown of DnaJA1, DnaJB12 and DnaJC8 decreased antioxidant defense ability and increased the degree of oxidative damage in the honeybee. Taken together, our results indicate that DnaJ genes play important roles under heat stress in the honeybee. Overexpression of DnaJ genes may protect honeybees from heat stress-induced injuries and increase their survival rate.