A Genetic Predictive Model for Canine Hip Dysplasia: Integration of Genome Wide Association Study (GWAS) and Candidate Gene Approaches

Canine hip dysplasia is one of the most prevalent developmental orthopedic diseases in dogs worldwide. Unfortunately, the success of eradication programs against this disease based on radiographic diagnosis is low. Adding the use of diagnostic genetic tools to the current phenotype-based approach might be beneficial. The aim of this study was to develop a genetic prognostic test for early diagnosis of hip dysplasia in Labrador Retrievers. To develop our DNA test, 775 Labrador Retrievers were recruited. For each dog, a blood sample and a ventrodorsal hip radiograph were taken. Dogs were divided into two groups according to their FCI hip score: control (A/B) and case (D/E). C dogs were not included in the sample. Genetic characterization combining a GWAS and a candidate gene strategy using SNPs allowed a case-control population association study. A mathematical model which included 7 SNPs was developed using logistic regression. The model showed a good accuracy (Area under the ROC curve = 0.85) and was validated in an independent population of 114 dogs. This prognostic genetic test represents a useful tool for choosing the most appropriate therapeutic approach once genetic predisposition to hip dysplasia is known. Therefore, it allows a more individualized management of the disease. It is also applicable during genetic selection processes, since breeders can benefit from the information given by this test as soon as a blood sample can be collected, and act accordingly. In the authors’ opinion, a shift towards genomic screening might importantly contribute to reducing canine hip dysplasia in the future. In conclusion, based on genetic and radiographic information from Labrador Retrievers with hip dysplasia, we developed an accurate predictive genetic test for early diagnosis of hip dysplasia in Labrador Retrievers. However, further research is warranted in order to evaluate the validity of this genetic test in other dog breeds.     

Long-Term Health Effects of Neutering Dogs: Comparison of Labrador Retrievers with Golden Retrievers

Our recent study on the effects of neutering (including spaying) in Golden Retrievers in markedly increasing the incidence of two joint disorders and three cancers prompted this study and a comparison of Golden and Labrador Retrievers. Veterinary hospital records were examined over a 13-year period for the effects of neutering during specified age ranges: before 6 mo., and during 6–11 mo., year 1 or years 2 through 8. The joint disorders examined were hip dysplasia, cranial cruciate ligament tear and elbow dysplasia. The cancers examined were lymphosarcoma, hemangiosarcoma, mast cell tumor, and mammary cancer. The results for the Golden Retriever were similar to the previous study, but there were notable differences between breeds. In Labrador Retrievers, where about 5 percent of gonadally intact males and females had one or more joint disorders, neutering at <6 mo. doubled the incidence of one or more joint disorders in both sexes. In male and female Golden Retrievers, with the same 5 percent rate of joint disorders in intact dogs, neutering at <6 mo. increased the incidence of a joint disorder to 4–5 times that of intact dogs. The incidence of one or more cancers in female Labrador Retrievers increased slightly above the 3 percent level of intact females with neutering. In contrast, in female Golden Retrievers, with the same 3 percent rate of one or more cancers in intact females, neutering at all periods through 8 years of age increased the rate of at least one of the cancers by 3–4 times. In male Golden and Labrador Retrievers neutering had relatively minor effects in increasing the occurrence of cancers. Comparisons of cancers in the two breeds suggest that the occurrence of cancers in female Golden Retrievers is a reflection of particular vulnerability to gonadal hormone removal.     

Genetic evaluation of hip score in UK Labrador Retrievers

