Skip to main content
  • Research article
  • Open access
  • Published:

Mating precedes selective immune priming which is maintained throughout bumblebee queen diapause



Understanding the mechanisms by which organisms adapt to unfavourable conditions is a fundamental question in ecology and evolutionary biology. One such mechanism is diapause, a period of dormancy typically found in nematodes, fish, crustaceans and insects. This state is a key life-history event characterised by arrested development, suppressed metabolism and increased stress tolerance and allows an organism to avoid prolonged periods of harsh and inhospitable environmental conditions. For some species, diapause is preceded by mating which can have a profound effect on female behaviour, physiology and key biological processes, including immunity. However, our understanding of how mating impacts long-term immunity and whether these effects persist throughout diapause is currently limited. To address this, we explored molecular changes in the haemolymph of the ecologically important pollinator, the buff-tailed bumblebee Bombus terrestris. B. terrestris queens mate prior to entering diapause, a non-feeding period of arrested development that can last 6–9 months. Using mass-spectrometry-based proteomics, we quantified changes in the pre-diapause queen haemolymph after mating, as well as the subsequent protein expression of mated queens during and post-diapause.


Our analysis identified distinct proteome profiles associated with diapause preparation, maintenance and termination. More specifically, mating pre-diapause was followed by an increase in the abundance of antimicrobial peptides, key effectors of the immune system. Furthermore, we identified the elevated abundance of these proteins to be maintained throughout diapause. This finding was in contrast to the general reduction observed in immune proteins during diapause suggestive of selective immune priming and expression during diapause. Diapause also affected the expression of proteins involved in cuticular maintenance, olfaction, as well as proteins of unknown function, which may have roles in diapause regulation.


Our results provide clear molecular evidence for the consequences and benefits of mating at the immune level as it precedes the selective increased abundance of antimicrobial peptides that are sustained throughout diapause. In addition, our results provide novel insights into the molecular mechanisms by which bumblebees prepare for, survive, and recover from diapause, insights that may have implications for our general understanding of these processes in other insect groups.


The preparation for, as well as maintenance of life within harsh climates, places enormous selective pressures on organisms, resulting in the evolution of adaptations to survive environmental extremes [1, 2]. One such adaptation is diapause, a dormancy-like state whereby organisms reduce or suppress metabolism and can arrest other important biological processes, such as development or reproduction [3,4,5]. The capacity for diapause is widespread across evolutionarily diverse taxa and has been suggested to contribute to their ability to exploit seasonal resources, as well as bridge periods of unfavourable environments (e.g. winter) [6]. For insects, within any given species diapause is restricted to one of several life stages, including different larval instars and pupal stages, as well as the adult stage of the life-cycle [6, 7]. Most species of diapausing insects can enter diapause in the absence of environmental cues but, more commonly, changes in environmental variables, such as photoperiod or temperature, lead to the preparation and programming of diapause [4, 8,9,10].

Regardless of the developmental stage of diapause entry, the dormancy-like state consists of three ecophysiological phases including pre-diapause (preparation and induction), diapause (initiation, maintenance and termination), and post-diapause, where development is reinitiated [8]. The expression of the diapause phenotype is regulated by ecdysteroids, maternal effects, and altered gene expression [3, 4] resulting in arrested development, suppressed metabolism, increased longevity and enhanced stress tolerance [3, 5, 11]. However, despite our increasing knowledge of certain physiological processes during diapause, such as reproduction and metabolism, our understanding of other key processes, such as immunity, remains limited [12, 13]. The reduced mobility and metabolic activity of diapausing individuals can place them at increased risk of pathogen infection, which may affect their survival rates and/or affect post-diapause fitness. This is coupled with constitutive immune expression being costly to generate and maintain [14,15,16,17], especially during periods of low nutritional intake [18, 19]. Diapausing organisms are also at greater risk of desiccation and require mechanisms to reduce water loss [20,21,22].

An additional consideration is the effect of mating on the preparation and maintenance of diapause, which in general can have profound long-term effects on female physiology affecting behaviour, immunity and longevity. In relation to diapause, both sexes of some insects undergo diapause with copulation occurring after termination of diapause, in others only females undergo diapause with mating occurring before diapause [23]. During mating, males transfer sperm into the female reproductive tract, where it is stored in a specialised organ for use post-diapause [24,25,26]. Males may also transfer additional compounds during copulation that modulate aspects of female physiology such as receptivity to remate [27, 28] or immune expression [29,30,31,32,33]. While mating can suppress aspects of the female immune response, in species where mating occurs before diapause this may incur a survival cost, and thus sexual selection might be expected to act on males to ensure that mating does not decrease, or perhaps even increases the likelihood of the female surviving diapause, to ensure paternity. However, at present, our understanding of the long-term effects of mating during diapause is limited.

An exemplary model to explore the prolonged effect of mating on immune expression during diapause is the bumblebee. Bumblebees are key ecological and commercial pollinators providing ecosystem services that maintain biodiversity [34] and improve agricultural crop yields [35]. Like other social Hymenoptera, caste differentiation is present [36], whereby the majority of females, known as workers, perform altruistic tasks, such as resource foraging, brood care, nest maintenance and defence [37], while generally a single reproductive individual, the queen, establishes the colony before transitioning primarily to egg-laying [38]. Despite morphological similarities between queens and workers in certain species, queens often display different behavioural and physiological repertoires in comparison to workers [36] and are the only colony members to undergo diapause [38]. Before diapause, virgin queens, known as gynes, emerge from the natal colony to find a mate. For most bumblebee species, queens mate with just one male [39, 40] before entering diapause, a non-feeding period that can last 6–9 months [41, 42]. Mating in the buff-tailed bumblebee, Bombus terrestris, results in an increased transcript expression of antimicrobial peptides, which coincides with increased resistance of mated queens to the pathogenic trypanosome, Crithidia bombi [43] demonstrating the immune benefit of mating. However, it is unknown whether this heightened immunity is maintained throughout bumblebee diapause and indeed, after, where queens are also susceptible to infection [44]. Previous studies have identified the role of the endocrine juvenile hormone (JH) [45] and genome-wide changes in transcriptional regulation associated with key metabolic processes, as well as nutrient storage and stress response [46], in the maintenance of bumblebee diapause. While transcriptomics provides valuable insights into gene expression, mRNA or transcript quantification is often a weak to modest proxy for the estimation of the abundance of proteins, which at the biochemical level perform functions directly related to the survival and fitness of an organism, especially during periods of stress [47].

Techniques, such as mass spectrometry-based proteomics, simultaneously identify and quantify thousands of proteins in a biological sample, which enables the assessment of environmental and/or physiological changes at the molecular phenotype level. In addition, RNA approaches are not conducive to the direct assessment of the protein constituents of biological fluids such as haemolymph. In this study, we characterised the haemolymph, a nutrient-rich plasma-like medium present in the haemocoel (open circulatory tract) of arthropods, which is responsible for the transport and delivery of nutrients to tissues [48, 49]. The haemolymph also contains proteins with roles in immunity [50,51,52,53], development [54,55,56], olfaction [48, 57] and neurotransmission [58, 59]. Profiling changes in haemolymph protein abundance has the potential to resolve global molecular signatures and biomarkers underlying key life stage events [59,60,61], including diapause [62]. Here, using mass spectrometry-based proteomics, we aim to establish a clear understanding of the molecular mechanisms underlying the ecophysiological stages of onset, maintenance and termination of diapause, by profiling the haemolymph proteome throughout diapause for the buff-tailed bumblebee, B. terrestris.

In our experiment, we specifically address three key questions. Firstly, how does the queen haemolymph proteome change as a result of mating, and are these changes maintained throughout and after diapause? Secondly, given that diapause represents a period of presumed metabolic inactivity, can changes in protein abundance, particularly of those associated with immune/defence systems and metabolic processes, provide insight into how the diapausing organism prepares for, survives and recovers from diapause. Finally, can the regulatory components of diapause be discerned and identified from queen haemolymph and do they concur with those identified in recent studies using transcriptomic analyses [46, 63] highlighting conserved genetic regulatory mechanisms across insect taxa. To address these questions, we quantified changes in the haemolymph proteome composition at six stages of the queen life-cycle, including two pre-diapause stages (virgin gynes and mated pre-diapause queens), two diapause stages (1 week and 14 weeks during diapause) and two post-diapause time-points (6 h and 48 h post-diapause). Our study reveals novel changes pre-, during and post-diapause in the queen haemolymph proteome highlighting the effect of mating on the queen molecular phenotype, as well as providing key insights into the molecular mechanisms underlying the regulation and expression of the diapause phenotype itself.


Haemolymph profiling reveals altered diapause phenotype

In total 129 unique multiple-peptide supported haemolymph proteins were identified from the virgin and mated queens collected across the six time-points examined (Additional file 1: Table S1) of which 79 were identified to be statistically significant differentially abundant (SSDA) following multivariate analysis (ANOVA, FDR < 0.05) (Additional file 2: Table S2a).

