Increased immune marker variance in a population of invasive birds

Immunity and parasites have been linked to the success of invasive species. Especially lower parasite burden in invasive populations has been suggested to enable a general downregulation of immune investment (Enemy Release and Evolution of Increased Competitive Ability Hypotheses). Simultaneously, keeping high immune competence towards potentially newly acquired parasites in the invasive range is essential to allow population growth. To investigate the variation of immune effectors of invasive species, we compared the mean and variance of multiple immune effectors in the context of parasite prevalence in an invasive and a native Egyptian goose (Alopochen aegyptiacus) population. Three of ten immune effectors measured showed higher variance in the invasive population. Mean levels were higher in the invasive population for three effectors but lower for eosinophil granulocytes. Parasite prevalence depended on the parasite taxa investigated. We suggest that variation of specific immune effectors, which may be important for invasion success, may lead to higher variance and enable invasive species to reduce the overall physiological cost of immunity while maintaining the ability to efficiently defend against novel parasites encountered.

down regulated inflammatory responses which are costly might lead to an increase in less costly responses (e.g. antibody-mediated immunity) 9 . However, in contrast to the revised-EICA, Cornet et al. 21 suggest that newly acquired local parasites might induce an equivalent immune response in invasive and native populations. Also, Brown and Shine 22 suggest that trade-offs between different immune effectors are more important for invasive species than the general downregulation of costly traits. Moreover, Møller and Cassey 23 propose that invasive bird species with strong immune response towards novel parasites are potentially more likely to become successful invaders. A better understanding of the role and complex interplay of parasites and immune effectors in the context of invasion is clearly necessary.
The main goals of this study were to compare parasite prevalences and immune effector levels between native and invasive Egyptian geese (Alopochen aegyptiacus) and shed light on their complex interplay in the frame of invasion. We tested (a) the prediction derived from the ERH that invasive populations are less infected than native populations and (b) the predicted differences in costly and less costly immune responses between native and invasive populations based on the revised-EICA hypothesis. Additionally, to test for potential differences in the strength of immune response between native and invasive birds, we did not restrict our analyses to investigating changes in mean levels of immune effectors. Instead, as an additional exploratory analysis, we studied how the variance in immune effectors differs between native and invasive populations, expecting changes in the variance of immune effectors in addition to changes in the mean levels. We chose a native and an invasive population of the Egyptian goose as one of the most successful invasive bird species in Europe. Its European population continues to increase 24,25 and it has been listed as an invasive species in Germany according to "EU Regulation No 1143/2014 on the prevention and management of the introduction and spread of invasive alien species" and the BfN script 47 in 2017 26 .

