Spermidine supplementation in honey bees: Autophagy and epigenetic modifications.

Polyamines (PAs), including putrescine (Put), spermidine (Spd), and spermine (Spm), are essential polycations with wide-ranging roles in cellular functions. PA levels decline with age, making exogenous PA supplementation, particularly Spd, an intriguing prospect. Previous research in honey bees demonstrated that millimolar Spd added to their diet increased lifespan and reinforced oxidative resilience. The present study is aimed to assess the anti-aging effects of spermidine supplementation at concentrations of 0.1 and 1 mM in honey bees, focusing on autophagy and associated epigenetic changes. Results showed a more pronounced effect at the lower Spd concentration, primarily in the abdomen. Spd induced site-specific histone 3 hypoacetylation at sites K18 and 27, hyperacetylation at K9, with no change at K14 in the entire body. Additionally, autophagy-related genes (ATG3, 5, 9, 13) and genes associated with epigenetic changes (HDAC1, HDAC3, SIRT1, KAT2A, KAT6B, P300, DNMT1A, DNMT1B) were upregulated in the abdomens of honey bees. In conclusion, our findings highlight profound epigenetic changes and autophagy promotion due to spermidine supplementation, contributing to increased honey bee longevity. Further research is needed to fully understand the precise mechanisms and the interplay between epigenetic alterations and autophagy in honey bees, underscoring the significance of autophagy as a geroprotective mechanism.

A field realistic model to assess the effects of pesticides residues and adulterants on honey bee gene expression.

While studies on the sublethal effects of chemical residues in beeswax on adult honey bees are increasing, the study protocols assessing the impacts on honey bee brood in realistic conditions still need to be investigated. Moreover, little is known about the residue's effect on gene expression in honey bee brood. This study reports the effects of chlorpyriphos-ethyl, acrinathrin and stearin worker pupae exposure through contaminated or adulterated beeswax on the gene expression of some key health indicators, using a novel in vivo realistic model. Larvae were reared in acrinathrin (12.5, 25, 10 and 100 ppb) and chlorpyriphos-ethyl (5, 10, 500 and 5000 ppb) contaminated or stearin adulterated beeswax (3, 4, 5, 6 and 9%) in newly formed colonies to reduce the influence of external factors. On day 11, mortality rates were assessed. Honey bee pupae were extracted from the comb after 19 days of rearing and were analysed for the gene expression profile of four genes involved in the immune response to pathogens and environmental stress factors (Imd, dorsal, domeless and defensin), and two genes involved in detoxifications mechanisms (CYP6AS14 and CYP9Q3). We found no linear relation between the increase in the pesticide concentrations and the brood mortality rates, unlike stearin where an increase in stearin percentage led to an exponential increase in brood mortality. The immune system of pupae raised in acrinathrin contaminated wax was triggered and the expression of CYP6AS14 was significantly upregulated (exposure to 12.5 and 25 ppb). Almost all expression levels of the tested immune and detoxification genes were down-regulated when pupae were exposed to chlorpyrifos-contaminated wax. The exposure to stearin triggered the immune system and detoxification system of the pupae. The identification of substance-specific response factors might ultimately serve to identify molecules that are safer for bees and the ecosystem's health.

Construction of a Full-Length Transcriptome of Western Honeybee Midgut Tissue and Improved Genome Annotation.