To determine whether samples clustered based on life-cycle stage, we performed a principal component analysis (PCA) on Z-score normalised label-free quantification (LFQ) intensities for the 79 SSDA proteins. The first two principal components explained 29.6 and 18.19% of the variance within the dataset, respectively. The first principal component clearly separated early, late diapause and 6 h post-diapause individuals from the three other time-points (Fig. 1) highlighting distinct differences in the haemolymph proteome profiles associated with diapause. As a complementary approach, we performed hierarchical clustering of Z-score normalised LFQ intensities for the 79 SSDA proteins and identified ten distinct clusters of proteins with similar expression profiles (Fig. 2; Additional file 3: Table S3). Two clusters (Clusters B and C) comprise proteins that represent a clear diapause phenotype with proteins of high or highest abundance in the haemolymph of early diapause, late diapause or 6 h post-diapause queens. An additional cluster (Cluster A; n = 6) consisted of proteins that had the highest abundance during late diapause. There were also two clusters (Clusters G and H) of proteins with reduced abundances during diapause with the largest identified cluster (Cluster G; 27/79 proteins) consisting of proteins with low abundances during late diapause and six hours post-diapause. By 48 h post-diapause, the majority of these proteins had increased at least two-fold in abundance with levels comparable to pre-diapause haemolymph. Five proteins were identified with elevated abundances post-mating (Cluster D). Cluster E (n = 6 proteins) and cluster F (n = 4 proteins) comprised proteins with the highest abundances in early diapausing queens and 48 h post-diapause, respectively. An additional cluster (Cluster I; n = 4 proteins) consisted of proteins that had relatively high abundances in gynes, pre- and early-diapausing queens but decreased abundance in the haemolymph of late diapausing queens and at 48 h post-diapause. The last cluster (Cluster J; n = 3 proteins) consisted of proteins elevated post-mating and maintained at high levels during diapause before a reduction post-diapause.

Fig. 1
figure 1

Label-free quantitative proteomic analysis of bumblebee queen haemolymph at six points within the life-cycle. Principal component analysis (PCA) of Z-score normalised label-free quantification (LFQ) intensity values for all statistically significant differentially abundant (SSDA) proteins (n = 79) identified (a) the first two components (PC1 and PC2) to explain 29.6 and 18.19% of the variance within the dataset, respectively. To further resolve differences in the haemolymph proteome, additional principal components (b) PC3, (c) PC4, and (d) PC5 were plotted individually against the first principal component. Each queen life stage is indicated by a unique colour

Fig. 2
figure 2

Distinct haemolymph proteome profile revealed during diapause. Two-way unsupervised hierarchical clustering of the median Z-score normalised label-free quantification (LFQ) intensity values of all statistically significant differentially abundant (SSDA) proteins (n = 79) for each life stage. Hierarchical clustering resolved 10 distinct clusters across the six queen life stages based on similar expression profiles of individual proteins. Differences in protein abundance are indicated by colour changes from low (green) to high (red) protein abundance representative of changes in Z-score normalised log2-fold transformed LFQ intensity values

Comparative analyses of consecutive stages of the queen life cycle

Pairwise t-tests (p < 0.05) were performed to identify differences in mean LFQ intensities between consecutive time-points which enabled the characterisation of stage-specific changes in the queen life-cycle, including changes in response to mating, diapause entry, as well as diapause termination (Additional file 2: Table S2b).

Comparison of virgin and mated pre-diapause queens: increased abundance of AMPs in mated queen haemolymph

We identified ten proteins with significant (two-sample t-test; p < 0.05) abundance changes between the haemolymph proteome of virgin and mated queens (Additional file 2: Table S2b). Three of these proteins were annotated as putative antimicrobial peptides (AMPs) with the greatest increase evident for hymenoptaecin (ADB29130.1), which was 183-fold more abundant (p < 10− 4; Fig. 3) in the haemolymph of mated queens. Increased abundances were also observed for defensin (17.9 fold increase; p < 10− 4) and abaecin (9.4 fold increase; p < 10− 3; Fig. 3). Post-mated queens also had increased abundance of a putative esterase FE4 (XP_003397300.2; 7.24 fold increase; p < 10− 2) and a putative kappa-theraphotoxin-like protein (XP_020718646.1; 3.48 fold increase; p < 0.05). Five proteins had reduced abundances within the mated queen haemolymph including a BMP endothelial regulator protein (XP_003393207.1), an ecdysteroid regulated 16 kDa (XP_003395667.1), a ferritin subunit (XP_012167824.1), a putative muscle protein (XP_020722973.1) and an odorant binding protein (XP_003397877.1) demonstrating that mating can affect the expression of proteins involved in diverse biological processes.

Fig. 3
figure 3

Elevated abundances of antimicrobial peptides in the queen haemolymph post-mating. Box plot representations of label-free quantification (LFQ) intensity values at six points throughout the bumblebee queen life-cycle for: (a) abaecin; (b) defensin; (c) hymenoptaecin; and (d) apidaecin. For each protein, the y-axis shows the LFQ intensity values representative of protein abundance. The x-axis represents the six sampled queen life stages. Each box is coloured to indicate a different life stage of the queen. Stages with a significant difference in abundance (two-sample t-test; p < 0.05) are indicated with asterisks

Comparison of mated pre-diapause and early diapause queens: early diapause phenotype characterised by an increase in cuticular and ejaculatory bulb-specific proteins

Sixteen proteins changed significantly in abundance between the haemolymph of mated pre-diapause and early diapause queens. There was a significant trend of elevated expression of proteins (n = 13) in the haemolymph of early diapause queens (binomial test, p = 0.02). Of the 13 proteins with elevated abundance, three were annotated as ejaculatory-bulb specific proteins (each had at least > 2 fold increase; p < 0.05; Fig. 4)a, while four were annotated as cuticular proteins, including one protein which was increased 478-fold within the haemolymph of early diapause queens in comparison to pre-diapause queens (Fig. 4)b. Additional proteins with significantly increased abundances included a serine protease inhibitor (XP_003401209.1), a neurofilament heavy polypeptide (XP_003394564.1), a heat shock protein (XP_003402976.1), a gamma interferon-inducible lysosomal thiol (GILT) reductase (XP_003397075.2), a putative scavenger receptor protein (XP_012163150.1) and a BMP-binding endothelial regulator (XP_003393256.1). We identified three proteins with significant reductions in the early diapause queen haemolymph, including a putative salivary protein (XP_012167732.1), a prion-like protein (XP_012175279.1) and a hexamerin (XP_003401781.1), which was reduced 23.9 fold.

Fig. 4
figure 4

Olfaction and structural-related protein changes in queen haemolymph in response to diapause. a Box plot representations of protein abundances at six points throughout the bumblebee queen life-cycle for six olfaction-related proteins including three odorant binding proteins and three ejaculatory bulb-specific proteins. For each protein, the y-axis shows the LFQ intensity values. The x-axis represents the queen life stages. Each box is coloured to indicate a different life stage of the queen. Stages with a significant difference in abundance (p < 0.05) are indicated with astrisks. b Line plot shows relative abundance changes between structural-related proteins within the queen haemolymph at six time-points. The y-axis shows the LFQ intensity values. The x-axis represents the individual queen life stages

Comparison of early and late diapause queens: selective reduction in immune expression in late diapause queen haemolymph

We identified 29 SSDA proteins between the haemolymph of early and late diapause queens (Additional file 2: Table S2b). The majority of these proteins (n = 21) had lower abundance within late diapause queens, which was a significant trend (binomial test, p < 0.03). Of the 21 proteins whose abundance was reduced at least 1.5 fold in the haemolymph of late diapause queens, 12 had putative roles in immunity, including two pathogen recognition proteins, two immune signalling proteins, five effector proteins, two regulatory proteins and one putative detoxification enzyme.

Comparison of late diapause and 6 h post-diapause queens: reduced abundance of detoxification and structural proteins in queen haemolymph post-diapause

We identified 14 SSDA proteins within the queen haemolymph between six hours post-diapause and late diapause (Additional file 2: Table S2b). Of the nine proteins significantly reduced at least 2 fold within the haemolymph of six hours post-diapause queen haemolymph, three proteins were annotated with putative roles in detoxification (XP_003397315.1: superoxide dismutase; XP_020719829.1: glutathione peroxidase; XP_003394648.1: cytochrome c). Four of the reduced proteins had putative structural roles with two annotated as muscle-associated proteins (XP_012172948.1; XP_020722973.1), one annotated as an actin-depolymerising protein (XP_012174797.1), and one protein annotated as a cuticular protein (XP_003394953.1). The final two proteins with significantly reduced expression were annotated as an odorant binding protein (XP_003397877.1) and pterin-4-alpha-carbinolamine dehydratase (XP_012170199.1). Five proteins had at least a 1.3 fold increase in abundance, including an endoribonuclease (XP_012167789.1), a glyoxylate reductase/hydroxypyruvate reductase (XP_012170485.1), an iron-binding ferritin (XP_012167818.1), a putative clotting factor (XP_012175299.1), as well as a hypothetical protein of unknown function (XP_003395010.1).

