pubmed.ncbi.nlm.nih.gov

Genomic analysis of morphometric traits in bighorn sheep using the Ovine Infinium® HD SNP BeadChip - PubMed

  • ️Mon Jan 01 2018

Genomic analysis of morphometric traits in bighorn sheep using the Ovine Infinium® HD SNP BeadChip

Joshua M Miller et al. PeerJ. 2018.

Abstract

Elucidating the genetic basis of fitness-related traits is a major goal of molecular ecology. Traits subject to sexual selection are particularly interesting, as non-random mate choice should deplete genetic variation and thereby their evolutionary benefits. We examined the genetic basis of three sexually selected morphometric traits in bighorn sheep (Ovis canadensis): horn length, horn base circumference, and body mass. These traits are of specific concern in bighorn sheep as artificial selection through trophy hunting opposes sexual selection. Specifically, horn size determines trophy status and, in most North American jurisdictions, if an individual can be legally harvested. Using between 7,994-9,552 phenotypic measures from the long-term individual-based study at Ram Mountain (Alberta, Canada), we first showed that all three traits are heritable (h2 = 0.15-0.23). We then conducted a genome-wide association study (GWAS) utilizing a set of 3,777 SNPs typed in 76 individuals using the Ovine Infinium® HD SNP BeadChip. We found suggestive association for body mass at a single locus (OAR9_91647990). The absence of strong associations with SNPs suggests that the traits are likely polygenic. These results represent a step forward for characterizing the genetic architecture of fitness related traits in sexually dimorphic ungulates.

Keywords: Animal model; GWAS; Ovis canadensis; SNP.

PubMed Disclaimer

Conflict of interest statement

David W Coltman is an Academic Editor for PeerJ.

Figures

Figure 1
Figure 1. Manhattan plots for morphological characteristics.

Horn length (A), horn base circumference (C), and body mass (E). The blue line represents the genome-wide significance threshold; the red line represents the threshold for suggestive association. Positions are relative to the domestic sheep genome assembly (version 3.1; Jiang et al., 2014). The green arrow indicates the suggestive locus for body mass. Next to each Manhattan plot is the corresponding QQ-plot (B, D, and F), with the genomic inflation factor (λ) and standard error indicated in the bottom right of each plot. The black line shows a 1:1 correspondence while the red line is a regression through the observed data.

Figure 2
Figure 2. Scatterplot of LD estimates versus inter-markers distance.

A non-linear least squares regression line is shown, with the round point indicating the half-length estimate.

Figure 3
Figure 3

Heat maps of expected percent power of a GWAS as a function of sample size and effect size for linkage disequilibrium (LD) estimates of 0.75 (A), 0.50 (B), and 0.25 (C). Light colors indicate higher power to detect associations (D). Dotted red lines correspond to the number of samples used in this study (N = 76).

Similar articles

Cited by

References

    1. Allais-Bonnet A, Grohs C, Medugorac I, Krebs S, Djari A, Graf A, Fritz S, Seichter D, Baur A, Russ I, Bouet S, Rothammer S, Wahlberg P, Esquerre D, Hoze C, Boussaha M, Weiss B, Thepot D, Fouilloux MN, Rossignol MN, Van Marle-Koster E, Hreidarsdottir GE, Barbey S, Dozias D, Cobo E, Reverse P, Catros O, Marchand JL, Soulas P, Roy P, Marquant-Leguienne B, Le Bourhis D, Clement L, Salas-Cortes L, Venot E, Pannetier M, Phocas F, Klopp C, Rocha D, Fouchet M, Journaux L, Bernard-Capel C, Ponsart C, Eggen A, Blum H, Gallard Y, Boichard D, Pailhoux E, Capitan A. Novel insights into the bovine polled phenotype and horn ontogenesis in bovidae. PLOS ONE. 2013;8(5):e63512. doi: 10.1371/journal.pone.0063512. - DOI - PMC - PubMed
    1. Andrews KR, Good JM, Miller MR, Luikart G, Hohenlohe PA. Harnessing the power of RADseq for ecological and evolutionary genomics. Nature Reviews Genetics. 2016;17:81–92. doi: 10.1038/nrg.2015.28. - DOI - PMC - PubMed
    1. Aulchenko YS, Ripke S, Isaacs A, Van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics. 2007;23:1294–1296. doi: 10.1093/bioinformatics/btm108. - DOI - PubMed
    1. Backström N, Qvarnström A, Gustafsson L, Ellegren H. Levels of linkage disequilibrium in a wild bird population. Biology Letters. 2006;2:435–438. doi: 10.1098/rsbl.2006.0507. - DOI - PMC - PubMed
    1. Barnett IJ, Lee S, Lin XH. Detecting rare variant effects using extreme phenotype sampling in sequencing association studies. Genetic Epidemiology. 2013;37:142–151. doi: 10.1002/gepi.21699. - DOI - PMC - PubMed

Grants and funding

Field work at RM has been supported by National Science and Engineering Research Council (NSERC) Discovery Grants, Alberta Conservation Association Grants in Biodiversity to Marco Festa-Bianchet. Alberta Fish & Wildlife provide logistic and financial support. The molecular work was supported by an NSERC Discovery Grant to David Coltman, as well as an Alberta Conservation Association Grant in Biodiversity, and an Alberta Sport, Recreation, Parks, and Wildlife Foundation Development Initiatives Program grant to Joshua Miller. Joshua Miller’s graduate research was supported by an NSERC Vanier scholarship, the Killam Foundation, and Alberta Innovates Technology Futures. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.