Quantitative genetics of growth and cryptic evolution of body size in an island population

While evolution occurs when selection acts on a heritable trait, empirical studies of natural systems have frequently reported phenotypic stasis under these conditions. We performed quantitative genetic analyses of weight and hindleg length in a free-living population of Soay sheep (Ovis aries) to test whether genetic constraints can explain previously reported stasis in body size despite evidence for strong positive directional selection. Genetic, maternal and environmental covariance structures were estimated across ontogeny using random regression animal models. Heritability increased with age for weight and hindleg length, though both measures of size were highly heritable across ontogeny. Genetic correlations among ages were generally strong and uniformly positive, and the covariance structures were also highly integrated across ontogeny. Consequently, we found no constraint to the evolution of larger size itself. Rather we expect size at all ages to increase in response to positive selection acting at any age. Consistent with expectation, predicted breeding values for age-specific size traits have increased over a twenty-year period, while maternal performance for offspring size has declined. Re-examination of the phenotypic data confirmed that sheep are not getting larger, but also showed that there are significant negative trends in size at all ages. The genetic evolution is therefore cryptic, with the response to selection presumably being masked at the phenotypic level by a plastic response to changing environmental conditions. Density-dependence, coupled with systematically increasing population size, may contribute to declining body size but is insufficient to completely explain it. Our results demonstrate that an increased understanding of the genetic basis of quantitative traits, and of how plasticity and microevolution can occur simultaneously, is necessary for developing predictive models of phenotypic change in nature.     

Natural Selection on Individual Variation in Tolerance of Gastrointestinal Nematode Infection

Hosts may mitigate the impact of parasites by two broad strategies: resistance, which limits parasite burden, and tolerance, which limits the fitness or health cost of increasing parasite burden. The degree and causes of variation in both resistance and tolerance are expected to influence host–parasite evolutionary and epidemiological dynamics and inform disease management, yet very little empirical work has addressed tolerance in wild vertebrates. Here, we applied random regression models to longitudinal data from an unmanaged population of Soay sheep to estimate individual tolerance, defined as the rate of decline in body weight with increasing burden of highly prevalent gastrointestinal nematode parasites. On average, individuals lost weight as parasite burden increased, but whereas some lost weight slowly as burden increased (exhibiting high tolerance), other individuals lost weight significantly more rapidly (exhibiting low tolerance). We then investigated associations between tolerance and fitness using selection gradients that accounted for selection on correlated traits, including body weight. We found evidence for positive phenotypic selection on tolerance: on average, individuals who lost weight more slowly with increasing parasite burden had higher lifetime breeding success. This variation did not have an additive genetic basis. These results reveal that selection on tolerance operates under natural conditions. They also support theoretical predictions for the erosion of additive genetic variance of traits under strong directional selection and fixation of genes conferring tolerance. Our findings provide the first evidence of selection on individual tolerance of infection in animals and suggest practical applications in animal and human disease management in the face of highly prevalent parasites.     

Repeated selection of morphometric traits in the Soay sheep on St Kilda

1. Long-term studies allow the outcomes of repeated selection events to be monitored. Here, we investigate phenotypic selection in successive winter mortality events in the Soay sheep of St Kilda, Scotland, between 1985 and 1996. Selection of three quantitative morphometric traits, body weight, hindleg length and incisor arcade breadth, was investigated in different sectors of the population.
2. Evidence from fitness differentials of positive directional selection for large size was repeatedly found in lambs and adult females. Selection in the opposing direction was only found in one year in lambs.
3. Selection gradients showed that in most years when significant selection occurred, body weight was the focus of direct selection, whereas selection of hindleg length and incisor breadth was indirect, arising from their correlation with body weight.
4. Selection was strongest in years of low over-winter survival and almost absent in years when survival was high. Intensity of selection was greatest in lambs, emphasizing the differences in selection pressure experienced by different sectors of the population, in addition to the temporal variation in selection pressure due to population density and environmental conditions.
5. Despite repeated positive selection of body weight, no evidence of a change in the population mean was found over the course of the study.     

Heterogeneity of genetic architecture of body size traits in a free‐living population