Comparison of 6 and 48 h post-diapause queens: recovery in queen reproductive and immune protein expression at 48 h post-diapause

The greatest number of SSDA proteins (n = 52) was identified between the haemolymph proteome of queens at 6 h and 48 h post-diapause with the majority (n = 32) increasing in abundance in the haemolymph 48 h post-diapause although this pattern was not significant (binomial test, p = 0.1263). The greatest increase in abundance was seen in a storage protein (XP_003401781.1: hexamerin), which had a 49.8 fold increase. Two other putative reproductive proteins (XP_020718283.1: membrane metalloendopeptidase-like 1; XP_012163499.1: vitellogenin) increased in abundance at least 5-fold in the haemolymph at 48 h post-diapause. The haemolymph proteome also had increased abundances of proteins with putative roles in immunity (n = 12), muscle-associated proteins (n = 2), olfaction (n = 2), venom (n = 3), as well as proteins of unknown function (n = 5). The remaining proteins were annotated with putative roles in cholesterol metabolism (XP_003395667.1: ecdysteroid-regulated 16 kDa), chitin binding (XP_003393340.1: peritrophin-1), neurotransmitter synthesis (XP_003401022.1: glutamate decarboxylase), as well as gene expression regulation (XP_020722259.1: trithorax).

Twenty proteins had significantly reduced expression in the queen haemolymph at 48 h post-diapause. The protein with the greatest reduction was a single domain Von Willebrand factor type C domain-containing protein (XP_003399812.1), which had a 22-fold reduction. Certain immune-associated proteins, including two AMPs (ADB29130.1: hymenoptaecin; ADB29128.1: abaecin), two putative serine protease inhibitors (XP_003398424.1: serpin-B3; XP_012169463.2: alaserpin), and a putative anti-viral protein (XP_003399869.1: protein son of sevenless) were reduced in comparison to the haemolymph of queens at six hours post-diapause. A sixth putative immune protein, a leucine-rich repeat domain-containing protein (XP_003394691.1) also had reduced abundance. Five proteins with putative structural roles were reduced, including three cuticular proteins, a neurofilament heavy polypeptide (XP_003394564.1) and a chitinase (XP_003394129.1). Similar to other time-points, putative olfaction-related proteins were also affected with two proteins (XP_003397885.1: OBP3; XP_012173545.1: PEBIII) reduced within the haemolymph at 48 h post-diapause. Other proteins reduced at least 2.6 fold included a glyoxylate reductase/hydroxypyruvate reductase (XP_012170485.1), a GILT protein (XP_003397075.2), an esterase (XP_003397300.2), a glutathione peroxidase (XP_003395541.1), a protein of unknown function (XP_003395010.1) and a putative BMP-binding endothelial regulator (XP_003393256.1). Interestingly, this last protein is coded for by a gene that forms part of a novel gene family consisting of five genes coding for proteins expressed at high abundance in the queen haemolymph (Additional file 4; Additional file 2: Table S2b).

Proteogenomic analysis of a novel haemolymph-associated protein family

It is estimated that more than 40% of genes in sequenced eukaryotic genomes do not have assigned functions [64]. In our analysis, we provide evidence for the expression of 33 hypothetical or uncharacterised proteins in the queen haemolymph, of which 16 changed in abundance over the life stages of the queen. Of most interest, we identified five novel proteins which shared high sequence similarity (minimum sequence similarity of 45.75%; E-value <1e-58) that were all expressed in the haemolymph proteome. These five proteins were coded for by four individual protein-coding genes suggestive of a previously undocumented gene family (Additional file 5: Table S4).

These four genes, hereafter known as the highly abundant haemolymph-associated protein (HAHP) family are located within a genomic region spanning approximately 100 kb on chromosome one (BG_1; NC_015762.1) in the B. terrestris genome (Additional file 4: Fig. S1). There is evidence of a fifth protein-coding gene (LOC110119444), but this predicted protein was not identified in the queen haemolymph in our analysis. Although one of the proteins was annotated as a BMP-binding endothelial regulator protein (XP_003393256.1), additional functional domain analysis did not identify any conserved domains present on any of the HAHP family members. At the nucleotide and protein level, we identified high sequence similarity with five putative orthologues identified in the genome of the Eastern bumblebee, B. impatiens. We performed homology searches against other insect predicted proteomes identifying single protein matches in each of the four sequenced honeybee (Apis) genomes. Similarly, functional domain analysis of honeybee species homologues only identified the presence of a predicted signal peptide domain. No homologous sequence was identified for the fruit fly Drosophila melanogaster.

Interestingly, one member of the HAHP family (XP_003393256.1; LOC100652150: HAHP1) displayed an expression profile that suggests a potential role in diapause (Additional file 4: Fig. S2). This protein has relatively high abundance levels in gynes and pre-diapause queens and significantly increases in early, mid diapause and 6 h post-diapause queens. Protein abundance returns to pre-diapause levels at 48 h post-diapause.

Weak conservation of protein and transcript expression profiles during diapause

In an attempt to explore the conservation of molecular mechanisms underlying diapause, we compared our proteomic dataset with a previously published transcriptomic dataset of bumblebee queens collected before, during and after diapause [46]. We identified 114 genes (88.4% of genes coding for haemolymph-associated proteins) expressed in the queen fat body at three time-points [46](Additional file 6: Table S5). Of this number, only three genes (LOC100643414: protein spaetzle; LOC100648549: cytochrome c; LOC100651094: glycine-rich cell wall structural protein) had conserved directional changes in transcript and protein abundance during diapause (Additional file 4). The majority differed in direction of expression profile at the transcript and protein level, which was a significant trend (binomial test, p < 0.005) and highlights a weak association between transcript and protein expression. We identified two Gene Ontology terms (GO:0032504, ‘multicellular organism reproduction’; GO:0005615, ‘extracellular space’) as significantly enriched in both the transcriptome and haemolymph proteome (Additional file 7: Table S6).


Diapause is a key component of insect life-cycles. Using bumblebee queens as a model system to investigate physiological adaptations to diapause, we show that: 1) mating is followed by significant increases in antimicrobial peptides, which are maintained at highly elevated levels throughout diapause; 2) in contrast to elevated AMP expression, other immune proteins reduce in abundance during diapause; 3) diapause is also characterised by increases in cuticular and muscular proteins, and 4) due to weak conservation in transcript and protein expression, transcriptomic approaches alone may not fully represent the biology of diapause phenotypes. Taken together these results provide novel insights into the molecular mechanisms involved in diapause preparation, maintenance and termination within an ecologically and economically important pollinator.

Dormancy strategies, such as diapause, have evolved to allow organisms avoid periods of harsh environmental conditions when resources may be low [3, 65]. A characteristic of diapause is an arrest in development and reproduction. These processes are metabolically and energetically demanding, which can result in physiological and evolutionary trade-offs with other biological processes [66]. While oogenesis places energetic demands on the adult female [67], the act of mating can result in behavioural and physiological changes in the female phenotype [66]. While the arrest of reproductive potential may free up resources for other costly physiological processes, such as immunity, this has not been previously investigated during diapause. We identified dynamic changes in protein profile and abundance during diapause (summarised as Fig. 5), especially in relation to immunity-associated proteins. Firstly, mating was followed by significant increases in the abundance of antimicrobial peptides, including a defensin, a hymenoptaecin and an abaecin in the haemolymph. Antimicrobial peptides are generally induced upon infection and display antibacterial and antifungal activities including inhibition of cell growth and lytic activity of foreign microorganisms [68,69,70,71]. While mating can result in expression of AMPs in the female reproductive tract possibly to restrict the movement of sexually transmitted diseases, copulation in insects generally results in reduction in female immune potential [66]. This immunosuppression is suggested to be due to resources being used to increase fecundity, which is regulated by products of the endocrine system, which promote reproductive development and suppress aspects of the immune system, including expression of AMPs. However, in females that undergo diapause, organs that synthesise and secrete hormones that regulate reproductive development, such as juvenile hormone (JH), are not fully developed [72]. The lack of development of the corpora allata may be the reason for both increased expression at the protein and transcript level in bumblebee queens. Furthermore, increased AMP expression in mated queens is associated with increased resistance to the trypanosome, Crithidia bombi, demonstrating the immune benefits of mating [43]. Here, not only do we provide evidence of elevated protein abundances post-mating but also show that these levels are sustained during diapause for a period up to 120 days (approximately one-third of a queen’s expected lifespan). Given how energetically costly activation and constitutive immune expression is for an organism [15, 73], the maintenance of elevated protein levels in the haemolymph would be suggested to directly increase diapause survival rather than as an indirect consequence of mating. Indeed, infections can occur during and immediately after diapause [44, 74, 75], which can affect survival and therefore, immune defences are required. However, selective expression of immune components may be sufficient to reduce infections.

