Discovery of SNP Molecular Markers and Candidate Genes Associated with Sacbrood Virus Resistance in Apis cerana cerana Larvae by Whole-Genome Resequencing.

Sacbrood virus (SBV) is a significant problem that impedes brood development in both eastern and western honeybees. Whole-genome sequencing has become an important tool in researching population genetic variations. Numerous studies have been conducted using multiple techniques to suppress SBV infection in honeybees, but the genetic markers and molecular mechanisms underlying SBV resistance have not been identified. To explore single nucleotide polymorphisms (SNPs), insertions, deletions (Indels), and genes at the DNA level related to SBV resistance, we conducted whole-genome resequencing on 90 Apis cerana cerana larvae raised in vitro and challenged with SBV. After filtering, a total of 337.47 gigabytes of clean data and 31,000,613 high-quality SNP loci were detected in three populations. We used ten databases to annotate 9359 predicted genes. By combining population differentiation index (FST) and nucleotide polymorphisms (π), we examined genome variants between resistant (R) and susceptible (S) larvae, focusing on site integrity (INT < 0.5) and minor allele frequency (MAF < 0.05). A selective sweep analysis with the top 1% and top 5% was used to identify significant regions. Two SNPs on the 15th chromosome with GenBank KZ288474.1_322717 (Guanine > Cytosine) and KZ288479.1_95621 (Cytosine > Thiamine) were found to be significantly associated with SBV resistance based on their associated allele frequencies after SNP validation. Each SNP was authenticated in 926 and 1022 samples, respectively. The enrichment and functional annotation pathways from significantly predicted genes to SBV resistance revealed immune response processes, signal transduction mechanisms, endocytosis, peroxisomes, phagosomes, and regulation of autophagy, which may be significant in SBV resistance. This study presents novel and useful SNP molecular markers that can be utilized as assisted molecular markers to select honeybees resistant to SBV for breeding and that can be used as a biocontrol technique to protect honeybees from SBV.

Degradation of an appetitive olfactory memory via devaluation of sugar reward is mediated by 5-HT signaling in the honey bee.

Conditioned taste aversion (CTA) learning induces the devaluation of a preferred food through its pairing with a stimulus inducing internal illness. In invertebrates, it is still unclear how this aversive learning impairs the memories of stimuli that had been associated with the appetitive food prior to its devaluation. Here we studied this phenomenon in the honey bee and characterized its neural underpinnings. We first trained bees to associate an odorant (conditioned stimulus, CS) with appetitive fructose solution (unconditioned stimulus, US) using a Pavlovian olfactory conditioning. We then subjected the bees that learned the association to a CTA training during which the antennal taste of fructose solution was contingent or not to the ingestion of quinine solution, which induces malaise a few hours after ingestion. Only the group experiencing contingent fructose stimulation and quinine-based malaise exhibited a decrease in responses to the fructose and a concomitant decrease in odor-specific retention in tests performed 23 h after the original odor conditioning. Furthermore, injection of dopamine- and serotonin-receptor antagonists after CTA learning revealed that this long-term decrease was mediated by serotonergic signaling as its blockade rescued both the responses to fructose and the odor-specific memory 23 h after conditioning. The impairment of a prior CS memory by subsequent CTA conditioning confirms that bees retrieve a devaluated US representation when presented with the CS. Our findings further highlight the importance of serotonergic signaling in aversive learning in the bee and uncover mechanisms underlying aversive memories induced by internal illness in invertebrates.

Population genomics of honey bees reveals a selection signature indispensable for royal jelly production.

