Validation of DNA test for hip dysplasia failed in Danish Labrador Retrievers

Canine hip dysplasia is characterized by poor hip joint conformation and laxity. The disease is a complex trait influenced by both genetics and environment. Diagnosis and quantification of hip dysplasia are performed by radiographic examination of the hip joint and the diagnosis is used for making breeding decisions in many breeds. A prognostic genetic test (the Dysgen test) based on seven associated SNPs has been developed in a study based on Spanish Labrador Retrievers. In our study this test has been evaluated in 39 Danish Labrador Retrievers with known radiographic hip score: 14 with hip dysplasia (grade D or E) and 25 without hip dysplasia (grade A or B). There was no significant correlation between the Dysgen test results and the radiographic hip status (P = 0.3203) in these dogs, indicating that Dysgen test results obtained for Danish Labrador Retrievers have no prognostic value.     

Accuracy of genotype imputation in Labrador Retrievers

The dog is a valuable model species for the genetic analysis of complex traits, and the use of genotype imputation in dogs will be an important tool for future studies. It is of particular interest to analyse the effect of factors like single nucleotide polymorphism (SNP) density of genotyping arrays and relatedness between dogs on imputation accuracy due to the acknowledged genetic and pedigree structure of dog breeds. In this study, we simulated different genotyping strategies based on data from 1179 Labrador Retriever dogs. The study involved 5826 SNPs on chromosome 1 representing the high density (HighD) array; the low‐density (LowD) array was simulated by masking different proportions of SNPs on the HighD array. The correlations between true and imputed genotypes for a realistic masking level of 87.5% ranged from 0.92 to 0.97, depending on the scenario used. A correlation of 0.92 was found for a likely scenario (10% of dogs genotyped using HighD, 87.5% of HighD SNPs masked in the LowD array), which indicates that genotype imputation in Labrador Retrievers can be a valuable tool to reduce experimental costs while increasing sample size. Furthermore, we show that genotype imputation can be performed successfully even without pedigree information and with low relatedness between dogs in the reference and validation sets. Based on these results, the impact of genotype imputation was evaluated in a genome‐wide association analysis and genomic prediction in Labrador Retrievers.     

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.     

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.     

Differential genetic regulation of canine hip dysplasia and osteoarthritis

Background

Canine hip dysplasia (HD) is a common polygenic trait characterized by hip malformation that results in osteoarthritis (OA). The condition in dogs is very similar to developmental dysplasia of the human hip which also leads to OA.

Methodology/Principal Findings

A total of 721 dogs, including both an association and linkage population, were genotyped. The association population included 8 pure breeds (Labrador retriever, Greyhounds, German Shepherd, Newfoundland, Golden retriever, Rottweiler, Border Collie and Bernese Mountain Dog). The linkage population included Labrador retrievers, Greyhounds, and their crosses. Of these, 366 dogs were genotyped at ∼22,000 single nucleotide polymorphism (SNP) loci and a targeted screen across 8 chromosomes with ∼3,300 SNPs was performed on 551 dogs (196 dogs were common to both sets). A mixed linear model approach was used to perform an association study on this combined association and linkage population. The study identified 4 susceptibility SNPs associated with HD and 2 SNPs associated with hip OA.

Conclusion/Significance

The identified SNPs included those near known genes (PTPRD, PARD3B, and COL15A1) reported to be associated with, or expressed in, OA in humans. This suggested that the canine model could provide a unique opportunity to identify genes underlying natural HD and hip OA, which are common and debilitating conditions in both dogs and humans.     

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.     

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.     

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.     

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.     

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.     

Multiple loci associated with canine hip dysplasia (CHD) in German shepherd dogs

Canine hip dysplasia (CHD) is the most common hereditary skeletal disorder in dogs. To identify common alleles associated with CHD, we developed 37 informative single nucleotide polymorphisms (SNPs) within 13 quantitative trait loci (QTL) previously identified for German shepherd dogs. These SNPs were genotyped in 95 German shepherd dogs affected by CHD and 95 breed, sex, and birth year-matched controls. A total of ten SNPs significant at a nominal P value of 0.05 were validated in 843 German shepherd dogs including 277 unaffected dogs and 566 CHD-affected dogs. Cases and controls were sampled from the whole German shepherd dog population in Germany in such a way that mean coancestry coefficients were below 0.1 % within cases and controls as well as among cases and controls. We identified nine SNPs significantly associated with CHD within five QTL on dog chromosomes (CFA) 3, 9, 26, 33, and 34. Genotype effects of these nine SNPs explained between 22 and 34 % of the phenotypic variance of hip dysplasia in German shepherd dogs. The strongest associated SNPs were located on CFA33 and 34 within the candidate genes PNCP, TRIO, and SLC6A3. Thus, the present study validated positional candidate genes within five QTL for CHD.     

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.     

