Skip to main content

Genetic and environmental variation in condition, cutaneous immunity, and haematocrit in house wrens



Life-history studies of wild bird populations often focus on the relationship between an individual’s condition and its capacity to mount an immune response, as measured by a commonly-employed assay of cutaneous immunity, the PHA skin test. In addition, haematocrit, the packed cell volume in relation to total blood volume, is often measured as an indicator of physiological performance. A multi-year study of a wild population of house wrens has recently revealed that those exhibiting the highest condition and strongest PHA responses as nestlings are most likely to be recruited to the breeding population and to breed through two years of age; in contrast, intermediate haematocrit values result in the highest recruitment to the population. Selection theory would predict, therefore, that most of the underlying genetic variation in these traits should be exhausted resulting in low heritability, although such traits may also exhibit low heritability because of increased residual variance. Here, we examine the genetic and environmental variation in condition, cutaneous immunity, and haematocrit using an animal model based on a pedigree of approximately 2,800 house wrens.


Environmental effects played a paramount role in shaping the expression of the fitness-related traits measured in this wild population, but two of them, condition and haematocrit, retained significant heritable variation. Condition was also positively correlated with both the PHA response and haematocrit, but in the absence of any significant genetic correlations, it appears that this covariance arises through parallel effects of the environment acting on this suite of traits.


The maintenance of genetic variation in different measures of condition appears to be a pervasive feature of wild bird populations, in contradiction of conventional selection theory. A major challenge in future studies will be to explain how such variation persists in the face of the directional selection acting on condition in house wrens and other species.


An animal’s condition is defined by the pool of resources available to allocate to the production and maintenance of traits that enhance fitness, and reflects an individual’s ability to acquire food, avoid predators, and resist disease [1]. Body condition is critically important to life history because a variety of traits affecting fitness are often condition dependent, such as those involved in mate attraction [1] or disease resistance [2]. A number of morphological and physiological indices for measuring condition have been devised, although the interpretation of these measures is often controversial and their validity frequently questioned [3]. In field studies of birds and mammals, an individual’s condition is often measured as body mass corrected for structural body size, which is often derived from the residuals of a regression of mass on body size [4]. Condition measured in this manner is assumed to reflect lipid reserves [5], a major form of energy storage in birds that is used to fuel the various processes that promote individual fitness (i.e., flight, migration, reproduction). Indeed, body-mass based measures appear to be tightly linked with fitness in a number of vertebrate species [6].

In studies of wild bird populations, the condition of both nestlings and breeding adults is often measured in the context of understanding life-history variation [7]. For example, all else being equal, individuals in better condition should be able to invest more in resource-limited traits associated with survival and/or reproduction. Such traits include an individual’s capacity to mount an immune response, which is an important component of maintenance [8]-[11]. This is particularly relevant to nestlings that, owing to their lack of mobility, are especially vulnerable to parasites and pathogens. A standard method employed for measuring cutaneous immune activity in wild birds is the phytohaemagglutinin (PHA) skin test, an assay that involves injecting a novel plant mitogen, PHA, into the wing prepatagium of the bird and measuring the ensuing swelling 24 h later. The swelling is taken as a measure of cutaneous immune responsiveness, which includes both innate and adaptive components of the immune system [12]-[15]. Although the conventional expectation has been that birds in good condition should be capable of mounting more robust immune responses [7], this expectation has been met in some studies [9],[16], but not in others [11],[17].

In addition to morphological measures of condition, haematocrit, the packed volume of red blood cells in relation to total blood volume, is often measured as a presumed physiological indicator of condition in field studies of wild birds [18],[19]. Because haematocrit is directly related to oxygen uptake, it is typically regarded as a measure of physiological performance [18]-[20], but its value as an indicator of condition or health state remains uncertain [11],[17],[21]. A growing body of evidence suggests, however, that physiological variation in haematocrit at different life-history stages is tied directly to such fundamental processes as reproduction, migration, and the acquisition of flight [22], and so there is reason to believe that haematocrit may be linked, at least under some circumstances, to fitness.

Over the past 10 years, we have routinely measured condition, PHA response, and haematocrit in a wild population of house wrens, Troglodytes aedon, in north-central Illinois, USA [23]-[26]. A multi-year study of a large subset of the nestlings measured over this time period revealed that those in the best condition and with the strongest PHA responses were more likely to be recruited to the breeding population and more likely to breed through two years of age [27]. The relationship between haematocrit and fitness was somewhat more complex, as intermediate haematocrit values resulted in the highest recruitment to our population, suggestive of stabilizing selection [27].

If, as our data on recruitment and subsequent breeding success would suggest, condition, immune response to PHA injection, and haematocrit in this population are subject to strong natural selection, we would anticipate that most of the underlying genetic variation in these traits would be exhausted [28]. However, a number of studies of wild bird populations have reported significant heritable variation in each of these traits [condition: [6],[29],[30]; PHA response: [10],[31],[32]; haematocrit: [33]]. This is, perhaps, most surprising in the case of condition, which is expected to reflect an individual’s current nutritional state and thus be determined primarily by environmental variation [6]. Here, we examine the genetic and environmental variation in condition, cutaneous immune responsiveness, and haematocrit over a three-year period (2004–2006) during which we used DNA profiling to determine the parentage of nestlings. Specifically, our objectives were to determine the degree of heritable variation in each of these traits and to measure the phenotypic and genetic correlations among these characters.


Study animals

House wrens are small (10–12 g), insectivorous songbirds that are sexually monomorphic in size and plumage. As obligate cavity-nesters, house wrens readily accept nestboxes within which to build their nests. Upon arrival at the study area in north-central Illinois after spring migration, females select a male that is defending a nest site and, after completing the nest, lay a clutch of 4–8 eggs. About half of the females that successfully rear their first brood produce a second [26],[34],[35]. Only female house wrens incubate the eggs and brood the nestlings, but both adults provision nestlings and fledglings. House wrens are usually socially monogamous, but social polygyny occurs and extra-pair fertilizations are common in the study population [23],[36]. Additional information concerning the breeding biology of house wrens is provided by Johnson [37].