In order to interpret the molecular mechanisms that modulating the organism variations and selection signatures to drive adaptive evolutionary changes are indispensable goals in the new evolutionary ecological genetics. Here, we identified the gene locus associated to royal jelly production through whole-genome sequencing of the DNA from eight populations of honeybees. The analysis of the samples was composed of 120 individuals and each pointed extremely opposite trait values for a given phenotype. We identified functional single nucleotide polymorphisms (SNPs) candidate that might be essential in regulating the phenotypic traits of honeybee populations. Moreover, selection signatures were investigated using pooling sequencing of eight distinct honeybee populations, and the results provided the evidence of signatures of recent selection among populations under different selection objectives. Furthermore, gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that selected genes were potentially involved in several biological processes and molecular functioning, which could directly or indirectly influence the production of royal jelly. Our findings can be used to understand the genomic signatures, as well as implicate a profound glance on genomic regions that control the production trait of royal jelly in honey bees.

Royal jelly-derived proteins enhance proliferation and migration of human epidermal keratinocytes in an in vitro scratch wound model.

BACKGROUND: Skin injury is inevitable in daily life. In recent years, with the increasing morbidity of diseases such as diabetes and metabolic disorders, chronic wounds have become a considerable challenge in clinical practice. Royal jelly, reported to have multifarious biological and physiological properties, has been used as a remedy for a variety of wounds since ancient times. However, the active components and mechanisms underlying the wound-healing properties of royal jelly are still largely unknown. METHODS: Water-soluble proteins of royal jelly were fractionated and investigated for the proliferative and migratory effects on human epidermal keratinocytes (HaCaT) in an in vitro wound healing model. The proteins present in bioactive fractions were characterised and quantified using Label-free protein quantification method. The potential functions of these proteins in biological systems were further analysed using bioinformatic tools. RESULTS: A protein fraction, mainly containing major royal jelly proteins 2 (MRJP2), MRJP3 and MRJP7, stimulated proliferative and migratory activities in HaCaT cells without visible cytotoxicity. It exerted the greatest effects on the growth of HaCaT cells in the first 48 h. Furthermore, when treated with this protein fraction, the closure rates of the in vitro scratch wound were significantly increased. Functional analysis indicated that MRJP2, MRJP3 and MRJP7 were associated with carbohydrate transport and metabolism. CONCLUSIONS: We fractionated the water-soluble proteins of royal jelly and identified one fraction (Fraction 2) that induced both proliferative and migratory effects on a human epidermal keratinocyte cell line. Major royal jelly proteins (MRJP2, MRJP3 and/or MRJP7) were speculated to possess potential wound-healing bioactivity. This is the first report that royal jelly may improve wound closure via MRJP-induced cellular proliferation and migration. These proteins may be valuable lead compounds for the development of novel wound healing medications. Our findings would facilitate better understanding of the wound repair mechanisms of royal jelly.

Brain transcriptome of honey bees (Apis mellifera) exhibiting impaired olfactory learning induced by a sublethal dose of imidacloprid.

Declines in honey bee populations represent a worldwide concern. The widespread use of neonicotinoid insecticides has been one of the factors linked to these declines. Sublethal doses of a neonicotinoid insecticide, imidacloprid, has been reported to cause olfactory learning deficits in honey bees via impairment of the target organ, the brain. In the present study, olfactory learning of honey bees was compared between controls and imidacloprid-treated bees. The brains of imidacloprid-treated and control bees were used for comparative transcriptome analysis by RNA-Seq to elucidate the effects of imidacloprid on honey bee learning capacity. The results showed that the learning performance of imidacloprid-treated bees was significantly impaired in comparison with control bees after chronic oral exposure to imidacloprid (0.02 ng/µl) for 11 days. Gene expression profiles between imidacloprid treatment and the control revealed that 131 genes were differentially expressed, of which 130 were downregulated in imidacloprid-treated bees. Validation of the RNA-Seq data using qRT-PCR showed that the results of qRT-PCR and RNA-Seq exhibited a high level of agreement. Gene ontology annotation indicated that the oxidation-reduction imbalance might exist in the brain of honey bees due to oxidative stress induced by imidacloprid exposure. KEGG and ingenuity pathway analysis revealed that transient receptor potential and Arrestin 2 in the phototransduction pathway were significantly downregulated in imidacloprid-treated bees, and that five downregulated genes have causal effects on behavioral response inhibition in imidacloprid-treated bees. Our results suggest that downregulation of brain genes involved in immune, detoxification and chemosensory responses may result in decreased olfactory learning capabilities in imidacloprid-treated bees.