Material and methods
Study species. The Egyptian goose is native to Africa and was introduced to Europe in the twentieth century 27 . Its native population is distributed on the sub-Saharan African continent and is one of the most common and wide spread African waterfowl species. Egg laying occurs throughout the whole year with a peak between late winter and early summer 28 .
The neozootic population invades Europe eastwards starting from the Netherlands, where they were introduced as ornamental species to parks 27 . It is now one of the most spreading neozootic bird species in Europe 24 . From the 1980s Egyptian geese also invade Germany where its population size increased rapidly 29,30 . The Egyptian goose is a resident (non-migratory / short distance migratory), monogamous, territorial bird species occurring as neozootic species in a variety of water habitats (e.g. streams, rivers, ponds, lakes,) in Europe 31 .
Sampling. Parasite prevalence and immunity of Egyptian geese from a native population in Namibia were investigated and compared to those of a currently spreading invasive population of the same species in Germany. In both regions, geese were sampled during ringing procedures (live trapped) or dissected after general pest control hunting (necropsy). Blood samples for immunological assays and serology exclusively originate from live trapped individuals whereas macro-parasite investigation was performed during necropsy. Micro-parasite investigation was performed in birds from both groups. Therefore, the resulting datasets are analysed separately (see method section: Statistical analysis) but a potential interplay between the different parasite prevalences and immune results is evaluated in the discussion.
Live trapping. Twenty-one Egyptian geese (9 male, 12 female) were live trapped in Namibia (22.35° S, 17.05° E) (native range) in February 2016. Additionally, data from a subset of 110 adult Egyptian geese from Germany (65 male, 45 female) investigated by Prüter et al. 32 were included. German geese were sampled in the Rhine and Mosel areas (50.4° N, 7.6° E) (invasive range) in 2015 (n = 78) and 2016 (n = 32) in different months (supplementary data Table S1). Sex and reproductive status were recorded. Reproductive status was defined as breeding (e.g. guiding gosling, showing territorial behavior with a partner, having an egg-laying active cloaca) or nonbreeding (e.g. not fulfilling criteria of breeding and/or being part of a non-family-flock). Numbers of breeding vs. non-breeding individuals can be found in Table S1. All Namibian birds were likely non-breeding individuals, which were sampled at a cattle feedlot were thousands of birds fly in to feed on the corn provided to the cattle. Blood was drawn from the vena metatarsalia plantaris superficialis using needles with a diameter of 0.06 mm for males and 0.04 mm for females. A fresh blood smear was prepared at capture and air dried. Blood samples were kept at 4-8 °C, centrifuged and sera were frozen in liquid nitrogen within eight hours after blood draw. Pharyngeal swabs were collected using sterile cotton swabs. Once field work finished, samples were transported to the Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany and sera, blood clot and pharyngeal swabs were kept frozen at -80 °C. Sampling in Germany was authorized by the Landesuntersuchungsamt Rheinland-Pfalz (G 15-20-005) and Landesamt für Natur, Umwelt und Verbraucherschutz Nordrhein-Westfalen (LANUV) (84-08.04.2015.A266). Permission to collect samples in Namibia was granted to GM and HK by the Ministry of Environment and Tourism (MET). Permission to export sample material from Namibia was granted by a MET export permit (No. 107513), and samples were transported to Germany in compliance with the Nagoya Protocol on Access and Benefit-sharing. All experimental procedures described in the materials and methods section were approved by the Internal Committee for Ethics  www.nature.com/scientificreports/ Egyptian geese (11 male, 16 female), which were shot in February 2016 during regular pest control in Central Namibia were dissected. One of twenty-seven Namibian birds was live trapped and sampled before death and is thus included in both groups (live trapped and necropsy). Geese from Germany were kept frozen at -20 °C after hunting until further analysis. Namibian geese were dissected immediately post mortem. During necropsy, ectoparasites, intestinal helminths and nasal leeches were collected. Additionally, pharyngeal swabs were taken for molecular analyses.
Parasitological and microbiological analysis. Both macro-parasites (ectoparasites, nasal leeches (Euhirundidae), intestinal helminths) and selected micro-parasites (blood parasites (Haematozoa), bacteria, viruses) of Egyptian geese from the two populations were characterized for community composition and prevalence (methods see Table S2). Hereafter we use the term "parasites" combining macro-and micro-parasites and only explicitly distinguish between the type of parasites when differences can be expected and/or occur. During necropsy, wing and breast feathers were macroscopically checked for the presence of ectoparasites. The upper beak was cut open and macroscopically investigated for the presence of nasal leeches. Intestinal helminths were extracted from the intestine of the birds and were determined to the family level based on morphology 33 . Additionally, blood smears of all live-trapped animals were investigated for the presence of blood parasites during white blood cell counts 34 .
To compare with previously determined bacterial prevalence of adult German Egyptian geese 32 (Table 2), the Namibian birds were screened for Mycoplasma spp. and Riemerella (R.) anatipestifer using conventional 16S rRNA-based PCR assays as described by Prüter et al. 32 . To verify the specificity of the Mycoplasma PCR assay, products with a clear band were further investigated by sequence analysis, following the procedure described by Prüter et al. 32 . Only samples with a clear sequencing result were designated positive.
The seroprevalence of antibodies (Ab) against Influenza A virus (IAV), Avian avulavirus 1 (AAvV-1) and West Nile virus (WNV) were determined 32 . For the detection of Abs against IAV, a commercial competitive enzyme linked immunosorbent assay (ELISA) was used following the manufacturer instructions (ID.vet, Grabels, France, Influenza A Antibody competition, FLUACA ver 0917DE). A commercial competitive ELISA for detection of Abs against AAvV-1 (Avian paramyxovirus 1; syn. Newcastle disease virus) was used according to the manufacturer protocol (ID.vet, Grabels, France, Newcastle Disease Competition, NDVC ver 0913 DE). Commercial competitive ELISA for Abs against Flaviviridae including WNV were applied following the manufacture protocol (ID. vet, West Nile Competition, WNC ver 1014-1P DE).
Immunological assays. Several eco-immunological tests were used to quantify both the cellular and humoral parts of the acquired and innate immune responses of Egyptian geese 35 . Most of the methods are speciesnon-specific and have been used in a wide variety of free-living avian species, including different waterfowl [36][37][38] . We quantified the amounts of different humoral (natural antibodies, complement, lysozyme and haptoglobin) and cellular (monocytes, heterophils, eosinophils and basophils) effectors of innate immunity. For adaptive immunity we measured the total immunoglobulin Y (IgY) concentration and the number of lymphocytes 36 . Sample sizes (n) for each assay were dependent on the total amount of serum available from each individual and therefore differ among the tests (Table 1).
Immunoglobulin Y. Total IgY, the avian equivalent to mammalian IgG, was measured using a sensitive ELISA with commercial anti-chicken antibodies 38,40 . Ninety six-well high-binding ELISA plates (82.1581.200, Sarstedt) were coated with 100 µl of diluted serum sample (2 samples per bird 1:16,000 diluted in carbonate-bicarbonate buffer) and incubated first for 1 h at 37 °C and then overnight at 4 °C. After incubation, the plates were washed with a 200 µl solution of phosphate buffer saline and PBS-Tween, before 100 µl of a solution of 1% gelatine in PBS-Tween was added. Plates were then incubated at 37 °C for 1 h, washed with PBS-Tween and 100 µl of polyclonal rabbit anti-chicken IgY conjugated with peroxidase (A-9046, Sigma) at 1:250 (v/v) was added. Following 2 h incubation at 37 °C, the plates were washed again with PBS-Tween three times. After washing, 100 µl of Table 1. Total sample sizes (total n), sample sizes grouped by sex (sex ratio (♂, ♀)) and year of sampling of blood and serum samples from Namibian (native) and German (invasive) Egyptian geese (Alopochen aegyptiacus) for each immunological effector grouped by the costs of immunity (low costs vs. high costs according to Klasing 39 and Lee and Klasing 9 .  www.nature.com/scientificreports/ revealing solution [peroxide diluted 1:1000 in ABTS (2,20-azino-bis-(3-ethylbenzthiazoline-6-sulphonic acid))] was added, and the plates were incubated for 1 h at 37 °C. The final absorbance was measured at 405 nm using a photometric microplate reader (μQuant Microplate Spectrophotometer, Biotek) and subsequently defined as total serum IgY levels 41 .
Lysozyme. To measure lysozyme concentration in serum, we used the lysoplate assay 37 : 25 μl serum were inoculated in the test holes of a 1% Noble agar gel (A5431, Sigma) containing 50 mg/100 ml lyophilized Micrococcus lysodeikticus (M3770, Sigma), a bacteria which is particularly sensitive to lysozyme concentration. Crystalline hen egg white lysozyme (L6876, Sigma) (concentration: 1, 1.25, 2.5, 5, 6.25, 10, 12.5, 20 and 25 µg/ml) was used to prepare a standard curve for each plate. Plates were incubated at room temperature (25-27 °C) for 20 h. During this period, as a result of bacterial lysis, a clear zone developed in the area of the gel surrounding the sample inoculation site. The diameters of the cleared zones are proportional to the log of the lysozyme concentration. This area was measured three times digitally using the software ImageJ (version 1.48, http://image j.nih.gov/ij/) and the mean was converted to a semi-logarithmic plot into hen egg lysozyme equivalents (HEL equivalents, expressed in μg/mL) according to the standard curve 42 .
Haemolysis-haemagglutination assay. The levels of the natural antibodies and complement were assessed by using a haemolysis-haemagglutination assay as described by 43 adjusted to the limited volume of serum. After pipetting 15 μl of serum into the first two columns of a U-shaped 96-well microtitre plate, 15 μl sterile PBS were added to columns 2-12. The content of the second column wells was serially diluted (1:2) until the 11th column, resulting in a dilution series for each sample from 1/1 to 1/1024. The last column of the plate was used as negative controls, containing PBS only. Fifteen μl of 1% rabbit red blood cells (supplied as 50% whole blood, 50% Alsever's solution, Envigo) suspension was added to all wells and incubated at 37 °C for 90 min. After incubation, in order to increase the visualisation of agglutination, the plates were tilted at a 45° angle at room temperature. Agglutination and lysis, which reflect the activity of the natural antibodies and the interaction between these antibodies and complement 43,44 , was recorded after 20 and 90 min, respectively. Haemagglutination is characterised by the appearance of clumped red blood cells, as a result of antibodies binding multiple antigens, while during haemolysis, the red blood cells are destroyed. Titres of the natural antibodies and complement were given as the log2 of the reciprocal of the highest dilution of serum showing positive haemagglutination or lysis, respectively 43,45 .
Haptoglobin. We measured haptoglobin concentrations with a commercial kit (TP801, Tri-Delta Diagnostics, Inc.) following the instructions of the manufacturer. Haptoglobin concentrations (mg/ml) in undiluted serum samples were calculated according to the standard curve on each plate 36 . www.nature.com/scientificreports/ White blood cell counts. To count leucocytes, blood smears were prepared, air-dried and stained using Giemsaand May-Grünwald staining. Smears were examined at 1,000 fold magnification with oil immersion and the relative number and types of leucocytes were assessed by counting 100 leucocytes. The number of white blood cells of different types was expressed per 10 4 erythrocytes 45 .

Statistical analyses. Parasite prevalence.
To investigate potential differences in the prevalence of parasites between native and invasive Egyptian geese, we used Fisher's exact tests because relatively low sample sizes of dissected animals did not allow to perform multivariate analysis.
Immunity. The means and variances of the different immune effectors were compared between the invasive and native Egyptian geese populations. To this end, we used an extension of commonly applied linear models. Linear models assume that the response variable y is a function of a linear combination of n predictor variables x with coefficients c 0 ,..,c n and an error ε: where ε is the so-called residual variance which captures all the variation in the response variable that is not explained by the predictors. In linear models ε is assumed to be normally and independently distributed around zero. An additional usual assumption is that the variance of this distribution is a constant σ 0 , i.e.: which corresponds to the assumption of normally distributed residuals with homogeneous variance. Thus, the estimated effects of the predictors c 1 ,..,c n only describe changes in the mean of the response variable, but not in the variance around that mean.
In our analysis, models were used in which the variance was allowed to be a linear function of some predictor variables z (which might be the same or different from the predictors x of the mean in Eq. 1), i.e.: Thus, we were able to estimate simultaneously the effect and respective p-values of predictors x upon the variation in the mean of the response variable (Eq. 1) and also the effects and respective p-values of predictors z upon the residual variation around that mean (Eq. 3).
In order to appropriately capture all the potential variation in the response variables we used linear mixedeffects models (LMMs), which in addition to fixed effect predictors in Eq. (1) also included a random effect as a predictor of the mean. However, for enhanced clarity random effects were omitted in the equations above.
Different immune effectors were used as response variables ( Table 1). As predictors for the mean sex (male vs. female), reproductive status (breeding vs. non-breeding) and invasion status (native vs. invasive) were included as fixed effects and month of sampling as a random effect. In this way, we control for potential confounding effects of breeding status between the two populations. As predictors for the variance, we included invasion status (native vs. invasive), sex and reproductive status, which allowed us to test our prediction that the variance in immune effects is higher among invasive individuals compared to native individuals.
Some of the immune effectors were transformed (see tables supplementary data S3, S4, S5) to ensure normality of residuals. For haptoglobin we were not able to perform a transformation that ensured normality, because of the high proportion of values below the detection threshold. To account for this, we performed a generalized linear mixed model (GLMM) with a binominal error distribution and with a binary response variable (haptoglobin being either above or below the detection threshold). In this model, it was necessary to constrain the error variance to a fixed value. Thus, for haptoglobin we were only able to test for a change in mean but not for a change in variance. In addition to analysing total leucocytes, different leucocyte subtypes were analysed separately.
The LMMs and GLMMs were implemented using the R package glmmTMB version 0.2.0 46 . To test for differences in residual variance we used the option disformula in the function glmmTMB. Potential collinearity of predictors was tested by calculating variance inflation factors using the R package car version 2.1-6 47 . All statistical analyses were performed using R version 3.3.2 48 .

Results
Parasite prevalences. Blood parasites and all groups of macro-parasites (ectoparasites, nasal leeches, intestinal cestodes and trematodes) were found at lower prevalence in the invasive than in the native population of Egyptian geese with the exception of intestinal nematodes. However, none of the observed trends reached statistical significance ( Table 2). The bacterium R. anatipestifer, was found in 67.0% of adult invasive Egyptian geese from Germany 32 but was not detected in native geese from Namibia. This difference was statistically significant ( Table 2). Mycoplasma spp., which was not detected in the German geese 32 , was detected in two individuals from Namibia. However, the result was not statistically significant ( Table 2). Seroprevalence of selected viruses (IAV, AAvV and WNV) were lower in geese from the invasive population but the difference was only statistically significant for antibodies against IAV (Table 2).

Immunity. Differences in mean immune function measures.
Of the four assays measuring 'low cost' immune effectors, mean haemolysis was significantly higher in the invasive population. Mean IgY, mean lysozyme and mean haemagglutination did not significantly differ between the two study populations (Fig. 1). Among the 'high cost' immune effectors, the mean number of leucocytes was not significantly different between the two study populations (Fig. 1). Detailed analyses of the means of differential white blood cells (Fig. 2) demonstrated (1) y i = c 0 + c 1 x 1,i + · · · + c n x n,i + ε i , ε i = N 0, σ 0 + σ 1 z 1,i + · · · + σ n z n,i .

Scientific Reports
| (2020) 10:21764 | https://doi.org/10.1038/s41598-020-78427-7 www.nature.com/scientificreports/  Tables S3, S4). For sample sizes of each immune measure for the two populations see Table 1. www.nature.com/scientificreports/ significantly higher mean numbers of heterophils and lymphocytes in the invasive population, whereas mean eosinophil concentration was significantly lower. No statistically significant differences in mean monocyte and basophil numbers were observed between the two groups. In contrast, haptoglobin abundance was significantly higher in the invasive (65%) than in the native (7.3%) population (Fig. 3). Males showed significantly higher mean levels of IgY, haptoglobin abundance, total leucocytes and lymphocytes than females. Mean lymphocytes were significantly higher in birds that were in breeding status than in non-breeding individuals. No statistically significant effects on the other means of measured immune effectors were associated with sex or reproductive status (supplementary data Tables S3, S4 and S5).
Variance of immune effector measures. Among the 'low cost' immune effectors the variance of IgY was significantly higher in the invasive population. However, the variance of lysozyme, haemagglutination, haemolysis did not differ significantly between the two groups (Fig. 1). For the 'high cost' immune effectors, the variance of total leucocytes did not differ significantly between the two groups (Fig. 1). The variance in the concentration of lymphocytes and heterophils was significantly higher in the invasive geese. Variance in basophil, eosinophil and monocyte counts was not significantly different between the groups (Fig. 2).
Males showed signifficantly higher variance of haemagglutinaion than females. The reproductive status of the birds significantly influenced the variance of three 'low cost' immune effectors. The variance of IgY was signifficantly higher in breeding than in non-breeding birds whereas the variance of lysozyme and haemagglutination was signifficantly lower in breeding birds than in non-breeding individuals (supplementary data Tables S3, S4 and S5).

Discussion
The results are not consistent with the revised-EICA. Four out of eleven immune effectors had higher mean values in invasive population when compared to their native conspecifics. Three of these immune effectors are considered energetically costly 9 , which is inconsistent with the predictions of the revised-EICA hypothesis which predicts a shift to low cost immune effectors. In line with the ERH, all but one of the parasites examined decreased in prevalence when comparing native and invasive populations of Egyptian geese, although this effect was not statistically significant. Only the bacterium R. anatipestifer was absent in the native population while the invasive population showed a 67% prevalence, which could indicate that the invading geese have encountered a novel pathogen during the invasion process 32 .
Most studies of invasive vertebrate hosts and their pathogens focus on helminths, which elicit a host Th2 response (e.g. promotes antibody production) and increase eosinophil numbers, and show reduced helminth infection and reduced eosinophil numbers in invasive populations compared to native 14,18,49,50 . Previous support for the revised-EICA related predictions have focused on macro-parasites, which might be biased toward the pathogens and immune effectors investigated. We observed low mean eosinophil numbers in invasive geese indicating a lower impact of macro-parasites on the immune system than in the native population (consistent with the ERH and the revised-EICA).
R. anatipestifer is a bacterium of relevance for domestic ducks and geese, leading to severe clinical symptoms 51,52 . The lack of clinical symptoms in the invasive Egyptian geese infected with R. anatipestifer may indicate that they are more resistant or tolerant than expected 21,53 or this waterfowl may carry R. anatipestifer as a commensal bacterium 32,54 . Overall, immune defences effective against micro-parasites, especially bacteria (heterophils, lymphocytes, haemolysis) were found to be higher on average in the invasive than in the native goose populations. Haptoglobin in particular, which is an energetically costly acute phase protein with bacteriostatic function 36 , was significantly more abundant in the invasive than in the native population of Egyptian geese which is inconsistent with the EICA-related predictions.
Studying a variety of micro-and macro-parasites and high numbers of different immune effector levels-as done in this study-indicate that presence or absence of support for existing hypotheses likely depends on the parasite and immune response investigated. This demonstrates that existing hypotheses predicting immune effector levels and parasite prevalence in invasive populations are likely too simplistic considering the loss and gain of parasites as well as the complexity of immune responses.
Immune effector variance. The variance for three out of ten immune effectors in the invasive population of Egyptian geese was higher than in the native population. In our opinion, existing hypotheses, such as the revised EICA, insufficiently cover the complex interplay of loss and gain of parasites in invasive populations. The higher variance of some immune effectors in the invasive population might indicate higher immune plasticity on the individual level (as measured here) which could reflect higher variation in immune reaction on the population level. Higher variation in immune reaction might enable invaders to efficiently defend against novel pathogens while simultaneously reducing the overall costs of immunity. The complexity of the immune system and the specificity of different immune effectors to specific parasites may be more important during invasion than previously acknowledged 55 . Furthermore, White and Perkins 56 suggested that higher plasticity of individual immune response in a founder population, the more likely the species will successfully invade. Ghalambor et al. 57 hypothesized that the period of persistence of invasive species (period after introduction and before rapid population growth starts) might be dependent on phenotypic plasticity. Changes in immune investment and variation of immune reactions in the different stages of invasion are likely (for the different phases of invasion see Duncan et al. 5 ), although the underlying cause of a variation of immune reactions might differ among invasive populations (e.g. genetically heterogeneous populations, diversity of outside stimuli). Variation of immunity, as opposed to generally reduced investment in immune function, might provide invading species with the necessary flexibility to colonize novel environments. Specifically, flexibility of immune reactions could enable invad-

Scientific Reports
| (2020) 10:21764 | https://doi.org/10.1038/s41598-020-78427-7 www.nature.com/scientificreports/ Figure 2. Differences in distributions of the differential white blood cells (nr per 10 4 erythrocytes) (a-e) between native and invasive Egyptian geese are shown; red = native; blue = invasive; P = p-values for the effects of population (native vs. invasive) on the mean (P mean ) and variance (P variance) of the respective immune measure from the GLMMs (see Table S5). For sample sizes of each immune measure for the two populations see Table 1. www.nature.com/scientificreports/ ers to balance reduced immune investment against their original parasites with increased defence against novel parasites. Individuals not confronted with novel parasites could reduce their immune investment (in line with the revised-EICA hypothesis), whereas infected individuals could selectively increase the appropriate immune defense. Here, a significantly higher mean value of heterophils, as effectors against bacterial infections, with evidence for higher variance in the invasive population indicate that pathogenic bacteria (as the here selected R. anatipestifer) might be particularly important to the invasive population. Depending on the immune effector and its role in the immune response against parasites and the stage of the invasion process we propose three scenarios consistent with our results (Fig. 4b-d). If the enemy release effect is dominant (e.g. early stages of invasion; invasion front), immune investment would be down-regulated compared to native populations (Fig. 4a) as suggested by the revised-EICA hypothesis (Fig. 4b). In this case, low parasite prevalence and reduced overall costs of immunity would be expected. During range expansion, members of invasive populations would likely increasingly face novel parasites. Invasive animals must have effective defences against diverse newly encountered parasites. Heterogeneous prevalence for parasite species would be expected with greater heterogeneity expected during the early stages of invasion. Some invading individuals would need to be able to defend themselves against the new parasites requiring higher immune investment than observed in native individuals 23 . Individuals of the invading population not infected may invest less in immunity compared to their native relatives. Variation in immune investment would be expected to be prominent during early stages of invasion and accordingly we would expect an increased variance in immune measures (Fig. 4c). In later phases, prevalence of parasites might increase further and thus, a need to defend against new parasites (Fig. 4d).
Eco-immunology of Egyptian geese. The Egyptian goose (Alopochen aegyptiacus) is a successful invasive bird species having and continuing to spread throughout Western Europe since the mid-twentieth century 25 . The impact of newly acquired pathogens on this population is most likely still ongoing and potentially increasing (middle phase of invasion, Fig. 4c). That the invasion is still on-going likely explains why increased immune effector variance is prominent. Our results indicate that among-individual variation in immunity is on average higher in invasive compared to native Egyptian geese but that these effects might be parasite dependent and therefore the effects differ among immune effectors.
In contrast to parasite infection, the higher bactericidal investment but less variance in bacterial related immunity, might indicate that population growth rates will slow in the future, as immune investment increases and consequently energy reallocates away from reproduction. After the invasion has progressed further or has been completed, the variance of all effectors is expected to decrease representing adaptation to the new environment (Fig. 4d).

Figure 3.
A barplot for the differences in haptoglobin abundance (percent of individuals with haptoglobin concentrations above the detection threshold) between native and invasive Egyptian geese is shown; red = native, blue = invasive; p-value for the effects of population (native vs. invasive) on the mean haptoglobin from the binomial distribution model from the GLMMs (see Table S4). (It was not possible to fit a Gaussian model for haptoglobin concentration. Thus, the assessment of the difference in variance was not possible). www.nature.com/scientificreports/ When interpreting our results, it is important to consider the limitations of our study. First, sample sizes, especially of dissected animals, are low and immune measures and macro-parasite data originate from different individuals and thus preclude a direct analysis of the interplay between immune marker level and infection status. Additionally, low sample sizes and sample size differences between the populations may impact the statistical power of the analysis. For future studies, we recommend avoiding this bias to allow statistical testing for the interplay of specific immune markers and the different parasitic infections, which could clarify the role of different parasites the related immune responses within the invasion process. Increasing the overall sample size and aligning the sample sizes of the compared populations could reveal more statistically significant differences in the variance and mean of immune effectors between native and invasive populations. Additionally, bigger sample sizes would allow to test for the interaction between sex and reproductive status in the context of immunity and invasion.
Second, the degree of genetic relatedness between the investigated invasive and native population is unknown. Invasive populations can have higher or lower genetic diversity than native conspecifics 58,59 . Variance in immunity might either compensate for potentially lower genotypic variation of invasive populations or reflect higher genetic diversity. Here we cannot distinguish if the observed variation is caused by one or the other. Studies of the genotypic diversity of the immune system (e.g. MHC and other immune genes) 60 of Egyptian geese from both populations could help to clarify the underlying mechanisms of the higher variance in immune effectors between the two populations in this study.
Third, local conditions (e.g. climate, habitats, community structures) of the two populations substantially differ. This differences in local conditions, which is a general aspect in invasion biology, clearly limits the interpretation of our results. Future studies of Egyptian geese populations (and those of other invasive species) should aim to increase the number of populations tested on a global scale to allow to control for local condition variables and their interplay with health indices of native and invasive populations 61 .

Conclusions
The interplay of parasites and host immunity in the frame of invasion biology is complex and still not clarified in detail. Its role likely depends on the characteristics of the parasites and the resulting costs of host's immune investment 9 . However, our study emphasizes that for unravelling the complexity of this interaction it could be very useful to measure both, mean and variance of different immune effectors in the light of parasitic infections. This could help to better understand contradictory findings in eco-immunological studies.
Contradictory effects of invasion on immune function have been reported for amphibians. Immune defences are weaker in invasive cane toads (Rhinella marina) that invade over long distances indicating a trade-off between dispersal and immune investment 22 . Cane toads on the invasion-front were found to have higher bactericidal and phagocytic activity than in more established populations 62 consistent with our observations in Egyptian geese. However, the overall immune investment in cane toads at the invasion front was lower compared to established cane toad populations 63 . In contrast, Cuban tree frogs (Osteopilus septentrionalis) on the invasion front exhibited reduced bactericidal ability compared to established frogs 64 . Studying flexibility of the immune system during www.nature.com/scientificreports/ range expansion of these amphibians might help to explain the contradictory findings. We expect that higher immune reaction flexibility would benefit individuals of invasive populations (1) in the early stages of invasion by reducing investment in the immune system and (2) in later stages of invasion by providing effective defence against novel pathogens. Intraspecific comparisons 21,65 will be essential to determine if flexible immune effector response is a general principle during biological invasions 36 allowing to allocate resources away from immunity whenever the parasitic burden is low. Common garden experiments combined with immune challenges or experimental infections would help to clarify the susceptibility to infection during range expansion of invasive species. Re-analysing data from previous studies based on the revised-EICA by investigating changes in the variance of immune effectors may also be a viable approach to determining if the results from Egyptian geese are relevant to other species. These approaches could help identify the drivers of successful invasion for neozootic species. Moreover, the interplay of local conditions and the strength of immune responses (variance) could add a new dimension to the field of invasion biology with predictive power regarding population dynamics of invasive populations in general.