Knowledge of the underlying genetic architecture of quantitative traits could aid in understanding how they evolve. In wild populations, it is still largely unknown whether complex traits are polygenic or influenced by few loci with major effect, due to often small sample sizes and low resolution of marker panels. Here, we examine the genetic architecture of five adult body size traits in a free‐living population of Soay sheep on St Kilda using 37 037 polymorphic SNPs. Two traits (jaw and weight) show classical signs of a polygenic trait: the proportion of variance explained by a chromosome was proportional to its length, multiple chromosomes and genomic regions explained significant amounts of phenotypic variance, but no SNPs were associated with trait variance when using GWAS. In comparison, genetic variance for leg length traits (foreleg, hindleg and metacarpal) was disproportionately explained by two SNPs on chromosomes 16 (s23172.1) and 19 (s74894.1), which each explained >10% of the additive genetic variance. After controlling for environmental differences, females heterozygous for s74894.1 produced more lambs and recruits during their lifetime than females homozygous for the common allele conferring long legs. We also demonstrate that alleles conferring shorter legs have likely entered the population through a historic admixture event with the Dunface sheep. In summary, we show that different proxies for body size can have very different genetic architecture and that dense SNP helps in understanding both the mode of selection and the evolutionary history at loci underlying quantitative traits in natural populations.     

Mapping quantitative trait Loci underlying fitness-related traits in a free-living sheep population

We searched for quantitative trait loci (QTL) underlying fitness-related traits in a free-living pedigree of 588 Soay sheep in which a genetic map using 251 markers with an average spacing of 15 cM had been established previously. Traits examined included birth date and weight, considered both as maternal and offspring traits, foreleg length, hindleg length, and body weight measured on animals in August and jaw length and metacarpal length measured on cleaned skeletal material. In some cases the data were split to consider different age classes separately, yielding a total of 15 traits studied. Genetic and environmental components of phenotypic variance were estimated for each trait and, for those traits showing nonzero heritability (N= 12), a QTL search was conducted by comparing a polygenic model with a model including a putative QTL. Support for a QTL at genome-wide significance was found on chromosome 11 for jaw length; suggestive QTL were found on chromosomes 2 and 5 (for birth date as a trait of the lamb), 8 (birth weight as a trait of the lamb), and 15 (adult hindleg length). We discuss the prospects for refining estimates of QTL position and effect size in the study population, and for QTL searches in free-living pedigrees in general.     

Natural selection on a measure of parasite resistance varies across ages and environmental conditions in a wild mammal

Parasites detrimentally affect host fitness, leading to expectations of positive selection on host parasite resistance. However, as immunity is costly, host fitness may be maximized at low, but nonzero, parasite infection intensities. These hypotheses are rarely tested on natural variation in free-living populations. We investigated selection on a measure of host parasite resistance in a naturally regulated Soay sheep population using a longitudinal data set and found negative correlations between parasite infection intensity and annual fitness in lambs, male yearlings and adult females. However, having accounted for confounding effects of body weight, the effect was only significant in lambs. Associations between fitness and parasite resistance were environment-dependent, being strong during low-mortality winters, but negligible during harsher high-mortality winters. There was no evidence for stabilizing selection. Our findings reveal processes that may shape variation in parasite resistance in natural populations and illustrate the importance of accounting for correlated traits in selection analysis.     

Estimating variance components and heritabilities in the wild: a case study using the ‘animal model’ approach

Using a genealogy containing over 1800 dams and nearly 400 sires (estimated by genetic paternity techniques), combined with maximum likelihood procedures and an ‘animal model’, we have estimated the heritabilities, genetic correlations and variance components of three morphometric traits in the Soay sheep (Ovis aries) on St Kilda, Scotland. This approach allows heritabilities to be estimated in natural populations that violate the assumptions of offspring–parent regression methods. Maternal (or paternal) effects can also be estimated under natural conditions. We demonstrate that all the traits, body weight, hind leg length and incisor arcade breadth, have low but significant heritabilities. Body weight, the trait that experiences the strongest selection, had the lowest heritability but the highest additive genetic coefficient of variation. An evolutionary response to selection is predicted. When maternal effects were not taken into consideration heritabilities were over‐estimated, although this effect was only significant in female offspring.     

Survival costs of reproduction are mediated by parasite infection in wild Soay sheep