Honeybees are an indispensable pollinator in nature with pivotal ecological, economic, and scientific value. However, a full-length transcriptome for Apis mellifera, assembled with the advanced third-generation nanopore sequencing technology, has yet to be reported. Here, nanopore sequencing of the midgut tissues of uninoculated and Nosema ceranae-inoculated A. mellifera workers was conducted, and the full-length transcriptome was then constructed and annotated based on high-quality long reads. Next followed improvement of sequences and annotations of the current reference genome of A. mellifera. A total of 5,942,745 and 6,664,923 raw reads were produced from midguts of workers at 7 days post-inoculation (dpi) with N. ceranae and 10 dpi, while 7,100,161 and 6,506,665 raw reads were generated from the midguts of corresponding uninoculated workers. After strict quality control, 6,928,170, 6,353,066, 5,745,048, and 6,416,987 clean reads were obtained, with a length distribution ranging from 1 kb to 10 kb. Additionally, 16,824, 17,708, 15,744, and 18,246 full-length transcripts were respectively detected, including 28,019 nonredundant ones. Among these, 43,666, 30,945, 41,771, 26,442, and 24,532 full-length transcripts could be annotated to the Nr, KOG, eggNOG, GO, and KEGG databases, respectively. Additionally, 501 novel genes (20,326 novel transcripts) were identified for the first time, among which 401 (20,255), 193 (13,365), 414 (19,186), 228 (12,093), and 202 (11,703) were respectively annotated to each of the aforementioned five databases. The expression and sequences of three randomly selected novel transcripts were confirmed by RT-PCR and Sanger sequencing. The 5' UTR of 2082 genes, the 3' UTR of 2029 genes, and both the 5' and 3' UTRs of 730 genes were extended. Moreover, 17,345 SSRs, 14,789 complete ORFs, 1224 long non-coding RNAs (lncRNAs), and 650 transcription factors (TFs) from 37 families were detected. Findings from this work not only refine the annotation of the A. mellifera reference genome, but also provide a valuable resource and basis for relevant molecular and -omics studies.

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.

The nuclear and mitochondrial genome assemblies of Tetragonisca angustula (Apidae: Meliponini), a tiny yet remarkable pollinator in the Neotropics.

BACKGROUND: The field of bee genomics has considerably advanced in recent years, however, the most diverse group of honey producers on the planet, the stingless bees, are still largely neglected. In fact, only eleven of the ~ 600 described stingless bee species have been sequenced, and only three using a long-read (LR) sequencing technology. Here, we sequenced the nuclear and mitochondrial genomes of the most common, widespread and broadly reared stingless bee in Brazil and other neotropical countries-Tetragonisca angustula (popularly known in Brazil as jataí). RESULTS: A total of 48.01 Gb of DNA data were generated, including 2.31 Gb of Pacific Bioscience HiFi reads and 45.70 Gb of Illumina short reads (SRs). Our preferred assembly comprised 683 contigs encompassing 284.49 Mb, 62.84 Mb of which (22.09%) corresponded to 445,793 repetitive elements. N50, L50 and complete BUSCOs reached 1.02 Mb, 91 contigs and 97.1%, respectively. We predicted that the genome of T. angustula comprises 17,459 protein-coding genes and 4,108 non-coding RNAs. The mitogenome consisted of 17,410 bp, and all 37 genes were found to be on the positive strand, an unusual feature among bees. A phylogenomic analysis of 26 hymenopteran species revealed that six odorant receptor orthogroups of T. angustula were found to be experiencing rapid evolution, four of them undergoing significant contractions. CONCLUSIONS: Here, we provided the first nuclear and mitochondrial genome assemblies for the ecologically and economically important T. angustula, the fourth stingless bee species to be sequenced with LR technology thus far. We demonstrated that even relatively small amounts of LR data in combination with sufficient SR data can yield high-quality genome assemblies for bees.

Behavioral and genetic correlates of heterogeneity in learning performance in individual honeybees, Apis mellifera.

Learning an olfactory discrimination task leads to heterogeneous results in honeybees with some bees performing very well and others at low rates. Here we investigated this behavioral heterogeneity and asked whether it was associated with particular gene expression patterns in the bee's brain. Bees were individually conditioned using a sequential conditioning protocol involving several phases of olfactory learning and retention tests. A cumulative score was used to differentiate the tested bees into high and low performers. The rate of CS+ odor learning was found to correlate most strongly with a cumulative performance score extracted from all learning and retention tests. Microarray analysis of gene expression in the mushroom body area of the brains of these bees identified a number of differentially expressed genes between high and low performers. These genes are associated with diverse biological functions, such as neurotransmission, memory formation, cargo trafficking and development.

Does coevolution in refugia drive mimicry in bumble bees? Insights from a South Asian mimicry group.

