The claim that human culture evolves through the differential adoption of cultural variants, in a manner analogous to the evolution of biological species, has been greeted with much resistance and confusion. Here we demonstrate that as compelling a case can now be made that cultural evolution has key Darwinian properties, as Darwin himself presented for biological evolution in The Origin of Species. Culture is shown to exhibit variation, competition, inheritance, and the accumulation of successive cultural modifications over time. Adaptation, convergence, and the loss or change of function can also be identified in culture. Just as Darwin knew nothing of genes or particulate inheritance, a case for Darwinian cultural evolution can be made irrespective of whether unitary cultural replicators exist or whether cultural transmission mechanisms are well understood.

Observed variations in rates of taxonomic diversification have been attributed to a range of factors including biological innovations, ecosystem restructuring, and environmental changes. Before inferring causality of any particular factor, however, it is critical to demonstrate that the observed variation in diversity is significantly greater than that expected from natural stochastic processes. Relative tests that assess whether observed asymmetry in species richness between sister taxa in monophyletic pairs is greater than would be expected under a symmetric model have been used widely in studies of rate heterogeneity and are particularly useful for groups in which paleontological data are problematic. Although one such test introduced by Slowinski and Guyer a decade ago has been applied to a wide range of clades and evolutionary questions, the statistical behavior of the test has not been examined extensively, particularly when used with Fisher's procedure for combining probabilities to analyze data from multiple independent taxon pairs. Here, certain pragmatic difficulties with the Slowinski-Guyer test are described, further details of the development of a recently introduced likelihood-based relative rates test are presented, and standard simulation procedures are used to assess the behavior of the two tests in a range of situations to determine: (1) the accuracy of the tests' nominal Type I error rate; (2) the statistical power of the tests; (3) the sensitivity of the tests to inclusion of taxon pairs with few species; (4) the behavior of the tests with datasets comprised of few taxon pairs; and (5) the sensitivity of the tests to certain violations of the null model assumptions. Our results indicate that in most biologically plausible scenarios, the likelihood-based test has superior statistical properties in terms of both Type I error rate and power, and we found no scenario in which the Slowinski-Guyer test was distinctly superior, although the degree of the discrepancy varies among the different scenarios. The Slowinski-Guyer test tends to be much more conservative (i.e., very disinclined to reject the null hypothesis) in datasets with many small pairs. In most situations, the performance of both the likelihood-based test and particularly the Slowinski-Guyer test improve when pairs with few species are excluded from the computation, although this is balanced against a decline in the tests' power and accuracy as fewer pairs are included in the dataset. The performance of both tests is quite poor when they are applied to datasets in which the taxon sizes do not conform to the distribution implied by the usual null model. Thus, results of analyses of taxonomic rate heterogeneity using the Slowinski-Guyer test can be misleading because the test's ability to reject the null hypothesis (equal rates) when true is often inaccurate and its ability to reject the null hypothesis when the alternative (unequal rates) is true is poor, particularly when small taxon pairs are included. Although not always perfect, the likelihood-based test provides a more accurate and powerful alternative as a relative rates test.

Sexual selection in the form of sperm competition is a major explanation for small size of male gametes. Can sexual selection in polyandrous species with reversed sex roles also lead to reduced female gamete size? Comparative studies show that egg size in birds tends to decrease as a lineage evolves social polyandry. Here, a quantitative genetic model predicts that female scrambles over mates lead to evolution of reduced female gamete size. Increased female mating success drives the evolution of smaller eggs, which take less time to produce, until balanced by lowered offspring survival. Mean egg size is usually reduced and polyandry increased by increasing sex ratio (male bias) and maximum possible number of mates. Polyandry also increases with the asynchrony (variance) in female breeding start. Opportunity for sexual selection increases with the maximum number of mates but decreases with increasing sex ratio. It is well known that parental investment can affect sexual selection. The model suggests that the influence is mutual: owing to a coevolutionary feedback loop, sexual selection in females also shapes initial parental investment by reducing egg size. Feedback between sexual selection and parental investment may be common.

