Genomic prediction in the wild: A case study in Soay sheep.

Genomic prediction, the technique whereby an individual's genetic component of their phenotype is estimated from its genome, has revolutionised animal and plant breeding and medical genetics. However, despite being first introduced nearly two decades ago, it has hardly been adopted by the evolutionary genetics community studying wild organisms. Here, genomic prediction is performed on eight traits in a wild population of Soay sheep. The population has been the focus of a >30 year evolutionary ecology study and there is already considerable understanding of the genetic architecture of the focal Mendelian and quantitative traits. We show that the accuracy of genomic prediction is high for all traits, but especially those with loci of large effect segregating. Five different methods are compared, and the two methods that can accommodate zero-effect and large-effect loci in the same model tend to perform best. If the accuracy of genomic prediction is similar in other wild populations, then there is a real opportunity for pedigree-free molecular quantitative genetics research to be enabled in many more wild populations; currently the literature is dominated by studies that have required decades of field data collection to generate sufficiently deep pedigrees. Finally, some of the potential applications of genomic prediction in wild populations are discussed.

A High-Density Linkage Map Reveals Sexual Dimorphism in Recombination Landscapes in Red Deer (Cervus elaphus).

High-density linkage maps are an important tool to gain insight into the genetic architecture of traits of evolutionary and economic interest, and provide a resource to characterize variation in recombination landscapes. Here, we used information from the cattle genome and the 50 K Cervine Illumina BeadChip to inform and refine a high-density linkage map in a wild population of red deer (Cervus elaphus). We constructed a predicted linkage map of 38,038 SNPs and a skeleton map of 10,835 SNPs across 34 linkage groups. We identified several chromosomal rearrangements in the deer lineage relative to sheep and cattle, including six chromosome fissions, one fusion, and two large inversions. Otherwise, our findings showed strong concordance with map orders in the cattle genome. The sex-averaged linkage map length was 2739.7 cM and the genome-wide autosomal recombination rate was 1.04 cM/Mb. The female autosomal map length was 1.21 longer than that of males (2767.4 cM vs. 2280.8 cM, respectively). Sex differences in map length were driven by high female recombination rates in peri-centromeric regions, a pattern that is unusual relative to other mammal species. This effect was more pronounced in fission chromosomes that would have had to produce new centromeres. We propose two hypotheses to explain this effect: (1) that this mechanism may have evolved to counteract centromeric drive associated with meiotic asymmetry in oocyte production; and/or (2) that sequence and structural characteristics suppressing recombination in close proximity to the centromere may not have evolved at neo-centromeres. Our study provides insight into how recombination landscapes vary and evolve in mammals, and will provide a valuable resource for studies of evolution, genetic improvement, and population management in red deer and related species.

Genomic analysis reveals depression due to both individual and maternal inbreeding in a free-living mammal population.

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.

Inbreeding depression across the lifespan in a wild mammal population.

Inbreeding depression is of major concern for the conservation of threatened species, and inbreeding avoidance is thought to be a key driver in the evolution of mating systems. However, the estimation of individual inbreeding coefficients in natural populations has been challenging, and, consequently, the full effect of inbreeding on fitness remains unclear. Genomic inbreeding coefficients may resolve the long-standing paucity of data on inbreeding depression in adult traits and total fitness. Here we investigate inbreeding depression in a range of life history traits and fitness in a wild population of red deer (Cervus elaphus) in Scotland using individual inbreeding coefficients derived from dense Single-Nucleotide Polymorphism (SNP) data (Fgrm). We find associations between[Formula: see text]and annual breeding success in both sexes, and between maternal inbreeding coefficient and offspring survival. We also confirm previous findings of inbreeding depression in birth weight and juvenile survival. In contrast, inbreeding coefficients calculated from a deep and comparatively complete pedigree detected inbreeding depression in juvenile survival, but not in any adult fitness component. The total effect of inbreeding on lifetime breeding success (LBS) was substantial in both sexes: for Fgrm = 0.125, a value resulting from a half-sib mating, LBS declined by 72% for females and 95% for males. Our results demonstrate that SNP-based estimates of inbreeding provide a powerful tool for evaluating inbreeding depression in natural populations, and suggest that, to date, the prevalence of inbreeding depression in adult traits may have been underestimated.

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.

Estimating quantitative genetic parameters in wild populations: a comparison of pedigree and genomic approaches.

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.