A trade‐off between current and future fitness potentially explains variation in life‐history strategies. A proposed mechanism behind this is parasite‐mediated reproductive costs: individuals that allocate more resources to reproduction have fewer to allocate to defence against parasites, reducing future fitness. We examined how reproduction influenced faecal egg counts (FEC) of strongyle nematodes using data collected between 1989 and 2008 from a wild population of Soay sheep in the St. Kilda archipelago, Scotland (741 individuals). Increased reproduction was associated with increased FEC during the lambing season: females that gave birth, and particularly those that weaned a lamb, had higher FEC than females that failed to reproduce. Structural equation modelling revealed future reproductive costs: a positive effect of reproduction on spring FEC and a negative effect on summer body weight were negatively associated with overwinter survival. Overall, we provide evidence that parasite resistance and body weight are important mediators of survival costs of reproduction.     

Heritable variation in resistance to gastro-intestinal nematodes in an unmanaged mammal population.

The impact of parasites on natural populations has received considerable attention from evolutionary biologists in recent years. Central to a number of theoretical developments during this period is the assumption of additive genetic variation in resistance to parasites. However, very few studies have estimated the heritability of parasite resistance under field conditions, and those that have are mainly restricted to birds and their ectoparasites. In this paper, to our knowledge, we show for the first time in a free-ranging mammal population, Soay sheep (Ovis aries) living on the islands of St Kilda, that there is significant heritable variation in resistance to gastrointestinal nematodes. This result is consistent with earlier studies on this population which have indicated locus-specific associations with parasite resistance. We discuss our results in the context of current studies examining heritable resistance to parasites in domestic sheep and the possible mechanisms of selective maintenance of genetic variation for resistance to gastrointestinal nematodes in the St Kilda Soay sheep population.     

Quantitative trait loci (QTL) mapping of resistance to strongyles and coccidia in the free-living Soay sheep (Ovis aries)

A genome-wide scan was performed to detect quantitative trait loci (QTL) for resistance to gastrointestinal parasites and ectoparasitic keds segregating in the free-living Soay sheep population on St. Kilda (UK). The mapping panel consisted of a single pedigree of 882 individuals of which 588 were genotyped. The Soay linkage map used for the scans comprised 251 markers covering the whole genome at average spacing of 15cM. The traits here investigated were the strongyle faecal egg count (FEC), the coccidia faecal oocyst count (FOC) and a count of keds (Melophagus ovinus). QTL mapping was performed by means of variance component analysis so that the genetic parameters of the study traits were also estimated and compared with previous studies in Soay and domestic sheep. Strongyle FEC and coccidia FOC showed moderate heritability (h(2)=0.26 and 0.22, respectively) in lambs but low heritability in adults (h(2)<0.10). Ked count appeared to have very low h(2) in both lambs and adults. Genome scans were performed for the traits with moderate heritability and two genomic regions reached the level of suggestive linkage for coccidia FOC in lambs (logarithm of the odds=2.68 and 2.21 on chromosomes 3 and X, respectively). We believe this is the first study to report a QTL search for parasite resistance in a free-living animal population and therefore may represent a useful reference for similar studies aimed at understanding the genetics of host-parasite co-evolution in the wild.     

Variation in the peacock’s train shows a genetic component

Female peafowl (Pavo cristatus) show a strong mating preference for males with elaborate trains. This, however, poses something of a paradox because intense directional selection should erode genetic variation in the males’ trains, so that females will no longer benefit by discriminating among males on the basis of these traits. This situation is known as the ‘lek paradox’, and leads to the theoretical expectation of low heritability in the peacock’s train. We used two independent breeding experiments, involving a total of 42 sires and 86 of their male offspring, to estimate the narrow sense heritabilities of male ornaments and other morphometric traits. Contrary to expectation, we found significant levels of heritability in a trait known to be used by females during mate choice (train length), while no significant heritabilities were evident for other, non-fitness related morphological traits (tarsus length, body weight or spur length). This study adds to the building body of evidence that high levels of additive genetic variance can exist in secondary sexual traits under directional selection, but further emphasizes the main problem of what maintains this variation.     

Genetic variation in a female genital trait evolved by sexual coevolution