Müllerian mimicry was proposed to be an example of a coevolved mutualism promoted by population isolation in glacial refugia. This, however, has not been well supported in butterfly models. Here, we use genomic data to test this theory while examining the population genetics behind mimetic diversification in a pair of co-mimetic bumble bees, Bombus breviceps Smith and Bombus trifasciatus Smith. In both lineages, populations were structured by geography but not as much by color pattern, suggesting sharing of color alleles across regions of restricted gene flow and formation of mimicry complexes in the absence of genetic differentiation. Demographic analyses showed mismatches between historical effective population size changes and glacial cycles, and niche modeling revealed only mild habitat retraction during glaciation. Moreover, mimetic subpopulations of the same color form in the two lineages only in some cases exhibit similar population history and genetic divergence. Therefore, the current study supports a more complex history in this comimicry than a simple refugium-coevolution model.

Identification of the cytochrome P450 gene AccCYP6A13 in Apis cerana cerana and its response to environmental stress.

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.

Biotic interactions promote local adaptation to soil in plants.

Although different ecological factors shape adaptative evolution in natural habitats, we know little about how their interactions impact local adaptation. Here we used eight generations of experimental evolution with outcrossing Brassica rapa plants as a model system, in eight treatment groups that varied in soil type, herbivory (with/without aphids), and pollination mode (hand- or bumblebee-pollination), to study how biotic interactions affect local adaptation to soil. First, we show that several plant traits evolved in response to biotic interactions in a soil-specific way. Second, using a reciprocal transplant experiment, we demonstrate that significant local adaptation to soil-type evolved in the "number of open flowers", a trait used as a fitness proxy, but only in plants that evolved with herbivory and bee pollination. Whole genome re-sequencing of experimental lines revealed that biotic interactions caused a 10-fold increase in the number of SNPs across the genome with significant allele frequency change, and that alleles with opposite allele frequency change in different soil types (antagonistic pleiotropy) were most common in plants with an evolutionary history of herbivory and bee pollination. Our results demonstrate that the interaction with mutualists and antagonists can facilitate local adaptation to soil type through antagonistic pleiotropy.

The first ITS1 profiling of honey samples from the Southeast Asian region Lombok, Bali and Banggi Island.

Southern Asian flowers offer honeybees a diversity of nectar. Based on its geographical origin, honey quality varies. Traditional methods are less authentic than DNA-based identification. The origin of honey is determined by pollen, polyphenolic, and macro-microorganisms. In this study, amplicon sequencing targets macro-microorganisms in eDNA using the ITS1 region to explore honey's geographical location and authentication. The variety of honey samples was investigated using ITS1 with Illumina sequencing. For all four honey samples, raw sequence reads showed 979,380 raw ITS1 amplicon reads and 375 ASVs up to the phylum level. The highest total number of 202 ASVs up to phylum level identified Bali honey with 211,189 reads, followed by Banggi honey with 309,207 a total number of 111 ASVs, and Lombok represents only 63 ASVs up to phylum level with several read 458,984. Based on Shannon and Chao1, honey samples from Bali (B2) and (B3) exhibited higher diversity than honey from Lombok (B1) and green honey from Sabah (B4), while the Simpson index showed that Banggi honey (B4) had higher diversity. Honey samples had significant variance in mycobiome taxonomic composition and abundance. Zygosaccharomyces and Aspergillus were the main genera found in Lombok honey, with percentages of 68.81% and 29.76% respectively. Bali honey samples (B2 and B3) were identified as having a significant amount of the genus Aureobasidium, accounting for 40.81% and 25% of the readings, respectively. The microbiome composition of Banggi honey (B4) showed a high presence of Zygosaccharomyces 45.17% and Aureobasidium 35.24%. The ITS1 analysis effectively distinguishes between honey samples of different origins and its potential as a discriminatory tool for honey origin and authentication purposes.

The black honey bee genome: insights on specific structural elements and a first step towards pangenomes.

