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.     

Identification of Quantitative Trait Loci (QTL) for Canine Hip Dysplasia and Canine Elbow Dysplasia in Bernese Mountain Dogs

A genome-wide association study for canine hip dysplasia (CHD) and canine elbow dysplasia (CED) using the Illumina canine high density bead chip had been performed for 174 Bernese mountain dogs. General and mixed linear model analysis identified two different regions with single nucleotide polymorphisms (SNPs) on dog chromosome (CFA) 14 significantly associated with CHD and a further significantly CHD-associated region on CFA37. For CED, four SNPs on CFA11 and 27 were significantly associated. The identified SNPs of four associated regions included nearby candidate genes. These possible positional candidates were the genes PON2 on CFA14 and FN1 on CFA37 for CHD and the genes LMNB1 on CFA11 and WNT10B on CFA27 for CED.     

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.     

Mapping quantitative trait loci for canine hip dysplasia in German Shepherd dogs

Canine hip dysplasia (CHD) is a common hereditary developmental disease of the coxofemoral joints. CHD is characterized by subluxation of the femoral head and deformation of the acetabulum leading to a painful osteoarthrosis. Analyses of mode of inheritance have shown the involvement of a major gene in expression of CHD in German Shepherd dogs. Thus, a whole genome scan for quantitative trait loci (QTL) was performed in German Shepherd dogs. For this purpose 11 paternal half-sib families, including a total of 459 purebred German Shepherd dogs with sires, dams, and offspring, were genotyped for 261 microsatellites. These markers were equidistantly distributed over all 38 autosomes and the X chromosome with an average marker distance of 11.7 cM. The mean observed heterozygosity of the marker set was 50%. The CHD status for the dogs was scored according to the official rules of the Fédération Cynologique Internationale. At the genome-wide level of significance at p < 0.05, QTL for CHD were located on nine different canine chromosomes: 1, 3, 4, 8, 9, 16, 19, 26, and 33. The minimal QTL regions containing the CHD genes spanned on average 5 Mb with a range between 1 and 8.2 Mb. Chromosome-wide level of significance at p < 0.05 was found for QTL on 19 chromosomes. Further analyses can now be performed to refine these map positions of QTL already identified in German Shepherd dogs. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

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.     

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.     

Complex segregation analysis of canine hip dysplasia in German shepherd dogs

Complex segregation analyses were carried out to clarify the mode of inheritance of canine hip dysplasia (CHD) in German shepherd dogs. Data were used from 8,567 animals examined for CHD from 20 families with three to four generations. The existence of a major gene in addition to polygenic gene effects was detected. In the present study, a mixed model with a dominant major gene effect seemed to be most probable for dichotomous encoding (0: dogs without signs of CHD; 1: dogs with borderline/slight to severe CHD). In addition, mixed major gene inheritance was shown for a binary trait where borderline was assigned to dogs scored free from CHD and for a trichotomously encoded trait (0: dogs without signs of CHD; 1: borderline CHD; 2: mild to severe CHD). Although only small frequencies were found for the unfavorable homozygotic genotype AA, the probability of the AB genotype was high in affected animals. Selection schemes to reduce the frequency of the allele A should therefore efficiently improve existing breeding programmes in German shepherd dogs.     

Quantitative Trait Loci (QTL) Study Identifies Novel Genomic Regions Associated to Chiari-Like Malformation in Griffon Bruxellois Dogs