Comparative transcriptome analysis of Apis mellifera antennae of workers performing different tasks.

Honey bee is a social insect. Its colony is mainly coordinated by the chemical signals such as pheromones produced by queen or brood. Correspondingly, the worker bee developed numerous complicated olfactory sensilla in antennae for detection of these colony chemical signals and nectar/pollen signals in foraging. With the normal development of new emerged workers, young adults (nurse bee) worked in colony at the first 2-3 weeks and then followed by the foraging activity outside of the hive, which give rise to great change of the surrounding chemical signals. However, the olfactory adaption mechanism of worker bee in these processes of behavioral development is still unclear. In this study, we conducted a comprehensive and quantitative analysis of gene expression in Apis mellifera antenna of newly emerged workers, nurses and foragers using transcriptome analysis. Meanwhile, we constructed experimental colonies to collect age-matched samples, which were used to determine whether task is the principal determinant of differential expression. RNA sequencing and quantitative real-time polymerase chain reaction revealed that 6 and 14 genes were closely associated with nurse and forager behaviors, respectively. Furthermore, a broad dynamic range of chemosensory gene families and candidate odorant degrading enzymes were analyzed at different behavior statuses. We firstly reported genes associated with nursing/foraging behavior from antennae and the variations of expression of genes belonging to various olfactory gene families at different development stages. These results not only could contribute to elucidating the relationship between olfactory and behavior-related changes, but also provide a new perspective into the molecular mechanism underlying honey bee division of labor.

Phylogenetic analysis and survey of Apis cerana strain of Sacbrood virus (AcSBV) in Taiwan suggests a recent introduction.

The Sacbrood virus (SBV) is widely distributed in European honey bees, Apis mellifera. AcSBV, a distinct SBV strain in Asian honey bees (A. cerana) causes larva death before pupation and often depopulates colonies, leading to collapse. It is the most severe disease in A. cerana beekeeping. AcSBV infects A. cerana in most natural habitats, yet occurrences were not reported in Taiwan before 2015 and were not a concern for local beekeepers. However, in 2016, A. cerana beekeepers in central Taiwan reported SBV-like symptoms. We screened samples of larvae using RT-PCR and surveyed asymptomatic apiaries in north Taiwan. Phylogenetic analyses suggested that AcSBV isolates from central Taiwan were introduced; all isolates had high similarity in sequences to AcSBV genomes identified in mainland China, Vietnam, and Korea and distinct differences to SBV sequence identified in Taiwan. The overall prevalence in symptomatic colonies was low. No latent infections were detected in asymptomatic colonies. The AcSBV epizootic may not yet have reached its highest potential.

Spore load and immune response of honey bees naturally infected by Nosema ceranae.

Nosema ceranae causes widespread infection in adult workers of European honey bees, Apis mellifera, and has often been linked to honey bee colony losses worldwide. Previous investigations of honey bee immune response to N. ceranae infection were largely based on laboratory experiment, however, little is known about the immune response of honey bees that are naturally infected by N. ceranae. Here, we compared the infection levels of N. ceranae in three different categories of adult bees (emergent bees, nurses, and foragers) and detected the host immune response to the N. ceranae infection under natural conditions. Our studies showed that the Nosema spore load and infection prevalence varied among the different types of adult workers, and both of them increased as honey bees aged: No infection was detected in emergent bees, nurses had a medium spore load and prevalence, while foragers were with the highest Nosema infection level and prevalence. Quantification of the mRNA levels of antimicrobial peptides (abaecin, apidaecin, defensin-1, defensin-2, and hymenoptaecin) and microbial recognition proteins (PGRP-S1, PGRP-S2, PGRP-S3, PGRP-LC, GNBP1-1, and GNBP1-2) confirmed the involvement of the Toll and/or Imd immune pathways in the host response to N. ceranae infection, and revealed an activation of host immune response by N. ceranae infection under natural conditions. Additionally, the levels of immune response were positively correlated with the Nosema spore loads in the infected bees. The information gained from this study will be relevant to the predictive modeling of honey bee disease dynamics for Nosema disease prevention and management.

