Monitoring Immune Modulation in Season Population: Identifying Effects and Markers Related to Apis mellifera ligustica Honey Bee Health

Honey bees play a significant role in ecology, producing biologically active substances used to promote human health. However, unlike humans, the molecular markers indicating honey bee health remain unknown. Unfortunately, numerous reports of honey bee collapse have been documented. To identify health markers, we analyzed ten defense system genes in Apis mellifera ligustica honey bees from winter (Owb) and spring (Fb for foragers and Nb for newly emerged) populations sampled in February and late April 2023, respectively. We focused on colonies free from SBV and DWV viruses. Molecular profiling revealed five molecular markers of honey bee health. Of these, two seasonal molecular markers-domeless and spz genes-were significantly downregulated in Owb compared to Nb and Fb honey bees. One task-related marker gene, apid-1, was identified as being downregulated in Owb and Nb compared to Fb honey bees. Two recommended general health markers, SOD and defensin-2, were upregulated in honey bees. These markers require further testing across various honey bee subspecies in different climatic regions. They can diagnose bee health without colony intervention, especially during low-temperature months like winter. Beekeepers can use this information to make timely adjustments to nutrients or heating to prevent seasonal losses.

Transcriptomic analysis to elucidate the response of Apis mellifera ligustica brain tissue to fluvalinate exposure

As a commonly used acaricide in apiculture, fluvalinate is used to kill Varroa mites, while it also damages the nervous system of honeybees. To date, the transcriptomic characteristics associated with fluvalinate-induced neuronal injury in the bee brain have not been reported. Here, we performed transcriptome sequencing on Apis mellifera ligustica (A. mellifera ligustica) brain tissues collected before and after fluvalinate treatment. A total of 546 differentially expressed genes (DEGs) were detected, and these DEGs mainly showed 4 different expression patterns. Further analysis revealed that DEGs with different expression patterns were mainly involved in lipid metabolism, amino acid metabolism, visual transduction, and neural response-related GO terms and KEGG pathways. Moreover, protein-protein interaction network analysis revealed five protein-coding DEGs as key genes, which may play important roles in the resistance to fluvalinate-induced honeybee brain nerve tissue damage. In summary, this study is the first to perform a detailed characterization and functional analysis of genes related to fluvalinate stimulation in honeybee brains.

First Characterization and Regulatory Function of piRNAs in the Apis mellifera Larval Response to Ascosphaera apis Invasion

piRNAs are a class of small non-coding RNAs that play essential roles in modulating gene expression and abundant biological processes. To decode the piRNA-regulated larval response of western honeybees (Apis mellifera) to Ascosphaera apis infection, the expression pattern of piRNAs in Apis mellifera ligustica larval guts after A. apis inoculation was analyzed based on previously obtained high-quality small RNA-seq datasets, followed by structural characterization, target prediction, regulatory network investigation, and functional dissection. Here, 504, 657, and 587 piRNAs were respectively identified in the 4-, 5-, and 6-day-old larval guts after inoculation with A. apis, with 411 ones shared. These piRNAs shared a similar length distribution and first base bias with mammal piRNAs. Additionally, 96, 103, and 143 DEpiRNAs were detected in the 4-, 5-, and 6-day-old comparison groups. Targets of the DEpiRNAs were engaged in diverse pathways such as the phosphatidylinositol signaling system, inositol phosphate metabolism, and Wnt signaling pathway. These targets were involved in three energy metabolism-related pathways, eight development-associated signaling pathways, and seven immune-relevant pathways such as the Jak-STAT signaling pathway. The expression trends of five randomly selected DEpiRNAs were verified using a combination of RT-PCR and RT-qPCR. The effective overexpression and knockdown of piR-ame-945760 in A. apis-infected larval guts were achieved by feeding a specific mimic and inhibitor. Furthermore, piR-ame-945760 negatively regulated the expression of two target immune mRNAs, SOCS5 and ARF1, in the larval gut during the A. apis infection. These findings indicated that the overall expression level of piRNAs was increased and the expression pattern of piRNAs in larval guts was altered due to the A. apis infection, DEpiRNAs were putative regulators in the A. apis-response of A. m. ligustica worker larvae. Our data provide not only a platform for the functional investigation of piRNAs in honeybees, especially in bee larvae, but also a foundation for illuminating the piRNA-involved mechanisms underlying the host response to the A. apis infection.

