Cytochrome P450 plays a crucial role in regulating insect growth, development, and resisting a variety of stresses. Insect metamorphosis and response to external stress are altered by deleting CYP450 genes. In this study, we identified and analyzed a novel gene of CYP450 family, AccCYP6A13, from Apis cerana cerana, and explored its role in the response of Apis cerana cerana to adverse external stressors. It was found that the expression of AccCYP6A13 was spatiotemporal specificity. The expression level increased with age and reached its highest value in the adult stage. The primarily expressiong location were legs, brain, and epidermis of honeybees. Stress conditions can affect the expression of AccCYP6A13 depending on treatment times. RNA interference experiments have shown that knocking down AccCYP6A13 reduces antioxidant activity and deactivates detoxification enzymes, resulting in oxidative damage accumulation and a decline in detoxification capability in bees, as well as inhibiting numerous antioxidant genes. Additionally, knockdown of the AccCYP6A13 gene in Apis cerana cerana resulted in increased sensitivity to pesticides and increased mortality when treated with neonicotinoid pesticides such as thiamethoxam. AccCYP6A13 overexpression in a prokaryotic system further confirmed its role in resistance to oxidative stress. To summarize, AccCYP6A13 may play an essential role in the normal development and response to environmental stress in Apis cerana cerana. Furthermore, this study contributed to the theoretical understanding of bee resistance biology.

Introduction: Guizhou Province, characterized by complex and diverse geographic and climatic environments, has rich genetic resources for the Chinese honeybee (Apis cerana cerana) and is one of the main bee-producing areas in China. However, research on the genetic diversity of Chinese honeybee in the Guizhou region is very limited, despite implications for conservation of biodiversity. Methods: In this study, we analyzed the genetic diversity, differentiation, and selection signals based on 116 Chinese honeybees from 12 regions in Guizhou Province using whole-genome sequencing. Results: We identified 1,400,430 high-quality SNPs across all samples. A population structure analysis revealed two independent genetic subgroups of Chinese honeybees in Guizhou, a Yunnan-Guizhou Plateau population in western Guizhou and a hilly-mountainous population in eastern Guizhou. The average nucleotide diversity (Pi) ranged from 0.00138 to 0.00161 and average expected heterozygosity (He) ranged from 0.2592 to 0.2604. The average genetic differentiation index (F ST) for Chinese honeybees in pairwise comparisons of 12 regions ranged from 0.0094 to 0.0293. There was clear genetic differentiation between the western plateau and the eastern hilly mountainous areas of Guizhou; however, F ST values between the eastern and western populations ranged from 0.0170 to 0.0293, indicating a low degree of differentiation. A genome-wide scan revealed a number of genes under selection in the Yunnan-Guizhou Plateau environment. These genes were related to growth and development, reproduction, and cold resistance, and several candidate genes involved in environmental adaptation were identified, including CTR, MAPK, MAST, HSF, and MKKK. Discussion: The results of the present study provide important theoretical bases for the conservation, evaluation, development, and utilization of genetic resources for Chinese honeybees in the Guizhou region and for further investigations of environmental adaptation and underlying mechanisms in the species.

In understudied regions of the world, beekeeper records can provide valuable insights into changes in pollinator population trends. We conducted a questionnaire survey of 116 beekeepers in a mountainous area of Western Nepal, where the native honeybee Apis cerana cerana is kept as a managed bee. We complemented the survey with field data on insect-crop visitation, a household income survey, and an interview with a local lead beekeeper. In total, 76% of beekeepers reported declines in honeybees, while 86% and 78% reported declines in honey yield and number of beehives, respectively. Honey yield per hive fell by 50% between 2012 and 2022, whilst the number of occupied hives decreased by 44%. Beekeepers ranked climate change and declining flower abundance as the most important drivers of the decline. This raises concern for the future food and economic security of this region, where honey sales contribute to 16% of total household income, and where Apis cerana cerana plays a major role in crop pollination, contributing more than 50% of all flower visits to apple, cucumber, and pumpkin. To mitigate further declines, we promote native habitat and wildflower preservation, and using well-insulated log hives to buffer bees against the increasingly extreme temperature fluctuations.

The honeybee, Apis cerana cerana (Ac), is an important pollinator and has adapted to the local ecological environment with relevant coloration. The cuticle coloration of the brown (br) mutant is brown instead of black in wild-type individuals. Therefore, this study aimed to identify and characterize the gene responsible for the br mutation. Genome resequencing with allele segregation measurement using Euclidean distance followed by Lowess regression analysis revealed that the color locus linked to the mutation was located on chromosome 11. A 2-base deletion on exon 4 was identified in the g7628 (yellow) gene after genome assembly and sequence cloning. In addition, the cuticle color of the abdomen of worker bees changed from black to brown when a defect was induced in the yellow gene using short interfering RNA (siRNA); however, the survival rate did not decrease significantly. These results indicate that the yellow gene participated in the body pigmentation, and its defect was responsible for the br mutation. This study promotes the understanding of the molecular basis of body coloration in honeybees, enriching the molecular mechanisms underlying insect pigmentation.