Hip dysplasia is an important and complex genetic disease in dogs with both genetic and environmental influences. Since the osteoarthritis that develops is irreversible the only way to improve welfare, through reducing the prevalence, is through genetic selection. This study aimed to evaluate the progress of selection against hip dysplasia, to quantify potential improvements in the response to selection via use of genetic information and increases in selection intensity, and to prepare for public provision of estimated breeding values (EBV) for hip dysplasia in the UK. Data consisted of 25,243 single records of hip scores of Labrador Retrievers between one and four years old, from radiographs evaluated between 2000 and 2007 as part of the British Veterinary Association (BVA) hip score scheme. A natural logarithm transformation was applied to improve normality and linear mixed models were evaluated using ASREML. Genetic correlations between left and right scores, and total hip scores at one, two and three years of age were found to be close to one, endorsing analysis of total hip score in dogs aged one to three as an appropriate approach. A heritability of 0.35±0.016 and small but significant litter effect (0.07±0.009) were estimated. The observed trends in both mean hip score and mean EBV over year of birth indicate that a small genetic improvement has been taking place, approximately equivalent to avoiding those dogs with the worst 15% of scores. Deterministic analysis supported by simulations showed that a 19% greater response could be achieved using EBV compared to phenotype through increases in accuracy alone. This study establishes that consistent but slow genetic improvement in the hip score of UK Labrador Retrievers has been achieved over the previous decade, and demonstrates that progress may be easily enhanced through the use of EBVs and more intense selection.     

Genome Wide Analysis Indicates Genes for Basement Membrane and Cartilage Matrix Proteins as Candidates for Hip Dysplasia in Labrador Retrievers

Hip dysplasia, an abnormal laxity of the hip joint, is seen in humans as well as dogs and is one of the most common skeletal disorders in dogs. Canine hip dysplasia is considered multifactorial and polygenic, and a variety of chromosomal regions have been associated with the disorder. We performed a genome-wide association study in Dutch Labrador Retrievers, comparing data of nearly 18,000 single nucleotide polymorphisms (SNPs) in 48 cases and 30 controls using two different statistical methods. An individual SNP analysis based on comparison of allele frequencies with a χ² statistic was used, as well as a simultaneous SNP analysis based on Bayesian variable selection. Significant association with canine hip dysplasia was observed on chromosome 8, as well as suggestive association on chromosomes 1, 5, 15, 20, 25 and 32. Next-generation DNA sequencing of the exons of genes of seven regions identified multiple associated alleles on chromosome 1, 5, 8, 20, 25 and 32 (p<0.001). Candidate genes located in the associated regions on chromosomes 1, 8 and 25 included LAMA2, LRR1 and COL6A3, respectively. The associated region on CFA20 contained candidate genes GDF15, COMP and CILP2. In conclusion, our study identified candidate genes that might affect susceptibility to canine hip dysplasia. These genes are involved in hypertrophic differentiation of chondrocytes and extracellular matrix integrity of basement membrane and cartilage. The functions of the genes are in agreement with the notion that disruptions in endochondral bone formation in combination with soft tissue defects are involved in the etiology of hip dysplasia.     

SNP-based heritability and genetic architecture of cranial cruciate ligament rupture in Labrador Retrievers

Cranial cruciate ligament rupture (CCLR) is one of the leading causes of pelvic limb lameness in dogs. About 6% of Labrador Retrievers suffer from this orthopedic problem. The aim of this study was to determine the heritability of CCLR in this breed using SNP array genotyping data. DNA samples were collected from CCLR-affected dogs (n = 190) and unaffected dogs over the age of 8 years (n = 143). All 333 dogs were genotyped directly or imputed up to approximately 710k SNPs on the Affymetrix Axiom CanineHD SNP array. Heritability of CCLR was calculated using multiple methodologies, including linear mixed models, Bayesian models and a model that incorporates LD. The covariates of sex and sterilization status were added to each analysis to assess their impact. Across the algorithms of these models, heritability ranged from 0.550 to 0.886, depending on covariate inclusion. The relatively high heritability for this disease indicates that a substantial genetic component contributes to CCLR in the Labrador Retriever.     

Higher medially-directed joint reaction forces are a characteristic of dysplastic hips: A comparative study using subject-specific musculoskeletal models

Acetabular dysplasia is a known cause of hip osteoarthritis. In addition to abnormal anatomy, changes in kinematics, joint reaction forces (JRFs), and muscle forces could cause tissue damage to the cartilage and labrum, and may contribute to pain and fatigue. The objective of this study was to compare lower extremity joint angles, moments, hip JRFs and muscle forces during gait between patients with symptomatic acetabular dysplasia and healthy controls. Marker trajectories and ground reaction forces were measured in 10 dysplasia patients and 10 typically developing control subjects. A musculoskeletal model was scaled in OpenSim to each subject and subject-specific hip joint centers were determined using reconstructions from CT images. Joint kinematics and moments were calculated using inverse kinematics and inverse dynamics, respectively. Muscle forces and hip JRFs were estimated with static optimization. Inter-group differences were tested for statistical significance (p ≤ 0.05) and large effect sizes (d ≥ 0.8). Results demonstrated that dysplasia patients had higher medially directed JRFs. Joint angles and moments were mostly similar between the groups, but large inter-group effect sizes suggested some restriction in range of motion by patients at the hip and ankle. Higher medially-directed JRFs and inter-group differences in hip muscle forces likely stem from lateralization of the hip joint center in dysplastic patients. Joint force differences, combined with reductions in range of motion at the hip and ankle may also indicate compensatory strategies by patients with dysplasia to maintain joint stability.     

Analysis of allele fidelity, polymorphic information content, and density of microsatellites in a genome-wide screening for hip dysplasia in a crossbreed pedigree

Recent advances in genomics resources and tools are facilitating quantitative trait locus mapping. We developed a crossbreed pedigree for mapping quantitative trait loci for hip dysplasia in dogs by crossing dysplastic Labrador Retrievers and normal Greyhounds. We show that one advantage to using a crossbreed pedigree is the increased marker informativeness in the backcross/F2 population relative to the founder populations. We also discuss three factors that affect the detection power in the context of this crossbreed pedigree: being able to detect and correct genotyping errors, increasing marker density for chromosomes with a sparse coverage, and adding individuals to the mapping population as soon as they become available.     

Quantitative Comparison of the Walk and Trot of Border Collies and Labrador Retrievers,Breeds with Different Performance Requirements

Introduction

We hypothesized that breed differences of Border Collies and Labrador Retrievers would be reflected in the temporospatial characteristics of the walk and trot.

Materials and Methods

Twenty healthy Border Collies and 20 healthy Labrador Retrievers made three passes across a pressure sensing walkway system that recorded quantitative temporospatial information at a walk and a trot. The following variables were measured for each dog: velocity, total pressure index percentage (TPI%), ratio of weight borne on the thoracic vs. pelvic limbs (T/P TPI%), stance time percentage (ST%), and thoracic limb stride length (TSrL).

Results

The mean T/P TPI% for Border Collies at a walk and at a trot were significantly lower than for Labrador Retrievers (p = 0.0007 and p = 0.0003). Border Collies had a significantly lower ST% than Labrador Retrievers for the thoracic limbs and pelvic limbs at a walk (p = 0.0058 and 0.0003) and the trot (p = 0.0280 and 0.0448). There was no relationship between ST% and TSrL in Border Collies and an inverse correlation between ST% and TSrL in Labrador Retrievers (p = 0.0002).

Discussion

Key quantitative gait differences were identified in Border Collies and Labrador Retrievers, which could potentially provide each breed with an advantage for their working function.     

Sequence analysis of three canine adipokine genes revealed an association between TNF polymorphisms and obesity in Labrador dogs

Obesity is an emerging health problem in purebred dogs. Due to their crucial role in energy homeostasis control, genes encoding adipokines are considered candidate genes, and their variants may be associated with predisposition to obesity. Searching for polymorphism was carried out in three adipokine genes (TNF, RETN and IL6). The study was performed on 260 dogs, including lean (n = 109), overweight (n = 88) and obese (n = 63) dogs. The largest cohort was represented by Labrador Retrievers (n = 136). Altogether, 24 novel polymorphisms were identified: 12 in TNF (including one missense SNP), eight in RETN (including one missense SNP) and four in IL6. Distributions of five common SNPs (two in TNF, two in RETN and one in IL6) were further analyzed with regard to body condition score. Two SNPs in the non‐coding parts of TNF (c.‐40A>C and c.233+14G>A) were associated with obesity in Labrador dogs. The obtained results showed that the studied adipokine genes are highly polymorphic and two polymorphisms in the TNF gene may be considered as markers predisposing Labrador dogs to obesity.     

Single nucleotide polymorphisms refine QTL intervals for hip joint laxity in dogs

