Exposure to low temperatures during honeybee development has been shown to impede brain development and affect cognitive function in adult bees. On the other hand, neuronal damage due to oxidative stress has been reported in many cases. Hence, biochemical parameters related to oxidative stress in honeybee pupae brain were determined. The levels of GSH in the pupal brain decreased after 24 h and 48 h of exposure to low temperatures; there were also reduced activities of SOD and CAT enzymes following 48 h of low-temperature treatment compared to the control group. Furthermore, analysis of transcriptome data post-24 h and -48 h low-temperature stress revealed the suppression of the glutathione metabolism and peroxisome pathways in pupal brains. Additionally, expression pattern clustering analysis and KEGG enrichment showed that 10 differentially expressed genes with down-regulated expression trends post-low-temperature treatment were significantly enriched in the peroxisome pathway, including PEX10, highlighting their connection to peroxisome function. RT-qPCR validation was conducted on 11 core enriched genes in pathways identified via GSEA, and all these genes exhibited a downregulated expression pattern, confirming the inhibition of glutathione metabolism and peroxisome function under low-temperature stress. The present study showed that exposing honeybee pupae to low temperatures suppressed both the glutathione metabolism and peroxisome pathways, resulting in increased oxidative stress. This research enhances our understanding of how the pupal brain reacts to cold stress and illuminates the neural damage associated with low temperatures during honeybee capped brood development.

Thermal condition affects the development and growth of ectotherms. The stenothermic honeybee brood, particularly the prepupae, are sensitive to low rearing temperature. The fat body plays important roles in energy reserve and metabolism during the honeybee brood development. To date, the fat body metabolic changes in prepupae responding to cold stress have not been completely understood. In this study, the ultra-performance liquid chromatography coupled with tandem mass spectrometry (UPLC-MS/MS)-based non-target metabolome was analyzed between the cold-treated (CT, 20 °C, 36 h) and control (CK, 35 °C) fat body in prepupal honeybees. The fat body metabolomic data showed that the levels of 1860 and 254 metabolites were significantly increased and decreased, respectively, in cold-stressed prepupae. These altered metabolites, glutamine, glutamic acid, pyroglutamic acid, and oxidized glutathione, were significantly enriched into glutamine metabolism and glutathione metabolism pathways. Furthermore, the expression levels of glutamine metabolism-related genes, glutaminase (GLS), glutamate dehydrogenase (GDH), and gamma-glutamyl transferase (GGT-1 and GGT-7), were significantly decreased in cold-exposed prepupae compared with the control groups. Meanwhile, the oxidized glutathione (GSSG), but not the reduced glutathione (GSH) content, was increased in the cold-exposed group compared with controls. Collectively, our data revealed the fat body metabolomic changes in larva-to-pupa transition when exposed to cold stress. Our data provided new insights into stenothermic honeybee sensitivity to cold, characterized by perturbation of glutamine metabolism and oxidative stress.

Cuticular proteins (CPs)—key components of the insect exoskeleton—not only regulate development but also serve as structural barriers that enhance resistance against environmental stressors. This study identified CP gene families in Apis mellifera and analyzed their expression patterns during the worker capped brood development stages from mature larva to pre-eclosion. Using a comprehensive genome-wide bioinformatic approach, we identified 85 CP genes in A. mellifera which comprise six families: CPR (n = 43), CPAPs (n = 27), CPF (n = 2), Tweedle (n = 2), CPLCP (n = 8) and Apidermin (n = 3). Analysis of CP gene evolutionary relationship revealed that each CP family forms a distinct, relatively independent clade. Domain and motif analyses confirmed that all CPR members harbor a conserved Chitin_Bind_4 domain, consistent with CPR family structures in other taxa. Additionally, CPAP members possess one or three Chitin-binding Peritrophin-A domain (CBM_14), CPF members possess a conserved Pupal cuticle protein C1 domain (Cuticle_3), and Tweedle members contain a conserved domain of unknown function (DUF243). In addition, the analysis found no conserved domain within the CPLCP and Apidermin families. RNA-seq data revealed dynamic expression patterns of AmCPs during pupal development, with each gene family displaying a relatively characteristic temporal profile. Quantitative PCR validation of eight highly expressed CPR genes at 9 days post-capping confirmed the RNA-seq results. This work provides a comprehensive bioinformatic characterization and transcriptional analysis of CP genes in A. mellifera, offering a foundation for future functional studies on cuticle formation and identifying candidate genes potentially involved in cuticle development in honeybees. This work relies on transcriptomic data and in silico analyses. All proposed biological roles are hypothetical and require experimental validation.

Melanin is an important component of the body color of honeybees, and its formation changes with the age of a capped brood of bees. However, up to now, the regulatory mechanism of melanin formation in honeybees remains unclear. To analyze the differential expression profile of microRNAs (miRNAs) in worker bees of Apis mellifera and to reveal the regulatory roles of differentially expressed miRNAs (DEmiRNAs) and mRNAs in the formation process of melanin during the capped brood stage, we used sRNA-seq technology and related software to analyze samples from four key developmental stages during the capped brood stage, when body color develops in Apis mellifera, namely, mature larvae (L0), pre-pupae (PP3), early pupae (P6) and mid-pupae (P9). A total of 1291 miRNAs were identified by bioinformatics. Three comparison groups were analyzed: L0 vs. PP3, PP3 vs. P6, and P6 vs. P9. A total of 171, 94, and 19 DEmiRNAs were identified in these groups, respectively, which regulate 1481, 690, and 182 differentially expressed target mRNAs (target DEmRNAs). The functional analysis of target DEmRNAs indicated that DEmiRNAs might regulate the formation of capped brood melanin in honeybees by activating expression changes in key genes in signaling pathways, such as the Wnt signaling pathway, melanogenesis, and the Toll and Imd signaling pathway, through activating miR-315-x, miR-8, ple, yellow family genes, wnt1, etc. Our research provides a theoretical basis for future analysis of the regulatory role of miRNAs in the formation of melanin in honeybees.