BACKGROUND: The honey bee reference genome, HAv3.1, was produced from a commercial line sample that was thought to have a largely dominant Apis mellifera ligustica genetic background. Apis mellifera mellifera, often referred to as the black bee, has a separate evolutionary history and is the original type in western and northern Europe. Growing interest in this subspecies for conservation and non-professional apicultural practices, together with the necessity of deciphering genome backgrounds in hybrids, triggered the necessity for a specific genome assembly. Moreover, having several high-quality genomes is becoming key for taking structural variations into account in pangenome analyses. RESULTS: Pacific Bioscience technology long reads were produced from a single haploid black bee drone. Scaffolding contigs into chromosomes was done using a high-density genetic map. This allowed for re-estimation of the recombination rate, which was over-estimated in some previous studies due to mis-assemblies, which resulted in spurious inversions in the older reference genomes. The sequence continuity obtained was very high and the only limit towards continuous chromosome-wide sequences seemed to be due to tandem repeat arrays that were usually longer than 10 kb and that belonged to two main families, the 371 and 91 bp repeats, causing problems in the assembly process due to high internal sequence similarity. Our assembly was used together with the reference genome to genotype two structural variants by a pangenome graph approach with Graphtyper2. Genotypes obtained were either correct or missing, when compared to an approach based on sequencing depth analysis, and genotyping rates were 89 and 76% for the two variants. CONCLUSIONS: Our new assembly for the Apis mellifera mellifera honey bee subspecies demonstrates the utility of multiple high-quality genomes for the genotyping of structural variants, with a test case on two insertions and deletions. It will therefore be an invaluable resource for future studies, for instance by including structural variants in GWAS. Having used a single haploid drone for sequencing allowed a refined analysis of very large tandem repeat arrays, raising the question of their function in the genome. High quality genome assemblies for multiple subspecies such as presented here, are crucial for emerging projects using pangenomes.

Variability and Number of Circulating Complementary Sex Determiner (Csd) Alleles in a Breeding Population of Italian Honeybees under Controlled Mating.

In Apis mellifera, csd is the primary gene involved in sex determination: haploid hemizygous eggs develop as drones, while females develop from eggs heterozygous for the csd gene. If diploid eggs are homozygous for the csd gene, diploid drones will develop, but will be eaten by worker bees before they are born. Therefore, high csd allelic diversity is a priority for colony survival and breeding. This study aims to investigate the variability of the hypervariable region (HVR) of the csd gene in bees sampled in an apiary under a selection scheme. To this end, an existing dataset of 100 whole-genome sequences was analyzed with a validated pipeline based on de novo assembly of sequences within the HVR region. In total, 102 allelic sequences were reconstructed and translated into amino acid sequences. Among these, 47 different alleles were identified, 44 of which had previously been observed, while 3 are novel alleles. The results show a high variability in the csd region in this breeding population of honeybees.

PABLOG: a Primer Analysis tool using a Bee-Like approach on Orthologous Genes.

RNA-seq data is currently generated in numerous non-model organisms that lack a reference genome. Nevertheless, the confirmation of gene expression levels using RT-qPCR remains necessary, and the existing techniques do not seamlessly interface with the omics pipeline workflow. Developing primers for many targets by utilising orthologous genes can be a laborious, imprecise, and subjective process, particularly for plant species that are not commonly studied and do not have a known genome. We have developed a primer design tool, named PABLOG, that analyses the alignments generated from long or short RNA-seq reads and a reference orthologous gene. PABLOG scans, much like a bee searching several flowers for pollen, and presents a sorted list of potential exon-exon junction locations, ranked according to their reliability. Through computational analysis across the whole genomes of several non-model species, we demonstrate that PABLOG performs more effectively than other methods in identifying exon-exon junctions since it generates significantly fewer false-positive results. Examination of candidate regions at the gene level, in conjunction with laboratory studies, shows that the suggested primers successfully amplified particular targets in non-model plants without any presence of genomic contamination. Our tool includes a consensus sequence feature that enables the complete process of primer design, from aligning with the target gene to determining amplification parameters. The utility can be accessed via the GitHub repository located at: https://github.com/tools4plant-omics/PABLOG.

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana.

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Early life exposure to queen mandibular pheromone mediates persistent transcriptional changes in the brain of honey bee foragers.