Oleic Acid Promotes the Biosynthesis of 10-Hydroxy-2-decenoic Acid via Species-Selective Remodeling of TAGs in Apis mellifera ligustica

This study aimed to assess the impact of oleic acid (OA) supplementation on the biosynthesis of 10-hydroxy-2-decenoic acid (10-HDA) in Apis mellifera ligustica. In experiment 1, varying concentrations of OA (2%, 4%, 6% and 8%) were added to an artificial diet for newly emerged bees reared in cages. Analysis of 10-HDA content and gene expression in the mandibular gland (MG) revealed that the 8% OA treatment had the greatest impact on promoting the synthesis of 10-HDA. Subsequent investigations utilized RNA-seq and lipidomics to characterize the molecular signature in the MG after feeding the 8% OA diet. Phosphatidylcholine (PC) and triacylglycerol (TAG) were found to be the predominant lipids in the MG of worker bees. A total of 154 TAGs were identified, with TAG (18:1-18:1-18:1) exhibiting the highest abundance, which increased by 1.5 times. The major TAG species contained palmitic acid (16:0) and oleic acid (18:1) in their structure, which was associated with fatty acid composition of diet. The increase in abundance of main TAGs may be attributed to the upregulation of glycerol-3-phosphate acyltransferase (Gpat) and glycerol kinase (GK) gene expression at the transcriptional level. The upregulation of differentially expressed genes (DEGs) related to carbohydrate metabolism may contribute to meeting the heightened metabolic demands of the MGs in worker bees. Royal jelly (RJ) samples from bee colonies fed with the 8% OA diet exhibited higher 10-HDA level than RJ collected from bee colonies fed with the artificial diet. These results indicate that 8% OA addition in the diet enhanced biosynthesis of 10-HDA in the mandibular gland, which was accompanied by significant and highly species-selective remodeling of TAGs.

Effect of fluvalinate on the expression profile of circular RNA in brain tissue of Apis mellifera ligustica workers

Fluvalinate is widely used in apiculture as an acaricide for removing Varroa mites, but there have been growing concerns about the negative effects of fluvalinate on honeybees in recent years. Previous research revealed changes in the miRNA and mRNA expression profiles of Apis mellifera ligustica brain tissues during fluvalinate exposure, as well as key genes and pathways. The role of circRNAs in this process, however, is unknown. The goal of this study was to discover the fluvalinate-induced changes in circular RNA (circRNA) expression profiles of brain tissue of A. mellifera ligustica workers. A total of 10,780 circRNAs were detected in A. mellifera ligustica brain tissue, of which eight were differentially expressed between at least two of the four time periods before and after fluvalinate administration, and six circRNAs were experimentally verified to be structurally correct, and their expression patterns were consistent with transcriptome sequencing results. Furthermore, ceRNA analysis revealed that five differentially expressed circRNAs (DECs) (novel_circ_012139, novel_circ_011690, novel_circ_002628, novel_circ_004765, and novel_circ_010008) were primarily involved in apoptosis-related functions by competitive binding with miRNAs. This study discovered changes in the circRNA expression profile of A. mellifera ligustica brain tissue caused by fluvalinate exposure, and it provides a useful reference for the biological function study of circRNAs in A. mellifera ligustica.

Fluvalinate-Induced Changes in MicroRNA Expression Profile of Apis mellifera ligustica Brain Tissue

Fluvalinate is a widely used and relatively safe acaricide for honeybees, but it still has a negative impact on honeybee colonies. Such negative effects may be related to fluvalinate-induced brain nerve tissue damage, but the detailed molecular regulatory mechanism of this phenomenon is still poorly understood. In this study, we analyzed the miRNA expression profile changes in the brain tissue of Apis mellifera ligustica by miRNA sequencing after fluvalinate treatment. A total of 1,350 miRNAs were expressed in Apis mellifera ligustica brain tissue, of which only 180 were previously known miRNAs in honeybees. Among all known and novel miRNAs, 15 were differentially expressed between at least two of the four time periods before and after fluvalinate administration. Further analysis revealed five significantly enriched KEGG pathways of the differentially expressed miRNA (DEM) potential target genes, namely, "Hippo signaling pathway-fly," "Phototransduction-fly," "Apoptosis-fly," "Wnt signaling pathway," and "Dorso-ventral axis formation," which indicates that differentially expressed miRNA function may be related to cell apoptosis and memory impairment in the fluvalinate-treated Apis mellifera ligustica brain. Ame-miR-3477-5p, ame-miR-375-3p, and miR-281-x were identified as key miRNAs. Overall, our research provides new insights into the roles of miRNAs in brain tissue during the process of fluvalinate-induced Apis mellifera ligustica poisoning.

