Genomic analysis reveals a polygenic architecture of antler morphology in wild red deer (Cervus elaphus).

Sexually selected traits show large variation and rapid evolution across the animal kingdom, yet genetic variation often persists within populations despite apparent directional selection. A key step in solving this long-standing paradox is to determine the genetic architecture of sexually selected traits to understand evolutionary drivers and constraints at the genomic level. Antlers are a form of sexual weaponry in male red deer (Cervus elaphus). On the island of Rum, Scotland, males with larger antlers have increased breeding success, yet there has been no evidence of any response to selection at the genetic level. To try and understand the mechanisms underlying this observation, we investigate the genetic architecture of ten antler traits and their principal components using genomic data from >38,000 SNPs. We estimate the heritabilities and genetic correlations of the antler traits using a genomic relatedness approach. We then use genome-wide association and haplotype-based regional heritability to identify regions of the genome underlying antler morphology, and an empirical Bayes approach to estimate the underlying distributions of allele effect sizes. We show that antler morphology is highly repeatable over an individual's lifetime, heritable and has a polygenic architecture and that almost all antler traits are positively genetically correlated with some loci identified as having pleiotropic effects. Our findings suggest that a large mutational target and genetic covariances among antler traits, in part maintained by pleiotropy, are likely to contribute to the maintenance of genetic variation in antler morphology in this population.

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.

Estimating selection on the act of inbreeding in a population with strong inbreeding depression.

Inbreeding depression is widely regarded as a driving force in the evolution of dispersal, mate choice and sperm selection. However, due to likely costs of inbreeding avoidance, which are poorly understood, it is unclear to what extent selection to avoid inbreeding is expected in nature. Moreover, there are currently very few empirical estimates of the strength of selection against the act of inbreeding (mating with a relative), as opposed to the fitness costs of being inbred. Here, we use data from the individual-based study of red deer on the Scottish island of Rum, a strongly polygynous system which harbours a large inbreeding load, to estimate selection against the act of inbreeding for each sex. We use pedigree and genomic estimates of relatedness between individuals and measure fitness using both lifetime breeding success (number of calves born) and lifetime reproductive success (number of calves surviving to independence), with the latter incorporating inbreeding depression in calf survival. We find for both sexes that the repeatability of the act of inbreeding was low (< 0.1), suggesting little among-individual variation for this trait on which selection can act. Using the genomic measures, there was significant selection against the act of inbreeding in males, but not in females, and there was considerable uncertainty in the estimate in both sexes. We discuss possible explanations for these patterns and their implications for understanding the evolution of inbreeding avoidance in natural populations.

Modelling wild-domestic interbreeding: how selection on a quantitative trait affects gene flow at a neutral locus.

Interbreeding between domesticated individuals and their wild conspecifics occurs in a range of species. The rate of gene flow into the wild population, estimated using genetic markers, is often smaller than the fraction of immigrants, as immigrants and their descendants generally have lower relative fitness. Here the difference between one-way migration rate and gene flow (effective migration rate) is explored using quantitative genetic simulations. We model a trait undergoing stabilizing selection in the recipient population, influenced by an infinite number of loci, each with small effect. Immigrants have suboptimal trait values, and differ in allele frequency at an unlinked neutral marker locus. We derive an analytical approximation for the effective migration rate, and show that in the limiting case of low migration rates, the ratio between effective and actual migration rate approximately equals the ratio between mean fitness of immigrants and in the admixed population. This ignores indirect selection on the marker locus within the admixed population due to covariance with the trait value, which may be substantial when the genetic difference between the populations is large. For higher migration rates, the gene flow/migration ratio rises with increasing migration rate, inflating the rate of gene flow.

Increased genetic marker density reveals high levels of admixture between red deer and introduced Japanese sika in Kintyre, Scotland.

Hybridization is a natural process at species range boundaries, but increasing numbers of species are hybridizing due to direct or indirect human activities. In such cases of anthropogenic hybridization, subsequent introgression can threaten the survival of native species. To date, many such systems have been studied with too few genetic markers to assess the level of threat resulting from advanced backcrossing. Here, we use 44,999 single nucleotide polymorphisms (SNPs) and the ADMIXTURE program to study two areas of Scotland where a panel of 22 diagnostic microsatellites previously identified introgression between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon). In Kintyre, we reclassify 26% of deer from the pure species categories to the hybrid category whereas in the NW Highlands we only reclassify 2%. As expected, the reclassified individuals are mostly advanced backcrosses. We also investigate the ability of marker panels selected on different posterior allele frequency criteria to find hybrids assigned by the full marker set and show that in our data, ancestry informative markers (i.e. those that are highly differentiated between the species, but not fixed) are better than diagnostic markers (those markers that are fixed between the species) because they are more evenly distributed in the genome. Diagnostic loci are concentrated on the X chromosome to the detriment of autosomal coverage.

A Genomic Region Containing REC8 and RNF212B Is Associated with Individual Recombination Rate Variation in a Wild Population of Red Deer (Cervus elaphus).