Behavioural regulation in insect societies remains a fundamental question in sociobiology. In hymenopteran societies, the queen plays a crucial role in regulating group behaviour by affecting individual behaviour and physiology through modulation of worker gene expression. Honey bee (Apis mellifera) queens signal their presence via queen mandibular pheromone (QMP). While QMP has been shown to influence behaviour and gene expression of young workers, we know little about how these changes translate in older workers. The effects of the queen pheromone could have prolonged molecular impacts on workers that depend on an early sensitive period. We demonstrate that removal of QMP impacts long-term gene expression in the brain and antennae in foragers that were treated early in life (1 day post emergence), but not when treated later in life. Genes important for division of labour, learning, chemosensory perception and ageing were among those differentially expressed in the antennae and brain tissues, suggesting that QMP influences diverse physiological and behavioural processes in workers. Surprisingly, removal of QMP did not have an impact on foraging behaviour. Overall, our study suggests a sensitive period early in the life of workers, where the presence or absence of a queen has potentially life-long effects on transcriptional activity.

Expression and function of the vitellogenin receptor in the hypopharyngeal glands of the honey bee Apis mellifera (Hymenoptera: Apidae) workers.

The vitellogenin receptor (VgR) is essential for the uptake and transport of the yolk precursor, vitellogenin (Vg). Vg is synthesized in the fat body, released in the hemolymph, and absorbed in the ovaries, via receptor-mediated endocytosis. Besides its important role in the reproductive pathway, Vg occurs in nonreproductive worker honey bee, suggesting its participation in other pathways. The objective was to verify if the VgR occurs in the hypopharyngeal glands of Apis mellifera workers and how Vg is internalized by these cells. VgR occurrence in the hypopharyngeal glands was evaluated by qPCR analyses of VgR and immunohistochemistry in workers with different tasks. The VgR gene is expressed in the hypopharyngeal glands of workers with higher transcript levels in nurse honey bees. VgR is more expressed in 11-day-old workers from queenright colonies, compared to orphan ones. Nurse workers with developed hypopharyngeal glands present higher VgR transcripts than those with poorly developed glands. The immunohistochemistry results showed the co-localization of Vg, VgR and clathrin (protein that plays a major role in the formation of coated vesicles in endocytosis) in the hypopharyngeal glands, suggesting receptor-mediated endocytosis. The results demonstrate that VgR performs the transport of Vg to the hypopharyngeal glands, supporting the Ovary Ground Plan Hypothesis and contributing to the understanding of the role of this gland in the social context of honey bees.

Transcriptional control of a metabolic switch regulating cellular methylation reactions is part of a common response to stress in divergent bee species.

Recent global declines in bee health have elevated the need for a more complete understanding of the cellular stress mechanisms employed by diverse bee species. We recently uncovered the biomarker lethal (2) essential for life [l(2)efl] genes as part of a shared transcriptional program in response to a number of cell stressors in the western honey bee (Apis mellifera). Here, we describe another shared stress-responsive gene, glycine N-methyltransferase (Gnmt), which is known as a key metabolic switch controlling cellular methylation reactions. We observed Gnmt induction by both abiotic and biotic stressors. We also found increased levels of the GNMT reaction product sarcosine in the midgut after stress, linking metabolic changes with the observed changes in gene regulation. Prior to this study, Gnmt upregulation had not been associated with cellular stress responses in other organisms. To determine whether this novel stress-responsive gene would behave similarly in other bee species, we first characterized the cellular response to endoplasmic reticulum (ER) stress in lab-reared adults of the solitary alfalfa leafcutting bee (Megachile rotundata) and compared this with age-matched honey bees. The novel stress gene Gnmt was induced in addition to a number of canonical gene targets induced in both bee species upon unfolded protein response (UPR) activation, suggesting that stress-induced regulation of cellular methylation reactions is a common feature of bees. Therefore, this study suggests that the honey bee can serve as an important model for bee biology more broadly, although studies on diverse bee species will be required to fully understand global declines in bee populations.

Preliminary Study on the Pathogenic Mechanism of Jujube Flower Disease in Honeybees (Apis mellifera ligustica) Based on Midgut Transcriptomics.

