Transcriptomic analysis of instinctive and learned reward-related behaviors in honey 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.

Comparative studies of genomic and epigenetic factors influencing transcriptional variation in two insect species

Different genes show different levels of expression variability. For example, highly expressed genes tend to exhibit less expression variability. Genes whose promoters have TATA box and initiator motifs tend to have increased expression variability. On the other hand, DNA methylation of transcriptional units, or gene body DNA methylation, is associated with reduced gene expression variability in many species. Interestingly, some insect lineages, most notably Diptera including the canonical model insect Drosophila melanogaster, have lost DNA methylation. Therefore, it is of interest to determine whether genomic features similarly influence gene expression variability in lineages with and without DNA methylation. We analyzed recently generated large-scale data sets in D. melanogaster and honey bee (Apis mellifera) to investigate these questions. Our analysis shows that increased gene expression levels are consistently associated with reduced expression variability in both species, while the presence of TATA box is consistently associated with increased gene expression variability. In contrast, initiator motifs and gene lengths have weak effects limited to some data sets. Importantly, we show that a sequence characteristics indicative of gene body DNA methylation is strongly and negatively associate with gene expression variability in honey bees, while it shows no such association in D. melanogaster. These results suggest the evolutionary loss of DNA methylation in some insect lineages has reshaped the molecular mechanisms concerning the regulation of gene expression variability.

Neural and Molecular Mechanisms of Biological Embedding of Social Interactions

Animals operate in complex environments, and salient social information is encoded in the nervous system and then processed to initiate adaptive behavior. This encoding involves biological embedding, the process by which social experience affects the brain to influence future behavior. Biological embedding is an important conceptual framework for understanding social decision-making in the brain, as it encompasses multiple levels of organization that regulate how information is encoded and used to modify behavior. The framework we emphasize here is that social stimuli provoke short-term changes in neural activity that lead to changes in gene expression on longer timescales. This process, simplified-neurons are for today and genes are for tomorrow-enables the assessment of the valence of a social interaction, an appropriate and rapid response, and subsequent modification of neural circuitry to change future behavioral inclinations in anticipation of environmental changes. We review recent research on the neural and molecular basis of biological embedding in the context of social interactions, with a special focus on the honeybee.

Sex-specific molecular specialization and activity rhythm-dependent gene expression in honey bee antennae

We performed an RNA-seq-based comparison of gene expression levels in the antennae of honey bee drones and time-trained foragers (workers) collected at different times of the day and different activity states. Interestingly, olfaction-related genes [i.e. odorant receptor (Or) genes, odorant binding protein (Obp) genes, carboxyl esterase (CEst) genes, etc.] showed stable gene expression differences between drone and worker antennae. Drone antennae showed higher expression of 24 Or genes, of which 21 belong to the clade X which comprises the receptor for the major queen pheromone compound 9-ODA. This high number of drone-biased Or genes suggests that more than previously thought play a role in sex-pheromone communication. In addition, we found higher expression levels for many non-olfaction-related genes including nitric oxide synthase (NOS), and the potassium channel Shaw In contrast, workers showed higher expression of 67 Or genes, which belong to different Or clades that are involved in pheromone communication as well as the perception of cuticular hydrocarbons and floral scents. Further, drone antennae showed higher expression of genes involved in energy metabolism, whereas worker antennae showed higher expression of genes involved in neuronal communication, consistent with earlier reports on peripheral olfactory plasticity. Finally, drones that perform mating flight in the afternoon (innate) and foragers that are trained to forage in the afternoon (adapted) showed similar daily changes in the expression of two major clock genes, period and cryptochrome2 Most of the other genes showing changes with time or onset of daily flight activity were specific to drones and foragers.

Neural activity mapping of bumble bee (Bombus ignitus) brains during foraging flight using immediate early genes

