Tachykinin signaling inhibits task-specific behavioral responsiveness in honeybee workers
Abstract
Behavioral specialization is key to the success of social insects and leads to division of labor among colony members. Response thresholds to task-specific stimuli are thought to proximally regulate behavioral specialization, but their neurobiological regulation is complex and not well understood. Here, we show that response thresholds to task-relevant stimuli correspond to the specialization of three behavioral phenotypes of honeybee workers in the well-studied and important Apis mellifera and Apis cerana. Quantitative neuropeptidome comparisons suggest two tachykinin-related peptides (TRP2 and TRP3) as candidates for the modification of these response thresholds. Based on our characterization of their receptor binding and downstream signaling, we confirm a functional role of tachykinin signaling in regulating specific responsiveness of honeybee workers: TRP2 injection and RNAi-mediated downregulation cause consistent, opposite effects on responsiveness to task-specific stimuli of each behaviorally specialized phenotype but not to stimuli that are unrelated to their tasks. Thus, our study demonstrates that TRP signaling regulates the degree of task-specific responsiveness of specialized honeybee workers and may control the context specificity of behavior in animals more generally.
Introduction
Behavioral responses of animals to external and internal stimuli have evolved to optimize survival and reproduction under average circumstances (Darwin, 1859). However, environmental and interindividual variability commonly causes deviations from the average, resulting in selection for context-specific and condition-dependent behavior (Dall et al., 2004; Bolnick et al., 2003; West-Eberhard, 1989). Such plasticity is limited and evolutionary constraint (Arnold, 1992) of behavior occurs in the form of behavioral syndromes: differences among individuals that manifest across different contexts (Sih et al., 2004). Advantages of behavioral plasticity and specificity have been documented in many systems and some neuroendocrine mechanisms have been identified (O’Connell and Hofmann, 2011; Kim et al., 2017). However, neural mechanisms that allow the sophistication of behavioral repertoires by regulating the context specificity of behavioral responses remain insufficiently understood.
Behavioral regulation is particularly important in social species in which social interactions provide a high diversity of behavioral context (Oliveira, 2009; Gronenberg and Riveros, 2009). Social evolution also allows individuals to restrict their behavioral repertoires through temporal or permanent behavioral specialization (Oster and Wilson, 1978). This specialization and the resulting division of labor are believed to be major contributors to the successful colony life of many social insects despite its potential costs (Jeanson and Weidenmüller, 2014). Advanced social evolution thus allows interindividual plasticity to replace individual behavioral plasticity. Nevertheless, the principal problem of behavioral plasticity across different contexts remains the same, and social insects can be constrained in their behavioral evolution by correlated selection responses across different behaviors or castes (Rueppell et al., 2006; Jandt et al., 2014).
Behavior often occurs in response to a specific stimulus exceeding an individual’s specific response threshold (Tinbergen, 1951; Mayr, 1974). Response thresholds depend on internal physiological states (Ricklefs and Wikelski, 2002) and translate the value of perceived stimuli into probabilities of behavioral responses and vary among individuals (Scheiner et al., 2004). In social insects, individual variation in response thresholds is linked to division of labor (Theraulaz et al., 1998; Beshers et al., 1999; Robinson, 1992; Schulz et al., 2002), and numerous studies have characterized this link across multiple levels of biological organization (Scheiner et al., 2004; Page et al., 2012; Johnson, 2010). The Western honeybee (Apis mellifera) is the most-studied social insect model. Many aspects of its division of labor are driven by a lifelong behavioral ontogeny, leading to age polyethism (Seeley, 1982). Young bees perform numerous inside tasks, most prominently brood care in form of alloparental nursing behavior, before transitioning to a mix of other in-hive tasks (Johnson, 2008). Similar to the highly specialized nursing stage, the final behavioral state of older bees as outside foragers is almost exclusive of other tasks (Seeley, 1982). Moreover, foragers often specialize on collecting only one of the principal food sources, pollen or nectar (Page et al., 2000). Specialized nurse bees (NBs), nectar foragers (NFs), and pollen foragers (PFs) are also found in the closely related Eastern honeybee, Apis cerana (Ji et al., 2020), thus offering a suitable parallel system for our study.
The behavioral specialists of NBs, NFs, and PFs exhibit pronounced differences in their responsiveness to task-related stimuli. Responsiveness to brood pheromones peaks at typical nursing age (Pankiw, 2004). In contrast, foragers have a lower response threshold to sugars and light than nurses (Değirmenci et al., 2018; Ben-Shahar et al., 2003). Among foragers, pollen specialists exhibit higher responsiveness to sucrose and pollen stimuli than NFs (Page et al., 1998; Pankiw and Page, 1999). Response thresholds can be quantified based on the honeybees’ reflexive extension of their proboscis in response to stimuli, such as sucrose (Scheiner et al., 2004). The spontaneous proboscis extension reflex (PER) to sucrose has been expanded to other stimuli that bees spontaneously respond to Nicholls and de Ibarra, 2013; Zhang et al., 2020 and conditioned stimuli to which no spontaneous responses occur (Giurfa and Sandoz, 2012).
Response thresholds can be modified by biogenic amines, and dopamine, 5-hydroxy-tryptamine, octopamine, and tyramine have been implicated in the regulation of different behaviors of worker bees although the mechanisms are not entirely clear (Schulz et al., 2002; Scheiner et al., 2006). Few studies have addressed the role of neuropeptides although they are a diverse group of neurotransmitters that can also act as neurohormones on distal targets to coordinate a wide range of internal states and behavioral processes (Nässel, 2002). Neuropeptides are intimately involved in food perception and social interaction of insects (Schoofs et al., 2017), two processes that are central to division of labor in social insects (Ament et al., 2010). Neuropeptides mediate pheromonal effects on physiology (Shankar et al., 2015; Gendron et al., 2014) and usually exhibit a high degree of specificity (Inagaki et al., 2014; Taghert and Nitabach, 2012). Therefore, neuropeptides are prime candidates for mediating the independent adjustment of socially relevant response thresholds and they have been implicated in honeybee worker specialization and division of labor (Ji et al., 2020; Brockmann et al., 2009).
More than 100 mature neuropeptides derived from 22 protein precursors have been identified in the Western honeybee, A. mellifera (Han et al., 2015; Boerjan et al., 2010). Several neuropeptides, including allatostatin, leucokinin, and tachykinin-related peptides (TRPs), may be involved in the control of social behavior of A. mellifera and the closely related A. cerana, such as aggression (Pratavieira et al., 2018), foraging (Ji et al., 2020; Brockmann et al., 2009), brood care (Han et al., 2015), and possibly a wide array of other behaviors (Pratavieira et al., 2014). However, most of these results are based on correlations between behavior and neuropeptide expression, and more detailed studies are needed to understand the causal roles of neuropeptides in the behavioral specialization among honeybee workers.
Here, we report the results of a comprehensive study to test the hypothesis that neuropeptides regulate the division of labor in honeybees. We initially compared response thresholds to task-relevant stimuli among behaviorally defined worker groups (NBs, NFs, PFs) of the two closely related honeybee species A. mellifera and A. cerana. These response thresholds were correlated with neuropeptide expression levels, especially TRPs, suggesting a role of TRPs in worker specialization. Based on these results, we demonstrated in a series of TRP injections and RNAi-mediated knockdown of the TRP and its receptor (TRPR) a causal role of this pathway in modulating different response thresholds in a task-specific manner. We characterized the TRP signaling pathway molecularly by studying the TRPR properties and the important downstream events of cAMP and calcium accumulation and extracellular-signal-regulated kinase (ERK) activation in cell culture. Finally, we confirmed this connection in vivo by demonstrating corresponding effects of TRP injections and RNAi-mediated knockdown of TRP and TRPR on ERK phosphorylation.
Results
The task-specific responsiveness of worker bees showed significant variations between behavioral phenotypes and the two honeybee species
Based on the response threshold model for division of labor, we initially studied the responsiveness of three different behavioral phenotypes to three different task-relevant stimuli in A. mellifera and A. cerana using the PER assay. In accordance with our predictions, we identified significant differences among worker behavioral phenotypes in the PER responsiveness to the task-specific stimuli, including sucrose solution, pollen, and larva (Figure 1, Figure 1—source data 1, Supplementary file 1).
The percentage of bees showing a PER increased with sucrose concentration across all experimental groups (Figure 1A). In both, Apis mellifera ligustica (AML) and Apis cerana cerana (ACC), the sucrose response scores (SRSs) of PFs were higher than the SRSs of NFs (AML: Z = 7.0, p≤0.001; ACC: Z = 6.1, p<0.001) and nurse bees (NBs) (AML: Z = 5.9, p<0.001; ACC: Z = 5.2, p<0.001), while no significant difference between NFs and NBs was observed in either species (Figure 1B). PFs were more responsive than NFs and NBs to all sucrose concentrations. The species comparison between AML and ACC showed significant higher sucrose responsiveness in PFs of AML than in PFs of ACC (Z = 2.361, p=0.018), specifically at sucrose concentrations of 0.3% (χ2 = 4.1, p=0.042), 1.0% (χ2 = 5.2, p=0.001), 3.0% (χ2 = 8.4, p=0.023), and 10.0% (χ2 = 5.3, p=0.021). NFs of AML and ACC showed no significant difference in overall SRS, but NFs of AML were more responsive than NFs of ACC at sucrose concentrations of 0.3% (χ2 = 4.5, p=0.035), 1.0% (χ2 = 4.5, p=0.033), and 3.0% (χ2 = 4.0, p=0.046). There was no significant difference between NBs of AML and ACC in sucrose responsiveness.
In AML, PFs were more responsive to pollen stimulation than NFs (χ2 = 14.9, p=0.002) and NBs (χ2 = 20.2, p<0.001), while there were no significantly statistical differences between NFs and NBs. Likewise, PFs of ACC were more sensitive than NFs (χ2 = 6.0, p=0.015) and NBs (χ2 = 7.8, p=0.001) without a statistically significant difference between NFs and NBs. PFs of AML showed a significant higher pollen responsiveness than of ACC (χ2 = 4.9, p=0.031), with no significant species differences in NFs and NBs (Figure 1C).
In larva responsiveness assay, NBs of AML showed increased responsiveness to larva stimulation compared to PFs (χ2 = 7.2, p=0.006) and NFs (χ2 = 10.3, p=0.001). Likewise, NBs of ACC were more sensitive than PFs (χ2 = 4.2, p=0.013) and NFs (χ2 = 6.1, p=0.002). NBs of AML were significantly more sensitive to larvae (χ2 = 4.3, p=0.027) than NBs of ACC, with no significant species differences in PFs and NFs (Figure 1D). In sum, our results indicated that behavioral specialization among the workers of both honeybee species correspond to differences in the workers’ responsiveness to task-relevant stimuli. Correlative evidence linking different neuropeptides to foraging specialization in A. mellifera (Pratavieira et al., 2018) motivated us to further investigate whether brain neuropeptides could regulate the response thresholds.
Quantitative neuropeptidomics revealed consistent correlations of TRPs with behavior
To identify neuropeptides that potentially influence behavioral specialization, we compared the whole brain neuropeptidomes of NBs, PFs, and NFs of AML and ACC with liquid chromatography with tandem mass spectrometry (LC-MS/MS). The comparisons revealed numerous differences among experimental groups, but only two tachykinin-related peptides (TRP2 and TRP3) showed consistent patterns relating to the task-specific responsiveness of the experimental groups.
Overall, 132 unique neuropeptides derived from 23 neuropeptide families were identified in the brain of AML worker bees (Supplementary file 2). In the brain of ACC worker bees, 116 unique neuropeptides derived from 22 neuropeptide families were identified (Supplementary file 3).
Quantitative comparison among the three behavioral phenotypes of AML showed that 40 neuropeptides derived from 16 neuropeptide families were differentially expressed the brain (Figure 2 and Figure 2—source data 1). Among 19 differentially expressed neuropeptides between PFs and NFs, nine neuropeptides were upregulated in PFs and 10 were upregulated in NFs. Among 24 differentially expressed neuropeptides between PFs and NBs, 18 were upregulated in PFs and six were upregulated in NBs. Moreover, 21 differentially expressed neuropeptides were found between NFs and NBs, with 14 upregulated in PFs and seven upregulated in NBs. In ACC, 18 neuropeptides were differentially expressed between PFs and NFs, with nine upregulated in each group. Between PFs and NBs, 27 neuropeptides showed different expression levels: 20 were upregulated in PFs and seven were upregulated in NBs (Figure 2—source data 2). Twenty-five neuropeptides were differentially expressed between NFs and NBs, with 19 upregulated in NFs and six in NBs. In the species comparison between AML and ACC, the number of differentially expressed neuropeptides in NBs, PFs, and NFs was 13, 10, and 11, of which 7, 6, and 6 were upregulated in AML, respectively (Figure 2—source data 3). Thus, this global comparison of brain neuropeptide levels revealed numerous quantitative and qualitative differences. However, consistent differences between NBs, NFs, and PFs of both species were only found for TRP2 and TRP3, prompting their further study.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/64830/elife-64830-fig2-data1-v2.docx
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/64830/elife-64830-fig2-data2-v2.docx
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/64830/elife-64830-fig2-data3-v2.docx
TRP signaling inhibited task-specific responsiveness
Based on our behavioral and peptidomics results, we hypothesized tachykinin signaling to affect response thresholds and tested this prediction with a combination of gain- and loss-of-function experiments. We focused on TRP2 because it showed a stronger binding affinity to the receptor TRPR than TRP3 in our biochemical work (see 'TRP/TRPR signaling was found to affect the Gαq and Gαs pathways and trigger the ERK cascade' section).
TRP2 injection decreased task-specific responsiveness
Task-specific responsiveness of the different behavioral phenotypes (PFs, NFs, and NBs) was decreased by injection of TRP2 in a task-specific manner (Figure 3 and Figure 3—source data 1).
Injection of the TRP2 peptide significantly reduced the SRS of PFs (Z = 2.2, p=0.031), significantly reducing PER responses to all sucrose concentrations used. Similarly, NFs injected with TRP2 displayed significantly lower SRS than control-injected NFs (Z = 2.3, p=0.019), significantly reducing PER responses to all sucrose concentrations except 0.1% (Figure 3A and B, Figure 3—source data 1). In contrast, TRP2-injected NBs did not show significant responsiveness changes to sucrose relative to controls. For pollen stimulation, PFs showed significantly decreased responsiveness to pollen loads after TRP2 injection (χ2 = 6.7, p=0.017), while no significant effects were observed in PFs and NFs (Figure 3C). In the larval responsiveness assay, injection of TRP2 only significantly affected the responsiveness of NBs (χ2 = 6.1, p=0.001) but not NFs or PFs (Figure 3D).
Downregulation of TRP or TRPR increased task-specific responsiveness
The function of TRP/TRPR signaling on task-specific responsiveness was further confirmed by RNAi-mediated downregulation of TRP or TRPR that increased task-specific responsiveness, opposite to the results of the TRP2 injection.
Knockdown efficiencies were close to 60% for TRP and TRPR mRNA levels at 24 hr post-injection of the corresponding dsRNA (Supplementary file 4). Therefore, subsequent PER assays were performed 24 hr after dsRNA injection. Relative to control injections, knockdown of either TRP or TRPR significantly increased the SRS of NFs (dsTRP: Z = 2.4, p=0.049; dsTRPR: Z = 2.6, p=0.025), specifically increasing the responses of NFs to sucrose at concentrations of 0.1% (dsTRP:χ2 = 3.9, p=0.039; dsTRPR: χ2 = 4.9, p=0.023), 0.3% (dsTRP: χ2 = 5.3, p=0.018; dsTRPR: χ2 = 4.3, p=0.030), 1.0% (dsTRP: χ2 = 7.0, p=0.007; dsTRPR: χ2 = 6.6, p=0.009), and 3.0% (dsTRP: χ2 = 6.0, p=0.012; dsTRPR: χ2 = 7.4, p=0.006) (Figure 4A and B, Figure 4—source data 1, Supplementary file 5). Knockdown of TRP and TRPR did not significantly change the overall SRS of PFs and NBs, although it significantly increased the responses of PFs to sucrose at concentrations of 0.1% (dsTRP: χ2 = 4.4, p=0.029; dsTRPR: χ2 = 6.1, p=0.011), 0.3% (dsTRP: χ2 = 5.2, p=0.018; dsTRPR: χ2 = 6.0, p=0.011), and 1.0% (dsTRP: χ2 = 5.0, p=0.020; dsTRPR: χ2 = 4.7, p=0.025). Responses to pollen stimulation after dsRNA injection indicated that knockdown of either TRP or TRPR specifically increased the pollen responsiveness of PFs (dsTRP: χ2 = 6.5, p=0.018; dsTRPR: χ2 = 6.4, p=0.010), whereas the effects on NFs and NBs were not significant (Figure 4C). The responsiveness of NBs to larvae was significantly increased after gene knockdown of either TRP (χ2 = 4.4, p=0.029) or TRPR (χ2 = 4.8, p=0.023), but NFs and PFs were not affected (Figure 4D).
The experimental up- and downregulation of TRP signaling resulted in complementary changes in the responsiveness of our experimental groups in a specific way: Only the particular, increased responsiveness to task-relevant stimuli of each behavioral group was affected, suggesting that TRP signaling generally regulates the degree of behavioral specialization by moderating the responsiveness to task-specific stimuli in behavioral specialists.
TRP/TRPR signaling was found to affect the Gαq and Gαs pathways and trigger the ERK cascade
The action of most insect neuropeptides is mediated by G-protein-coupled receptors (GPCRs) that activate cAMP- and Ca2+-dependent pathways, such as ERK signaling to affect diverse biological processes (Hauser et al., 2006; He et al., 2014; Werry et al., 2005). However, the details of TRP signaling vary among insects (Birse et al., 2006; Poels et al., 2007) and have not been well studied in honeybees. To further support our behavioral studies by providing a plausible biochemical action of TRP signaling, we thus confirmed that the honeybee TRPR was localized in the cell membrane and specifically activated by TRP, triggering intracellular cAMP accumulation, Ca2+ mobilization, and ERK phosphorylation by dually coupling Gαs and Gαq signaling pathways. Because TRP2 and TRP3 are very similar and TRP2 displayed slightly higher binding affinity to their common receptor (see below), we used TRP2 only in a few of the below experiments.
The honeybee TRPR gene was cloned and expressed in the human embryonic kidney cells (HEK293) and the insect Spodoptera frugiperda pupal ovary cells (Sf21). Significant cell surface expression was observed by fluorescence microscopy (Figure 5A and Figure 5—figure supplement 1), revealing that the honeybee TRPR was exclusively localized in the cell membrane. Competitive binding assays with labeled TRP2 and TRP3 confirmed the predicted high affinity of the TRPR for both, although it was higher for TRP2 than for TRP3 (Figure 5—figure supplement 1, Figure 5—source data 1).
-
Figure 5—source data 1
- https://cdn.elifesciences.org/articles/64830/elife-64830-fig5-data1-v2.xlsx
-
Figure 5—source data 2
- https://cdn.elifesciences.org/articles/64830/elife-64830-fig5-data2-v2.xlsx
-
Figure 5—source data 3
- https://cdn.elifesciences.org/articles/64830/elife-64830-fig5-data3-v2.xlsx
-
Figure 5—source data 4
- https://cdn.elifesciences.org/articles/64830/elife-64830-fig5-data4-v2.xlsx
-
Figure 5—source data 5
- https://cdn.elifesciences.org/articles/64830/elife-64830-fig5-data5-v2.xlsx
The detected accumulation of intracellular cAMP concentration only in HEK293 cells transformed with TRPR (Figure 5B, Figure 5—source data 2) confirmed that TRP2 and TRP3 can activate TRPR and trigger cAMP signaling. This effect was confirmed for both cell types in further dose-response experiments and compared to other neuropeptides, including short neuropeptide F (SNF), pigment-dispersing hormone (PDH), and corazonin (CRZ), which did not induce any detectable cAMP accumulation (Figure 5—figure supplement 2, Figure 5—source data 2). Correspondingly, selective inhibition/activation experiments implicated the Gαs and Gαq (but not Gαi) subunits in this signaling mechanism (Figure 5—figure supplement 3, Figure 5—source data 3). Furthermore, intracellular Ca2+ mobilization was also found to result from TRP2 or TRP3 stimulation, dependent on Gαq (Figure 5C, Figure 5—figure supplement 4, and Figure 5—source data 4).
The final experiment of our in vitro studies demonstrated ERK phosphorylation in response to TRP/TRPR signaling by Western blot analysis. Treatment with different concentrations of TRP2 induced a transient, dose-dependent phosphorylation of ERK in both HEK293 (EC50 = 68.04 nM) and Sf21 (EC50 = 1.68 nM) cells (Figure 5D, Figure 5—figure supplement 5, Figure 5—source data 5). Further time-dependent analysis indicated that TRP2 elicited transient phosphorylation of ERK with maximal phosphorylation at 2 min and near basal levels by 90 min (Figure 5—figure supplement 5, Figure 5—source data 5). Moreover, specific inhibitors were used to elucidate TRP/TRPR signaling-mediated ERK activation in both HEK293 and Sf21 cells. Treatment with MEK inhibitor U0126, PKA inhibitor H89, and PKC inhibitor Go6983, respectively, led to a significant inhibition of TRP/TRPR-mediated ERK activation, whereas Gαi inhibitor pertussis toxin (PTX) had no effect, demonstrating that honeybee TRP/TRPR signaling dually coupled to Gαs and Gαq proteins to activate the ERK signaling pathway (Figure 5—figure supplement 5, Figure 5—source data 5).
Regulation of ERK activity by TRP/TRPR signaling was confirmed in vivo
To complement our finding that TRP/TRPR signaling activates ERK phosphorylation in cell culture, we used our in vivo manipulations of TRP signaling to confirm the link between TRP and ERK signaling in living honeybee workers. Western blot results confirmed that TRP/TRPR signaling triggers ERK signaling in vivo. The level of phosphorylated ERK significantly increased after injection of TRP2 peptide into NBs, PFs, and NFs (Figure 6A) and decreased after knockdown of the TRP or TRPR transcripts (Figure 6B, Figure 6—source data 1).
Discussion
Behavioral plasticity plays a central role in animal adaptation and modulating behavioral responsiveness to different stimuli and contexts is key to individual fitness. The success of social insects is partly due to their efficient division of labor, a form of behavioral plasticity among instead of within individuals. In this study, we demonstrated that the responsiveness to task-relevant stimuli correlates with behavioral specialization in two different honeybee species. Through parallel characterization of the neuropeptidome, we identified two tachykinin-related peptides (TRP2 and TRP3) as putative mechanism to adjust task-specific response thresholds and thus proximally guide division of labor. Subsequently, we characterized the molecular action of TRP2 and TRP3 in cell culture by verifying their binding to their membrane-bound receptor and demonstrating activation of multiple downstream signaling mechanisms. Finally, we verified causal involvement of TRP signaling in modulating task-specific behavioral response thresholds through complementary outcomes of TRP2 injection and RNAi-mediated knockdown of TRP and TRPR: while injection decreased task-specific responses, downregulation of TRP or TRPR increased the same specific responses. Thus, we present a mechanism that tunes the behavioral responsiveness of animals to specific stimuli compared to others. We use behaviorally specialized honeybee workers as models but hypothesize that this function of TRP signaling could be more widely conserved to adjust the context specificity of behavioral responses in animals.
Among all the signaling molecules in the nervous system, neuropeptides represent the largest and most diverse category and are crucial in orchestrating various biological processes and behavioral actions (Burbach, 2011; Hökfelt et al., 2000). Thus, we quantitatively compared the entire neuropeptidome among three behavioral worker phenotypes of AML and ACC without an a priori assumption. In addition to discovering several new neuropeptides from the ACC and AML brain, we identified TRP2 and TRP3 as candidates. TRPs have been associated with the modulation of appetitive olfactory sensation (Ki et al., 2015; Winther et al., 2006; Gui et al., 2017), foraging (Pratavieira et al., 2018), sex pheromone perception (Shankar et al., 2015), and aggression (Asahina et al., 2014). Particularly in honeybees, TRP is preferentially expressed in the mushroom body and some neurons scattered in the antennal and optic lobes (Takeuchi et al., 2004), and some expression has also been found in antennae (Jain and Brockmann, 2020). These expression patterns are consistent with our hypothesis that TRP signaling may be modulating behavioral responsiveness due to various inputs. Although specific knockdown of TRPs in Drosophila has revealed spatially variable effects (Nässel et al., 2019), the changes in response thresholds that we found are most likely due to systemic manipulation of TRP or TRPR levels. Global TRP expression in the honeybee worker brain also increases during the transition from nursing to foraging, further implicating it in the regulation of honeybee social behavior (Pratavieira et al., 2014; Takeuchi et al., 2003). The spatiotemporal expression of TRPR remains to be explored in A. mellifera and A. cerana, but the specific functions that have been linked to specific neurons in Drosophila (Ki et al., 2015) suggests that this information is critical for more detailed mechanistic models of the behavioral effects of tachykinin in honeybees.
In our study, only expression of TRP2 and TRP3 varied consistently among behavioral phenotypes of AML and ACC. In both species, TRP2 and TRP3 were most abundant in the brain of NFs, followed by PFs, and finally NBs. This is consistent with the very specific responsiveness of NBs to brood stimuli observed in our PER experiments, while the responsiveness of PFs and NFs was successively less specific: PFs responded specifically to two stimuli, while NFs did not show specifically strong responses to any stimuli. Moreover, the comparison between AML and ACC indicated higher TRP2 and TRP3 abundance in ACC in each behavioral phenotype, commensurate with the less specific PER responsiveness in ACC compared to AML. A few other neuropeptides, such as apidaecins, diuretic hormone, and prohormone-3, showed somewhat similar expression patterns in both species, but none of these was as tightly correlated to behavioral responsiveness and none has previously been connected with behavioral regulation in insects or other animals. Therefore, the TRPs were chosen as candidates of the control of honeybee division of labor for subsequent functional tests and molecular characterization.
The action of most insect neuropeptides is mediated by binding to GPCRs and often involves cAMP and Ca2+ as second messengers (Hauser et al., 2006). The TRPR is activated by TRPs triggering intracellular cAMP accumulation and Ca2+ mobilization in fruit flies and silkworms (Bombyx mori) (He et al., 2014; Birse et al., 2006), while no cAMP responses were discovered in stable flies (Stomoxys calcitrans) (Poels et al., 2007). The results of our peptide-based binding assays functionally confirmed that the honeybee TRPR is indeed the receptor for TRP2 and TRP3. The subsequent functional assays revealed that TRP signaling results in a dose-dependent increase in both intracellular cAMP and Ca2+. Together, these results indicate that TRPs can activate TRPR and trigger second messengers to regulate downstream functions. TRP2 displayed a higher affinity to TRPR and induced higher cAMP and Ca2+ signaling than TRP3, leading us to focus on TRP2 in the later in vivo experiments. Moreover, TRP signaling is sensitive to Gαs activation and is significantly blocked by Gαq and PKA inhibitors, suggesting both Gαs and Gαq are involved in TRP signaling in honeybees. Many GPCRs are able to induce mitogen-activated protein kinase cascades via cooperation of Gαs, Gαq, and Gαi signals, leading to the phosphorylation of ERK1/2, which plays critical roles in diverse biological processes (Rozengurt, 2007). Our results indicate that honeybee TRP signaling mediates phosphorylation of ERK1/2 in a dose- and time-dependent manner in both HEK293 and Sf21 cells. In addition, ERK1/2 activation was significantly inhibited by the PKA, PKC, and MEK inhibitors, which is in line with the observation of intracellular cAMP accumulation and Ca2+ mobilization. Thus, honeybees seem to be very similar to silkworms with regard to the involvement of the Gαs/cAMP/PKA and Gαq/Ca2+/PKC signaling pathways in the regulation of TRP-induced ERK1/2 activation (He et al., 2014). Taken together, our results demonstrate that the honeybee TRPR is specifically activated by TRPs, eliciting intracellular cAMP accumulation, Ca2+ mobilization, and ERK phosphorylation by dually coupling Gαs and Gαq signaling pathways.
Our in vitro and in vivo demonstrations that TRP signaling activates the ERK pathway indicated that the basic signaling mechanisms are conforming to the general patterns in insects (Hauser et al., 2006; He et al., 2014; Werry et al., 2005). Thus, our behavioral results may extend to other species. ERK links TRP signaling also to the insulin/insulin-like signaling (IIS) pathway. IIS is controlled by neuropeptides through ERK in Drosophila (Lee et al., 2008), and this connection in honeybees ties TRP back to the age-based division of labor among workers: IIS influences the timing of the behavioral maturation of honeybee workers and brain AmIlp1 is significantly higher expressed in foragers than nurses (Ament et al., 2008), consistent with our finding that TRPs are higher in foragers than nurses. Numerous other physiological changes accompany the transition from in-hive nurses to foragers (Scheiner et al., 2006; Robinson, 1987; Wang et al., 2012; Toth and Robinson, 2005, and our results integrate TRPs as the most important neuropeptides into the regulation of the behavioral ontogeny of honeybee workers and potential feedback loops to the modulation of behavioral response thresholds. The specialization of NFs and PFs has also been linked to IIS (Hunt et al., 2007; Wang et al., 2009) and explained by differences in sucrose response thresholds (Pankiw and Page, 2000). Our findings here may connect the differences in response thresholds and IIS mechanistically through the TRP and ERK signaling pathways.
The PER paradigm is well suited to test behavioral response thresholds and has been used for over 50 years in honeybees (Giurfa and Sandoz, 2012). Consistent with previous studies, we found PFs to be more responsive to sucrose than NFs and NBs in A. mellifera (Page et al., 1998). Moreover, we found corresponding differences between these behavioral groups in the closely related A. cerana. The PF’s responsiveness to low sucrose concentrations might also make them more responsive to pollen, but the causation of the PER to pollen is unclear (Grüter et al., 2008) and other components of pollen may functionally distinguish pollen from sucrose responsiveness (Nicholls and de Ibarra, 2013). Our results support the view that pollen and sucrose are distinct stimuli: While our experimental manipulations of TRP signaling altered the responsiveness of PFs to pollen and sucrose, only responsiveness to sucrose was affected in NFs and only responsiveness to larvae was affected in NBs. The functional significance of the PER in response to larvae is currently unclear, but we show that it is specific to nurses and it has previously been linked to brood provisioning (Zhang et al., 2020). Thus, our diverse PER results in two species comprehensively support the hypothesis that task-specific response thresholds guide behavioral specialization, leading to division of labor among honeybee workers (Theraulaz et al., 1998; Beshers et al., 1999; Robinson, 1992).
TRPs may adjust specific sensory neural circuits, potentially acting in concert with other neuromodulators (Kahsai et al., 2010; Jung et al., 2013). However, we have currently no evidence to support the hypothesis of different molecular TRP actions in different stimulus-response pathways (cf Ki et al., 2015). Thus, we have to restrict our conclusion to the simple hypothesis that TRP signaling decreases task-specific response thresholds of behavioral specialists without affecting non-task-related thresholds: It decreases pollen and sucrose responsiveness in PFs, sucrose responsiveness in NFs, and responsiveness to larvae in NBs. TRP signaling in honeybees may thus be a general regulator of how task-specific stimuli are weighted relative to others and consequently how specialized behavioral specialists are (Figure 7). This effect translates to different degrees of division of labor in social insect colonies and may control the context specificity of behavioral responses in animals more generally (Sih et al., 2004).
Although AML and ACC are close relatives with similar basic biology, some behavioral differences have evolved since their speciation (Oldroyd and Wongsiri, 2009). AML and ACC share the age-based division of labor, with younger bees specializing on nursing before maturing to foraging activities (Hepburn and Radloff, 2011) and ACC foragers also specialize in nectar or pollen collection (Ji et al., 2020; Rueppell et al., 2008) similar to AML (Page et al., 2000). Accordingly, we found the main differences of stimulus responsiveness and TRPs' expression among worker phenotypes conserved. However, ACC exhibited less responses to the task-specific stimuli than AML. Consistent PER differences in AML and ACC between NFs and PFs and a generally lower responsiveness of ACC have been identified before (Yang et al., 2013), but the biological interpretation has remained unclear. It is possible that the species differences arise due to methodological bias, favoring AML performance in PER assays. However, our study offers the alternative explanation that ACC workers are less specialized than AML workers due to higher TRP signaling. Lower innate specialization may accompany better learning of ACC (Qin et al., 2012), facilitating its more opportunistic worker task allocation and resource exploitation than AML (Tan et al., 2008). These alternative life history strategies are plausible, given the typical differences in colony size and habitat (Jung et al., 2013; Sih et al., 2004; Ruttner, 1988). All three worker phenotypes of ACC exhibited higher levels of TRPs than their AML counterparts, but functional verification at the level of colony phenotypes will be required to unambiguously link TRP signaling to such interspecific differences in life history.
Materials and methods
Honeybee sources and sampling
Request a detailed protocolTwo honeybee species, AML and ACC, were maintained in the apiary of the Institute of Apicultural Research at the Chinese Academy of Agricultural Sciences in Beijing. Three colonies of each species with mated queens of identical age were selected as experimental colonies, and before experiments the colonies were equalized in terms of adult bee population, brood combs, and food storage. Frames containing old pupae (1–2 days before emergence) were put into an incubator (34°C and 80% relative humidity) for eclosion. Newly emerged worker bees were paint-marked on their thoraxes and placed back into their parent colonies. Ten days later, marked bees that had their head and thorax in open brood cells while contracting their abdomen for more than 10 s were collected as NBs. Twenty days after eclosion, marked bees were collected during early morning (between 8:00 a.m. and 10:00 a.m.) in good weather conditions during the blooming period of black locusts (Robinia pseudoacacia L.) as forager bees. The entrance to the hives were blocked to facilitate collecting. Bees flying into the hive with pollen loads were collected as PFs, returning foragers without pollen loads were collected as NFs. The experimental design of six groups (three behavioral phenotypes in two species) was used to compare responsiveness to task-specific stimuli ('Comparative PER experiments' section) and to relate these phenotypes to differences in the brain neuropeptidome ('Quantitative comparisons of brain neuropeptidomes' section).
Comparative PER experiments
Request a detailed protocolTo investigate the responsiveness of different worker bee behavioral phenotypes (NBs, PFs, and NFs of AML and ACC) to different stimulus modalities (sucrose solution, pollen, and larva), series of PER experiments were performed. One hundred bees of each behavioral phenotype were collected from each experimental colony in the morning, transferred to the laboratory and narcotized on ice, then harnessed using a previously described protocol (Wang and Tan, 2014). All harnessed bees were fed to satiation with 50% sucrose solution and placed in a dark incubator (20°C and 65% relative humidity) overnight. After 24 hr, all surviving bees were assayed for their PER following the methodology of Page et al., 1998. Each stimulus was assessed independently with a new set of bees.
To investigate the sucrose responsiveness, bees were assayed using an ascending order of sucrose concentrations: 0.1%, 0.3%, 1%, 3%, 10%, and 30% (weight/weight). A small droplet of each solution was touched to the bees’ antennae for 3 s and the extension of the proboscis was monitored during this time. The interval between each sucrose solution trial was 5 min to exclude sensitization or habituation effects. The total number of PER responses after stimulation with the six different sucrose concentrations was combined into an SRS of a bee (Pankiw et al., 2001; Scheiner et al., 2002; Scheiner et al., 2003). The SRSs of the three behavioral phenotypes in the same species were compared using Kruskal-Wallis tests with Bonferroni correction. Pairwise Mann-Whitney U tests were used for comparing the same phenotype from two honeybee species. The sucrose responsiveness for specific sucrose concentrations was further compared between different groups with independent chi-square tests.
To test pollen stimulation, fresh pollen loads that had been removed from the leg of randomly selected PFs of the test group were used: AML were tested with pollen collected by AML foragers and ACC with pollen collected by ACC foragers. These loads contained a mixture of different pollen, predominated by black locust (R. pseudoacacia). As a control for mechanical stimulation, each bee had both antennae first touched with a piece of filter paper and spontaneous responders were excluded. Subsequently, both antennae of each bee were gently touched with a pollen load and PER responses were recorded. The pollen responsiveness was compared with independent chi-square tests between different groups.
To test responsiveness to larva, 1-day-old larvae from each honeybee species were collected, briefly rinsed in distilled water to remove royal jelly residue and dried on a filter paper. As before, both antennae of bees were touched with a piece of filter paper first and spontaneous responders were excluded, then PERs in response to a larva touching the antennae were recorded. The responsiveness to larvae was compared with independent chi-square tests between different groups. Statistical analyses were conducted using SPSS Statistics 20.0 (IBM, USA).
Quantitative comparisons of brain neuropeptidomes
Request a detailed protocolTo explore brain neuropeptide functions in behavioral regulation, a label-free quantitative strategy was employed to compare neuropeptidomic variations between behavioral phenotypes and the two honeybee species. Three independent biological replicate samples (120 bees per sample) of NBs, PFs, and NFs of both AML and ACC (18 samples total) were collected and immediately frozen in liquid nitrogen. Individual brains were carefully dissected from the head capsule while remaining chilled on ice, and the dissected brains were frozen at −80°C until neuropeptide extraction.
The brains were homogenized at 4°C by using a 90:9:1 solution of methanol, H2O, and acetic acid. The homogenates were centrifuged at 12,000 g for 10 min at 4°C. The resulting supernatant containing the neuropeptides was collected and dried. The extracted neuropeptide samples were dissolved in 0.1% formic acid in distilled water, and the peptide concentration was quantified using a Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, USA). LC-MS/MS analysis was performed on an Easy-nLC 1200 (Thermo Fisher Scientific) coupled Q-Exactive HF mass spectrometer (Thermo Fisher Scientific). Buffer A (0.1% formic acid in water) and buffer B (0.1% formic acid in acetonitrile) were used as mobile phase buffers. Neuropeptides were separated using the following gradients: from 3% to 8% buffer B in 5 min, from 8% to 20% buffer B in 80 min, from 20% to 30% buffer B in 20 min, from 30% to 90% buffer B in 5 min, and remaining at 90% buffer B for 10 min. The eluted neuropeptides were injected into the mass spectrometer via a nano-ESI source (Thermo Fisher Scientific). Ion signals were collected in a data-dependent mode and run with the following settings: full scan resolution at 70,000, automatic gain control (AGC) target 3 × 106; maximum inject time (MIT) 20 ms; scan range m/z 300–1800; MS/MS scans resolution at 17,500; AGC target 1 × 105; MIT 60 ms; isolation window 2 m/z; normalized collision energy 27; loop count 10; and dynamic exclusion:charge exclusion: unassigned, 1, 8, >8; peptide match preferred; exclude isotopes on; dynamic exclusion: 30 s. Raw data were retrieved using Xcalibur 3.0 software (Thermo Fisher Scientific).
The extracted MS/MS spectra were searched against a composite database of A. mellifera (23,491 protein sequences, downloaded from NCBI on July 2018) or A. cerana (20,934 protein sequences, downloaded from NCBI on July 2018) using in-house PEAKS 8.5 software (Bioinformatics Solutions, Canada). Amidation (A, −0.98) and pyro-glu from Q (P, −17.03) were selected as variable modifications. The other parameters used were: parent ion mass tolerance, 20.0 ppm; fragment ion mass tolerance, 0.05 Da; enzyme, none; maximum allowed variable PTM per peptide, 2. A fusion target-decoy approach was used for the estimation of the false discovery rate and controlled at ≤1.0% (−10 log p≥20.0) both at protein and peptide levels. Neuropeptide identifications were only used if ≥2 spectra were identified in at least two of the three replicates of each sample type.
Quantitative comparison of brain neuropeptidomes was performed by the label-free approach in PEAKS Q module. Feature detection was performed separately on each sample by using the expectation-maximization algorithm. The features of the same peptide from different samples were reliably aligned together using a high-performance retention time alignment algorithm (Lin et al., 2013). Peptide features were considered significantly different between experimental groups if pairwise p<0.01 and fold change ≥1.5. A heat map of differentially expressed proteins was created by Gene cluster 3.0 using the unsupervised hierarchical clustering, and the result was visualized using Java Tree view software. The LC−MS/MS data and search results are deposited in ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository with the dataset identifier PXD018713.
Characterization of honeybee TRP signaling pathway
To characterize honeybee TRP signaling pathway, the TRPR gene was first cloned and expressed in human and insect cell lines to identify its cellular location and verify its binding to TRPs. Additionally, these cells were used to test whether TRP/TRPR signaling triggers intracellular cAMP accumulation, Ca2+ mobilization, and ERK phosphorylation.
TRPR gene clone and expression
Request a detailed protocolTo amplify the full-length sequence encoding TRPR of A. mellifera, primers were designed using Primer Premier 5.0 software (PREMIER Biosoft, USA) based on the sequence from GenBank KT232312. The coding sequence of TRPR was amplified and cloned into FLAG-tag expression vectors (pCMV-FLAG and pBmIE1-FLAG) and EGFP-tag expression vectors (pEGFP-N1 and pBmIE1-EGFP). The primers used are documented in Supplementary file 6. All constructs were sequenced to verify the correct sequence, orientation, and reading frame of the inserts.
The human embryonic kidney cell line HEK293 and the insect S. frugiperda pupal ovary cell line Sf21 were used for honeybee TRPR expression. HEK293 cells (RRID:CVCL_0045) were purchased from American Type Culture Collection (ATCC, CRL-1573, the identity has been authenticated using STR profiling) and cultured in DMEM medium (Gibco, USA) supplemented with 10% fetal bovine serum (FBS). Cells were routinely tested for mycoplasma contamination every 6 months. Sf21 cells (RRID: CVCL_0518) were purchased from Thermo Fisher Scientific and were cultured in TC100 medium (Gibco) supplemented with heat-inactivated 10% FBS. Cells were routinely tested for mycoplasma contamination every 6 months. Transfection of HEK293 cells was performed using Lipo6000 transfection reagent (Beyotime, China), while transfection of Sf21 cells was performed using LipoInsect transfection reagent (Beyotime), according to the manufacturer’s instructions.
Cellular location of TRPR
Request a detailed protocolTo confirm the location of the honeybee TRPR, receptor surface expression assays were performed. HEK293 or Sf21 cells expressing TRPR-EGFP were seeded onto poly-L-lysine-coated glass coverslips and allowed to attach overnight under normal growth conditions. After 24 hr, cells were incubated with the membrane probe DiI (Beyotime) and the nucleic acid probe Hoechst 33342 (Beyotime) at 37°C for 10 min, then fixed with 4% paraformaldehyde for 15 min. Cells transfected with empty EGFP-tag expression vectors were used as a control. The cells were imaged using a Leica SP8 (Leica Microsystems, Germany) confocal microscope equipped with an HC PL APO CS2 63×/1.40 oil objective. Images were acquired with the sequence program in the Leica LAS X software.
Binding of TRPs to TRPR
Request a detailed protocolTo confirm the direct binding of the honeybee TRPs to TRPR, competitive binding experiments were performed using synthesized TAMRA-TRP2 (TAMRA-ALMGFQGVRa) and TAMRA-TRP3 (TAMRA-APMGFQGMRa), with TAMRA labeled at the N-terminus. The neuropeptides used as ligands here and in later sections were commercially synthesized by SynPeptide Co, Ltd (China). All peptides were purified by reverse-phase high-performance liquid chromatography with a purity >98%, lyophilized, and diluted to the desired concentrations for subsequent experiments. The peptide sequences were verified by us using a Q-Exactive HF mass spectrometer (Thermo Fisher Scientific).
HEK293 and Sf21 cells expressing FLAG-TRPR were first seeded onto poly-L-lysine-coated 96-well plates and cultured overnight. On the next day, cells were washed once with phosphate-buffered saline (PBS), then incubated with 25 ml TAMRA-TRP2 or TAMRA-TRP3 (10 nM) in the presence of increasing concentrations of unlabeled TRP2 and TRP3 in a final volume of 100 ml of binding buffer (PBS containing 0.2% bovine serum albumin). Cells were incubated at room temperature for 2 hr. Fluorescence intensity was measured with a fluorescence spectrometer microplate reader (Tecan Infinite 200 PRO, Tecan, Switzerland) after washing twice with binding buffer. The cells transfected with empty FLAG-tag expression vectors were used as a control. The binding displacement curves were analyzed by GraphPad Prism 8.0 (GraphPad Software, USA) using the non-linear logistic regression method.
TRP/TRPR signaling targets: cAMP, Ca2+, and ERK
Request a detailed protocolTo test whether TRP/TRPR signaling affects cAMP accumulation, intracellular cAMP was measured after incubation of HEK293 and Sf21 cells expressing FLAG-TRPR and pCRE-Luc with TRP2 and TRP3. After seeding in a 96-well plate overnight, HEK293 or Sf21 cells co-transfected with pFLAG-TRPR and pCRE-Luc were grown to about 90% confluence. After washing once with PBS, cells were incubated with the neuropeptides TRP2, TRP3, SNF, PDH, and CRZ in serum-free medium for 4 hr at 37°C for HEK293 cells, and at 28°C for Sf21 cells. Cells transfected with empty EGFP-tag expression vectors were used as a control. Luciferase activity was detected by a luciferase assay system (Promega, USA). Fluorescence intensity was measured with a Tecan fluorescence spectrometer. When characterizing the TRP-mediated cAMP accumulation, cells were pretreated with Gαi inhibitor PTX, Gαs activator cholera toxin (CTX), Gαq inhibitor YM-254890, and PKA inhibitor H89 before stimulation with TRP2.
To test whether TRP signaling also affects intracellular Ca2+ concentrations, intracellular Ca2+ was measured after incubation of HEK293 and Sf21 cells expressing FLAG-TRPR with TRP2 or TRP3. Cells were detached by a non-enzymatic cell dissociation solution (Sigma-Aldrich, USA), washed twice with PBS, and resuspended at a density of 5 × 106 cells/ml in HEPES buffered saline (Macklin, China). Cells were then incubated with 3 μM Fura-2 AM (MedChemExpress, USA) for 30 min at 37°C for HEK293 cells, and at 28°C for Sf21 cells. Intracellular Ca2+ flux was measured using excitation wavelengths alternating at 340 and 380 nm with emission measured at 510 nm in a Tecan fluorescence spectrometer. When characterizing the detailed TRP-mediated intracellular Ca2+ mobilization, cells were pretreated with Gαq inhibitor YM-254890 and PLC inhibitor U73122 before stimulation with TRP2.
To assess whether TRP signaling mediates ERK1/2 signaling, ERK1/2 phosphorylation was measured by Western blot analysis after incubation of HEK293 and Sf21 cells expressing FLAG-TRPR with TRP2. Cells were seeded in 24-well plates and starved for 4 hr in serum-free medium to reduce background ERK1/2 activation and eliminate the effects of the change of medium. After incubation with TRP2, cells were lysed by RIPA buffer (Beyotime) at 4°C for 30 min. Protein concentration was determined according to the Bradford method using BSA as the standard and the absorption was measured at 595 nm (spectrophotometer DU800, Beckman Coulter, USA), then all the samples were kept in −80°C for further use. For Western blot, equal amounts of total cell lysate (20 μg/lane) were fractionated by SDS-PAGE (10%) and transferred to a PVDF membrane (Millipore, USA) using an iBlot dry blotting system (Invitrogen, USA). The membranes were blocked for 2 hr at room temperature and then incubated with rabbit monoclonal anti-pERK1/2 antibody (Cell Signaling Technology, USA) and anti-rabbit horseradish peroxidase-conjugated secondary antibody (Cell Signaling Technology) according to the manufacturers’ protocols. Antibody reactive bands were visualized using Pierce ECL western blotting substrate (Thermo Fisher Scientific, USA) followed by photographic film exposure. Total ERK1/2 was assessed as a loading control after p-ERK1/2 chemiluminescence detection. Quantification analyses were performed using Gel-Pro Analyzer 4.0 software (Media Cybernetics, USA).
To explore the detailed TRP-mediated ERK1/2 signaling, cells were pretreated with Gαi inhibitor PTX, MEK inhibitor U0126, PKA inhibitor H89, and PKC inhibitor Go6983 before stimulation with TRP2.
Effects of TRP2 injection on task-specific responsiveness
Request a detailed protocolTo confirm the function of TRP on task-specific responsiveness, NBs, PFs, and NFs of AML were injected with TRP2 and tested for their PER response to sucrose solution, pollen, and larva. About 150 bees of each behavioral phenotype were collected in the morning, then harnessed, fed and placed in a dark incubator as described in 'Comparative PER experiments' section. After 24 hr, all surviving bees were evenly divided into two groups and injected with 1 μl TRP2 solution (1 μg/μl, synthesized TRP2 dissolved in ddH2O) or 1 μl of ddH2O into the head of honeybees via the central ocellus using a glass capillary needle coupled to a microinjector. Bees injected with ddH2O were used as control. All injected bees were put back to the dark incubator and 1 hr after injection all surviving bees were assayed for their PER to stimulations of sucrose solution, pollen, and larva as described in 'Comparative PER experiments' section. Each experiment was performed with a new set of bees containing about 55 individuals per experimental and control group.
The average SRSs of the TRP2 injection group and the ddH2O injection group were compared separately for each of the three behavioral phenotypes (NBs, PFs, and NFs) using pairwise Mann-Whitney U tests. The sucrose responsiveness was further compared between different groups at each specific sucrose concentration with independent chi-square tests. The responsiveness to pollen and larvae was compared between TRP2 injection group and ddH2O injection group with independent chi-square tests for each behavioral phenotype separately. All statistical analyses were performed with SPSS Statistics 20.0 (IBM).
Effects of RNAi-mediated downregulation of TRP or TRPR on responsiveness
Request a detailed protocolTo further confirm the hypothesized effects of TRP/TRPR signaling on task-specific responsiveness, RNAi-mediated downregulation of TRP and TRPR was performed on NBs, PFs, and NFs of AML and then their PER to sucrose solution, pollen, and larva was compared to controls.
Before evaluating the behavioral effects of transcript knockdown of TRP or TRPR, preliminary experiments were performed to test the dsRNA-mediated knockdown efficiencies of TRP and TRPR. The dsRNAs of the TRP and TRPR genes were prepared using the T7 RiboMAX Express RNAi system (Promega). The primers used are listed in Supplementary file 6. Sixty bees were randomly collected from each of the three AML colonies. Bees were harnessed, fed with sucrose, and put into the dark incubator (20°C and 65% relative humidity) to acclimatize to the experimental conditions. After 30 min, dsRNA (200 ng/bee for TRP, 2 μg/bee for TRPR) was microinjected into the head of honeybees via the central ocellus using a glass capillary needle-coupled microinjector. dsRNA of green fluorescent protein gene (dsGFP, 2 μg/bee) was used as control in all RNAi experiments. All harnessed bees were fed with 50% sucrose solution every 12 hr. At 0, 12, 24, and 48 hr after injection, a group of six individual bees were collected from each injection group (dsTRP, dsTRPR, and dsGFP) for comparing TRP and TRPR expression. Individual brains were carefully dissected and frozen at −80°C until RNA extraction. Three independent replicate groups per condition were collected and qRT-PCR was performed to calculate the RNAi efficiency. Total RNA was isolated using TRIzol reagent (Takara, Japan). Total RNA quantification was performed by NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific), and the quality of RNA was evaluated by 1.0% denaturing agarose gel electrophoresis. Reverse transcription was performed using a PrimeScript RT reagent kit (Takara), according to the manufacturer’s instructions. Gene-specific mRNA levels were assessed by qPCR using TB Green Fast qPCR Mix (Takara) on a LightCycler 480II instrument (Roche, Switzerland). The β-actin gene was used as a reference gene. After verifying amplification efficiency of the selected genes and β-actin (from 96.8% to 100.5%), the differences in gene expression levels were calculated using the 2−ΔΔCt method. Pairwise differences in gene expression were considered significant at p<0.05, using one-way ANOVA (SPSS Statistics 20.0). The primers used for qPCR are shown in Supplementary file 6.
After determination of knockdown efficiencies (see 'Results' section), 24 hr post-injection was chosen as the time point to study the PER effects of dsRNA-mediated knockdown of TRP and TRPR. About 200 bees of each behavioral phenotype (NBs, PFs, and NFs of AML) were collected in the morning, harnessed, and remained in a dark incubator to acclimatize. After 30 min, all surviving bees of each behavioral phenotype were evenly divided into three groups, injected with dsTRP, dsTRPR, and dsGFP, and kept as described above. After 24 hr, all surviving bees were assayed for their PER to stimulations of sucrose solution, pollen, or larvae as described in 'Comparative PER experiments' section. Each stimulus was assessed with a new set of bees containing about 55 individuals for each treatment group (dsTRP, dsTRPR, and dsGFP). The SRSs of the TRP-knockdown, TRPR-knockdown, and control groups were compared using Kruskal-Wallis tests with Bonferroni correction for each behavioral phenotype separately. The sucrose responsiveness was further compared between the different groups at the same sucrose concentration with independent chi-square tests. The responsiveness to pollen and larvae was compared between the TRP knockdown, TRPR knockdown, and control groups using independent chi-square tests for each behavioral phenotype separately. All statistical analyses were performed with SPSS Statistics 20.0 (IBM).
Effects of TRP2 injection and RNAi-mediated downregulation of TRP and TRPR on ERK signaling in honeybee workers
Request a detailed protocolTo test whether manipulating TRP/TRPR signaling has effect on honeybee ERK signaling, a group of 10 individual worker bees were collected from each injection group (ddH2O, TRP2, dsTRP, dsTRPR, and dsGFP) to compare ERK phosphorylation levels. Three independent replicate groups per condition were collected and Western blot analyses were performed: Honeybee brains were carefully dissected and frozen at −80°C until protein extraction. Brain protein extractions were carried out according to our previously described method with some modifications. Briefly, the larvae were homogenized with lysis buffer (LB, 8 M urea, 2 M thiourea, 4% CHAPS, 20 mM Tris-base, 30 mM dithiothreitol). The mixture was homogenized for 30 min on ice and sonicated 20 s per 5 min during this time, then centrifuged at 12,000 g and 4°C for 10 min. Ice-cold acetone was added to the collected supernatants, and then the mixture was kept on ice for 30 min for protein precipitation. Subsequently, the mixture was centrifuged at 12,000 g and 4°C for 10 min. The supernatant was discarded and the pellets were resolved in LB and kept at −20°C for further use. Western blot analyses were performed as described in 'TRP/TRPR signaling targets: cAMP, Ca2+, and ERK' section.
Data availability
Original data have been deposited to ProteomeXchange Consortium with the dataset identifier PXD018713 under http://proteomecentral.proteomexchange.org or are provided as supplementary data files.
-
PRIDEID PXD018713. Honeybee brain neuropeptides LC-MSMS.
References
-
Nutritional regulation of division of labor in honey bees: toward a systems biology perspectiveWiley Interdisciplinary Reviews: Systems Biology and Medicine 2:566–576.https://doi.org/10.1002/wsbm.73
-
Constraints on phenotypic evolutionThe American Naturalist 140:S85–S107.https://doi.org/10.1086/285398
-
cGMP-dependent changes in phototaxis: a possible role for the foraging gene in honey bee division of laborJournal of Experimental Biology 206:2507–2515.https://doi.org/10.1242/jeb.00442
-
BookResponse thresholds and division of labor in insect coloniesIn: Detrain C, Deneubourg J. L, Pasteels J. M, editors. Information Processing in Social Insects. Basel: Springer. pp. 115–139.https://doi.org/10.1007/978-3-0348-8739-7
-
The ecology of individuals: incidence and implications of individual specializationThe American Naturalist 161:1–28.https://doi.org/10.1086/343878
-
Open source clustering softwareBioinformatics 20:1453–1454.https://doi.org/10.1093/bioinformatics/bth078
-
BookSocial brains and behavior: past and presentIn: Gadau J, Fewell J, editors. Organization of Insect Societies: From Genome to Sociocomplexity. Camridge MA: Harvard University Press. pp. 377–401.
-
Role of a tachykinin-related peptide and its receptor in modulating the olfactory sensitivity in the oriental fruit fly, bactrocera dorsalis (Hendel)Insect Biochemistry and Molecular Biology 80:71–78.https://doi.org/10.1016/j.ibmb.2016.12.002
-
Neuropeptides--an overviewNeuropharmacology 39:1337–1356.https://doi.org/10.1016/S0028-3908(00)00010-1
-
Behavioral genomics of honeybee foraging and nest defenseNaturwissenschaften 94:247–267.https://doi.org/10.1007/s00114-006-0183-1
-
Sex-specific molecular specialization and activity rhythm-dependent gene expression in honey bee antennaeThe Journal of Experimental Biology 223:jeb217406.https://doi.org/10.1242/jeb.217406
-
Behavioural syndromes and social insects: personality at multiple levelsBiological Reviews 89:48–67.https://doi.org/10.1111/brv.12042
-
Within-nest temporal polyethism in the honey beeBehavioral Ecology and Sociobiology 62:777–784.https://doi.org/10.1007/s00265-007-0503-2
-
Division of labor in honeybees: form, function, and proximate mechanismsBehavioral Ecology and Sociobiology 64:305–316.https://doi.org/10.1007/s00265-009-0874-7
-
Neuromodulation of innate behaviors in DrosophilaAnnual Review of Neuroscience 40:327–348.https://doi.org/10.1146/annurev-neuro-072116-031558
-
Behavior programs and evolutionary strategies: natural selection sometimes favors a genetically "closed" behavior program, sometimes an "open" oneAmerican Scientist 62:650–659.
-
BookAsian Honey Bees: Biology, Conservation, and Human InteractionsCambridge, MA: Harvard University Press.
-
Social behavior in context: hormonal modulation of behavioral plasticity and social competenceIntegrative and Comparative Biology 49:423–440.https://doi.org/10.1093/icb/icp055
-
BookCaste and Ecology in the Social InsectsPrinceton NJ: Princeton University Press.
-
Genes, hormones, and circuits: an integrative approach to study the evolution of social behaviorFrontiers in Neuroendocrinology 32:320–335.https://doi.org/10.1016/j.yfrne.2010.12.004
-
The effect of genotype on response thresholds to sucrose and foraging behavior of honey bees (Apis mellifera L.)Journal of Comparative Physiology. A, Sensory, Neural, and Behavioral Physiology 182:489–500.https://doi.org/10.1007/s003590050196
-
Genetic dissection of honeybee (Apis mellifera L.) foraging behaviorJournal of Heredity 91:474–479.https://doi.org/10.1093/jhered/91.6.474
-
Genetics of reproduction and regulation of honeybee (Apis mellifera L.) social behaviorAnnual Review of Genetics 46:97–119.https://doi.org/10.1146/annurev-genet-110711-155610
-
Modulation of sucrose response thresholds in honey bees (Apis mellifera L.): influence of genotype, feeding, and foraging experienceJournal of Comparative Physiology A: Sensory, Neural, and Behavioral Physiology 187:293–301.https://doi.org/10.1007/s003590100201
-
The effect of genotype, age, sex, and caste on response thresholds to sucrose and foraging behavior of honey bees (Apis mellifera L.)Journal of Comparative Physiology A: Sensory, Neural, and Behavioral Physiology 185:207–213.https://doi.org/10.1007/s003590050379
-
Response thresholds to sucrose predict foraging division of labor in honeybeesBehavioral Ecology and Sociobiology 47:265–267.https://doi.org/10.1007/s002650050664
-
MALDI imaging analysis of neuropeptides in the africanized honeybee (Apis mellifera) brain: effect of ontogenyJournal of Proteome Research 13:3054–3064.https://doi.org/10.1021/pr500224b
-
MALDI imaging analysis of neuropeptides in africanized honeybee ( Apis mellifera) Brain: effect of aggressivenessJournal of Proteome Research 17:2358–2369.https://doi.org/10.1021/acs.jproteome.8b00098
-
Comparison of learning and memory of Apis cerana and Apis melliferaJournal of Comparative Physiology A 198:777–786.https://doi.org/10.1007/s00359-012-0747-9
-
The physiology/life-history nexusTrends in Ecology & Evolution 17:462–468.https://doi.org/10.1016/S0169-5347(02)02578-8
-
Regulation of honey bee age polyethism by juvenile hormoneBehavioral Ecology and Sociobiology 20:329–338.https://doi.org/10.1007/BF00300679
-
Regulation of division of labor in insect societiesAnnual Review of Entomology 37:637–665.https://doi.org/10.1146/annurev.en.37.010192.003225
-
Mitogenic signaling pathways induced by G protein-coupled receptorsJournal of Cellular Physiology 213:589–602.https://doi.org/10.1002/jcp.21246
-
Behavioural pharmacology of octopamine, tyramine and dopamine in honey beesBehavioural Brain Research 136:545–553.https://doi.org/10.1016/S0166-4328(02)00205-X
-
Aminergic control and modulation of honeybee behaviourCurrent Neuropharmacology 4:259–276.https://doi.org/10.2174/157015906778520791
-
Neuropeptides as regulators of behavior in insectsAnnual Review of Entomology 62:35–52.https://doi.org/10.1146/annurev-ento-031616-035500
-
A role for octopamine in honey bee division of laborBrain, Behavior and Evolution 60:350–359.https://doi.org/10.1159/000067788
-
Adaptive significance of the age polyethism schedule in honeybee coloniesBehavioral Ecology and Sociobiology 11:287–293.https://doi.org/10.1007/BF00299306
-
Behavioral syndromes: an ecological and evolutionary overviewTrends in Ecology & Evolution 19:372–378.https://doi.org/10.1016/j.tree.2004.04.009
-
Prepro-tachykinin gene expression in the brain of the honeybee Apis melliferaCell and Tissue Research 316:281–293.https://doi.org/10.1007/s00441-004-0865-y
-
Dancing to different tunes: heterospecific deciphering of the honeybee waggle danceNaturwissenschaften 95:1165–1168.https://doi.org/10.1007/s00114-008-0437-1
-
Response threshold reinforcements and division of labour in insect societiesProceedings of the Royal Society of London. Series B: Biological Sciences 265:327–332.https://doi.org/10.1098/rspb.1998.0299
-
Worker nutrition and division of labour in honeybeesAnimal Behaviour 69:427–435.https://doi.org/10.1016/j.anbehav.2004.03.017
-
"Ins and outs" of seven-transmembrane receptor signalling to ERKTrends in Endocrinology and Metabolism 16:26–33.https://doi.org/10.1016/j.tem.2004.11.008
-
Phenotypic plasticity and the origins of diversityAnnual Review of Ecology and Systematics 20:249–278.https://doi.org/10.1146/annurev.es.20.110189.001341
-
Tachykinin-related peptides modulate odor perception and locomotor activity in DrosophilaMolecular and Cellular Neurosciences 31:399–406.https://doi.org/10.1016/j.mcn.2005.10.010
Article and author information
Author details
Funding
National Natural Science Foundation of China (31970428)
- Bin Han
Agricultural Science and Technology Innovation Program (CAAS-ASTIP-2015-IAR)
- Jianke Li
National Project for Upgrading the Beekeeping Industry of China
- Jianke Li
Earmarked Fund for Modern Agro-industry Technology Research System (CARS-44)
- Jianke Li
University of North Carolina at Greensboro
- Olav Rueppell
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We appreciate Dr. Huipeng Yang at Institute of Apicultural Research for kind gifts of vectors.
Copyright
© 2021, Han et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,597
- views
-
- 237
- downloads
-
- 14
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
- Genetics and Genomics
Evolutionary arms races can arise at the contact surfaces between host and viral proteins, producing dynamic spaces in which genetic variants are continually pursued. However, the sampling of genetic variation must be balanced with the need to maintain protein function. A striking case is given by protein kinase R (PKR), a member of the mammalian innate immune system. PKR detects viral replication within the host cell and halts protein synthesis to prevent viral replication by phosphorylating eIF2α, a component of the translation initiation machinery. PKR is targeted by many viral antagonists, including poxvirus pseudosubstrate antagonists that mimic the natural substrate, eIF2α, and inhibit PKR activity. Remarkably, PKR has several rapidly evolving residues at this interface, suggesting it is engaging in an evolutionary arms race, despite the surface’s critical role in phosphorylating eIF2α. To systematically explore the evolutionary opportunities available at this dynamic interface, we generated and characterized a library of 426 SNP-accessible nonsynonymous variants of human PKR for their ability to escape inhibition by the model pseudosubstrate inhibitor K3, encoded by the vaccinia virus gene K3L. We identified key sites in the PKR kinase domain that harbor K3-resistant variants, as well as critical sites where variation leads to loss of function. We find K3-resistant variants are readily available throughout the interface and are enriched at sites under positive selection. Moreover, variants beneficial against K3 were also beneficial against an enhanced variant of K3, indicating resilience to viral adaptation. Overall, we find that the eIF2α-binding surface of PKR is highly malleable, potentiating its evolutionary ability to combat viral inhibition.
-
- Ecology
- Evolutionary Biology
Seasonal polyphenism enables organisms to adapt to environmental challenges by increasing phenotypic diversity. Cacopsylla chinensis exhibits remarkable seasonal polyphenism, specifically in the form of summer-form and winter-form, which have distinct morphological phenotypes. Previous research has shown that low temperature and the temperature receptor CcTRPM regulate the transition from summer-form to winter-form in C. chinensis by impacting cuticle content and thickness. However, the underling neuroendocrine regulatory mechanism remains largely unknown. Bursicon, also known as the tanning hormone, is responsible for the hardening and darkening of the insect cuticle. In this study, we report for the first time on the novel function of Bursicon and its receptor in the transition from summer-form to winter-form in C. chinensis. Firstly, we identified CcBurs-α and CcBurs-β as two typical subunits of Bursicon in C. chinensis, which were regulated by low temperature (10 °C) and CcTRPM. Subsequently, CcBurs-α and CcBurs-β formed a heterodimer that mediated the transition from summer-form to winter-form by influencing the cuticle chitin contents and cuticle thickness. Furthermore, we demonstrated that CcBurs-R acts as the Bursicon receptor and plays a critical role in the up-stream signaling of the chitin biosynthesis pathway, regulating the transition from summer-form to winter-form. Finally, we discovered that miR-6012 directly targets CcBurs-R, contributing to the regulation of Bursicon signaling in the seasonal polyphenism of C. chinensis. In summary, these findings reveal the novel function of the neuroendocrine regulatory mechanism underlying seasonal polyphenism and provide critical insights into the insect Bursicon and its receptor.