Learning processes potentially play a role in speciation but are often ignored in speciation models. Learning may, for instance, play a role when a new niche is being colonized, because the learning of niche features may cause niche-specific assortative mating and a tendency to produce young in this niche. Several animal species learn about their environmental features that may be important in finding or attracting mates. We use a gene-culture coevolutionary model to look into the effect of such learning on the colonization of new niches and on the genetic divergence between groups using different niches, which are steps necessary in achieving speciation. We assume that density is regulated separately in each of the two niches and that the viability of an individual depends on its genotype as well as on which niche it exploits. Our results show that genetic adaptation to the new niche is enhanced by a high female fecundity and a low viability selection against heterozygotes. Furthermore, when initial colonization (without genetic adaptation) fails, genetic divergence is more difficult when the mating preference is stronger. In contrast, when colonization without genetic adaptation is successful, a stronger mating preference makes genetic divergence easier. An increase in the number of egg-laying mistakes by females can have a positive or negative effect on the success of genetic adaptation depending on other parameters. We show that genetic divergence can be prevented by a niche shift, which can occur only if viabilities in the two niches are asymmetrical.

The genetic basis of fluctuating asymmetry (FA), a measure of random deviations from perfect bilateral symmetry, has been the subject of much recent work. In this paper we compare two perspectives on the quantitative genetic analysis of FA and directional asymmetry (DA). We call these two approaches the character-state model and the environmental responsiveness model. In the former approach, the right and left sides are viewed as separate traits whose genetic coupling is manifested by the genetic correlation. This model leads to the relationship, h²_DA = h²[(1 − r_A)/(1 − r_P)], where h² is the heritability of each component trait (assumed to be the same), r_A and r_P are the genetic and phenotypic correlations between traits, respectively. Simulation shows that, under this model, the heritability of FA is considerably less than that of DA, except when heritabilities are very close to zero. The environmental responsiveness model permits genetic variance in FA even when the genetic correlation between traits is 1. Simulation shows that under this model the heritability of FA can be uncoupled from that of DA. The additive and nonadditive components of the component (right and left) traits, their DA and FA values are estimated using a diallel cross of seven inbred lines of the sand cricket, Gryllus firmus. Four leg measurements were made and both the individual DA and FA values and the compound measures DASUM and CFA estimated. The heritabilities of the compound measures are slightly larger than the individual estimates. Dominance variance is observed in the individual traits but predicted to be an even smaller component of the phenotypic variance than the additive genetic variance. The estimated values confirm this, although a previous study has demonstrated that dominance variance is present. Because the heritabilities of FA are generally larger than those of DA, which never exceed 0.02, the environmental responsiveness model is more consistent with the data than the character-state model. A review of other data suggests that both sources of variation might be found in some species.

Both genetic differentiation and phenotypic plasticity might be expected to affect the location of geographic range limits. Co-gradient variation (CoGV), plasticity that is congruent with genetic differentiation, may enhance performance at range margins, whereas its opposite, counter-gradient variation (CnGV) may hinder performance. Here we report findings of reciprocal transplant experiments intended to tease apart the roles of differentiation and plasticity in producing phenotypic variation across a geographic border between two plant subspecies. Clarkia xantiana ssp. xantiana and C. xantiana ssp. parviflora are California-endemic annuals that replace each other along a west–east gradient of declining precipitation. We analyzed variation in floral traits, phenological traits, and vegetative morphological and developmental traits by sowing seeds of 18 populations (six of ssp. xantiana and 12 of ssp. parviflora) at three sites (one in each subspecies' exclusive range and one in the subspecies' contact zone), in two growing seasons (an exceptionally wet El Niño winter and a much drier La Niña winter). Significant genetic differences between subspecies appeared in 11 of 12 traits, and differences were of the same sign as in nature. These findings are consistent with the hypothesis that selection is responsible for subspecies differences. Geographic variation within subspecies over part of the spatial gradient mirrored between-subspecies differences present at a larger scale. All traits showed significant plasticity in response to spatial and temporal environmental variation. Plasticity patterns ranged from spatial and temporal CoGV (e.g., in node of first flower), to spatial CnGV (e.g., in flowering time), to patterns that were neither CoGV nor CnGV (the majority of traits). Instances of CoGV may reflect adaptive plasticity and may serve to increase performance under year-to-year environmental variation and at sites near the subspecies border. However, the presence of spatial CnGV in some critical traits suggests that subspecies ranges may also be constrained by patterns of plasticity.