Chiari-like malformation (CM) is a developmental abnormality of the craniocervical junction that is common in the Griffon Bruxellois (GB) breed with an estimated prevalence of 65%. This disease is characterized by overcrowding of the neural parenchyma at the craniocervical junction and disturbance of cerebrospinal fluid (CSF) flow. The most common clinical sign is pain either as a direct consequence of CM or neuropathic pain as a consequence of secondary syringomyelia. The etiology of CM remains unknown but genetic factors play an important role. To investigate the genetic complexity of the disease, a quantitative trait locus (QTL) approach was adopted. A total of 14 quantitative skull and atlas measurements were taken and were tested for association to CM. Six traits were found to be associated to CM and were subjected to a whole-genome association study using the Illumina canine high density bead chip in 74 GB dogs (50 affected and 24 controls). Linear and mixed regression analyses identified associated single nucleotide polymorphisms (SNPs) on 5 Canis Familiaris Autosomes (CFAs): CFA2, CFA9, CFA12, CFA14 and CFA24. A reconstructed haplotype of 0.53 Mb on CFA2 strongly associated to the height of the cranial fossa (diameter F) and an haplotype of 2.5 Mb on CFA14 associated to both the height of the rostral part of the caudal cranial fossa (AE) and the height of the brain (FG) were significantly associated to CM after 10 000 permutations strengthening their candidacy for this disease (P = 0.0421, P = 0.0094 respectively). The CFA2 QTL harbours the Sall-1 gene which is an excellent candidate since its orthologue in humans is mutated in Townes-Brocks syndrome which has previously been associated to Chiari malformation I. Our study demonstrates the implication of multiple traits in the etiology of CM and has successfully identified two new QTL associated to CM and a potential candidate gene.     

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.     

Genetics of canine anal furunculosis in the German shepherd dog

Canine anal furunculosis (AF) is characterised by ulceration and fistulation of perianal tissue and is a disease that particularly affects German shepherd dogs (GSDs). There are some similarities between AF and perianal Crohn’s disease (CD) in man. An immune-mediated aetiopathogenesis for AF has been suggested due to tissue pathology, a major histocompatibility complex (MHC) association and clinical response to ciclosporin. Genome-wide association studies (GWAS) can be conducted in dogs with fewer markers and individuals than would be required in a human study. A discovery GWAS was performed on 21 affected and 25 control GSDs from the UK. No SNPs reached genome-wide significance levels at this stage. However, 127 nominally associated SNPs were genotyped in further 76 cases and 191 controls from the UK and Finland. Sequencing of these regions was undertaken to discover novel genetic variation. Association testing of these variants in the UK and Finnish cohorts revealed nine significantly associated SNPs, six of which cause non-synonymous changes in protein sequence. The ADAMTS16 and CTNND2 gene regions were most significantly associated with disease. Members of the butyrophilin protein family, important in intestinal inflammatory regulation, were also associated with disease, but their independence from the MHC region remains to be established. The CTNND2 gene region is also interesting as this locus was implicated in human ulcerative colitis and CD, albeit at a different candidate gene: DAP. We suggest that this represents a common association between inflammatory bowel disease-related conditions in both species and believe that future studies will strengthen this link.     

Congenital Sensorineural Deafness in Dalmatian Dogs Associated with Quantitative Trait Loci

A genome-wide association study (GWAS) was performed for 235 Dalmatian dogs using the canine Illumina high density bead chip to identify quantitative trait loci (QTL) associated with canine congenital sensorineural deafness (CCSD). Data analysis was performed for all Dalmatian dogs and in addition, separately for brown-eyed and blue-eyed dogs because of the significant influence of eye colour on CCSD in Dalmatian dogs. Mixed linear model analysis (MLM) revealed seven QTL with experiment-wide significant associations (-log₁₀P>5.0) for CCSD in all Dalmatian dogs. Six QTL with experiment-wide significant associations for CCSD were found in brown-eyed Dalmatian dogs and in blue-eyed Dalmatian dogs, four experiment-wide significant QTL were detected. The experiment-wide CCSD-associated SNPs explained 82% of the phenotypic variance of CCSD. Five CCSD-loci on dog chromosomes (CFA) 6, 14, 27, 29 and 31 were in close vicinity of genes shown as causative for hearing loss in human and/or mouse.     

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.     

The Shepherds’ Tale: A Genome-Wide Study across 9 Dog Breeds Implicates Two Loci in the Regulation of Fructosamine Serum Concentration in Belgian Shepherds