Differential physiological effects of neonicotinoid insecticides on honey bees: A comparison between Apis mellifera and Apis cerana.

Acute toxicities (LD50s) of imidacloprid and clothianidin to Apis mellifera and A. cerana were investigated. Changing patterns of immune-related gene expressions and the activities of four enzymes between the two bee species were compared and analyzed after exposure to sublethal doses of insecticides. Results indicated that A. cerana was more sensitive to imidacloprid and clothianidin than A. mellifera. The acute oral LD50 values of imidacloprid and clothianidin for A. mellifera were 8.6 and 2.0ng/bee, respectively, whereas the corresponding values for A. cerana were 2.7 and 0.5ng/bee. The two bee species possessed distinct abilities to mount innate immune response against neonicotinoids. After 48h of imidacloprid treatment, carboxylesterase (CCE), prophenol oxidase (PPO), and acetylcholinesterase (AChE) activities were significantly downregulated in A. mellifera but were upregulated in A. cerana. Glutathione-S-transferase (GST) activity was significantly elevated in A. mellifera at 48h after exposure to imidacloprid, but no significant change was observed in A. cerana. AChE was downregulated in both bee species at three different time points during clothianidin exposure, and GST activities were upregulated in both species exposed to clothianidin. Different patterns of immune-related gene expression and enzymatic activities implied distinct detoxification and immune responses of A. cerana and A. mellifera to imidacloprid and clothianidin.

Ranatensin-HL: A Bombesin-Related Tridecapeptide from the Skin Secretion of the Broad-Folded Frog, Hylarana latouchii.

Bombesin-related peptides are a family of peptides whose prototype was discovered in amphibian skin and which exhibit a wide range of biological activities. Since the initial isolation of bombesin from Bombina bombina skin, diverse forms of bombesin-related peptides have been found in the skins across Anura. In this study, a novel bombesin-related peptide of the ranatensin subfamily, named ranatensin-HL, was structurally-characterised from the skin secretion of the broad-folded frog, Hylarana latouchii, through combination of molecular cloning and mass spectrometric methodologies. It is composed of 13 amino acid residues, pGlu-RAGNQWAIGHFM-NH2, and resembles an N-terminally extended form of Xenopus neuromedin B. Ranatensin-HL and its C-terminal decapeptide (ranatensin-HL-10) were chemically synthesised and subjected to in vitro smooth muscle assays in which they were found to display moderate stimulatory effects on rat urinary bladder and uterus smooth muscles with EC50 values in the range of 1-10 nM. The prepro-ranatensin-HL was highly homological to a bombesin-like peptide from Rana catesbeiana at both nucleotide and amino acid levels, which might provide a clue for the taxonomic classification of ranid frogs in the future.

Identification of genes related to high royal jelly production in the honey bee (Apis mellifera) using microarray analysis.

China is the largest royal jelly producer and exporter in the world, and high royal jelly-yielding strains have been bred in the country for approximately three decades. However, information on the molecular mechanism underlying high royal jelly production is scarce. Here, a cDNA microarray was used to screen and identify differentially expressed genes (DEGs) to obtain an overview on the changes in gene expression levels between high and low royal jelly producing bees. We developed a honey bee gene chip that covered 11,689 genes, and this chip was hybridised with cDNA generated from RNA isolated from heads of nursing bees. A total of 369 DEGs were identified between high and low royal jelly producing bees. Amongst these DEGs, 201 (54.47%) genes were up-regulated, whereas 168 (45.53%) were down-regulated in high royal jelly-yielding bees. Gene ontology (GO) analyses showed that they are mainly involved in four key biological processes, and pathway analyses revealed that they belong to a total of 46 biological pathways. These results provide a genetic basis for further studies on the molecular mechanisms involved in high royal jelly production.

