Figures
Abstract
The ability to control infections is a key trait for migrants that must be balanced against other costly features of the migratory life. In this study we explored the links between migration and disease ecology by examining natural variation in parasite exposure and immunity in several populations of Lesser Black-backed Gulls (Larus fuscus) with different migratory strategies. We found higher activity of natural antibodies in long distance migrants from the nominate subspecies L.f.fuscus. Circulating levels of IgY showed large variation at the population level, while immune parameters associated with antimicrobial activity showed extensive variation at the individual level irrespective of population or migratory strategy. Pathogen prevalence showed large geographical variation. However, the seroprevalence of one of the gull-specific subtypes of avian influenza (H16) was associated to the migratory strategy, with lower prevalence among the long-distance migrants, suggesting that migration may play a role in disease dynamics of certain pathogens at the population level.
Citation: Arriero E, Müller I, Juvaste R, Martínez FJ, Bertolero A (2015) Variation in Immune Parameters and Disease Prevalence among Lesser Black-Backed Gulls (Larus fuscus sp.) with Different Migratory Strategies. PLoS ONE 10(2): e0118279. https://doi.org/10.1371/journal.pone.0118279
Academic Editor: Paul A. Bartell, Pennsylvania State University, UNITED STATES
Received: July 23, 2014; Accepted: January 12, 2015; Published: February 13, 2015
Copyright: © 2015 Arriero et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This study was supported by the Max Planck Institute for Ornithology (Wikelski group), an EU-funded SYNTHESYS grant to EA, and by the Spanish Ministry of Economy and Competitiveness (project BFU2011-25957 to EA). During writing, EA was supported by the Wissenschaftskolleg zu Berlin. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Migration has always attracted much attention as a fascinating life strategy involving displacements which can span several continents [1]. At the same time, migration still poses many questions regarding the considerable energetic demands and sophisticated physiological adaptations and mechanisms of navigation [1–3]. More recently, migration has attracted attention as an important factor presumably involved in spreading zoonotic diseases [4].
Migrants are exposed to a wide variety of environments, including all their biotic and abiotic components. Pathogens are an important biotic component of the environment and represent one of the major challenges for migrants [5,6]. Several phenomena imply that migratory animals may suffer higher impact from parasitism than non-migrants. Those include higher encounter rate with parasites from different areas (e.g. [7,8]), immunomodulation due to physiological demands incurred in migration [9,10], relapse of chronic infections prior or during the trip [11,12], or increased risk of infection due to host aggregation during migration or at the wintering sites [13].
However, despite the exposure to a broader range of pathogens and high susceptibility to infection presumably associated with long-distance migration, at the population level, the prevalence of infection may be lower in migrants than in residents, and some recent studies advocate for a role of migration in keeping populations healthy [14]. Several non-excluding hypotheses have been proposed to explain the role of migration in disease ecology. On one hand, from the perspective of the selection pressure imposed by parasitism, animals are expected to match their investment in immune defence to the parasitic threats they face, and in that context, migrant birds are expected to invest highly in immune defense [14,16]. Under this scenario, we would expect to find differences in parasite prevalence associated with migratory strategy and also associations between mechanisms of disease resistance and migration. At the interspecific level, a comparative study across 263 bird species showed lower prevalence of hemosporidian parasites in migrant species than in residents, although no explicit association between immunity and migration was tested in that study [18]. On the other hand, lower parasite prevalence in migrants has been proposed to be mediated by mechanisms such as migratory escape from infected areas [19], or migratory culling of sick individuals due to the exacerbated physiological demands of the trip and the infection [6]. Evidences from parasite prevalence in monarch butterflies support this role of migration in host-parasite dynamics [16,20]. This alternative scenario does not entail predictions about different immunity in migrants and residents. Thus, the main two mechanisms proposed underlying associations between migration and disease dynamics are elimination of infected individuals from the population, or selection for higher investment in mechanisms of disease resistance.
In the last decades, several animal species have shown changes in their migratory habits [21]. Changes in migratory strategy may be followed by changes in physiological adaptations and different host-pathogen dynamics, which may result in higher or lower disease risk [14,17]. In migratory birds with wide-ranging dietary preferences, the cost-benefit balance is shifted towards resident life when food is available year round in human-made settings like rubbish dumps [1]. Migration imposes tremendous physiological costs that may drive numerous trade-offs and physiological adaptations, including the investment in immune defence [6]. Therefore, changes towards a resident life style may alter for example migration-mediated immunomodulation [17,22].
On the other hand, variation in migratory strategy can result from the selection pressures that lead to species divergence [23]. When a taxonomic group diversifies in different migratory strategies, phenotypic or behavioural changes may be followed by changes in physiological adaptations [24].
The aim of this study is to examine variation in immune parameters and exposure to pathogens in a taxonomic group exhibiting wide variation in migratory strategies to test for associations between migration and parasite prevalence and immunity. The Lesser Black-backed Gull is a suitable model, because the taxonomic group includes several subspecies with different migratory habits, and gulls are a good example of recent changes in migratory behavior driven by human environmental modification. The Lesser Black-backed gull group (Larus fuscus) comprises three subspecies L.f. fuscus, L.f. intermedius and L.f. graellsii [25,26]. Although the subspecies grade into each other and have no clear limits, taxon designation is based on geographic location of the breeding colony and phenotypic differences [26]. The subspecies of Lesser Black-backed Gulls differ in dorsal plumage colour, ranging from black in the nominate taxa (from Northern Europe) to slate-grey in graellsii (from South-western Europe) [27]. The migratory strategy also varies within the Lesser Black-backed group: the nominate L.f. fuscus migrates long distances mainly to tropical Africa where it overwinters in large inland lakes, while intermedius and graellsii adopted a more resident life strategy and are medium or short distance migrants, wintering in South-western Europe and North-western Africa [1]. Despite the marked phenotypic differences in plumage color and migratory habits, molecular analyses show little genetic variation and poor phylogeographical structure within the Lesser Black-backed group [25].
Material and Methods
Study populations and sampling
The study was conducted during the breeding season of 2009 in five localities distributed along Europe, and including the only breeding colony in the Mediterranean. The five sampling sites include the subspecies graellsii (population from Northwest Spain (Sisargas)), intermedius (populations from the Netherlands (Moerdjik) and Northeast Spain (Ebro delta)), and fuscus (populations from Finland, Häme and Kokkola) (Fig. 1). Populations from Finland (and in general from Scandinavia) are considered long-distant migrants, as their main wintering sites are located in lakes from East Africa [28]. Lesser black-backed gulls from The Netherlands have been reported to winter predominantly along the coasts of South-western Europe, and are considered short-distance migrants [29]. Although there are, to our knowledge, no published data about the migratory strategy of the two colonies from the Iberian peninsula (colonies from Sisargas and Ebro Delta), ringing data, and especially observations of color rings outside the breeding season provide support for a short-distance migratory strategy in these colonies (20 out of the 25 gulls ringed in Sisargas for this study have been reported at least once wintering in the Mediterranean coast or Atlantic coast, and ringing data from Ebro delta showed that the Mediterranean area is the main wintering area for individuals of that colony, with only 6.3% of birds reported wintering in the Atlantic coast. In each locality ca. 25 adult breeding birds were captured at the nest with walk-in traps and blood samples and oropharyngeal and cloacal swabs were obtained according to standard sampling procedures. Heparinized whole blood samples were obtained from the wing vein (Vena ulnaris) and centrifuged (1000 G, 10 min) in order to obtain blood plasma. Plasma samples were stored at -20°C until analyses were performed. Swab samples (Virocult, Medical Wire and Equipment Co Ltd, Corsham, UK) were stored and shipped at 4°C and arrived in the laboratory within less than three days after sampling. All samples were obtained thanks to collaboration with local groups that were already monitoring breeding populations and catching adult birds. Procedures were approved by Finnish National Animal Experiment Board (ESLH-2009-03944/Ym-23). Sampling in the different locations was approved by local and regional authorities; Galicia: Dirección Xeral de Conservación da Naturaleza (Xunta de Galicia), Ebro Delta: Servei de Protecció i Gestió de la Fauna (Generalitat de Catalunya), The Netherlands: Vogeltrekstation, Finland: National Animal Experiment Board. Researchers involved in the design and sampling are certified in the use of animals for research purposes according to the current European legislation.
(Sisargas (Galicia, Spain); Delta (Ebro delta, Spain); Moerdjik (The Netherlands); Häme (Finland); Kokkola (Finland). Häme and Kokkola are long-distance migrants. Different symbols denote subspecies: circles L f.fuscus, triangles Lf. Intermedius, squares Lf.graellsii
Immune parameters
Natural antibodies and complement activity
The activity of natural antibodies and complement proteins were used as an estimation of constitutive immune defenses. Natural antibodies have been classified as constitutive components of both the innate and adaptive immune defense [30,31], while the activity of complement proteins is considered a constitutive innate defense [32]. Natural antibodies and complement cascade provide a first line of defense against pathogens. We estimated the activity of both components through the hemolysis-agglutination assay developed by Matson et al.[31]. In brief, plasma samples were serially diluted twofold with 0.01M phosphate buffered saline (PBS Sigma P3813) in 96-well round bottom plates (Corning Costar #3795) and incubated at 37°C for 90 min with a 1% rabbit blood cell suspension. A negative control of PBS was added in all the plates. Natural antibody (Nab) titers and complement activity were scored as –log2 of the highest dilution exhibiting agglutination (Nab) or hemolysis (complement).
Total IgY
Immunoglobulin levels were determined by a direct enzyme-linked immunosorbent assay (ELISA) optimized for this study following Martinez et al. [33]. Plates were coated with plasma diluted 1:4000 in 0.1M carbonate-bicarbonate buffer. A peroxidase-conjugated polyclonal rabbit anti-chicken IgY (Sigma A-9046) was used and colour reactions were developed by incubation with the substrate ABTS (Sigma A-1888). Absorbance was measured at 405nm with a plate reader (VersaMax ELISA Microplate Reader, Molecular Devices, Inc., Sunnyvale, CA, USA).
Haptoglobin concentration
Haptoglobin is an acute phase protein involved in binding free haem released into the circulation [34]. Haptoglobins play an important role as immunomodulators, and plasma circulating levels may increase several fold in the face of inflammatory stimulus, such as those caused by infections [34]. To estimate haptoglobin concentration in plasma we measured pseudoperoxidase activity of haemoglobin, which is directly proportional to the concentration of haptoglobin in the sample. We used a commercial haptoglobin test kit (TP-801, Tridelta Development Limited, Maynooth, Co. Kildare, Ireland), and plates were read at 630nm. OD values were ln-transformed before the analyses.
Lysozyme concentration
Plasma lysozyme levels were estimated following the assay adapted by Millet et al. [35]. We prepared a suspension of 50ml 1% agar (Difco, Detroit, MI) with 25mg dried Micrococcus lysodeikticus (Sigma M3770). Then 300μl of this suspension and 15μl of plasma sample from the birds were added in duplicates to a 96-well plate. We obtained a standard curve of lysozyme concentration by adding the bacterial suspension to serial dilutions of a standard lysozyme (Sigma L-6876) solution. Plates were incubated overnight at room temperature and absorbance was measured after 18hrs at 650nm with a microplate reader. Results were analyzed in ln-transformed OD units. Higher OD values indicate lower lysozyme activity.
Leucocyte profile
Leucocyte profile provides information on circulating immune cells and can be used as an indicator of current infections [36]. Relative numbers of the different leucocyte types were obtained examining blood smears stained with Giemsa under a light microscope at 1000x magnification. We obtained relative numbers of the different leucocyte types (heterophils, eosinophils, basophils, lymphocytes, monocytes) up to a total of 100 leucocytes. We calculated the H/L ratio as an indirect index of physiological stress [37].
Pathogen prevalence
Infections by blood parasites from the genus Plasmodium or Haemoproteus were assessed by amplifying a fraction (390bp) of the cytochrome-b gene of the parasite, using the primer pair Palu-F/Palu-R [38]. PCR products were visualized in 1.5% agarose gel and amplicons of positive samples were sequenced. Infections were confirmed by light microscopy with blood smears stained with Giemsa.
The infection status of viral infections was determined using swab and serum samples. Oropharyngeal and cloacal swabs were tested for avian influenza by real-time RT-PCR targeting the viral M-gene (Super-ScriptIII One-Step RT-PCR Kit with Platinum Taq polymerase, Invitrogen). To detect prior infections with avian influenza viruses (AIV) and avian paramyxoviruses (APMV) serum samples were analyzed by ELISA and hemagglutination inhibition assays (HI-test). Antibodies directed against the AIV nucleocapsid protein were detected using a competitive NP-ELISA (cELISA; ID Screen, Influenza A NP-Antibody Competition; ID.VET, Montpellier, France). The results were evaluated as percent inhibition (PI = 1 - (OD sample/mean OD negative) x 100). Samples with a PI <46% were classified as positive, samples with a PI>46% were classified as negative [39]. Positive samples were subsequently analysed by hemagglutination inhibition assays (HI-test) following standard procedures [39], and using a panel of four different reference antigens including the subtypes H5, H7, H13 and H16. All serum samples were additionally tested by HI-test for antibodies against three avian paramyxovirus serotypes APMV-1/Clone 30, APMV-4/duck/HongKong/D3/75 and APMV-6/duck/HongKong/199/77, following standard procedures [40].
Statistical analyses
Variation in immunological parameters was analyzed using linear mixed effects models, with the different immunological parameters as the dependent variable, and sex and infection status as fixed factors, and population, subspecies and migratory strategy as random factors with hierarchical structure (populations within subspecies and subspecies within migratory strategy). Models were subsequently simplified based on AIC values and likelihood ratio test of model comparison. Because several studies have shown sexual differences in immune parameters, [41], in the prevalence of Haemoproteus infections [42], or in the immune response to viral infections [43], sex was included in all statistical models as a fixed effect, and subsequently removed when the fit of the model was not improved. Variation in pathogen prevalence was analyzed with generalized linear mixed effects models with the infectious status (1/0) as the dependent variable with binomial distribution and the same hierarchical structure of random factors. A total of 106 samples were available for this study (48 males, 58 females), although not all the assays could be performed with all the samples due to technical issues, and therefore the sample size slightly varies in the statistical analyses. Means ± standard deviations are reported through the text. Analyses were performed using R 3.1.1 (http://www.R-project.org).
Results
Immune parameters
Migratory strategy explained 33% of the variation in the concentration of natural antibodies, with long distance migrants showing higher titers than short distance migrants (agglutination score long distance: 6.50 ± 1.20, short distance: 4.96±1.25, Fig. 2). Population identity explained 9% of the variation in natural antibodies, while subspecies explained less than 1% of the variation (Table 1). The best-fit model also included prevalence of Haemoproteus as a significant explanatory variable of individual variation in natural antibodies, with infected birds showing lower levels than uninfected birds (t = -3.181, P = 0.002, non-infected: 5.66 ± 1.35, infected: 4.61±1.19). Total IgY levels showed significant variation among populations (19% of variance explained), but migratory strategy explained less than 1% of the variance (Table 1). The model included sex (t = 1.752, P = 0.083), Haemoproteus prevalence (t = 3.247, P = 0.002) and seroprevalence of avian influenza (t = 2.312, P = 0.023) as explanatory variables, with males and birds which tested positive for the parasitic infection showing higher IgY concentrations. We also found a significant interaction between sex and prevalence of Haemoproteus (t = -2.126, P = 0.036), showing that in females, but not in males there was a significant increase in IgY levels associated with infection. The interaction between avian influenza and Haemoproteus prevalence was also significant (t = -3.622, P < 0.001), showing a much higher increase in IgY levels when infected by haemoparasites among birds which showed no previous exposure to avian influenza.
Mean, SE and SD.
Immune parameters mainly correlated with antimicrobial activity such as complement and lysozyme activity did not show significant differences among populations, and migratory strategy explained less than 1% of the variation observed in these parameters (Table 1). We did not find association between pathogen prevalence and complement or lysozyme activity (i.e. none of the variables of pathogen prevalence were retained in the final models). However, lysozyme activity was higher in females than in males (t = -2.695, P = 0.008; females: 1.08 ± 0.14, males: 1.02 ± 0.07).
Haptoglobin levels were not associated with migratory strategy, but 8% of the variation was explained by differences among populations (Table 1). Birds suffering from Haemoproteus sp. infections had lower levels of haptoglobin (t = -2.055, P = 0.043, infected: 0.08 ± 0.07; non-infected: 0.13 ± 0.14). Leucocyte profiles and H/L ratio did not show any pattern with migratory strategy or subspecies, and we did not find significant variation among populations (Table 1). No association between leucocyte profiles and infectious status was found.
Natural antibody titers were positively and significantly correlated with complement activity (t = 5.139, P < 0.001). Lysozyme activity was negatively correlated with IgY (t = -2.071, P = 0.041).
Pathogen prevalence
Total prevalence of protozoan blood parasites was 34%. Haemoparasites were identified as belonging to the genus Haemoproteus, except for two individuals infected by Parahaemoproteus. Variation in Haemoproteus prevalence was explained by a model that included population id, but migratory strategy and subspecies explained less than 1% of the variance and were subsequently removed as explanatory factors (Table 2). The other variables retained in the best-fit model were the activity of natural antibodies and plasma concentrations of haptoglobin, with infected birds showing lower concentration of natural antibodies (z = -2.557, P = 0.011), and a trend for lower concentration of haptoglobin (z = -1.817, P = 0.069; see mean values above).
Exposure to the different viral infections tested varied greatly among populations (Table 2). All samples were tested negative for the presence of AIV-RNA with RT-PCR; ongoing infections with influenza A viruses at the time of sampling could thus be excluded. Overall seroprevalence of avian influenza (AI) was high (68.3%), with significant variation among populations (z = 3.481, P< 0.001, Table 2). Migratory strategy and subspecies did not explain differences in seroprevalence and were not retained in the final model. However, we found a significantly lower proportion of males than females with anti-influenza antibodies (z = -2.66, P = 0.007, males: 0.54, females: 0.83). H16 was the most common AI-subtype detected, and its prevalence among the AI seropositives was explained by migratory strategy (z = 3.009, P = 0.003), with long distance migrant populations showing lower prevalence (Table 2, Fig. 3). Seroprevalence of H13 and H5 did not differ significantly among populations, and H7 was only detected in an individual from Islas Sisargas (Spain). Presence of antibodies against several hemagglutinin subtypes in the same individual was relatively common, but we did not find significant differences among populations.
Mean, SE and SD.
Seroprevalence of Avian paramyxovirus serotypes differed significantly among populations (Table 2). The proportion of individuals with detected antibodies against Newcastle Disease Virus (APMV-1) was 45%, and population id was the only explanatory factor retained in the statistical model (z = 1.989, P = 0.045). Total seroprevalence of APMV-6 was 30%, also with significant differences among populations (Table 2). The best-fit model for APMV-6 included population id (z = -1.743, P = 0.081) and plasma levels of IgY (z = 1.757, P = 0.079) as explanatory variables, with slightly higher IgY levels in birds seropositive to APMV-6 (negative: 2.65, positive: 2.76). Subspecies and migratory strategy explained less than 1% of the variance in seroprevalence of APMV-6 and APMV-1, and were not retained in the final models. We did not detect exposure to APMV-4. In total, 17% of the birds had antibodies against multiple viral strains, but there were no significant differences among populations (P = 0.65).
Discussion
Our study shows significant geographical variation in pathogen prevalence and immune parameters in Lesser Black-backed Gull populations with different migratory strategies. We found higher activity of natural antibodies in gulls of the long-distance migrant nominate subspecies Larus fuscus fuscus, suggesting that long-distance migrants maintain higher constitutive immunity than short-distance migrants.
Since the end of the 60’s the Western populations of Lesser Black-backed Gulls adopted an omnivorous diet and a resident life style [1]. Avoiding tropical areas and less energetically demanding migratory habits [29] seem to be associated with lower investment in constitutive antibody-mediated immunity in the subspecies intermedius and graellsii. Several studies have highlighted the important role of the environment in explaining geographical variation in immune parameters at intra- or interspecific level [44–47].
The taxonomic group constituted by Larus fuscus has been reported to have little genetic variation [48]. Our study shows little differentiation in physiological parameters involved in immune defense among subspecies, which suggests a minor role of genetic effects explaining the variation found in immunity and disease susceptibility in our study [49,50]. The different subspecies of Larus fuscus have undergone some degree of phenotypic differentiation in morphological characters, (e.g. plumage coloration) and in migratory habits [51]. However, only the nominate subspecies L. fuscus fuscus exhibits differences in physiological parameters involved in the immune response. This may have important implications in the susceptibility to pathogens and the ability to control infections [52], although not all pathogens had lower prevalence in the Finnish population, and variation in pathogen prevalence was not explained at the taxonomic level of subspecies.
We found large geographical differences in the prevalence of some of the major avian diseases in Lesser Black-backed Gulls. Furthermore, short-distance migratory strategy was associated with higher seroprevalence of one of the gull-specific subtypes of avian influenza, H16. Although we do not have serologic evidence of exposure to highly pathogenic strains of avian influenza virus in the gulls of our study, infections with low pathogenic strains of avian influenza have been reported to impact migratory performance in other species [53,54]. Thus, infections with low pathogenic strains of avian influenza may have higher clinical and ecological consequences than previously thought, or may have exacerbated effects due to other concurrent infections, highlighting migratory culling as a likely explanation for the lower seroprevalence of H16 found in long-distance migrants. Further support for migratory culling is provided by the fact that seroprevalence indicates the proportion of individuals which have been infected and survived an infection by that particular serotype of avian influenza.
Our results show that over 50% of the birds tested had been exposed to different hemagglutinin subtypes of AIV, and a large proportion of them had also been exposed to several avian paramyxovirus serotypes, including New Castle Disease Virus (APMV-1). Prolonged use of habitats may explain high prevalence of pathogens with oral-fecal mode of transmission [55,56]. This may explain the important spatial variation found in our study in exposure to avian viruses, and the remarkably high seroprevalence in some of the colonies tested (Table 2).
Anseriformes (ducks, geese, swans) and Charadriiformes (gulls, terns, shorebirds) are considered the main natural reservoirs for avian influenza viruses, and important vectors of several other infectious diseases (e.g.[57–59]). In particular, gulls are considered key in facilitating intercontinental exchanges of AIVs [60] and may have a relevant epidemiological role because of their close contact with human settlements [61].
Despite the fact that none of the birds tested in our study were shedding influenza viruses at the time of sampling, we found high seroprevalence rates (52–88%), indicating that a large proportion of the gulls had been into contact with the virus and had mounted an immune response [62]. Our results are comparable to previous results obtained in other gull species (e.g. [63,64]), and indicate that the Lesser Black backed Gulls of the five populations sampled have been exposed to several hemagglutinin subtypes, with H13 and H16 being the most prevalent. H13 and H16 are considered gull-specific subtypes [65–67], and their maintenance in gull populations is attributed to high incidence of LPAIV among first year birds [68]. In our study, adult gulls were sampled in the early breeding season (i.e. before hatching), suggesting that the high seroprevalence of H13 and H16 is an indication of perpetuation of these serotypes within the populations.
Interestingly, Lesser Black backed Gulls had antibodies against the hemagglutinin subtypes H5 and H7, which are potential precursors of highly pathogenic avian influenza and are considered notifiable diseases with special regulatory measures in most countries. However, while individuals from all colonies sampled were seropositive for H5 (with seroprevalence ranging between 13–47%), only a single individual from the colony in North-Western Spain tested positive for antibodies against H7.
We found an overall seroprevalence of APMV-1 of 45%. Although velogenic strains of APVM-1 (NDV) are the causation agent of Newcastle Disease, one of the most devastating diseases of poultry, the strains of the virus commonly found in waterfowl and sea birds are considered avirulent [55]. These strains grow mainly in the intestinal tract and cause loss of appetite and dehydration [69]. However, for migrating birds, digestive problems may have severe consequences, by depleting the required body stores for the trip [6], or impairing immune capacity [70].
Despite the previous general perception of very low incidence of haematozoan infections in marine and coastal birds [12,71,72], data provided in recent studies suggest that the incidence of haematozoan infections in marine birds is higher than previously thought (rev. Quillfeldt et al. [73]). Our results show an overall haematozoan prevalence of 34% in gulls of the Larus fuscus taxonomic group. Infections by avian haematozoa are known to impair several physiological functions and significantly affect flying performance [12,74]. Therefore, it seems likely that migratory culling plays a role in maintaining low prevalence of infections among some of the long-distance migratory populations [14]. However, our results show extensive geographical variation in the prevalence of hemoparasites, not associated to migratory habits.
Conclusions
Our study explored the links between animal migration and disease risk, and found large differences among populations in parameters of anti-parasite defence and in the prevalence of avian pathogenic diseases. Our results offer partial support for a role of migration in disease dynamics of Lesser black-backed gulls. We found higher activity of constitutive immunity in long-distance migrants, which is in line with previous studies suggesting that long-distance migration selects for higher investment in immunity [15]. In addition, we found higher seroprevalence of one of the most common serotypes of avian influenza in populations with short-distance migratory strategy, supporting the idea that adopting a resident or short-distance migratory strategy may result in higher disease prevalence at the population level [14]. However, we found large geographical variation in other immune parameters and prevalence of other viral and protozoan infections, suggesting important variation in disease dynamics at the population level not explained by differences in migratory strategy.
Supporting Information
S1 Table. Data set used in the study.
Delta = Delta Ebro (Spain); Häme = Finland; Kokkola = Finland; Moerdjik = The Netherlands; Sisargas = Galicia (Spain); Agglu = agglutination score; Hp = haptoglobin concentration; Lysozyme = lysozyme activity (OD); H = heterophils; L = lymphocytes; E = eosinophils; B = basophils; M = monocytes; Haemoproteus = infection status; FluA = prevalence of antibodies against influenza A; PMV1–6 = prevalence of antibodies against avian paramixovirus 1, 4 or 6.
https://doi.org/10.1371/journal.pone.0118279.s001
(DOCX)
Acknowledgments
This study was possible thanks to the collaboration of several groups who provided samples or logistics to access the breeding colonies. We thank T. Muuse, A. Velando, J. Morales and C. Perez, for assistance in the field, E. Fricke, C. Illing, J. A. Tendero and C. Remacha for helping with laboratory assays. S. Merino (MNCN-CSIC) provided support and advice on haemoparasite infections, and T. Harder and C. Grund (Friedrich- Loeffler-Institute) on serological and virological analyses.
Author Contributions
Conceived and designed the experiments: EA. Performed the experiments: EA IM RJ FJM AB. Analyzed the data: EA IM. Contributed reagents/materials/analysis tools: EA IM RJ FJM AB. Wrote the paper: EA IM RJ FJM AB.
References
- 1.
Alerstam T. Bird Migration. Cambridge University Press; 1993. pmid:25590127
- 2. Wikelski M, Tarlow EM, Raim A, Diehl RH, Larkin RP, Visser GH. Costs of migration in free-flying songbirds. Nature 2003; 423: 704. pmid:12802324
- 3.
McWilliams SR, Karasov WH. Migration takes guts. In: Greenberg R and Marra P, editors. Birds of two worlds: the ecology and evolution of migration. Smithsonian Institution Press, Washington, D.C.; 2005. pp. 67–78.
- 4. Altizer S, Bartel R, Han BA. Animal migration and infectious disease risk. Science 2011; 331: 296–302. pmid:21252339
- 5.
Ricklefs RE, Fallon SM, Latta SC. Migrants and their parasites. In: Greenberg R and Marra P, editors. Birds of two worlds: the ecology and evolution of migration. Smithsonian Institution Press, Washington, D.C.; 2005. pp. 210–221.
- 6. Klaassen M, Hoye BJ, Nolet BA, Buttemer WA. Ecophysiology of avian migration in the face of current global hazards. Philos Trans R Soc Lond B Biol Sci 2012; 367: 1719–1732. pmid:22566678
- 7. Waldenstrom J, Bensch S, Kiboi S, Hasselquist D, Ottosson U. Cross-species infection of blood parasites between resident and migratory songbirds in Africa. Mol Ecol 2002; 11: 1545–1554. pmid:12144673
- 8. Ricklefs RE, Fallon SM. Diversification and host switching in avian malaria parasites. Proc Biol Sci 2002; 269: 885–892. pmid:12028770
- 9. Owen J, Moore F. Seasonal differences in immunological condition of three species of thrushes. Condor 2006; 108: 389–398.
- 10. Buehler DM, Piersma T, Matson K, Tieleman BI. Seasonal redistribution of immune function in a migrant shorebird: annual‐cycle effects override adjustments to thermal regime. Am Nat 2008; 172: 783–796. pmid:18999941
- 11. Gylfe A, Bergström S, Lundström J, Olsen B. Reactivation of Borrelia infection in birds. Nature 2000; 403: 724–725. pmid:10693792
- 12.
Valkiunas G. Avian malaria parasites and other haemosporidia. CRC press; 2004. pmid:25057686
- 13. Krauss S, Stallknecht DE, Negovetich NJ, Niles LJ, Webby RJ, et al. Coincident ruddy turnstone migration and horseshoe crab spawning creates an ecological “hot spot” for influenza viruses. Proc Biol Sci 2010; 277: 3373–3379. pmid:20630885
- 14. Altizer S, Bartel R, Han BA. Animal migration and infectious disease risk. Science 2011; 331: 296–302. pmid:21252339
- 15. Møller A, Erritzøe J. Host immune defence and migration in birds. Evol Ecol 1998; 12: 945–953.
- 16. Altizer SM. Migratory behaviour and host—parasite co-evolution in natural populations of monarch butterflies infected with a protozoan parasite. Evol Ecol Res 2001; 3: 611–632.
- 17. Buehler DM, Tieleman BI, Piersma T. How do migratory species stay healthy over the annual cycle? A conceptual model for immune function and for resistance to disease. Integr Comp Biol 2010; 50: 346–357. pmid:21558209
- 18. Arriero E, Møller AP. Host ecology and life-history traits associated with blood parasite species richness in birds. J Evol Biol 2008; 21:1504–1513. pmid:18823451
- 19. Loehle C. Social barriers to pathogen transmission in wild animal populations. Ecology 1995; 76: 326–335.
- 20. Bradley CA, Altizer S. Parasites hinder monarch butterfly flight: implications for disease spread in migratory hosts. Ecol Lett 2005; 8: 290–300.
- 21.
Wilcove D. No way home: the decline of the world’s great animal migrations. Washington D.C: Island Press; 2008. pmid:25506952
- 22. Owen JFM. Swainson’s thrushes in migratory disposition exhibit reduced immune function. J Ethol 2008; 26: 383–388.
- 23.
Mayr E. Animal species and their evolution. Harvard University Press; 1963. pmid:14191062
- 24. Garland TJ, Carter PA. Evolutionary physiology. Annu Rev Physiol 1994; 56: 579–621. pmid:8010752
- 25. Liebers D, Helbig AJ. Phylogeography and colonization history of Lesser Black-backed Gulls (Larus fuscus) as revealed by mtDNA sequences. J Evol Biol 2002; 15: 1021–1033. .
- 26. Liebers D, de Knijff P, Helbig AJ. The herring gull complex is not a ring species. Proc Biol Sci 2004; 271: 893–901. pmid:15255043
- 27.
Malling OK, Larsson H. Gulls of Europe, Asia and North America. Helm Identification guides B, London; 2003. pmid:25057689
- 28. Bustnes JO, Moe B, Helberg M, Phillips RA. Short communication Rapid long-distance migration in Norwegian Lesser Black-backed Gulls Larus fuscus fuscus along their eastern flyway. Ibis 2013; 155: 402–406.
- 29. Klaassen RHG, Ens BJ, Shamoun-Baranes J, Exo KM, Bairlein F. Migration strategy of a flight generalist, the Lesser Black-backed Gull Larus fuscus. Behav Ecol 2011; 23: 58–68.
- 30. Ochsenbein A, Zinkernagel R. Natural antibodies and complement link innate and acquired immunity. Immunol Today 2000; 5699: 624–630.
- 31. Matson KD, Ricklefs RE, Klasing KC. A hemolysis-hemagglutination assay for characterizing constitutive innate humoral immunity in wild and domestic birds. Dev Comp Immunol 2005; 29: 275–286. pmid:15572075
- 32. Lee KA. Linking immune defenses and life history at the levels of the individual and the species. Integr Comp Biol 2006; 46: 1000–1015. pmid:21672803
- 33. Martinez J, Tomás G, Merino S, Arriero E, Moreno J. Detection of serum immunoglobulins in wild birds by direct ELISA: a methodological study to validate the technique in different species using antichicken antibodies. Funct Ecol 2003; 17: 700–706.
- 34. Quaye IK. Haptoglobin, inflammation and disease. Trans R Soc Trop Med Hyg 2008; 102: 735–742. pmid:18486167
- 35. Millet S, Bennett J, Lee KA, Hau M, Klasing KC. Quantifying and comparing constitutive immunity across avian species. Dev Comp Immunol 2007; 31: 188–201. pmid:16870251
- 36.
Campbell T, Ellis CK. Avian and exotic animal hematology and cytology. 3rd ed. Ames: Willey, Blackwell; 2007. https://doi.org/10.1093/jxb/erm028 pmid:25506957
- 37. Gross WB, Siegel HS. Evaluation of the heterophil/lymphocyte ratio as a measure of stress in chickens. Avian Dis 1982; 27: 972–979.
- 38. Martínez J, Martínez-de-la-Puente J, Herrero J, Del Cerro S, Lobato E, et al. A restriction site to differentiate Plasmodium and Haemoproteus infections in birds: on the inefficiency of general primers for detection of mixed infections. Parasitology 2009; 136: 713–722. pmid:19416554
- 39.
CEC. Commission decision 2006/437/EC Aproving a diagnostic manual for avian influenza as provided for in Council Directive 2005/94/EC. Off J Eur Communities; 2006.
- 40.
O.I.E. Newcastle disease. Manual of diagnostic tests and vaccines for terrestrial animals. Paris: World organisation for animal health, O.I.E. 2008; pp. 576–589. https://doi.org/10.1111/j.1440-1681.2008.05117.x pmid:19704473
- 41.
O’Neal DM, Keterson ED. Life-history evolution, hormones, and avian immune function. In: Demas GE, Nelson R, editor. Ecoimmunology. Oxford University Press. 2012; pp.7–44.
- 42. McCurdy DG, Shutler D, Mullie A, Forbes MR. Sex-biased parasitism of avian hosts: relations to blood parasite taxon and mating system. Oikos 1998; 82: 303–312.
- 43.
Robinson DP, Klein SL. Sex differences in immune response to viruses. In: Demas G.E., Nelson R, editor. Immunoecology. Oxford University Press. 2012; pp. 497–529.
- 44. Grasman KA, Armstrong M, Hammersley DL, Scanlon PF, Fox GA. Geographic variation in blood plasma protein concentrations of young herring gulls (Larus argentatus) and Caspian terns (Sterna caspia) from the Great Lakes and Lake Winnipeg. Comp Biochem Physiol C Toxicol Pharmacol 2000; 125: 365–375. pmid:11790357
- 45. Ardia D. The ability to mount multiple immune responses simultaneously varies across the range of the tree swallow. Ecography 2007; 30: 23–30.
- 46. Barbosa A, Merino S, Benzal J, Martinez J, García-Fraile S. Geographic variation in the immunoglobulin levels in pygoscelid penguins. Polar Biol 2006; 30: 219–225.
- 47. Buehler DM, Tieleman IB, Piersma T. Age and environment affect constitutive immune function in Red Knots (Calidris canutus). J Ornithol 2009; 150: 815–825.
- 48. Liebers D, Helbig AJ. Phylogeography and colonization history of Lesser Black-backed Gulls (Larus fuscus) as revealed by mtDNA sequences. J Evol Biol 2002;15: 1021–1033.
- 49. Bonneaud C, Pérez-Tris J, Federici P, Chastel O, Sorci G. Major histocompatibility alleles associated with local resistance to malaria in a passerine. Evolution 2006; 60: 383–389. pmid:16610328
- 50. Turner AK, Begon M, Jackson J A, Bradley J E, Paterson S. Genetic diversity in cytokines associated with immune variation and resistance to multiple pathogens in a natural rodent population. PLoS Genet 2011; 7: e1002343. pmid:22039363
- 51.
Cramp S. Handbook of the birds of Europe, the Middle East and North Africa. Oxford: Oxford University Press; 1983. pmid:25121320
- 52.
Murphy K, Travers P, Walport M. Janeway’s Immunobiology. 7th ed. New York: Garland Science; 2008. https://doi.org/10.1080/13803390802438454 pmid:19044189
- 53. Van Gils JA, Munster VJ, Radersma R, Liefhebber D, Fouchier RAM, et al. Hampered foraging and migratory performance in swans infected with low-pathogenic avian influenza A virus. PLoS One 2007; 2: e184. pmid:17264886
- 54. Latorre-Margalef N, Gunnarsson G, Munster VJ, Fouchier RAM, Osterhaus ADME, et al. Effects of influenza A virus infection on migrating mallard ducks. Proc Biol Sci 2009; 276: 1029–1036. pmid:19129127
- 55.
Alexander D. Newcastle disease. Kluwer Academic Publisher. Boston, MA, USA; 1988. pmid:25144059
- 56. Munster VJ, Baas C, Lexmond P, Waldenström J, Wallensten A, et al. Spatial, temporal, and species variation in prevalence of influenza A viruses in wild migratory birds. PLoS Pathog 2007; 3: e61. pmid:17500589
- 57. Hubalek Z. An annotated checklist of pathogenic microorganisms. J Wildl Dis 2004; 40: 639–659. pmid:15650082
- 58. Olsen B, Munster VJ, Wallensten A, Waldenström J, Osterhaus ADME, et al. Global patterns of influenza a virus in wild birds. Science 2006; 312: 384–388. pmid:16627734
- 59. Krauss S, Obert CA, Franks J, Walker D, Jones K, et al. Influenza in migratory birds and evidence of limited intercontinental virus exchange. PLoS Pathog 2007; 3: e167. pmid:17997603
- 60. Wille M, Robertson GJ, Whitney H, Ojkic D, Lang AS. Reassortment of American and Eurasian genes in an influenza A virus isolated from a great black-backed gull (Larus marinus), a species demonstrated to move between these regions. Arch Virol 2011; 156: 107–115. pmid:21053031
- 61. Arnal A, Vittecoq M, Pearce-Duvet J, Gauthier-Clerc M, Boulinier T, et al. Laridae: A neglected reservoir that could play a major role in avian influenza virus epidemiological dynamics. Crit Rev Microbiol 2014; 7828: 1–12.
- 62. Brown JD, Luttrell MP, Berghaus RD, Kistler W, Keeler SP, et al. Prevalence of antibodies to type a influenza virus in wild avian species using two serologic assays. J Wildl Dis 2010; 46: 896–911. pmid:20688695
- 63. Lewis NS, Javakhishvili Z, Russell CA, Machablishvili A, Lexmond P, et al. Avian influenza virus surveillance in wild birds in Georgia: 2009–2011. PLoS One 2013; 8: e58534. pmid:23516501
- 64. Verhagen JH, Munster VJ, Majoor F, Lexmond P, Vuong O, et al. Avian influenza a virus in wild birds in highly urbanized areas. PLoS One 2012; 7: e38256. pmid:22761671
- 65. Hinshaw VS, Air GM, Gibbs AJ, Graves L, Prescott B, Karunakaran D. Antigenic and genetic characterization of a novel hemagglutinin subtype of influenza A viruses from gulls. J Virol 1982; 42: 865–872. pmid:7097861
- 66. Fouchier RAM, Munster V, Wallensten A, Bestebroer TM, Herfst S, et al. Characterization of a novel influenza A virus hemagglutinin subtype (H16) obtained from black-headed gulls. J Virol 2005; 79: 2814–2822. pmid:15709000
- 67. Obenauer JC, Denson J, Mehta PK, Su X, Mukatira S, et al. Large-scale sequence analysis of avian influenza isolates. Science 2006; 311: 1576–1580. pmid:16439620
- 68. Verhagen JH, Majoor F, Lexmond P, Vuong O, Kasemir G, et al. Epidemiology of influenza A virus among Black- headed gulls, the Netherlands 2006–2010. Emerg Infect Dis 2014; 20: 2008–2011.
- 69.
McFerran J, McCracken R. Newcastle disease. Kluwer Academic Publisher; 1988. pmid:25144059
- 70. Klasing KC. Nutritional modulation of resistance to infectious diseases. Poult Sci 1998; 77: 1119–1125. pmid:9706075
- 71. Mendes L, Piersma T, Lecoq M, Spaans B, Ricklefs R. Disease‐limited distributions? Contrasts in the prevalence of avian malaria in shorebird species using marine and freshwater habitats. Oikos 2005; 109: 396–404.
- 72. Yohannes E, Križanauskienė A, Valcu M, Bensch S, Kempenaers B. Prevalence of malaria and related haemosporidian parasites in two shorebird species with different winter habitat distribution. J Ornithol 2009; 150: 287–291.
- 73. Quillfeldt P, Arriero E, Martínez J, Masello JF, Merino S. Prevalence of blood parasites in seabirds—a review. Front Zool 2011; 8: 26. pmid:22035144
- 74.
Atkinson C. Haemoproteus. In: Atkinson C, Thomas N, Hunter B, editors. Parasitic Diseases of Wild Birds. Wiley-Blackwell. 2008; pp. 13–34.