The effects of beta-cypermethrin, chlorbenzuron, chlorothalonil, and pendimethalin on Apis mellifera ligustica and Apis cerana cerana larvae reared in vitro

BACKGROUND

Declines in bee populations and diversity have drawn international attention. The long-term use of chemical pesticides has affected bee behavior and physiology. This study aimed to investigate the effects of chronic exposure to four commonly used chemical pesticides (beta-cypermethrin, chlorbenzuron, chlorothalonil and pendimethalin) on the growth of Apis mellifera ligustica and Apis cerana cerana larvae reared in vitro.

RESULTS

Pesticide type and concentration were the main factors affecting honeybee fitness. Beta-cypermethrin and chlorbenzuron had chronic toxic effects on bee larvae. They reduced the fitness of A. m. ligustica and A. c. cerana even at low doses of 323.5 ng g^-1 for beta-cypermethrin and 62.6 ng g^-1 for chlorbenzuron in bee bread. The effects were positively associated with the dietary amounts of pesticides. By contrast, chlorothalonil and pendimethalin exposure did not affect bee larvae despite changes in enzyme activities. Caution is still needed with chlorothalonil, which led to a decrease in harvest adult bee numbers at a high dose (6937.2 ng g^-1 ). Furthermore, a difference in pesticide resistance was observed, suggesting that A. m. ligustica may tolerate toxic effects better than A. c. cerana.

CONCLUSION

This study sheds new light on chronic toxicity in bee larvae exposed to residues in bee bread. The results could guide the scientific and rational use of chemical pesticides to reduce the potential risks to A. m. ligustica and A. c. cerana. © 2021 Society of Chemical Industry.

Characterization of the mitochondrial genome of a high royal jelly-producing honeybee strain (Apis mellifera ligustica)

The complete mitochondrial genome (mitogenome) of a high royal jelly-producing honeybee strain selected from Italian bees (Apis mellifera ligustica) in China was determined. The entire mitogenome was 16,467 bp in length with an A + T content of 84.99%. It contained the typical 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a control region. The order and orientation of these genes were identical to those of other A. mellifera subspecies sequenced to date. Phylogenetic analysis supported that the high royal jelly-producing honeybee strain was sister to formerly sequenced Italian bees.

Non-targeted lipidomics and transcriptomics analysis reveal the molecular underpinnings of mandibular gland development in Apis mellifera ligustica

The mandibular gland is an important exocrine gland of worker bees, which mainly secretes fatty acids and pheromones. Lipids have important roles in energy storage, membrane structure stabilization, and signaling. However, molecular underpinnings of mandibular gland development and lipid remodeling at the different physiological stages of worker bees is still lacking. In this study, we used scanning and transmission electron microscopy to reveal the morphological changes in secretory cells, and liquid chromatography-mass spectrometry and RNA-seq to investigate the lipidome and gene transcripts during development. The morphology of secretory cells was flat in newly emerged workers, becoming vacuolated and turgid when they were activated in nurse bees and foragers. Transport vesicles became denser from newly emerged bees to 21-day worker bees. Concentrations of 10-HDA reached a maximum within 15d workers and changes in genes expression were consistent with 10-HDA content. Non-targeted lipidomics analysis of newly emerged, 6d, and 15d worker bees revealed that PC and TAG were the main lipids in mandibular gland, and lipids dramatically altered across developmental stages. TAG 54:4 was increased most strongly at 6d and 15d worker bees, meanwhile, the abundances of TAG 54:1 and TAG 54:2 were decreased sharply. Further, transcriptomics analysis showed that differentially expressed genes were significantly enriched in key nutrient metabolic pathways, particularly lipid metabolism, in 6d and 15d bees. This multi-omic perspective provides a unique resource and deeper insight into bee mandibular gland development and baseline data for further study of the mandibular gland in worker bees.

Heavy metal bioaccumulation in honey bee matrix, an indicator to assess the contamination level in terrestrial environments