We conducted this study during the 2004–2006 breeding seasons on the Mackinaw study site (40°40′N, 88°53′W) in a second-growth, deciduous forest bordering the Mackinaw River. Nestboxes (N =700) of uniform construction [38] were spaced 30 m apart along north–south transects separated by 60 m (Figure one in [39]). All nestboxes were mounted on 1.5-m metal poles that had been greased or under which 48.3-cm diameter aluminium baffles had been mounted to discourage nest predators. The nestboxes used in this study (N =301) were located in three semi-isolated neighbourhoods within the study area [23]. All activities complied with the Illinois State University Institutional Animal Care and Use Committee (Protocol No.17-2003) and United States Geological Survey banding permit 09211.

Field procedures

We visited all nestboxes every 1–3 days, noting their contents, the behaviours of any wrens present, and male identity as revealed by unique combinations of coloured, Darvic leg rings. We trapped and individually ringed any previously uncaught females when incubating or males when provisioning nestlings. To obtain DNA for paternity analyses, we collected a blood sample from adults upon capture and from all nestlings. On brood-day 11 or 12 (brood-day 0 is the day the first egg of a clutch hatches), we ringed nestlings with a numbered aluminium band, weighed them on an electronic balance (Acculab, Pocket Pro 250-B) to the nearest 0.1 g, and measured their right tarsus length with dial callipers to the nearest 0.1 mm. At the same time, we obtained a blood sample (≈75 μL) from the brachial vein in heparinized microcapillary tubes that were stored on ice in coolers in the field. Blood was taken to the laboratory later the same day to be centrifuged at 1,610 g for 60 s (Hematastat II, Separation Technologies) to separate cellular and plasma components. We measured haematocrit as the percentage of whole blood constituted by red blood cells, using the mean of three measurements. Red blood cells were stored in lysis buffer at 4°C and plasma at −20°C until further analysis.

We induced a cutaneous immune response in nestlings on brood-day 11 by injecting 50 μL of sterile phosphate buffered saline (PBS) containing phytohaemagglutinin (Sigma Aldrich, St. Louis, MO, USA) at a concentration of 5 mg mL−1 [40] into their prepatagium (wing-web). Phytohaemagglutinin (PHA), a plant-derived mitogen, induces a measurable tissue swelling as a result of responses from both the innate and adaptive axes of the immune system [12]; but see [15]. We used the change in wing-web thickness (mean of three successive readings before and 24 h after PHA injection) as a measure of cutaneous immune activity. Wing-web thickness was measured with a Mitutoyo thickness gauge (no. 547–500, Mitutoyo America Corp., Aurora, IL, USA). The repeatability of haematocrit and the PHA response was 0.99 and 0.95, respectively [41].

Determination of parentage

We isolated DNA from blood samples using a high-salt extraction protocol following Bruford et al. [42]. Polymerase chain reaction (PCR) amplifications were carried out in 15-μL volumes containing 200 μM dNTPs, 2.5 mM MgCl2, 1X PCR Buffer II (Applied Biosystems), and 0.133 μM forward and reverse primers. We used a thermal profile that followed the touchdown protocol described in Johnson et al. [43]. Forward primers were fluorescently labelled, and PCR products were analysed using a Beckman Coulter CEQ 8000 Genetic Analysis System. We typed all samples at three loci: TA-C3 (B)2 [44], Mcyμ4 [45], and LTMR6 [46]. When more resolution was needed, two additional loci, TA-A5-15 and TA-B4-2 [44], were used. We analysed allele data using Cervus 2.0. No locus deviated significantly from Hardy-Weinberg equilibrium, with the exception of TA-A5-15, whose null allele frequency was estimated at 0.094.

Attendant females matched nestlings at all loci with rare exceptions (i.e., six cases in which there was a mismatch at a single locus, attributable to mutation) and were assumed to be the genetic parent in all cases. At some nests, we were unable to obtain blood samples from attending adults, in which case, maternity of the attending female was assumed because intraspecific brood parasitism does not occur in our population [36]; paternity of nestlings in these cases was scored as unknown. Nestlings that matched attendant males at all loci were assigned within-pair paternity, and those that failed to match attendant males at two or more loci were designated as extra-pair. A few nestlings failed to match attendant males at one locus, so we re-typed at the anomalous locus to prevent typing error. If they still failed to match the attendant male, we typed them at additional loci to attempt to resolve the anomaly. In some cases, they failed to match the attendant male at one or more of these additional loci and were therefore designated as extra-pair nestlings. If the probability of false assignment [47] based on matching loci was lower than 0.005, then nestlings were designated as within-pair nestlings. In some cases, the paternity of nestlings could not be resolved as a result of loss of sample or when repeated attempts to genotype failed to produce reliable allele data for one or more loci because of the poor quality of the DNA sample. Exclusion probabilities for each locus were Mcyμ4, 0.811; LTMR6, 0.691; TA-C3(B)2, 0.841; TA-B4-2, 0.641; TA-A5-15, 0.464. For the three-locus set and five-locus set exclusion, probabilities were 0.991 and 0.998, respectively. Overall, the probability of false assignment for nestlings designated as within-pair was <0.008.

Based on these methods, we were able to construct a pedigree comprised of 2,768 house wrens, in which we were able to identify both parents of 1,661 individuals and only the mother of an additional 552 individuals. The pedigree spanned three generations, and included 229 paternal and 265 maternal identities. Not all phenotypic traits were measured on all individuals; sample sizes and mean values (±SE) for all traits are shown in Table 1. Unlike cross-fostering studies, our ability to separate nest effects (i.e., the effect of common environment) from genetic effects depends on the extent of extra-pair paternity and mate-switching across different reproductive episodes. As noted earlier, extra-pair paternity is common in our house-wren population (35-40% of broods; [23],[36]), and, on our study area, nearly 60% of pairs switch mates across broods produced during the same breeding season [34]. Birds returning to the study area in subsequent years (as was the case with many of the birds in our pedigree) invariably pair with a different mate from the one with which they have paired the previous year (pers. obs.). As a consequence of the collective force of these factors, a high proportion of females in our pedigree produced young with more than one male (39%), and an even higher proportion of males sired young with more than one female (49%) (Figure 1). Thus, our pedigree is well suited to separating common nest effects and genetic effects.

Table 1 Mean trait values (and associated sample sizes) for nestling house wrens included in the animal model analysis
Figure 1

Relative frequency of male house wrens males siring young with different numbers of females (solid bars) and the relative frequency of female house wrens producing young with different numbers of males (open bars).