Insulin receptors (InR) are an integral component of the insulin/insulin-like growth factor signaling pathway, which plays a vital role in insect development, lifespan, reproduction, and olfactory sensitivity. However, whether InR participate in the peripheral olfactory system of insects remains unclear. Recently, we found that 2-heptanone (2-HT) affects AcerInR expression, the gene for an InR protein, in Apis cerana cerana. We then examined the spatiotemporal expression profile of the gene in A. cerana cerana. The mRNA of AcerInR was primarily expressed in the antennae, wings, and legs of forager bees, which are probable chemosensory tissues. The results of fluorescence competitive binding assays, combined with site-directed mutagenesis, demonstrated that AcerOBP6 and AcerOBP14 exhibit strong binding affinities to 2-HT. Furthermore, after foragers were fed with double-stranded AcerInR, the expression levels of AcerOBP6 and AcerOBP14 decreased significantly, as did the electroantennogram responsiveness to 2-HT and some other odorants. In conclusion, our findings provide a foundation for understanding the involvement of AcerInR in the odor perception of A. cerana cerana. Moreover, they offer novel insights into the olfactory recognition mechanism in insects.

Long non-coding RNAs (lncRNAs) are crucial modulators in a variety of biological processes, such as gene expression, development, and immune defense. However, little is known about the function of lncRNAs in the development of Asian honey bee (Apis cerana) larval guts. Here, on the basis of our previously obtained deep-sequencing data from the 4-, 5-, and 6-day-old larval guts of A. cerana workers (Ac4, Ac5, and Ac6 groups), an in-depth transcriptome-wide investigation was conducted to decipher the expression pattern, regulatory manners, and potential roles of lncRNAs during the developmental process of A. cerana worker larval guts, followed by the verification of the relative expression of differentially expressed lncRNAs (DElncRNAs) and the targeting relationships within a competing endogenous RNA (ceRNA) axis. In the Ac4 vs. Ac5 and Ac5 vs. Ac6 comparison groups, 527 and 498 DElncRNAs were identified, respectively, which is suggestive of the dynamic expression of lncRNAs during the developmental process of larval guts. A cis-acting analysis showed that 330 and 393 neighboring genes of the aforementioned DElncRNAs were respectively involved in 29 and 32 functional terms, such as cellular processes and metabolic processes; these neighboring genes were also respectively engaged in 246 and 246 pathways such as the Hedgehog signaling pathway and the Wnt signaling pathway. Additionally, it was found that 79 and 76 DElncRNAs as potential antisense lncRNAs may, respectively, interact with 72 and 60 sense-strand mRNAs. An investigation of competing endogenous RNA (ceRNA) networks suggested that 75 (155) DElncRNAs in the Ac4 vs. Ac5 (Ac5 vs. Ac6) comparison group could target 7 (5) DEmiRNAs and further bind to 334 (248) DEmRNAs, which can be annotated to 33 (29) functional terms and 186 (210) pathways, including 12 (16) cellular- and humoral-immune pathways (lysosome pathway, necroptosis, MAPK signaling pathway, etc.) and 11 (10) development-associated signaling pathways (Wnt, Hippo, AMPK, etc.). The RT-qPCR detection of five randomly selected DElncRNAs confirmed the reliability of the used sequencing data. Moreover, the results of a dual-luciferase reporter assay were indicative of the binding relationship between MSTRG.11294.1 and miR-6001-y and between miR-6001-y and ncbi_107992440. These results demonstrate that DElncRNAs are likely to modulate the developmental process of larval guts via the regulation of the source genes\' transcription, interaction with mRNAs, and ceRNA networks. Our findings not only yield new insights into the developmental mechanism underlying A. cerana larval guts, but also provide a candidate ceRNA axis for further functional dissection.