Most theoretical treatments of the evolutionary ecology of offspring size assume a simple and direct effect of investment per offspring on offspring fitness. In this paper I experimentally determine the relationship between seed mass and several main fitness components of the oak Quercus ilex, to estimate phenotypic selection acting on seed mass during the early life cycle and to discover any potential selective conflicts occurring between different stages from dispersal to establishment. I found a positive effect of acorn size on most fitness components related to seedling establishment. Large size increased germination rate and seedling survival, accelerated germination timing, and enhanced seedling growth. Nevertheless, there was also a direct negative effect of acorn size on survival to predation, because large acorns were highly preferred by the main postdispersal seed predators at the study site, wild boars and wood mice. Because of the low probability of escape from predation, the fitness of large acorns estimated on this component was significantly lower than the fitness of smaller acorns. Therefore, seed size affected fitness in two different ways, yielding opposing and conflicting selective forces. These findings suggest that the general assumption that offspring fitness is a fixed positive function of seed size needs to be reconsidered for some systems. The existence of conflicting selection might explain the occurrence of an optimal seed size in some plant species without invoking a seed number-size trade-off.

Pigmentation in the freshwater isopod Asellus aquaticus (Crustacea) differed between habitats in two Swedish lakes. In both lakes, isopods had lighter pigmentation in stands of submerged vegetation, consisting of stoneworts (Chara spp.), than in nearby stands of reed (Phragmites australis). Experimental crossings of light and dark isopods in a common environment showed that pigmentation had a genetic basis and that genetic variance was additive. Environmental effects of diet or chromatophore adjustment to the background had minor influence on pigmentation, as shown by laboratory rearing of isopods on stonewort or reed substrates, as well as analyses of stable isotope ratios for isopods collected in the field. In both study lakes, the average phenotype became lighter with time (across generations) in recently established stonewort stands. Taken together, these results indicate that altered phenotype pigmentation result from evolutionary responses to local differences in natural selection. Based on the assumption of two generations per year, the evolutionary rate of change in pigmentation was 0.08 standard deviations per generation (haldanes) over 20 generations in one lake and 0.22 haldanes over two generations in the other lake. This genetic change occurred during an episode of population growth in a novel habitat, a situation known to promote adaptive evolution. In addition, stonewort stands constitute large and persistent patches, characteristics that tend to preserve local adaptations produced by natural selection. Results from studies on selective forces behind the adaptive divergence suggest that selective predation from visually oriented predators is a possible selective agent. We found no indications of phenotype-specific movements between habitats. Mating within stonewort stands was random with respect to pigmentation, but on a whole-lake scale it is likely that mating is assortative, as a result of local differences in phenotype distribution.

Repeated patterns among biological communities suggest similar evolutionary and ecological forces are acting on the communities. Conversely, the lack of such patterns suggests that similar forces are absent or additional ones are present. Coevolution between a seed predator, the red crossbill (Loxia curvirostra complex), and lodgepole pine (Pinus contorta var. latifolia) exemplifies the ecological and evolutionary predictions for coevolving systems. In the absence of another seed predator and preemptive competitor (pine squirrels Tamiasciurus hudsonicus), natural selection by crossbills results in the evolution of larger cones with thicker distal scales, while relaxation of selection by squirrels results in the evolution of cones with more seeds and a greater ratio of seed mass to cone mass. However, in one range, the Little Rocky Mountains, distal scale thickness has diverged as expected but cone size has not. In these mountains seed predation by lodgepole pine cone borer moths (Eucosma recissoriana) was about 10 times greater than in other ranges lacking squirrels. We quantified moth predation and cone traits and found that moths select for smaller cones with fewer seeds. Thus, selection by moths in the Little Rocky Mountains counters both selection by crossbills for large cone size and relaxation of selection by squirrels favoring more seeds per cone and accounts for the relatively small and few-seeded cones in these mountains. It is also apparent that selection by crossbills changes seed defenses in a manner that favors seed predation by moths, whereas selection by squirrels likely reduces such predation. These results demonstrate the importance of considering the evolutionary consequences of community context in locally evolved (coevolved) traits and interactions.

Population differentiation often reflects a balance between divergent natural selection and the opportunity for homogenizing gene flow to erode the effects of selection. However, during ecological speciation, trait divergence results in reproductive isolation and becomes a cause, rather than a consequence, of reductions in gene flow. To assess both the causes and the reproductive consequences of morphological differentiation, we examined morphological divergence and sexual isolation among 17 populations of Timema cristinae walking-sticks. Individuals from populations adapted to using Adenostoma as a host plant tended to exhibit smaller overall body size, wide heads, and short legs relative to individuals using Ceonothus as a host. However, there was also significant variation in morphology among populations within host-plant species. Mean trait values for each single population could be reliably predicted based upon host-plant used and the potential for homogenizing gene flow, inferred from the size of the neighboring population using the alternate host and mitochondrial DNA estimates of gene flow. Morphology did not influence the probability of copulation in between-population mating trials. Thus, morphological divergence is facilitated by reductions in gene flow, but does not cause reductions in gene flow via the evolution of sexual isolation. Combined with rearing data indicating that size and shape have a partial genetic basis, evidence for parallel origins of the host-associated forms, and inferences from functional morphology, these results indicate that morphological divergence in T. cristinae reflects a balance between the effects of host-specific natural selection and gene flow. Our findings illustrate how data on mating preferences can help determine the causal associations between trait divergence and levels of gene flow.