Hip laxity is one characteristic of canine hip dysplasia (CHD), an inheritable disease that leads to hip osteoarthritis. Using a genome-wide screen with 250 microsatellites in a crossbreed pedigree of 159 dysplastic Labrador retrievers and unaffected greyhounds, we previously identified putative ( P  <   0.01) QTL on canine chromosomes 11 and 29 (CFA11 and CFA29). To refine these QTL locations, we have genotyped 257 dogs including 105 Labrador retrievers, seven greyhounds, four generations of their crossbreed offspring and three German shepherds for 111 and 171 SNPs on CFA11 and CFA29 respectively. The distraction index (DI, a measure of maximum hip laxity) was used as an intermediate phenotype that predicts whether a hip joint will or will not develop osteoarthritis. Using a multipoint linkage analysis, significant evidence (95% posterior probability) was found for QTL contributing to hip laxity in the 16.2–21 cM region on CFA11 that explained 15–18% of the total variance in DI. Evidence for an independent QTL on CFA29 was weaker than that on CFA11. Identification of the causative mutation(s) will lead to better understanding of biochemical pathways in both dogs and humans with hip laxity and dysplasia.     

Factors affecting gestation duration in the bitch

A retrospective analysis was performed to determine the effects of age, breed, parity, and litter size on the duration of gestation in the bitch. Bitches at two locations were monitored from breeding to whelping. A total of 764 litters whelped from 308 bitches (36 large hounds, 34 Golden Retrievers, 23 German Shepherd Dogs (GSD), and 215 Labrador Retrievers). By breed, the number of whelpings was 152, 72, 58, and 482 for the hounds, Golden Retrievers, German Shepherd Dogs, and Labrador Retrievers, respectively. Whelping was predicted to be 57 d from the first day of cytologic diestrus in the hounds or 65 d from the initial progesterone rise in the other breeds. The average gestation duration (calculated as 8 d prior to Day 1 of cytologic diestrus in hounds or measured from the initial progesterone rise in other breeds) by breed (days +/- S.D.) was 66.0 +/- 2.8, 64.7 +/- 1.5, 63.6 +/- 2.1, and 62.9 +/- 1.3 for the hounds, Golden Retrievers, German Shepherd Dogs, and Labrador Retrievers, respectively. The relationship of age, breed, parity, and litter size with the difference in gestation duration was evaluated using log linear modeling. Age or parity had no effect on gestation duration. Compared to Labrador Retrievers, the German Shepherd Dogs, Golden Retrievers and hounds were more likely to have a longer gestation duration; three, four and nearly eight times as likely, respectively. Bitches whelping four or fewer pups were significantly more likely to have a longer gestation duration than those whelping five or more pups; the prolongation averaging 1 d.     

Hip biomechanics in orthopaedic clinical practice

Radiographic and clinical studies, coupled with biomechanical assessment of the hip, are important tools for predicting the development of osteoarthitis of the hip. In order to better understand the treatment of hip dysplasia, it is necessary to determine the contact stress in the hip joint. In this study, a three-dimensional mathematical model was used to determine hip joint contact stress. Because of the discrepancy in the results of analyses of different radiographic indicators of hip dysplasia, the calculation of hip joint contact stress is proposed for a more accurate assessment of the severity of hip dysplasia.     

Retrospective Analysis for Genetic Improvement of Hip Joints of Cohort Labrador Retrievers in the United States: 1970–2007

Background

Canine Hip Dysplasia (CHD) is a common inherited disease that affects dog wellbeing and causes a heavy financial and emotional burden to dog owners and breeders due to secondary hip osteoarthritis. The Orthopedic Foundation for Animals (OFA) initiated a program in the 1960''s to radiograph hip and elbow joints and release the OFA scores to the public for breeding dogs against CHD. Over last four decades, more than one million radiographic scores have been released.

Methodology/Principal Findings

