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.

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.

Comprehensive analysis of the microbiome in Apis cerana honey highlights honey as a potential source for the isolation of beneficial bacterial strains.

Background: Honey is a nutritious food made by bees from nectar and sweet deposits of flowering plants and has been used for centuries as a natural remedy for wound healing and other bacterial infections due to its antibacterial properties. Honey contains a diverse community of bacteria, especially probiotic bacteria, that greatly affect the health of bees and their consumers. Therefore, understanding the microorganisms in honey can help to ensure the quality of honey and lead to the identification of potential probiotic bacteria. Methods: Herein, the bacteria community in honey produced by Apis cerana was investigated by applying the next-generation sequencing (NGS) method for the V3-V4 hypervariable regions of the bacterial 16S rRNA gene. In addition, lactic acid bacteria (LAB) in the honey sample were also isolated and screened for in vitro antimicrobial activity. Results: The results showed that the microbiota of A. cerana honey consisted of two major bacterial phyla, Firmicutes (50%; Clostridia, 48.2%) and Proteobacteria (49%; Gammaproteobacteria, 47.7%). Among the 67 identified bacterial genera, the three most predominant genera were beneficial obligate anaerobic bacteria, Lachnospiraceae (48.14%), followed by Gilliamella (26.80%), and Enterobacter (10.16%). Remarkably, among the identified LAB, Lactobacillus kunkeei was found to be the most abundant species. Interestingly, the isolated L. kunkeei strains exhibited antimicrobial activity against some pathogenic bacteria in honeybees, including Klebsiella spp., Escherichia coli, Enterococcus faecalis, Pseudomonas aeruginosa and Staphylococcus aureus. This underscores the potential candidacy of L. kunkeei for developing probiotics for medical use. Taken together, our results provided new insights into the microbiota community in the A. cerana honey in Hanoi, Vietnam, highlighting evidence that honey can be an unexplored source for isolating bacterial strains with potential probiotic applications in honeybees and humans.

Molecular identification of Hymenopteran insects collected by using Malaise traps from Hazarganji Chiltan National Park Quetta, Pakistan.

The order Hymenoptera holds great significance for humans, particularly in tropical and subtropical regions, due to its role as a pollinator of wild and cultivated flowering plants, parasites of destructive insects and honey producers. Despite this importance, limited attention has been given to the genetic diversity and molecular identification of Hymenopteran insects in most protected areas. This study provides insights into the first DNA barcode of Hymenopteran insects collected from Hazarganji Chiltan National Park (HCNP) and contributes to the global reference library of DNA barcodes. A total of 784 insect specimens were collected using Malaise traps, out of which 538 (68.62%) specimens were morphologically identified as Hymenopteran insects. The highest abundance of species of Hymenoptera (133/538, 24.72%) was observed during August and least in November (16/538, 2.97%). Genomic DNA extraction was performed individually from 90/538 (16.73%) morphologically identified specimens using the standard phenol-chloroform method, which were subjected separately to the PCR for their molecular confirmation via the amplification of cytochrome c oxidase subunit 1 (cox1) gene. The BLAST analyses of obtained sequences showed 91.64% to 100% identities with related sequences and clustered phylogenetically with their corresponding sequences that were reported from Australia, Bulgaria, Canada, Finland, Germany, India, Israel, and Pakistan. Additionally, total of 13 barcode index numbers (BINs) were assigned by Barcode of Life Data Systems (BOLD), out of which 12 were un-unique and one was unique (BOLD: AEU1239) which was assigned for Anthidium punctatum. This indicates the potential geographical variation of Hymenopteran population in HCNP. Further comprehensive studies are needed to molecularly confirm the existing insect species in HCNP and evaluate their impacts on the environment, both as beneficial (for example, pollination, honey producers and natural enemies) and detrimental (for example, venomous stings, crop damage, and pathogens transmission).

Odorant Receptors Expressing and Antennal Lobes Architecture Are Linked to Caste Dimorphism in Asian Honeybee, Apis cerana (Hymenoptera: Apidae).

