Sexual selection, through intra-male competition or female choice, is assumed to be a source of strong and sustained directional selection in the wild. In the presence of such strong directional selection, alleles enhancing a particular trait are predicted to become fixed within a population, leading to a decrease in the underlying genetic variation. However, there is often considerable genetic variation underlying sexually selected traits in wild populations, and consequently, this phenomenon has become a long-discussed issue in the field of evolutionary biology. In wild Soay sheep, large horns confer an advantage in strong intra-sexual competition, yet males show an inherited polymorphism for horn type and have substantial genetic variation in their horn size. Here we show that most genetic variation in this trait is maintained by a trade-off between natural and sexual selection at a single gene, relaxin-like receptor 2 (RXFP2). We found that an allele conferring larger horns, Ho(+), is associated with higher reproductive success, whereas a smaller horn allele, Ho(P), confers increased survival, resulting in a net effect of overdominance (that is, heterozygote advantage) for fitness at RXFP2. The nature of this trade-off is simple relative to commonly proposed explanations for the maintenance of sexually selected traits, such as genic capture ('good genes') and sexually antagonistic selection. Our results demonstrate that by identifying the genetic architecture of trait variation, we can determine the principal mechanisms maintaining genetic variation in traits under strong selection and explain apparently counter-evolutionary observations.

The estimation of quantitative genetic parameters in wild populations is generally limited by the accuracy and completeness of the available pedigree information. Using relatedness at genomewide markers can potentially remove this limitation and lead to less biased and more precise estimates. We estimated heritability, maternal genetic effects and genetic correlations for body size traits in an unmanaged long-term study population of Soay sheep on St Kilda using three increasingly complete and accurate estimates of relatedness: (i) Pedigree 1, using observation-derived maternal links and microsatellite-derived paternal links; (ii) Pedigree 2, using SNP-derived assignment of both maternity and paternity; and (iii) whole-genome relatedness at 37 037 autosomal SNPs. In initial analyses, heritability estimates were strikingly similar for all three methods, while standard errors were systematically lower in analyses based on Pedigree 2 and genomic relatedness. Genetic correlations were generally strong, differed little between the three estimates of relatedness and the standard errors declined only very slightly with improved relatedness information. When partitioning maternal effects into separate genetic and environmental components, maternal genetic effects found in juvenile traits increased substantially across the three relatedness estimates. Heritability declined compared to parallel models where only a maternal environment effect was fitted, suggesting that maternal genetic effects are confounded with direct genetic effects and that more accurate estimates of relatedness were better able to separate maternal genetic effects from direct genetic effects. We found that the heritability captured by SNP markers asymptoted at about half the SNPs available, suggesting that denser marker panels are not necessarily required for precise and unbiased heritability estimates. Finally, we present guidelines for the use of genomic relatedness in future quantitative genetics studies in natural populations.

Meiotic recombination breaks down linkage disequilibrium (LD) and forms new haplotypes, meaning that it is an important driver of diversity in eukaryotic genomes. Understanding the causes of variation in recombination rate is important in interpreting and predicting evolutionary phenomena and in understanding the potential of a population to respond to selection. However, despite attention in model systems, there remains little data on how recombination rate varies at the individual level in natural populations. Here we used extensive pedigree and high-density SNP information in a wild population of Soay sheep (Ovis aries) to investigate the genetic architecture of individual autosomal recombination rates. Individual rates were high relative to other mammal systems and were higher in males than in females (autosomal map lengths of 3748 and 2860 cM, respectively). The heritability of autosomal recombination rate was low but significant in both sexes (h2 = 0.16 and 0.12 in females and males, respectively). In females, 46.7% of the heritable variation was explained by a subtelomeric region on chromosome 6; a genome-wide association study showed the strongest associations at locus RNF212, with further associations observed at a nearby ∼374-kb region of complete LD containing three additional candidate loci, CPLX1, GAK, and PCGF3. A second region on chromosome 7 containing REC8 and RNF212B explained 26.2% of the heritable variation in recombination rate in both sexes. Comparative analyses with 40 other sheep breeds showed that haplotypes associated with recombination rates are both old and globally distributed. Both regions have been implicated in rate variation in mice, cattle, and humans, suggesting a common genetic architecture of recombination rate variation in mammals.