Fig. 5
figure 5

Summary overview of protein expression profiles in relation to life stage and colony life cycle. Broad trends in the abundance of protein categories in unmated and mated queens prior to (green), during (blue) and after (red) diapause. Proteins associated with the immune system, cuticle, storage and reproduction, regulation and stress displayed characteristic expression profiles across diapause suggestive in roles in the preparation for, maintenance of and recovery from diapause. L and H represent low and high relative protein abundances, respectively

In contrast to AMP expression, we identified a reduction in proteins with putative immunological roles. Periods of reduced nutritional intake, such as hibernation and diapause, can reduce resources available for certain immune expression and activity. Even in non-immune challenged individuals, metabolic trade-offs are evident between different facets of the immune response, most notably between the humoral and cellular response [76, 77], with different nutritional demands suggested for the expression of each arm [78]. The reduction in immune proteins observed in our study may also be self-regulatory, as suggested by the increase in serine protease inhibitors, proteins which reduce and inactivate serine proteases, during queen diapause. Increased expression of AMPs, such as drosomycin and cecropin, has been seen in diapausing D. melanogaster adults in response to bacterial infection [13]. In the same study, a peptidoglycan recognition protein (PGRP) was not elevated in diapause samples in comparison to post-diapause samples providing additional evidence that aspects of the host immune system display differences in expression during diapause. Overall, we suggest that sentinel and responsive components of the immune system are depressed during diapause when the metabolic costs of driving these systems may not be sustainable. This view is supported by the fact that the end products of these recognition and signalling pathways are AMPs, which are highly elevated throughout diapause. As such, an anticipatory or prophylactic defence system may have evolved in bumblebee queens to protect the diapausing queen from microbial attack without the need to first detect and respond to pathogen presence.

While infections can reduce rates of diapause survival, diapause also represents a period where organisms are at risk of desiccation and adaptations are required to reduce water loss. One important barrier to water loss is the cuticle [79]. During the early diapause stage, cuticular proteins had elevated abundances in the queen haemolymph. Previously, a cuticular protein identified in the head tissues of early diapause Culex adult females [80] was suggested to function in cuticle reinforcement during this period. Increased protein expression in the haemolymph of both Bombus and Culex suggests a putative conserved role of cuticular proteins in diapause. While the exact function of these proteins are uncharacterised, the expression of cuticular proteins, which function in cuticular development during metamorphosis, is regulated by ecdysteroids, such as JH [81]. Consequently, low levels of JH [45, 46] evident during adult diapause for queens and other insects may facilitate the elevated expression of cuticular proteins providing an important defence against desiccation.

The diapause phenotype is tightly regulated by hormones affecting numerous biological processes of the diapausing individuals [82]. Our study has determined that such processes may involve the expression of olfaction-related proteins, most notably odorant binding proteins (OBPs) and ejaculatory bulb-specific proteins (PEBs), which display complex expression changes in queen haemolymph. OBPs and PEBs represent families of small hydrophilic molecules, which have been associated with the transportation of odour molecules in insect haemolymph [83, 84]. While the association of these proteins with olfactory organs, such as antennae [85], has suggested a role in olfaction, additional roles in other processes have been identified, including reproduction and immunity. For example, ejaculatory bulb-specific proteins, which contain chemosensory protein domains, have been characterised as major components of ejaculatory bulbs in the fruit fly D. melanogaster [86, 87]. While the exact role of these proteins in B. terrestris queens will require experimental validation, these types of proteins have also been shown to increase in abundance in flies in response to viral infection, although abundance did not coincide with a reduction in pathology suggesting indirect roles in viral response, such as tissue repair or behavioural response to infection [88]. Interestingly, this study also identified suppression of PEBIII expression in S2 cells exposed to the moulting steroid hormone, 20-hydroxyecdysone (20E), an ecdysteroid also involved in eliciting diapause termination within a number of insect species [88]. Similarly, the expression of OBPs in the developing antennae of the tobacco hawkmoth Manduca sexta is regulated by declining levels of the ecdysteroid 20E [89]. We identified the abundance of three OBPs to be affected during diapause although there was no general pattern of their response to diapause with only an OBP19 having elevated expression during early and late diapause suggestive of a potential role in diapause maintenance. While 20E ecdysteroid levels were not examined in the present study, low ecdysteroid levels have been found in bumblebee queens during diapause [46]. Thus, changes in their abundance in the haemolymph may be a consequence of changes in ecdysteroid levels associated with insect diapause. Further work will be required to elucidate the specific roles performed by olfaction-associated haemolymph proteins during insect diapause.

The convergence of diapause phenotypes across phylogenetically diverse taxa has led to the suggestion that a conserved genetic toolkit may regulate diapause [46, 63, 90]. To explore conserved changes at the transcript and protein level in queen diapause, we compared changes in the haemolymph proteome to transcriptional changes identified by a previous study in the fat body of bumblebee queens [46]. The insect fat body has an important role in storage and energy usage, as well as playing a central role in biosynthesis and metabolic activity [91], and also synthesises most of the haemolymph proteins [92]. Due to this role in haemolymph protein synthesis, it was unsurprising to find that the majority of haemolymph-associated proteins were transcribed in the fat body across three time-points analysed by Amsalem et al. [46]. However, we identified only weak overlap in expression profiles between the fat body transcriptome and haemolymph proteome. While transcriptomes provide important insights into genome-wide transcriptional regulation, here we show that these signals may not translate at the protein level, which is an important finding when considering candidate genes involved in the expression of the diapause phenotype. While these differences may be explained by biological and technical differences between the studies [93], a similarly modest correlation between mRNA transcript and protein levels has been described within model organisms, such as mice [94] and humans [95]. The limited overlap between the transcriptome and proteome highlights that in certain cases gene expression may not be an appropriate proxy for protein products expressed during diapause.

Given that the haemolymph represents a non-cellular fraction, proteomics has an added advantage of providing direct information on spatial and temporal protein expression and potential function not amenable by RNA-directed methods. For example, our results supplement recent genomic resources for bumblebee species [96,97,98,99], through the provision of translational and spatial evidence for 129 multiple-peptide supported haemolymph-associated proteins, including seven uncharacterised proteins devoid of known functional motifs and several members of a putative expanded gene family in queen haemolymph. Previous work by Sadd et al. [92] described 17 haemolymph-associated proteins of unknown function in post-diapause queen haemolymph (HAP1–14 and three putative members of the HAHP family). Three of these HAHP members, including an additional fourth, expressed in the haemolymph throughout the queen life cycle were identified in this study. In addition, the abundance of one member (HAHP1) changed during diapause suggesting a potential role in diapause maintenance. The HAHPs share sequence similarity with homologous proteins annotated as BMP endothelial regulator proteins and kielin/chordin-like proteins, which is suggestive of a role in development. However, their exact role in bumblebees will require further genomic, functional and expression studies to elucidate the biological processes associated with this diverse, haemolymph profuse and potentially important group of proteins.


Here we provide a novel insight into the molecular phenotype underlying insect diapause, identifying distinct haemolymph profiles associated with each of the key ecophysiological stages of diapause. Our data indicate that mating precedes a specific increase in antimicrobial peptides, which persists throughout queen diapause suggesting a key role for these proteins in the protection and survival of queens. We provide new perspectives on the molecular mechanisms underlying diapause, highlighting the increase in cuticular proteins to reduce desiccation, the reduction in non-AMP immune components, as well as highlighting dynamic changes in aspects of olfactory proteins throughout diapause. We also demonstrate a lag period in recovery post-diapause. Finally, our results highlight the importance of proteomics in understanding the molecular mechanisms underlying diapause regulation and maintenance. Future research in this area will benefit through increased tissue-specific profiling, as well as functional characterisation of haemolymph-associated proteins, to provide a greater understanding into the molecular mechanisms underlying diapause in bumblebees.


Animal handling and haemolymph collections

Four week old queen-right Bombus terrestris audax colonies were obtained from a commercial supplier (Koppert, the Netherlands) and maintained in an environmentally controlled room (27 ± 1 °C, 45% relative humidity (RH)) under red light illumination. Colonies were provided with pollen and sugar water (ApiInvert) ad libitum and maintained for approximately six weeks until new sexuals were produced.