Silencing the Honey Bee (Apis mellifera) Naked Cuticle Gene (nkd) Improves Host Immune Function and Reduces Nosema ceranae Infections.

Nosema ceranae is a new and emerging microsporidian parasite of European honey bees, Apis mellifera, that has been implicated in colony losses worldwide. RNA interference (RNAi), a posttranscriptional gene silencing mechanism, has emerged as a potent and specific strategy for controlling infections of parasites and pathogens in honey bees. While previous studies have focused on the silencing of parasite/pathogen virulence factors, we explore here the possibility of silencing a host factor as a mechanism for reducing parasite load. Specifically, we used an RNAi strategy to reduce the expression of a honey bee gene, naked cuticle (nkd), which is a negative regulator of host immune function. Our studies found that nkd mRNA levels in adult bees were upregulated by N. ceranae infection (and thus, the parasite may use this mechanism to suppress host immune function) and that ingestion of double-stranded RNA (dsRNA) specific to nkd efficiently silenced its expression. Furthermore, we found that RNAi-mediated knockdown of nkd transcripts in Nosema-infected bees resulted in upregulation of the expression of several immune genes (Abaecin, Apidaecin, Defensin-1, and PGRP-S2), reduction of Nosema spore loads, and extension of honey bee life span. The results of our studies clearly indicate that silencing the host nkd gene can activate honey bee immune responses, suppress the reproduction of N. ceranae, and improve the overall health of honey bees. This study represents a novel host-derived therapeutic for honey bee disease treatment that merits further exploration. IMPORTANCE: Given the critical role of honey bees in the pollination of agricultural crops, it is urgent to develop strategies to prevent the colony decline induced by the infection of parasites/pathogens. Targeting parasites and pathogens directly by RNAi has been proven to be useful for controlling infections in honey bees, but little is known about the disease impacts of RNAi silencing of host factors. Here, we demonstrate that knocking down the honey bee immune repressor-encoding nkd gene can suppress the reproduction of N. ceranae and improve the overall health of honey bees, which highlights the potential role of host-derived and RNAi-based therapeutics in controlling the infections in honey bees. The information obtained from this study will have positive implications for honey bee disease management practices.

The ability to cause infection in a pathogenic fungus uncovers a new biological feature of honey bee viruses.

We demonstrated that honey bee viruses including Deformed wing virus (DWV), Black queen cell virus (BQCV) and Israeli acute paralysis virus (IAPV) could infect and replicate in the fungal pathogen Ascosphaera apis that causes honey bee chalkbrood disease, revealing a novel biological feature of honey bee viruses. The phylogenetic analysis show that viruses of fungal and honey bee origins form two clusters in the phylogenetic trees distinctly and that host range of honey bee viruses is dynamic. Further studies are warranted to investigate the impact of the viruses on the fitness of their fungal host and phenotypic effects the virus-fungus combination has on honey bee hosts.

Comparison of Brain Gene Expression Profiles Associated with Auto-Grooming Behavior between Apis cerana and Apis mellifera Infested by Varroa destructor.

The grooming behavior of honeybees serves as a crucial auto-protective mechanism against Varroa mite infestations. Compared to Apis mellifera, Apis cerana demonstrates more effective grooming behavior in removing Varroa mites from the bodies of infested bees. However, the underlying mechanisms regulating grooming behavior remain elusive. In this study, we evaluated the efficacy of the auto-grooming behavior between A. cerana and A. mellifera and employed RNA-sequencing technology to identify differentially expressed genes (DEGs) in bee brains with varying degrees of grooming behavior intensity. We observed that A. cerana exhibited a higher frequency of mite removal between day 5 and day 15 compared to A. mellifera, with day-9 bees showing the highest frequency of mite removal in A. cerana. RNA-sequencing results revealed the differential expression of the HTR2A and SLC17A8 genes in A. cerana and the CCKAR and TpnC47D genes in A. mellifera. Subsequent homology analysis identified the HTR2A gene and SLC17A8 gene of A. cerana as homologous to the HTR2A gene and SLC17A7 gene of A. mellifera. These DEGs are annotated in the neuroactive ligand-receptor interaction pathway, the glutamatergic synaptic pathway, and the calcium signaling pathway. Moreover, CCKAR, TpnC47D, HTR2A, and SLC17A7 may be closely related to the auto-grooming behavior of A. mellifera, conferring resistance against Varroa infestation. Our results further explain the relationship between honeybee grooming behavior and brain function at the molecular level and provide a reference basis for further studies of the mechanism of honeybee grooming behavior.