Tetraploid inheritance has two extremes: disomic in allotetraploids and tetrasomic in autotetraploids. The possibility of mixed, or intermediate, inheritance models has generally been neglected. These could well apply to newly formed hybrids or to diploidizing (auto)tetraploids. We present a simple likelihood-based approach that is able to incorporate disomic, tetrasomic, and intermediate inheritance models and estimates the double-reduction rate. Our model shows that inheritance of microsatellite markers in natural tetraploids of Rorippa amphibia and R. sylvestris is tetrasomic, confirming their autotetraploid origin. However, in F 1 hybrids inheritance was intermediate to disomic and tetrasomic inheritance. Apparently, in meiosis, chromosomes paired preferentially with the homolog from the same parental species, but not strictly so. Detected double-reduction rates were low. We tested the general applicability of our model, using published segregation data. In two cases, an intermediate inheritance model gave a better fit to the data than the tetrasomic model advocated by the authors. The existence of inheritance intermediate to disomic and tetrasomic has important implications for linkage mapping and population genetics and hence breeding programs of tetraploids. Methods that have been developed for either disomic or tetrasomic tetraploids may not be generally applicable, particularly in systems where hybridization is common.

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.

The degree to which changes in lifespan are coupled to changes in senescence in different physiological systems and phenotypic traits is a central question in biogerontology. It is underpinned by deeper biological questions about whether or not senescence is a synchronised process, or whether levels of synchrony depend on species or environmental context. Understanding how natural selection shapes patterns of synchrony in senescence across physiological systems and phenotypic traits demands the longitudinal study of many phenotypes under natural conditions. Here, we examine the patterns of age-related variation in late adulthood in a wild population of Soay sheep (Ovis aries) that have been the subject of individual-based monitoring for thirty years. We examined twenty different phenotypic traits in both males and females, encompassing vital rates (survival and fecundity), maternal reproductive performance (offspring birth weight, birth date and survival), male rutting behaviour, home range measures, parasite burdens, and body mass. We initially quantified age-related variation in each trait having controlled for annual variation in the environment, among-individual variation and selective disappearance effects. We then standardised our age-specific trait means and tested whether age trajectories could be meaningfully grouped according to sex or the type of trait. Whilst most traits showed age-related declines in later life, we found striking levels of asynchrony both within and between the sexes. Of particular note, female fecundity and reproductive performance declined with age, but male annual reproductive success did not. We also discovered that whilst home range size and quality decline with age in females, home range size increases with age in males. Our findings highlight the complexity of phenotypic ageing under natural conditions and, along with emerging data from other wild populations and laboratory models, suggest that the long-standing hypothesis within evolutionary biology that fitness-related traits should senesce in a synchronous manner is seriously flawed.

A 25-year study of wild sheep shows that individuals vary in how quickly they lose weight as parasite infections increase, and that those who lose the least weight when heavily infected produce more offspring.

There is ample evidence for inbreeding depression manifested as a reduction in fitness or fitness‐related traits in the focal individual. In many organisms, fitness is not only affected by genes carried by the individual, but also by genes carried by their parents, for example if receiving parental care. While maternal effects have been described in many systems, the extent to which inbreeding affects fitness directly through the focal individual, or indirectly through the inbreeding coefficients of its parents, has rarely been examined jointly. The Soay sheep study population is an excellent system in which to test for both effects, as lambs receive extended maternal care. Here, we tested for both maternal and individual inbreeding depression in three fitness‐related traits (birthweight and weight and hindleg length at 4 months of age) and three fitness components (first‐year survival, adult annual survival and annual breeding success), using either pedigree‐derived inbreeding or genomic estimators calculated using ~37 000 SNP markers. We found evidence for inbreeding depression in 4‐month hindleg and weight, first‐year survival in males, and annual survival and breeding success in adults. Maternal inbreeding was found to depress both birthweight and 4‐month weight. We detected more instances of significant inbreeding depression using genomic estimators than the pedigree, which is partly explained through the increased sample sizes available. In conclusion, our results highlight that cross‐generational inbreeding effects warrant further exploration in species with parental care and that modern genomic tools can be used successfully instead of, or alongside, pedigrees in natural populations.