For the purpose of characterising changes in the haemolymph proteome over life cycle stages related to diapause, six time-points were examined: 1) virgin queens (gynes); 2) pre-diapause mated queens; 3) early diapause queens; 4) late diapause queens; and post-diapause queens collected at 5) 6 h post-diapause; and 6) 48 h post-diapause. For each time-point, we collected haemolymph from four individual queens. In brief, each natal colony was checked twice daily for the presence of newly emerged gynes. Callow gynes detected within a 24 h period were transferred from the natal colony to a separate wooden box (22 × 15 × 10 cm) where they were supplied with sugar water (50%) and pollen ad libitum. Gynes were maintained in groups in these nurseries until 5 days post-eclosion, at which point we collected haemolymph from four individual virgin queens. On the same day, remaining sexually mature gynes were mated. Matings were performed under standard laboratory conditions [100] and queens were mated with males from different colonies to avoid full-sibling inbreeding [101]. Three days post-mating, we collected haemolymph from four newly mated queens. The remaining queens were transferred to individual 50 ml centrifuge tubes, which contained 5 ml autoclaved sand. Queens were entered into diapause via incubation of animals at 4 °C (RH of 70%) in the dark. Diapausing queens (n = 4) were sampled at seven days into diapause to examine proteomic changes associated with early diapause. We also collected a second time-point (n = 4) at 14 weeks into diapause (late diapause). Post-diapause, queens were transferred to a temperature-controlled room at 27 ± 1 °C (RH of 45%). Each queen was maintained in an individual plastic box (dimensions: (L) 11.5 cm; (W) 9 cm; (H) 7 cm) lined with autoclaved clean sand. We collected haemolymph from individual bees at six (n = 4) and 48 (n = 4) hours post-diapause, respectively.

For each haemolymph collection, queens were anaesthetized using carbon dioxide (CO2) (BOC, Ireland) and an incision was performed between the third and fourth sternite. Using this newly made incision, we collected haemolymph using a microcapillary tube and transferred the collected neat haemolymph into a fresh 1.5 ml collection tube containing 100 μl autoclaved phosphate buffered saline (PBS) (Sigma, Ireland). From each bee, we collected approximately 30 μl of neat haemolymph. We centrifuged this solution at 8000 x g at 4 °C for 5 min to pellet down haemocytes and other cellular debris. We transferred the supernatant to a fresh 1.5 ml collection tube and stored at − 80 °C for later use.

Mass spectrometry and data analysis

Haemolymph samples were removed from − 80 °C storage, thawed on ice and biological impurities removed from the samples using the 2-D clean-up kit (GE Healthcare, UK), following the manufacturer’s instructions. The resultant protein pellets were resuspended in 300 μl 50 mM ammonium bicarbonate (Ambic), 1 mM calcium chloride (CaCl2). Protein abundance was then quantified using the Qubit fluorometer with the Qubit protein assay kit (Invitrogen, Ireland). We removed 50 μg from each sample and reduced each with 4 μl of 200 mM dithiothreitol (DTT) (Fisher Scientific, Ireland). We then heated reduced samples to 95 °C for 15 min followed by acetylation with 1 M iodoacetamide (IAA) and incubated at 25 °C for 45 min. The acetylation reaction was stopped by adding 20 μl of 200 mM DTT followed by incubation of the samples at 25 °C for 45 min. We digested 10 μg of reduced and acetylated proteins with 0.5 μl of 0.5 μg/μl trypsin (Medical Supply Company, Ireland) through incubation at 37 °C overnight. The samples were then speed vacuumed and stored at 4 °C once dry.

Dried peptide samples were resolubilised in 0.1% formic acid (Sigma Aldrich, Ireland) and subjected to liquid chromatography tandem mass spectrometry (LC-MS/MS) using a Thermo Scientific LTQ Orbitrap XL mass spectrometer connected to a Dionex Ultimate 3000 (RSLCnano) chromatography system. Peptides were separated by an increasing acetonitrile gradient on a Biobasic C18 Picofrit™ column (100 mm length, 75 mm ID), using a 60 min reverse phase gradient at a flow rate of 300 nL min− 1. All data were acquired with the mass spectrometer operating in automatic data-dependent switching mode. A high-resolution MS scan (300–2000 Da) was performed using the Orbitrap to select the seven most intense ions prior to MS/MS. For each scan performed, the twenty most intense ions were selected for MS/MS analysis.

Protein identification from the MS/MS data was performed using the Andromeda search engine in MaxQuant (version;; [102]) to correlate the data against the protein reference sequences derived from the B. terrestris genome [96] obtained from the National Centre for Biotechnology Information (NCBI) repository (17,508 entries, downloaded September 2018). The following search parameters were used: first search peptide tolerance of 20 ppm (parts per million), second search peptide tolerance of 4.5 ppm (standard parameters for Orbitrap data) with cysteine carbamidomethylation as a fixed modification and N-acetylation of protein, oxidation of methionine as variable modifications and a maximum of two missed cleavage sites allowed, all standard. False discovery rates (FDR) were set to 1% for both peptides and proteins and the FDR was estimated following searches against a target-decoy database. Peptides with a minimum length of seven amino acids were considered for identification and proteins were only considered identified when more than one unique peptide for each protein was observed or when a single peptide was identified across multiple samples by multiple MS/MS spectra. The MS proteomics data and MaxQuant search output files have been deposited to the ProteomeXchange Consortium [103] via the PRIDE partner repository with the dataset identifier PXD014384.

Label-free quantification (LFQ) intensities were calculated using the MaxLFQ algorithm [104] from razor and unique peptides with a minimum ratio count of two peptides across samples, log2-transformed and replicates were assigned to their respective samples. Only proteins found in all four replicates in at least one life stage group were retained in the analysis. Where proteins were absent in one or more replicates of the other groups a data-imputation step was conducted to replace non-detected values with values that simulate signals of low abundant proteins. These imputed values were chosen randomly from the normal distribution of all LFQ intensities for each sample and specified by a downshift of 1.8 times the mean standard deviation (SD) of all measured values and a width of 0.1 times this SD.

Statistical analysis to identify differentially abundant haemolymph proteins

Using Perseus (v. [105]), we performed an analysis of variance (ANOVA) for multiple samples across the six life stage groups using a permutation-based false discovery rate (FDR) of 5% and below to indicate statistically significant differentially abundant (SSDA) proteins. The null hypothesis being tested was that there is no difference in haemolymph proteome prior to, in preparation for, during and after diapause. As a complementary measure to the analyses performed using Perseus, we performed a one-way ANOVA in R (v.3.5.1) using the LFQ intensity values. We identified a significant overlap in proteins detected as differentially abundant by the two approaches (binomial test, p < 10− 10, Additional file 4, Additional file 2: Table S2c) providing independent support for the findings of the analysis by Perseus. As the Perseus software package was specifically developed for the loading, imputation and analysis of proteomic data, we report only these results. Full details on the cross-validation are provided in Additional file 4 and the results are provided in Additional file 8: Table S7. A principal component analysis was performed using Z-score normalised LFQ intensity values for SSDA proteins and output visualised using ggfortify (v.0.4.5 [106]). Hierarchical clustering was performed using Perseus to produce a heatmap of protein abundance using Z-score normalised LFQ intensity values for all SSDA proteins by clustering both proteins using Euclidean distance and complete linkage.

To identify proteomic changes between each consecutive time-point of pre-, mid, and post-diapause samples, Student’s t-tests (p < 0.05) were performed by Perseus with the log2 difference between mean LFQ intensities obtained and used to determine relative abundance differences in individual proteins. Box and line plots for SSDA proteins were generated for specific protein groups to visualise their expression profile before, during and after diapause using ggplot2 (v. 3.1.0 [107]). Gene Ontology term enrichment for SSDA proteins was performed using a Fisher’s exact test (Benjamini Hochberg (BH) adjusted p < 0.05) with the ‘weight01’ algorithm (node size = 20) as part of the topGO package (v. 2.34.0 [108]).

Analysis of a novel haemolymph-associated gene family

For each hypothetical protein identified within the haemolymph, we performed individual BLASTp [109] searches against the NCBI non-redundant (nr) database to identify putative homologues across taxonomically diverse groups to allow for inference of functional annotation. During these searches, we identified certain hypothetical proteins to share close sequence similarity (E-value threshold: 1e-40) to other hypothetical proteins suggestive that the proteins may be isoforms or coded for paralogous genes. To explore this further, we investigated the genomic coordinates of underlying genes to identify the proximity of putative paralogous genes to each other. Genomic coordinates of underlying protein-coding genes were visualised using arrow plots generated by gggenes (v.0.3.2).

Transcript expression of haemolymph-associated proteins during diapause

Using data generated by Amsalem et al. [46], we quantified transcript expression, as well as differential expression of genes coding for haemolymph-associated proteins in the fat bodies of queens collected pre-, during and post-diapause. Additional information on the methods related to transcript quantification and differential expression is provided in Additional file 4. Scripts used for transcript quantification, differential expression and Gene Ontology enrichment analyses are modifications of the scripts published by Colgan et al. [110]. The modified scripts used in the present study are archived at

Availability of data and materials

Data are archived at PRIDE (PXD014384) and will be publicly available upon manuscript acceptance. Analytical scripts for gene expression analysis, Gene Ontology term enrichment and data visualisation are achieved on a public repository at



Ammonium bicarbonate


Antimicrobial peptide


Analysis of variance


Bone morphogenetic protein




False Discovery Rate


Gamma-interferon-inducible lysosomal thiol reductase


Gene Ontology


Highly abundant haemolymph-associated protein




Juvenile hormone


Liquid chromatography tandem mass spectrometry