Climate oscillations of the Quaternary drove the repeated expansion and contraction of ecosystems. Alpine organisms were probably isolated in sky island refugia during warm interglacials, such as now, and expanded their range by migrating down-slope during glacial periods. We used population genetic and phylogenetic approaches to infer how paleoclimatic events influenced the distribution of genetic variation in the predominantly alpine butterfly Parnassius smintheus. We sequenced a 789 bp region of cytochrome oxidase I for 385 individuals from 20 locations throughout the Rocky Mountains, ranging from southern Colorado to northern Montana. Analyses revealed at lease two centers of diversity in the northern and southern Rocky Mountains and strong population structure. Nested clade analysis suggested that the species experienced repeated cycles of population expansion and fragmentation. The estimated ages of these events, assuming a molecular clock, corresponded with paleoclimatic data on habitat expansion and contraction over the past 400,000 years. We propose that alpine butterflies persisted in an archipelago of isolated sky islands during interglacials and that populations expanded and became more connected during cold glacial periods. An archipelago model implies that the effects of genetic drift and selection varied among populations, depending on their latitude, area, and local environment. Alpine organisms are sensitive indicators of climate change and their history can be used to predict how high-elevation ecosystems might respond to further climate warming.

Genetic changes conferring adaptation to a new environment may induce a fitness cost in the previous environment. Although this prediction has been verified in laboratory conditions, few studies have tried to document this cost directly in natural populations. Here, we evaluated the pleiotropic effects of insecticide resistance on putative fitness components of the mosquito Culex pipiens. Experiments using different larval densities were performed during the summer in two natural breeding sites. Two loci that possess alleles conferring organophosphate (OP) resistance were considered: ace-1 coding for an acetylcholinesterase (AChE1, the OP target) and Ester, a “super locus” including two closely linked loci coding for esterases A and B. Resistance ace-1 alleles coding for a modified AChE1 were associated with a longer development time and shorter wing length. The pleiotropic effects of two resistance alleles Ester¹ and Ester⁴ coding for the overproduced esterases A1 and A4-B4, respectively, were more variable. Both A1 and A4-B4 reduced wing length, although only A1 was associated with a longer preimaginal stage. The fluctuating asymmetry (FA) of the wing did not respond to the presence or to the interaction of resistance alleles at the two loci at any of the density levels tested. Conversely, the FA of one wing section decreased when larval density increased. This may be the consequence of selection against less developmentally stable individuals. The results are discussed in relation to the local evolution of insecticide resistance genes.

In the conditional strategy model, divergence in reproductive phenotypes depends on whether the individual's condition is above or below a genetically determined threshold. The relative contribution of the genetic and environmental components that lead to the expression of a reproductive tactic by an individual is not well understood. In the present field study, we determined when condition diverged between males that develop the mature parr phenotype and those that do not in Atlantic salmon (Salmo salar). We also investigated the uniformity of the threshold value in the population. We sampled mature parr and immature males at age one, of the same population at six different sites for four consecutive years. Our study provides an example of the interaction of genotype and environment on the expression of a reproductive tactic. Size was significantly greater for future mature parr than for future immature males as early as 20 days after hatching (emergence), suggesting that there may be a parental effect component in the tactic adopted, since no exogenous feeding takes place before this time. Size advantage at emergence was maintained through the next spring at age one to different degrees depending on the year, thus suggesting the presence of an environmental component of tactic expression. Our results support the contention that within the conditional strategy, the environment faced by a male and his condition at the moment of reproduction consistently predicts neither the environment faced by his offspring nor the fitness they will obtain by expressing the same tactic as their father. Furthermore, higher mean size at a site did not always translate into a higher proportion of mature parr, therefore supporting the hypothesis that thresholds vary across habitats within the same population.