Diabetes mellitus is a serious health problem in both dogs and humans. Certain dog breeds show high prevalence of the disease, whereas other breeds are at low risk. Fructosamine and glycated haemoglobin (HbA1c) are two major biomarkers of glycaemia, where serum concentrations reflect glucose turnover over the past few weeks to months. In this study, we searched for genetic factors influencing variation in serum fructosamine concentration in healthy dogs using data from nine dog breeds. Considering all breeds together, we did not find any genome-wide significant associations to fructosamine serum concentration. However, by performing breed-specific analyses we revealed an association on chromosome 3 (p_corrected ≈ 1:68 × 10^-6) in Belgian shepherd dogs of the Malinois subtype. The associated region and its close neighbourhood harbours interesting candidate genes such as LETM1 and GAPDH that are important in glucose metabolism and have previously been implicated in the aetiology of diabetes mellitus. To further explore the genetics of this breed specificity, we screened the genome for reduced heterozygosity stretches private to the Belgian shepherd breed. This revealed a region with reduced heterozygosity that shows a statistically significant interaction (p = 0.025) with the association region on chromosome 3. This region also harbours some interesting candidate genes and regulatory regions but the exact mechanisms underlying the interaction are still unknown. Nevertheless, this finding provides a plausible explanation for breed-specific genetic effects for complex traits in dogs. Shepherd breeds are at low risk of developing diabetes mellitus. The findings in Belgian shepherds could be connected to a protective mechanism against the disease. Further insight into the regulation of glucose metabolism could improve diagnostic and therapeutic methods for diabetes mellitus.     

Genetic variants at newly identified lipid loci are associated with coronary heart disease in a Chinese Han population

Background

Recent genome-wide association studies (GWAS) have mapped several novel loci influencing blood lipid levels in Caucasians. We sought to explore whether the genetic variants at newly identified lipid-associated loci were associated with CHD susceptibility in a Chinese Han population.

Methodology/Principal Findings

We conducted a two-stage case-control study in a Chinese Han population. The first-stage, consisting of 1,376 CHD cases and 1,376 sex and age- frequency matched controls, examined 5 novel lipid-associated single-nucleotide polymorphisms (SNPs) identified from GWAS among Caucasians in relation to CHD risk in Chinese. We then validated significant SNPs in the second-stage, consisting of 1,269 cases and 2,745 controls. We also tested associations between SNPs within the five novel loci and blood lipid levels in 4,121 controls. We identified two novel SNPs (rs599839 in CELSR2-PSRC1-SORT1 and rs16996148 in NCAN-CILP2) that were significantly associated with reduced CHD risk in Chinese (odds ratios (95% confidence intervals) in the dominant model 0.76 (0.61-0.90; P = 0.001), 0.67 (0.57-0.77; P = 3.4×10⁻⁸), respectively). Multiple linear regression analyses using dominant model showed that rs599839 was significantly associated with decreased LDL levels (P = 0.022) and rs16996148 was significantly associated with increased LDL and HDL levels (P = 2.9×10⁻⁴ and 0.001, respectively).

Conclusions/Significance

We identified two novel SNPs (rs599839 and rs16996148) at newly identified lipid-associated loci that were significantly associated with CHD susceptibility in a Chinese Han population.     

Utilizing the Dog Genome in the Search for Novel Candidate Genes Involved in Glioma Development—Genome Wide Association Mapping followed by Targeted Massive Parallel Sequencing Identifies a Strongly Associated Locus

Gliomas are the most common form of malignant primary brain tumors in humans and second most common in dogs, occurring with similar frequencies in both species. Dogs are valuable spontaneous models of human complex diseases including cancers and may provide insight into disease susceptibility and oncogenesis. Several brachycephalic breeds such as Boxer, Bulldog and Boston Terrier have an elevated risk of developing glioma, but others, including Pug and Pekingese, are not at higher risk. To identify glioma-associated genetic susceptibility factors, an across-breed genome-wide association study (GWAS) was performed on 39 dog glioma cases and 141 controls from 25 dog breeds, identifying a genome-wide significant locus on canine chromosome (CFA) 26 (p = 2.8 x 10⁻⁸). Targeted re-sequencing of the 3.4 Mb candidate region was performed, followed by genotyping of the 56 SNVs that best fit the association pattern between the re-sequenced cases and controls. We identified three candidate genes that were highly associated with glioma susceptibility: CAMKK2, P2RX7 and DENR. CAMKK2 showed reduced expression in both canine and human brain tumors, and a non-synonymous variant in P2RX7, previously demonstrated to have a 50% decrease in receptor function, was also associated with disease. Thus, one or more of these genes appear to affect glioma susceptibility.     