Label-free quantification


National Centre for Biotechnology Information


Odorant binding protein


Phosphate buffered saline


Principal component analysis


Ejaculatory bulb-specific protein


Peptidoglycan recognition protein


Parts per million


Relative humidity


Standard Deviation


Statistically significant differentially abundant


  1. Tehei M, Franzetti B, Madern D, Ginzburg M, Ginzburg BZ, Giudici-Orticoni M-T, et al. Adaptation to extreme environments: macromolecular dynamics in bacteria compared in vivo by neutron scattering. EMBO Rep. 2004;5:66–70.

    Article  CAS  PubMed  Google Scholar 

  2. Rothschild LJ, Mancinelli RL. Life in extreme environments. Nature. 2001;409:1092–101.

    Article  CAS  PubMed  Google Scholar 

  3. Denlinger DL. Regulation of diapause. Annu Rev Entomol. 2002;47:93–122.

    Article  CAS  PubMed  Google Scholar 

  4. Hand SC, Denlinger DL, Podrabsky JE, Roy R. Mechanisms of animal diapause: recent developments from nematodes, crustaceans, insects, and fish. Am J Physiol Regul Integr Comp Physiol. 2016;310:R1193–211.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Tatar M, Yin C. Slow aging during insect reproductive diapause: why butterflies, grasshoppers and flies are like worms. Exp Gerontol. 2001;36:723–38.

    Article  CAS  PubMed  Google Scholar 

  6. Danks HV. Insect dormancy: an ecological perspective. Entomological Society of Canada; 1987.

    Google Scholar 

  7. Tauber MJ, Tauber CA. Masaki S. Seasonal adaptations of insects: Oxford University Press; 1986.

    Google Scholar 

  8. Kostál V. Eco-physiological phases of insect diapause. J Insect Physiol. 2006;52:113–27.

    Article  PubMed  CAS  Google Scholar 

  9. Sim C, Denlinger DL. Insulin signaling and the regulation of insect diapause. Front Physiol. 2013;4:189.

  10. Emerson KJ, Bradshaw WE, Holzapfel CM. Complications of complexity: integrating environmental, genetic and hormonal control of insect diapause. Trends Genet. 2009;25:217–25.

    Article  CAS  PubMed  Google Scholar 

  11. Flatt T, Tu M-P, Tatar M. Hormonal pleiotropy and the juvenile hormone regulation of Drosophila development and life history. Bioessays. 2005;27:999–1010.

    Article  CAS  PubMed  Google Scholar 

  12. Nakamura A, Miyado K, Takezawa Y, Ohnami N, Sato M, Ono C, et al. Innate immune system still works at diapause, a physiological state of dormancy in insects. Biochem Biophys Res Commun. 2011;410:351–7.

    Article  CAS  PubMed  Google Scholar 

  13. Kubrak OI, Kučerová L, Theopold U, Nässel DR. The sleeping beauty: how reproductive diapause affects hormone signaling, metabolism, immune response and somatic maintenance in Drosophila melanogaster. PLoS One. 2014;9:e113051. .

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Demas GE, Chefer V, Talan MI, Nelson RJ. Metabolic costs of mounting an antigen-stimulated immune response in adult and aged C57BL/6J mice. Am J Phys. 1997;273(5 Pt 2):R1631–7.

    Article  CAS  Google Scholar 

  15. Moret Y, Schmid-Hempel P. Survival for immunity: the price of immune system activation for bumblebee workers. Science. 2000;290:1166–8.

    Article  CAS  PubMed  Google Scholar 

  16. Ardia DR, Gantz JE, Brent C, Schneider, Strebel S. Costs of immunity in insects: an induced immune response increases metabolic rate and decreases antimicrobial activity: Energetic costs of immunity. Funct Ecol. 2012;26:732–9.

    Article  Google Scholar 

  17. Poulsen M, Bot AN, Nielsen MG, Boomsma JJ. Experimental evidence for the costs and hygienic significance of the antibiotic metapleural gland secretion in leaf-cutting ants. Behav Ecol Sociobiol. 2002;52:151–7.

    Article  Google Scholar 

  18. Lord GM, Matarese G, Howard JK, Baker RJ, Bloom SR, Lechler RI. Leptin modulates the T-cell immune response and reverses starvation-induced immunosuppression. Nature. 1998;394:897–901.

    Article  CAS  PubMed  Google Scholar 

  19. Chandra RK. Nutrition and the immune system: an introduction. Am J Clin Nutr. 1997;66:460S–3S.

    Article  CAS  PubMed  Google Scholar 

  20. Benoit JB, Denlinger DL. Suppression of water loss during adult diapause in the northern house mosquito, Culex pipiens. J Exp Biol. 2007;210(Pt 2):217–26.

    Article  PubMed  Google Scholar 

  21. Sota T, Mogi M. Survival time and resistance to desiccation of diapause and non-diapause eggs of temperate Aedes (Stegomyia) mosquitoes. Entomologia Experimentalis et Applicata. 1992;63:155–61.

    Article  Google Scholar 

  22. Hand SC, Moore DS, Patil Y. Challenges during diapause and anhydrobiosis: mitochondrial bioenergetics and desiccation tolerance. IUBMB Life. 2018;70:1251–9.

    Article  CAS  PubMed  Google Scholar 

  23. Pener MP. Environmental cues, endocrine factors, and reproductive diapause in male insects. Chronobiol Int. 1992;9:102–13.

    Article  CAS  PubMed  Google Scholar 

  24. Pitnick S, Marrow T, Spicer GS. Evolution of multiple kinds of female sperm-storage organs in Drosophila. Evol. 1999;53:1804–22.

    Article  PubMed  Google Scholar 

  25. Presgraves DC, Baker RH, Wilkinson GS. Coevolution of sperm and female reproductive tract morphology in stalk--eyed flies. Proc R Soc Lond Ser B Biol Sci. 1999;266:1041–7.

    Article  PubMed Central  Google Scholar 

  26. Baer B, Eubel H, Taylor NL, O’Toole N, Millar AH. Insights into female sperm storage from the spermathecal fluid proteome of the honeybee Apis mellifera. Genome Biol. 2009;10:R67.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Gillott C. Male accessory gland secretions: modulators of female reproductive physiology and behavior. Annu Rev Entomol. 2003;48:163–84.

    Article  CAS  PubMed  Google Scholar 

  28. Poiani A. Complexity of seminal fluid: a review. Behav Ecol Sociobiol. 2006;60:289–310.

    Article  Google Scholar 

  29. Lung O, Kuo L, Wolfner MF. Drosophila males transfer antibacterial proteins from their accessory gland and ejaculatory duct to their mates. J Insect Physiol. 2001;47:617–22.

    Article  CAS  PubMed  Google Scholar 

  30. Markow TA, Ankney PF. Insemination reaction in Drosophila found in species whose males contribute material to oocytes before fertilization. Evol. 1988;42:1097–101.

    PubMed  Google Scholar 

  31. Lawniczak MKN, Begun DJ. A genome-wide analysis of courting and mating responses in Drosophila melanogaster females. Genome. 2004;47:900–10.

    Article  CAS  PubMed  Google Scholar 

  32. McGraw LA, Gibson G, Clark AG, Wolfner MF. Genes regulated by mating, sperm, or seminal proteins in mated female Drosophila melanogaster. Curr Biol. 2004;14:1509–14.

    Article  CAS  PubMed  Google Scholar 

  33. Peng J, Zipperlen P, Kubli E. Drosophila sex-peptide stimulates female innate immune system after mating via the toll and Imd pathways. Curr Biol. 2005;15:1690–4.

    Article  CAS  PubMed  Google Scholar 

  34. Garibaldi LA, Steffan-Dewenter I, Winfree R, Aizen MA, Bommarco R, Cunningham SA, et al. Wild pollinators enhance fruit set of crops regardless of honey bee abundance. Science. 2013;339:1608–11.

    Article  CAS  PubMed  Google Scholar 

  35. Kleijn D, Winfree R, Bartomeus I, Carvalheiro LG, Henry M, Isaacs R, et al. Delivery of crop pollination services is an insufficient argument for wild pollinator conservation. Nat Commun. 2015;6:7414.

  36. Wheeler DE. Developmental and physiological determinants of caste in social Hymenoptera: evolutionary implications. Am Nat. 1986;128:13–34.

    Article  Google Scholar 

  37. Wilson EO. The insect societies: Belknap Press; 1971.

  38. Alford DV. Bumblebees. Davis-Poynter; 1975.

    Google Scholar 

  39. Schmid-Hempel R, Schmid-Hempel P. Female mating frequencies in Bombus spp. from Central Europe. Insect Soc. 2000;47:36–41.

    Article  Google Scholar 

  40. Estoup A, Scholl A, Pouvreau A, Solignac M. Monoandry and polyandry in bumble bees (Hymenoptera; Bombinae) as evidenced by highly variable microsatellites. Mol Ecol. 1995;4:89–94.

    Article  CAS  PubMed  Google Scholar 

  41. Alford DV. A study of the hibernation of bumblebees (Hymenoptera:Bombidae) in southern England. J Anim Ecol. 1969;38:149–70.

    Article  Google Scholar 

  42. Amsalem E, Grozinger CM, Padilla M, Hefetz A. The physiological and genomic bases of bumble bee social behaviour. In: Advances in Insect Physiology. Elsevier; 2015. p. 37–93.

    Google Scholar 

  43. Barribeau SM, Schmid-Hempel P. Sexual healing: mating induces a protective immune response in bumblebees. J Evol Biol. 2017;30:202–9.

    Article  PubMed  Google Scholar 

  44. Brown MJF, Schmid-Hempel R, Schmid-Hempel P. Strong context-dependent virulence in a host--parasite system: reconciling genetic evidence with theory. J Anim Ecol. 2003;72:994–1002.

    Article  Google Scholar 

  45. Larrere M, Lavenseau L, Tasei JN, Couillaud F. Juvenile hormone biosynthesis and diapause termination in Bombus terrestris. Invertebr Reprod Dev. 1993;23:7–14.

    Article  CAS  Google Scholar 

  46. Amsalem E, Galbraith DA, Cnaani J, Teal PEA, Grozinger CM. Conservation and modification of genetic and physiological toolkits underpinning diapause in bumble bee queens. Mol Ecol. 2015;24:5596–615.

    Article  CAS  PubMed  Google Scholar 

  47. Feder ME, Walser J-C. The biological limitations of transcriptomics in elucidating stress and stress responses. J Evol Biol. 2005;18:901–10.

    Article  CAS  PubMed  Google Scholar 

  48. Kanost MR, Kawooya JK, Law JH, Ryan RO, Van Heusden MC, Ziegler R. Insect haemolymph proteins. Advances in Insect Physiology. 1990;:299–396. doi:

    Google Scholar 

  49. Paskewitz SM, Shi L. The hemolymph proteome of Anopheles gambiae. Insect Biochem Mol Biol. 2005;35:815–24.

    Article  CAS  PubMed  Google Scholar 

  50. Vierstraete E, Verleyen P, Baggerman G, D’Hertog W, Van den Bergh G, Arckens L, et al. A proteomic approach for the analysis of instantly released wound and immune proteins in Drosophila melanogaster hemolymph. Proc Natl Acad Sci U S A. 2004;101:470–5.

    Article  CAS  Google Scholar 

  51. Randolt K, Gimple O, Geissendörfer J, Reinders J, Prusko C, Mueller MJ, et al. Immune-related proteins induced in the hemolymph after aseptic and septic injury differ in honey bee worker larvae and adults. Arch Insect Biochem Physiol. 2008;69:155–67.

    Article  CAS  PubMed  Google Scholar 

  52. Chan QWT, Melathopoulos AP, Pernal SF, Foster LJ. The innate immune and systemic response in honey bees to a bacterial pathogen, Paenibacillus larvae. BMC Genomics. 2009;10:387.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Guarna MM, Melathopoulos AP, Huxter E, Iovinella I, Parker R, Stoynov N, et al. A search for protein biomarkers links olfactory signal transduction to social immunity. BMC Genomics. 2015;16:63.

  54. Hou Y, Zou Y, Wang F, Gong J, Zhong X, Xia Q, et al. Comparative analysis of proteome maps of silkworm hemolymph during different developmental stages. Proteome Sci. 2010;8:45.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Li X-H, Wu X-F, Yue W-F, Liu J-M, Li G-L, Miao Y-G. Proteomic analysis of the silkworm (Bombyx mori L.) hemolymph during developmental stage. J Proteome Res. 2006;5:2809–14.

    Article  CAS  PubMed  Google Scholar 

  56. Furusawa T, Rakwal R, Nam HW, Hirano M, Shibato J, Kim YS, et al. Systematic investigation of the hemolymph proteome of Manduca sexta at the fifth instar larvae stage using one- and two-dimensional proteomics platforms. J Proteome Res. 2008;7:938–59.

    Article  CAS  PubMed  Google Scholar 

  57. Chan QWT, Howes CG, Foster LJ. Quantitative comparison of caste differences in honeybee hemolymph. Mol Cell Proteomics. 2006;5:2252–62.

    Article  CAS  Google Scholar 

  58. Lu Y-X, Xu W-H. Proteomic and phosphoproteomic analysis at diapause initiation in the cotton bollworm, Helicoverpa armigera. J Proteome Res. 2010;9:5053–64.

    Article  CAS  PubMed  Google Scholar 

  59. Zhang Q, Lu Y-X, Xu W-H. Proteomic and metabolomic profiles of larval hemolymph associated with diapause in the cotton bollworm, Helicoverpa armigera. BMC Genomics. 2013;14:751.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Mansingh A. Changes in the free amino acids of the haemolymph of Antheraea pernyi during induction and termination of diapause. J Insect Physiol. 1967;13:1645–55.

    Article  CAS  Google Scholar 

  61. Isabel G, Gourdoux L, Moreau R. Changes of biogenic amine levels in haemolymph during diapausing and non-diapausing status in Pieris brassicae L. Comp Biochem Physiol A Mol Integr Physiol. 2001;128:117–27.

    Article  CAS  Google Scholar 

  62. Pérez-Hedo M, Sánchez-López I, Eizaguirre M. Comparative analysis of hemolymph proteome maps in diapausing and non-diapausing larvae of Sesamia nonagrioides. Proteome Sci. 2012;10:58.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Poelchau MF, Reynolds JA, Elsik CG, Denlinger DL, Armbruster PA. Deep sequencing reveals complex mechanisms of diapause preparation in the invasive mosquito, Aedes albopictus. Proc Biol Sci. 2013;280:20130143.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Galperin MY, Koonin EV. “Conserved hypothetical” proteins: prioritization of targets for experimental study. Nucleic Acids Res. 2004;32:5452–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Bouma HR, Carey HV, Kroese FGM. Hibernation: the immune system at rest? J Leukoc Biol. 2010;88:619–24.

    Article  CAS  PubMed  Google Scholar 

  66. Schwenke RA, Lazzaro BP. Wolfner MF. Reproduction–immunity trade-offs in insects. 2016.

    Article  CAS  PubMed  Google Scholar 

  67. Wheeler D. The role of nourishment in oogenesis. Annu Rev Entomol. 1996;41:407–31.

    Article  CAS  PubMed  Google Scholar 

  68. Bulet P, Hetru C, Dimarcq JL, Hoffmann D. Antimicrobial peptides in insects; structure and function. Dev Comp Immunol. 1999;23:329–44.

    Article  CAS  Google Scholar 

  69. Hoffmann JA, Reichhart J-M. Drosophila innate immunity: an evolutionary perspective. Nat Immunol. 2002;3:121–6.

    Article  CAS  PubMed  Google Scholar 

  70. De Gregorio E, Spellman PT, Tzou P, Rubin GM, Lemaitre B. The Toll and Imd pathways are the major regulators of the immune response in Drosophila. EMBO J. 2002;21:2568–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Rees JA, Moniatte M, Bulet P. Novel antibacterial peptides isolated from a European bumblebee, Bombus pascuorum (Hymenoptera, Apoidea). Insect Biochem Mol Biol. 1997;27:413–22.

    Article  CAS  PubMed  Google Scholar 

  72. Williams CM. The juvenile hormone. I. Endocrine activity of the corpora allata of the adult Cecropia silkworm. Biol Bull. 1959;116:323–38.

    Article  Google Scholar 

  73. Schmid-Hempel P. Variation in immune defence as a question of evolutionary ecology. Proc Biol Sci. 2003;270:357–66.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Colgan TJ, Carolan JC, Sumner S, Blaxter ML, Brown MJF. Infection by the castrating parasitic nematode Sphaerularia bombi changes gene expression in Bombus terrestris bumblebee queens. Insect Mol Biol. 2019.

  75. Poinar GO, Van Der Laan PA. Morphology and life history of Sphaerularia bombi. Nematologica. 1972;18:239–52.

    Article  Google Scholar 

  76. Cotter SC, Kruuk LEB, Wilson K. Costs of resistance: genetic correlations and potential trade-offs in an insect immune system. J Evol Biol. 2003;17:421–9.

    Article  Google Scholar 

  77. Lochmiller RL, Deerenberg C. Trade-offs in evolutionary immunology: just what is the cost of immunity? Oikos. 2000;88:87–98.

    Article  Google Scholar 

  78. Cotter SC, Simpson SJ, Raubenheimer D, Wilson K. Macronutrient balance mediates trade-offs between immune function and life history traits. Funct Ecol. 2011;25:186–98.

    Article  Google Scholar 

  79. Gibbs AG. Water-proofing properties of cuticular lipids. Integr Comp Biol. 1998;38:471–82.

    Article  CAS  Google Scholar 

  80. Li A, Denlinger DL. Pupal cuticle protein is abundant during early adult diapause in the mosquito Culex pipiens. J Med Entomol. 2009;46:1382–6.

    Article  CAS  PubMed  Google Scholar 

  81. Andersen SO, Hojrup P, Roepstorff P. Insect cuticular proteins. Insect Biochem Mol Biol. 1995;25:153–76.

    Article  CAS  PubMed  Google Scholar 

  82. Denlinger DL. Hormonal control of diapause. Endocrinology II. 1985;:353–412. doi:

    Chapter  Google Scholar 

  83. Pelosi P, Maida R. Odorant-binding proteins in insects. Comp Biochem Physiol B Biochem Mol Biol. 1995;111:503–14.

    Article  CAS  Google Scholar 

  84. Foret S, Maleszka R. Function and evolution of a gene family encoding odorant binding-like proteins in a social insect, the honey bee (Apis mellifera). Genome Res. 2006;16:1404–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Zhou J-J. Odorant-binding proteins in insects. Vitamins & Hormones. 2010;:241–72. doi:

    Google Scholar 

  86. Lung O, Wolfner MF. Identification and characterization of the major Drosophila melanogaster mating plug protein. Insect Biochem Mol Biol. 2001;31:543–51.

    Article  CAS  PubMed  Google Scholar 

  87. Ludwig MZ, Uspensky II, Ivanov AI, Kopantseva MR, Dianov CM, Tamarina NA, et al. Genetic control and expression of the major ejaculatory bulb protein (PEB-me) in Drosophila melanogaster. Biochem Genet. 1991;29:215–39.

    Article  CAS  PubMed  Google Scholar 

  88. Sabatier L, Jouanguy E, Dostert C, Zachary D, Dimarcq J-L, Bulet P, et al. Pherokine-2 and-3: Two Drosophila molecules related to pheromone/odor-binding proteins induced by viral and bacterial infections. Eur J Biochem. 2003;270:3398–407.

    Article  CAS  PubMed  Google Scholar 

  89. Vogt RG, Rybczynski R, Cruz M, Lerner MR. Ecdysteroid regulation of olfactory protein expression in the developing antenna of the tobacco hawk moth, Manduca sexta. J Neurobiol. 1993;24:581–97.

    Article  CAS  PubMed  Google Scholar 

  90. Yocum GD, Rinehart JP, Horvath DP, Kemp WP, Bosch J, Alroobi R, et al. Key molecular processes of the diapause to post-diapause quiescence transition in the alfalfa leafcutting bee Megachile rotundata identified by comparative transcriptome analysis. Physiol Entomol. 2015;40:103–12.

    Article  CAS  Google Scholar 

  91. Law JH, Wells MA. Insects as biochemical models. J Biol Chem. 1989;264:16335–8.

    CAS  PubMed  Google Scholar 

  92. Arrese EL, Soulages JL. Insect fat body: energy, metabolism, and regulation. Annu Rev Entomol. 2010;55:207–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  93. Haider S, Pal R. Integrated analysis of transcriptomic and proteomic data. Curr Genomics. 2013;14:91–110.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  94. Ghazalpour A, Bennett B, Petyuk VA, Orozco L, Hagopian R, Mungrue IN, et al. Comparative analysis of proteome and transcriptome variation in mouse. PLoS Genet. 2011;7:e1001393.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Wilhelm M, Schlegl J, Hahne H, Gholami AM, Lieberenz M, Savitski MM, et al. Mass-spectrometry-based draft of the human proteome. Nature. 2014;509:582–7.

    Article  CAS  PubMed  Google Scholar 

  96. Sadd BM, Barribeau SM, Bloch G, de Graaf DC, Dearden P, Elsik CG, et al. The genomes of two key bumblebee species with primitive eusocial organization. Genome Biol. 2015;16:76.

  97. Lin G, Huang Z, Wang L, Chen Z, Zhang T, Gillman LN, et al. Evolutionary rates of bumblebee genomes are faster at lower elevations. Mol Biol Evol. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  98. Barribeau SM, Sadd BM, du Plessis L, Brown MJF, Buechel SD, Cappelle K, et al. A depauperate immune repertoire precedes evolution of sociality in bees. Genome Biol. 2015;16:83.

  99. Colgan TJ, Carolan JC, Bridgett SJ, Sumner S, Blaxter ML, Brown MJ. Polyphenism in social insects: insights from a transcriptome-wide analysis of gene expression in the life stages of the key pollinator, Bombus terrestris. BMC Genomics. 2011;12:623.

  100. Sauter A, Brown MJF. To copulate or not? The importance of female status and behavioural variation in predicting copulation in a bumblebee. Anim Behav. 2001;62:221–6.

    Article  Google Scholar 

  101. Beekman M, Stratum P, Veerman A. Selection for non-diapause in the bumblebee Bombus terrestris, with notes on the effect of inbreeding. Entomol Exp Appl. 1999;93:69–75.

    Article  Google Scholar 

  102. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26:1367–72.

    Article  CAS  PubMed  Google Scholar 

  103. Côté RG, Griss J, Dianes JA, Wang R, Wright JC, van den Toorn HWP, et al. The PRoteomics IDEntification (PRIDE) converter 2 framework: an improved suite of tools to facilitate data submission to the PRIDE database and the ProteomeXchange consortium. Mol Cell Proteomics. 2012;11:1682–9.

    Article  CAS  Google Scholar 

  104. Cox J, Hein MY, Luber CA, Paron I, Nagaraj N, Mann M. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol Cell Proteomics. 2014;13:2513–26.

    Article  CAS  Google Scholar 

  105. Tyanova S, Cox J. Perseus: a bioinformatics platform for integrative analysis of proteomics data in cancer research. Methods Mol Biol. 1711;2018:133–48.

    Google Scholar 

  106. Tang Y. Horikoshi M. Li W ggfortify: unified interface to visualize statistical result of popular R packages RJ. 2016;8:474–85.

    Article  Google Scholar 

  107. Wickham H. ggplot2: Elegant graphics for data analysis. Springer Science & Business Media; 2009.

  108. Alexa A, Rahnenfuhrer J. topGO: Enrichment analysis for Gene Ontology. R package version 2.28. 0. Cranio. 2016.

  109. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    Article  CAS  PubMed  Google Scholar 

  110. Colgan TJ, Fletcher IK, Arce AN, Gill RJ, Rodrigues AR, Stolle E, et al. Caste- and pesticide-specific effects of neonicotinoid pesticide exposure on gene expression in bumblebees. Mol Ecol. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We acknowledge Ms. Heather Cuddy for help with the maintenance of animals. We also thank Mr. Peter Stafford and Ms. Alison Boyce for technical assistance throughout the project. We thank Prof. Christine Griffin (Maynooth University) for beneficial comments on an earlier draft and Mr. Pedro Busana for assistance in the generation of Fig. 5. We would also like to extend our gratitude to the Editor and two anonymous reviewers that provided helpful comments, which greatly improved this manuscript.