Honey bees and bumble bees belong to the same family (Apidae) and their workers exhibit a division of labor, but the style of division of labor differs between species. The molecular and neural bases of the species-specific social behaviors of Apidae workers have not been analyzed. Here, we focused on two immediate early genes, hormone receptor 38 (HR38) and early growth response gene-1 (Egr1), and late-upregulated ecdysone receptor (EcR), all of which are upregulated by foraging flight and expressed preferentially in the small-type Kenyon cells of the mushroom bodies (MBs) in the honey bee brain. Gene expression analyses in Bombus ignitus revealed that HR38 and Egr1, but not EcR, exhibited an immediate early response during awakening from CO₂ anesthesia. Both premature mRNA for HR38 and mature mRNA for Egr1 were induced during foraging flight, and mRNAs for HR38 and Egr1 were sparsely detected inside the whole MB calyces. In contrast, EcR expression was higher in forager brains than in nurse bees and was expressed preferentially in the small-type Kenyon cells inside the MBs. Our findings suggest that Kenyon cells are active during foraging flight and that the function of late-upregulated EcR in the brain is conserved among these Apidae species.

Transcriptomic analysis to elucidate the response of honeybees (Hymenoptera: Apidae) to amitraz treatment

Amitraz is an acaricide that is widely used in apiculture. Several studies have reported that in honeybees (Apis mellifera Linnaeus; Hymenoptera: Apidae), amitraz affects learning, memory, behavior, immunity, and various other physiological processes. Despite this, few studies have explored the molecular mechanisms underlying the action of amitraz on honeybees. Here, we investigated the transcriptome of honeybees after exposure to 9.4 mg/L amitraz for 10 d, a subchronic dose. Overall, 279 differentially expressed genes (DEGs) were identified (237 upregulated, 42 downregulated). Several, including Pla2, LOC725381, LOC413324, LOC724386, LOC100577456, LOC551785, and P4504c3, were validated by quantitative PCR. According to gene ontology, DEGs were mainly involved in metabolism, biosynthesis, and translation. Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed that amitraz treatment affected the relaxin signaling pathway, platelet activation, and protein digestion and absorption.

Enantioselective Olfactory Effects of the Neonicotinoid Dinotefuran on Honey Bees (Apis mellifera L.)

Sublethal exposure to neonicotinoids affects honey bee olfaction, but few studies have investigated the sublethal effects of the enantioselective neonicotinoid dinotefuran on honey bee olfaction. This study assessed the sublethal olfactory toxicity of dinotefuran enantiomers to honey bees. Compared to R-dinotefuran, S-dinotefuran had higher acute oral toxicity, sucrose sensitivity effects, octopamine concentrations, lower learning ability, and memory effects on honey bees. High-throughput circular RNA sequencing of the honey bee brain revealed that R-dinotefuran caused more gene regulatory changes than S-dinotefuran. Gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses demonstrated that the SERCA, Kca, and Maxik genes may be related to the enantioselective effects of dinotefuran isomers on honey bee olfaction. These results indicated that the current ecotoxicological safety knowledge about chiral dinotefuran effects on honey bees should be amended.

Mass Spectrometric Quantification of Arousal Associated Neurochemical Changes in Single Honey Bee Brains and Brain Regions

Honey bee foragers show a strong diurnal rhythm of foraging activity, and such behavioral changes are likely under the control of specific neuromodulators. To identify and quantify neuromodulators involved in regulating rest and arousal in honey bees, we established a mass spectrometric method for quantifying 14 different neurochemicals and precursor molecules. We measured forager type and brain region specific differences in amine levels from individual honey bee brains and brain regions. The observed differences in amine levels between resting and aroused foragers resemble findings in other species indicating a conserved molecular mechanism by glutamate and GABA in regulating arousal. Subesophageal ganglion specific changes in the histaminergic system and global increases in aspartate during arousal suggest a possible role of histamine and aspartate in feeding and arousal, respectively. More aminergic systems were significantly affected due to arousal in nectar foragers than in pollen foragers, implying that forager phenotypes differ not only in their food preference but also in their neuromodulatory signaling systems (brain states). Finally, we found that neurotransmitter precursors were better at distinguishing brain states in the central brain, while their end products correlated with arousal associated changes in sensory regions like the optic and antennal lobes.

Immediate early genes in social insects: a tool to identify brain regions involved in complex behaviors and molecular processes underlying neuroplasticity