This study uses a combined methodological approach including phylogenetic, phylogeographic, and demographic analyses to understand the evolutionary history of the northern leopard frog, Rana pipiens. We tested hypotheses concerning how (or if) known geological events and key features of the species biology influenced the contemporary geographic and genetic distribution of R. pipiens. We assayed mitochondrial DNA variation from 389 individuals within 35 populations located throughout the species range. Our a priori expectations for patterns and processes influencing the current genetic structure of R. pipiens were supported by the data. However, our analyses revealed specific aspects of R. pipiens evolutionary history that were unexpected. The phylogenetic analysis indicated that R. pipiens is split into populations containing discrete eastern or western haplotypes, with the Mississippi River and Great Lakes region dividing the geographic ranges. Nested clade analysis indicated that the biological process most often invoked to explain the pattern of haplotype position is restricted gene flow with isolation by distance. Demographic analyses showed evidence of both historical bottlenecks and population expansions. Surprisingly, the genetic evidence indicated that the western haplotypes had significantly reduced levels of genetic diversity relative to the eastern haplotypes and that major range expansions occurred in both regions well before the most recent glacial retreat. This study provides a detailed history of how a widespread terrestrial vertebrate responded to episodic Pleistocene glacial events in North America. Moreover, this study illustrates how complementary methods of data analysis can be used to disentangle recent and ancient effects on the genetic structure of a species.

Prior studies have shown that macrogeographic gradients in temperature associated with latitude and altitude can lead to countergradient patterns of variation in a number of taxa: individuals from colder environments are known to grow or develop faster than their conspecifics from warmer environments when placed in a common setting. In this study, I hypothesized that countergradient variation also is important at microgeographic scales. The wood frog, Rana sylvatica, breeds in open-canopied, temporary wetlands as well as those heavily shaded by vegetation. Shading leads to cooler thermal environments that are associated with embryonic development rates as much as 50% slower than those in unshaded wetlands. Wetlands with contrasting canopy environments are often found within tens or hundreds of meters of each other. In a common garden experiment, embryos from nearby natural wetlands displayed countergradient variation: individuals collected from shaded wetlands developed up to 12% faster than those collected from relatively unshaded wetlands. The results of this study suggest that the concept of countergradient variation may be extended to small scales of space. In addition, the rate and scale of vegetation dynamics (the agent of wetland shading) imply that divergence in development among residents of nearby wetlands may be relatively rapid, on the order of decades.

Colonization of novel environments creates new selection pressures. Sexually selected traits are affected by the physical and social environment and should be especially susceptible to change, but this has rarely been studied. In southern California, dark-eyed juncos, (Junco hyemalis) naturally breed in mixed-coniferous temperate forests, typically from 1500 m to 3000 m in elevation. In the early 1980s, a small population became established in a coastal habitat, the University of California, San Diego campus, which has a mild, Mediterranean climate. I show that a sexually and socially selected signaling trait—the amount of white in the tail—has declined by approximately 22% as compared to mountain juncos. I address three main factors that could explain the difference between mountain and coastal juncos: phenotypic plasticity, genetic drift, and selection. Results indicate that the first two can be ruled out as the sole cause of the plumage change, which implies that selection contributed to the genetic differentiation from the mountain population. The estimated rate of evolution is about 0.2 haldanes, comparable with rates of change in systems where individuals have been artificially introduced into new environments (e.g., guppies and Drosophila). This is the first study to demonstrate evolution of a sexually selected trait after only several generations resulting from a natural invasion into a novel environment.

Diapause, the temporary cessation of development at an early life-history stage, is widespread among animals and plants. The range of taxa exhibiting various forms of diapause indicates its enormous ecological significance and highlights its value as a model for examining life-history trait evolution. However, despite the impact of diapause on species ecology, there is little understanding of its adaptive value in many groups. Furthermore, the relative roles of phylogeny and ecology in determining the contemporary expression of the trait remain unresolved. Delayed implantation (DI) is a type of diapause found in several orders of mammals. It is particularly prevalent in the Mustelidae, with mustelids making up more than half of all mammals known to exhibit DI. This taxon is thus ideal for examining life-history predictors of DI and investigating the mode of evolution. Both maximum likelihood and maximum parsimony methods of ancestral state reconstruction indicated DI to be plesiomorphic in the mustelids, although multiple state changes are required to explain its contemporary distribution. After controlling for phylogeny, species with and without DI could be discriminated using just three variables: longevity, maximum latitude of the geographical distribution, and a term describing maternal investment. Our analyses supported the hypothesis that DI is more prevalent in seasonal climates. We also showed that longer-lived species are more likely to exhibit DI, suggesting a time cost to the trait. We found no correlate for the highly variable duration of DI, which remains unexplained. Although ecological factors can predict the distribution of DI in modern mustelids, phylogenetic constraint is likely to play an important role.