The Association Study between Twenty One Polymorphisms in Seven Candidate Genes and Coronary Heart Diseases in Chinese Han Population

Previous genome-wide association studies (GWAS) in multiple populations identified several genetic loci for coronary heart diseases (CHD). Here we utilized a 2-stage candidate gene association strategy in Chinese Han population to shed light on the putative association between several metabolic-related candidate genes and CHD. At the 1^st stage, 190 patients with CHD and 190 controls were genotyped through the MassARRAY platform. At the 2^nd stage, a larger sample including 400 patients and 392 controls was genotyped by the High Resolution Melt (HRM) method to confirm or rule out the associations with CHD. MLXIP expression level was quantified by the real time PCR in 65 peripheral blood samples. From the 21 studied single nucleotide polymorphisms (SNPs) of seven candidate genes: MLXIPL, MLXIP, MLX, ADIPOR1, VDR, SREBF1 and NR1H3, only one tag SNP rs4758685 (T→C) was found to be statistically associated with CHD (P-value = 0.02, Odds ratio (OR) of 0.83). After adjustment for the age, sex, lipid levels and diabetes, the association remained significant (P-value = 0.03). After adjustment for the hypertension, P-value became 0.20 although there was a significant difference in the allele distribution between the CHD patients with hypertension and the controls (P-value = 0.04, 406 vs 582). In conclusion, among the 21 tested SNPs, we identified a novel association between rs4758685 of MLXIP gene and CHD. The C allele of common variant rs4758685 interacted with hypertension, and was found to be protective against CHD in both allelic and genotypic models in Chinese Han population.     

The SNPs of Melanocortin 4 Receptor (MC4R) Associated with Body Weight in Beagle Dogs

Melanocortin 4 receptor (MC4R), which is associated with inherited human obesity, is involoved in food intake and body weight of mammals. To study the relationships between MC4R gene polymorphism and body weight in Beagle dogs, we detected and compared the nucleotide sequence of the whole coding region and 3′- and 5′- flanking regions of the dog MC4R gene (1214 bp). In 120 Beagle dogs, two SNPs (A420C, C895T) were identified and their relation with body weight was analyzed with RFLP-PCR method. The results showed that the SNP at A420C was significantly associated with canine body weight trait when it changed amino acid 101 of the MC4R protein from asparagine to threonine,while canine body weight variations were significant in female dogs when MC4R nonsense mutation at C895T. It suggested that the two SNPs might affect the MC4R gene’s function which was relative to body weight in Beagle dogs. Therefore, MC4R was a candidate gene for selecting different size dogs with the MC4R SNPs (A420C, C895T) being potentially valuable as a genetic marker.     

Linked genetic variants on chromosome 10 control ear morphology and body mass among dog breeds

Background

The domestic dog is a rich resource for mapping the genetic components of phenotypic variation due to its unique population history involving strong artificial selection. Genome-wide association studies have revealed a number of chromosomal regions where genetic variation associates with morphological characters that typify dog breeds. A region on chromosome 10 is among those with the highest levels of genetic differentiation between dog breeds and is associated with body mass and ear morphology, a common motif of animal domestication. We characterised variation in this region to uncover haplotype structure and identify candidate functional variants.

Results

We first identified SNPs that strongly associate with body mass and ear type by comparing sequence variation in a 3 Mb region between 19 breeds with a variety of phenotypes. We next genotyped a subset of 123 candidate SNPs in 288 samples from 46 breeds to identify the variants most highly associated with phenotype and infer haplotype structure. A cluster of SNPs that associate strongly with the drop ear phenotype is located within a narrow interval downstream of the gene MSRB3, which is involved in human hearing. These SNPs are in strong genetic linkage with another set of variants that correlate with body mass within the gene HMGA2, which affects human height. In addition we find evidence that this region has been under selection during dog domestication, and identify a cluster of SNPs within MSRB3 that are highly differentiated between dogs and wolves.

Conclusions

We characterise genetically linked variants that potentially influence ear type and body mass in dog breeds, both key traits that have been modified by selective breeding that may also be important for domestication. The finding that variants on long haplotypes have effects on more than one trait suggests that genetic linkage can be an important determinant of the phenotypic response to selection in domestic animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1702-2) contains supplementary material, which is available to authorized users.