Social insects show complex behaviors and master cognitive tasks. The underlying neuronal mechanisms, however, are in most cases only poorly understood due to challenges in monitoring brain activity in freely moving animals. Immediate early genes (IEGs) that get rapidly and transiently expressed following neuronal stimulation provide a powerful tool for detecting behavior-related neuronal activity in vertebrates. In social insects, like honey bees, and in insects in general, this approach is not yet routinely established, even though these genes are highly conserved. First studies revealed a vast potential of using IEGs as neuronal activity markers to analyze the localization, function, and plasticity of neuronal circuits underlying complex social behaviors. We summarize the current knowledge on IEGs in social insects and provide ideas for future research directions.

Time-restricted foraging under natural light/dark condition shifts the molecular clock in the honey bee, Apis mellifera

Honey bees have a remarkable sense of time and individual honey bee foragers are capable of adjusting their foraging activity with respect to the time of food availability. Although, there is compelling experimental evidence that foraging behavior is guided by the circadian clock, nothing is known about the underlying molecular mechanisms. Here we present for the first time a study that explores whether time-restricted foraging under natural light-dark (LD) condition affects the molecular clock in honey bees. Food was presented in an enclosed flight chamber (12 m × 4 m × 4 m) either for 2 hours in the morning or 2 hours in the afternoon for several consecutive days and daily cycling of the two major clock genes, cryptochrome2 (cry2) and period (per), were analyzed for three different parts of the nervous system involved in feeding-related behaviors: brain, subesophageal ganglion (SEG), and the antennae with olfactory sensory neurons. We found that morning and afternoon trained foragers showed significant phase differences in the cycling of both clock genes in all three tissues. In addition, the phase differences were more pronounced when the feeder was scented with the common plant odor, linalool. Together our findings suggest that foraging time may function as a Zeitgeber that might have the capability to modulate the light entrained molecular clock.

Honey bee neurogenomic responses to affiliative and agonistic social interactions

Social interactions can be divided into two categories, affiliative and agonistic. How neurogenomic responses reflect these opposing valences is a central question in the biological embedding of experience. To address this question, we exposed honey bees to a queen larva, which evokes nursing, an affiliative alloparenting interaction, and measured the transcriptomic response of the mushroom body brain region at different times after exposure. Hundreds of genes were differentially expressed at distinct time points, revealing a dynamic temporal patterning of the response. Comparing these results to our previously published research on agonistic aggressive interactions, we found both shared and unique transcriptomic responses to each interaction. The commonly responding gene set was enriched for nuclear receptor signaling, the set specific to nursing was enriched for olfaction and neuron differentiation, and the set enriched for aggression was enriched for cytoskeleton, metabolism, and chromosome organization. Whole brain histone profiling after the affiliative interaction revealed few changes in chromatin accessibility, suggesting that the transcriptomic changes derive from already accessible areas of the genome. Although only one stimulus of each type was studied, we suggest that elements of the observed transcriptomic responses reflect molecular encoding of stimulus valence, thus priming individuals for future encounters. This hypothesis is supported by behavioral analyses showing that bees responding to either the affiliative or agonistic stimulus exhibited a higher probability of repeating the same behavior but a lower probability of performing the opposite behavior. These findings add to our understanding of the biological embedding at the molecular level.

Honey bee foraging induces upregulation of early growth response protein 1, hormone receptor 38 and candidate downstream genes of the ecdysteroid signalling pathway

In honey bees, continuous foraging at an artificial feeder induced a sustained upregulation of the immediate early genes early growth response protein 1 (Egr-1) and hormone receptor 38 (Hr38). This gene expression response was accompanied by an upregulation of several Egr-1 candidate downstream genes: ecdysone receptor (EcR), dopamine/ecdysteroid receptor (DopEcR), dopamine decarboxylase and dopamine receptor 2. Hr38, EcR and DopEcR are components of the ecdysteroid signalling pathway, which is highly probably involved in learning and memory processes in honey bees and other insects. Time-trained foragers still showed an upregulation of Egr-1 when the feeder was presented at an earlier time of the day, suggesting that the genomic response is more dependent on the food reward than training time. However, presentation of the feeder at the training time without food was still capable of inducing a transient increase in Egr-1 expression. Thus, learnt feeder cues, or even training time, probably affect Egr-1 expression. In contrast, whole brain Egr-1 expression changes did not differ between dancing and nondancing foragers. On the basis of our results we propose that food reward induced continuous foraging ultimately elicits a genomic response involving Egr-1 and Hr38 and their downstream genes. Furthermore this genomic response is highly probably involved in foraging-related learning and memory responses.

