- Research
- Open access
- Published:
Identification of genetic variants associated with anterior cruciate ligament rupture and AKC standard coat color in the Labrador Retriever
BMC Genomic Data volume 24, Article number: 60 (2023)
Abstract
Canine anterior cruciate ligament (ACL) rupture is a common complex disease. Prevalence of ACL rupture is breed dependent. In an epidemiological study, yellow coat color was associated with increased risk of ACL rupture in the Labrador Retriever. ACL rupture risk variants may be linked to coat color through genetic selection or through linkage with coat color genes. To investigate these associations, Labrador Retrievers were phenotyped as ACL rupture case or controls and for coat color and were single nucleotide polymorphism (SNP) genotyped. After filtering, ~ 697 K SNPs were analyzed using GEMMA and mvBIMBAM for multivariate association. Functional annotation clustering analysis with DAVID was performed on candidate genes. A large 8 Mb region on chromosome 5 that included ACSF3, as well as 32 additional SNPs, met genome-wide significance at P < 6.07E-7 or Log10(BF) = 3.0 for GEMMA and mvBIMBAM, respectively. On chromosome 23, SNPs were located within or near PCCB and MSL2. On chromosome 30, a SNP was located within IGDCC3. SNPs associated with coat color were also located within ADAM9, FAM109B, SULT1C4, RTDR1, BCR, and RGS7. DZIP1L was associated with ACL rupture. Several significant SNPs on chromosomes 2, 3, 7, 24, and 26 were located within uncharacterized regions or long non-coding RNA sequences. This study validates associations with the previous ACL rupture candidate genes ACSF3 and DZIP1L and identifies novel candidate genes. These variants could act as targets for treatment or as factors in disease prediction modeling. The study highlighted the importance of regulatory SNPs in the disease, as several significant SNPs were located within non-coding regions.
Introduction
The cruciate ligaments within the knee provide stability and oppose rotational and translational forces. In humans, the anterior cruciate ligament (ACL) is prone to rupture mostly through non-contact injury, especially in young female athletes [1]. In dogs, the cranial cruciate ligament is anatomically equivalent to the ACL and is also vulnerable to non-contact injury [2]. ACL rupture is the most common cause for canine lameness [3]. The underlying disease mechanism for most cases is described as the progressive tearing of ligament fibers in the presence of stifle synovitis [4, 5]. A combination of genetic and environmental factors defines an individual’s risk of complex disease development. Established genetic risk factors include breed, joint conformation, and joint immune responses, while age, obesity, and other extrinsic factors also have roles in disease development [6,7,8,9]. Prevalence is breed dependent; ACL rupture is more common in large breed dogs, such as Newfoundlands, Rottweilers, and Labrador Retrievers [9]. Heritability estimates for canine ACL rupture range from 0.27 to 0.89 [10,11,12,13]. High prevalence in a breed is a consequence of intense selection leading to a concentration of risk loci [14]. Linkage disequilibrium (LD) is extensive in dogs, which enhances the statistical power of genome-wide association studies (GWAS) using single nucleotide polymorphism (SNP) markers as an approach for variant discovery [15]. ACL rupture in the Labrador Retriever is a polygenic complex disease whereby many loci with small and moderate effects influence disease risk [12, 16]. Several candidate genetic variants have been identified [12, 17,18,19], but given its polygenic nature, it is likely that additional variants remain to be discovered.
Labrador Retrievers have three AKC recognized coat colors: black, chocolate, and yellow. Risk of ACL rupture is increased in yellow Labrador Retrievers compared to black and chocolate dogs [20]. Inheritance of coat color is controlled by two genes: MC1R and TYRP1. Yellow coat color is inherited in an autosomal recessive manner and is caused by a loss of function mutation in MC1R [21, 22]. In dogs that are homozygous for MC1R mutations, pheomelanin is produced, creating a yellow coat color, while dogs with at least one wildtype MC1R allele produce eumelanin resulting in a black or chocolate phenotype. In these dogs, color is determined by the TYRP1 gene. Chocolate color is autosomal recessive to black and is seen with any of three TYRP1 mutations [23]. Through LD, selection for coat color may have inadvertently selected risk variants for other phenotypes. For example, research suggests that chocolate Labrador Retrievers have shorter lifespans and are affected with skin and ear disease at higher rates than black or yellow Labrador Retrievers [24]. Behavioral differences have been associated with coat color as well [25]. Alternatively, genes that regulate coat color may have additional unknown direct biologic effects on ACL rupture risk. Studies have shown that MC1R has an important role in regulating inflammatory pathways [26,27,28], and variants associate with osteoarthritis [29].
High LD and diversity within and between dog breeds enhance use of GWAS as an approach for causal variant discovery. Joint analysis of multiple phenotypes is a valuable technique that can increase statistical power to detect small and moderate associations expected with complex traits, such as ACL rupture [18, 30], and detect associations with ACL rupture, coat color, or both phenotypes in this study. Such knowledge will advance analysis of candidate variants that could improve genomic prediction of disease and further mechanistic study of influential biological pathways.
Materials and methods
Research approach
A within-breed GWAS in the Labrador Retriever was performed by analyzing SNPs from case and control dogs. All procedures were performed in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health and the American Veterinary Medical Association and with approval from the Animal Care Committee of the University of Wisconsin-Madison (protocols V1070, V5463). Informed consent of each owner was obtained before participation in the study. Phenotypes included ACL rupture case/control status and coat color. Contact with the owners of control dogs was maintained to ensure accurate phenotyping and coding was updated if a dog became affected. Preparation of the manuscript conformed with the ARRIVE guidelines.
Inclusion and exclusion criteria for phenotyping
Most Labrador Retrievers exhibit signs of ACL rupture by 8 years old [31]. Lameness, knee instability, and ligament damage are usually confirmed during surgical treatment. Dogs were excluded as a case if contact injury was diagnosed from the clinical history. Phenotype negative controls were dogs ≥ 8 years old with normal orthopaedic exams and normal knee radiographs.
Clinical examination of each dog was performed, including knee palpation for instability. Bilateral weight-bearing radiographs were reviewed and graded for stifle effusion and osteophytosis [32]. Dog age, sex, weight, neuter status, and coat color were also recorded [18].
DNA isolation, SNP genotyping, and quality control
DNA was isolated from buffy coat leukocytes from EDTA blood or cheek swab saliva. DNA was isolated using standard reagents (Blood – Puregene Cell Core Kit, Qiagen, Germantown, MD; saliva – DNA Genotek prepIT-L2P, Ottawa, ON, Canada). Genotyping was performed using the Illumina Canine HD BeadChip which has ~ 220,000 SNPs mapped to CanFam3.1. SNPs were then imputed to the Thermofisher Axiom HD Canine 770 K array using Beagle 5.0 [19, 33] and a multi-breed reference panel that consisted of 646 purebred dogs of various breeds that were genotyped using the Axiom Canine HD 770 K array, including 96 Labrador Retrievers. Before imputation of the ACL rupture SNP data, we validated our method using a group of 22 Labrador Retrievers with whole genome sequence (WGS) data. Illumina SNPs were extracted from WGS data and imputed to denser set of Axiom SNPs using Beagle 5.0 with the reference panel, a window size of 3 cM with 1 cM overlap, and effective population size of 100. Imputation accuracy was 96%. Using PLINK1.9 [34], SNPs with a minor allele frequency (MAF) < 0.005 (to analyze potential rare variants) and a call rate < 95%, and dogs with a call rate of < 95% were removed. SNPs with deviations from Hardy-Weinberg proportions at P < 1E-7 were also filtered out.
Multivariate GWAS
ACL rupture case-control and coat color phenotypes were used for multivariate association using the linear mixed model approach implemented in GEMMA [35]. P-values are calculated to measure support for each model compared to the null (no association). The value for genome-wide significance was calculated using a Bonferroni correction for the number of haplotype blocks in the genome using PLINK [34] and was determined to be P < 6.07E-7. Additionally, a Bayesian statistical model was used with mvBIMBAM [36] to identify associations between genotypes and phenotypes (directly associated, indirectly associated, or unassociated with one or both phenotypes). Bayes Factors were evaluated for evidence of association with the multivariate phenotype. SNPs with a Log10 Bayes Factor > 3 were considered to have moderate evidence of association, and SNPs with a Log10 Bayes Factor > 6 were considered to have strong evidence of association. Marginal posterior probabilities of associated SNPs were evaluated to determine which phenotypes may be influencing the association. Sex, age, neuter status, and weight were included as covariates in both models. For mvBIMBAM, case-control phenotypes were residuals of multiple logistic regression against the covariates.
Regions with evidence of association with either GEMMA or mvBIMBAM analysis were evaluated using the CanFam3.1 Broad Improved Canine Annotation catalog in the UCSC Table Browser to identify associated genes. LocusZoom [37] and triangle heat map plotting with the ‘gpart’ R package and the Big-LD algorithm [38, 39] were also used to investigate haplotype structure and candidate genes in selected regions.
Functional annotation clustering
A list of genes that were within a ±25 kb flanking region of each significant SNP was created using the canFam3.1 Broad Improved Canine Annotation catalog in the UCSC Table Browser. Functional annotation clustering was then performed using DAVID [40, 41]. Functional clusters with significant P-values were evaluated for biological relevance to the ACL rupture and coat color phenotypes.
Results
Phenotype data were collected from 367 dogs. After quality control, the final dataset included 696,846 SNPs from 148 cases and 219 controls. There were 41 cases and 96 controls with a black coat color, 33 cases and 51 controls with a chocolate coat color, and 73 cases and 73 controls with a yellow coat color. Black dogs had a decreased risk of ACL rupture (odds ratio = 0.50, P = 0.0024), while yellow dogs had an increased risk (odds ratio = 1.99, P = 0.0017). Chocolate color did not influence ACL rupture risk (Table 1).
The multivariate GWAS using GEMMA provided evidence of association (P < 6.02E-7) for 337 SNPs (Figs. 1 and 2; Table 2, Supplementary File S1). Of these, 331 were located on chromosome 5, including the most significant SNP (P = 3.95E-45), located within the TCF25 gene, which is a transcription factor that is important in embryonic development. Significant SNPs on chromosome 5 spanned a large region of > 8 Mb that included more than 100 genes. The 6 other significant SNPs were located on chromosomes 7, 23, 24, and 30. SNPs were within or near non-coding or uncharacterized regions on chromosome 7 and 24 (Table 2). On chromosome 23, the association was near the PCCB and MSL2 genes. On chromosome 30, the association was within the IGDCC3 gene. To better discern nearby genes and haplotype structure in the large region chromosome 5, a LocusZoom plot was built around the most significant SNP (P = 3.95E-45) at chr5:63697949 (Fig. 3). There were many SNPs in LD within the ANKRD11 gene. There were also SNPs in LD within MC1R, PIK3CD, DPEP1, SPG7, ACSF3, and CBFA2T3 (Fig. 3). Additionally, a triangle plot was built spanning a ~ 3 Mb region (62.3-65.3 Mb) which revealed several areas of LD, including a haplotype block that contained TCF25 and MC1R.
The multivariate analysis using mvBIMBAM revealed a total of 589 SNPs directly associated with either coat color, ACL rupture, or both (Fig. 1, Supplementary File S2). 237 of these SNPs displayed strong evidence of association (Log10 BF≥6.0). Of the associated SNPs, 556 were on chromosome 5, and were within or near 56 genes. The SNP with the highest significance (Log10 BF = 31.30) was located on chromosome 5 at the same position as the most significant SNP revealed in GEMMA analysis (Table 3) and was associated with coat color. There were nine additional loci associated with coat color found on chromosomes 2, 3, 7, 10, 16, and 26 (Table 3). The associations on chromosomes 2 and 3 were within noncoding regions. On chromosome 7, an association was identified within RGS7. Chromosome 10 associations were within FAM109B and SULT1C4. Three SNPs were associated ADAM9 on chromosome 16. On chromosome 26, there were associations within lncRNA sequences, and the RTDR1 and BCR genes. We identified one DZIP1L SNP on chromosome 23 that was associated with ACL rupture. Additionally, there was one SNP on chromosome 24 that was associated with both color and ACL rupture that was within a noncoding region. Functional annotation clustering with DAVID did not identify significant pathway enrichment.
Discussion
Multivariate GWAS improves power to detect loci with weaker associations with disease risk. Multivariate analysis of ACL rupture case and control Labrador Retrievers of differing coat color revealed associations with several ACL rupture candidate genes after GEMMA and mvBIMBAM analysis, consistent with a polygenic disease. This study validated ACSF3 and DZIP1L [19] as candidate genes, discovered other novel candidate loci, and revealed association with many genes in a chromosome 5 locus.
Considering coat color as a potential risk factor in the development of ACL rupture is a novel approach. Although the exact mechanisms remain unknown, coat color has been found to influence the risk of other disease processes and behavior patterns in the dog [24, 25, 42]. Since ACL rupture is a complex disease, it is possible that coat color genes may have small effects on disease risk, or that other genes that influence risk are inherited together with coat color genes through LD.
The region of the genome that displayed the strongest association was on chromosome 5. There were numerous SNPs that met genome wide significance in both GEMMA and mvBIMBAM analysis. This locus contained many associated genes with proximity to the MC1R gene. In both multivariate analyses, the most significant SNP was within the TCF25 gene. This gene encodes a transcription factor that acts as a transcriptional repressor and is important during embryonic development. It is not currently known to have effects on ACL rupture or coat color, or to directly influence biological pathways that precede disease development. Interestingly, earlier research [43] and our analysis suggests that TCF25 is in LD with MC1R. Mutations in MC1R control yellow versus chocolate or black coat color in Labrador Retrievers [21, 22]. MC1R is also known to be expressed within articular cartilage [44], and melanocortins promote anti-inflammatory states within joints [26, 28]. Nonfunctional MC1 receptors have been linked to the development of osteoarthritis in mice [29]. Knee osteoarthritis is an inflammatory process that is also associated with ACL rupture in dogs [4]. LD between TCF25 and MC1R may indicate that the effect of this locus is mediated by MC1R. Defective MC1 receptors may promote an inflammatory state within the knee that leads to cruciate ligament fiber rupture. It should be noted that no genome-wide associations were found within the TYRP1 gene, which also controls coat color in Labrador Retrievers, so it is unlikely that the association near MC1R is an unintended consequence of including coat color as an additional phenotype.
There were many other significant SNP associations within the haplotype block on chromosome 5 that were located within or near other candidate ACL rupture genes. Interestingly, one of these genes, ACSF3, has been previously associated with ACL rupture [19]. ACSF3 is differentially expressed in ligament [19] and may influence risk of rheumatoid arthritis [45]. ANKRD11 is another interesting gene in LD with the most significant SNP. This gene regulates cell proliferation and apoptosis. Maladaptive responses to injury in joints could lead to increased inflammation or structural changes that compromise joint homeostasis. Other genes in LD in this locus, including DPEP1, LZIC, SLC7A5, PIK3CD, RPL13, and TNFRSF25, have been associated with OA or are involved with processes that could link them to ACL rupture pathogenesis [46,47,48,49,50,51,52,53,54]. PIK3CD and TNFRSF25 regulate lymphocyte development, which could play a role in ACL rupture since lymphoplasmacytic inflammation is seen within affected knees [55,56,57]. CLSTN1 is associated with weight and obesity, which are risk factors for ACL rupture development [58].
In both analyses, significant SNPs were located within multiple non-coding regions or lncRNA sequences. These sequences are thought to have regulatory effects on transcription, translation, or post-translational portions of protein production [59, 60]. Effects of regulatory SNPs in the risk of complex canine diseases are not understood. In complex diseases such as ACL rupture, many small-effect SNPs combine to influence disease risk. It is plausible that regulatory mechanisms play a role in gene expression and ultimately influence disease risk. Gene regulation effects on protein production and cell signaling have been implicated in changes to the extracellular matrix, ligament, and cartilage homeostasis, as well as ACL rupture pathophysiology [61]. Other than those on chromosome 5 and the associations within the uncharacterized or non-coding regions, GEMMA analysis identified several additional ACL rupture candidate genes. On chromosome 23, there was a significant SNP located near the MSL2 gene. MSL2 is responsible for histone acetylation and gene activation, as well as cellular responses to damage, such as apoptosis [62], and could modulate inflammatory responses leading to ligament fiber rupture [4, 5].
Due to the different statistical approaches of GEMMA and mvBIMBAM programs, it is not surprising that unique SNPs were discovered in each analysis. mvBIMBAM revealed significant associations with 26 SNPs that were not identified by GEMMA, most being associated with coat color. A chromosome 10 locus included the genes FAM109B and SULT1C4. FAM109B plays a role in endocytic trafficking [63]. It is expressed in many tissues including skin melanocytes and keratinocytes [64], which could explain association with coat color. RGS7 on chromosome 7, SULT1C4 on chromosome 10, ADAM9 on chromosome 16, and RTDR1 and BCR on chromosome 26 have no known association with coat color.
DZIP1L on chromosome 23 was associated with ACL rupture only, which validates earlier work [19]. DZIP1L encodes a protein found in the transition zone of cilia, and mutations have been associated with autosomal recessive polycystic kidney disease, as well as craniofacial deformities and polydactyly [65]. Cilia are found in chondrocytes, fibroblasts, and other connective tissue cells. Primary cilia of connective tissues are within the extracellular matrix and transduce chemical stimuli, mechanical stimuli, or respond to growth factors to control homeostasis, fibroblast migration, and cell cycling [66]. DZIP1L mutations resulting in ciliary dysfunction could influence ligament mechanotransduction and cause dysfunction in fibroblast homeostasis.
Through the inclusion of coat color as an additional phenotype, several genes were surprisingly associated with Labrador Retriever coat color. There are two genes that are known to control coat color in Labrador Retrievers [21, 23] as well as many other genes that play a role in coloring patterns, fur length, hair structure, and other related characteristics in dogs [41, 67,68,69]. Besides these, there have been numerous other loci with effects on coat color in other animal models [70] that have not been studied extensively in dogs. Therefore, it is plausible that some of the coat color gene associations in the current study may influence pathways that determine phenotypic appearance of dogs in ways that are not currently understood. Future research is warranted in this area to determine effects on skin, hair follicles, strands, or related structures that may explain the associations found in this study. There were many unique additional observations in this study. Differences in SNP discovery across studies are likely due to a difference in selection of subject breeds, study design, and analytical approaches [12, 17,18,19, 71].
In conclusion, GWAS using multivariate linear mixed model and Bayesian model approaches has identified several novel variants associated with ACL rupture and coat color in the Labrador Retriever. The study has also identified two variants, ACSF3 and DZIP1L, that were validated from previous studies [19], suggesting that these genes merit additional investigation. Associated genes in this study have effects on bone and cartilage pathology, inflammatory pathways, metabolism, development, and gene expression and regulation, supporting the complexity of ACL rupture [12]. Although coat color has been linked to many other disease states, this is the first study to examine the relationship between coat color and ACL rupture in the Labrador Retriever. Canine ACL rupture is an important model for human ACL rupture [2], which is also a heritable disease [72] that leads to fatigue injury of ACL fibers and eventual non-contact rupture [2, 73]. Results from this study have translational value for the development of treatment and prevention strategies in both species, particularly polygenic risk score prediction of disease risk in both species.
Data Availability
The SNP data set used for this analysis cannot be fully shared publicly because of restrictions relating to development of a commercial genetic screening test for ACL rupture in the dog at the University of Wisconsin-Madison. ACL rupture phenotypes are proprietary. The datasets generated and analyzed in this study are available in the Dryad repository at https://doi.org/10.5061/dryad.brv15dvfw.
References
Agel J, Rockwood T, Klossner D. Collegiate ACL injury rates across 15 sports: National collegiate athletic association injury surveillance system data update (2004–2005 through 2012–2013). Clin J Sport Med. 2016;26:518–23.
Binversie EE, Walczak BE, Cone SG, Baker LA, Scerpella TA, Muir P. Canine ACL rupture: a spontaneous large animal model of human ACL rupture. BMC Musculoskelet Disord. 2022;23:116.
Engdahl K, Hanson J, Bergström A, Bonnett B, Höglund O, Emanuelson U. The epidemiology of stifle joint Disease in an insured Swedish dog population. Vet Rec. 2021;189:e197.
Bleedorn J, Greuel E, Manley P, Schaefer S, Markel M, Holzman G, Muir P. Synovitis in dogs with stable stifle joints and incipient cranial cruciate ligament rupture: a cross-sectional study. Vet Surg. 2011;40:531–43.
Doring AK, Junginger J, Hewicker-Trautwein M. Cruciate ligament degeneration and stifle joint synovitis in 56 dogs with intact cranial cruciate ligaments: correlation of histological findings and numbers and phenotypes of inflammatory cells with age, body weight and breed. Vet Immunol Immunopathol. 2018;196:5–13.
Comerford E, Smith K, Hayashi K. Update on the aetiopathogenesis of canine cranial cruciate ligament Disease. Vet Comp Orthop Traumatol. 2011;24:91–8.
Healey E, Murphy RJ, Hayward JJ, Castelhano M, Boyko AR, Hayashi K, Krotscheck U, Todhunter RJ. Genetic mapping of distal femoral, stifle, and tibial radiographic morphology in dogs with cranial cruciate ligament Disease. PLoS ONE. 2019;14:e0223094.
Sellon DC, Marcellin-Little DJ. Risk factors for cranial cruciate ligament rupture in dogs participating in canine agility. BMC Vet Res. 2022;18:39.
Witsberger T, Villamil J, Schultz L, Hahn A, Cook J. Prevalence of and risk factors for hip dysplasia and cranial cruciate ligament deficiency in dogs. J Am Vet Med Assoc. 2008;232:1818–24.
Nielen AL, Janss LL, Knol BW. Heritability estimations for diseases, coat color, body weight, and height in a birth cohort of Boxers. Am J Vet Res. 2001;62:1198–2006.
Wilke V, Conzemius M, Kinghorn B, Macrossan P, Cai W, Rothschild M. Inheritance of rupture of the cranial cruciate ligament in Newfoundlands. J Am Vet Med Assoc. 2006;228:61–4.
Baker L, Kirkpatrick B, Rosa G, Gianola D, Valente B, Sumner J, Baltzer W, Hao Z, Binversie EE, Volstad N, Piazza A, Sample SJ, Muir P. Genome-wide association analysis in dogs implicates 99 loci as risk variants for anterior cruciate ligament rupture. PLoS ONE. 2017;12:e0173810.
Cook SR, Conzemius MG, McCue ME, Ekenstedt KJ. SNP-based heritability and genetic architecture of cranial cruciate ligament rupture in Labrador Retrievers. Anim Genet. 2020;51:824–8.
Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, et al. Genome sequence, comparative analysis, and haplotype structure of the domestic dog. Nature. 2005;438:803–19.
Karlsson EK, Lindblad-Toh K. Leader of the pack: gene mapping in dogs and other model organisms. Nat Rev Genet. 2008;9:713–25.
Baker LA, Momen M, Chan K, Bollig N, Lopes FB, Rosa GJM, Todhunter RJ, Binversie EE, Sample SJ, Muir P. Bayesian and machine learning models for genomic prediction of anterior cruciate ligament rupture in the canine model. (Bethesda). 2020;G3:10:2619–28.
Baird AE, Carter SD, Innes JF, Ollier W, Short A. Genome-wide association study identifies genomic regions of association for cruciate ligament rupture in Newfoundland dogs. Anim Genet. 2014;45:542–9.
Baker L, Rosa G, Hao A, Piazza A, Hoffman C, Binversie E, Sample SJ, Muir P. Multivariate genome-wide association analysis identifies novel and relevant variants associated with anterior cruciate ligament rupture risk in the dog model. BMC Genet. 2018;19:39.
Baker LA, Momen M, McNally R, Berres ME, Binversie EE, Sample SJ, Muir P. Biologically enhanced genome-wide association study provides further evidence for candidate loci and discovers novel loci that influence risk of anterior cruciate ligament rupture in a dog model. Front Genet. 2021;2:593515.
Terhaar H, Muir P, Baker L, Binversie E, Chi J, Sample SJ. Contribution of habitual activity to cruciate ligament rupture in Labrador Retrievers. Vet Comp Orthop Traumatol. 2020;33:82–8.
Everts RE, Rothuizen J, van Oost BA. Identification of a premature stop codon in the melanocyte-stimulating hormone receptor gene (MC1R) in Labrador and Golden retrievers with yellow coat colour. Anim Genet. 2000;31:194–9.
Newton JM, Wilkie AL, He L, Jordan SA, Metallinos DL, Holmes NG, Jackson IJ, Barsh GS. Melanocortin 1 receptor variation in the domestic dog. Mamm Genome. 2000;11:24–30.
Schmutz SM, Berryere TG, Goldfinch AD. TYRP1 and MC1R genotypes and their effects on coat color in dogs. Mamm Genome. 2002;13:380–7.
McGreevy P, Wilson B, Mansfield C, Brodbelt D, Church D, Dhand N. Soares Magalhães RJ, O’Neill DG. Labrador retrievers under primary veterinary care in the UK: demography, mortality, and disorders. Canine Genet Epidemiol. 2018;5:8.
van Rooy D, Wade CM. Association between coat colour and the behaviour of Australian Labrador retrievers. Canine Genet Epidemiol. 2019;6:10.
Kang L, McIntyre KW, Gillooly KM, Yang Y, Haycock J, Roberts S, Khanna A, Herpin TF, Yu G, Wu X, Morton GC, Tuerdi H, Koplowitz B, Walker SG, Wardwell-Swanson J, Macor JE, Lawrence RM, Carlson KE. A selective small molecule agonist of the melanocortin-1 receptor inhibits lipopolysaccharide-induced cytokine accumulation and leukocyte infiltration in mice. J Leukoc Biol. 2006;80:897–904.
Carter DW, Sood RF, Seaton ME, Muffley LA, Honari S, Hocking AM, Arbabi SA, Gibran NS. MC1R gene polymorphisms are associated with dysfunctional immune responses and wound Infection after burn injury. J Surg Res. 2018;231:448–52.
Montero-Melendez T, Nagano A, Chelala C, Filer A, Buckley CD, Perretti M. Therapeutic senescence via GPCR activation in synovial fibroblasts facilitates resolution of arthritis. Nat Commun. 2020;11:745.
Lorenz J, Seebach E, Hackmayer G, Greth C, Bauer R, Kleinschmidt K, Bettenworth D, Böhm M, Grifka J, Grässel S. Melanocortin 1 receptor-signaling deficiency results in an articular cartilage phenotype and accelerates pathogenesis of surgically induced murine osteoarthritis. PLoS ONE. 2014;9:e105858.
Robinson MR, Wray NR, Visscher PM. Explaining additional genetic variation in complex traits. Trends Genet. 2014;30:124–32.
Reif U, Probst CW. Comparison of tibial plateau angles in normal and cranial cruciate deficient stifles of Labrador retrievers. Vet Surg. 2003;32:385–9.
Chuang C, Ramaker MA, Kaur S, Csomos RA, Kroner KT, Bleedorn JA, Schaefer SL, Muir P. Radiographic risk factors for contralateral rupture in dogs with unilateral cranial cruciate ligament rupture. PLoS ONE. 2014;9:e106389.
Browning BL, Zhou Y, Browning SR. A one-penny imputed genome from next-generation reference panels. Am J Human Genet. 2018;103:338–48.
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.
Zhou X, Stephens M. Efficient algorithms for multivariate linear mixed models in genome-wide association studies. Nat Methods. 2014;11:407–9.
Stephens M. A unified framework for association analysis with multiple related phenotypes. PLoS ONE. 2013;8:e65245.
Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics. 2010;26:2336–7.
Kim SA, Brossard M, Roshandel D, Paterson AD, Bull SB, Yoo YJ. Gpart: human genome partitioning and visualization of high-density SNP data by identifying haplotype blocks. Bioinformatics. 2019;35:4419–21.
R Core Team. (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/.
Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Brancalion L, Haase B, Wade CM. Canine coat pigmentation genetics: a review. Anim Genet. 2022;53:3–34.
Ishida K, Sato JJ, Kinoshita G, Hosoda T, Kryukov AP, Suzuki H. Evolutionary history of the sable (Martes zibellina brachyura) on Hokkaido inferred from mitochondrial cytb and nuclear Mc1r and Tcf25 gene sequences. Acta Theriol. 2013;58:13–24.
Grässel S, Opolka A, Anders S, Straub RH, Grifka J, Luger TA, Böhm M. The melanocortin system in articular chondrocytes: melanocortin receptors, pro-opiomelanocortin, precursor proteases, and a regulatory effect of alpha-melanocyte-stimulating hormone on proinflammatory cytokines and extracellular matrix components. Arthritis Rheum. 2009;60:3017–27.
Julia A, Absher D, Lopez-Lasanta M, Palau N, Pluma A, Waite Jones L, Glossop JR, Farrell WE, Myers RM, Marsal S. Epigenome-wide association study of rheumatoid arthritis identifies differentially methylated loci in B cells. Hum Mol Genet. 2017;26:2803–11.
Tachmazidou I, Hatzikotoulas K, Southam L, Esparza-Gordillo J, Haberland V, Zheng J, Johnson T, Koprulu M, Zengini E, Steinberg J, Wilkinson JM, Bhatnagar S, Hoffman JD, Buchan N, Süveges D, arcOGEN Consortium, Yerges-Armstrong L, Smith GD, Gaunt TR, Scott RA, McCarthy LC, Zeggini E. Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data. Nat Genet. 2019;51:230–6.
Zhang T, Yang Y, Yin X, Ni J, Dong Z, Li C, Lu W. Two loss-of-function ANKRD11 variants in Chinese patients with short stature and a possible molecular pathway. Am J Med Genet A. 2021;185:710–8.
Zhu M, Chen M, Zuscik M, Wu Q, Wang YJ, Rosier RN, O’Keefe RJ, Chen D. Inhibition of beta-catenin signaling in articular chondrocytes results in articular cartilage destruction. Arthritis Rheum. 2008;58:2053–64.
Yoon BR, Oh YJ, Kang SW, Lee EB, Lee WW. Role of SLC7A5 in metabolic reprogramming of human monocyte/macrophage immune responses. Front Immunol. 2018;9:53.
Xu J, Jiang C, Cai Y, Guo Y, Wang X, Zhang J, Xu J, Xu K, Zhu W, Wang S, Zhang F, Geng M, Han Y, Ning Q, Xu P, Meng L, Lu S. Intervening upregulated SLC7A5 could mitigate inflammatory mediator by mTOR-P70S6K signal in rheumatoid arthritis synoviocytes. Arthritis Res Ther. 2020;22:200.
Bartok B, Hammaker D, Firestein GS. Phosphoinositide 3-kinase δ regulates migration and invasion of synoviocytes in rheumatoid arthritis. J Immunol. 2014;192:2063–70.
Bartok B, Boyle DL, Liu Y, Ren P, Ball ST, Bugbee WD, Rommel C, Firestein GS. PI3 kinase δ is a key regulator of synoviocyte function in rheumatoid arthritis. Am J Pathol. 2012;180:1906–16.
Le Caignec C, Ory B, Lamoureux F, O’Donohue MF, Orgebin E, Lindenbaum P, Téletchéa S, Saby M, Hurst A, Nelson K, Gilbert SR, Wilnai Y, Zeitlin L, Segev E, Tesfaye R, Nizon M, Cogne B, Bezieau S, Geoffroy L, Hamel A, Mayrargue E, de Courtivron B, Decock-Giraudaud A, Charrier C, Pichon O, Retière C, Redon R, Pepler A, McWalter K, Da Costa L, Toutain A, Gleizes PE, Baud’huin M, Isidor B. RPL13 variants cause spondyloepimetaphyseal dysplasia with severe short stature. Am J Hum Genet. 2019;105:1040–7.
Collins FL, Williams JO, Bloom AC, Singh RK, Jordan L, Stone MD, Stone MD, McCabe LR, Wang ECY, Williams AS. CCL3 and MMP-9 are induced by TL1A during death receptor 3 (TNFRSF25)-dependent osteoclast function and systemic bone loss. Bone. 2017;97:94–104.
Nunes-Santos CJ, Uzel G, Rosenzweig SD. PI3K pathway defects leading to immunodeficiency and immune dysregulation. J Allergy Clin Immunol. 2019;143:1676–87.
Hashiramoto A, Konishi Y, Murayama K, Kawasaki H, Yoshida K, Tsumiyama K, Tanaka K, Mizuhara M, Shiotsuki T, Kitamura H, Komai K, Kimura T, Yagita H, Shiozawa K, Shiozawa S. A variant of death-receptor 3 associated with rheumatoid arthritis interferes with apoptosis-induction of T cell. J Biol Chem. 2018;293:1933–43.
Valatas V, Kolios G, Bamias G. TL1A (TNFSF15) and DR3 (TNFRSF25): a co-stimulatory system of cytokines with diverse functions in gut mucosal immunity. Front Immunol. 2019;10:583.
Hohos NM, Smith AK, Kilaru V, Park HJ, Hausman DB, Bailey LB, Lewis RD, Phillips BG, Meagher RB. CD4+ and CD8+ T-cell-specific DNA cytosine methylation differences associated with obesity. Obes (Silver Spring). 2018;26:1312–21.
Zhang F, Lupski JR. Non-coding genetic variants in human disease. Hum Mol Genet. 2015;24:R102–10.
Le Beguec C, Wucher V, Lagoutte L, Cadieu E, Botherel N, Hedan B, De Brito C, Guillory AS, André C, Derrien T, Hitte C. Characterisation and functional predictions of canine long non-coding RNAs. Sci Rep. 2018;8:13444.
Rahim M, Mannion S, Klug B, Hobbs H, van der Merwe W, Posthumus M, Collins M, September AV. Modulators of the extracellular matrix and risk of anterior cruciate ligament ruptures. J Sci Med Sport. 2017;20:152–8.
Lai Z, Moravcova S, Canitrot Y, Andrzejewski LP, Walshe DM, Rea S. Msl2 is a novel component of the vertebrate DNA damage response. PLoS ONE. 2013;8:e68549.
Elkin SR, Lakoduk AM, Schmid SL. Endocytic pathways and endosomal trafficking: a primer. Wien Med Wochenschr. 2016;166:196–204.
Karlsson M, Zhang C, Mear L, Zhong W, Digre A, Katona B, Sjöstedt E, Butler L, Odeberg J, Dusart P, Edfors F, Oksvold P, von Feilitzen K, Zwahlen M, Arif M, Altay O, Li X, Ozcan M, Mardinoglu A, Fagerberg L, Mulder J, Luo Y, Ponten F, Uhlén M, Lindskog C. A single-cell type transcriptomics map of human tissues. Sci Adv. 2021;7:eabh2169.
Lu H, Galeano MCR, Ott E, Kaeslin G, Kausalya PJ, Kramer C, Ortiz-Brüchle N, Hilger N, Metzis V, Hiersche M, Tay SY, Tunningley R, Vij S, Courtney AD, Whittle B, Wühl E, Vester U, Hartleben B, Neuber S, Frank V, Little MH, Epting D, Papathanasiou P, Perkins AC, Wright GD, Hunziker W, Gee HY, Otto EA, Zerres K, Hildebrandt F, Roy S, Wicking C, Bergmann C. Mutations in DZIP1L, which encodes a ciliary-transition-zone protein, cause autosomal recessive polycystic Kidney Disease. Nat Genet. 2017;49:1025–34.
Satir P, Christensen ST. Overview of structure and function of mammalian cilia. Annu Rev Physiol. 2007;69:377–400.
Cadieu E, Neff MW, Quignon P, Walsh K, Chase K, Parker HG, Vonholdt BM, Rhue A, Boyko A, Byers A, Wong A, Mosher DS, Elkahloun AG, Spady TC, André C, Lark KG, Cargill M, Bustamante CD, Wayne RK, Ostrander EA. Coat variation in the domestic dog is governed by variants in three genes. Science. 2009;326:150–3.
Schmutz SM, Berryere TG. Genes affecting coat colour and pattern in domestic dogs: a review. Anim Genet. 2007;38:539–49.
Parker HG, Chase K, Cadieu E, Lark KG, Ostrander EA. An insertion in the RSPO2 gene correlates with improper coat in the Portuguese water dog. J Hered. 2010;101:612–7.
Chintala S, Li W, Lamoreux ML, Ito S, Wakamatsu K, Sviderskaya EV, Bennett DC, Park YM, Gahl WA, Huizing M, Spritz RA, Ben S, Novak EK, Tan J, Swank RT. Slc7a11 gene controls production of pheomelanin pigment and proliferation of cultured cells. Proc Natl Acad Sci U S A. 2005;102:10964–9.
Hayward JJ, Castelhano MG, Oliveira KC, Corey E, Balkman C, Baxter TL, et al. Complex disease and phenotype mapping in the domestic dog. Nat Commun. 2016;7:10460.
Magnusson K, Turkiewicz A, Hughes V, Frobell R, Englund M. High genetic contribution to anterior cruciate ligament rupture: heritability ~ 69. Br J Sports Med. 2021;55:385–9.
Beaulieu ML, Ashton-Miller JA, Wojtys EM. Loading mechanisms of the anterior cruciate ligament. Sports Biomech. 2023;22:1–29.
Acknowledgements
The authors would like to thank the faculty, residents and students throughout the University of Wisconsin-Madison UW Veterinary Care Hospital and many individual breeders and dog owners for their help in recruitment of Labrador Retrievers for this study.
Funding
EEB was supported by the Comparative Biomedical Sciences Training Grant (T32OD010423) funded by the National Institutes of Health. LAB was supported by the National Institutes of Health (T32OD10999) and by the Computation and Informatics in Biology and Medicine Training Program Grant (NLM5T15LM007359) funded by the National Library of Medicine. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author information
Authors and Affiliations
Contributions
LAB and PM conceived and designed the study. BTL, LAB, MM, HT, EEB, SJS, and PM undertook the investigation and generated the data. BTL, MM, and LAB undertook the data analysis and prepared the figures. PM contributed to funding acquisition. PM provided project administration. PM and SJS provided the resources for the study. PM supervised the study. BTL, LAB, and PM wrote the main manuscript text. BTL, LAB, MM, SJS, and PM reviewed and edited the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
All procedures were performed in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health and the American Veterinary Medical Association and with approval from the Animal Care Committee of the University of Wisconsin-Madison (protocols V1070, V5463). Informed consent of each owner was obtained before participation in the study.
Consent for publication
Not applicable.
Competing interests
The authors have declared that no competing interests exist.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
12863_2023_1164_MOESM1_ESM.xlsx
Supplementary File S1. SNP GWAS associations with the multivariate phenotype using the GEMMA algorithm [35].
12863_2023_1164_MOESM2_ESM.xlsx
Supplementary File S2. SNP GWAS associations with ACL rupture and coat color phenotypes using the mvBIMBAM algorithm [36].
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Lee, B., Baker, L., Momen, M. et al. Identification of genetic variants associated with anterior cruciate ligament rupture and AKC standard coat color in the Labrador Retriever. BMC Genom Data 24, 60 (2023). https://doi.org/10.1186/s12863-023-01164-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12863-023-01164-z