The Young, the Weak and the Sick: Evidence of Natural Selection by Predation

It is assumed that predators mainly prey on substandard individuals, but even though some studies partially support this idea, evidence with large sample sizes, exhaustive analysis of prey and robust analysis is lacking. We gathered data from a culling program of yellow-legged gulls killed by two methods: by the use of raptors or by shooting at random. We compared both data sets to assess whether birds of prey killed randomly or by relying on specific individual features of the prey. We carried out a meticulous post-mortem examination of individuals, and analysing multiple prey characteristics simultaneously we show that raptors did not hunt randomly, but rather preferentially predate on juveniles, sick gulls, and individuals with poor muscle condition. Strikingly, gulls with an unusually good muscle condition were also predated more than expected, supporting the mass-dependent predation risk theory. This article provides a reliable example of how natural selection may operate in the wild and proves that predators mainly prey on substandard individuals.


Introduction
Predation is an important selective force in evolution [1][2][3][4][5][6][7] and is generally assumed to select against substandard individuals, i.e. the young, senescent, sick, or individuals in poor physical condition [8][9]. Although some studies support this hypothesis and have contributed substantially to the understanding of selection by predation [10][11][12][13][14][15][16][17][18][19][20], (but see [18]), most of them are based on opportunistic observations or rely on some specific traits of the prey [15,[18][19][20] or parasite load [12,13]. Kenward [10] and Temple [11] investigated morphological traits of individuals together with their healthy state, however, sample sizes were rather small (less than 30 prey) and the effects of different traits were not analysed simultaneously, so the contribution of each trait on differential predation it is difficult to evaluate. Additionally, all evidence typically comes from the typical predator-prey system, where traits may have coevolved in parallel, and thus predation upon substandard individuals could be an opportunistic foraging strategy rather than a response to substandard features of the prey. To fully understand the role of predation as a selective force, it is also necessary to collect evidence of predation outside the typical predator-prey system, gather information on large sample sizes, and investigate multiple traits of prey simultaneously.
Populations of large gulls have increased substantially over the last century and some species are currently perceived as a pest by wildlife managers [21][22][23][24] but see [25]. As a consequence, many conservation agencies have set up culling programs to control gull populations, which typically consist of the systematic removal of large numbers of eggs, chicks or breeding adults. Even if the efficacy of these culling programs is still under debate [25][26][27][28], many conservation agencies still control gulls by culling. The Local Government of the Balearic Islands (Spain) began a gull culling programme on a refuse tip in the island of Mallorca as a part of the population control of yellow-legged gulls -Larus michahellisin the Balearic archipelago. From 2003 to 2007 birds were culled by two methods: by shooting or by the use of trained birds of prey (peregrine falcon -Falco peregrinus-, saker falcon -F. cherrugand Harris's hawk -Parabuteo unicinctus-). We gathered data from this culling program and examined killed birds to determine 1) the sex and age of the individual, 2) individual body condition, assessed from muscle condition, and 3) any sign of parasitism (internal and external), infection, malformation or chronic disease (e.g. aspergillosis). We used these data to investigate multiple prey traits simultaneously and to assess whether birds of prey killed randomly or by relying on specific individual prey features.

Results
We examined 506 gulls that had been shot and 122 gulls removed by raptors. The age structure in the shooting sample was similar to that observed at the dump over five available censuses (Breslow-Day homogeneity test of odds ratio, x 2 4~3 .311, P = 0.507; Mantel-Haenszel, odds-ratio -log transformed-95% Confidence Intervals: [20.066; 0.102], see Table 1). Post-mortem examination of gulls (Table 2) showed a low prevalence of external parasites; however, half of the individuals examined had internal parasites, mainly cestodes. About 7% of the gulls showed infection by Salmonella and 4% by Aspergillus. A few individuals also showed some kind of congenital malformation (bill deformity) and others had alteration of internal organs ( Table 2), but prevalence of most veterinarian findings were low. Individuals that had been shot showed different characteristics than those killed by raptors ( Table 3). The Multiple Component Analysis (MCA) scores plot also showed differences between groups, with more healthy adults and a higher average muscle condition within the group of individuals shot, and more juveniles, gulls in poor condition or showing some signs of illness in the group of individuals killed by birds of prey (Fig. 1). To test for the significance of these differences we used logistic regressions with predation by raptors as the response variable, being ''killed by raptor'' = 1 and ''killed by shooting'' = 0. In this way we retained shooting as the intercept of the regression to check for differential predation. The overdispersion value of the saturated model was 1.09, indicating a good fit of the data. The best ranked model (based on AIC values, Table S1) included a negative effect of age and a positive effect of both sickness and poor muscle condition on the probability of being predated by birds of prey (Table 4); this model did not include a gender effect. Strikingly, individuals in unusually good condition were also predated more frequently than expected by chance (Table 4). A model including an effect of sex explained the data equally well and was statistically equivalent to the previous model (Table S1) but the effect was not significant (z = 21.293, P = 0.196). Note that all models with lower AIC values unequivocally showed that age, muscle condition and sickness were clues for differential predation by birds of prey (see also Fig. 2). When these three variables were tested separately, results showed that muscle condition was the main factor affecting predation, this variable alone explaining 71% of the total variance.

Discussion
Natural selection of certain prey traits (e.g. morphological traits) has repeatedly been shown to be driven by predation [10][11][12][13][14][15][16][17][18][19][20]. Our paper could not address such particular issue, but on the other hand, the exhaustive analysis of a large number of prey, combined with the simultaneous analysis of a variety of traits, give us some general insights into how predation may operate in the wild. Here we show that predators did not kill individuals at random, but rather selected their prey on the basis of several, not always related traits. Our results indicated that age, muscle condition and sickness influence the probability of being predated, with juveniles, sick gulls, and individuals with poor muscle condition being killed  preferentially; thus strongly supporting the hypothesis that predators prey primarily on substandard individuals [8,9]. Natural selection acts on many characters simultaneously [29] and accordingly, our results dealing with natural selection by predation would suggest that two or more traits may affect fitness in an interactive way (i.e. correlational selection) [7,29,30]. Nevertheless, even if our sample sizes were relatively large, they still lack sufficient power to simultaneously test interactive effects between many characters.
Our study also shows that not only individuals with severe diseases but also those with mild diseases are predated preferentially, indicating that subtle changes in behaviour or condition may have been sufficient to increase susceptibility to predation. This was also found by Miller et al (2000) who showed that prion infection in deer increased the rate of predation of deer by mountain lions (Puma concolor) nearly fourfold, even if few of the deer killed were recorded as ''noticeably ill'' by field observers prior to their deaths [31].
Our results clearly support the mass-dependent predation risk (MDPR) theory, which predicts that birds should keep their mass as low as possible to reduce their likelihood of being killed by predators [32,33]. To date, empirical evidence for this theory comes only from small passerines with body masses between 10-150gr [34,35 and references therein] and recently from one mammal [36]. Additionally, results showed preferential predation on those individuals with poor muscle condition, suggesting that stabilizing selection [37] could be operating on traits linked to body mass.
This article provides a reliable, robust example of how natural selection by predation operates in the wild and strongly supports the paradigm that predators kill substandard individuals. Since gulls are an occasional prey of falcons and hawks, results probably indicate the ability of predators to detect substandard individuals in the wild rather than showing an optimal foraging strategy or possible coevolution within a natural predator-prey system.

Materials and Methods
From 2003 to 2007, we carried out 5 gull censuses at the landfill and estimated the proportion of birds in each age class. We identified four age classes by plumage features: juveniles (from 0 to 1 year old), 2-year-olds (from 1 to 2 years old), 3-year-olds (from 2 to 3 years old), and adults (.3 years old); these data were used to assess whether shooting was performed randomly regarding to age, by means of a goodness-of-fit test. We cannot exclude that shooting was biased in relation to veterinarian findings (health state) and muscular condition; however if a bias existed it would rend our comparison more conservative as individuals that were ill or in poorer condition should be shot preferentially [38].
Post-mortem examination of culled individuals was carried out at the Fundació Natura Parc-COFIB wildlife recovery centre in Removed birds were examined to determine the sex and age of the individual, the individual nutritional state, assessed from fat layers and muscle condition, and any sign of infection, malformation or disease. We identified four age classes by plumage features: juveniles (from 0 to 1 year old), 2 years old (from 1 to 2 years old), 3 years old (from 2 to 3 years old), and adults (.3 years old). We determined three levels of body condition: normal, low and high, depending on the layers of muscle mass). Veterinarians determined if the illness detected was either severe or mild. doi:10.1371/journal.pone.0009774.t003 Mallorca. Within each age-by-sex class, we sorted individuals as healthy, mildly sick or severely sick depending on the veterinarian diagnosis on infections, mechanical or internal dysfunctions; individuals were also classified according to three levels of body condition depending on the pectoral musculature mass: normal, low or high (extremely large muscle mass). In subsequent analyses, sickness was treated in two ways, one by separating mild and severe sickness and a second one by grouping all sick individuals into a single group. Gulls breed in springtime and differences may exist in the number of individuals of different sex or age visiting the dump throughout the year, so we defined two periods: one encompassing the breeding season, from March to July, and a second one including the non-breeding period, from August to February. The number of gulls captured by falcons and those killed by shooting were not equally distributed throughout the year, with comparatively more gulls shot during the breeding period. As a consequence, we randomly balanced the sample size between both periods to assure the comparability of the two data sets.
Individuals were shot at random (see above), with the same probability for all birds at the dump to be shot and shooting did not occur in a spatially segregated manner. Hence, comparison between the birds shot and those killed by falcons should reveal any predator preferences. We first used a MCA analysis including age, muscle condition, and sickness to visualize patterns of differentiation among gulls depending on the type of capture and the correlation among traits. We then used logistic regressions to test for differences between the two groups. We assessed the goodness-of-fit of the saturated model by estimating the overdispersion parameter (a value close to 1 indicated a good fit of the data). Model selection was based on Akaike's Information Criterion (AIC); the model with the lowest AIC value was considered as the best compromise between model deviance and model parameters [39]. We also calculated the AIC weight as a measure of relative plausibility of each model. To avoid model over-parameterization, only additive models were considered and individual age was treated as a continuous covariate. All statistical analyses were done using the software R (www.r-project.org) and SPSS (version 16.0).

Ethics Statement
Animals were killed as a part of a culling program that the Local Government began for conservation issues as well as public health concerns. All animal work has been conducted according to relevant national and international guidelines and permits were Table 4. Estimates of the factors affecting predation from the best ranked model, which included age as a continuous covariate, muscle condition as a factor with three levels (normal, low and high), and sickness, as a factor with two levels (healthy and sick).  provided by Conselleria de Medi Ambient (Govern Balear). Authors were not responsible nor executed the culling programme.