Two FGF4 retrogenes on chromosomes 12 (12-FGF4RG) and 18 (18-FGF4RG) contribute to short-limbed phenotypes in dogs. 12-FGF4RG has also been associated with intervertebral disc disease (IVDD). Both of these retrogenes were found to be widespread among dog breeds with allele frequencies ranging from 0.02 to 1; however, their additive contribution to disease is unknown. Surgical cases of IVDD (n = 569) were evaluated for age of onset, disc calcification, and genotypes for the FGF4 retrogenes. Multivariable linear regression analysis identified the presence of one or two copies of 12-FGF4RG associated with significantly younger age at first surgery in a dominant manner. 18-FGF4RG had only a minor effect in dogs with one copy. Multivariable logistic regression showed that 12-FGF4RG had an additive effect on radiographic disc calcification, while 18-FGF4RG had no effect. Multivariable logistic regression using mixed breed cases and controls identified only 12-FGF4RG as highly associated with disc herniation in a dominant manner (Odds Ratio, OR, 18.42, 95% Confidence Interval (CI) 7.44 to 50.26; p < 0.001). The relative risk for disc surgery associated with 12-FGF4RG varied from 5.5 to 15.1 within segregating breeds and mixed breeds. The FGF4 retrogene on CFA12 acts in a dominant manner to decrease the age of onset and increase the overall risk of disc disease in dogs. Other modifiers of risk may be present within certain breeds, including the FGF4 retrogene on CFA18.

Primary ciliary dyskinesia (PCD) is a hereditary defect of motile cilia in humans and several domestic animal species. Typical clinical findings are chronic recurrent infections of the respiratory tract and fertility problems. We analyzed an Alaskan Malamute family, in which two out of six puppies were affected by PCD. The parents were unaffected suggesting autosomal recessive inheritance. Linkage and homozygosity mapping defined critical intervals comprising ~118 Mb. Whole genome sequencing of one case and comparison to 601 control genomes identified a disease associated frameshift variant, c.43delA, in the NME5 gene encoding a sparsely characterized protein associated with ciliary function. Nme5^-/- knockout mice exhibit doming of the skull, hydrocephalus and sperm flagellar defects. The genotypes at NME5:c.43delA showed the expected co-segregation with the phenotype in the Alaskan Malamute family. An additional unrelated Alaskan Malamute with PCD and hydrocephalus that became available later in the study was also homozygous mutant at the NME5:c.43delA variant. The mutant allele was not present in more than 1000 control dogs from different breeds. Immunohistochemistry demonstrated absence of the NME5 protein from nasal epithelia of an affected dog. We therefore propose NME5:c.43delA as the most likely candidate causative variant for PCD in Alaskan Malamutes. These findings enable genetic testing to avoid the unintentional breeding of affected dogs in the future. Furthermore, the results of this study identify NME5 as a novel candidate gene for unsolved human PCD and/or hydrocephalus cases. Motile cilia are required for clearing mucous, infectious agents and inhaled dust from the airways. Primary ciliary dyskinesia (PCD) is a hereditary defect of motile cilia. Clinical findings may include recurrent airway infections, fertility problems, and sometimes hydrocephalus. We analyzed an Alaskan Malamute family, in which two out of six puppies were affected by an autosomal recessive form of PCD. Whole genome sequencing of an affected dog identified a one base pair deletion in the NME5 gene, c.43delA, leading to an early frame-shift and premature stop codon. Later in the study, we became aware of a previously published Alaskan Malamute with PCD involving respiratory infections and hydrocephalus. We observed perfect concordance of the NME5 genotypes with the PCD phenotype in all three affected Alaskan Malamutes and more than 1000 controls. The fact that the third case, which had no documented close relationship to the initial two cases, was homozygous for the same rare mutant NME5 allele, strongly supports our hypothesis that NME5:c.43delA causes the PCD phenotype. We confirmed absence of NME5 protein expression in nasal epithelium of an affected dog. Our results enable genetic testing in dogs and identify NME5 as novel candidate gene for unsolved human PCD cases.

Two FGF4 retrogenes (FGF4L1 on chromosome 18 and FGF4L2 on chromosome 12) have been identified to cause dwarfism across many dog breeds. Some breeds are nearly homozygous for both retrogenes (e.g., Dachshunds) and others are homozygous for just one (e.g., Beagles and Scottish Terriers). Since most breeds do not segregate both of these retrogenes, it is challenging to evaluate their individual effects on long bone length and body size. We identified two dog breeds selected for hunting ability, the Alpine Dachsbracke and the Schweizer Niederlaufhund, that segregate both of these retrogenes. Using individual measurements of height at the shoulder, back length, head width, thorax depth and width, and thoracic limb measurements, we evaluated the combined effects of FGF4 retrogenes within these breeds. We applied multivariable linear regression analysis to determine the effects of retrogene copy numbers on the measurements. Copy numbers of both retrogenes had significant effects reducing height at the shoulders and antebrachial length, with FGF4L1 having a much greater effect than FGF4L2. FGF4L1 alone influenced the degree of carpal valgus and FGF4L2 alone increased head width. Neither retrogene had an effect on thorax width or depth. Selectively breeding dogs with FGF4L1 and without FGF4L2 would likely lead to a reduction in the FGF4L2-related risk of intervertebral disc herniation while maintaining the reduction in leg length resulting from FGF4L1.