Toward the most ideal case–control design with related and unrelated dogs in whole‐exome sequencing studies

With the recent development of whole‐exome sequencing enrichment designs for the dog, a novel tool for disease‐association studies became available. The aim of disease‐association studies is to identify one or a very limited number of putative causal variants or genes from the large pool of genetic variation. To maximize the efficiency of these studies and to provide some directions of what to expect, we evaluated the effect on variant reduction for various combinations of cases and controls for both dominant and recessive types of inheritance assuming variable degrees of penetrance and detectance. In this study, variant data of 14 dogs (13 Labrador Retrievers and one Dogue de Bordeaux), obtained by whole‐exome sequencing, were analyzed. In the filtering process, we found that unrelated dogs from the same breed share up to 70% of their variants, which is likely a consequence of the breeding history of the dog. For the designs tested with unrelated dogs, combining two cases and two controls gave the best result. These results were improved further by adding closely related dogs. Reduced penetrance and/or detectance has a drastic effect on the efficiency and is likely to have a profound effect on the sample size needed to elucidate the causal variant. Overall, we demonstrated that sequencing a small number of dogs results in a marked reduction of variants that are likely sufficient to pinpoint causal variants or genes.     

Polymorphism and methylation of the <Emphasis Type="Italic">MC4R</Emphasis> gene in obese and non-obese dogs

The dog is considered to be a useful biomedical model for human diseases and disorders, including obesity. One of the numerous genes associated with human polygenic obesity is MC4R, encoding the melanocortin 4 receptor. The aim of our study was to analyze polymorphisms and methylation of the canine MC4R in relation to adiposity. Altogether 270 dogs representing four breeds predisposed to obesity: Labrador Retriever (n?=?187), Golden Retriever (n?=?38), Beagle (n?=?28) and Cocker Spaniel (n?=?17), were studied. The dogs were classified into three groups: lean, overweight and obese, according to the 5-point Body Condition Score (BCS) scale. In the cohort of Labradors a complete phenotypic data (age, sex, neutering status, body weight and BCS) were collected for 127 dogs. The entire coding sequence as well as 5′ and 3′-flanking regions of the studied gene were sequenced and six polymorphic sites were reported. Genotype frequencies differed considerably between breeds and Labrador Retrievers appeared to be the less polymorphic. Moreover, distribution of some polymorphic variants differed significantly (P?<?0.05) between small cohorts with diverse BCS in Golden Retrievers (c.777T>C, c.868C>T and c.*33C>G) and Beagles (c.-435T>C and c.637G>T). On the contrary, in Labradors no association between the studied polymorphisms and BCS or body weight was observed. Methylation analysis, using bisulfite DNA conversion followed by Sanger sequencing, was carried out for 12 dogs with BCS?=?3 and 12 dogs with BCS?=?5. Two intragenic CpG islands, containing 19 cytosines, were analyzed and the methylation profile did not differ significantly between lean and obese animals. We conclude that an association of the MC4R gene polymorphism with dog obesity or body weight is unlikely, in spite of the fact that some associations were found in small cohorts of Beagles and Golden Retrievers. Also methylation level of this gene is not related with dog adiposity.     

A splice site variant in the SUV39H2 gene in Greyhounds with nasal parakeratosis

Hereditary nasal parakeratosis (HNPK), described in the Labrador Retriever breed, is a monogenic autosomal recessive disorder that causes crusts and fissures on the nasal planum of otherwise healthy dogs. Our group previously showed that this genodermatosis may be caused by a missense variant located in the SUV39H2 gene encoding a histone 3 lysine 9 methyltransferase, a chromatin modifying enzyme with a potential role in keratinocyte differentiation. In the present study, we investigated a litter of Greyhounds in which six out of eight puppies were affected with parakeratotic lesions restricted to the nasal planum. Clinically and histologically, the lesions were comparable to HNPK in Labrador Retrievers. Whole genome sequencing of one affected Greyhound revealed a 4‐bp deletion at the 5′‐end of intron 4 of the SUV39H2 gene that was absent in 188 control dog and three wolf genomes. The variant was predicted to disrupt the 5′‐splice site with subsequent loss of SUV39H2 function. The six affected puppies were homozygous for the variant, whereas the two non‐affected littermates were heterozygous. Genotyping of a larger cohort of Greyhounds revealed that the variant is segregating in the breed and that this breed might benefit from genetic testing to avoid carrier × carrier matings.     

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.     

A COL11A2 Mutation in Labrador Retrievers with Mild Disproportionate Dwarfism

We describe a mild form of disproportionate dwarfism in Labrador Retrievers, which is not associated with any obvious health problems such as secondary arthrosis. We designate this phenotype as skeletal dysplasia 2 (SD2). It is inherited as a monogenic autosomal recessive trait with incomplete penetrance primarily in working lines of the Labrador Retriever breed. Using 23 cases and 37 controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 4.44 Mb interval on chromosome 12. We re-sequenced the genome of one affected dog at 30x coverage and detected 92 non-synonymous variants in the critical interval. Only two of these variants, located in the lymphotoxin A (LTA) and collagen alpha-2(XI) chain gene (COL11A2), respectively, were perfectly associated with the trait. Previously described COL11A2 variants in humans or mice lead to skeletal dysplasias and/or deafness. The dog variant associated with disproportionate dwarfism, COL11A2:c.143G>C or p.R48P, probably has only a minor effect on collagen XI function, which might explain the comparatively mild phenotype seen in our study. The identification of this candidate causative mutation thus widens the known phenotypic spectrum of COL11A2 mutations. We speculate that non-pathogenic COL11A2 variants might even contribute to the heritable variation in height.     

Identification and Validation of Quantitative Trait Loci (QTL) for Canine Hip Dysplasia (CHD) in German Shepherd Dogs

Canine hip dysplasia (CHD) is the most common hereditary skeletal disorder in dogs. To identify common alleles associated with CHD, we genotyped 96 German Shepherd Dogs affected by mild, moderate and severe CHD and 96 breed, sex, age and birth year matched controls using the Affymetrix canine high density SNP chip. A mixed linear model analysis identified five SNPs associated with CHD scores on dog chromosomes (CFA) 19, 24, 26 and 34. These five SNPs were validated in a by sex, age, birth year and coancestry stratified sample of 843 German Shepherd Dogs including 277 unaffected dogs and 566 CHD-affected dogs. Mean coancestry coefficients among and within cases and controls were <0.1%. Genotype effects of these SNPs explained 20–32% of the phenotypic variance of CHD in German Shepherd Dogs employed for validation. Genome-wide significance in the validation data set could be shown for each one CHD-associated SNP on CFA24, 26 and 34. These SNPs are located within or in close proximity of genes involved in bone formation and related through a joint network. The present study validated positional candidate genes within two previously known quantitative trait loci (QTL) and a novel QTL for CHD in German Shepherd Dogs.     

An approach to canine behavioural genetics employing guide dogs for the blind

The purpose of this study was to attempt to find related variables of the canine genome with behavioural traits of dogs maintained and tested in a guide dog facility which provided a relatively uniform environment. The study involved 81 Labrador Retrievers that were being trained as guide dogs. Each dog was taken on walk-out sessions in which the trainer weekly recorded observations that were related to behavioural traits. The records were subjected to key-word analysis of 14 behaviour-related words. A factor analysis on the appearance rate of the 14 key words or phrases resulted in the extraction of six factors that accounted for 67.4% of the variance. Factor 1, referred to as aggressiveness, was significantly related to the success or failure of the dog in qualifying as a guide dog, and was also related to the variable of litter identification. Factor 2, referred to as distraction, was related to the variable of trainer. Factor 3, activity level, was related to the variable of sex, and was significantly related to the polymorphisms of c.471T>C in the solute carrier family 1 ( neuronal/epithelial high affinity glutamate transporter ) member 2 gene and c.216G>A in the catechol-O-methyltransferase gene. The involvement of polymorphisms c.471T>C and c.216G>A in behavioural patterns related to activity level is similar to comparable genetic studies in other mammalian species. These results contribute to a greater understanding of the role of these genes in behaviour.