This publication has emanated from research conducted with the financial support of Science Foundation Ireland ( grant EEEOBF131 to MJFB. This funding was used to support sample collection and mass-spectrometry analysis of haemolymph samples.

Author information

Authors and Affiliations



TJC, MJFB and JCC conceived and designed the experiment. TJC, SF and JCC performed the experiments. TJC and JCC analysed the data. All authors wrote, read and have approved the final manuscript.

Corresponding author

Correspondence to Thomas J. Colgan.

Ethics declarations

Ethics approval and consent to participate

This work did not require any ethical approval and complies with the Animals Scientific Procedures Act 1986.

Consent for publication

Not applicable to this present study.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

The 129 proteins identified from the haemolymph of unmated and mated queens prior to, during and after diapause.

Additional file 2: Table S2.

Multivariate analysis (ANOVA FDR 5%) of all identified proteins from the haemolymph of unmated and mated queens prior to, during and after diapause.

Additional file 3: Table S3.

All ANOVA statistically significant proteins (FDR < 0.05) used for hierarchical clustering.

Additional file 4: Figure S1.

Genomic coordinates of members of a putative expanded haemolymph-associated gene family within the buff-tailed bumblebee genome. Figure S2. One HAHP family member increases in abundance during diapause. Figure S3. Gene Ontology term enrichment within significantly differentially abundant proteins in the bumblebee queen haemolymph proteome.

Additional file 5: Table S4.

Genomic coordinates for members of a novel haemolymph-associated gene family.

Additional file 6: Table S5

. Gene expression of the haemolymph-associated proteins in the bumblebee queen fat body.

Additional file 7: Table S6.

Enriched Gene Ontology terms in genes coding for haemolymph-associated differentially expressed proteins during the bumblebee queen life cycle stages.

Additional file 8: Table S7.

Pairwise t-tests were performed for each protein between consecutive life-cycle stages.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Colgan, T.J., Finlay, S., Brown, M.J.F. et al. Mating precedes selective immune priming which is maintained throughout bumblebee queen diapause. BMC Genomics 20, 959 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: