Behavioural and neurogenomic responses of host workers to social parasite invasion in a social insect

The strong coevolutionary arms race between social parasites and their hosts has dramatically shaped the life-history traits of both parties. One of the main strategies exhibited by hosts in response to parasitism is reproduction by host workers. We lack a mechanistic understanding of how these defence strategies unfold and, specifically, whether hosts exhibit more subtle strategies to reduce the costs of parasitism from the outset. Here we test the hypothesis that there are both behavioural and neurogenomic signatures of worker responses to parasitism, prior to overt expression in the form of egg-laying; we test this using the social parasite—social host system of the paper wasps Polistes sulcifer-Polistes dominula. We characterized individual workers’ position within the social interaction network of queenright and host colonies immediately after parasite usurpation, weeks before the workers’ reproductive rebellion is evident. Parasitism influenced network centrality measures, with workers in parasitized colonies showing increased connectedness and centrality compared to those in unparasitized ones. Next, we quantified brain gene expression levels for five genes related to physiological and behavioural phenotypes in Polistes wasps. The gene Imaginal disc growth factor (Idgf4), thought to be responsive to changes in the social environment, was significantly down-regulated in workers from parasitized colonies; this may be an indication that parasitized workers are anticipating a shift toward a less worker-like phenotype in preparation for their reproductive rebellion. Our results provide the first evidence of early behavioural and neurogenomic responses of host workers toward the presence of an inquiline social parasite in a social insect.


Introduction
For social insect colonies, the threat of invasion by social parasites is part of life. From beetles to butterflies, many species have evolved the ability to exploit the socially acquired resources of insect colonies, e.g. their brood care, nest building and food supplies (Wilson 1971). Many of the social parasites that specialise in exploiting the social Hymenoptera (bees, wasps and ants) evolved from a social ancestor (Lowe et al. 2002): social parasites are widespread in ants (Hölldobler and Wilson 1990), allodapine bees (Smith et al. 2013), honeybees and bumblebees (Lhomme and Hines 2018) and wasps (Cervo 2006;Cini et al. 2019b). As parasites strive (in evolutionary terms) to be better exploiters of their social hosts, hosts respond by evolving strategies to reduce the chance and costs of invasion, leading to a coevolutionary arms race between the two species, and dramatically shaping their morphologies, behaviours and other life-history traits (Dawkins and Krebs 1979). Social parasite and host systems, therefore, present opportunities to study the ultimate and proximate bases of co-evolution.
Several lines of defence against the social parasite can be used by host species (Grüter et al. 2018). A first line of host defence is to avoid being conquered, for example through increased recognition or defensive strategies (Brandt et al. 2005;Ortolani and Cervo 2010): host species might evolve strategies that make it more difficult for the social parasite to integrate into the colony, for instance by evolving more complex chemical colony signatures which are more difficult to be mimicked by the social parasites (e.g. Bruschini et al. 2010;Lorenzi et al. 2014;Martin et al. 2011). In the event of an invasion, a second line of defence kicks in: host workers often launch a rebellion. Rarely can hosts expel their parasite; instead, they have evolved strategies to reduce the fitness of the parasite. Indeed, incomplete control by the social parasite over host worker reproduction has been highlighted in a few species of wasps and bumblebees Reed and Akre 1983) and two main strategies have been reported. Workers may show an overt response, as seen in Temnothorax ants, where workers enslaved by the social parasite Protomognathus americanus kill up to two thirds of the parasite's pupae (Achenbach and Foitzik 2009). On the other hand, another common response is more covert, through worker reproduction . From the workers' point of view this strategy is clearly beneficial, as they can increase their direct fitness through the production of both males and females (as workers can mate and produce reproductive females (Strassmann et al. 2004)). By contrast, from the parasite's perspective, investment by workers in reproduction should be minimized in order to maximize worker investment in rearing the parasite's brood . Indeed, reproductive workers might be costly as they might contribute less to colony efficiency (e.g. Heinze 2008; Monnin et al. 2003).
However, both these strategies are 'last resorts' and only kick in several weeks after the parasite first invades. It seems remarkable that there has not been selection for host colonies to respond more immediately to the invasion of a parasite, as an early form of 'damage limitation'. Signatures of early response may be apparent through more subtle changes in behaviour and also in the expression of genes in the brain which are the front line response to changes in the environment (Robinson et al. 2005;Rubenstein and Hofmann 2015).
Here we explore whether there are signatures of early response to parasitism through subtle changes in behaviour and/or shifts in brain gene expression in a host paper wasp and its congeneric social parasite. Polistes sulcifer is an obligate inquiline (workerless) social parasite that invades nests of the social species Polistes dominula through violent fights Turillazzi et al. 1990). The parasite replaces the host queen both behaviourally and chemically, achieving the top social rank in the dominance hierarchy (Cervo 2006;Dapporto et al. 2004). After usurpation, parasitized colonies produce parasite offspring, but also a small number of host offspring (Cervo 2006;Dapporto et al. 2004). Host workers are thus largely fooled (or suppressed) by the social parasite: they remain on the nest and rear parasite offspring, preferentially directing their provisioning effort toward the parasite brood as opposed to host brood . However, recent results show that workers of parasitized colonies show signs of reproductive rebellion: they have more developed ovaries than workers from unparasitized colonies, and they lay eggs more rapidly when given the opportunity ). However, this response is not evident until 6 weeks after parasitism , despite the fact that workers are able to develop their ovaries in only six days (Monnin et al. 2009) (Fig. 1a). The most parsimonious explanation, therefore, is that on usurpation host workers cannot detect the presence of the parasite, or if they can, the parasite is able to effectively suppress any ovary activation in host workers; but, after a considerable period of time (several weeks), the workers are able to overcome these issues and develop their ovaries. An outstanding question is whether workers show subtle signs of response to the parasite from the outset of usurpation, suggesting that perhaps they can detect and/or are imperfectly suppressed by the parasite; or whether workers are indeed effectively deceived/suppressed on usurpation, and remain so until the parasite is well established. Subtle signs of response to the parasite may be apparent in their behavioural social networks within the colony, and/or changes in brain gene expression (neurogenomic signatures), as the front-line response to changes in the social environment.
We tested these two hypotheses by comparing worker behaviour and brain gene expression in parasitized and non-parasitized P. dominula colonies across different stages of parasitism ( Fig. 1). First, we characterized colony-level dominance networks in the early phase after parasite usurpation, within 2 weeks of usurpation, well before ovarian development is apparent  to determine whether there are any subtle changes in worker social behaviour as an early response to parasitism (Prediction 1). Dominance interactions are crucial to the regulation of worker activities in primitive insect societies (Naug 2008) and dominance interaction networks have been shown to reflect changes in social structure following major events such as changes in colony composition [e.g. queen removal, (Bhadra and Jordán 2013;Bhadra et al. 2009)]. Moreover, changes in individual network properties correlate with changes in behavioural phenotype in response to colony-level changes in social structure (Bhadra and Jordán 2013). We characterized the network positions of individual host workers using centrality measures that characterized both the local (first order) and "extended" (second to n-order) networks for each individual (Pinter-Wollman et al. 2014;Wey et al. 2008); we then compared social network features in workers from parasitized and unparasitized colonies. If parasitism has an effect on workers' social behaviour soon after usurpation, we expect to find differences in the dominance network metrics of workers from parasitized and unparasitized colonies. Specifically, as dominance behaviour and reproductive skew are positively correlated in this species (Dapporto et al. 2010;Reeve et al. 1991), we predict workers from parasitized colonies to have higher centrality (i.e. engaging in more dominance interactions) than in unparasitized colonies (Fig. 1b).
To determine whether workers show any neurogenomic changes in response to the early stages of parasitism, we quantified brain gene expression for a set of five candidate genes in workers from parasitized and unparasitized colonies. We chose genes known to indicate response to changes in the social environment and reproductive status in Polistes Ferreira et al. 2013;Manfredini et al. 2018;Sumner et al. 2006;Toth et al. 2014;Toth et al. 2007) (see Table 2 for details). If workers detect a change in their colony (putatively, the presence of the parasite herself) soon after usurpation, we expect the gene expression of workers in parasitized and unparasitized colonies to differ (Prediction 2). In particular, we expect to find some downregulation of worker-biased genes and upregulation of queen-biased genes in parasitized workers (Table 2), indicating an attempt to instigate a brain transcriptomic shift toward a reproductive phenotype (Fig. 1c).

Prediction 1: Positions of workers differ in social networks of parasitized and unparasitized colonies
In Spring 2012, eight pre-emergence multiple foundresses colonies (range 2-4 foundresses) of P. dominula were collected in the vicinity of Florence (Tuscany, Italy). Each nest was fixed to the ceiling of a glass cage (15 cm × 15 cm × 15 cm) and maintained under controlled laboratory conditions with ad libitum sugar, water and paper material for nest building. All colony members and any newly emerged individual were individually marked with enamel paint dots (Testor Corporation) on the thorax and the wings before each observation. P. sulcifer females were collected in Spring from their overwintering sites at 2050 m altitude on Monti Sibillini (Umbria-Marche, Italy), transferred to the laboratory and kept in glass cages (15 cm × 15 cm × 15 cm) with ad libitum water and sugar, under controlled and constant low temperature inside a fridge, until it was time for their "activation" (temperature and exposure to direct natural and artificial sunlight following the protocol used by Ortolani et al. 2008). In the last week of May, when usurpation usually occurs in the field (Cervo and Turillazzi 1996;Ortolani et al. 2008;Turillazzi et al. 1990), four randomly chosen host colonies were subjected to parasitism by introducing a social parasite female inside the glass cage of a putative host nest (following the protocol in   (Fig. 1a). Parasitism was considered successful only when the parasite exhibited typical behaviours of a dominant social parasite (i.e. stroking the abdomen on the nest cell rims and dominating host individuals (Turillazzi et al. 1990)). Parasitized and unparasitized colonies were reared under the same conditions for two weeks, until the end of the experiment. Unparasitized and parasitized colonies did not differ in nest size (number of cells) nor in colony size (number of adult colony members at the beginning of the experiment) (Mann Whitney test, nest size U = 4, p = 0.314, n = 4; colony size: U = 7, p = 0.881, n = 4; Unparasitized colonies; mean ± st.dev: 54.80 ± 3.95 cells and 13.00 ± 2.94 adults; Parasitized colonies, mean ± st. dev: 63.53 ± 9.68 cells and 10.60 ± 2.22 adults). All colonies appeared healthy, with brood of all stages in the comb at the start of the experiment.
Each colony (n = 8) was videotaped for 1 h between 10 am-3 pm (when wasps are most active) on the day before usurpation, on day 1, day 7 and day 14 after usurpation (Fig. 1a); a total of 4 h of video observations were recorded for each of the four parasitized colonies and for each of the four unparasitized colonies.
All Ritualized Dominance Interactions (hereafter RDI) that occurred between any two individuals were recorded, noting the identities of the actor (i.e. the wasp that started the interaction) and the recipient. RDI includes domination acts, in which one wasp (the dominant individual) climbs on and antennates another wasp (the subordinate), and trophallaxis, i.e. liquid food exchange by mouth-to-mouth contact (Cini and Dapporto 2009;Pardi 1948) (Fig. 1b). Both interactions are directional, with dominant individuals (the actor) performing and not receiving RDI, and subordinate individuals receiving and not performing RDI. During trophallaxis the liquid is regurgitated by the subordinate individual toward the dominant one. To maintain the directionality of dominance, for trophallaxis the actor was considered the individual which received the liquid (i.e. the dominant individual) and the recipient the individual who gave the liquid (the subordinate one) [as in (Cini and Dapporto 2009)].
In animal social networks, usually each node represents an individual, and each link between two nodes ('edge') represents the existence and direction (from actor to recipient) of a dyadic interaction, and (in weighted networks) indicates the strength of the interaction (i.e. a value proportional to the frequency) by the thickness of the link. In the present study, we built directional weighted networks of parasitized and non-parasitized wasp colonies as video observations allowed us measuring the frequency of interactions, and RDI are directional interactions. Social networks for each colony at each time point (Fig. 1) were computed with UCINET 6 (v 6.665) (Borgatti et al. 2002).
For each individual and for each network we computed five standard node-centrality measures: weighted node degree, closeness centrality, eigenvector centrality, betweenness and clustering coefficient (Table 1). Tie direction is a critical aspect of a dominance network, as more dominant individuals deliver, but do not receive, RDIs, while the opposite is true for subordinates. For this reason, we separately computed Incoming (i.e. received) and Outgoing (i.e. performed) measures for the three directional network metrics weighted node degree, closeness centrality and eigenvector centrality. For example: outgoing weighted degree considers only the interactions started at a specific node and directed toward another node (e.g. by the actor of RDI); incoming closeness considers only the incoming ties (e.g. by the recipient of RDI).
Centrality measures are often positively correlated in social networks (Valente et al. 2008;Wey et al. 2008). For this reason, we applied a Principal Component Analysis (Varimax rotation) to reduce the number of variables into smaller number of uncorrelated principal components, extracting principal components with eigenvalues greater than 1, and then using these as social network metrics in our study (see results section and Supplementary material Table 1 The network metrics used in the study to describe worker position in the colonial social network Network metrics. Considered and their definition (Krause et al. 2009;Wey et al. 2008) Weighted node degree. Number of direct links (or 'edges') that a focal individual has with other individuals, weighted by the strength of each links (number of interactions) Closeness centrality. Average distance (measured as the number of edges) from a given starting node to all other nodes in the network; this measure captures how well connected an individual is to all the other individuals in the nest Eigenvector centrality. Assigns relative scores to all nodes in the network based on the concept that connections to high-scoring nodes (those with more edges) contribute more to the score of the node in question than equal connections to low-scoring nodes; this method captures the notion that a degree extends beyond just the first-order layer of connections that a focal node has Betweenness centrality. Measures how often a node appears on shortest paths between nodes in the network. For every pair of vertices in a connected graph, there exists at least one path between the vertices such that the sum of the weights of the edges (for weighted graphs) is minimized and the betweenness centrality for each node is the number of these shortest paths that pass through the node Clustering coefficient. Measures the density of a focal node neighborhood, i.e. the degree to which an individual's immediate neighbours are connected (Krause et al. 2009)

Table 2
The set of target genes investigated, with information about putative function and predictions Gene Putative gene function Evidence for putative function in social insects

Prediction
Our results

Arrestin (ARR )
Modulation of both visual and olfactory perception (Dolph et al., 1993;Merrill et al., 2002) Reproductive phenotype: linked to division of labour (upregulated in foragers) and to defence in bees (upregulated in defensive families) (Hunt, 2007;Hunt et al., 2007;Sasagawa et al., 2003) Social Environment: upregulated in social parasites during the phases of colony usurpation and integration  Downregulated in parasitized workers No significant differences, but significant interaction between parasitism and colony size. It suggests an effect of parasitism on ARR responsiveness to social environment Apolipophorin-III (ApoLp-III) General metabolic processes and lipid transport (Weers and Ryan, 2006) Reproductive phenotype: Apolipoprotein was down-regulated in adult queens vs. workers of the tropical paper wasp P. canadensis (Sumner et al., 2006) Social Environment: upregulated in social parasitic wasps during usurpation  Downregulated in parasitized workers No significant differences Heat Shock Protein 70 kDa (HSP70) Involved in response to environmental stress, in particular to thermal stress (King and MacRae, 2015) Associated with caste in P. metricus  and upregulated in queens compared to workers in the host species  Upregulated in parasitized workers Dataset3.csv). Comparisons made with each single network measure gave an overall similar pattern (data not shown).
We assessed the importance of parasitism on individual centrality in the dominance network running three nodelevel regression analyses (implemented in Ucinet, using 10,000 random permutations and 32,767 random number seeds), one for each centrality measure: outgoing centrality PC1, incoming centrality PC2, and clustering coefficient PC3. We used each centrality measure as dependent variable, and parasitism status and colony of origin as independent variables.

Prediction 2: Worker brain gene expression differs in parasitized and unparasitized colonies
Data examined here come from samples collected during a previous experiment, which were reared in a similar way to the above behavioural experiments . Wasp samples were collected from Central Italy in June 2005. Unparasitized colonies were collected from the field and parasitized colonies were obtained either by field collection or laboratory usurpations (as above). Colonies were brought to the lab and reared under controlled laboratory conditions (as above). For each colony, behavioural observations confirmed that parasites were dominant and successfully took over the nest . Workers were collected within the first two weeks of usurpation. Forty-two workers were sampled, of which 26 were from 14 unparasitized colonies and 18 were from 9 parasitized colonies. Two workers were analysed per colony, except for two colonies in which only one worker was analysed. Wasps were killed directly into RNAlater (AMBION) and stored at − 20 °C until qPCR analyses (performed in 2009). Nest size (number of cells) did not differ between unparasitized and parasitized colonies (nest size, Mann Whitney U = 53.5, p = 0.537; n = 14 vs 9); colony size (i.e. number of colony members) was slightly larger in unparasitized colonies (Mann Whitney U = 29.5, p = 0.036; n = 14 vs 9), so we included colony size in the model to account for this (see below).
To select the genes to test this prediction, we initially referred to the only study published at the time, which identified several genes differently expressed between the workers and the queen of Polistes canadensis (Sumner et al. 2006). At the time, neither genome nor transcriptome data for P. dominula nor P. sulcifer were available; accordingly, we tested different primers specific to P. canadensis on synthesized cDNAs from our two species of interest. Since genes associated with worker behaviour appear to be conserved across species, using genes identified in a congeneric species like P. canadensis is an appropriate approach. A first selection of candidate genes was made on the basis of the amplification of a single band observed on agarose gel (Supplementary Material). Then we isolated, purified and sequenced each single band to verify that the amplified genes were indeed homologous between the Polistes species (Supplementary Material S2). Finally, we aligned all sequences and measured the percentage of identity using CLUSTALW and nucleotide BLAST bioinformatic tools (Supplementary Material S2). Five candidate genes Apolipophorin III, ApoLp-III; Arrestin, ARR ; Heat shock protein 70KDa, HSP70; Imaginal disc growth factor, Idgf4; Major Royal Jelly Protein 1, MRJP1, (Table 2) were selected for further quantification of gene expression by q-PCR, using Tubulin alpha1 chain, TUBA1C as a control. On the basis of previous studies , four out of the five candidate genes were expected to be downregulated in workers from parasitized colonies compared to workers from unparasitized colonies (ARR, ApoLpIII, Idgf4, MRJP1) while one gene (HSP70) was expected to be upregulated (see Table 2).
RNA extraction from P. dominula and P. sulcifer and quantification of q-PCR product followed the protocol in . Briefly, total RNAs were isolated with TRIZOL (Invitrogen) from individual dissected brains and DNase treated (Turbo, PROMEGA) (Fig. 1c). Quantification was conducted using a Nanodrop and samples were checked for RNA degradation on an agarose gel. First strand cDNAs were synthesized with 1 μg of total RNA and reverse transcribed with Retrotranscriptase kit (AppliedBiosystem). cDNAs were amplified by quantitative PCR using the LightCycler-DNA Master SYBR Green I kit (AB) and the following primers: HSP70-5′-CTG TTC TTA GCG GCA ATG GTC-3′ and 5′-TGG ACA AAT CTA CTG GCA AGGAG-3′, Idgf4-5′-TGC CAT ACT TTT TCG ATG GGT-3′ and 5′-GCG TTC CAC CGA TAA TAG CTG-3′, TUBA1C-5′-AGC ACC ATC GAA TCG TAA GGA-3′ and 5′-ACA TTC CGA CTG CGC ATT TAT-3′, ApoLp-III-5′-AGC ACG AAT TTC GCC AAT AAA-3′ and 5′-CGC ATC TGG AAT TTG AGC ATT-3′, MRJP1-5′-GCT AGA CCA ACG ACA CCA TCATC-3′ and 5′-TTC GCG AAA AAT CCA AAA GC-3′, ARR-5′-TCG CAG TGT TGC ACC ACA A-3′ and 5′-TCA CTT TCT CTC AGG GCA AACTT-3′. Since primer sets were chosen to target gene regions where both species share > 90% similarity (Supplementary Material S2), comparable amplification efficiencies were obtained in both species. Gene specific standard curves were established using serially diluted cDNAs from a mixture of both host and social parasite brain RNAs. Q-PCR was performed on a LightCycler STEP-ONE (AB). Three technical replicates were run for both samples and standards. Cycle Threshold (CT) values were removed if amplification baseline was 15 < or > 35 and if melting curve patterns were not sigmoidal. Selected CT were submitted for quality control using Q-base software, and normalized against tubulin. Absolute quantitative means were derived from triplicate CT values obtained from the standard curve, normalized to the value of the tubulin gene. Stability of tubulin gene was validated using GeNorm software. Primers' specificity for each gene was checked by RT-PCR (PCR masterMix, PROMEGA) on both parasite and host cDNAs. Unique band amplification was confirmed on an agarose gel, after which PCR products were purified (purification kit, QUIAGEN), sequenced (BigDye kit, AB) and homology between parasite and host sequences was checked with CLUSTAL-W. After data quality control, sample sizes of workers from unparasitized and parasitized colonies for each of the five genes (see below and Table 2) were as follows: ApoLp-III: 26 vs 18; ARR : 11 vs 11; HSP70: 26 vs 18; Idgf4: 26 vs 18; MRJP1: 12 vs 12.
To account for the non-independence of data (i.e. workers belonging to the same colony), we used a generalized estimating equations (GEE) approach, which extends the generalized linear model to allow for analysis of correlated observations, as is the case for clustered data (Burton et al. 1998). This approach is robust against misspecification of the error structure model and has relaxed assumptions on data distributions (Hubbard et al. 2010;Overall and Tonidandel 2004).
We selected the best model on the basis of the "quasilikelihood under the independence model" criterion QIC, by choosing the model parameters that resulted in the smallest QIC (Pan and Wall 2002). We ran a GEE for each gene, with the following parameters: tweedie probability distribution, log link function, independent working matrix correlation. We ran a full factorial model which tested the effect of parasitism (fixed effect), colony size (covariate) and their interaction, by controlling for colony of origin (included as the subject variable). GEE were performed in SPSS 20.0 (IBM 2012).
All workers from unparasitized and parasitized colonies used in both behavioural and genomic experiments were dissected to confirm that they did not have developed ovaries nor yellow bodies (as sign of previous ovarian development), thus confirming they were not reproductively active (Cini 2014;Cini et al. 2013). During dissection workers were also checked for the presence of Strepsiptera endoparasites, which are known to be common in P. dominula (Hughes et al. 2004) and dramatically influence the wasp's behaviour and gene expression (Geffre et al. 2017). None of the workers included in the behavioural and transcriptomic analysis was found to be infected by Strepsiptera.

Results
None of the workers showed any signs of ovarian development in this study: all workers had threadlike ovarioles with no visible yolk nor yellow bodies.

Prediction 1: positions of workers differ in social networks of parasitized and unparasitized colonies
Colony size increased in both parasitized and unparasitized colonies during the experimental period (Wilcoxon test, W = 36, p = 0.008, n = 8). At the end of the experiment colony size did not differ between parasitized and unparasitized colonies (Mann-Whitney test, U = 5.5, p = 0.571, n = 4) (Supplementary Material S3).
Principal component analysis extracted three components with Eigenvalues greater than 1: PC1 (hereafter referred to as 'outgoing centrality') was mainly represented by outdegree, outcloseness, out eigenvector and betweenness, loadings range 0.604-0.852); PC2 (hereafter referred to as 'incoming centrality') was mainly represented by indegree, incloseness and in eigenvector, loadings range 0.717-0.808); finally PC3 was mainly correlated with clustering coefficient (loading 0.980). These three PCs performed well in condensing the social network metrics of individuals (variance explained 67.97%) (Supplementary Material S3).
Social network metrics differed between workers on parasitized and unparasitized colonies, suggesting workers were responding at the behavioural level to the social parasite during the early stages of parasitism. Early evidence of response was apparent at day 1 after usurpation, at which time PC3 (mainly clustering coefficient) differed significantly with workers from parasitized colonies showing higher clustering centrality than workers from unparasitized colonies; this difference was not apparent by day 7 and day 14 (Table 3, Supplementary Material S4). Later evidence of response was apparent at day 14, where workers on parasitized colonies showed significantly greater outgoing centrality (PC1) and significantly higher incoming centrality (PC2) than workers from unparasitized colonies (Table 3, Supplementary Material S4); there was no significant difference before usurpation, at day 1 and at day 7. Colony of worker origin had no significant effect for any of the social network measures (Table 3, Supplementary Material S4) (Fig. 2). Overall, these data demonstrate that workers from parasitized colonies show significant modifications in their social behaviour within the dominance interaction network compared to workers from unparasitized colonies. In particular, they showed (at specific time-points) more interconnected networks and higher number of individuals performed and received dominance interactions.

