BACKGROUND
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  <cit> .
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  <cit>  or disease resistance  <cit> .
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  <cit> .
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  <cit> .
condition measured in this manner is assumed to reflect lipid reserves  <cit> , a major form of energy storage in birds that is used to fuel the various processes that promote individual fitness .
indeed, body-mass based measures appear to be tightly linked with fitness in a number of vertebrate species  <cit> .
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  <cit> .
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  <cit> .
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  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  <cit> .
although the conventional expectation has been that birds in good condition should be capable of mounting more robust immune responses  <cit> , this expectation has been met in some studies  <cit> , but not in others  <cit> .
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  <cit> .
because haematocrit is directly related to oxygen uptake, it is typically regarded as a measure of physiological performance  <cit> , but its value as an indicator of condition or health state remains uncertain  <cit> .
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  <cit> , 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  <cit> .
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  <cit> .
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  <cit> .
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  <cit> .
however, a number of studies of wild bird populations have reported significant heritable variation in each of these traits .
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  <cit> .
here, we examine the genetic and environmental variation in condition, cutaneous immune responsiveness, and haematocrit over a three-year period  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.
methods
study animals
house wrens are small , 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– <dig> eggs.
about half of the females that successfully rear their first brood produce a second  <cit> .
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  <cit> .
additional information concerning the breeding biology of house wrens is provided by johnson  <cit> .
we conducted this study during the 2004– <dig> breeding seasons on the mackinaw study site  in a second-growth, deciduous forest bordering the mackinaw river.
nestboxes  of uniform construction  <cit>  were spaced 30 m apart along north–south transects separated by 60 m .
all nestboxes were mounted on  <dig> -m metal poles that had been greased or under which  <dig> -cm diameter aluminium baffles had been mounted to discourage nest predators.
the nestboxes used in this study  were located in three semi-isolated neighbourhoods within the study area  <cit> .
all activities complied with the illinois state university institutional animal care and use committee  and united states geological survey banding permit  <dig>
field procedures
we visited all nestboxes every 1– <dig> 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  <dig> or  <dig> , we ringed nestlings with a numbered aluminium band, weighed them on an electronic balance  to the nearest  <dig>  g, and measured their right tarsus length with dial callipers to the nearest  <dig>  mm.
at the same time, we obtained a blood sample  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  <dig>  g for 60 s  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  <dig> by injecting 50 μl of sterile phosphate buffered saline  containing phytohaemagglutinin  at a concentration of 5 mg ml− <dig>  <cit>  into their prepatagium .
phytohaemagglutinin , 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  <cit> ; but see  <cit> .
we used the change in wing-web thickness  as a measure of cutaneous immune activity.
wing-web thickness was measured with a mitutoyo thickness gauge .
the repeatability of haematocrit and the pha response was  <dig>  and  <dig> , respectively  <cit> .
determination of parentage
we isolated dna from blood samples using a high-salt extraction protocol following bruford et al.
<cit> .
polymerase chain reaction  amplifications were carried out in 15-μl volumes containing 200 μm dntps,  <dig>  mm mgcl <dig>  1x pcr buffer ii , and  <dig>  μm forward and reverse primers.
we used a thermal profile that followed the touchdown protocol described in johnson et al.
<cit> .
forward primers were fluorescently labelled, and pcr products were analysed using a beckman coulter ceq  <dig> genetic analysis system.
we typed all samples at three loci: ta-c <dig>  <dig>  <cit> , mcyμ <dig>  <cit> , and ltmr <dig>  <cit> .
when more resolution was needed, two additional loci, ta-a5- <dig> and ta-b4- <dig>  <cit> , were used.
we analysed allele data using cervus  <dig> .
no locus deviated significantly from hardy-weinberg equilibrium, with the exception of ta-a5- <dig>  whose null allele frequency was estimated at  <dig> .
attendant females matched nestlings at all loci with rare exceptions  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  <cit> ; 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  <cit>  based on matching loci was lower than  <dig> , 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μ <dig>   <dig> ; ltmr <dig>   <dig> ; ta-c <dig>   <dig> ; ta-b4- <dig>   <dig> ; ta-a5- <dig>   <dig> .
for the three-locus set and five-locus set exclusion, probabilities were  <dig>  and  <dig> , respectively.
overall, the probability of false assignment for nestlings designated as within-pair was < <dig> .
based on these methods, we were able to construct a pedigree comprised of  <dig>  house wrens, in which we were able to identify both parents of  <dig>  individuals and only the mother of an additional  <dig> individuals.
the pedigree spanned three generations, and included  <dig> paternal and  <dig> maternal identities.
not all phenotypic traits were measured on all individuals; sample sizes and mean values  for all traits are shown in table  <dig>  unlike cross-fostering studies, our ability to separate nest effects  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 , and, on our study area, nearly 60% of pairs switch mates across broods produced during the same breeding season  <cit> .
birds returning to the study area in subsequent years  invariably pair with a different mate from the one with which they have paired the previous year .
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 , and an even higher proportion of males sired young with more than one female  .
thus, our pedigree is well suited to separating common nest effects and genetic effects.table  <dig>
mean trait values  for nestling house wrens included in the animal model analysis
trait
n
mean
sdrelative frequency of male house wrens males siring young with different numbers of females  and the relative frequency of female house wrens producing young with different numbers of males .
determination of sex
sex was determined by amplifying sex-specific introns of the chd- <dig> gene.
polymerase chain reactions  were carried out with sexing primers 1237l and 1272h  <cit>  using a touchdown protocol as described in johnson et al.
<cit> , 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– <dig> 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 .
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 variance provided pedigree data are available  <cit> .
we first fitted univariate models to estimate the additive variance  for each trait .
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 .
similarly, we assume that residuals from the model are uncorrelated and normally distributed, with the residual variance  arising from unmodelled environmental effects and/or measurement error.
fixed effects included year , and sex  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  from our model of mass can be viewed as pertaining to body condition .
for each trait we determined heritability  as the ratio of va to vp  and the corresponding proportion of variance explained by nest effects .
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  that we denote p. we then added individual genetic merit and nest effects back into the model to partition p into additive genetic , nest  and residual  matrices.
statistical inference was based on likelihood-ratio tests  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: <dig> mix of χ <dig> and χ <dig>  <cit> .
for comparing multivariate models with differing random effects, we  set the degrees of freedom to equal the number of additional variance parameters estimated in the more complex model.
RESULTS
fixed-effect estimates for nestling house wren traits are shown in table  <dig> and their associated heritabilities are provided in table  <dig>  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 .
notwithstanding this large amount of environmental variation, both nestling haematocrit and condition  were significantly heritable, with additive genetic variance accounting for about 15% of the phenotypic variation in each case .
nestling condition differed across years, and male nestlings were of significantly higher body condition than female nestlings .
there were no significant year effects or sex differences in nestling haematocrit .table  <dig>
univariate models showing fixed effect estimates for nestling house wren traits
trait/effect
level
coefficient ± se
df
num
df
den
f
p
year1
year1
year1
year1
tarsus2
tarsus3
1year effects for  <dig> and  <dig> are shown relative to the predicted mean in  <dig>
estimated variance components of nestling house wren traits  under the animal model, heritabilities , and nest effects
trait/component
estimate ± se
full model lnl
reduced lnl
χ
2
p
vp
va
vnest
vresidual
h2
n2
vp
va
vnest
vresidual
h2
n2
vp
va
vnest
vresidual
h2
n2
vp
va
vnest
vresidual
h2
n2
neither tarsus length nor pha response were significantly heritable .
there was a significant difference between male and female nestlings in tarsus length, with males exhibiting longer tarsi .
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 .
our multivariate modelling with no random effects included showed that p contained significant among-trait covariance .
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 .
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  versus variance only matrix  lrt; χ26 =  <dig> , p =  <dig> ) and r .
among broods , nestling body condition was positively associated with both tarsus length and haematocrit, but was not correlated with pha response .
however, there was no evidence of significant covariance structure in g , and pairwise genetic correlations were characterised by large standard errors in all cases .table  <dig>
phenotypic covariance structure  and its additive genetic , nest  and residual  components
matrix
trait
tarsus
pha
hematocrit
condition
p
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
g
n
<dig>
<dig>
<dig>
<dig>
r
<dig>
<dig>
<dig>
<dig>
p was estimated from a multivariate model with no random effects while individual identification and nest were then included to decompose the covariance matrix into g, n and r. values below the diagonals are covariances  and values above the diagonals are the corresponding correlations.
bold values denote nominally significant parameter at p < <dig>  based |estimate| ≥2se.
discussion
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  <cit> .
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  <cit> , and other cross-fostering studies showing an absence of nest-of-origin effects on nestling haematocrit  <cit> .
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  <cit>  or spatial variation in the availability of insect prey  <cit> .
parasites might also contribute to the observed nest effects on these traits  <cit> , as nestlings in our house wren population are frequently infested with ectoparasitic mites  <cit> ; however, in an earlier study of our population, nestling mass was not affected by ectoparasite numbers  <cit>  and in a different house wren population, ectoparasites had no detectable effect on nestling tarsus and only a slight effect on mass  <cit> .
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  <dig>   <cit> .
however, christe et al.,  <cit>  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  <cit> , so our estimate of heritability of tarsus length in house wrens  is relatively conservative.
alternatively, low heritability of tarsus may be due to high residual variance rather than low additive genetic variation  <cit> .
a long-term study of savannah sparrows , 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   <cit> .
consistent with a number of other studies employing the animal model  <cit> , 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  <cit>  and that nestling condition in our house wren population is positively correlated with recruitment to the breeding population and subsequent reproductive success  <cit> .
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  <cit> .
in addition, selection may be acting primarily on the environmental component of variation in condition  <cit> .
this occurs if the positive correlation found between condition and fitness reflects a shared dependence on environmental factors rather than a causal relationship  <cit> .
pha response was not significantly heritable, consistent with other studies employing the animal model to estimate additive genetic variance  <cit> , as well as those relying on parent-offspring resemblance  <cit> .
a prominent exception to this pattern is a study of common kestrels  employing the animal model that revealed a surprising high heritability of pha response in fledglings   <cit> ; 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  <cit> , these have relied on full-sib comparisons that likely produce inflated estimates because of dominance variance and maternal effects  <cit> .
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  <cit> .
additional environmental factors might also include food availability, temperature, and parasite load  <cit> , and we also cannot rule out the possibility of parental effects arising from differential provisioning of broods  <cit>  or transgenerational priming of immunity arising from the transfer of maternal antibodies in eggs  <cit> .
given the complex nature of the pha response, which involves both innate and adaptive components of the immune system  <cit> , martin et al.
<cit>  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  <cit> .
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  <cit> .
a number of environmental factors can influence nestling haematocrit, among them parasite load  <cit> , food availability  <cit> , and temperature  <cit> .
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  <cit> .
in our house wren population, nestlings are frequently infested with two species of ectoparasitic mites  <cit>  and their density varies greatly among nests  <cit> .
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  <cit> .
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  <cit> .
although we did not measure provisioning of nestlings in the present study, provisioning rates within and among pairs can vary widely  <cit> , 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  <cit>  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  <cit>  or parent-offspring resemblance  <cit>  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 .
in support of this possibility, a study of a captive population of zebra finches  based on a pedigree of similar size as the one considered here  <cit> , reported a heritability of  <dig>  ±  <dig>   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  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  <cit> .
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  <cit> , 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  <cit> .
although oxygen transport generally increases with increasing haematocrit, increases in blood viscosity beyond a certain threshold may actually hinder oxygen transport  <cit> .
indeed, experimental studies in mice have shown that oxygen uptake and physical endurance are maximized at intermediate values of haematocrit  <cit> .
CONCLUSIONS
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  <cit> .
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 .
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  <cit> .
availability of supporting data
the data supporting the results of this article are available in the dryad digital repository doi: <dig> /dryad.jk2m <dig>  <cit>  .
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.
we thank the 2004– <dig> wren crews for field assistance and the parklands foundation  and the sears and butler families for the use of their properties.
financial support was provided by nsf grants gk12- <dig>  ibn- <dig>  ios- <dig> and ios-1118160; nih grant r15hd076308-01; a visiting professorship from the leverhulme trust ; the school of biological sciences, illinois state university; a bbsrc david phillips research fellowship ; and student-research grants from the beta lambda chapter of the phi sigma biological sciences honor society .

 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. 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 or disease resistance. 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. 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. condition measured in this manner is assumed to reflect lipid reserves, a major form of energy storage in birds that is used to fuel the various processes that promote individual fitness. indeed, body-mass based measures appear to be tightly linked with fitness in a number of vertebrate species. 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. 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. 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 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. although the conventional expectation has been that birds in good condition should be capable of mounting more robust immune responses, this expectation has been met in some studies, but not in others. 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. because haematocrit is directly related to oxygen uptake, it is typically regarded as a measure of physiological performance, but its value as an indicator of condition or health state remains uncertain. 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, 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. 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. 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. 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. however, a number of studies of wild bird populations have reported significant heritable variation in each of these traits. 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. here, we examine the genetic and environmental variation in condition, cutaneous immune responsiveness, and haematocrit over a three-year period 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. methods study animals house wrens are small, 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– eggs. about half of the females that successfully rear their first brood produce a second. 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. additional information concerning the breeding biology of house wrens is provided by johnson. we conducted this study during the 2004– breeding seasons on the mackinaw study site in a second-growth, deciduous forest bordering the mackinaw river. nestboxes of uniform construction were spaced 30 m apart along north–south transects separated by 60 m. all nestboxes were mounted on -m metal poles that had been greased or under which -cm diameter aluminium baffles had been mounted to discourage nest predators. the nestboxes used in this study were located in three semi-isolated neighbourhoods within the study area. all activities complied with the illinois state university institutional animal care and use committee and united states geological survey banding permit field procedures we visited all nestboxes every 1– 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 or, we ringed nestlings with a numbered aluminium band, weighed them on an electronic balance to the nearest  g, and measured their right tarsus length with dial callipers to the nearest  mm. at the same time, we obtained a blood sample 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  g for 60 s 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 by injecting 50 μl of sterile phosphate buffered saline containing phytohaemagglutinin at a concentration of 5 mg ml− into their prepatagium. phytohaemagglutinin, 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; but see. we used the change in wing-web thickness as a measure of cutaneous immune activity. wing-web thickness was measured with a mitutoyo thickness gauge. the repeatability of haematocrit and the pha response was and, respectively. determination of parentage we isolated dna from blood samples using a high-salt extraction protocol following bruford et al.. polymerase chain reaction amplifications were carried out in 15-μl volumes containing 200 μm dntps,  mm mgcl 1x pcr buffer ii, and  μm forward and reverse primers. we used a thermal profile that followed the touchdown protocol described in johnson et al.. forward primers were fluorescently labelled, and pcr products were analysed using a beckman coulter ceq genetic analysis system. we typed all samples at three loci: ta-c, mcyμ, and ltmr. when more resolution was needed, two additional loci, ta-a5- and ta-b4-, were used. we analysed allele data using cervus. no locus deviated significantly from hardy-weinberg equilibrium, with the exception of ta-a5- whose null allele frequency was estimated at. attendant females matched nestlings at all loci with rare exceptions 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; 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 based on matching loci was lower than, 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μ; ltmr; ta-c; ta-b4-; ta-a5-. for the three-locus set and five-locus set exclusion, probabilities were and, respectively. overall, the probability of false assignment for nestlings designated as within-pair was <. based on these methods, we were able to construct a pedigree comprised of house wrens, in which we were able to identify both parents of individuals and only the mother of an additional individuals. the pedigree spanned three generations, and included paternal and maternal identities. not all phenotypic traits were measured on all individuals; sample sizes and mean values for all traits are shown in table  unlike cross-fostering studies, our ability to separate nest effects 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, and, on our study area, nearly 60% of pairs switch mates across broods produced during the same breeding season. birds returning to the study area in subsequent years invariably pair with a different mate from the one with which they have paired the previous year. 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, and an even higher proportion of males sired young with more than one female. thus, our pedigree is well suited to separating common nest effects and genetic effects.table mean trait values for nestling house wrens included in the animal model analysis trait n mean sdrelative frequency of male house wrens males siring young with different numbers of females and the relative frequency of female house wrens producing young with different numbers of males. determination of sex sex was determined by amplifying sex-specific introns of the chd- gene. polymerase chain reactions were carried out with sexing primers 1237l and 1272h using a touchdown protocol as described in johnson et al., 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– 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. 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 variance provided pedigree data are available. we first fitted univariate models to estimate the additive variance for each trait. 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. similarly, we assume that residuals from the model are uncorrelated and normally distributed, with the residual variance arising from unmodelled environmental effects and/or measurement error. fixed effects included year, and sex 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 from our model of mass can be viewed as pertaining to body condition. for each trait we determined heritability as the ratio of va to vp and the corresponding proportion of variance explained by nest effects. 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 that we denote p. we then added individual genetic merit and nest effects back into the model to partition p into additive genetic, nest and residual matrices. statistical inference was based on likelihood-ratio tests 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: mix of χ and χ. for comparing multivariate models with differing random effects, we set the degrees of freedom to equal the number of additional variance parameters estimated in the more complex model.  fixed-effect estimates for nestling house wren traits are shown in table  and their associated heritabilities are provided in table  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. notwithstanding this large amount of environmental variation, both nestling haematocrit and condition were significantly heritable, with additive genetic variance accounting for about 15% of the phenotypic variation in each case. nestling condition differed across years, and male nestlings were of significantly higher body condition than female nestlings. there were no significant year effects or sex differences in nestling haematocrit.table univariate models showing fixed effect estimates for nestling house wren traits trait/effect level coefficient ± se df num df den f p year1 year1 year1 year1 tarsus2 tarsus3 1year effects for and are shown relative to the predicted mean in estimated variance components of nestling house wren traits under the animal model, heritabilities, and nest effects trait/component estimate ± se full model lnl reduced lnl χ 2 p vp va vnest vresidual h2 n2 vp va vnest vresidual h2 n2 vp va vnest vresidual h2 n2 vp va vnest vresidual h2 n2 neither tarsus length nor pha response were significantly heritable. there was a significant difference between male and female nestlings in tarsus length, with males exhibiting longer tarsi. 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. our multivariate modelling with no random effects included showed that p contained significant among-trait covariance. 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. 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 versus variance only matrix lrt; χ26 = , p =  ) and r. among broods, nestling body condition was positively associated with both tarsus length and haematocrit, but was not correlated with pha response. however, there was no evidence of significant covariance structure in g, and pairwise genetic correlations were characterised by large standard errors in all cases.table phenotypic covariance structure and its additive genetic, nest and residual components matrix trait tarsus pha hematocrit condition p g n r p was estimated from a multivariate model with no random effects while individual identification and nest were then included to decompose the covariance matrix into g, n and r. values below the diagonals are covariances and values above the diagonals are the corresponding correlations. bold values denote nominally significant parameter at p < based |estimate| ≥2se. discussion 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. 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, and other cross-fostering studies showing an absence of nest-of-origin effects on nestling haematocrit. 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 or spatial variation in the availability of insect prey. parasites might also contribute to the observed nest effects on these traits, as nestlings in our house wren population are frequently infested with ectoparasitic mites; however, in an earlier study of our population, nestling mass was not affected by ectoparasite numbers and in a different house wren population, ectoparasites had no detectable effect on nestling tarsus and only a slight effect on mass. 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. however, christe et al., 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, so our estimate of heritability of tarsus length in house wrens is relatively conservative. alternatively, low heritability of tarsus may be due to high residual variance rather than low additive genetic variation. a long-term study of savannah sparrows, 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. consistent with a number of other studies employing the animal model, 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 and that nestling condition in our house wren population is positively correlated with recruitment to the breeding population and subsequent reproductive success. 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. in addition, selection may be acting primarily on the environmental component of variation in condition. this occurs if the positive correlation found between condition and fitness reflects a shared dependence on environmental factors rather than a causal relationship. pha response was not significantly heritable, consistent with other studies employing the animal model to estimate additive genetic variance, as well as those relying on parent-offspring resemblance. a prominent exception to this pattern is a study of common kestrels employing the animal model that revealed a surprising high heritability of pha response in fledglings; 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, these have relied on full-sib comparisons that likely produce inflated estimates because of dominance variance and maternal effects. 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. additional environmental factors might also include food availability, temperature, and parasite load, and we also cannot rule out the possibility of parental effects arising from differential provisioning of broods or transgenerational priming of immunity arising from the transfer of maternal antibodies in eggs. given the complex nature of the pha response, which involves both innate and adaptive components of the immune system, martin et al. 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. 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. a number of environmental factors can influence nestling haematocrit, among them parasite load, food availability, and temperature. 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. in our house wren population, nestlings are frequently infested with two species of ectoparasitic mites and their density varies greatly among nests. 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. 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. although we did not measure provisioning of nestlings in the present study, provisioning rates within and among pairs can vary widely, 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 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 or parent-offspring resemblance 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. in support of this possibility, a study of a captive population of zebra finches based on a pedigree of similar size as the one considered here, reported a heritability of  ±  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 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. 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, 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. although oxygen transport generally increases with increasing haematocrit, increases in blood viscosity beyond a certain threshold may actually hinder oxygen transport. indeed, experimental studies in mice have shown that oxygen uptake and physical endurance are maximized at intermediate values of haematocrit.  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. 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. 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. availability of supporting data the data supporting the results of this article are available in the dryad digital repository doi: /dryad.jk2m. 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. we thank the 2004– wren crews for field assistance and the parklands foundation and the sears and butler families for the use of their properties. financial support was provided by nsf grants gk12- ibn- ios- and ios-1118160; nih grant r15hd076308-01; a visiting professorship from the leverhulme trust; the school of biological sciences, illinois state university; a bbsrc david phillips research fellowship; and student-research grants from the beta lambda chapter of the phi sigma biological sciences honor society.
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. in addition, haematocrit, the packed cell volume in relation to total blood volume, is often measured as an indicator of physiological performance. here, we examine the genetic and environmental variation in condition, cutaneous immunity, and haematocrit using an animal model based on a pedigree of approximately house wrens. 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. 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. 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. animal modelconditionhaematocritheritabilitygenetic variationimmunitylife-history theoryphatroglodytes aedonissue-copyright-statement© the author 2014 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. 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.
 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 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. keywords animal modelconditionhaematocritheritabilitygenetic variationimmunitylife-history theoryphatroglodytes aedonissue-copyright-statement© the author 2014 