In the alpine valley of Valtellina there are two Robertsonian chromosomal races of house mouse, the Poschiavo (POS: 2n = 24–26) characterized by metacentric 8.12 and acrocentrics 2 and 10 and the Upper Valtellina (UV: 2n = 22–24) characterized by metacentrics 2.8 and 10.12. The races inhabit separate villages in the valley except in Sommacologna and Sondalo, where they both occur together with hybrids. A total of 179 mice from 16 villages were typed at 13 microsatellite loci. Seven of these loci were localized close to the centromeres of chromosomes 10 and 12, with the prediction that these regions on the race-specific chromosomes would be the most likely to experience a barrier to gene flow. The remaining six loci were localized at the telomeres of chromosomes 10 and 12 and at the centromeres of chromosomes that do not differ between the races. Substantial differences in allelic frequencies were found between the villages with POS and UV races at five of the loci at the centromeres of chromosomes 10 and 12 but at none of the other loci. These differences were not found to distinguish the two races in Sommacologna and Sondalo. Therefore, the centromeric regions of race-specific chromosomes do appear to experience a barrier to gene flow, although this can break down under intense interbreeding between the races. These results are considered in the context of Harrison's (1990) concept of the semipermeability of hybrid zones to gene exchange and in relation to parapatric speciation.

Evolutionary biologists have often suggested that ecology is important in speciation, in that natural selection may drive adaptive divergence between lineages that inhabit different environments. I suggest that it is the tendency of lineages to maintain their ancestral ecological niche (phylogenetic niche conservatism) and their failure to adapt to new environments which frequently isolates incipient species and begins the process of speciation. Niche conservatism may be an important and widespread component of allopatric speciation but is largely unstudied. The perspective outlined here suggests roles for key microevolutionary processes (i.e., natural selection, adaptation) that are strikingly different from those proposed in previous literature on ecology and speciation. Yet, this perspective is complementary to the traditional view because it focuses on a different temporal stage of the speciation process.

Little is known about what controls effective sizes and migration rates among parasite populations. Such data are important given the medical, veterinary, and economic (e.g., fisheries) impacts of many parasites. The autogenic-allogenic hypothesis, which describes ecological patterns of parasite distribution, provided the foundation on which we studied the effects of life cycles on the distribution of genetic variation within and among parasite populations. The hypothesis states that parasites cycling only in freshwater hosts (autogenic life cycle) will be more limited in their dispersal ability among aquatic habitats than parasites cycling through freshwater and terrestrial hosts (allogenic life cycle). By extending this hypothesis to the level of intraspecific genetic variation, we examined the effects of host dispersal on parasite gene flow. Our a priori prediction was that for a given geographic range, autogenic parasites would have lower gene flow among subpopulations. We compared intraspecific mitochondrial DNA variation for three described species of trematodes that infect salmonid fishes. As predicted, autogenic species had much more highly structured populations and much lower gene flow among subpopulations than an allogenic species sampled from the same locations. In addition, a cryptic species was identified for one of the autogenic trematodes. These results show how variation in life cycles can shape parasite evolution by predisposing them to vastly different genetic structures. Thus, we propose that knowledge of parasite life cycles will help predict important evolutionary processes such as speciation, coevolution, and the spread of drug resistance.

We tested the effect of detrimental environmental conditions during growth on the heritability of chick body size in a wild population of blue tits (Parus caeruleus) highly parasitized by blowfly larvae. During nine years, we experimentally induced deparasitized broods, whereas unmanipulated control broods remained heavily infested by two species of Protocalliphora ectoparasites. The heritability of tarsus length was significantly higher in deparasitized broods than control broods, due in part to a very low common brood environment effect in deparasitized broods. We also found evidence for significant genotype-by-environment interactions, which further reflected the effect of the ecological conditions on the expression of additive genetic effects and could represent an additional constraint on the evolution of tarsus length. To our knowledge, this experiment provides the first evidence of host quantitative genetics being influenced by parasitism, and illustrates the potential for parasitism to constrain an evolutionary response to selection.