Insects heavily rely on the olfactory system for food, mating, and predator evasion. However, the caste-related olfactory differences in Apis cerana, a eusocial insect, remain unclear. To explore the peripheral and primary center of the olfactory system link to the caste dimorphism in A. cerana, transcriptome and immunohistochemistry studies on the odorant receptors (ORs) and architecture of antennal lobes (ALs) were performed on different castes. Through transcriptomesis, we found more olfactory receptor genes in queens and workers than in drones, which were further validated by RT-qPCR, indicating caste dimorphism. Meanwhile, ALs structure, including volume, surface area, and the number of glomeruli, demonstrated a close association with caste dimorphism. Particularly, drones had more macroglomeruli possibly for pheromone recognition. Interestingly, we found that the number of ORs and glomeruli ratio was nearly 1:1. Also, the ORs expression distribution pattern was very similar to the distribution of glomeruli volume. Our results suggest the existence of concurrent plasticity in both the peripheral olfactory system and ALs among different castes of A. cerana, highlighting the role of the olfactory system in labor division in insects.

Unique spatially and temporary-regulated/sex-specific expression of a long ncRNA, Nb-1, suggesting its pleiotropic functions associated with honey bee lifecycle.

Honey bees are social insects, and each colony member has unique morphological and physiological traits associated with their social tasks. Previously, we identified a long non-coding RNA from honey bees, termed Nb-1, whose expression in the brain decreases associated with the age-polyethism of workers and is detected in some neurosecretory cells and octopaminergic neurons, suggesting its role in the regulation of worker labor transition. Herein, we investigated its spatially and temporary-regulated/sex-specific expression. Nb-1 was expressed as an abundant maternal RNA during oogenesis and embryogenesis in both sexes. In addition, Nb-1 was expressed preferentially in the proliferating neuroblasts of the mushroom bodies (a higher-order center of the insect brain) in the pupal brains, suggesting its role in embryogenesis and mushroom body development. On the contrary, Nb-1 was expressed in a drone-specific manner in the pupal and adult retina, suggesting its role in the drone visual development and/or sense. Subcellular localization of Nb-1 in the brain during development differed depending on the cell type. Considering that Nb-1 is conserved only in Apidae, our findings suggest that Nb-1 potentially has pleiotropic functions in the expression of multiple developmental, behavioral, and physiological traits, which are closely associated with the honey bee lifecycle.

The dual coding of a single sex pheromone receptor in Asian honeybee Apis cerana.

In Asian honeybees, virgin queens typically only mate during a single nuptial flight before founding a colony. This behavior is controlled by the queen-released mandibular pheromone (QMP). 9-oxo-(E)-2-decenoic acid (9-ODA), a key QMP component, acts as sex pheromone and attracts drones. However, how the queens prevent additional mating remains elusive. Here, we show that the secondary QMP component methyl p-hydroxybenzoate (HOB) released by mated queens inhibits male attraction to 9-ODA. Results from electrophysiology and in situ hybridization assay indicated that HOB alone significantly reduces the spontaneous spike activity of 9-ODA-sensitive neurons, and AcerOr11 is specifically expressed in sensilla placodea from the drone's antennae, which are the sensilla that narrowly respond to both 9-ODA and HOB. Deorphanization of AcerOr11 in Xenopus oocyte system showed 9-ODA induces robust inward (regular) currents, while HOB induces inverse currents in a dose-dependent manner. This suggests that HOB potentially acts as an inverse agonist against AcerOr11.

The Effects of Artificial Diets on the Expression of Molecular Marker Genes Related to Honey Bee Health.

Honey bees are commonly used to study metabolic processes, yet the molecular mechanisms underlying nutrient transformation, particularly proteins and their effects on development, health, and diseases, still evoke varying opinions among researchers. To address this gap, we investigated the digestibility and transformation of water-soluble proteins from four artificial diets in long-lived honey bee populations (Apis mellifera ligustica), alongside their impact on metabolism and DWV relative expression ratio, using transcriptomic and protein quantification methods. Diet 2, characterized by its high protein content and digestibility, was selected for further analysis from the other studied diets. Subsequently, machine learning was employed to identify six diet-related molecular markers: SOD1, Trxr1, defensin2, JHAMT, TOR1, and vg. The expression levels of these markers were found to resemble those of honey bees who were fed with Diet 2 and bee bread, renowned as the best natural food. Notably, honey bees exhibiting chalkbrood symptoms (Control-N) responded differently to the diet, underscoring the unique nutritional effects on health-deficient bees. Additionally, we proposed a molecular model to elucidate the transition of long-lived honey bees from diapause to development, induced by nutrition. These findings carry implications for nutritional research and beekeeping, underscoring the vital role of honey bees in agriculture.