Determination of sex

Sex was determined by amplifying sex-specific introns of the CHD-1 gene. Polymerase chain reactions (PCR) were carried out with sexing primers 1237L and 1272H [48] using a touchdown protocol as described in Johnson et al. [43], and the products electrophoresed through 2% agarose gels and stained with ethidium bromide. DNA isolated from adult house wrens of known sex was included in all sets of PCR runs as controls, and their PCR products were always electrophoresed with those of nestling samples. Because some nestlings did not survive to brood-day 11–12 and some DNA samples did not successfully amplify, not all nestlings in each nest were sexed.

Genetic analyses

Data were analysed using animal models fitted in ASReml (version 3). An animal model is a form of linear mixed-effect model in which an individual’s genetic merit is included as a random effect allowing estimation of the additive genetic (co)variance provided pedigree data are available [49],[50]. We first fitted univariate models to estimate the additive variance (VA) for each trait (tarsus, PHA, haematocrit, mass). Nest identity was included as an additional random effect to minimise the potential for upward bias in genetic variance from common environment effects. Nest effects are assumed to be normally distributed with a mean of zero and a variance to be estimated (Vnest). Similarly, we assume that residuals from the model are uncorrelated and normally distributed, with the residual variance (VR) arising from unmodelled environmental effects and/or measurement error. Fixed effects included year (2004, 2005, 2006), and sex (unknown, male, female) as three-level factors for all traits. In addition, we included a cubic regression of tarsus length as a fixed effect in the model of mass. By correcting for the average mass-tarsus length relationship, variance components (and derived parameters) from our model of mass can be viewed as pertaining to body condition (i.e. mass “corrected” for skeletal size). For each trait we determined heritability (h2) as the ratio of VA to VP (with VP estimated as the sum of the variance components; i.e., conditional on fixed effects) and the corresponding proportion of variance explained by nest effects (denoted n2). To test the genetic basis of hypothesised relationships among traits, we then fitted a multivariate model of all traits simultaneously with fixed effects as described above but no random effects. Note that the cubic function of tarsus length was included as a fixed effect on mass only, not the other response variables. This yielded an estimate of the phenotypic variance-covariance matrix (conditional on fixed effects) that we denote P. We then added individual genetic merit and nest effects back into the model to partition P into additive genetic (G), nest (N) and residual (R) matrices. Statistical inference was based on likelihood-ratio tests (LRT) for random effects and Wald F-statistics for fixed effects. For LRT of single random-effect terms in the univariate models, we assume the test statistic to be asymptotically distributed as a 50:50 mix of χ20 and χ21 [51]. For comparing multivariate models with differing random effects, we (conservatively) set the degrees of freedom to equal the number of additional (co)variance parameters estimated in the more complex model.


Fixed-effect estimates for nestling house wren traits are shown in Table 2 and their associated heritabilities are provided in Table 3. Nest effects accounted for a major portion of the variation in all traits, ranging from 21% of the variation in tarsus length to 54% of the variation in the PHA response (Table 3). Notwithstanding this large amount of environmental variation, both nestling haematocrit and condition (i.e., size-adjusted mass) were significantly heritable, with additive genetic variance accounting for about 15% of the phenotypic variation in each case (Table 3). Nestling condition differed across years, and male nestlings were of significantly higher body condition than female nestlings (Table 2). There were no significant year effects or sex differences in nestling haematocrit (Table 2).

Table 2 Univariate models showing fixed effect estimates for nestling house wren traits
Table 3 Estimated variance components of nestling house wren traits (conditioned on fixed effects) under the animal model, heritabilities (h 2 ), and nest effects (n 2 )

Neither tarsus length nor PHA response were significantly heritable (Table 3). There was a significant difference between male and female nestlings in tarsus length, with males exhibiting longer tarsi (Table 2). Tarsus length did not vary across the three years of the study. There was no sex difference in the PHA response, but PHA response varied across the three years of the study (Table 2).

Our multivariate modelling with no random effects included showed that P contained significant among-trait covariance (LRT comparison of model with full P matrix to one with diagonal elements only, i.e., all covariance terms set to zero; χ26 = 136, P < 0.001). This result was driven by significant positive phenotypic correlations between nestling body condition and each of the other three traits, tarsus length, PHA response, and haematocrit (Table 4). There was also a small, but significant correlation between tarsus length and PHA response. Partitioning the P matrix provided support for statistically significant covariance in N (full covariance matrix (full) versus variance only matrix (diagonal) LRT; χ26 = 22.4, P = 0.001) and R (full versus diagonal matrix LRT; χ26 = 38.5, P < 0.001). Among broods (i.e., nests), nestling body condition was positively associated with both tarsus length and haematocrit, but was not correlated with PHA response (Table 4). However, there was no evidence of significant covariance structure in G (full versus diagonal matrix LRT; χ26 = 3.24, P = 0.778), and pairwise genetic correlations were characterised by large standard errors in all cases (Table 4).

Table 4 Phenotypic covariance structure (P) and its additive genetic (G), nest (N) and residual (R) components


In keeping with a number of other studies of wild bird populations, the majority of variation in condition, cutaneous immune response, and haematocrit in our house wren population could be accounted for by nest effects, attesting to the paramount role of the early environment in influencing the expression of these traits. This is, perhaps, not surprising given that these traits often are subject to strong natural selection, which might be expected to exhaust any genetic variation underlying these traits [28],[52]-[54]. Nevertheless, two of the traits measured in this study, condition and haematocrit, were significantly heritable, while our measure of cutaneous immune responsiveness was not, in seeming contrast to some studies showing heritable variation in the PHA response [10],[31],[32], and other cross-fostering studies showing an absence of nest-of-origin effects on nestling haematocrit [17],[21]. Below, we explore the proximate underpinnings of the environmental and genetic factors influencing these important fitness-related traits, and the phenotypic and genetic relationships among them.