The pedigrees in the OFA database consisted of 258,851 Labrador retrievers, the major breed scored by the OFA (25% of total records). Of these, 154,352 dogs had an OFA hip score reported between 1970 and 2007. The rest of the dogs (104,499) were the ancestors of the 154,352 dogs to link the pedigree relationships. The OFA hip score is based on a 7-point scale with the best ranked as 1 (excellent) and the worst hip dysplasia as 7. A mixed linear model was used to estimate the effects of age, sex, and test year period and to predict the breeding value for each dog. Additive genetic and residual variances were estimated using the average information restricted maximum likelihood procedure. The analysis also provided an inbreeding coefficient for each dog. The hip scores averaged 1.93 (±SD = 0.59) and the heritability was 0.21. A steady genetic improvement has accrued over the four decades. The breeding values decreased (improved) linearly. By the end of 2005, the total genetic improvement was 0.1 units, which is equivalent to 17% of the total phenotypic standard deviation.

Conclusion/Significance

A steady genetic improvement has been achieved through the selection based on the raw phenotype released by the OFA. As the heritability of the hip score was on the low end (0.21) of reported ranges, we propose that selection based on breeding values will result in more rapid genetic improvement than breeding based on phenotypic selection alone.     

BMI and Severity of Clinical and Radiographic Signs of Hip Osteoarthritis

Obesity might be involved in the pathogenesis of osteoarthritis (OA) not only via increased mechanical loading, but also via an inflammatory component possibly causing increased pain and functional disability. The study aim was to examine the relationship between BMI and clinical symptoms as well as radiographic severity of OA in patients scheduled for primary total hip arthroplasty (THA). We conducted a cross‐sectional study of 855 patients scheduled for a first THA for primary OA at a single centre between November 2001 and December 2006. The primary outcome was clinical and radiographic severity of OA, which was evaluated in four BMI categories (18.5–24.9, 25–29.9, 30–34.9 and ≥35 kg/m²). We used the Harris Hip Score (HHS) and the Western Ontario McMaster Universities Osteoarthritis Index (WOMAC) to assess pain and function. The severity of radiographic hip joint damage was evaluated using the Kellgren–Lawrence classification. Multivariate analyses were performed to adjust for potential confounders. In patients scheduled for THA, increasing BMI was associated with significantly higher levels of pain and functional disability on both HHS (P for trend <0.001) and WOMAC (P for trend <0.001). However, the degree of radiographic joint damage remained similar across BMI categories. These findings emphasize the need to further investigate the potential pathogenic role of obesity and low‐grade inflammation in OA and underscore the importance of obesity prevention to avoid early prosthetic replacement of the hip.     

A chromosomal duplication that includes the canine microsatellite INRA21 in Labrador Retrievers

INRA21 is one of the canine microsatellites recommended for parentage verification by the International Society for Animal Genetics. In Labrador Retrievers, abnormal peak patterns such as three-peak patterns during capillary electrophoresis were frequently observed at INRA21. Pedigree analysis indicated that the abnormal peak patterns were due to inheritable causes, and semiquantitative multiplex (SQM) PCR analysis showed that the abnormal peak patterns were caused by chromosomal duplication. Walking SQM-PCR analysis revealed that the size of the duplicated segment was approximately 1.58 Mb. Genotypes of microsatellites within the duplicated segment indicated that the duplication was an identical-by-descent mutation. This duplication is probably carried by more than half of the dogs in the Japanese population of Labrador Retrievers. The abnormal peak patterns at INRA21 were also observed in German Shorthaired Pointers and Flat-Coated Retrievers. Genotyping analysis of the microsatellites within the duplicated segment in Labrador Retrievers suggested that the abnormal peak patterns observed in the two breeds were due to the duplication inherited from the same ancestor as the duplication of Labrador Retrievers. This study urges attention to the use of INRA21 and shows an example of copy number polymorphisms that are characteristic to dog breeds or lineages.     

Quantitative trait loci for hip dysplasia in a crossbreed canine pedigree