Understanding the patterns of genetic variation of traits subject to sexual selection is fundamental for explaining its evolutionary dynamics and potential for sexual coevolution. The signa of female Lepidoptera are sclerotized structures located on the inner surface of the genital receptacle that receives the spermatophore during copulation (the corpus bursae), whose main function is tearing the spermatophore envelope. Comparative data indicate that the evolution of signa has been influenced by sexually antagonistic coevolution with spermatophore envelopes. We looked for additive genetic variation in the size and shape of signa in females of the butterfly Callophrys xami (Lycaenidae) from two localities (BG and FC) in Mexico City. We also looked for genetic variation in female body size and in the size of corpus bursae. There were significant between-population differences in female body size, signa width and three signa shape traits. We found significant extranuclear maternal effects in one component of signa shape in the BG population, and in body weight, signa length and in one uniform component of signa shape in the FC population. Extranuclear maternal contributions could permit the evolution of female adaptations even if these reduce male fitness. We found additive genetic variation in signa length and width only in one population (BG); heritability estimates were high: 0.96 and 0.8, respectively. The existence of additive genetic variation in signa size could be, at least in part, a result of relaxed sexually antagonistic selection pressures due to the low level of polyandry exhibited by this species. Our results imply that there is currently potential for further sexual coevolution in this trait.     

Possibilities to breed for resistance to nematode parasite infections in small ruminants in tropical production systems

Gastrointestinal nematode parasitism is the most important disease affecting livestock production systems in developing countries, particularly small ruminant production systems. Of particular importance are infections with the strongyle Haemonchus contortus. Integrated disease control strategies are required, including improved management, nutrition and wise use of anthelmintic chemicals. Increasingly, selection of sheep or goats for improved nematode resistance is viewed as a valuable option to complement other control measures. Breeding for resistance is possible because of the existence of extensive genetic variation in resistance, both within and between breeds of sheep and goats. Such breeding schemes are most likely to be based on choice of appropriate breeds adapted to the local environmental conditions, followed by phenotypic selection for resistance. Goal and selection objective traits are likely to include performance (e.g. growth rate) under conditions of parasite challenge, faecal egg count (FEC) and measures of anaemia. With current technologies, genetic markers are likely to be too expensive and logistically difficult to incorporate into breeding schemes in tropical or developing countries. Genotype by environment interactions may be expected, particularly when comparing animals in environments that differ in the extent of parasite challenge or differ in the quality of available nutrition. However, there is no reason to expect antagonistic genetic relationships between performance and resistance, and selection indices should be readily constructed that improve both performance and resistance. If FEC is decreased, then pasture contamination should also decrease, leading to additional benefits for all sheep grazing the same pasture. Finally, breeding for nematode resistance should lead to lasting and sustained improvements in resistance or tolerance. There is no empirical evidence to suggest that nematodes will evolve rapidly in response to resistant hosts, and mathematical models based on genetic and biological principles also suggest that resistance should be sustainable.     

Environmental and genetic influences on queen and worker body size in the social wasp Vespula maculifrons

Many social insects exhibit morphologically distinct worker and queen castes that perform different functions. These functional differences may generate unique selection regimes operating on body size. For example, queens may be under directional selection for large body size, whereas directional selection on worker body size may be limited. Such contrasting selection pressures may differentially affect levels of genetic variation associated with size variation in the two castes. This study sought to determine if genetic effects underlying phenotypic differences varied between the worker and queen castes of the social wasp Vespula maculifrons. We predicted that directional selection would remove genetic variation associated with size differences in the queen caste, whereas a lack of directional selection would tend to maintain genetic variation associated with size differences in the worker caste. We thus (1) calculated broad and narrow sense heritabilities for several morphological traits, (2) examined whether some paternal genotypes produced more morphologically diverse offspring than others, and (3) determined whether trait size variation was associated with genetic variation within colonies. We found that few morphological traits were significantly heritable, indicating that little genetic variance for those traits existed within our study population. We also found that some patrilines produced more morphologically variable offspring than others, suggesting a role of genotype in phenotypic plasticity. And finally, no significant correlations between genetic diversity arising from multiple mating by queens within colonies and trait variation in either caste were found. Overall, our findings indicate a weak effect of genotype on both worker and queen body size variation and are suggestive of a large environmental influence on morphological trait size. Moreover, our results do not indicate that levels of genetic variation underlying size variation differ substantially between castes in this species.     

Live fast, die young: trade-offs between fitness components and sexually antagonistic selection on weaponry in Soay sheep