Long non-coding RNAs (lncRNAs) play an essential part in controlling gene expression and a variety of biological processes such as immune defense and stress-response. However, whether and how lncRNAs regulate responses of Apis cerana larvae to Ascosphaera apis invasion has remained unclear until now. Here, the identification and structural analysis of lncRNAs in the guts of A. cerana worker larvae were conducted, and the expression profile of larval lncRNAs during the A. apis infection process was then analyzed, followed by an investigation of the regulatory roles of differentially expressed lncRNAs (DElncRNAs) in the host response. In total, 76 sense lncRNAs, 836 antisense lncRNAs, 184 intron lncRNAs, 362 bidirectional lncRNAs, and 2181 intron lncRNAs were discovered in the larval guts. Additionally, 30 known and 9 novel lncRNAs were potential precursors for 36 and 11 miRNAs, respectively. In the three comparison groups, 386, 351, and 272 DElncRNAs were respectively identified, indicating the change in the overall expression pattern of host lncRNAs following the A. apis invasion. Analysis of cis-acting effect showed that DElncRNAs in the 4-, 5-, and 6-day-old comparison groups putatively regulated 55, 30, and 20 up- and down-stream genes, respectively, which were involved in a series of crucial functional terms and pathways, such as MAPK signaling pathway, and cell process. Analysis showed that 31, 8, and 11 DElncRNAs as potential antisense lncRNAs may interact with 26, 8, and 9 sense-strand mRNAs. Moreover, investigation of the competing endogenous RNA (ceRNA) network indicated that 148, 283, and 257 DElncRNAs were putatively regulated. The expression of target genes by targeting corresponding DEmiRNAs included those associated with antioxidant enzymes and immune responses. These results suggested that DElncRNAs played a potential part in the larval guts responding to the A. apis infection through a cis-acting manner and ceRNA mechanisms. Our findings deepen our understanding of interactions between A. cerana larvae and A. apis and offer a basis for clarifying the DElncRNA-mediated mechanisms underlying the host response to fungal invasion.

Environmental pollution has gained negative attention in recent years. The pesticides and heavy metals are top list of environmental toxicants directly endangering the survival and development of Apis cerana cerana. Cyclin-dependent kinases (CDKs) are heteromeric serine/threonine kinases that participate in cell cycle regulation and have a vital role in pesticide and heavy metal stress in Apis cerana cerana. In this experiment, we filtered out CDK8 gene from Apis cerana cerana (AccCDK8) and investigated its functions of pesticide and heavy metals resistance. Sequence analysis indicated that AccCDK8 is highly homologous to multiple CDK8s and contains a highly conserved CDK active site sequence. Phylogenetic analysis showed that AmCDK8 and AccCDK8 were closely related evolutionarily in Apis mellifera. Transcriptome analysis revealed that AccCDK8 expression was differentially affected after exposure to pesticide and heavy metal stresses. This indicates that AccCDK8 has a significant role in the resistance of Apis cerana cerana to pesticide and heavy metal stresses. It has implications for studying the function of CDK in other insects in response to stress.

Heavy metals and pesticides represent prominent sources of pollution in the natural habitat of Apis cerana cerana, potentially endangering their health through the induction of oxidative stress reactions. This study aimed to address this issue by isolating AccCDK2-like and AccCINP-like proteins from Apis cerana cerana and investigating their functional roles in honey bee resistance against pesticide and heavy metal stresses. Bioinformatics analysis revealed significant homology of these proteins with those found in other species. Functional studies confirmed their participation in interaction with each other, alongside demonstrating distinct patterns of expression and localization. Specifically, AccCDK2-like exhibited higher expression levels in prepupae and muscle tissues, while AccCINP-like showed maximal expression in brown pupae and abdomen. Furthermore, the expression levels of these proteins were found to be modulated in response to pesticide and heavy metal stresses. Notably, overexpression of AccCDK2-like and AccCINP-like led to a noticeable alteration in E. coli\'s ability to withstand external stresses. Additionally, silencing of the AccCDK2-like and AccCINP-like genes resulted in a significant reduction in antioxidant enzyme activity and the expression levels of genes related to antioxidant function. Consequently, the mortality rate of Apis cerana cerana under pesticide and heavy metal stresses conspicuously increased. Hence, our findings suggest that AccCDK2-like and AccCINP-like proteins potentially play a crucial role in the response of Apis cerana cerana to pesticide and heavy metal stress, likely by modulating the antioxidant pathway.

BACKGROUND: Environmental stress can induce oxidative stress in Apis cerana cerana, leading to cellular oxidative damage, reduced vitality, and even death. Currently, owing to an incomplete understanding of the molecular mechanisms by which A. cerana cerana resists oxidative damage, there is no available method to mitigate the risk of this type of damage. Cyclin plays an important role in cell stress resistance. The aim of this study was to explore the in vivo protection of cyclin H against oxidative damage induced by abiotic stress in A. cerana cerana and clarify the mechanism of action. We isolated and identified the AccCyclin H gene in A. cerana cerana and analysed its responses to different exogenous stresses.
RESULTS: The results showed that different oxidative stressors can induce or inhibit the expression of AccCyclin H. After RNA-interference-mediated AccCyclin H silencing, the activity of antioxidant-related genes and related enzymes was inhibited, and trehalose metabolism was reduced. AccCyclin H gene silencing reduced A. cerana cerana high-temperature tolerance. Exogenous trehalose supplementation enhanced the total antioxidant capacity of A. cerana cerana, reduced the accumulation of oxidants, and improved the viability of A. cerana cerana under high-temperature stress.
CONCLUSIONS: Our findings suggest that trehalose can alleviate adverse stress and that AccCyclin H may participate in oxidative stress reactions by regulating trehalose metabolism. © 2023 Society of Chemical Industry.