Genetic architecture and polygenic risk score prediction of degenerative suspensory ligament desmitis (DSLD) in the Peruvian Horse.
- Journal Article
Summary
The research aims to understand the genetic factors that affect horses’ risk for developing degenerative suspensory ligament desmitis (DSLD), a degenerative tendon disease. Using a variety of genetic analysis methods on a pool of 183 Peruvian Horses, the team identified associated genomic regions, estimated the trait’s heritability, and created a polygenic risk score model to predict DSLD risk.
Overview of the Study
- The researchers gathered a reference population of 183 Peruvian Horses, identified as either having DSLD (the cases) or not (the controls). The Peruvian Horse was chosen due to the breed’s increased predisposition to DSLD.
- Using this population, they performed different forms of genetic analysis: a genome-wide association study (GWAS), a regional window variance analysis using local genome partitioning, a signatures of selection (SOS) analysis, and computed polygenic risk scores (PRS) for predicting DSLD risk.
- Besides this, the researchers estimated the trait heritability from pedigrees, using a larger population pool of 1,927 Peruvian horses.
Findings of the Study
- The trait’s heritability was calculated to be 0.22 ± 0.08. Heritability, in this context, refers to the extent that genetic variation can explain the observed variation in DSLD occurrence among the Peruvian Horse population.
- After setting a threshold for genome-wide significance, GWAS identified 151 risk single nucleotide polymorphisms (SNPs, variations in a single base pair in a DNA sequence) associated with DSLD.
- There were multiple regions across the genome showing enriched local heritability, particularly on chromosomes 1, 2, 6, 10, 13, 16, 18, 22, and the X chromosome. “Enriched local heritability” implies these chromosomal regions contribute more to the heritability of DSLD.
- The SOS analysis revealed 66 genes that had selection signatures in DSLD cases, which did not appear in the control group. A selection signature indicates these genes were subjected to selection, suggesting they contribute to disease susceptibility.
- Pathways enriched in DSLD cases were related to proteoglycan metabolism, extracellular matrix homeostasis, and certain signal transduction pathways, including the “hedgehog” signaling pathway. These are biological process and actions that, when disrupted, may negatively impact tendon health.
Polygenic Risk Score
- The team developed a predictive model using PRS, fitting the top 1% of SNPs through a Bayesian Ridge Regression model.
- This model showed strong predictive performance for DSLD risk, evidenced by a high mean R-value on both probit and logit liability scales. The higher the R-value, the more accurately the model predicts the outcome (DSLD risk).
In conclusion, the study reveals that DSLD in Peruvian Horses has a moderate heritability and polygenic architecture. The researchers were able to pinpoint multiple SNP associations in novel tendinopathy candidate genes influencing disease risk, opening avenues for further genetic study and potential preventative measures.
Cite This Article
Publication
Researcher Affiliations
- Department of Surgical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, WI, United States.
- Department of Surgical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, WI, United States.
- Department of Surgical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, WI, United States.
- Department of Surgical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, WI, United States.
- Department of Surgical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, WI, United States.
- Department of Veterinary Integrative Biosciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX, United States.
- Department of Veterinary Integrative Biosciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX, United States.
- Department of Animal and Dairy Sciences, College of Agriculture and Life Sciences, University of Wisconsin-Madison, Madison, WI, United States.
- Department of Surgical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, WI, United States.
- Department of Surgical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, WI, United States.
Conflict of Interest Statement
References
- Banfield RE, Hall LO, Bowyer KW, Kegelmeyer WP. A comparison of decision tree ensemble creation techniques.. IEEE Trans Pattern Anal Mach Intell 2007 Jan;29(1):173-80.
- Brisset M, Metay C, Carlier RY, Badosa C, Marques C, Schalkwijk J, vanVlijmen-Willems I, Jimenez-Mallebrera C, Keren B, Jobic V, Laforêt P, Malfatti E. Biallelic mutations in Tenascin-X cause classical-like Ehlers-Danlos syndrome with slowly progressive muscular weakness.. Neuromuscul Disord 2020 Oct;30(10):833-838.
- Broder C, Arnold P, Vadon-Le Goff S, Konerding MA, Bahr K, Müller S, Overall CM, Bond JS, Koudelka T, Tholey A, Hulmes DJ, Moali C, Becker-Pauly C. Metalloproteases meprin α and meprin β are C- and N-procollagen proteinases important for collagen assembly and tensile strength.. Proc Natl Acad Sci U S A 2013 Aug 27;110(35):14219-24.
- Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering.. Am J Hum Genet 2007 Nov;81(5):1084-97.
- Bryant FB, Yarnold PR. Principal-components analysis and exploratory and confirmatory factor analysis. in Reading and understanding multivariate statistics, 99–136.
- Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets.. Gigascience 2015;4:7.
- Chen T, Guestrin C. Xgboost: A scalable tree boosting system. Proc. 22nd Acm Sigkdd Int. Conf. Knowl. Discov. Data Min. 785–794.
- Cutter AD, Payseur BA. Genomic signatures of selection at linked sites: unifying the disparity among species.. Nat Rev Genet 2013 Apr;14(4):262-74.
- Dai G, Li Y, Liu J, Zhang C, Chen M, Lu P, Rui Y. Higher BMP Expression in Tendon Stem/Progenitor Cells Contributes to the Increased Heterotopic Ossification in Achilles Tendon With Aging.. Front Cell Dev Biol 2020;8:570605.
- Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, Handsaker RE, Lunter G, Marth GT, Sherry ST, McVean G, Durbin R. The variant call format and VCFtools.. Bioinformatics 2011 Aug 1;27(15):2156-8.
- Dempster ER, Lerner IM. Heritability of Threshold Characters.. Genetics 1950 Mar;35(2):212-36.
- Fenwick SA, Curry V, Harrall RL, Hazleman BL, Hackney R, Riley GP. Expression of transforming growth factor-beta isoforms and their receptors in chronic tendinosis.. J Anat 2001 Sep;199(Pt 3):231-40.
- Freedman BR, Adu-Berchie K, Barnum C, Fryhofer GW, Salka NS, Shetye S, Soslowsky LJ. Nonsurgical treatment reduces tendon inflammation and elevates tendon markers in early healing.. J Orthop Res 2022 Oct;40(10):2308-2319.
- Friedman J, Hastie T, Tibshirani R. Regularization Paths for Generalized Linear Models via Coordinate Descent.. J Stat Softw 2010;33(1):1-22.
- Friedman JH. Stochastic gradient boosting. Comput. Statistics Data Analysis 38 (4), 367–378.
- Giral H, Landmesser U, Kratzer A. Into the Wild: GWAS Exploration of Non-coding RNAs.. Front Cardiovasc Med 2018;5:181.
- de Simoni Gouveia JJ, da Silva MV, Paiva SR, de Oliveira SM. Identification of selection signatures in livestock species.. Genet Mol Biol 2014 Jun;37(2):330-42.
- Greenwell B, Boehmke B, Cunningham J, Developers GBM, Greenwell MB. Package ‘gbm. R. Package Version 2 (5).
- Hadfield JD. MCMC methods for multi-response generalized linear mixed models: the MCMCglmm R package. J. Stat. Softw. 33 (2), 1–22.
- Halper J, Kim B, Khan A, Yoon JH, Mueller PO. Degenerative suspensory ligament desmitis as a systemic disorder characterized by proteoglycan accumulation.. BMC Vet Res 2006 Apr 12;2:12.
- Haythorn A, Young M, Stanton J, Zhang J, Mueller POE, Halper J. Differential gene expression in skin RNA of horses affected with degenerative suspensory ligament desmitis.. J Orthop Surg Res 2020 Oct 7;15(1):460.
- Huang LZ, Li YJ, Xie XF, Zhang JJ, Cheng CY, Yamashiro K, Chen LJ, Ma XY, Cheung CM, Wang YS, Zhang CF, Bai YJ, Hou J, Chen XL, Qi Y, Li SS, Sun YY, Mei JP, Cheng Y, Yu WZ, Hu XB, Zhuang FF, Fan L, Lu Y, Sun XH, Zhu XJ, Shen DF, Chan CC, Zhao MW, Yoshimura N, Pang CP, Wong TY, Khor CC, Zhang K, Zhou P, Li XX. Whole-exome sequencing implicates UBE3D in age-related macular degeneration in East Asian populations.. Nat Commun 2015 Apr 15;6:6687.
- Karlsson EK, Sigurdsson S, Ivansson E, Thomas R, Elvers I, Wright J, Howald C, Tonomura N, Perloski M, Swofford R, Biagi T, Fryc S, Anderson N, Courtay-Cahen C, Youell L, Ricketts SL, Mandlebaum S, Rivera P, von Euler H, Kisseberth WC, London CA, Lander ES, Couto G, Comstock K, Starkey MP, Modiano JF, Breen M, Lindblad-Toh K. Genome-wide analyses implicate 33 loci in heritable dog osteosarcoma, including regulatory variants near CDKN2A/B.. Genome Biol 2013 Dec 12;14(12):R132.
- Kemper KE, Saxton SJ, Bolormaa S, Hayes BJ, Goddard ME. Selection for complex traits leaves little or no classic signatures of selection.. BMC Genomics 2014 Mar 28;15(1):246.
- Kendal AR, Layton T, Al-Mossawi H, Appleton L, Dakin S, Brown R, Loizou C, Rogers M, Sharp R, Carr A. Multi-omic single cell analysis resolves novel stromal cell populations in healthy and diseased human tendon.. Sci Rep 2020 Sep 3;10(1):13939.
- Kim B, Yoon JH, Zhang J, Eric Mueller PO, Halper J. Glycan profiling of a defect in decorin glycosylation in equine systemic proteoglycan accumulation, a potential model of progeroid form of Ehlers-Danlos syndrome.. Arch Biochem Biophys 2010 Sep 15;501(2):221-31.
- Kinitz R, Heyne E, Koch LG, Britton SL, Thierbach M, Wildemann B. The Effect of Age and Intrinsic Aerobic Exercise Capacity on the Expression of Inflammation and Remodeling Markers in Rat Achilles Tendons.. Int J Mol Sci 2021 Dec 22;23(1).
- Kuhn M. A short introduction to the caret package. R. Found. Stat. Comput. 1, 1–10.
- Lee SH, Goddard ME, Wray NR, Visscher PM. A better coefficient of determination for genetic profile analysis.. Genet Epidemiol 2012 Apr;36(3):214-24.
- Leegwater PA, Vos-Loohuis M, Ducro BJ, Boegheim IJ, van Steenbeek FG, Nijman IJ, Monroe GR, Bastiaansen JW, Dibbits BW, van de Goor LH, Hellinga I, Back W, Schurink A. Dwarfism with joint laxity in Friesian horses is associated with a splice site mutation in B4GALT7.. BMC Genomics 2016 Oct 28;17(1):839.
- Liu Y, Deng XH, Zhang X, Cong T, Chen D, Hall AJ, Ying L, Rodeo SA. The Role of Indian Hedgehog Signaling in Tendon Response to Subacromial Impingement: Evaluation Using a Mouse Model.. Am J Sports Med 2022 Feb;50(2):362-370.
- Loh PR, Danecek P, Palamara PF, Fuchsberger C, A Reshef Y, K Finucane H, Schoenherr S, Forer L, McCarthy S, Abecasis GR, Durbin R, L Price A. Reference-based phasing using the Haplotype Reference Consortium panel.. Nat Genet 2016 Nov;48(11):1443-1448.
- Ma Y, Ding X, Qanbari S, Weigend S, Zhang Q, Simianer H. Properties of different selection signature statistics and a new strategy for combining them.. Heredity (Edinb) 2015 Nov;115(5):426-36.
- Mero JL, Pool R. Twenty cases of degenerative suspensory ligament desmitis in Peruvian Paso horses. Abstract for AAEP 48, 329–334.
- Mero JL, Scarlett JM. Diagnostic criteria for degenerative suspensory ligament desmitis in Peruvian Paso horses. J. Equine Veterinary Sci. 25 (5), 224–228.
- Metzger J, Distl O. Genetics of Equine Orthopedic Disease.. Vet Clin North Am Equine Pract 2020 Aug;36(2):289-301.
- Meuwissen TH, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps.. Genetics 2001 Apr;157(4):1819-29.
- Momen M, Brounts SH, Binversie EE, Sample SJ, Rosa GJM, Davis BW, Muir P. Selection signature analyses and genome-wide association reveal genomic hotspot regions that reflect differences between breeds of horse with contrasting risk of degenerative suspensory ligament desmitis.. G3 (Bethesda) 2022 Sep 30;12(10).
- Momen M, Mehrgardi AA, Sheikhi A, Kranis A, Tusell L, Morota G, Rosa GJM, Gianola D. Predictive ability of genome-assisted statistical models under various forms of gene action.. Sci Rep 2018 Aug 17;8(1):12309.
- Montesinos-López OA, Gonzalez HN, Montesinos-López A, Daza-Torres M, Lillemo M, Montesinos-López JC, Crossa J. Comparing gradient boosting machine and Bayesian threshold BLUP for genome-based prediction of categorical traits in wheat breeding.. Plant Genome 2022 Sep;15(3):e20214.
- Nakagawa S. A farewell to Bonferroni: the problems of low statistical power and publication bias. Behav. Ecol. 15 (6), 1044–1045.
- Ogutu JO, Schulz-Streeck T, Piepho HP. Genomic selection using regularized linear regression models: ridge regression, lasso, elastic net and their extensions.. BMC Proc 2012 May 21;6 Suppl 2(Suppl 2):S10.
- Oppong RF, Boutin T, Campbell A, McIntosh AM, Porteous D, Hayward C, Haley CS, Navarro P, Knott S. SNP and Haplotype Regional Heritability Mapping (SNHap-RHM): Joint Mapping of Common and Rare Variation Affecting Complex Traits.. Front Genet 2021;12:791712.
- Park T, Casella G. The bayesian lasso. J. Am. Stat. Assoc. 103 (482), 681–686.
- Patron J, Serra-Cayuela A, Han B, Li C, Wishart DS. Assessing the performance of genome-wide association studies for predicting disease risk.. PLoS One 2019;14(12):e0220215.
- Perdry H, Dandine-Roulland C. Genetic Data Handling (QC, GRM, LD, PCA) and linear mixed models. Gaston R package .
- Pérez P, de Los Campos G, Crossa J, Gianola D. Genomic-Enabled Prediction Based on Molecular Markers and Pedigree Using the Bayesian Linear Regression Package in R.. Plant Genome 2010;3(2):106-116.
- Pérez P, de los Campos G. Genome-wide regression and prediction with the BGLR statistical package.. Genetics 2014 Oct;198(2):483-95.
- Petersen JL, Mickelson JR, Rendahl AK, Valberg SJ, Andersson LS, Axelsson J, Bailey E, Bannasch D, Binns MM, Borges AS, Brama P, da Câmara Machado A, Capomaccio S, Cappelli K, Cothran EG, Distl O, Fox-Clipsham L, Graves KT, Guérin G, Haase B, Hasegawa T, Hemmann K, Hill EW, Leeb T, Lindgren G, Lohi H, Lopes MS, McGivney BA, Mikko S, Orr N, Penedo MC, Piercy RJ, Raekallio M, Rieder S, Røed KH, Swinburne J, Tozaki T, Vaudin M, Wade CM, McCue ME. Genome-wide analysis reveals selection for important traits in domestic horse breeds.. PLoS Genet 2013;9(1):e1003211.
- Plaas A, Sandy JD, Liu H, Diaz MA, Schenkman D, Magnus RP, Bolam-Bretl C, Kopesky PW, Wang VM, Galante JO. Biochemical identification and immunolocalizaton of aggrecan, ADAMTS5 and inter-alpha-trypsin-inhibitor in equine degenerative suspensory ligament desmitis.. J Orthop Res 2011 Jun;29(6):900-6.
- nR Core Team. R: A language and environment for statistical computing. Vienna, A: R Foundation for Statistical Computing.
- Rai MF, Cai L, Tycksen ED, Chamberlain A, Keener J. RNA-Seq analysis reveals sex-dependent transcriptomic profiles of human subacromial bursa stratified by tear etiology.. J Orthop Res 2022 Dec;40(12):2713-2727.
- Raudvere U, Kolberg L, Kuzmin I, Arak T, Adler P, Peterson H, Vilo J. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update).. Nucleic Acids Res 2019 Jul 2;47(W1):W191-W198.
- Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, Müller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves.. BMC Bioinformatics 2011 Mar 17;12:77.
- Salminen A. Aryl hydrocarbon receptor (AhR) reveals evidence of antagonistic pleiotropy in the regulation of the aging process.. Cell Mol Life Sci 2022 Aug 20;79(9):489.
- Sargolzaei M, Iwaisaki H, Colleau JJ. CFC: A tool for monitoring genetic diversity. Proc. 8th World Congr. Genet. Appl. Livest. Prod. CD-ROM Commun. 27–28, 13–18.
- Schweitzer R, Chyung JH, Murtaugh LC, Brent AE, Rosen V, Olson EN, Lassar A, Tabin CJ. Analysis of the tendon cell fate using Scleraxis, a specific marker for tendons and ligaments.. Development 2001 Oct;128(19):3855-66.
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks.. Genome Res 2003 Nov;13(11):2498-504.
- Shirali M, Pong-Wong R, Navarro P, Knott S, Hayward C, Vitart V, Rudan I, Campbell H, Hastie ND, Wright AF, Haley CS. Regional heritability mapping method helps explain missing heritability of blood lipid traits in isolated populations.. Heredity (Edinb) 2016 Mar;116(3):333-8.
- Strong DI. The use of a whole genome scan to find a genetic marker for degenerative suspensory ligament desmitis in the Peruvian Paso horse. Master’s thesis University of Kentucky.
- Szpiech ZA, Hernandez RD. selscan: an efficient multithreaded program to perform EHH-based scans for positive selection.. Mol Biol Evol 2014 Oct;31(10):2824-7.
- Thomopoulos S, Parks WC, Rifkin DB, Derwin KA. Mechanisms of tendon injury and repair.. J Orthop Res 2015 Jun;33(6):832-9.
- Tibshirani R. Regression shrinkage and selection via the lasso. J. R. Stat. Soc. Ser. B Methodol. 58 (1), 267–288.
- VanRaden PM. Efficient methods to compute genomic predictions.. J Dairy Sci 2008 Nov;91(11):4414-23.
- Verity R, Collins C, Card DC, Schaal SM, Wang L, Lotterhos KE. minotaur: A platform for the analysis and visualization of multivariate results from genome scans with R Shiny.. Mol Ecol Resour 2017 Jan;17(1):33-43.
- Wang Z, Kim SS, Hutton WC, Yoon ST. E-cadherin upregulates expression of matrix macromolecules aggrecan and collagen II in the intervertebral disc cells through activation of the intracellular BMP-Smad1/5 pathway.. J Orthop Res 2012 Nov;30(11):1746-52.
- Xu B, Ye Y, Nie L. An improved random forest classifier for image classification. 2012 IEEE International Conference on Information and Automation 795–800.
- Young M, Moshood O, Zhang J, Sarbacher CA, Mueller POE, Halper J. Does BMP2 play a role in the pathogenesis of equine degenerative suspensory ligament desmitis?. BMC Res Notes 2018 Sep 18;11(1):672.