Expression of honey bee (Apis mellifera) sterol homeostasis genes in food jelly producing glands of workers.

Adult workers of Western honey bees (Apis mellifera L.) acquire sterols from their pollen diet. These food sterols are transported by the hemolymph to peripheral tissues such as the mandibular and the hypopharyngeal glands in the worker bees' heads that secrete food jelly which is fed to developing larvae. As sterols are obligatory components of biological membranes and essential precursors for molting hormone synthesis in insects, they are indispensable to normal larval development. Thus, the study of sterol delivery to larvae is important for a full understanding of honey bee larval nutrition and development. Whereas hypopharyngeal glands only require sterols for their membrane integrity, mandibular glands add sterols, primarily 24-methylenecholesterol, to its secretion. For this, sterols must be transported through the glandular epithelial cells. We have analyzed for the first time in A. mellifera the expression of genes which are involved in intracellular movement of sterols. Mandibular and hypopharyngeal glands were dissected from newly emerged bees, 6-day-old nurse bees that feed larvae and 26-day-old forager bees. The expression of seven genes involved in intracellular sterol metabolism was measured with quantitative real-time PCR. Relative transcript abundance of sterol metabolism genes was significantly influenced by the age of workers and specific genes but not by gland type. Newly emerged bees had significantly more transcripts for six out of seven genes than older bees indicating that the bulk of the proteins needed for sterol metabolism are produced directly after emergence.

A rapid return to normal: temporal gene expression patterns following cold exposure in the bumble bee Bombus impatiens.

Bumble bees are common in cooler climates and many species likely experience periodic exposure to very cold temperatures, but little is known about the temporal dynamics of cold response mechanisms following chill exposure, especially how persistent effects of cold exposure may facilitate tolerance of future events. To investigate molecular processes involved in the temporal response by bumble bees to acute cold exposure, we compared mRNA transcript abundance in Bombus impatiens workers exposed to 0°C for 75 min (inducing chill coma) and control bees maintained at a constant ambient temperature (28°C). We sequenced the 3' end of mRNA transcripts (TagSeq) to quantify gene expression in thoracic tissue of bees at several time points (0, 10, 30, 120 and 720 min) following cold exposure. Significant differences from control bees were only detectable within 30 min after the treatment, with most occurring at the 10 min recovery time point. Genes associated with gluconeogenesis and glycolysis were most notably upregulated, while genes related to lipid and purine metabolism were downregulated. The observed patterns of expression indicate a rapid recovery after chill coma, suggesting an acute differential transcriptional response during recovery from chill coma and return to baseline expression levels within an hour, with no long-term gene expression markers of this cold exposure. Our work highlights the functions and pathways important for acute cold recovery, provides an estimated time frame for recovery from cold exposure in bumble bees, and suggests that cold hardening may be less important for these heterothermic insects.

Range-wide genetic analysis of an endangered bumble bee (Bombus affinis, Hymenoptera: Apidae) reveals population structure, isolation by distance, and low colony abundance.

Declines in bumble bee species range and abundances are documented across multiple continents and have prompted the need for research to aid species recovery and conservation. The rusty patched bumble bee (Bombus affinis) is the first federally listed bumble bee species in North America. We conducted a range-wide population genetics study of B. affinis from across all extant conservation units to inform conservation efforts. To understand the species' vulnerability and help establish recovery targets, we examined population structure, patterns of genetic diversity, and population differentiation. Additionally, we conducted a site-level analysis of colony abundance to inform prioritizing areas for conservation, translocation, and other recovery actions. We find substantial evidence of population structuring along an east-to-west gradient. Putative populations show evidence of isolation by distance, high inbreeding coefficients, and a range-wide male diploidy rate of ~15%. Our results suggest the Appalachians represent a genetically distinct cluster with high levels of private alleles and substantial differentiation from the rest of the extant range. Site-level analyses suggest low colony abundance estimates for B. affinis compared to similar datasets of stable, co-occurring species. These results lend genetic support to trends from observational studies, suggesting that B. affinis has undergone a recent decline and exhibit substantial spatial structure. The low colony abundances observed here suggest caution in overinterpreting the stability of populations even where B. affinis is reliably detected interannually. These results help delineate informed management units, provide context for the potential risks of translocation programs, and help set clear recovery targets for this and other threatened bumble bee species.