Large cells suppress the reproduction of Varroa destructor.

BACKGROUND: The parasitic mite, Varroa destructor has posed a threat to the health and survival of European honey bees, Apis mellifera worldwide. There is a prevailing belief that small comb cells could provide a management tool against Varroa mites. However, the hypothesis that smaller cells can impede Varroa reproduction has not been fully tested. Here, we tested this hypothesis under laboratory conditions by using two distinct Varroa in vitro rearing systems: one involved gelatin capsules of different sizes, specifically size 00 (0.95 mL) versus size 1 (0.48 mL), and the second consisted of brood comb cells drawn on 3D printed foundations with varying cell sizes, ranging from 5.0 mm to 7.0 mm at 0.5 mm intervals. RESULTS: The results showed that mother mites in size 00 cells had significantly lower fecundity and fertility compared to those in size 1 cells. Interestingly, the reproductive suppression in larger cells could be reversed by adding an extra worker larva. Similarly, gonopore size of mother mites was smaller in size 00 cells, but restored with another host larva. Furthermore, both the fecundity and fertility of mother mites decreased linearly with the size of brood comb cells. CONCLUSIONS: Our results suggest that the reproduction of V. destructor is hindered by larger cells, possibly because larger brood cells disperse or weaken host volatile chemical cues that are crucial for Varroa reproduction. The insights derived from this study are expected to hold significant implications for the implementation of Varroa management programs. © 2024 Society of Chemical Industry.

Response to comment on "Food Wanting is Mediated by Transient Activation of Dopaminergic Signaling in the Honeybee Brain".

In a technical comment, Barron et al. (1) criticized the work of Huang et al. (2) putting the accent on the quantification of dopamine levels via high-performance liquid chromatography (HPLC), yet also including data interpretation through alternative hypotheses aimed at invalidating the original ones proposed by Huang et al. We thank the authors of this technical comment, which allows us to clarify technical aspects of our work that may have been unclear, and for promoting discussion around the conclusions of our work. Below we provide answers to the points raised in their comment.

Survey Results of Honey Bee Colony Losses in Winter in China (2009-2021).

There is growing concern that massive loss of honey bees can cause serious negative effects on biodiversity and ecosystems. Surveys of colony losses have been performed worldwide to monitor the dynamic changes and health status of honey bee colonies. Here, we present the results of surveys regarding winter colony losses from 21 provinces in China from 2009 to 2021, with a total of 1,744,324 colonies managed by 13,704 beekeepers. The total colony losses were low (9.84%; 95% Confidence Interval (CI): 9.60-10.08%) but varied among years, provinces, and scales of apiaries. As little is known about the overwintering mortality of Apis cerana, in this study, we surveyed and compared the loss rates between Apis mellifera and A. cerana in China. We found colonies of A. mellifera suffered significantly lower losses than A. cerana in China. Larger apiaries resulted in higher losses in A. mellifera, whereas the opposite was observed in A. cerana. Furthermore, we used generalized linear mixed-effects models (GLMMs) to evaluate the effects of potential risk factors on winter colony losses and found that the operation size, species, migration, migration×species interaction, and queen problems were significantly related to the loss rates. New queens can increase their colony overwintering survival. Migratory beekeepers and large operations reported lower loss rates.

Differential Brain Expression Patterns of microRNAs Related to Olfactory Performance in Honey Bees (Apis mellifera).

