- Open Access
Evidence for a chemical arms race between cuckoo wasps of the genus Hedychrum and their distantly related host apoid wasps
BMC Ecology and Evolution volume 22, Article number: 138 (2022)
Brood parasites can exert strong selection pressure on their hosts. Many brood parasites escape their detection by mimicking sensory cues of their hosts. However, there is little evidence whether or not the hosts are able to escape the parasites’ mimicry by changing these cues. We addressed this question by analyzing cuticular hydrocarbon (CHC) profiles of Cerceris and Philanthus wasps and their brood parasites, cuckoo wasps mimicking the CHC profiles of their hosts. Some of these hosts use hydrocarbons to preserve their prey against fungal infestation and thus, they cannot significantly change their CHC composition in response to chemical mimicry by Hedychrum brood parasites.
We found that the CHC overlap between brood parasites and their hosts was lower in case of host wasps not preserving their prey than in case of prey-preserving host wasps, whose CHC evolution is constrained. Furthermore, the CHC profiles in non-preserving host wasps is more strongly diversified in females than in males, thus in the sex that is chemically mimicked by brood parasites.
Our results provide evidence for a chemical arms race between those hosts that are liberated from stabilizing selection on their chemical template and their parasites.
Coevolution between interacting species is considered one of the major forces generating biological diversity . Coevolution plays an important role in the organization of communities, by shaping, for example, both symbiotic and parasitic interactions and reciprocal specialization among free-living taxa [2, 3], and in promoting the appearance of key innovations [4, 5]. Yet, demonstrating that a given trait has evolved in response to another trait from a different species and that the evolved trait caused counter-adaptations in return is not straightforward. It is particularly difficult to rule out effects of external selection pressures on a particular trait of interest .
Among antagonistic interactions, interspecific brood parasitism is a widespread strategy of kleptoparasites and parasitoids (subsequently collectively referred to as brood parasites). Brood parasites exploit their host’s provision for offspring (in case of kleptoparasites) or the immature stages of the host’s offspring itself (in case of parasitoids) to raise their own offspring [6,7,8,9]. One of the best studied examples of coevolution between brood parasites and their hosts is that between cuckoos (Aves: Cuculidae) and their hosts. Here, hosts of cuckoos evolved various strategies to compromise optical mimicry of host eggs by cuckoos, such as increased egg color polymorphism [6, 10]. For such antagonistic coevolutionary relationships, Van Valen’s Red Queen Hypothesis provides a theoretical framework that explains the evolutionary arms race between hosts and their parasites .
Olfactory communication plays a fundamental role for insects . As virtually all insects are covered by cuticular hydrocarbons (CHCs), this diverse class of semiochemicals has gained particular importance for communicating important information, such as species identity, age, sex, reproductive status, and—in eusocial species—caste and colony membership . It therefore comes as no surprise that some brood parasitic insects evolved the ability to chemically deceive their hosts in order to conceal their detrimental action.
Insect brood parasites are known to adopt one of at least three different strategies to successfully deceive their host’s olfactory recognition system : (1) the brood parasite synthesizes a low total amount and/or a simplified blend of hydrocarbons on its cuticle and thus provides few recognition cues to the host [15,16,17]; this is called the chemical insignificance strategy. (2) The brood parasite physically adopts the host’s CHC profile by grooming; this strategy is referred to as chemical camouflage [18, 19]. (3) The brood parasite synthesizes de novo a CHC profile that is very similar to that of its host ; this strategy is known as chemical mimicry. These three strategies are not mutually exclusive, and some brood parasites apply them simultaneously. For example, the butterfly Phengaris (Maculinea) rebeli, whose larvae develop in ant nests, applies both a chemical mimicry and a chemical camouflage strategy to deceive its host .
Hosts can counteract the chemical deception strategies of their brood parasites by evolving an improved ability to recognize and to discriminate chemical cues and/or by expressing chemical phenotypes that differ from those of their brood parasites (e.g., via negative frequency-dependent selection; [21, 22]. However, the evolutionary response by hosts to chemical mimicry by their brood parasites has received comparatively little attention so far.
Despite the importance of chemical mimicry for insect brood parasites, very few studies have provided evidence for this strategy to result in a coevolutionary arms race between brood parasites and their hosts. Examples include slave-making ants [23,24,25], bumble bees , vespid wasps , and cuckoo wasps and their hosts [21, 28]. In most of the studied cases, the brood parasitic species apply a chemical camouflage strategy (i.e., ants, bumble bees, vespid wasps), often in association with the insignificance strategy. Furthermore, most of these cases involve host-brood parasite species pairs that are phylogenetically closely related, a phenomenon known as Emery’s rule . Clear evidence for chemical mimicry is currently only known from host-brood parasite species pairs that are phylogenetically distantly related. Here, the brood parasites have been forced to evolve a chemical profile that matches that of their hosts [20, 21].
To shed light on chemical adaptations in host-brood parasite species pairs, we investigated coevolutionary patterns in CHC profiles across species, using cuckoo wasps of the genus Hedychrum (Hymenoptera: Chrysididae: Elampini) and their distantly related hosts, apoid wasps of the genera Cerceris and Philanthus (Hymenoptera: Philanthidae) [30,31,32,33], as models. Females of Cerceris and Philanthus species dig brood cells in the ground and provision their larvae with paralyzed prey, which consists either of beetles (Coleoptera-hunting wasps: COLw) or of stinging wasps and bees (Hymenoptera-hunting wasps: HYMw).
What makes the above mentioned hosts of Hedychrum cuckoo wasps particularly interesting is the fact that females of HYMw (but not of COLw) are known to embalm their prey with a secretion from their postpharyngeal gland to delay/prevent fungal infestation of the prey [34,35,36,37]. These secretions consist primarily of unsaturated long-chain hydrocarbons (e.g., alkenes), which form a hydrophobic oily layer on the prey that prevents water condensation and impairs mold development on the prey [38, 39]. Moreover, the hydrocarbon composition of the postpharyngeal gland content of female HYMw strongly matches the CHC composition of the wasps’ own cuticle , similar to what was previously observed in ants . As the CHC composition of the postpharyngeal gland content is highly similar among species, the specific composition likely represents an adaptation to preserve prey and the composition is thus likely under stabilizing selection . In contrast, females of COLw do not embalm their prey with a secretion of CHCs, because their prey is less susceptible to fungal growth . Females of HYMw consequently cannot alter their CHC profile without losing the ability to preserve their prey. In contrast, females of COLw do not seem to be constrained by stabilizing selection on their CHC profile and thus have the ability to shift their CHC profiles to avoid chemical mimicry by their brood parasites .
In this study, we report results from testing three hypotheses: (1) CHC profile mimicry is less precisely accomplished by cuckoo wasps that are brood parasites of COLw than by those that are brood parasites of HYMw, because COLw have more abilities to alter their CHC profiles in response to chemical mimicry than HYMw have. (2) Only female cuckoo wasps apply a chemical mimicry strategy, because only females benefit from not leaving detectable traces of their intrusion in the host nest. (3) Possible chemical counter-adaptations by the hosts should primarily affect females, whose nests are intruded by cuckoo wasps. This is because only female hosts are under selective parasite pressure to evolve counter-adaptations to chemical mimicry by cuckoo wasps. To be better able to interpret possible adaptations and counteradaptations in the CHC profiles of hosts and parasites, we considered the phylogenetic history of the investigated species by inferring the phylogeny for the genus Hedychrum and considering a published phylogenetic tree of the hosts.
Molecular phylogeny and evolutionary relationships
Our inferred cuckoo wasp phylogeny shows that the species parasitizing COLw (Hedychrum chalybaeum, Hedychrum niemelai, and Hedychrum nobile) and the species parasitizing HYMw (Hedychrum gerstaeckeri, Hedychrum longicolle, and Hedychrum rutilans) each form a natural group. In contrast, the host species that hunt Hymenoptera are paraphyletic with respect to the host species that hunt Coleoptera  (Fig. 1).
Cuticular hydrocarbon profile composition
In total, we identified 112 different CHCs across the two sexes of the twelve species (hosts and brood parasites) (N = 277; Additional file 2: Tables S4,S5). Among the 112 CHCs, twelve are linear alkanes, 30 are alkenes, nine are alkadienes, 39 are monomethyl- and 22 are dimethyl-branched alkanes. Linear alkanes constitute 30–45% of the total CHC amount in each species (except in females of H. niemelai, whose alkanes make up 64% of the total CHC amount). Unsaturated compounds comprise up to 65% of a given CHC profile of Hymenoptera-hunting wasps/hosts (HYMw), whereas the proportion of methyl-branched alkanes in HYMw is less than 3% of the total CHCs (Table 1). We found no dimethyl-branched alkanes in any HYMw. In contrast, methyl-branched alkanes range from 20% (males of Cerceris interrupta) to more than 60% of the CHC profile of COLw (Fig. 2).
The CHC profiles of female HYMw are very similar to each other across species and are predominantly composed of alkenes (60–65%) and linear alkanes (34–37%). Males of HYMw species are similar in their CHC substance class composition to their female conspecifics (Fig. 2), but they systematically differ from conspecific females by a high relative abundance of alkenes that exhibit their double bonds at positions different from those of alkenes of the females. In contrast, females and males of COLw are more dissimilar in their composition of alkenes (0–30% and 8–38% in females and males, respectively) and methyl-branched alkanes (31–61% and 20–55% in females and males, respectively). Females of COLw possess a larger proportion of methyl-branched alkanes in their CHC profiles than their conspecific males (50% vs. 38%, t(54) = 3.03, p = 0.049; Table 3; Additional file 1: Fig. S1c). In general, both the CHC profiles of female and male HYMw possess less diverse CHC profiles (with an average number of 24.7 ± 3.49 and 21.9 ± 4.46 CHC compounds in the profiles of females and males, respectively) than CHC profiles of female and male COLw (average number of 30.0 ± 6.17 and 34.5 ± 3.65 CHC compounds in the profiles of females and males, respectively). The lower number of compounds in HYMw than in COLw is due to the very low number of methyl-branched compounds that HYMw express (on average fewer than six compounds, Additional file 2: Tables S4,S5). Among COLw, only Cerceris arenaria females possess more compounds in their CHC profiles than males (35.9 ± 1.13 and 31 ± 1.26 in females and males, respectively; t(17) = 8.8427, p < 0.001). In the two other COLw species (i.e., Cerceris interrupta and Cerceris quinquefasciata), males possess a larger number of CHC compounds in their CHC profiles than their conspecific females (Table 3; Additional file 1: Fig. S1a,b).
The profiles of cuckoo wasps contain all CHC substance classes found in their hosts. Species that parasitize HYMw feature larger proportions of alkenes compared to species parasitizing COLw (t(108) = − 4.99, p < 0.001). Species that parasitize COLw tend to show larger proportions of methyl-branched alkanes than HYMw except for H. chalybaeum (t(88) = − 2.88, p = 0.059). The latter species exhibits larger proportions of alkenes (> 60% of its CHC profile) similar to Hedychrum species parasitizing HYMw. Female brood parasites of COLw produce high relative amounts (14–48%) of methyl-branched alkanes except for H. chalybaeum (< 0.5%), while conspecific males produce high relative amounts (51–62%) of alkenes, similar in their proportion to the proportion of alkenes in the profiles of female and male brood parasites of HYMw (Fig. 2 and Table 1). Female and male brood parasites of HYMw differ from each other in the type of synthesized alkenes: females synthesize primarily alkenes with double bonds at position 7 and position 9, while conspecific males synthesize mainly alkenes with double bonds at position 9 and position 11.
The average CHC chain length in COLw in both sexes is approximately two carbon atoms longer than that in HYMw and in all investigated brood parasites (Fig. 3). Both females and males of COLw species synthesize significantly less relative amounts of hydrocarbons with chain lengths of fewer than 24 carbon atoms than HYMw (t(97) = − 16.7, p < 0.001). In general, females of hosts and parasites synthesize CHCs with a longer mean chain length than their males (t(259) = 3.58, p = 0.007). When investigating different guilds (COLw, HYMw, and brood parasites) separately, this pattern was found when comparing males and females of COLw hosts (t(54) = 3.23, p = 0.029) and when comparing males and females of brood parasites (t(96) = 7.99, p < 0.001). This pattern was not found when comparing males and females of HYMw (t(50) = 1.47, p = 0.7).
Chemical profile overlap between brood parasites and hosts
NMDS revealed a clear separation of COLw and HYMw based on CHC profile data, with the two host types forming distinct groups with little overlap of their chemical space (Fig. 4a, ANOSIM R: 0.943, mean rank distance within groups: 685, mean rank distance between groups: 1890, p = 0.0029, 9999 permutations). This separation mainly results from differences in the relative amounts of methyl-branched compounds and of unsaturated compounds in both types of hosts (i.e., COLw and HYMw). Hedychrum nobile and—to a lesser extent—H. niemelai produce large amounts of methyl-branched compounds and are therefore separated from the other cuckoo wasps by occupying a chemical space that falls between that of their hosts and that of HYMw. In contrast, all HYMw and their brood parasites occupy a comparatively small range of chemical space and exhibit overlapping chemical signatures (i.e., they are quite similar to each other, Fig. 4a). An R value close or identical to 1 indicates no overlap between groups. Indeed, the R statistic of ANOSIM is larger than 0.99 in all COLw-brood parasites comparisons, whereas it ranges between 0.66 and 0.93 when comparing HYMw and their brood parasites. A different pattern was observed when the CHC profiles of male brood parasites were plotted with those of their female hosts: males of cuckoo wasps that parasitize COLw species are chemically similar among each other and separated from their COLw hosts (Fig. 4b, ANOSIM among male brood parasites of COLw hosts, R: 0.4873, p = 0.0029; ANOSIM between male brood parasites of COLw and their female COLw hosts, R: 1, p = 0.0029), whereas female brood parasites of COLw are very distinct from each other (ANOSIM R: 0.9686, p = 0.0029). With the exception of H. rutilans, male cuckoo wasps overlap more in their chemical profiles among each other than their females do (ANOSIM between female brood parasites of COLw and HYMw, R: 0.4524, p = 0.0029; ANOSIM between male brood parasites of COLw and HYMw, R: 0.7184, p = 0.0029). The overlap between female COLw and their female brood parasites is larger than that between female COLw and male parasites (ANOSIM between female COLw and their female brood parasites, R: 0.5873, p = 0.0029; ANOSIM between female COLw hosts and their male brood parasites, R: 0.9349, p = 0.0029). However, when analyzing HYMw, the chemical space overlap between female brood parasites of HYMw and their female hosts is similar than that between male brood parasites of HYMw and female hosts (ANOSIM between female HYMw and their female brood parasites, R: 0.4297, p = 0.0029; ANOSIM between female HYMw hosts and their male brood parasites, R: 0.4651, p = 0.0029).
Female host/female brood parasite chemical distances were almost always smaller than female host/male brood parasite chemical distances (Bray–Curtis dissimilarity index calculated across all compounds present). Exceptions were comparisons between the two Philanthus species and H. rutilans (Fig. 5). We did not find a difference between the chemical distances of female HYMw hosts to their female brood parasites and between the chemical distances of female COLw hosts to their female parasites (t(820) = − 1.36, p = 0.77). Chemical distances of COLw female hosts to their female brood parasites were smaller than the chemical distances to their male brood parasites (0.58 vs. 0.68, t(643) = − 10, p < 0.001). However, chemical distances of HYMw female hosts to their female brood parasites were similar to those of their male brood parasites (0.57 vs. 0.58; t(1,191) = 1.84, p = 0.5) (Fig. 5).
Intra- and interspecific CHC variability of HYMw and COLw hosts
We observed a pattern indicating differences in inter- and intraspecific variability of CHC profiles in the NMDS plot (Fig. 4a). To further assess this pattern, we calculated Bray–Curtis dissimilarities between and within species. Interspecific differences were significantly larger between females of COLw hosts than between females of HYMw hosts (average 0.57 ± 0.09 vs. 0.47 ± 0.17; Bray–Curtis distance; t(897) = − 11.39, p < 0.001). In contrast, intraspecific variability was smaller between females of COLw hosts than between females of HYMw hosts (on average, Bray Curtis distances between samples within COLw species were 0.17 ± 0.08; those within HYMw species were 0.32 ± 0.19; t(334) = 10.83, p < 0.001). Hence, in contrast to females of COLw host species, females of HYMw host species show a larger intraspecific profile variance and thus a much larger dissimilarity between individuals belonging to the same species. This trend was not observed in conspecific males (Fig. 6). Males of COLw and HYMw show a smaller intraspecific variability relative to conspecific females. To test whether these results are not an artifact of sample size differences among species, we repeated the analyses using only five samples per species (the maximum number of samples consistently available across all analyzed species) and obtained comparable results (see Additional file 1: Fig. S2). In addition, we analyzed randomly picked CHC profiles of two females of a HYMw host species and of two females of COLw host species, with each pair of samples coming from the same population (same species and same locality of collection). We found that intraspecific differences in HYMw females are significantly larger than in COLw females [t(1324) = 22.2, p < 0.001]. By contrast, when repeating the analysis using the CHC profiles of two HYMw males and of two COLw males, no such differences were observed (see Additional file 1: Fig. S3).
The results of our study confirm the existence of clear CHC profile composition differences between COLw and HYMw, first suggested by Wurdack and colleagues . HYMw species exhibit profiles that consist primarily of unsaturated compounds and that are more similar among each other than those of COLw species. In contrast, COLw species possess CHC profiles with a large number of mono- and dimethyl-branched alkanes, and these compounds also represent a substantial fraction of the total CHC amount. The compositions of the CHC profiles of all investigated species are species- and sex-specific and heritable due to their genetic basis of their biosynthetic pathways [13, 42]. Although we know that CHC profiles are plastic in some species, we can rule out plasticity as an explanation for the CHC differences to a large extent as we collected all individuals from the same area with very similar environmental conditions.
CHC diversification in COLw has led to less chemical overlap between COLw and their brood parasites
Our first hypothesis stated that brood parasites of COLw species should chemically mimic their hosts’ CHC profile less precisely than brood parasites of HYMw. This is because COLw can change their CHC profile to a large extent since it is less constrained due to the abandonment of prey embalming, whereas HYMw seemingly cannot alter their CHC composition significantly at all, the latter thereby representing an easy target to chemically mimic. Brood parasites of HYMw should thus show a strong overlap in their CHC profile composition with that of their host, since HYMw are constrained to produce a relatively fixed and alkene-enriched CHC profile to cope with prey embalming. The abandonment of prey embalming behavior, which according to  likely evolved later in COLw, may have liberated COLw from producing alkene-enriched CHC profiles by relaxation from stabilizing selection.
We showed that COLw hosts differ in a number of features from HYMw hosts, all of which could be attributed to adaptations to escape mimicry by cuckoo wasps: (a) the biosynthesis of a larger number and higher amounts of methyl-branched CHCs, (b) the production of CHC compounds of higher chain length and reduction of shorter chain-length compounds, and (c) smaller intraspecific variation of CHC profiles in comparison to HYMw hosts (Fig. 6a). Moreover, the high between species dissimilarity observed in COLw could have arisen to avoid overlapping CHC profiles with their parasites but higher interspecific differences between COLw (Fig. 6a) to avoid overlapping CHC profiles with their parasites. As expected, there is a larger CHC profile overlap between HYMw and their female brood parasites but there is no overlap of CHC profiles between COLw and their female brood parasites. In contrast to our expectation, however, Bray–Curtis dissimilarities between HYMw and their brood parasites were not smaller than those between COLw female hosts and their brood parasites.
Female brood parasites of HYMw produce similar amounts of the same compound classes as their hosts . However, differences in the double bond positions of alkenes in hosts and brood parasites contribute to larger chemical distances between the CHC profiles of female HYMw and their female brood parasites and can be interpreted as a strategy applied by HYMw to escape chemical mimicry. In addition, CHC profiles of female brood parasites of HYMw (i.e., H. gerstaeckeri and H. rutilans) contain low amounts (~ 5%) of compounds (e.g., methyl-branched compounds) that are absent in the CHC profile of their hosts and which contributes to the relatively large chemical distances observed between female HYMw and their female brood parasites. However, whether CHC compounds occurring at low abundances on the cuticle of HYMw brood parasites can be recognized by their hosts is unknown, especially when considering that these cues may interfere with scents of the environment (e.g., nest material, prey items provided).
Interestingly, although there is no overlap in the chemical space of the CHC profiles between COLw and their associated brood parasites, female brood parasites of COLw hosts (with the exception of H. chalybaeum) often produce the same methyl-branched compounds as their hosts (or at least methyl-branched compounds with the same branching position). The similarities in methyl-branched compounds between hosts and brood parasites are mainly responsible for the relatively short chemical distances between COLw hosts and their cuckoo wasps. A diversification of the CHC profiles with methyl-branched compounds may have resulted from selection on the hosts. In fact, some studies suggest that CHC profiles with more methyl-branched alkanes prevail in highly parasitized populations of other Hymenoptera: for example, parasitized ant colonies of Formica fusca show a higher diversity of dimethyl-branched alkanes than non-parasitized ones, and this increase in compound diversity correlates with increased recognition abilities in the host populations . Likewise, the proportion of methyl-branched hydrocarbons in CHC profiles of the paper wasp Polistes biglumis is larger in populations highly parasitized by the social parasite Polistes atrimandibularis .
Apart from the diversification of the CHC profile in COLw, we found that the mean chain length of CHCs in COLw profiles is significantly larger than in the profiles of HYMw and those of all investigated brood parasites. Thus, COLw and their brood parasites differ in mean chain length from each other, while HYMw and their parasites do not. Several studies showed that shifts in chain-length in CHC profiles can be due to changes in climatic conditions (e.g., [45,46,47]). However, since both types of hosts inhabit similar habitats (warm/dry conditions) and even co-occur in mixed nest aggregations—at least in the regions where we studied the species—the differences in the CHC chain lengths between the two groups of hosts are unlikely the result of differential adaptation to climatic conditions. It appears more likely that HYMw species require alkenes of a specific chain length, because this chain length allows maintaining a semifluid texture in the secretions that are spread on their prey against fungal infestation . Thus, it is possible that the elimination of short-chain compounds in COLw may represent an additional escape strategy from chemical mimicry, as this increases the qualitative differences between the CHC profiles of hosts and parasites. A shift of the mean chain length could be beneficial if it allows the hosts to better discriminate between otherwise similar CHC profiles and thus allows them to detect differences between the CHC profiles of a potential brood parasite and their own. Moreover, an increase of the mean chain length of CHC would likely not require the evolution of new enzymatic pathways .
As mentioned above, the CHC profile similarity between HYMw and their brood parasites was lower than expected. We hypothesize that HYMw host counteract chemical mimicry of their brood parasites by applying an alternative evasion strategy. Chemically, HYMw are restrained to maintain the same proportion of alkenes in their CHCs because of the advantages alkenes confer to prey preservation. However, the CHC profiles of females of HYMw have a larger within-species variability than those of the females of COLw. In case of HYMw, negative frequency-dependent selection may favor the existence of rare host chemotypes. Highly variable CHC profiles could result in some individuals in a population exhibiting chemotypes that are likely not being well mimicked by a brood parasite. A similar phenomenon (i.e., increasing within clutch variation) has been observed in hosts of avian cuckoos : females of some host species of cuckoos have evolved the ability to lay eggs with different color hues. It has been suggested that hosts of cuckoo that are in an evolutionary arms race with their brood parasites should increase inter-clutch variation (i.e., differences among eggs within a population, eggs laid by different females), but at the same time should keep intra-clutch color variation low  in order to make mimicry more difficult for the brood parasites. Similar adaptations for increasing within species variation in chemical signals have been observed in insects. For example, within population variation in the relative proportion of different CHCs is higher in highly parasitized populations of Polistes biglumis than in non-parasitized populations of this species, probably because of negative frequency-dependent selection of rare phenotypes . Similarly, when the CHCs of Temnothorax longispinosus ants are compared within and between colonies in populations with and without the slave-making ant Protomognathus americanus, the CHC profiles of the host ant species are more variable in parasitized populations . We have shown that two randomly chosen CHC profiles of a HYMw female host coming from the same population are more dissimilar to each other than two randomly chosen CHCs of a COLw host (Additional file 1: Fig. S3). CHC variation in HYMw females is largely quantitative. For example, in two out of nine female individuals of Philanthus triangulum in our analysis, the main compound being produced was (Z)-9-C27:1, whereas in the remaining individuals, (Z)-9-C25:1 was more abundant (see Fig. 4a). This polymorphic variation in the production of one or the other alkene in females of P. triangulum within the same population had already been reported before [28, 50]. It would be interesting to test the hypothesis that this chemical polymorphism is an adaptation to partly escape chemical mimicry by H. rutilans.
CHC profile diversification is more pronounced in females than in males of COLw hosts
We hypothesized that females of COLw hosts responded stronger to selection pressure exerted by the brood parasites on their CHC profile and thus show a greater CHC profile diversity than males of COLw hosts. If COLw species are able to evolve CHC profiles that differ from that of their brood parasites, we should see this pattern more pronounced in females than in males, assuming that CHC profiles evolve sex-specifically. While CHC profiles of the analyzed hosts do not show striking sex-specific qualitative differences, irrespective of whether we analyzed COLw or HYMw, we found quantitative differences between the CHC profiles of females and of males. Specifically, females of COLw are characterized by possessing larger proportions of methyl-branched alkanes than their conspecific males on their cuticle. In the COLw species C. arenaria, females furthermore synthesize a larger number of methyl-branched compounds than conspecific males. Male-exclusive methyl-branched compounds that contribute to CHC diversity in male COLw of C. interrupta and C. quinquefasciata (Additional file 1: Fig. S1a,b) were not produced in large amounts and made up less than 3% of the total CHC profile in males. Overall, the CHC profile differences between males and females are consistent with the idea with female COLw hosts escaping chemical mimicry of their brood parasites.
CHC profiles of female cuckoo wasps parasitizing COLw hosts are more similar to those of their hosts than to those of conspecific males
We hypothesized that only cuckoo wasp females have been selected to mimic the CHC profile of their female hosts. In contrast, cuckoo wasp males should not have been selected to mimic the CHC profile of their female hosts due to the lack of interaction with their host or their host’s nest. This pattern should be particularly recognizable in COLw and their brood parasites, because COLw are free from selection to synthesize alkene-enriched CHC profiles due to the abandonment of prey embalming behavior . As predicted, we found that cuckoo wasp females parasitizing COLw species (with the notable exception of H. chalybaeum) largely synthesize the same type of compounds as their female COLw hosts (especially the same homologous series of methyl-branched alkanes). In contrast, male cuckoo wasps, irrespective of what type of host they develop from (i.e., COLw or HYMw), show alkene rich CHC profiles that are compositionally similar to the CHC profiles of HYMw. Interestingly, a phylogenetic study of the host species gives evidence that HYMw species are ancestral to COLw . Thus, males of all cuckoo wasps seem to still synthesize a more ancestral CHC profile and are evolutionarily less affected by the shift towards higher quantities of methyl-branched hydrocarbons in COLw than are cuckoo wasp females.
Brood parasitism can exert strong selection pressure on the host to improve its ability to detect brood parasites. The brood parasites, in turn, are counter-selected for traits that enable them to remain undetected by their host. The switch from using Hymenoptera to using Coleoptera as prey in Philanthidae liberated COLw from producing CHC profiles enriched with unsaturated CHCs, a requirement for embalming and preserving Hymenoptera prey. We demonstrated that this relaxation likely allowed COLw species to evolve distinct and species-specific CHC profiles dominated by methyl-branched alkanes as a strategy to escape chemical mimicry by their Hedychrum brood parasites. Several lines of evidence are consistent with our conclusions: (1) COLw CHC profiles are conspicuously distinct from the alkene-enriched CHC profiles of HYMw necessary for prey preservation and very distinct between each species, with species-specific types of methyl-branched alkanes; (2) the CHC profile overlap between COLw and their brood parasites is smaller than between HYMw and their brood parasites, although the Bray–Curtis dissimilarities of the CHC profiles were not significantly smaller in HYMw compared to COLw; (3) female cuckoo wasps show a CHC profile that resembles that of their COLw host species (except for H. chalybaeum), whereas the CHC profile of their male conspecifics resembles the host’s likely ancestral CHC profile of HYMw. In addition, COLw species show CHC profiles with longer chain lengths than HYMw and all Hedychrum species, which may constitute an additional strategy to escape chemical mimicry by their brood parasites. In summary, the observed patterns suggest that a diversification of CHC profiles with methyl-branched alkanes has evolved in COLw host species and that their cuckoo wasps “followed” their hosts in chemical space, consistent with an arms race hypothesis. Furthermore, our data suggest that HYMw may counteract the chemical mimicry of cuckoo wasps by exhibiting a larger intraspecific CHC profile variability than COLw. Although we provide evidence for a coevolutionary arms race on the phenotypic level, future studies should consider searching for footprints of this arms race also at the genetic level.
Collection and origin of the insect samples
We studied the CHC profiles of 5–14 individual wasps per species and sex of twelve species (five species of the genus Hedychrum, five species of the genus Cerceris, and two species of the genus Philanthus, Table 2). To avoid considering individuals from populations that have experienced very different climatic conditions in the past, we collected all specimens in the Upper Rhein valley between Freiburg and Germersheim (150 km). These samples were collected over the course of multiple years, as some species are considered rare or difficult to find in larger numbers (e.g., Hedychrum chalybaeum). As sampling over multiple years could potentially increase variation (compared to sampling in just one year), we intentionally also collected the more common species over the course of multiple several years. We screened for possible differences in CHC profile composition between samples collected at different sample locations and in different years but did not find any significant differences (data not shown). Besides the five target species of Hedychrum (see above), we included an additional species of this genus (i.e., H. longicolle collected in Almarail, Spain in 2011) in the phylogenetic analysis for reasons unrelated to the present study. All samples were collected with an insect net and were identified by ON using a reference collection and DNA barcoding. Each wasp was placed in a glass vial (1.5 mL) and transported to the lab, where the specimen was killed by freezing. All specimens were stored at − 20 °C until extracting their CHCs and (only in a subset of samples) their DNA. Note that CHCs were extracted prior to any molecular work. Voucher specimens of samples whose DNA we studied are deposited at the Zoological Research Museum Alexander Koenig in Bonn or in the personal collections of ON and TS (Additional file 2: Table S1). Species of Hedychrum are known to be very host specific, and the currently known host associations of the studied species is summarized in Table 3.
Frozen insects were allowed to thaw for about five minutes and were subsequently submerged in pure n-hexane (Merck, Darmstadt, Germany) to extract the CHCs. After 10 min, the CHC extract was transferred into another glass vial, concentrated with a gentle stream of CO2 until approximately 80–100 µL of the solvent remained. The CHC extract was stored at − 20 °C, and the insect was stored separately in 100% ethanol to preserve its DNA.
A HP 6890 gas chromatograph (GC) coupled with a HP 5973 Mass Selective Detector (MS) (Hewlett Packard, Waldbronn, Germany) and an Agilent 7890/5975 GC/MS System were used for analyzing the extracts. The GCs (split/splitless injector in splitless mode for 1 min, injected volume: 1 µl at 300 °C injector temperature) were equipped with a DB-5 Fused Silica capillary column (30 m × 0.25 mm ID, df = 0.25 µm2, J&W Scientific, Folsom, USA). Helium was used as carrier gas with a constant flow of 1 mL/min. We applied the following temperature program: start temperature at 60 °C, with an increase of 5 °C/min until 300 °C, and isotherm at 300 °C for 10 min. An ionization voltage of 70 eV (source temperature: 230 °C) was set for the acquisition of the mass spectra by electron ionization. Ion masses of 40–600 units were recorded.
In order to identify the double-bond position of alkenes, 1–5 extracts of each sex and species (depending on the available amount of CHC extract) were pooled and used for dimethyl disulfide (DMDS) derivatization following the protocol provided by Carlson and colleagues . The double bond positions of alkadienes remained unidentified due to the low quantity of this substance class in the CHC extracts. Alkadienes were characterized according to their retention indices.
We analyzed the cuticular hydrocarbon composition of cuckoo wasps and their hosts with the software AMDIS version 2.71 (Automated Mass Spectral Deconvolution and Identification System, http://chemdata.nist.gov/mass-spc/amdis/). AMDIS uses mass spectra similarities and retention indices to select target compounds. Specifically, we first created a mass spectral library (which contains more than 600 identified mass spectra of common hydrocarbons and their retention indices). We used the protocol by Carlson and colleagues , in which the elution patterns of CHC are described to confirm the identified methyl-branched alkanes in this library. The parameters used in AMDIS were as follows: component width = 22, adjacent peak subtraction = 2, resolution = medium, sensitivity = low, and shape requirements = medium.
After CHC had been identified, we calculated their relative abundance by dividing the total ion count of each peak relative to the total ions count of CHC in a profile occurring in the range of C21–C33. To ensure that a compound occurred in the majority of samples of a given species and did not represent an artifact of CHC extract concentration differences or of the sensitivity of the GC/MS, we set a threshold for the consideration of any compound within a group. This minimum threshold requirement was met if the compound occurred in at least 50% of all specimens (per sex and species) and if the mean relative quantitative abundance was at least 0.1%.
We estimated the total number of compounds and the number of compounds per compound class (i.e., linear alkanes, unsaturated compounds, and methyl-branched alkanes). Additionally, we calculated the mean chain length of the CHCs in a given CHC extract by summing up the relative amount of each peak within the range C21 to C33, weighted by its retention index. This value indicates the retention index at which half of the relative amount of the CHC profile occurs.
Statistical analyses of chemical data
CHC profiles were compared and differences were visualized using multivariate methods using the software R version 3.02 . We conducted Non-metric Multidimensional Scaling (NMDS) using Bray–Curtis dissimilarity to visualize CHC profile similarity in two-dimensional graphs [54, 55]. All inferred stress values fell below 0.15. We assessed the degree of similarity and overlap between CHC profiles of the different groups using an analysis of similarity (ANOSIM; , a non-parametrical test that operates on a ranked dissimilarity matrix. The obtained test statistic R (− 1 < R < 1) indicates the degree of similarity between and within groups. An R value close to 1 indicates complete separation between the tested groups and values close to zero indicates more similarity between the groups (greater overlap, less separation). Negative values of R are uncommon and do not have a biological meaning.
We used Welch corrected t-tests [57, 58] to compare various traits between COLw and HYMw and their brood parasites. Specifically, we compared Bray–Curtis dissimilarities between hosts and their brood parasites, the proportion, number, and diversity of CHC compounds, and the mean chain length. P values from the ANOSIM tests and the t-tests were adjusted using the Holm-Bonferroni correction  to account for multiple testing.
For most of the analyses, we used the R package vegan . R packages used for plotting were ade4 , ape , flagme , phytools , and xcms .
We used muscle tissue from one or two individuals per species for DNA extraction (see below) after CHC extraction.
We used a set of degenerated oligonucleotide primer pairs  to amplify twelve single-copy protein-coding nuclear genes in cuckoo wasps (same as studied by ; Additional file 2: Table S2). Polymerase chain reactions (PCR) and bidirectional direct Sanger sequencing followed the protocol given by .
Forward and reverse DNA strands were assembled to contigs and trimmed (to exclude binding sites of PCR primers) in Geneious (version 6.1; . All contigs were aligned with the LINS-i algorithm of MAFFT version 7.123 . Intronic and exonic regions were annotated manually by aligning a reference DNA sequence of each target gene from the 1KITE transcriptome of the cuckoo wasp Chrysis terminata and identifying canonical splice sites in the genomic DNA sequencing (i.e., the dinucleotide pair GC-AG) (see . We manually removed uninformative and ambiguously aligned sites from the intronic DNA sequences. Finally, we concatenated all exons and introns to a supermatrix and defined three partitions: (1) 1st and 2nd codon position of exons, (2) 3rd codon position of exons, (3) introns.
We inferred the phylogenetic relationships of the studied Hedychrum species by applying the maximum likelihood optimality criterion. Specifically, we first searched for the best fitting substitution model for each partition using Modelfinder  implemented in IQ-TREE (version 1.5.5; ) and applying the corrected Akaike information criterion (AICc) to choose between models. The partition-specific substitution models are listed in Additional file 2: Table S3. The selected substitution models and the empirically inferred substitution parameters were subsequently used to infer a phylogenetic tree within IQ-TREE. Statistical branch support was estimated from 1000 non-parametric bootstrap replicates. We additionally conducted a phylogenetic analysis in a Bayesian framework with the software MrBayes version 3.1.2 [72, 73]. For this purpose, we repeated the model selection step in Modelfinder, but restricting the tested models to those implemented in MrBayes. We started two parallel runs, both with a random starting tree, over 107 generations. We sampled trees every 105 generation and discarded the first 106 generations of both runs as burn-in. We used the remaining trees from all runs in order to calculate a 50% majority rule consensus tree. Convergence was assessed with the software Tracer (version 1.6; ).
To infer the phylogenetic relationship of Cerceris and Philanthus species, we use the phylogeny published by Wurdack and colleagues  and trimmed it to the species analyzed in the present study.
Availability of data and materials
Voucher specimens are deposited in the repositories listed in Additional file 2: Table S1. Accession numbers for the genes used for the molecular phylogenies are listed in Additional file 2: Table S6. Relative compositions of the CHC compounds of all specimens are listed in Additional file 2: Table S7.
Gas chromatography/mass spectrometry
Laine A-L. Role of coevolution in generating biological diversity: spatially divergent selection trajectories. J Exp Bot. 2009;60:2957–70.
Thompson JN. The coevolutionary process. Chicago: University of Chicago Press; 1994.
Thompson JN. The coevolving web of life (American society of naturalists presidential address). Am Nat. 2009;173:125–40.
Thompson JN. The role of coevolution. Science. 2012;335:410–1.
Meyer JR, Dobias DT, Weitz JS, Barrick JE, Quick RT, Lenski RE. Repeatability and contingency in the evolution of a key innovation in phage lambda. Science. 2012;335:428–32.
Rothstein SI. A model system for coevolution: avian brood parasitism. Ann Rev Ecol Syst. 1990;21:481–508.
Davies NB, Bourke AF, Brooke MD. Cuckoos and parasitic ants: interspecific brood parasitism as an evolutionary arms race. TREE. 1989;4:274–8.
Winfree R. Cuckoos, cowbirds and the persistence of brood parasitism. TREE. 1999;14:338–43.
Spottiswoode CN, Kilner RM, Davies NB. Brood parasitism. In: Royle NJ, Smiseth PT, Kölliker M, editors. The evolution of parental care. Oxford: Oxford University Press; 2012. p. 226–356.
Spottiswoode CN, Stevens M. How to evade a coevolving brood parasite: egg discrimination versus egg variability as host defences. Proc Roy Soc B. 2011;278:3566–73.
van Valen L. A new evolutionary law. Evol Theory. 1973;1:1–30.
Roitberg BD. Chemical communication. In: Córdoba-Aguilar A, González-Tokman D, González-Santoyo I, editors. Insect behavior: from mechanisms to ecological and evolutionary consequences. Oxford: Oxford University Press; 2018. p. 145–57.
Blomquist GJ, Bagnères A-G. Insect hydrocarbons: biology, biochemistry, and chemical ecology. Cambridge: Cambridge University Press; 2010.
Lenoir A, Dettorre P, Errard C, Hefetz A. Chemical ecology and social parasitism in ants. Ann Rev Entomol. 2001;46:573–99.
Lambardi D, Dani FR, Turillazzi S, Boomsma JJ. Chemical mimicry in an incipient leaf-cutting ant social parasite. Behav Ecol Sociobiol. 2007;61:843–51.
Kroiss J, Schmitt T, Strohm E. Low level of cuticular hydrocarbons in a parasitoid of a solitary digger wasp and its potential for concealment. Entomol Sci. 2009;12:9–16.
Polidori C, Geyer M, Schmitt T. Do Sphecodes cuckoo bees use chemical insignificance to invade the nests of their social Lasioglossum bee hosts? Apidologie. 2020;51:147–62.
Dettorre P, Mondy N, Lenoir A, Errard C. Blending in with the crowd: Social parasites integrate into their host colonies using a flexible chemical signature. Proc Roy Soc B. 2002;2002(269):1911–8.
von Beeren C, Schulz S, Hashim R, Witte V. Acquisition of chemical recognition cues facilitates integration into ant societies. BMC Ecol. 2011;11:30.
Akino T, Knapp JJ, Thomas JA, Elmes GW. Chemical mimicry and host specificity in the butterfly Maculinea rebeli, a social parasite of Myrmica ant colonies. Proc Roy Soc B. 1999;266:1419–26.
Wurdack M, Herbertz S, Dowling D, Kroiss J, Strohm E, Baur H, Niehuis O, Schmitt T. Striking cuticular hydrocarbon dimorphism in the mason wasp Odynerus spinipes and its possible evolutionary cause (Hymenoptera: Chrysididae, Vespidae). Proc Roy Soc B. 2015;282:20151777.
Jongepier E, Foitzik S. Ant recognition cue diversity is higher in the presence of slavemaker ants. Behav Ecol. 2016;27:304–11.
Brandt M, Heinze J, Schmitt T, Foitzik S. A chemical level in the coevolutionary arms race between an ant social parasite and its hosts. J Evol Biol. 2005;18:576–86.
Errard C, Ruano F, Richard F-J, Lenoir A, Tinaut A, Hefetz A. Co-evolution-driven cuticular hydrocarbon variation between the slave-making ant Rossomyrmex minuchae and its host Proformica longiseta (Hymenoptera: Formicidae). Chemoecology. 2006;16:235–40.
Guillem RM, Drijfhout F, Martin SJ. Chemical deception among ant social parasites. Curr Zool. 2014;60:62–75.
Martin SJ, Carruthers JM, Williams PH, Drijfhout FP. Host specific social parasites (Psithyrus) indicate chemical recognition system in bumblebees. J Chem Ecol. 2010;36:855–63.
Lorenzi MC. The result of an arms race: the chemical strategies of Polistes social parasites. Ann Zool Fenn. 2006;43:550–63.
Strohm E, Kroiss J, Herzner G, Laurien-Kehnen C, Boland W, Schreier P, Schmitt T. A cuckoo in wolves ́clothing? Chemical mimicry in a specialized cuckoo wasp of the European beewolf (Hymenoptera, Chrysididae and Crabronidae). Front Zool. 2008;5:2.
Emery C. Über den Ursprung der dulotischen, parasitischen und myrmekophilen Ameisen. Biol Zent Bl. 1909;29:352–62.
Peters RS, Krogmann L, Mayer C, Donath A, Gunkel S, Meusemann K, Kozlov A, Podsiadlowski L, Petersen M, Lanfear R, Diez PA, Heraty J, Kjer KM, Klopfstein S, Meier R, Polidori C, Schmitt T, Liu S, Zhou X, Wappler T, Rust J, Misof B, Niehuis O. Evolutionary history of the Hymenoptera. Curr Biol. 2017;27:1013–8.
Sann M, Niehuis O, Peters RS, Mayer C, Kozlov A, Podsiadlowski L, Bank S, Meusemann K, Misof B, Bleidorn C, Ohl M. Phylogenomic analysis of Apoidea sheds new light on the sister group of bees. BMC Evol Biol. 2018;18:1–15.
Sann M, Meusemann K, Niehuis O, Escalona HE, Mokrousov M, Ohl M, Pauli T, Schmid-Egger C. Reanalysis of the apoid wasp phylogeny with additional taxa and sequence data confirms the placement of Ammoplanidae as sister to bees. Syst Entomol. 2021;46:558–69.
Pauli T, Meusemann K, Kukowka S, Sann M, Donath A, Mayer C, Oeyen JP, Ballesteros Y, Berg A, van den Berghe E, Escalona HE, Guglielmino A, Niehuis M, Olmi M, Podsiadlowski L, Polidori C, de Rond J, Rosa P, Schmitt T, Strumia F, Wurdack M, Liu S, Zhou X, Misof B, Peters RS, Niehuis O. Analysis of RNA-seq, DNA target enrichment, and Sanger nucleotide sequence data resolves deep splits in the phylogeny of cuckoo wasps (Hymenoptera: Chrysididae). Insect Syst Div. 2021;5:1–14.
Strohm E, Linsenmair KE. Females of the European beewolf preserve their honeybee prey against competing fungi. Ecol Entomol. 2001;26:198–203.
Herzner G, Kaltenpoth M, Poettinger T, Weiss K, Koedam D, Kroiss J, Strohm E. Morphology, chemistry and function of the postpharyngeal gland in the South American digger wasps Trachypus boharti and Trachypus elongatus. PLoS ONE. 2013;8: e82780.
Weiss K, Strohm E, Kaltenpoth M, Herzner G. Comparative morphology of the postpharyngeal gland in the Philanthinae (Hymenoptera, Crabronidae) and the evolution of an antimicrobial brood protection mechanism. BMC Evol Biol. 2015;15:1.
Wurdack M, Polidori C, Keller A, Feldhaar H, Schmitt T. Release from prey preservation behavior via prey switch allowed diversification of cuticular hydrocarbon profiles in digger wasps. Evolution. 2017;71:2562–71.
Herzner G, Strohm E. Fighting fungi with physics: food wrapping by a solitary wasp prevents water condensation. Curr Biol. 2007;17:R46-47.
Herzner G, Strohm E. Food wrapping by females of the European Beewolf, Philanthus triangulum, retards water loss of larval provisions. Physiol Entomol. 2008;33:101–9.
Strohm E, Kaltenpoth M, Herzner G. Is the postpharyngeal gland of a solitary digger wasp homologous to ants? Evidence from chemistry and physiology. Insectes Soc. 2010;57:285–91.
Bagnères AG, Morgan ED. The postpharyngeal glands and the cuticle of Formicidae contain the same characteristic hydrocarbons. Experientia. 1991;47:106–11.
Walsh J, Pontieri L, d’Ettorre P, Linksvayer TA. Ant cuticular hydrocarbons are heritable and associated with variation in colony productivity. Proc Roy Soc B. 2020;287:20201029.
Martin SJ, Helanterä H, Drijfhout FP. Is parasite pressure a driver of chemical cue diversity in ants? Proc Roy Soc B. 2011;278:496–503.
Lorenzi MC, Azzani L, Bagnéres A-G. Evolutionary consequences of deception: complexity and informational content of colony signature are favored by social parasitism. Curr Zool. 2014;60:137–48.
Gibbs AG, Chippindale AK, Rose MI. Physiological mechanisms of evolved desiccation resistance in Drosophila melanogaster. J Exp Biol. 1997;200:1821–32.
Gibbs A, Pomonis JG. Physical properties of insect cuticular hydrocarbons: the effects of chain length, methyl-branching and unsaturation. Comp Biochem Physiol B. 1995;112:243–9.
Rouault J-D, Marican C, Wicker-Thomas C, Jallon J-M. Relations between cuticular hydrocarbon (HC) polymorphism, resistance against desiccation and breeding temperature; a model for HC evolution in D. melanogaster and D. simulans. Genetica. 2004;120:195–212.
Blomquist GJ. Biosynthesis of cuticular hydrocarbons. In: Blomquist GJ, Bagnères A-G, editors. Insect hydrocarbons: biology, biochemistry, and chemical ecology. Cambridge: Cambridge University Press; 2010. p. 35–52.
Øien IJ, Moksnes A, Røskaft E. Evolution of variation in egg color and marking pattern in European passerines: adaptations in a coevolutionary arms race with the cuckoo, Cuculus canorus. Behav Ecol. 1995;6:166–74.
Kroiss J, Bordon S, Strohm E. Hydrocarbons in the nest material of a solitary digger wasp represent a kairomone for a specialized cuckoo wasp. Anim Behav. 2008;76:1555–63.
Carlson DA, Roan CS, Yost RA, Hector J. Dimethyl disulfide derivatives of long chain alkenes, alkadienes, and alkatrienes for gas chromatography/mass spectrometry. Anal Chem. 1989;61:1564–71.
Carlson DA, Bernier UR, Sutton BD. Elution patterns from capillary GC for methyl-branched alkanes. J Chem Ecol. 1998;24:1845–65.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria, 2013. URL http://www.R-project.org/.
Kruskal JB. Multidimensional scaling by optimizing goodness of fit to a nonmetric hypothesis. Psychometrika. 1964;29:1–27.
Kruskal JB. Nonmetric multidimensional scaling: a numerical method. Psychometrika. 1964;29:115–29.
Clarke KR. Non-parametric multivariate analyses of changes in community structure. Aust J Ecol. 1993;18:117–43.
Welch BL. The significance of the difference between two means when the population variances are unequal. Biometrika. 1938;29:350–62.
Ruxton GD. The unequal variance t-test is an underused alternative to Student’s t-test and the Mann-Whitney U test. Behav Ecol. 2006;17:688–90.
Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6:65–70.
Oksanen J, Blanchet FG & others. Vegan: Community Ecology Package. R-package version 2.0-10, 2013.
Dray S, Dufour A-B. The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007;22:1–20.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.
Robinson M. Flagme: Fragment-level analysis of GCMS-based metabolomics data. R package version 1.14, 2013.
Revell LJ. Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–23.
Smith CA, Want EJ, O’Maille G, Abagyan R, Siuzdak G. XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem. 2006;78:779–87.
Hartig G, Peters RS, Borner J, Etzbauer C, Misof B, Niehuis O. Oligonucleotide primers for targeted amplification of single-copy nuclear genes in apocritan Hymenoptera. PLoS ONE. 2012;7: e39826.
Pauli T, Castillo-Cajas RF, Rosa P, Kukowka S, Berg A, van den Berghe E, Fornoff F, Hopfenmüller S, Niehuis M, Peters RS, Staab M, Strumia F, Tischendorf S, Schmitt T, Niehuis O. Phylogenetic analysis of cuckoo wasps (Hymenoptera: Chrysididae) reveals a partially artificial classification at the genus level and a species-rich clade of bee parasitoids. Syst Entomol. 2019;44:322–35.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Meintjes P, Drummond A. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Kalyaanamoorthy S, Minh BQ, Wong TK, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9.
Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2014;32:268–74.
Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754–5.
Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6. Available from http://beast.bio.ed.ac.uk/Tracer, 2014.
Alfken JD. Verzeichnis der Goldwespen (Chrysiden) Nordwestdeutschlands. Abh Naturwiss Verein Bremen. 1915;23:291–5.
Berland L, Bernard F. Hyménoptères vespiformes III. (Cleptidae, Chrysidae, Trigonalidae). Faune de France, vol. 34. Paris: Paul Lechevalier; 1938.
Simon-Thomas RT, Simon-Thomas AMJ. Some observations on the behavior of females of Philanthus triangulum (F.) (Hymenoptera, Sphecidae). Tijdschr Entomol. 1972;115:123–39.
Mingo Pérez E. Las especies españolas de Hedychrum LATR., 1806 (Hym., Chrysididae). Eos. 1981;55:143–54.
Morgan D. Cuckoo-Wasps. (Hymenoptera, Chrysididae). Handbooks for the Identification of British Insects, vol. 6. St. Albans: Royal Entomological Society; 1984.
Veenendaal RL. Het verborgen ei van Hedychrum rutilans (Hymenoptera: Chrysididae). Entomolog Ber. 1987;47:169–71.
Bitsch J, Barbier Y, Gayubo SF, Schmidt K, Ohl M. Hymenopteres Sphecidae d’Europe Occidentale (Vol 2). Faune de France, 82. Paris: Fédération Française des Sociétés de Science Naturel; 1997.
Saure C. Beobachtungen und Anmerkungen zur Wirtsbindung einiger Goldwespenarten im nordostdeutschen Raum (Hymenoptera: Chrysididae: Chrysidinae). Bembix. 1998;10:15–8.
van der Smissen J. Die Wildbienen und Wespen Schleswig-Holsteins—Rote Liste. Band I-III. Flintbeck: Landesamt für Natur und Umwelt des Landes Schleswig-Holstein; 2001.
Reder G, Burger R. Nachweise der Goldwespe Hedychrum chalybaeum (Dahlbom, 1854) in Rheinland-Pfalz (Hymenoptera: Chrysididae). Fauna Flora Rheinl-Pfalz. 2009;11:851–6.
Paukkunen J, Berg A, Soon V, Ødegaard F, Rosa P. An illustrated key to the cuckoo wasps (Hymenoptera, Chrysididae) of the Nordic and Baltic countries, with description of a new species. ZooKeys. 2015;548:1–116.
Wiesbauer H, Rosa P, Zettel H. Die Goldwespen Mitteleuropas—Biologie, Lebensräume, Artenporträts. Stuttgart: Eugen Ulmer Verlag; 2020.
Open Access funding enabled and organized by Projekt DEAL. This publication was supported by the Open Access Publication Fund of the University of Würzburg. O.N. and T.S. acknowledge the German Research Foundation (DFG) for supporting this study (NI1387/1-1; NI1387/5-1; SCHM 2645/1-1; SCHM 2645/2-1). R.C. acknowledges the German Academic Exchange Service (DAAD) for a scholarship to finish her PhD. The study was also supported by the Ministerio de Economía, Industria y Competitividad (España) in the project CGL2017- 83046-P (to C.P. and T.S.).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
Boxplots of the number of a) total CHC compounds and b) methyl-branched alkanes and the c) proportion of methyl-branched alkanes in the profiles of Hymenoptera preying (HYMw) and Coleoptera preying (COLw) apoid wasps, hosts to Hedychrum species in this study. Females are depicted in dark grey and males in light gray colors. Figure S2. Intra- and interspecific variability of cuticular hydrocarbon profiles in a) female and b) male individuals of all host species. Bray–Curtis dissimilarities were calculated between all individuals of species hunting the same type of prey (“between species”) and between individuals of the same species (species name indicated). In this case, however, in comparison to the figure presented in the text, only a maximum of five individuals were randomly selected for each group so that the number of specimens used in each group remains the same (5). Figure S3. Bray–Curtis dissimilarities between two randomly selected CHC profiles of samples collected from the same population of a species belonging to HYMw or to COLw in females (a) and males (b). Pairwise distances between two individuals coming from the same species were randomly selected in 1000 simulations in each case.
Additional file 2: Table S1.
Provenance and details of collection and storage of all specimens used in the analyses. Table S2. Characteristics of the primers used and number of nucleotides in the exonic regions of the nuclear genes used for the phylogenetic analyses. Numbers before the semicolon indicate number of exonic regions, and the numbers following the semicolon indicate the total number of nucleotides used per gene. Table S3. Substitution models chosen by Modelfinder for the respective phylogenetic analysis in IQTree and MrBayes. Table S4. Mean relative abundance ± standard deviation for each peak or mixed peak included in the NMDS of the females of all species analyzed. Table S5. Mean relative abundance ± standard deviation for each peak or mixed peak included in the NMDS of the males of all species analyzed. Table S6. Accession numbers for the genes used for the molecular phylogeny in Fig. 1. Table S7. Relative composition of the CHC compounds of all specimens. First row shows retention index with compound name abbreviation. First column shows the specimens name, sex, and number
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Castillo, R., Wurdack, M., Pauli, T. et al. Evidence for a chemical arms race between cuckoo wasps of the genus Hedychrum and their distantly related host apoid wasps. BMC Ecol Evo 22, 138 (2022). https://doi.org/10.1186/s12862-022-02093-8
- Chemical mimicry
- Cuticular hydrocarbons
- Evolutionary arms race