Males are predicted to compete for reproductive opportunities, with sexual selection driving the evolution of large body size and weaponry through the advantage they confer for access to females. Few studies have explored potential trade-offs of investment in secondary sexual traits between different components of fitness or tested for sexually antagonistic selection pressures. These factors may provide explanations for observed polymorphisms in both form and quality of secondary sexual traits. We report here an analysis of selection on horn phenotype in a feral population of Soay sheep (Ovis aries) on the island of Hirta, St. Kilda, Scotland. Soay sheep display a phenotypic polymorphism for horn type with males growing either normal or reduced (scurred) horns, and females growing either normal, scurred, or no (polled) horns; further variation in size exists within horn morphs. We show that horn phenotype and the size of the trait displayed is subject to different selection pressures in males and females, generating sexually antagonistic selection. Furthermore, there was evidence of a trade-off between breeding success and longevity in normal-horned males, with both the normal horn type and larger horn size being associated with greater annual breeding success but reduced longevity. Therefore, selection through lifetime breeding success was not found to act upon horn phenotype in males. In females, a negative association of annual breeding success within the normal-horned phenotype did not result in a significant difference in lifetime fitness when compared to scurred individuals, as no significant difference in longevity was found. However, increased horn size within this group was negatively associated with breeding success and longevity. Females without horns (polled) suffered reduced longevity and thus reduced lifetime breeding success relative the other horn morphs. Our results therefore suggest that trade-offs between different components of fitness and antagonistic selection between the sexes may maintain genetic variation for secondary sexual traits within a population.     

The Impact of Environmental Heterogeneity on Genetic Architecture in a Wild Population of Soay Sheep

This work demonstrates that environmental conditions experienced by individuals can shape their development and affect the stability of genetic associations. The implication of this observation is that the environmental response may influence the evolution of traits in the wild. Here, we examined how the genetic architecture of a suite of sexually dimorphic traits changed as a function of environmental conditions in an unmanaged population of Soay sheep (Ovis aries) on the island of Hirta, St. Kilda, northwest Scotland. We examined the stability of phenotypic, genetic, and environmental (residual) covariance in males during the first year of life between horn length, body weight, and parasite load in environments of different quality. We then examined the same covariance structures across environments within and between the adult sexes. We found significant genotype-by-environment interactions for lamb male body weight and parasite load, leading to a change in the genetic correlation among environments. Horn length was genetically correlated with body weight in males but not females and the genetic correlation among traits within and between the sexes was dependent upon the environmental conditions experienced during adulthood. Genetic correlations were smaller in more favorable environmental conditions, suggesting that in good environments, loci are expressed that have sex-specific effects. The reduction in genetic correlation between the sexes may allow independent evolutionary trajectories for each sex. This study demonstrates that the genetic architecture of traits is not stable under temporally varying environments and highlights the fact that evolutionary processes may depend largely upon ecological conditions.     

Metazoan-protozoan parasite co-infections and host body weight in St Kilda Soay sheep

For hundreds of years, the unmanaged Soay sheep population on St Kilda has survived despite enduring presumably deleterious co-infections of helminth, protozoan and arthropod parasites and intermittent periods of starvation. Important parasite taxa in young Soay sheep are strongyles (Trichostrongylus axei, Trichostrongylus vitrinus and Teladorsagia circumcincta), coccidia (11 Eimeria species) and keds (Melophagus ovinus) and in older animals, Teladorsagia circumcincta. In this research, associations between the intensity of different parasite taxa were investigated. Secondly, the intensities of different parasite taxa were tested for associations with variation in host weight, which is itself a determinant of over-winter survival in the host population. In lambs, the intensity of strongyle eggs was positively correlated with that of Nematodirus spp. eggs, while in yearlings and adults strongyle eggs and coccidia oocysts were positively correlated. In lambs and yearlings, of the parasite taxa tested, only strongyle eggs were significantly and negatively associated with host weight. However, in adult hosts, strongyles and coccidia were independently and negatively associated with host weight. These results are consistent with the idea that strongyles and coccidia are exerting independent selection on Soay sheep.     

Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population