Canine hip dysplasia is a common developmental inherited trait characterized by hip laxity, subluxation or incongruity of the femoral head and acetabulum in affected hips. The inheritance pattern is complex and the mutations contributing to trait expression are unknown. In the study reported here, 240 microsatellite markers distributed in 38 autosomes and the X chromosome were genotyped on 152 dogs from three generations of a crossbred pedigree based on trait-free Greyhound and dysplastic Labrador Retriever founders. Interval mapping was undertaken to map the QTL underlying the quantitative dysplastic traits of maximum passive hip laxity (the distraction index), the dorsolateral subluxation score, and the Norberg angle. Permutation testing was used to derive the chromosome-wide level of significance at p < 0.05 for each QTL. Chromosomes 4, 9, 10, 11 (p < 0.01), 16, 20, 22, 25, 29 (p < 0.01), 30, 35, and 37 harbor putative QTL for one or more traits. Successful detection of QTL was due to the crossbreed pedigree, multiple-trait measurements, control of environmental background, and marked advancement in canine mapping tools.     

Pre-treatment systemic immune-inflammation index is a useful prognostic indicator in patients with breast cancer undergoing neoadjuvant chemotherapy

The systemic immune-inflammation index (SII = N × P/L) based on neutrophil (N), platelet (P) and lymphocyte (L) counts is used to predict the survival of patients with malignant tumours and can fully reflect the balance between host inflammatory and immune status. This study is conducted to explore the potential prognostic significance of SII in patients with breast cancer undergoing neoadjuvant chemotherapy (NACT). A total of 262 patients with breast cancer received NACT were enrolled in this study. According to the receiver operating characteristic curve, the optimal cut-off value of SII was divided into two groups: low SII group (<602 × 10⁹/L) and high SII group (≥602 × 10⁹/L). The associations between breast cancer and clinicopathological variables by SII were determined by chi-squared test or Fisher's exact test. The Kaplan-Meier plots and log-rank test were used to determine clinical outcomes of disease-free survival (DFS) and overall survival (OS). The prognostic value of SII was analysed by univariate and multivariate Cox proportional hazards regression models. The toxicity of NACT was accessed by National Cancer Institute Common Toxicity Criteria (NCICTC). According to univariate and multivariate Cox regression survival analyses, the results showed that the value of SII had prognostic significance for DFS and OS. The patients with low SII value had longer DFS and OS than those with high SII value (31.11 vs 40.76 months, HR: 1.075, 95% CI: 0.718-1.610, P = .006; 44.47 vs 53.68 months, HR: 1.051, 95% CI: 0.707-1.564, P = .005, respectively). The incidence of DFS and OS in breast cancer patients with low SII value was higher than that in those patients with high SII value in 3-, 5- and 10-year rates. The common toxicities after NACT were haematological and gastrointestinal reaction, and there were no differences by SII for the assessment of side effects of neoadjuvant chemotherapy. Meanwhile, the results also proved that breast cancer patients with low SII value and high Miller and Payne grade (MPG) survived longer than those breast cancer with high SII value and low MPG grade. In patients without lymph vessel invasion, these breast cancer patients with low SII value had better prognosis and lower recurrence rates than those with high SII value. Pre-treatment SII with the advantage of reproducible, convenient and non-invasive was a useful prognostic indicator for breast cancer patients undergoing neoadjuvant chemotherapy and is a promising biomarker for breast cancer on treatment strategy decisions.     

The association between leptin,interleukin-6, and hip radiographic osteoarthritis in older people: a cross-sectional study

Introduction  

The associations between leptin, interleukin (IL)-6, and hip radiographic osteoarthritis (OA) have not been reported, and their roles in obesity-related hip OA are unclear. The aim of this study was to describe the associations between leptin, IL-6, and hip radiographic osteoarthritis (ROA) in older adults.     

Genome-wide association analysis for lethal brachycephalic-like facial dysmorphia in Labrador Retrievers

A GWAS was performed for inborn X-linked facial dysmorphia with severe growth retardation in Labrador Retrievers. This lethal condition was mapped on the X chromosome at 17–21 Mb and supported by eight SNPs in complete LD. Dams of affected male puppies were heterozygous for the significantly associated SNPs and male affected puppies carried the associated alleles hemizygously. In the near vicinity to the associated region, RPS6KA3 was identified as a candidate gene causing facial dysmorphia in humans and mice known as Coffin–Lowry syndrome. Haplotype analysis showed significant association with the phenotypes of all 18 animals under study. This haplotype was validated through normal male progeny from a dam with the not-associated haplotype on both X chromosomes but male affected full-sibs with the associated haplotype.