Deep evolutionary conservation of autism-related genes

E. O. Wilson proposed in Sociobiology that similarities between human and animal societies reflect common mechanistic and evolutionary roots. When introduced in 1975, this controversial hypothesis was beyond science's ability to test. We used genomic analyses to determine whether superficial behavioral similarities in humans and the highly social honey bee reflect common molecular mechanisms. Here, we report that gene expression signatures for individual bees unresponsive to various salient social stimuli are significantly enriched for autism spectrum disorder-related genes. These signatures occur in the mushroom bodies, a high-level integration center of the insect brain. Furthermore, our finding of enrichment was unique to autism spectrum disorders; brain gene expression signatures from other honey bee behaviors do not show this enrichment, nor do datasets from other human behavioral and health conditions. These results demonstrate deep conservation for genes associated with a human social pathology and individual differences in insect social behavior, thus providing an example of how comparative genomics can be used to test sociobiological theory.

Short reads from honey bee (Apis sp.) sequencing projects reflect microbial associate diversity

High throughput (or ‘next generation’) sequencing has transformed most areas of biological research and is now a standard method that underpins empirical study of organismal biology, and (through comparison of genomes), reveals patterns of evolution. For projects focused on animals, these sequencing methods do not discriminate between the primary target of sequencing (the animal genome) and ‘contaminating’ material, such as associated microbes. A common first step is to filter out these contaminants to allow better assembly of the animal genome or transcriptome. Here, we aimed to assess if these ‘contaminations’ provide information with regard to biologically important microorganisms associated with the individual. To achieve this, we examined whether the short read data from Apis retrieved elements of its well established microbiome. To this end, we screened almost 1,000 short read libraries of honey bee (Apis sp.) DNA sequencing project for the presence of microbial sequences, and find sequences from known honey bee microbial associates in at least 11% of them. Further to this, we screened ∼500 Apis RNA sequencing libraries for evidence of viral infections, which were found to be present in about half of them. We then used the data to reconstruct draft genomes of three Apis associated bacteria, as well as several viral strains de novo. We conclude that ‘contamination’ in short read sequencing libraries can provide useful genomic information on microbial taxa known to be associated with the target organisms, and may even lead to the discovery of novel associations. Finally, we demonstrate that RNAseq samples from experiments commonly carry uneven viral loads across libraries. We note variation in viral presence and load may be a confounding feature of differential gene expression analyses, and as such it should be incorporated as a random factor in analyses.

Identification of two novel biomarkers of rectal carcinoma progression and prognosis via co-expression network analysis

mRNA expression profiles provide important insights on a diversity of biological processes involved in rectal carcinoma (RC). Our aim was to comprehensively map complex interactions between the mRNA expression patterns and the clinical traits of RC. We employed the integrated analysis of five microarray datasets and The Cancer Genome Atlas rectal adenocarcinoma database to identify 2118 consensual differentially expressed genes (DEGs) in RC and adjacent normal tissue samples, and then applied weighted gene co-expression network analysis to parse DEGs and eight clinical traits in 66 eligible RC samples. A total of 16 co-expressed gene modules were identified. The green-yellow and salmon modules were most appropriate to the pathological stage (R = 0.36) and the overall survival (HR =13.534, P = 0.014), respectively. A diagnostic model of the five pathological stage hub genes (SCG3, SYP, CDK5R2, AP3B2, and RUNDC3A) provided a powerful classification accuracy between localized RC and non-localized RC. We also found increased Secretogranin III (SCG3) expression with higher pathological stage and poorer prognosis in the test and validation set. The increased Homer scaffolding protein 2 (HOMER2) expression with the favorable survival prediction efficiency significantly correlated with the markedly reduced overall survival of RC patients and the higher pathological stage during the test and validation set. Our findings indicate that the SCG3 and HOMER2 mRNA levels should be further evaluated as predictors of pathological stage and survival in patients with RC.