In dogs, symmetrical lupoid onychodystrophy (SLO) results in nail loss and an abnormal regrowth of the claws. In Bearded Collies, an autoimmune nature has been suggested because certain dog leukocyte antigen (DLA) class II haplotypes are associated with the condition. A genome-wide association study of the Bearded Collie revealed two regions of association that conferred risk for disease: one on canine chromosome (CFA) 12 that encompasses the DLA genes, and one on CFA17. Case-control association was employed on whole genome sequencing data to uncover putative causative variants in SLO within the CFA12 and CFA17 associated regions. Genotype imputation was then employed to refine variants of interest. Although no SLO-associated protein-coding variants were identified on CFA17, multiple variants, many with predicted damaging effects, were identified within potential candidate genes on CFA12. Furthermore, many potentially damaging alleles were fully correlated with the presence of DLA class II risk haplotypes for SLO, suggesting that the variants may reflect DLA class II haplotype association with disease or vice versa. Strong linkage disequilibrium in the region precluded the ability to isolate and assess the individual or combined effect of variants on disease development. Nonetheless, all were predictive of risk for SLO and, with judicious assessment, their application in selective breeding may prove useful to reduce the incidence of SLO in the breed.

Symmetrical lupoid onychodystrophy (SLO) is characterized by inflammation of the nail bed and nail sloughing that causes affected dogs considerable pain. Disease etiology remains unclear, although an autoimmune component is suspected. A genome-wide association study on Bearded Collies revealed regions of association on canine chromosomes (CFA) 12 and 17. The large region of association on CFA12 likely consists of two smaller linked regions, both of which are also linked to the dog leukocyte antigen (DLA) class II genes. Dogs homozygous for the alternate allele at the top CFA12 SNP also carried two DLA class II risk haplotypes for SLO, and this locus explained most of the increased risk for disease seen throughout the CFA12 region of association. A stronger peak was seen on CFA17 when analysis was done solely on dogs that carried DLA class II risk haplotypes for SLO. The majority of SLO dogs carried a homozygous alternate genotype on CFA12 and at least one CFA17 risk haplotype. Our findings offer progress toward uncovering the genetic basis of SLO. While the contribution of the CFA17 region remains unclear, both CFA12 and CFA17 regions are significantly associated with SLO disease expression in the Bearded Collie and contain potential candidate genes for this disease.

A frameshift deletion variant in the Wnt pathway gene dishevelled 2 (DVL2) is associated with a truncated, kinked tail ("screw tail") in English Bulldogs, French Bulldogs and Boston Terriers. These breeds are also characterized by distinctive morphological traits, including a wide head, flat face and short-limbed dwarfism, which are characteristic of Robinow syndrome in humans, caused by defects in genes such as DVL1 and DVL3. Based on these phenotypic and genetic similarities, it has previously been hypothesized that the canine DVL2 variant results in a syndromic phenotype called the Robinow-like syndrome. In our study, we investigated the distribution of the DVL2 variant in 1954 dogs from 15 breeds, identifying breeds with allele variation and enabling the dissection of the genotype-phenotype correlation for the first time. With CT examinations in American Staffordshire Terriers, we confirmed that the DVL2 allele is associated with caudal vertebral malformations and a brachycephalic phenotype. We also hypothesize that the variant may be linked to additional health conditions, including brachycephalic obstructive airway syndrome and congenital heart defects. Altogether, our study strengthens the role of DVL2 as one of the contributors to the "bulldog type" morphology and features on the spectrum of human Robinow syndrome.

Idiopathic epilepsy (IE) has been known to be inherited in the Belgian Tervuren for many decades. Risk genotypes for IE in this breed have recently been identified on Canis familiaris chromosomes (CFA) 14 and 37. In the current study, the allele frequencies of these loci were analyzed to determine whether dog breeders had employed a purposeful selection against IE, leading to a reduction in risk-associated allele frequency within the breed over time. The allele frequencies of two generational groupings of Belgian Tervuren with and without IE were compared. Allele frequencies for risk-associated alleles on CFA14 were unchanged between 1985 and 2015, whereas those on CFA37 increased during that time in the control population (p < 0.05). In contrast, dogs with IE showed a decrease (p < 0.05) in the IE risk-associated allele frequency at the CFA37 locus. Seizure prevalence in the Belgian Tervuren appears to be increasing. These results suggest that, despite awareness that IE is inherited, selection against IE has not been successful.