Recombination is a fundamental feature of sexual reproduction, ensuring proper disjunction, preventing mutation accumulation and generating new allelic combinations upon which selection can act. However it is also mutagenic, and breaks up favorable allelic combinations previously built up by selection. Identifying the genetic drivers of recombination rate variation is a key step in understanding the causes and consequences of this variation, how loci associated with recombination are evolving and how they affect the potential of a population to respond to selection. However, to date, few studies have examined the genetic architecture of recombination rate variation in natural populations. Here, we use pedigree data from ∼ 2,600 individuals genotyped at ∼ 38,000 SNPs to investigate the genetic architecture of individual autosomal recombination rate in a wild population of red deer (Cervus elaphus). Female red deer exhibited a higher mean and phenotypic variance in autosomal crossover counts (ACC). Animal models fitting genomic relatedness matrices showed that ACC was heritable in females ([Formula: see text] = 0.12) but not in males. A regional heritability mapping approach showed that almost all heritable variation in female ACC was explained by a genomic region on deer linkage group 12 containing the candidate loci REC8 and RNF212B, with an additional region on linkage group 32 containing TOP2B approaching genome-wide significance. The REC8/RNF212B region and its paralogue RNF212 have been associated with recombination in cattle, mice, humans and sheep. Our findings suggest that mammalian recombination rates have a relatively conserved genetic architecture in both domesticated and wild systems, and provide a foundation for understanding the association between recombination loci and individual fitness within this population.

Pedigree reconstruction from SNP data: parentage assignment, sibship clustering and beyond.

Data on hundreds or thousands of single nucleotide polymorphisms (SNPs) provide detailed information about the relationships between individuals, but currently few tools can turn this information into a multigenerational pedigree. I present the r package sequoia, which assigns parents, clusters half-siblings sharing an unsampled parent and assigns grandparents to half-sibships. Assignments are made after consideration of the likelihoods of all possible first-, second- and third-degree relationships between the focal individuals, as well as the traditional alternative of being unrelated. This careful exploration of the local likelihood surface is implemented in a fast, heuristic hill-climbing algorithm. Distinction between the various categories of second-degree relatives is possible when likelihoods are calculated conditional on at least one parent of each focal individual. Performance was tested on simulated data sets with realistic genotyping error rate and missingness, based on three different large pedigrees (N = 1000-2000). This included a complex pedigree with overlapping generations, occasional close inbreeding and some unknown birth years. Parentage assignment was highly accurate down to about 100 independent SNPs (error rate <0.1%) and fast (<1 min) as most pairs can be excluded from being parent-offspring based on opposite homozygosity. For full pedigree reconstruction, 40% of parents were assumed nongenotyped. Reconstruction resulted in low error rates (<0.3%), high assignment rates (>99%) in limited computation time (typically <1 h) when at least 200 independent SNPs were used. In three empirical data sets, relatedness estimated from the inferred pedigree was strongly correlated to genomic relatedness.

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.

Comparison of non-Gaussian quantitative genetic models for migration and stabilizing selection.

The balance between stabilizing selection and migration of maladapted individuals has formerly been modeled using a variety of quantitative genetic models of increasing complexity, including models based on a constant expressed genetic variance and models based on normality. The infinitesimal model can accommodate nonnormality and a nonconstant genetic variance as a result of linkage disequilibrium. It can be seen as a parsimonious one-parameter model that approximates the underlying genetic details well when a large number of loci are involved. Here, the performance of this model is compared to several more realistic explicit multilocus models, with either two, several or a large number of alleles per locus with unequal effect sizes. Predictions for the deviation of the population mean from the optimum are highly similar across the different models, so that the non-Gaussian infinitesimal model forms a good approximation. It does, however, generally estimate a higher genetic variance than the multilocus models, with the difference decreasing with an increasing number of loci. The difference between multilocus models depends more strongly on the effective number of loci, accounting for relative contributions of loci to the variance, than on the number of alleles per locus.

Varying disease-mediated selection at different life-history stages of Atlantic salmon in fresh water.

Laboratory studies on associations between disease resistance and susceptibility and major histocompatibility (MH) genes in Atlantic salmon Salmo salar have shown the importance of immunogenetics in understanding the capacity of populations to fight specific diseases. However, the occurrence and virulence of pathogens may vary spatially and temporally in the wild, making it more complicated to predict the overall effect that MH genes exert on fitness of natural populations and over several life-history stages. Here we show that MH variability is a significant determinant of salmon survival in fresh water, by comparing observed and expected genotype frequencies at MH and control microsatellite loci at parr and migrant stages in the wild. We found that additive allelic effects at immunogenetic loci were more likely to determine survival than dominance deviation, and that selection on certain MH alleles varied with life stage, possibly owing to varying pathogen prevalence and/or virulence over time. Our results highlight the importance of preserving genetic diversity (particularly at MH loci) in wild populations, so that they have the best chance of adapting to new and increased disease challenges as a result of projected climate warming and increasing aquaculture.