Comparing modern identification methods for wild bees: Metabarcoding and image-based morphological taxonomic assignment.

With the decline of bee populations worldwide, studies determining current wild bee distributions and diversity are increasingly important. Wild bee identification is often completed by experienced taxonomists or by genetic analysis. The current study was designed to compare two methods of identification including: (1) morphological identification by experienced taxonomists using images of field-collected wild bees and (2) genetic analysis of composite bee legs (multiple taxa) using metabarcoding. Bees were collected from conservation grasslands in eastern Iowa in summer 2019 and identified to the lowest taxonomic unit using both methods. Sanger sequencing of individual wild bee legs was used as a positive control for metabarcoding. Morphological identification of bees using images resulted in 36 unique taxa among 22 genera, and >80% of Bombus specimens were identified to species. Metabarcoding was limited to genus-level assignments among 18 genera but resolved some morphologically similar genera. Metabarcoding did not consistently detect all genera in the composite samples, including kleptoparasitic bees. Sanger sequencing showed similar presence or absence detection results as metabarcoding but provided species-level identifications for cryptic species (i.e., Lasioglossum). Genus-specific detections were more frequent with morphological identification than metabarcoding, but certain genera such as Ceratina and Halictus were identified equally well with metabarcoding and morphology. Genera with proportionately less tissue in a composite sample were less likely to be detected using metabarcoding. Image-based methods were limited by image quality and visible morphological features, while genetic methods were limited by databases, primers, and amplification at target loci. This study shows how an image-based identification method compares with genetic techniques, and how in combination, the methods provide valuable genus- and species-level information for wild bees while preserving tissue for other analyses. These methods could be improved and transferred to a field setting to advance our understanding of wild bee distributions and to expedite conservation research.

Cloning, functional expression, and pharmacological characterization of inwardly rectifying potassium channels (Kir) from Apis mellifera.

Potassium channels belong to the super family of ion channels and play a fundamental role in cell excitability. Kir channels are potassium channels with an inwardly rectifying property. They play a role in setting the resting membrane potential of many excitable cells including neurons. Although putative Kir channel family genes can be found in the Apis mellifera genome, their functional expression, biophysical properties, and sensitivity to small molecules with insecticidal activity remain to be investigated. We cloned six Kir channel isoforms from Apis mellifera that derive from two Kir genes, AmKir1 and AmKir2, which are present in the Apis mellifera genome. We studied the tissue distribution, the electrophysiological and pharmacological characteristics of three isoforms that expressed functional currents (AmKir1.1, AmKir2.2, and AmKir2.3). AmKir1.1, AmKir2.2, and AmKir2.3 isoforms exhibited distinct characteristics when expressed in Xenopus oocytes. AmKir1.1 exhibited the largest potassium currents and was impermeable to cesium whereas AmKir2.2 and AmKir2.3 exhibited smaller currents but allowed cesium to permeate. AmKir1 exhibited faster opening kinetics than AmKir2. Pharmacological experiments revealed that both AmKir1.1 and AmKir2.2 are blocked by the divalent ion barium, with IC50 values of 10-5 and 10-6 M, respectively. The concentrations of VU041, a small molecule with insecticidal properties required to achieve a 50% current blockade for all three channels were higher than those needed to block Kir channels in other arthropods, such as the aphid Aphis gossypii and the mosquito Aedes aegypti. From this, we conclude that Apis mellifera AmKir channels exhibit lower sensitivity to VU041.

Population-wide modelling reveals prospects of marker-assisted selection for parasitic mite resistance in honey bees.

In 2019, a joint eight-variant model was published in which eight single nucleotide polymorphisms (SNPs) in seven Apis mellifera genes were associated with Varroa destructor drone brood resistance (DBR, i.e. mite non-reproduction in drone brood). As this model was derived from only one Darwinian Black Bee Box colony, it could not directly be applied on a population-overarching scale in the northern part of Belgium (Flanders), where beekeepers prefer the carnica subspecies. To determine whether these eight SNPs remained associated with the DBR trait on a Flemish colony-broad scope, we performed population-wide modelling through sampling of various A. mellifera carnica colonies, DBR scoring of Varroa-infested drone brood and variant genotyping. Novel eight-variant modelling was performed and the classification performance of the eight SNPs was evaluated. Besides, we built a reduced three-variant model retaining only three genetic variants and found that this model classified 76% of the phenotyped drones correctly. To examine the spread of beneficial alleles and predict the DBR probability distribution in Flanders, we determined the allelic frequencies of the three variants in 292 A. mellifera carnica queens. As such, this research reveals prospects of marker-assisted selection for Varroa drone brood resistance in honeybees.