Understanding the genetic architecture of phenotypic variation in natural populations is a fundamental goal of evolutionary genetics. Wild Soay sheep (Ovis aries) have an inherited polymorphism for horn morphology in both sexes, controlled by a single autosomal locus, Horns. The majority of males have large normal horns, but a small number have vestigial, deformed horns, known as scurs; females have either normal horns, scurs or no horns (polled). Given that scurred males and polled females have reduced fitness within each sex, it is counterintuitive that the polymorphism persists within the population. Therefore, identifying the genetic basis of horn type will provide a vital foundation for understanding why the different morphs are maintained in the face of natural selection. We conducted a genome-wide association study using ～36000 single nucleotide polymorphisms (SNPs) and determined the main candidate for Horns as RXFP2, an autosomal gene with a known involvement in determining primary sex characters in humans and mice. Evidence from additional SNPs in and around RXFP2 supports a new model of horn-type inheritance in Soay sheep, and for the first time, sheep with the same horn phenotype but different underlying genotypes can be identified. In addition, RXFP2 was shown to be an additive quantitative trait locus (QTL) for horn size in normal-horned males, accounting for up to 76% of additive genetic variation in this trait. This finding contrasts markedly from genome-wide association studies of quantitative traits in humans and some model species, where it is often observed that mapped loci only explain a modest proportion of the overall genetic variation.     

Natural selection on the genetical component of variance in body condition in a wild bird population

Although there is substantial evidence that skeletal measures of body size are heritable in wild animal populations, it is frequently assumed that the nonskeletal component of body weight (or ‘condition’) is determined primarily by environmental factors, in particular nutritional state. We tested this assumption by quantifying the genetic and environmental components of variance in fledgling body condition index (=relative body weight) in a natural population of collared flycatchers (Ficedula albicollis), and compared the strength of natural selection on individual breeding values with that on phenotypic values. A mixed model analysis of the components of variance, based on an ‘animal model’ and using 18 years of data on 17 717 nestlings, revealed a significant additive genetic component of variance in body condition, which corresponded to a narrow sense heritability (h²) of 0.30 (SE=0.03). Nongenetic contributions to variation in body condition were large, but there was no evidence of dominance variance nor of contributions from early maternal or common environment effects (pre‐manipulation environment) in condition at fledging. Comparison of pre‐ and post‐selection samples revealed virtually identical h² of body condition index, despite the fact that there was a significant decrease (35%) in the levels of additive genetic variance from fledging to breeding. The similar h² in the two samples occurred because the environmental component of variance was also reduced by selection, suggesting that natural selection was acting on both genotypic and environmental variation. The effects of selection on genetic variance were confirmed by calculation of the selection differentials for both phenotypic values and best linear unbiased predictor (BLUP) estimates of breeding values: there was positive directional selection on condition index both at the phenotypic and the genotypic level. The significant h² of body condition index is consistent with data from human and rodent populations showing significant additive genetic variance in relative body mass and adiposity, but contrasts with the common assumption in ecology that body condition reflects an individual’s nongenetic nutritional state. Furthermore, the substantial reduction in the additive genetic component of variance in body condition index suggests that selection on environmental deviations cannot alone explain the maintenance of additive genetic variation in heritable traits, but that other mechanisms are needed to explain the moderate to high heritabilities of traits under consistent and strong directional selection.     

SELECTION ON MOTHERS AND OFFSPRING: WHOSE PHENOTYPE IS IT AND DOES IT MATTER?

Reproductive and early life-history traits can be considered aspects of either offspring or maternal phenotype, and their evolution will therefore depend on selection operating through offspring and maternal components of fitness. Furthermore, selection at these levels may be antagonistic, with optimal offspring and maternal fitness occurring at different phenotypic values. We examined selection regimes on the correlated traits of birth weight, birth date, and litter size in Soay sheep (Ovis aries) using data from a long-term study of a free-living population on the archipelago of St. Kilda, Scotland. We tested the hypothesis that selective constraints on the evolution of the multivariate phenotype arise through antagonistic selection, either acting at offspring and maternal levels, or on correlated aspects of phenotype. All three traits were found to be under selection through variance in short-term and lifetime measures of fitness. Analysis of lifetime fitness revealed strong positive directional selection on birth weight and weaker selection for increased birth date at both levels. However, there was also evidence for stabilizing selection on these traits at the maternal level, with reduced fitness at high phenotypic values indicating lower phenotypic optima for mothers than for offspring. Additionally, antagonistic selection was found on litter size. From the offspring's point of view it is better to be born a singleton, whereas maternal fitness increases with average litter size. The decreased fitness of twins is caused by their reduced birth weight; therefore, this antagonistic selection likely results from trade-offs between litter size and birth weight that have different optimal resolutions with respect to offspring and maternal fitness. Our results highlight how selection regimes may vary depending on the assignment of reproductive and early life-history traits to either offspring or maternal phenotype.