Nest effects accounted for the largest portion of phenotypic variation in both tarsus length and condition, but unlike tarsus length, condition was significantly heritable. Although we can only speculate as to the sources of environmental variation for these two traits, likely candidates include differential provisioning among broods [39],[55] or spatial variation in the availability of insect prey [16]. Parasites might also contribute to the observed nest effects on these traits [33], as nestlings in our house wren population are frequently infested with ectoparasitic mites [56]; however, in an earlier study of our population, nestling mass was not affected by ectoparasite numbers [56] and in a different house wren population, ectoparasites had no detectable effect on nestling tarsus and only a slight effect on mass [57]. There was also significant variation among years in nestling condition, which might also reflect temporal variation in prey availability arising from annual variation in rainfall and summer temperatures.

The absence of significant genetic variation in tarsus length was surprising, given that the heritability of tarsus length in birds typically hovers around 0.5 [58]. However, Christe et al., [33] also found that tarsus length was not significantly heritable in house martins, Delichon urbicum. Because environmental factors such as parasites and food availability can profoundly influence variation in tarsus length, heritability of tarsus may vary according to the environment in which it is measured. Moreover, heritability estimates based on the animal model are typically about 24% lower than those based on other methods such as parent-offspring resemblance [58], so our estimate of heritability of tarsus length in house wrens (0.11) is relatively conservative. Alternatively, low heritability of tarsus may be due to high residual variance rather than low additive genetic variation [59]. A long-term study of Savannah sparrows (Passerculus sandwichensis), for example, has shown that heritability can vary with the age at which traits are measured, and residual variation was generally much higher in traits expressed early in development (e.g., tarsus length in nestlings) [59].

Consistent with a number of other studies employing the animal model [6],[11],[29],[30], condition was significantly heritable, with additive genetic variation accounting for about 14% of the phenotypic variation in this trait. This result might be considered surprising, given that body condition often is subject to strong directional selection [2],[6],[60] and that nestling condition in our house wren population is positively correlated with recruitment to the breeding population and subsequent reproductive success [27]. However, assuming that condition is influenced by the expression of many loci, including those involved in the acquisition and efficient use of resources, this alone might sustain genetic variation in the face of strong directional selection [1]. In addition, selection may be acting primarily on the environmental component of variation in condition [6]. This occurs if the positive correlation found between condition and fitness reflects a shared dependence on environmental factors rather than a causal relationship [61].

PHA response was not significantly heritable, consistent with other studies employing the animal model to estimate additive genetic variance [11],[62], as well as those relying on parent-offspring resemblance [63],[64]. A prominent exception to this pattern is a study of common kestrels (Falco tinnunculus) employing the animal model that revealed a surprising high heritability of PHA response in fledglings (h 2 = 0.47-0.55) [65]; as with tarsus length discussed above, such differences could be due to age-related variation in the magnitude of residual variance. Although a number of cross-fostering studies also have reported significant heritable variation in the PHA response [10],[31],[32], these have relied on full-sib comparisons that likely produce inflated estimates because of dominance variance and maternal effects [62].

Nest effects accounted for over half the phenotypic variation in nestling PHA response, a result consistent with previous studies in this population showing that hatching date and time of injection can influence PHA response of nestling house wrens [24],[25],[41],[66]. Additional environmental factors might also include food availability, temperature, and parasite load [12], and we also cannot rule out the possibility of parental effects arising from differential provisioning of broods [39],[55] or transgenerational priming of immunity arising from the transfer of maternal antibodies in eggs [67]. Given the complex nature of the PHA response, which involves both innate and adaptive components of the immune system [12],[13], Martin et al. [12] have cautioned against interpreting larger swellings as indicative of stronger cell-mediated immunocompetence. In our house wren population, however, a multi-year study has revealed that nestlings with the strongest response to PHA injection have the highest likelihood of recruitment to the breeding population, and those with the strongest PHA responses are more likely to breed through two years of age [27]. Thus, the PHA response appears to be the focus of significant directional selection in our population and intimately associated with offspring fitness.

Nest effects had a major influence on nestling haematocrit, consistent with studies of other wild bird populations [17],[20],[21]. A number of environmental factors can influence nestling haematocrit, among them parasite load [19],[20], food availability [16],[68], and temperature [69]. Parasite load appears to be an especially important factor in other species, but ectoparasite loads in another population of house wrens did not affect nestling haematocrit [57]. In our house wren population, nestlings are frequently infested with two species of ectoparasitic mites [56] and their density varies greatly among nests [56]. Although a relationship between the density of mites and nestling haematocrit has not been established in our population, variability in the level of infestation could account for at least a portion of the nest effects detected in the present study, as it has been in other species [70]. In addition to variation in ectoparasite load, parental effects mediated through differences among broods in levels of parental provisioning of nestlings could further contribute to nest effects on nestling haematocrit [71]. Although we did not measure provisioning of nestlings in the present study, provisioning rates within and among pairs can vary widely [39],[55], and thus differences among broods in the amount of food provided to nestlings could explain some of the nest effects on nestling haematocrit.

In addition to the strong environmental effect on nestling haematocrit, haematocrit was significantly heritable, in contrast to cross-fostering studies of other species [17],[21] that have shown significant effects of nest-of-rearing on nestling haematocrit, but no effect of nest-of-origin. Additive genetic variation accounted for approximately 15% of the phenotypic variation in this trait. Studies that have employed cross-fostering [17],[21] or parent-offspring resemblance [20] to detect genetic effects on haematocrit have relied on much smaller sample sizes than the one on which our pedigree was based, and, thus, the absence of a genetic effect in these studies may have been due to a lack of statistical power (but see [33] for evidence of a significant nest-of-origin effect in house martins). In support of this possibility, a study of a captive population of zebra finches (Taeniopygia guttata) based on a pedigree of similar size as the one considered here [11], reported a heritability of 0.38 ± 0.077 (P <0.001) as determined by the animal model.

There was a significant positive phenotypic correlation, but no genetic correlation, between nestling haematocrit and condition. This suggests that the phenotypic correlation arises because of parallel effects of the common nest environment (including parental effects) on the two traits, and not because of any underlying pleiotropic effect. The absence of genetic correlations in this study should be regarded with some caution, however, as our power to detect such correlations given our sample sizes is probably limited. This is reflected in the large standard errors associated with our rG estimates. Increased precision could be obtained through larger sample sizes and/or the use of cross-fostering designs that better facilitate statistical separation of additive genetic from nest effects [72]. Both nestling condition and haematocrit appear to be closely related to fitness in our population, predictive of both recruitment to the breeding population and subsequent reproductive success [27], and would thus seem to represent different axes of condition, one reflecting fat stores and the other oxygen-uptake capacity. The relationship between haematocrit and fitness is more complex than the one between condition and fitness, however, as intermediate haematocrit values result in the highest recruitment to our population, suggestive of stabilizing selection [27]. Although oxygen transport generally increases with increasing haematocrit, increases in blood viscosity beyond a certain threshold may actually hinder oxygen transport [73]. Indeed, experimental studies in mice have shown that oxygen uptake and physical endurance are maximized at intermediate values of haematocrit [73].