The most significant risk factor for organisms living in an environment contaminated by heavy metals is the metal bioavailability. Therefore, an efficient ecotoxicological approach to metal contamination is the measure of bioaccumulation level in target organisms. In this work, we characterized the heavy metal bioaccumulation in honey bees, Apis mellifera ligustica, collected at 35 sites from Umbria (Central Italy). The comparison of our data with selected Italian investigations revealed metal bioaccumulation in honey bee matrix of the same order of magnitude, with Cd showing a higher variability. To generalize the results, we developed a Honeybee Contamination Index (HCI) based on metal bioaccumulation in honey bees. An application of the HCI to the present dataset revealed cases of low (sixteen sites), intermediate (eighteen sites), and high (one site) metal contaminations. The comparison of HCI values from the Umbrian dataset with values calculated for other Italian and European metadata showed that most of the Umbrian sites fell in the portion of low and intermediate contamination conditions. HCI represented a reliable tool that provided a piece of concise information on metal contamination in terrestrial environments. Parallel to this effort, we have determined, the metal concentrations in the airborne particulate matter (PM₁₀) at three regional background-monitoring stations in Umbria. These stations are representative of the average air quality of the areas of the investigated apiaries. A comparative analysis of metal enrichment factors in PM₁₀, and honey bees suggested that the contamination in the bees was related to the PM₁₀ values only to a minor extent. On the other side, a clear enrichment of metals such as Cd, Mn, Zn, and Cu in the honey bees appeared to depend on very local conditions and was probably related to the use of pesticides and fertilizers, and the resuspension of the locally contaminated soils and agriculture residues.

A comprehensive transcriptome data of normal and Nosema ceranae-stressed midguts of Apis mellifera ligustica workers

Honeybees are pivotal pollinators of crops and wild flora, and of great importance in supporting critical ecosystem balance. Nosema ceranae, a unicellular fungal parasite that infects midgut epithelial cells of honeybees, can dramatically reduce honeybee population and productivity. Here, midguts of Apis mellifera ligustica workers at 7 d and 10 d post inoculation (dpi) with sucrose solution (Ac7CK and Ac10CK) and midguts at 7 dpi and 10 dpi with sucrose solution containing N. ceranae spores (Ac7T and Ac10T) were sequenced using strand-specific cDNA library construction and next-generation sequencing. A total of 1956129858 raw reads were gained in this article, and after quality control, 1946489304 high-quality clean reads with a mean Q30 of 93.82% were obtained. The rRNA-removed clean reads were then aligned to the reference genome of Apis mellifera with TopHat2. For more insight please see “Genome-wide identification of long non-coding RNAs and their regulatory networks involved in Apis mellifera ligustica response to Nosema ceranae infection” [1]. Raw data were deposited in NCBI Sequence Read Archive (SRA) database under the BioProject number PRJNA406998. These data can be used for comparative analysis to identify differentially expressed coding RNAs and non-coding RNAs involved in A. m. ligustica responses to N. ceranae stress, and for investigation of molecular mechanisms regulating host N. ceranae -response.

The complete mitochondrial genome of a Buckfast bee, Apis mellifera (Insecta: Hymenoptera: Apidae) in Northern Ireland

We analyzed the complete mitochondrial genome of the 'Buckfast bee', Apis mellifera, collected from North Ireland, UK. It consisted of a circular molecule of 16,353 bp. The genome contained 13 protein-coding, 22 tRNA, and 2 rRNA genes, along with one A + T-rich control region. The average AT content was 84.9%. The genes ATP8 and ATP6 shared 19 nucleotides. A phylogenetic analysis, suggested that the matriline 'Buckfast bee' has remained most closely related to the A. mellifera ligustica race from which it originated in 1917, despite being cross-bred with many other A. mellifera races over the past 100 years.

The complete mitochondrial genome of the yellow coloured honeybee Apis mellifera (Insecta: Hymenoptera: Apidae) of New Zealand

The complete mitochondrial genome of the yellow coloured honeybee Apis mellifera from North Island, New Zealand was analyzed using next-generation sequencing. The mitochondrial genome was a 16,349bp circular molecule and was predicted to contain 13 protein-coding genes (PCGs), 22 tRNA genes and two rRNA genes. The initiation codon ATA was found in two genes, ATG in four genes, ATT in six genes, and ATC in one gene, while the termination codon TAA was observed in all the PCGs. Phylogenetic analysis using the sequence of 23 closely related taxa suggested a sister relationship with the Italian strain A. mellifera ligustica.