Background

Primary hypoadrenocorticism (or Addison's disease, AD) is an autoimmune disease that results in destruction of the adrenal cortex and consequent adrenal insufficiency. The disease has been described in purebred and mixed breed dogs, although some breeds, including the Bearded Collie, are at increased risk for AD. Candidate gene approaches have yielded few associations that appear to be breed-specific. A single other genome-wide association study reported no significant regions of association for AD in Standard Poodles. The present study aimed to identify genomic regions of association for canine AD in Bearded Collies.

Results

Our study consists of the first genome-wide association analysis to identify a genome-wide significant region of association with canine AD (CFA18). Peaks of suggestive association were also noted on chromosomes 11, 16 and 29. Logistic regression analysis supported an additive effect of risk genotypes at these smaller effect loci on the probability of disease associated with carrying a risk genotype on CFA18. Potential candidate genes involved in adrenal steroidogenesis, regulation of immune responses and/or inflammation were identified within the associated regions of chromosomes 11 and 16. The gene-poor regions of chromosomes 18 and 29 may, however, harbor regulatory sequences that can modulate gene expression and contribute to disease susceptibility.

Conclusion

Our findings support the polygenic and complex nature of canine AD and identified a strongly associated locus on CFA18 that, when combined with three other smaller effect loci, was predictive of disease. The results offer progress in the identification of susceptibility loci for canine AD in the Bearded Collie. Further studies are needed to confirm association with the suggested candidate genes and identify actual causative mutations involved with AD susceptibility in this breed.

Maternally skewed transmission of traits has been associated with genomic imprinting and oocyte-derived mRNA. We report canine congenital eye malformations, caused by an amino acid deletion (K12del) near the N terminus of retinol-binding protein (RBP4). The disease is only expressed when both dam and offspring are deletion homozygotes. RBP carries vitamin A (retinol) from hepatic stores to peripheral tissues, including the placenta and developing eye, where it is required to synthesize retinoic acid. Gestational vitamin A deficiency is a known risk factor for ocular birth defects. The K12del mutation disrupts RBP folding in vivo, decreasing its secretion from hepatocytes to serum. The maternal penetrance effect arises from an impairment in the sequential transfer of retinol across the placenta, via RBP encoded by maternal and fetal genomes. Our results demonstrate a mode of recessive maternal inheritance, with a physiological basis, and they extend previous observations on dominant-negative RBP4 alleles in humans.

An idiopathic epilepsy (IE) risk haplotype on canine chromosome (CFA) 14 has been reported to interact with the CFA37 common risk haplotype in the Belgian shepherd (BS). Additional IE cases and control dogs were genotyped for the risk haplotypes to validate these previous findings. In the new cohort, the interaction between the two regions significantly elevated IE risk. When the haplotypes were analyzed individually, particular haplotypes on both CFA14 (ACTG) and 37 (GG) were associated with elevated IE risk, though only the CFA37 AA was significantly associated (p < 0.003) with reduced risk in the new cohort. However, the CFA14 ACTG risk was statistically significant when the new and previous cohort data were combined. The frequency of the ACTG haplotype was four-fold higher in BS dogs than in other breeds. Whole genome sequence analysis revealed that a 3-base pair predicted disruptive insertion in the RAPGEF5 gene, which is adjacent to the CFA14 risk haplotype. RAPGEF5 is involved in the Wnt-β-catenin signaling pathway that is crucial for normal brain function. Although this risk variant does not fully predict the likelihood of a BS developing IE, the association with a variant in a candidate gene may provide insight into the genetic control of canine IE.

Background

Primary hypoadrenocorticism (Addison's disease, AD) and symmetrical lupoid onychodystrophy (SLO) are two clinical conditions with an autoimmune etiology that occur in multiple dog breeds. In man, autoimmunity is associated with polymorphisms in immune-related genes that result in a reduced threshold for, or defective regulation of, T cell activation. The major histocompatibility complex (MHC) class II genes encode molecules that participate in these functions, and polymorphisms within these genes have been associated with autoimmune conditions in dogs and humans. Bearded collies have a relatively high prevalence of autoimmune diseases, particularly AD and SLO. Our study assessed the relationship between particular MHC (dog leukocyte antigen, DLA) class II haplotypes and the two autoimmune diseases most common in this breed. Moreover, five unrelated breeds at increased risk for AD were studied for comparative purposes and analyzed in the context of extant literature.

Results

A single DLA class II three-locus haplotype, determined by sequence-based typing, was associated with increased risk for AD (DLA-DRB1*009:01/DQA1*001:01/DQB1*008:02) in bearded collies. Comparative analysis with the five additional breeds showed limited allele sharing, with DQA1*001:01 and DQB1*002:01 being the only alleles observed in all breeds. A distinct three-locus risk haplotype (DLA-DRB1*001:01/DQA1*001:01/DQB1*002:01) was associated with AD in the West Highland white terrier and Leonberger. Two different risk haplotypes were associated with increased risk for SLO in the bearded collie (DLA-DRB1*018:01/DQA1*001:01/DQB1*002:01 and DLA-DRB1*018:01/DQA1*001:01/ DQB1*008:02).

Conclusion

Two-locus DQ haplotypes composed of DLA-DQA1*001:01 in association with DLA-DQB1*002:01 or DLA-DQB1*008:02 make up the four risk haplotypes identified in the present study and are also found in other risk haplotypes previously associated with diabetes mellitus and hypothyroidism across different dog breeds. Our findings build upon previously published data to suggest that this two-locus (DQ) model serves as a good indicator for susceptibility to multiple organ-specific autoimmune diseases in the canine population. However, it is also clear that additional loci are necessary for actual disease expression. Investigation of affected and unaffected dogs carrying these predisposing DQ haplotype signatures may allow for the identification of those additional genetic components that determine autoimmune disease expression and organ specificity.

Domestic dogs can suffer from hearing losses that can have profound impacts on working ability and quality of life. We have identified a type of adult-onset hearing loss in Border Collies that appears to have a genetic cause, with an earlier age of onset (3-5 years) than typically expected for aging dogs (8-10 years). Studying this complex trait within pure breeds of dog may greatly increase our ability to identify genomic regions associated with risk of hearing impairment in dogs and in humans. We performed a genome-wide association study (GWAS) to detect loci underlying adult-onset deafness in a sample of 20 affected and 28 control Border Collies. We identified a region on canine chromosome 6 that demonstrates extended support for association surrounding SNP Chr6.25819273 (p-value = 1.09 × 10(-13)). To further localize disease-associated variants, targeted next-generation sequencing (NGS) of one affected and two unaffected dogs was performed. Through additional validation based on targeted genotyping of additional cases (n = 23 total) and controls (n = 101 total) and an independent replication cohort of 16 cases and 265 controls, we identified variants in USP31 that were strongly associated with adult-onset deafness in Border Collies, suggesting the involvement of the NF-κB pathway. We found additional support for involvement of RBBP6, which is critical for cochlear development. These findings highlight the utility of GWAS-guided fine-mapping of genetic loci using targeted NGS to study hereditary disorders of the domestic dog that may be analogous to human disorders.

Retinitis pigmentosa (RP) is the leading cause of blindness with nearly two million people affected worldwide. Many genes have been implicated in RP, yet in 30-80% of the RP patients the genetic cause remains unknown. A similar phenotype, progressive retinal atrophy (PRA), affects many dog breeds including the Miniature Schnauzer. We performed clinical, genetic and functional experiments to identify the genetic cause of PRA in the breed. The age of onset and pattern of disease progression suggested that at least two forms of PRA, types 1 and 2 respectively, affect the breed, which was confirmed by genome-wide association study that implicated two distinct genomic loci in chromosomes 15 and X, respectively. Whole-genome sequencing revealed a fully segregating recessive regulatory variant in type 1 PRA. The associated variant has a very recent origin based on haplotype analysis and lies within a regulatory site with the predicted binding site of HAND1::TCF3 transcription factor complex. Luciferase assays suggested that mutated regulatory sequence increases expression. Case-control retinal expression comparison of six best HAND1::TCF3 target genes were analyzed with quantitative reverse-transcriptase PCR assay and indicated overexpression of EDN2 and COL9A2 in the affected retina. Defects in both EDN2 and COL9A2 have been previously associated with retinal degeneration. In summary, our study describes two genetically different forms of PRA and identifies a fully penetrant variant in type 1 form with a possible regulatory effect. This would be among the first reports of a regulatory variant in retinal degeneration in any species, and establishes a new spontaneous dog model to improve our understanding of retinal biology and gene regulation while the affected breed will benefit from a reliable genetic testing.