Combined pesticides in field doses weaken honey bee (Apis cerana F.) flight ability and analyses of transcriptomics and metabolomics.

Imidacloprid, chlorpyrifos, and glyphosate rank among the most extensively employed pesticides worldwide. The effects of these pesticides and their combined on the flight capability of Apis cerana, and the potential underlying mechanisms remain uncertain. To investigate these effects, we carried out flight mill, transcriptome, and metabolome experiments. Our findings reveal that individual acute oral treatments with pesticides, specifically 20 µL of 10 ng/g imidacloprid (0.2 ng per bee), 30 ng/g chlorpyrifos (0.6 ng per bee), and 60 ng/g glyphosate (1.2 ng per bee), did not impact the flight capability of the bees. However, when bees were exposed to a combination of two or three pesticides, a notable reduction in flight duration and distance was observed. In the transcriptomic and metabolomic analyses, we identified 307 transcripts and 17 metabolites that exhibited differential expression following exposure to combined pesticides, primarily associated with metabolic pathways involved in energy regulation. Our results illuminate the intricate effects and potential hazards posed by combined pesticide exposures on bee behavior. These findings offer valuable insights into the synergistic potential of pesticide combinations and their capacity to impair bee behavior. Understanding these complex interactions is essential for comprehending the broader consequences of pesticide formulations on honey bee populations.

Serial founder effects slow range expansion in an invasive social insect.

Invasive populations often experience founder effects: a loss of genetic diversity relative to the source population, due to a small number of founders. Even where these founder effects do not impact colonization success, theory predicts they might affect the rate at which invasive populations expand. This is because secondary founder effects are generated at advancing population edges, further reducing local genetic diversity and elevating genetic load. We show that in an expanding invasive population of the Asian honey bee (Apis cerana), genetic diversity is indeed lowest at range edges, including at the complementary sex determiner, csd, a locus that is homozygous-lethal. Consistent with lower local csd diversity, range edge colonies had lower brood viability than colonies in the range centre. Further, simulations of a newly-founded and expanding honey bee population corroborate the spatial patterns in mean colony fitness observed in our empirical data and show that such genetic load at range edges will slow the rate of population expansion.

Roles of DNA Methylation in Color Alternation of Eastern Honey Bees (Apis cerana) Induced by the Royal Jelly of Western Honey Bees (Apis mellifera).

Honey bees have a very interesting phenomenon where the larval diets of two different honey bee species are exchanged, resulting in altered phenotypes, namely, a honey bee nutritional crossbreed. This is a classical epigenetic process, but its underlying mechanisms remain unclear. This study aims to investigate the contribution of DNA methylation to the phenotypic alternation of a Apis mellifera-Apis cerana nutritional crossbreed. We used a full nutritional crossbreed technique to rear A. cerana queens by feeding their larvae with A. mellifera royal-jelly-based diets in an incubator. Subsequently, we compared genome-wide methylation sequencing, body color, GC ratio, and the DMRs between the nutritional crossbreed, A. cerana queens (NQs), and control, A. cerana queens (CQs). Our results showed that the NQ's body color shifted to yellow compared to the black control queens. Genome methylation sequencing revealed that NQs had a much higher ratio of mCG than that of CQs. A total of 1020 DMGs were identified, of which 20 DMGs were enriched into key pathways for melanin synthesis, including tryptophan, tyrosine, dopamine, and phenylalanine KEGG pathways. Three key differentially methylated genes [OGDH, ALDH(NAD+) and ALDH7] showed a clear, altered DNA methylation in multiple CpG islands in NQs compared to CQs. Consequently, these findings revealed that DNA methylation participates in A. cerana-A. mellifera nutritional crossbreeding as an important epigenetic modification. This study serves as a model of cross-kingdom epigenetic mechanisms in insect body color induced by environmental factors.