Cutaneous lupus erythematosus (CLE) in humans encompasses multiple subtypes that exhibit a wide array of skin lesions and, in some cases, are associated with the development of systemic lupus erythematosus (SLE). We investigated dogs with exfoliative cutaneous lupus erythematosus (ECLE), a dog-specific form of chronic CLE that is inherited as a monogenic autosomal recessive trait. A genome-wide association study (GWAS) with 14 cases and 29 controls confirmed a previously published result that the causative variant maps to chromosome 18. Autozygosity mapping refined the ECLE locus to a 493 kb critical interval. Filtering of whole genome sequence data from two cases against 654 controls revealed a single private protein-changing variant in this critical interval, UNC93B1:c.1438C>A or p.Pro480Thr. The homozygous mutant genotype was exclusively observed in 23 ECLE affected German Shorthaired Pointers and an ECLE affected Vizsla, but absent from 845 controls. UNC93B1 is a transmembrane protein located in the endoplasmic reticulum and endolysosomes, which is required for correct trafficking of several Toll-like receptors (TLRs). The p.Pro480Thr variant is predicted to affect the C-terminal tail of the UNC93B1 that has recently been shown to restrict TLR7 mediated autoimmunity via an interaction with syndecan binding protein (SDCBP). The functional knowledge on UNC93B1 strongly suggests that p.Pro480Thr is causing ECLE in dogs. These dogs therefore represent an interesting spontaneous model for human lupus erythematosus. Our results warrant further investigations of whether genetic variants affecting the C-terminus of UNC93B1 might be involved in specific subsets of CLE or SLE cases in humans and other species.

Two transcribed retrocopies of the fibroblast growth factor 4 (FGF4) gene have previously been described in the domestic dog. An FGF4 retrocopy on chr18 is associated with disproportionate dwarfism, while an FGF4 retrocopy on chr12 is associated with both disproportionate dwarfism and intervertebral disc disease (IVDD). In this study, whole-genome sequencing data were queried to identify other FGF4 retrocopies that could be contributing to phenotypic diversity in canids. Additionally, dogs with surgically confirmed IVDD were assayed for novel FGF4 retrocopies. Five additional and distinct FGF4 retrocopies were identified in canids including a copy unique to red wolves (Canis rufus). The FGF4 retrocopies identified in domestic dogs were identical to domestic dog FGF4 haplotypes, which are distinct from modern wolf FGF4 haplotypes, indicating that these retrotransposition events likely occurred after domestication. The identification of multiple, full length FGF4 retrocopies with open reading frames in canids indicates that gene retrotransposition events occur much more frequently than previously thought and provide a mechanism for continued genetic and phenotypic diversity in canids.

Gene retrocopies arise from the reverse transcription and insertion into the genome of processed mRNA transcripts. Although many retrocopies have acquired mutations that render them functionally inactive, most mammals retain active LINE-1 sequences capable of producing new retrocopies. New retrocopies, referred to as retro copy number variants (retroCNVs), may not be identified by standard variant calling techniques in high-throughput sequencing data. Although multiple functional FGF4 retroCNVs have been associated with skeletal dysplasias in dogs, the full landscape of canid retroCNVs has not been characterized. Here, retroCNV discovery was performed on a whole-genome sequencing data set of 293 canids from 76 breeds. We identified retroCNV parent genes via the presence of mRNA-specific 30-mers, and then identified retroCNV insertion sites through discordant read analysis. In total, we resolved insertion sites for 1911 retroCNVs from 1179 parent genes, 1236 of which appeared identical to their parent genes. Dogs had on average 54.1 total retroCNVs and 1.4 private retroCNVs. We found evidence of expression in testes for 12% (14/113) of the retroCNVs identified in six Golden Retrievers, including four chimeric transcripts, and 97 retroCNVs also had significantly elevated F_ST across dog breeds, possibly indicating selection. We applied our approach to a subset of human genomes and detected an average of 4.2 retroCNVs per sample, highlighting a 13-fold relative increase of retroCNV frequency in dogs. Particularly in canids, retroCNVs are a largely unexplored source of genetic variation that can contribute to genome plasticity and that should be considered when investigating traits and diseases.