Discussion
In the strong coevolutionary arms race between social parasites and their hosts, hosts show a wide suite of defensive strategies, from avoidance strategies to reduce the chance of parasitism, to rebellion after parasitism has occurred (Brandt et al. 2005;Cervo 2006;Cini et al. 2019b;Grüter et al. 2018). The timing of host defence affects host fitness and a fast response to parasitism may help limit the costs incurred by the host due to the invasion. When parasitized by the social parasite P. sulcifer, P. dominula host workers show a covert reaction which consists of ovarian development, such that they are ready to reproduce if the opportunity arises ). Here we highlighted that well before this physiological response is evident, significant modifications in worker social behaviour and brain gene expression occur. To our knowledge, this is the first evidence of early behavioural and neurogenomic responses by host workers towards the presence of an inquiline social parasite in a social insect.
Polistes workers from parasitized colonies were thought to not exhibit changes in their behaviours after being invaded by a social wasps, continuing to perform typical worker task, such as foraging and nest activities (Cervo 2006). Our timeseries dissection of social interactions, however, reveals subtle changes in their social behaviour from within a day of usurpation. On the day after the usurpation by the social parasite, workers showed significant changes in their clustering coefficient: workers from usurped colonies exhibited denser (more interconnected) social networks than workers from unparasitized colonies. This could indicate that workers in parasitized colonies detect the intrusion of the social parasites from the outset of usurpation. The effect disappears by day 7. It is unclear at this stage whether the change in worker social networks on day 1 after usurpation is a response to the threat of the parasite specifically (i.e. they detect her as an alien intruder), or simply a response to a social perturbation that could equally arise within an unparasitized colony. In the former case, the response could be a reflection that the parasite has not yet become chemically integrated into the host colony: this is only achieved a few days after usurpation (Dapporto et al. 2004;Sledge et al. 2001;Turillazzi et al. 2000). The higher levels of connectedness amongst nest mates could be a form of information transfer and/or an attempt to reinforce the colony odour to defend against the successful take-over of the intruder. Conversely, the response on day 1 could simply reflect social disruption that occurs when a queen is lost and replaced, which may  Table 3 for details) Fig. 3 a Comparison in brain gene expression level for the target gene Idgf4 between workers from unparasitized (blue squares) and parasitized (red circles) colonies. Box plot reports the median value, the interquartile range and outliers are indicated with black circles. Bold font indicates significant differences for that target gene (see text and Table 3 for details); (b), (c) and (d) scatterplot of gene expression level vs colony size for workers from unparasitized (blue squares) and parasitized (red circles) colonies for the target genes HSP70, Idgf4 and ARR disrupt the dominance hierarchy (Jandt et al. 2014;Pardi 1948;Strassmann et al. 2004), resulting in an increase in RDI among colony members as they try to re-establish their dominance position in the new network. We lack social network analyses of the behavioural dynamics that take place after queen loss and replacement in Polistes; such studies would help determine whether the changes in centrality observed in parasitized workers is typical or not of queen succession. Equally, the parasite herself elicits high rates of RDI as she is particularly active during the first days of usurpation as she tries to dominate all individuals on the nest (Cervo 2006;Turillazzi et al. 1990): workers may be responding to general social perturbation within their worker ranks, rather than specific detection of the parasite per se. Future studies are required to directly assess whether such an increase in workers clustering is specific to the social disruption due to the parasite presence or if it is a generic response to any disturbance of colonial stability. These arguments cannot, however, explain the result that workers from parasitized colonies also showed differences in their centrality within the colonial dominance network compared to unparasitized colonies two weeks after parasitism. We found that workers from parasitized colonies, by day 14 after usurpation, were significantly more central to the social network than workers from unparasitized colonies. In particular, workers from parasitized colonies showed a higher outgoing centrality (i.e. instigated more dominance interactions) and also a higher incoming centrality (i.e. received more dominance interactions) in the colonial dominance network compared to workers of unparasitized colonies. This finding is intriguing, as (a) dominance behaviour and reproduction are positively correlated in this species (Dapporto et al. 2010;Reeve et al. 1991), (b) a dominance hierarchy is established among workers, which regulate the access to reproduction in case the breeder is to be replaced (Jandt et al. 2014;Pardi 1948) and finally (c) more central individuals in dominance networks are typically the reproductive ones in wasps, as shown in Ropalidia species (Bhadra and Jordán 2013;Bhadra et al. 2009). Taking these factors into account, the finding that workers in parasitized colonies have higher outgoing and incoming centrality (i.e. perform and receive more dominance interactions) than workers from unparasitized colonies might suggest that the presence of the parasite leads to a social re-organization within the colony, with workers trying to escalate the dominance hierarchy to ensure future access to reproduction. The fact that workers in parasitized colonies showed higher levels of both performed and received dominance interactions might suggest that a clear dominance hierarchy was not yet achieved among workers. Our findings indeed parallel what has been found in Ropalidia wasps, where a reorganization of the dominance network, with increased centrality, has been shown after experimental breeder removal (Bhadra and Jordán 2013).
We also found some evidence that workers respond to parasitism from their brain gene expression. First, the gene Idgf4 was downregulated in parasitized workers. Idgf4 is known to be involved in cuticle formation, organization of proteins in the chitin-matrix, and in adult moulting (Gu et al. 2019); this gene is putatively involved in responses to changes in the social environment in Polistes wasps . With respect to reproduction, in P. dominula Idgf4 is upregulated in workers relative to queens in unparasitized colonies . A downregulation in the expression of Idgf4 in parasitized workers, therefore, may suggest that parasitized workers downregulate worker-biased genes, possibly instigating the switch from helper to reproductive. With respect to social environment, Idgf4 is upregulated in the social parasites when they shift from a solitary to social environment, upon taking over the colony and integrate into the host colony . The finding that workers downregulate Idgf4 expression in response to parasitism may indicate that, in the presence of the parasite, workers experience a change in the social environment. This is in accordance with the behavioural data, which shows behavioural changes in the colonial social network and signature of a social rearrangement.
The second gene whose expression pattern was influenced by parasite presence in our experiment was Arrestin; this has several roles of potential importance here. Firstly, the ARR gene family is involved in visual and olfactory reception (Dolph et al. 1993;Merrill et al. 2002), it is linked to division of labour in bees, being upregulated in foragers and more defensive bees Sasagawa et al. 2003). In Polistes wasps, ARR was found to be associated with reproductive phenotype ) and, specifically in the parasite species, was found responsive to changes in the social environment . While there was not an overall significant effect of parasitism on ARR expression, there was a significant crossover effect with colony size (i.e. a significant interaction). In workers from unparasitized colonies the expression of ARR increased with colony size, this effect was not detected in workers from parasitized colonies. The increased expression of ARR in bigger colonies (as found in workers from unparasitized colonies) thus makes sense under this perspective, as bigger colonies represent more crowded social environments, where more visual and chemical stimuli for social recognition are indeed to be processed. The finding that such a relationship does not hold for parasitized colonies is intriguing and suggests that, perhaps, social environment is disrupted in parasitized colonies such that workers do not (or cannot) respond to colony-size effects.
The neurogenomic responses of workers to social parasites reported here offer a tantalising first-look at the interactions between social parasites and their hosts, and support the behavioural responses we also report. However, we investigated only five candidate genes, and so our results represent a first step at understanding the response to parasitism. We hope these findings will encourage future work in this area, especially taking a genome-wide perspective (e.g. using RNAseq to sequence whole brain transcriptomes) along a temporal gradient from usurpation to the onset of evident, physiological modification in the parasite workers (6 weeks) . Matched analyses of genomic and behavioural changes over such a time course will also enable the critical question of whether the observed responses of parasitized workers are due to detection of the parasites (and thus a co-evolved response), or simply a response to conspecific colony perturbations. Genomic analyses will help reveal the nature of any co-evolved mechanism, e.g. in olfactory proteins, such as chemosensory and odorant binding molecules Pelosi et al. 2005Pelosi et al. , 2006Pelosi et al. , 2014Zhu et al. 2019, but also see Xiao et al. 2019 for recent findings on the controversial role of odorant binding proteins) that may allow the host workers to detect the parasite.
The idea that workers might be able detect the social parasite would then raise interesting and important questions about the cues that workers use to do so. Workers could recognize the presence of the parasite by two main sensory modalities: smell, by detecting cues associated with the parasite itself or with its brood (eggs, larvae and pupae) Cini et al. 2014;Lorenzi 2006) or sight (Cini et al. 2019), by recognising the facial pattern of the parasite, a major point of difference in parasite and host morphology Ortolani et al. 2010). The timing of workers response (day 1 and day 14) suggests that visual (facial recognition) and chemical cues (recognition of the odour of the adult) might be the cues informing workers of parasite presence. Indeed, facial pattern is evidently a static cue, which does not change with time and can be assessed from the beginning of the host-parasite interaction [and indeed it plays a role in parasite infection, reducing host aggressiveness toward social parasites (Ortolani et al. 2010)]. Similarly, the odour of the social parasite is expected to play a role in the very first days after the usurpation since P. sulcifer is able to chemically integrate into the host colony within 3-4 days after usurpation. Therefore, we would expect a behavioural reaction to occur rapidly after usurpation if parasite odour was the detection cue. The idea that parasite brood odour is the key cue is not supported by our findings. Indeed, eggs and larvae of the social parasites are unlikely to allow detection, as the odour of parasite eggs is very similar to that of host eggs, and larvae escape detection by bearing a meaningless odour Dani et al. 2004). The remaining developmental stages of parasite brood, i.e. pupae and newly emerged brood, will not be present in the colony until, on average, 15 and 25 days ). This suggests that chemical cues from parasite brood are unlikely to allow detection of the parasite, leaving open the possibility that social parasites are detected because of some visual or odorant cue of the social parasite itself.
Overall, our study shows that host workers respond at the behavioural and brain transcriptomic level to social parasitism, from the outset. We cannot yet determine which of these responses are a result of a co-evolutionary adaptation by hosts to minimise costs of the parasite, or if they are a typical response to within-colony perturbations that any unparasitized colony could have experienced. A co-evolved response should carry strong selection pressures because an immediate and overt reaction to the social parasite would be risky for the workers. Indeed, P. sulcifer is bigger than the host, has a greater fighting ability and can easily defeat up to 11 workers (Cervo and Turillazzi 1996;. A subtle response, as observed here, therefore may be adaptive. Our results provide the first evidence of early behavioural and neurogenomic responses of host workers toward the presence of an inquiline social parasite in a social insect, highlighting the importance of considering the timing of defensive responses when analysing the ongoing co-evolutionary arms race between hosts and their social parasites. 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://creat iveco mmons .org/licen ses/by/4.0/.