MicroRNAs (miRNAs) play a vital role in the nerve regulation of honey bees (Apis mellifera). This study aims to investigate the differences in expression of miRNAs in a honey bee's brain for olfactory learning tasks and to explore their potential role in a honey bee's olfactory learning and memory. In this study, 12 day old honey bees with strong and weak olfactory performances were utilized to investigate the influence of miRNAs on olfactory learning behavior. The honey bee brains were dissected, and a small RNA-seq technique was used for high-throughput sequencing. The data analysis of the miRNA sequences revealed that 14 differentially expressed miRNAs (DEmiRNAs) between the two groups, strong (S) and weak (W), for olfactory performance in honey bees were identified, which included seven up-regulated and seven down-regulated. The qPCR verification results of the 14 miRNAs showed that four miRNAs (miR-184-3p, miR-276-3p, miR-87-3p, and miR-124-3p) were significantly associated with olfactory learning and memory. The target genes of these DEmiRNAs were subjected to the GO database annotation and KEGG pathway enrichment analyses. The functional annotation and pathway analysis showed that the neuroactive ligand-receptor interaction pathway, oxidative phosphorylation, biosynthesis of amino acids, pentose phosphate pathway, carbon metabolism, and terpenoid backbone biosynthesis may be a great important pathway related to olfactory learning and memory in honey bees. Our findings together further explained the relationship between olfactory performance and the brain function of honey bees at the molecular level and provides a basis for further study on miRNAs related to olfactory learning and memory in honey bees.

Understanding of Waggle Dance in the Honey Bee (Apis mellifera) from the Perspective of Long Non-Coding RNA.

The ethological study of dance behaviour has yielded some findings since Karl Von Frisch discovered and interpreted the 'dance language' in the honey bee. However, the function and role of long non-coding RNAs on dance behaviour are hardly known until now. In this study, the differential expression patterns of lncRNAs in the brains of waggling dancers and non-dancing bees were analysed by RNA sequencing. Furthermore, lncRNA-mRNA association analysis was constructed to decipher the waggle dance. The results of RNA sequencing indicated that a total of 2877 lncRNAs and 9647 mRNAs were detected from honey bee brains. Further comparison analysis displayed that two lncRNAs, MSTRG.6803.3 and XR_003305156.1, may be involved in the waggle dance. The lncRNA-mRNA association analysis showed that target genes of differentially expressed lncRNAs in the brains between waggling dancers and non-dancing bees were mainly annotated in biological processes related to metabolic process, signalling and response to stimulus and in molecular function associated with signal transducer activity, molecular transducer activity and binding. Nitrogen metabolism was likely implicated in the modulation of the waggle dance. Our findings contribute to further understanding the occurrence and development of waggle dance.

Pharmacological effects and mechanisms of bee venom and its main components: Recent progress and perspective.

Bee venom (BV), a type of defensive venom, has been confirmed to have favorable activities, such as anti-tumor, neuroprotective, anti-inflammatory, analgesic, anti-infectivity effects, etc. This study reviewed the recent progress on the pharmacological effects and mechanisms of BV and its main components against cancer, neurological disorders, inflammatory diseases, pain, microbial diseases, liver, kidney, lung and muscle injury, and other diseases in literature during the years 2018-2021. The related target proteins of BV and its main components against the diseases include Akt, mTOR, JNK, Wnt-5α, HIF-1α, NF-κB, JAK2, Nrf2, BDNF, Smad2/3, AMPK, and so on, which are referring to PI3K/Akt/mTOR, MAPK, Wnt/ß-catenin, HIF-1α, NF-κB, JAK/STAT, Nrf2/HO-1, TrkB/CREB/BDNF, TGF-ß/Smad2/3, and AMPK signaling pathways, etc. Further, with the reported targets, the potential effects and mechanisms on diseases were bioinformatically predicted via Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, disease ontology semantic and enrichment (DOSE) and protein-protein interaction (PPI) analyses. This review provides new insights into the therapeutic effects and mechanisms of BV and its main components on diseases.