Identification of chemosensory genes in the stingless bee Tetragonisca fiebrigi.

Reception of chemical information from the environment is crucial for insects' survival and reproduction. The chemosensory reception mainly occurs by the antennae and mouth parts of the insect, when the stimulus contacts the chemoreceptors located within the sensilla. Chemosensory receptor genes have been well-studied in some social hymenopterans such as ants, honeybees, and wasps. However, although stingless bees are the most representative group of eusocial bees, little is known about their odorant, gustatory, and ionotropic receptor genes. Here, we analyze the transcriptome of the proboscis and antennae of the stingless bee Tetragonisca fiebrigi. We identified and annotated 9 gustatory and 15 ionotropic receptors. Regarding the odorant receptors, we identified 204, and we were able to annotate 161 of them. In addition, we compared the chemosensory receptor genes of T. fiebrigi with those annotated for other species of Hymenoptera. We found that T. fiebrigi showed the largest number of odorant receptors compared with other bees. Genetic expansions were identified in the subfamilies 9-exon, which was also expanded in ants and paper wasps; in G02A, including receptors potentially mediating social behavior; and in GUnC, which has been related to pollen and nectar scent detection. Our study provides the first report of chemosensory receptor genes in T. fiebrigi and represents a resource for future molecular and physiological research in this and other stingless bee species.

The inheritance of alternative nest architectural traditions in stingless bees.

The transmission of complex behavior and culture in humans has long been attributed to advanced forms of social learning,1,2 which play a crucial role in our technological advancement.3 While similar phenomena of behavioral traditions and cultural inheritance have been observed in animals,1,2,4,5,6 including in primates,7 whales,8 birds,9 and even insects,10 the underlying mechanisms enabling the persistence of such animal traditions, particularly in insects, are less well understood. This study introduces pioneering evidence of enduring architectural traditions in the stingless bee Scaptotrigona depilis, which are maintained without any evidence for social learning. We demonstrate that S. depilis exhibits two distinct nest architectures, comprising either helicoidal or flat, stacked horizontal combs, which are transmitted across generations through stigmergy11,12,13,14,15,16,17-an environmental feedback mechanism whereby the presence of the existing comb structures guides subsequent construction behaviors-thereby leading to a form of environmental inheritance.18,19,20 Cross-fostering experiments further show that genetic factors or prior experience does not drive the observed variation in nest architecture. Moreover, the experimental introduction of corkscrew dislocations within the combs prompted helicoidal building, confirming the use of stigmergic building rules. At a theoretical level, we establish that the long-term equilibrium of building in the helicoidal pattern fits with the expectations of a two-state Markov chain model. Overall, our findings provide compelling evidence for the persistence of behavioral traditions in an insect, based on a simple mechanism of environmental inheritance and stigmergic interactions, without requiring any sophisticated learning mechanism, thereby expanding our understanding of how traditions can be maintained in non-human species.

Effects of microRNA-305 knockdown on brain gene expression associated with division of labor in honey bee colonies (Apis mellifera).

Division of labor in honey bee colonies is based on the behavioral maturation of adult workers that involves a transition from working in the hive to foraging. This behavioral maturation is associated with distinct task-related transcriptomic profiles in the brain and abdominal fat body that are related to multiple regulatory factors including juvenile hormone (JH) and queen mandibular pheromone (QMP). A prominent physiological feature associated with behavioral maturation is a loss of abdominal lipid mass as bees transition to foraging. We used transcriptomic and physiological analyses to study whether microRNAs (miRNAs) are involved in the regulation of division of labor. We first identified two miRNAs that showed patterns of expression associated with behavioral maturation, ame-miR-305-5p and ame-miR-375-3p. We then downregulated the expression of these two miRNAs with sequence-specific antagomirs. Neither ame-miR-305-5p nor ame-miR-375-3p knockdown in the abdomen affected abdominal lipid mass on their own. Similarly, knockdown of ame-miR-305-5p in combination with JH or QMP also did not affect lipid mass. By contrast, ame-miR-305-5p knockdown in the abdomen caused substantial changes in gene expression in the brain. Brain gene expression changes included genes encoding transcription factors previously implicated in behavioral maturation. The results of these functional genomic experiments extend previous correlative associations of microRNAs with honey bee division of labor and point to specific roles for ame-miR-305-5p.