Exploring poisonous mechanism of honeybee, Apis mellifera ligustica Spinola, caused by pyrethroids

As the important intracellular secondary messengers, calcium channel is the target of many neurotoxic pesticides as calcium homeostasis in the neuroplasm play important role in neuronal functions and behavior in insects. This study investigated the effect of deltamethrin (DM) on calcium channel in the brain nerve cells of adult workers of Apis mellifera ligustica Spinola that were cultured in vitro. The results showed that the intracellular calcium concentration was significantly elevated even with a very low concentration of the DM (3.125×10^-2mg/L). Further testing revealed that T-type voltage-gated calcium channels (VGCCs), except for sodium channels, was one of the target of DM on toxicity of Apis mellifera, while DM has no significant effect on the L-type VGCCs, N-methyl-d-aspartate receptor-gated calcium channels and calcium store. These results suggesting that the DM may act on T-type VGCCs in brain cells of honeybees and result in behavioral abnormalities including swarming, feeding, learning, and acquisition.

In-Depth N-Glycosylation Reveals Species-Specific Modifications and Functions of the Royal Jelly Protein from Western (Apis mellifera) and Eastern Honeybees (Apis cerana)

Royal jelly (RJ), secreted by honeybee workers, plays diverse roles as nutrients and defense agents for honeybee biology and human health. Despite being reported to be glycoproteins, the glycosylation characterization and functionality of RJ proteins in different honeybee species are largely unknown. An in-depth N-glycoproteome analysis and functional assay of RJ produced by Apis mellifera lingustica (Aml) and Apis cerana cerana (Acc) were conducted. RJ produced by Aml yielded 80 nonredundant N-glycoproteins carrying 190 glycosites, of which 23 novel proteins harboring 35 glycosites were identified. For Acc, all 43 proteins glycosylated at 138 glycosites were reported for the first time. Proteins with distinct N-glycoproteomic characteristics in terms of glycoprotein species, number of N-glycosylated sites, glycosylation motif, abundance level of glycoproteins, and N-glycosites were observed in this two RJ samples. The fact that the low inhibitory efficiency of N-glycosylated major royal jelly protein 2 (MRJP2) against Paenibacillus larvae (P. larvae) and the absence of antibacterial related glycosylated apidaecin, hymenoptaecin, and peritrophic matrix in the Aml RJ compared to Acc reveal the mechanism for why the Aml larvae are susceptible to P. larvae, the causative agent of a fatal brood disease (American foulbrood, AFB). The observed antihypertension activity of N-glycosylated MRJP1 in two RJ samples and a stronger activity found in Acc than in Aml reveal that specific RJ protein and modification are potentially useful for the treatment of hypertensive disease for humans. Our data gain novel understanding that the western and eastern bees have evolved species-specific strategies of glycosylation to fine-tune protein activity for optimizing molecular function as nutrients and immune agents for the good of honeybee and influence on the health promoting activity for human as well. This serves as a valuable resource for the targeted probing of the biological functions of RJ proteins for honeybee and medical communities.

Determination of exposure levels of honey bees foraging on flowers of mature citrus trees previously treated with imidacloprid

BACKGROUND

Field and tunnel cage studies were undertaken to determine the extent to which honey bees foraging on citrus blossoms were exposed to imidacloprid and its metabolites when citrus trees were treated with soil applications of the insecticide. Residues were measured by LC/MS/MS in nectar and pollen samples from trees treated up to 232 days prior to bloom.

RESULTS

Imidacloprid, imidacloprid olefin and 5-hydroxy imidacloprid were detected in nectar and pollen sampled from the flowers of citrus trees treated with imidacloprid up to 232 days prior to bloom. In tunnel studies, where foraging was restricted exclusively to citrus, imidacloprid residues in nectar extracted from flowers and from bee crops were similar (<10 ng mL(-1) ) and below the established no-observed-effects limit; however, the residue levels were about threefold higher in nectar sampled from comb. Concentrations of imidacloprid in nectar were higher in trees treated with higher application rates.

CONCLUSIONS

Imidacloprid and its metabolites were detected in the nectar and pollen of citrus trees treated up to 232 days prior to the onset of bloom. However, based on published bioassay data, the imidacloprid concentrations in the floral nectar did not surpass levels that would compromise foraging activity under normal use conditions for imidacloprid. Further research is needed to assess the impact of elevated levels of imidacloprid within stored nectar in the comb.