Honeybees are prone to poisoning, also known as jujube flower disease, after collecting nectar from jujube flowers, resulting in the tumultuous demise of foragers. The prevalence of jujube flower disease has become one of the main factors affecting the development of the jujube and beekeeping industries in Northern China. However, the pathogenic mechanisms underlying jujube flower disease in honeybees are poorly understood. Herein, we first conducted morphological observations of the midgut using HE-staining and found that jujube flower disease-affected honeybees displayed midgut damage with peritrophic membrane detachment. Jujube flower disease was found to increase the activity of chitinase and carboxylesterase (CarE) and decrease the activity of superoxide dismutase (SOD), catalase (CAT), glutathione S-transferase (GST), and the content of CYP450 in the honeybee midgut. Transcriptomic data identified 119 differentially expressed genes in the midgut of diseased and healthy honeybees, including CYP6a13, CYP6a17, CYP304a1, CYP6a14, AADC, and AGXT2, which are associated with oxidoreductase activity and vitamin binding. In summary, collecting jujube flower nectar could reduce antioxidant and detoxification capacities of the honeybee midgut and, in more severe cases, damage the intestinal structure, suggesting that intestinal damage might be the main cause of honeybee death due to jujube nectar. This study provides new insights into the pathogenesis of jujube flower disease in honeybees.

Regulatory roles of long non-coding RNAs in short-term heat stress in adult worker bees.

Long non-coding RNAs (lncRNAs) are crucial modulators of post-transcriptional gene expression regulation, cell fate determination, and disease development. However, lncRNA functions during short-term heat stress in adult worker bees are poorly understood. Here, we performed deep sequencing and bioinformatic analyses of honeybee lncRNAs. RNA interference was performed by using siRNA targeting the most highly expressed lncRNA. The silencing effect on lncRNA and the relative expression levels of seven heat shock protein (HSP) genes, were subsequently examined. Overall, 7,842 lncRNAs and 115 differentially expressed lncRNAs (DELs) were identified in adult worker bees following heat stress exposure. Structural analysis revealed that the overall expression abundance, length of transcripts, exon number, and open reading frames of lncRNAs were lower than those of mRNAs. GO analysis revealed that the target genes were mainly involved in "metabolism," "protein folding," "response to stress," and "signal transduction" pathways. KEGG analysis indicated that the "protein processing in endoplasmic reticulum" and "longevity regulating pathway-multiple species" pathways were most enriched. Quantitative real-time polymerase chain reaction (qRT-PCR) detection of the selected DELs confirmed the reliability of the sequencing data. Moreover, the siRNA experiment indicated that feeding siRNA yielded a silencing efficiency of 77.51% for lncRNA MSTRG.9645.5. Upon silencing this lncRNA, the expression levels of three HSP genes were significantly downregulated (p < 0.05), whereas those of three other HSP genes were significantly upregulated (p < 0.05). Our results provide a new perspective for understanding the regulatory mechanisms of lncRNAs in adult worker bees under short-term heat stress.

Telomere length is longer following diapause in two solitary bee species.

The mechanisms that underlie senescence are not well understood in insects. Telomeres are conserved repetitive sequences at chromosome ends that protect DNA during replication. In many vertebrates, telomeres shorten during cell division and in response to stress and are often used as a cellular marker of senescence. However, little is known about telomere dynamics across the lifespan in invertebrates. We measured telomere length in larvae, prepupae, pupae, and adults of two species of solitary bees, Osmia lignaria and Megachile rotundata. Contrary to our predictions, telomere length was longer in later developmental stages in both O. lignaria and M. rotundata. Longer telomeres occurred after emergence from diapause, which is a physiological state with increased tolerance to stress. In O. lignaria, telomeres were longer in adults when they emerged following diapause. In M. rotundata, telomeres were longer in the pupal stage and subsequent adult stage, which occurs after prepupal diapause. In both species, telomere length did not change during the 8 months of diapause. Telomere length did not differ by mass similarly across species or sex. We also did not see a difference in telomere length after adult O. lignaria were exposed to a nutritional stress, nor did length change during their adult lifespan. Taken together, these results suggest that telomere dynamics in solitary bees differ from what is commonly reported in vertebrates and suggest that insect diapause may influence telomere dynamics.