Environmental effects played a paramount role in shaping the expression of the fitness-related traits measured in this wild population of house wrens, but two of them, condition and haematocrit, retained significant heritable variation. It is becoming increasingly evident that the maintenance of genetic variation in condition measures, in particular, appears to be a pervasive feature of wild bird populations, in apparent contradiction of conventional selection theory [6],[30],[74]. A major challenge in future studies will be to explain how such variation persists in the face of the directional selection acting on condition in house wrens and other species. Condition was also positively correlated with both the PHA response and haematocrit, but in the absence of any significant genetic correlations, it appears that this covariance arises through parallel effects of the environment on this suite of traits (but see [9]). This would seem to bolster the case for the utility of metrics designed to capture the multi-dimensionality of condition based on multiple physiological and morphological measures including mass, haematocrit, and immune responsiveness [3].

Availability of supporting data

The data supporting the results of this article are available in the Dryad digital repository doi:10.5061/dryad.jk2m0 [75] (


  1. 1.

    Rowe L, Houle D: The lek paradox and the capture of genetic variance by condition dependent traits. Proc R Soc Lond B. 1996, 263: 1415-1421. 10.1098/rspb.1996.0207.

    Article  Google Scholar 

  2. 2.

    Wilcoxen TE, Boughton RK, Schoech SJ: Selection on innate immunity and body condition in Florida scrub-jays throughout an epidemic. Biol Lett. 2010, 6: 552-554. 10.1098/rsbl.2009.1078.

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Milot E, Cohen AA, Vézina F, Buehler DM, Matson KD, Piersma T: A novel integrative method for measuring body condition in ecological studies based on physiological dysregulation. Methods Ecol Evol. 2014, 5: 146-155. 10.1111/2041-210X.12145.

    Article  Google Scholar 

  4. 4.

    Schulte-Hostedde AI, Zinner B, Millar JS, Hickling GJ: Restitution of mass-size residuals: validating body condition indices. Ecology. 2005, 86: 155-163. 10.1890/04-0232.

    Article  Google Scholar 

  5. 5.

    Ardia DR: Super size me: an experimental test of the factors affecting lipid content and the ability of residual body mass to predict lipid stores in nestling European starlings. Funct Ecol. 2005, 19: 414-420. 10.1111/j.1365-2435.2005.00997.x.

    Article  Google Scholar 

  6. 6.

    Merilä J, Kruuk LEB, Sheldon BC: Natural selection on the genetical component of variance in body condition in a wild bird population. J Evol Biol. 2001, 14: 918-929. 10.1046/j.1420-9101.2001.00353.x.

    Article  Google Scholar 

  7. 7.

    Norris K, Evans MR: Ecological immunology: life history trade-offs and immune defense in birds. Behav Ecol. 2000, 11: 19-26. 10.1093/beheco/11.1.19.

    Article  Google Scholar 

  8. 8.

    Tella JL, Scheuerlein A, Ricklefs RE: Is cell-mediated immunity related to the evolution of life-history strategies in birds?. Proc R Soc Lond B. 2002, 269: 1059-1066. 10.1098/rspb.2001.1951.

    Article  Google Scholar 

  9. 9.

    Gleeson DJ, Blows MW, Owens IPF: Genetic covariance between indices of body condition and immunocompetence in a passerine bird. BMC Evol Biol. 2005, 5: 61-10.1186/1471-2148-5-61.

    PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Cichoń M, Sendecka J, Gustafsson L: Genetic and environmental variation in immune response of collared flycatcher nestlings. J Evol Biol. 2006, 19: 1701-1706. 10.1111/j.1420-9101.2006.01110.x.

    PubMed  Article  Google Scholar 

  11. 11.

    Birkhead TR, Pellatt EJ, Matthews IM, Roddis NJ, Hunter FM, McPhie F, Castillo-Juarez H: Genic capture and the genetic basis of sexually selected traits in the zebra finch. Evolution. 2006, 60: 2389-2398. 10.1111/j.0014-3820.2006.tb01873.x.

    PubMed  Article  Google Scholar 

  12. 12.

    Martin LB, Han P, Lewittes J, Kuhlman JR, Lasing KC, Wikelski M: Phytohemagglutinin-induced skin swelling in birds: histological support for a classic immunoecological technique. Funct Ecol. 2006, 20: 290-299. 10.1111/j.1365-2435.2006.01094.x.

    Article  Google Scholar 

  13. 13.

    Tella JL, Lemus JA, Carrete M, Blanco G: The PHA test reflects acquired T-cell mediated immunocompetence in birds. PLoS One. 2008, 3: e3295-10.1371/journal.pone.0003295.

    PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Vinkler M, Bainová H, Albrecht T: Functional analysis of the skin-swelling response to phytohaemagglutinin. Funct Ecol. 2010, 24: 1081-1086. 10.1111/j.1365-2435.2010.01711.x.

    Article  Google Scholar 

  15. 15.

    Vinkler M, Svobodová J, Gabrielová B, Bainová H, Bryjová A: Cytokine expression in phytohaemagglutinin-induced skin inflammation in a galliform bird. J Avian Biol. 2014, 45: 43-50. 10.1111/j.1600-048X.2011.05860.x.

    Article  Google Scholar 

  16. 16.

    Hoi-Leitner M, Romero-Pujante M, Hoi H, Pavlova A: Food availability and immune capacity in serin (Serinus serinus) nestlings. Behav Ecol Sociobiol. 2001, 49: 333-339. 10.1007/s002650000310.

    Article  Google Scholar 

  17. 17.

    Morrison ES, Ardia DR, Clotfelter ED: Cross-fostering reveals sources of variation in innate immunity and hematocrit in nestling tree swallows Tachycineta bicolor . J Avian Biol. 2009, 40: 573-578. 10.1111/j.1600-048X.2009.04910.x.

    Article  Google Scholar 

  18. 18.

    Saino N, Cuervo JJ, Krivacek M, de Lope F, Møller AP: Experimental manipulation of tail ornament size affects the hematocrit of male barn swallows (Hirundo rustica). Oecologia. 1997, 110: 186-190. 10.1007/s004420050148.

    Article  Google Scholar 

  19. 19.

    Potti J, Moreno J, Merino S, Frías O, Rodríguez R: Environmental and genetic variation in the haematocrit of fledgling pied flycatchers Ficedula hypoleuca . Oecologia. 1999, 120: 1-8. 10.1007/s004420050826.

    Article  Google Scholar 

  20. 20.

    Simon A, Thomas DW, Bourgault P, Blondel J, Perret P, Lambrechts MM: Between-population differences in nestling size and hematocrit level in blue tits (Parus caeruleus): a cross-fostering test for genetic and environmental effects. Can J Zool. 2005, 83: 694-701. 10.1139/z05-059.

    Article  Google Scholar 

  21. 21.

    Potti J: Variation in the hematocrit of a passerine bird across life stages is mainly of environmental origin. J Avian Biol. 2007, 38: 726-730. 10.1111/j.2007.0908-8857.04073.x.

    Article  Google Scholar 

  22. 22.

    Williams TD: Physiological Adaptations for Breeding in Birds. 2012, Princeton University Press, Princeton

    Google Scholar 

  23. 23.

    Forsman AM, Vogel LA, Sakaluk SK, Johnson BG, Masters BS, Johnson LS, Thompson CF: Female house wrens (Troglodytes aedon) increase the size, but not immunocompetence, of their offspring through extra-pair mating. Mol Ecol. 2008, 17: 3697-3706. 10.1111/j.1365-294X.2008.03860.x.

    PubMed  Article  Google Scholar 

  24. 24.

    Grana SC, Sakaluk SK, Bowden RM, Doellman MA, Vogel LA, Thompson CF: Reproductive allocation in female house wrens is not influenced by experimentally altered male attractiveness. Behav Ecol Sociobiol. 2012, 66: 1247-1258. 10.1007/s00265-012-1378-4.

    Article  Google Scholar 

  25. 25.

    Sutherland JL, Thompson CF, Sakaluk SK: No effect of carotenoid supplementation on PHA response or body condition of nestling house wrens. Physiol Biochem Zool. 2012, 85: 21-28. 10.1086/663353.

    PubMed  CAS  Article  Google Scholar 

  26. 26.

    Bowers EK, Smith RA, Hodges CJ, Zimmerman LM, Thompson CF, Sakaluk SK: Sex-biased terminal investment in offspring induced by maternal immune challenge in the house wren (Troglodytes aedon). Proc R Soc B. 2012, 279: 2891-2898. 10.1098/rspb.2012.0443.

    PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Bowers EK, Hodges CJ, Forsman AM, Vogel LA, Masters BS, Johnson BGP, Johnson LS, Thompson CF, Sakaluk SK: Neonatal body condition, immune responsiveness, and hematocrit predict longevity in a wild bird population. Ecology. 2014, 95: 3027-3034. 10.1890/14-0418.1.

    PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Falconer DS: Introduction to Quantitative Genetics. 1981, Longman, New York

    Google Scholar 

  29. 29.

    Jensen H, Sæther B-E, Ringby TH, Tufto J, Griffith SC, Ellegren H: Sexual variation in heritability and genetic correlations of morphological traits in house sparrow (Passer domesticus). J Evol Biol. 2003, 16: 1296-1307. 10.1046/j.1420-9101.2003.00614.x.

    PubMed  CAS  Article  Google Scholar 

  30. 30.

    Gienapp P, Merilä J: Genetic and environmental effects on a condition-dependent trait: feather growth in Siberian jays. J Evol Biol. 2010, 23: 715-723. 10.1111/j.1420-9101.2010.01949.x.

    PubMed  CAS  Article  Google Scholar 

  31. 31.

    Ardia DR: Cross-fostering reveals an effect of spleen size and nest temperatures on immune responses in nestling European starlings. Oecologia. 2005, 145: 327-334. 10.1007/s00442-005-0120-6.

    PubMed  Article  Google Scholar 

  32. 32.

    Drobniak SM, Wiejaczka D, Arct A, Dubiec A, Gustafsson L, Cichoń M: Sex-specific heritability of cell-mediated immune response in the blue tit nestlings (Cyanistes caeruleus). J Evol Biol. 2010, 23: 1286-1292. 10.1111/j.1420-9101.2010.01993.x.

    PubMed  Article  Google Scholar 

  33. 33.

    Christe P, Møller AP, Saino N, de Lope F: Genetic and environmental components of phenotypic variation in immune response and body size of a colonial bird, Delichon urbica (the house martin). Heredity. 2000, 85: 75-83. 10.1046/j.1365-2540.2000.00732.x.

    PubMed  Article  Google Scholar 

  34. 34.

    Drilling NE, Thompson CF: Mate switching in multibrooded house wrens. Auk. 1991, 108: 60-70.

    Google Scholar 

  35. 35.

    Bowers EK, Sakaluk SK, Thompson CF: Experimentally increased egg production constrains future reproduction of female house wrens. Anim Behav. 2012, 83: 495-500. 10.1016/j.anbehav.2011.11.026.

    Article  Google Scholar 

  36. 36.

    Soukup SS, Thompson CF: Social mating system affects the frequency of extra-pair paternity in house wrens. Anim Behav. 1997, 54: 1089-1105. 10.1006/anbe.1997.0556.

    PubMed  Article  Google Scholar 

  37. 37.

    Johnson LS: House Wren (Troglodytes aedon). The Birds of North America No. 380. Edited by: Poole A, Gill F. 1998, The Birds of North America, Inc, Philadelphia

    Google Scholar 

  38. 38.

    Lambrechts MM, Adriaensen F, Ardia DR, Artemyev AV, Atiénzar F, Bańbura J, Barba E, Bouvier J–C, Camprodon J, Cooper CB, Dawson RD, Eens M, Eeva T, Faivre B, Garamszegi LZ, Goodenough AE, Gosler AG, Grégoire A, Griffith SC, Gustafsson L, Johnson LS, Kania W, Keišs O, Llambías PE, Mainwaring MC, Mänd R, Massa B, Mazgajski TD, Møller AP, Moreno J: The design of artificial nestboxes for the study of secondary hole-nesting birds: a review of methodological inconsistencies and potential biases. Acta Ornithol. 2010, 45: 1-26. 10.3161/000164510X516047.

    Article  Google Scholar 

  39. 39.

    DeMory ML, Thompson CF, Sakaluk SK: Male quality influences male provisioning in house wrens independent of attractiveness. Behav Ecol. 2010, 21: 1156-1164. 10.1093/beheco/arq123.

    Article  Google Scholar 

  40. 40.

    Smits JE, Bortolotti GR, Tella JL: Simplifying the phytohaemagglutinin skin-testing technique in studies of avian immunocompetence. Funct Ecol. 1999, 13: 567-572. 10.1046/j.1365-2435.1999.00338.x.

    Article  Google Scholar 

  41. 41.

    Forsman AM, Sakaluk SK, Thompson CF, Vogel LA: Cutaneous immune activity, but not innate immune responsiveness, covaries with mass and environment in nestling house wrens (Troglodytes aedon). Physiol Biochem Zool. 2010, 83: 512-518. 10.1086/649894.

    PubMed  Article  Google Scholar 

  42. 42.

    Bruford MW, Hanotte O, Brookfield JFY, Burke T: Single-Locus and Multilocus DNA Fingerprinting. Molecular Genetic Analysis of Populations. Edited by: Hoelzel AR. 1992, Oxford University Press, Oxford, 225-269.

    Google Scholar 

  43. 43.

    Johnson LS, Hicks BG, Masters BS: Increased cuckoldry as a cost of breeding late for male house wrens (Troglodytes aedon). Behav Ecol. 2002, 13: 670-675. 10.1093/beheco/13.5.670.

    Article  Google Scholar 

  44. 44.

    Cabe PR, Marshall KE: Microsatellite loci from the house wren (Troglodytes aedon). Mol Ecol Notes. 2001, 1: 155-156. 10.1046/j.1471-8278.2001.00057.x.

    CAS  Article  Google Scholar 

  45. 45.

    Double MC, Dawson D, Burke T, Cockburn A: Finding the fathers in the least faithful bird: a microsatellite-based genotyping system for the superb fairy-wren Malurus cyaneus . Mol Ecol. 1997, 6: 691-693. 10.1046/j.1365-294X.1997.00228.x.

    CAS  Article  Google Scholar 

  46. 46.

    McDonald DB, Potts WK: Cooperative display and relatedness among males in a lek-mating bird. Science. 1994, 266: 1030-1032. 10.1126/science.7973654.

    PubMed  CAS  Article  Google Scholar 

  47. 47.

    Masters BS, Hicks BG, Johnson LS, Erb LA: Genotype and extra-pair paternity in the house wren: a rare-male effect?. Proc R Soc Lond B. 2003, 270: 1393-1397. 10.1098/rspb.2003.2380.

    Article  Google Scholar 

  48. 48.

    Kahn NW, John J, St. Quinn TW: Chromosome-specific intron size differences in the avian CHD gene provides an efficient method for sex identification in birds. Auk. 1998, 115: 1074-1078. 10.2307/4089527.

    Article  Google Scholar 

  49. 49.

    Kruuk LEB: Estimating genetic parameters in natural populations using the ‘animal model’. Phil Trans R Soc Lond B. 2004, 359: 873-890. 10.1098/rstb.2003.1437.

    Article  Google Scholar 

  50. 50.

    Wilson AJ, Réale D, Clements MN, Morrissey MM, Postma E, Walling CA, Kruuk LEB, Nussey DH: An ecologist’s guide to the animal model. J Anim Ecol. 2010, 79: 13-26. 10.1111/j.1365-2656.2009.01639.x.

    PubMed  Article  Google Scholar 

  51. 51.

    Visscher PM: A note on the asymptotic distribution of likelihood ratio tests to test variance components. Twin Res Hum Genet. 2006, 9: 490-495. 10.1375/twin.9.4.490.

    PubMed  Article  Google Scholar 

  52. 52.

    Price T, Schluter D: On the low heritability of life-history traits. Evolution. 1991, 45: 853-861. 10.2307/2409693.

    Article  Google Scholar 

  53. 53.

    Houle D: Comparing evolvability and variability of quantitative traits. Genetics. 1992, 130: 195-204.

    PubMed  CAS  PubMed Central  Google Scholar 

  54. 54.

    McCleery RH, Pettifor RA, Armbruster P, Meyer K, Sheldon BC, Perrins CM: Components of variance underlying fitness in a natural population of great tits Parus major . Am Nat. 2004, 164: E62-E72. 10.1086/422660.

    PubMed  CAS  Article  Google Scholar 

  55. 55.

    Barnett CA, Thompson CF, Sakaluk SK: Aggressiveness, boldness and parental food provisioning in male house wrens (Troglodytes aedon). Ethology. 2012, 118: 984-993. 10.1111/j.1439-0310.2012.02092.x.

    Article  Google Scholar 

  56. 56.

    Pacejka AJ, Gratton CM, Thompson CF: Do potentially virulent mites affect house wren (Troglodytes aedon) reproductive success?. Ecology. 1998, 79: 1797-1806. 10.1890/0012-9658(1998)079[1797:DPVMAH]2.0.CO;2.

    Article  Google Scholar 

  57. 57.

    Johnson LS, Albrecht DJ: Effects of haematophagous ectoparasites on nestling house wrens, Troglodytes aedon: who pays the cost of parasitism?. Oikos. 1993, 66: 255-262. 10.2307/3544812.

    Article  Google Scholar 

  58. 58.

    Postma E: Four Decades of Estimating Heritabilities in Wild Vertebrate Populations: Improved Methods, More Data, Better Estimates?. Quantitative Genetics in the Wild. Edited by: Charmantier A, Garant D, Kruuk LEB. 2014, Oxford University Press, Oxford, 16-33. 10.1093/acprof:oso/9780199674237.003.0002.

    Google Scholar 

  59. 59.

    Wheelwright NT, Keller LF, Postma E: The effect of trait type and strength of selection on heritability and evolvability in an island bird population. Evolution. 2014, 68: 3325-3336. 10.1111/evo.12499.

    PubMed  Article  Google Scholar 

  60. 60.

    Merilä J, Kruuk LEB, Sheldon BC: Cryptic evolution in a wild bird population. Nature. 2001, 412: 76-79. 10.1038/35083580.

    PubMed  Article  Google Scholar 

  61. 61.

    Morrissey MB, Kruuk LEB, Wilson AJ: The danger of applying the breeder’s equation outside of the context of artificial selection. J Evol Biol. 2010, 23: 2277-2288. 10.1111/j.1420-9101.2010.02084.x.

    PubMed  CAS  Article  Google Scholar 

  62. 62.

    Pitala N, Gustafsson L, Sendecka J, Brommer JE: Nestling immune response to phytohaemagglutinin is not heritable in collared flycatchers. Biol Lett. 2007, 3: 418-421. 10.1098/rsbl.2007.0135.

    PubMed  PubMed Central  Article  Google Scholar 

  63. 63.

    Cucco M, Malacarne G, Ottonelli R, Patrone M: Repeatability of cell-mediated and innate immunity, and other fitness-related traits, in the grey partridge. Can J Zool. 2006, 84: 72-79. 10.1139/z05-179.

    Article  Google Scholar 

  64. 64.

    Bonato M, Evans MR, Hasselquist D, Sherley RB, Cloete SWP, Cherry MI: Ostrich chick humoral immune responses and growth rate are predicted by parental immune responses and paternal coloration. Behav Ecol Sociobiol. 2013, 67: 1891-1901. 10.1007/s00265-013-1597-3.

    Article  Google Scholar 

  65. 65.

    Kim S-Y, Fargallo JA, Vergara O, Martínez-Padilla J: Multivariate heredity of melanin-based coloration, body mass and immunity. Heredity. 2013, 111: 139-146. 10.1038/hdy.2013.29.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  66. 66.

    Thompson CF, Sakaluk SK, Masters BS, Johnson BPG, Vogel LA, Forsman AM, Johnson LS: Condition-dependent sex difference in nestling house wren (Troglodytes aedon) response to phytohaemagglutinin injection. Can J Zool. 2014, 92: 1-7. 10.1139/cjz-2013-0140.

    CAS  Article  Google Scholar 

  67. 67.

    Grindstaff JL, Hasselquist D, Nilsson J-Å, Sandell M, Smith HG, Stjerman M: Transgenerational priming of immunity: maternal exposure to a bacterial antigen enhances offspring humoral immunity. Proc R Soc B. 2006, 273: 2551-2557. 10.1098/rspb.2006.3608.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  68. 68.

    Merino S, Potti J: Growth, nutrition, and blow fly parasitism in nestling pied flycatchers. Can J Zool. 1998, 76: 936-941. 10.1139/z98-013.

    Article  Google Scholar 

  69. 69.

    Rehder NB, Bird DM: Annual profiles of blood packed cell volumes of captive American kestrels. Can J Zool. 1983, 61: 2550-2555. 10.1139/z83-337.

    Article  Google Scholar 

  70. 70.

    Richner H, Oppliger A, Christe P: Effect of an ectoparasite on reproduction in great tits. J Anim Ecol. 1993, 62: 703-710. 10.2307/5390.

    Article  Google Scholar 

  71. 71.

    Santangeli A, Hakkarainen H, Laaksonen T, Korpimäki E: Home range size is determined by habitat composition but feeding rate by food availability in male Tengmalm’s owls. Anim Behav. 2012, 83: 1115-1123. 10.1016/j.anbehav.2012.02.002.

    Article  Google Scholar 

  72. 72.

    Kruuk LEB, Hadfield JD: How to separate genetic and environmental causes of similarity between relatives. J Evol Biol. 2007, 20: 1890-1903. 10.1111/j.1420-9101.2007.01377.x.

    PubMed  CAS  Article  Google Scholar 

  73. 73.

    Schuler B, Arras M, Keller S, Rettich A, Lundby C, Vogel J, Gassmann M: Optimal hematocrit for maximal exercise performance in acute and chronic erythropoietin-treated mice. Proc Natl Acad Sci U S A. 2010, 107: 419-423. 10.1073/pnas.0912924107.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  74. 74.

    Kruuk LEB, Merilä J, Sheldon BC: Phenotypic selection on a heritable size trait revisited. Am Nat. 2001, 158: 557-581. 10.1086/323585.

    PubMed  CAS  Article  Google Scholar 

  75. 75.

    Sakaluk SK, Wilson AJ, Bowers EK, Johnson LS, Masters BS, Johnson BG, Vogel LA, Forsman AM, Thompson CF: Data from: Genetic and environmental variation in condition, cutaneous immunity, and haematocrit in house wrens. Dryad. doi:10.5061/dryad.jk2m0.

Download references


We thank the 2004–2006 Wren Crews for field assistance and the ParkLands Foundation (Merwin Preserve) and the Sears and Butler families for the use of their properties. Financial support was provided by NSF grants GK12-0086354, IBN-0316580, IOS-0718140 and IOS-1118160; NIH grant R15HD076308-01; a visiting professorship from the Leverhulme Trust (SKS); the School of Biological Sciences, Illinois State University; a BBSRC David Phillips Research Fellowship (AJW); and student-research grants from the Beta Lambda Chapter of the Phi Sigma Biological Sciences Honor Society (AMF).

Author information



Corresponding author

Correspondence to Scott K Sakaluk.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SKS, CFT, EKB, and AJW conceived the research theme, SKS, CFT, and AMF conducted the field work and collected the blood samples, LAV supervised the immunological work, LSJ, BSM, and BGPJ performed the molecular analyses of sex and parentage of nestlings, AJW and SKS performed the statistical analyses, and SKS wrote the original draft of the manuscript. All authors contributed to the manuscript writing and approved its final version.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

This article is published under license to BioMed Central Ltd. 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 work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sakaluk, S.K., Wilson, A.J., Bowers, E.K. et al. Genetic and environmental variation in condition, cutaneous immunity, and haematocrit in house wrens. BMC Evol Biol 14, 242 (2014).

Download citation


  • Animal model
  • Condition
  • Haematocrit
  • Heritability
  • Genetic variation
  • Immunity
  • Life-history theory
  • PHA
  • Troglodytes aedon