A RAPGEF6 variant constitutes a major risk factor for laryngeal paralysis in dogs

Insights into the genetics of body size in the Bull Terrier

The Bull Terrier (Miniature) and Bull Terrier are two varieties of a dog breed historically divided by size. We identify variety-associated chromosomal regions identified using stratified genome-wide association analysis of 69 Bull Terriers (Miniature) and 33 Bull Terriers. Next, we assess the significance of possible functional variants for body size using height (N = 1458) and weight (N = 1282) of Dog10K individuals with breed-representative metrics available. Variants significant for size across breeds that are consistent with size alleles observed in four Bull Terriers and four Bull Terriers (Miniature) represented in Dog10K are highlighted. From five identified regions, two include genes already known to influence canine body size and a third contains a potential new height gene (ARFGEF3). Near LCORL, the most highly associated variant for height in Bull Terriers was chr3:91734656A>G (pAcross-breed height = 2.459 × 10-99) and for weight it was chr3:91706639G>A (pAcross-breed weight = 9.762 × 10-85). All Bull Terriers (including Miniature) were monomorphic for the derived allele at the known size variant in LCORL (chr3:91872822A>del). In the first exon of IGF2BP2, the derived allele at chr34:18694869-71ins>del significantly reduces both height and weight in Bull Terriers and across breeds (Dog10K breed representative height and weight) (pAcross-breed height = 1.65 × 10-9; pAcross-breed weight = 1.79 × 10-8). The derived allele of the missense variant in ARFGEF3 chr1:30793904G>A, XP_038382065.1 p.V243I significantly reduces breed representative height but not weight (pAcross-breed height = 0.01). The effects on the variants assessed are limited to small variants identified in the Dog10K resource using breed-representative sizes.

iDog: a multi-omics resource for canids study

iDog (https://ngdc.cncb.ac.cn/idog/) is a comprehensive public resource for domestic dogs (Canis lupus familiaris) and wild canids, designed to integrate multi-omics data and provide data services for the worldwide canine research community. Notably, iDog 2.0 features a 15-fold increase in genomic samples, including 29.55 million single nucleotide polymorphisms (SNPs) and 16.54 million insertions/deletions (InDels) from 1929 modern samples and 29.09 million SNPs from 111 ancient Canis samples. Additionally, 43487 breed-specific SNPs and 530 disease/trait-associated variants have been identified and integrated. The platform also includes data from 141 BioProjects involving gene expression analyses and a single-cell transcriptome module containing data from 105 057 Beagle hippocampus cells. iDog 2.0 also includes an epignome module that evaluates DNA methylation patterns across 547 samples and chromatin accessibility across 87 samples for the analysis of gene expression regulation. Additionally, it provies phenotypic data for 897 dog diseases, 3207 genotype-to-phenotype (G2P) pairs, and 349 dog disease-associated genes, along with two newly constructed ontologies for breed and disease standardization. Finally, 13 new analytical tools have been added. Given these enhancements, the updated iDog 2.0 is an invaluable resource for the global cannie research community.

Genomic analyses identify 15 susceptibility loci and reveal HDAC2, SOX2-OT, and IGF2BP2 in a naturally-occurring canine model of gastric cancer

Gastric cancer (GC) is the fifth most common human cancer worldwide, but the genetic etiology is largely unknown. We performed a Bayesian genome-wide association study and selection analyses in a naturally-occurring canine model of GC, the Belgian Tervuren and Sheepdog breeds, to elucidate underlying genetic risk factors. We identified 15 loci with over 90% predictive accuracy for the GC phenotype. Variant filtering revealed germline putative regulatory variants for the EPAS1 (HIF2A) and PTEN genes and a coding variant in CD101. Although closely related to Tervuren and Sheepdogs, Belgian Malinois rarely develop GC. Across-breed analyses uncovered protective haplotypes under selection in Malinois at SOX2-OT and IGF2BP2. Among Tervuren and Sheepdogs, HDAC2 putative regulatory variants were present at comparatively high frequency and were associated with GC. Here, we describe a complex genetic architecture governing GC in a dog model, including genes such as PDZRN3, that have not been associated with human GC.

Genome Wide Association Study with Imputed Whole Genome Sequence Data Identifies a 431 kb Risk Haplotype on CFA18 for Congenital Laryngeal Paralysis in Alaskan Sled Dogs

Congenital laryngeal paralysis (CLP) is an inherited disorder that affects the ability of the dog to exercise and precludes it from functioning as a working sled dog. Though CLP is known to occur in Alaskan sled dogs (ASDs) since 1986, the genetic mutation underlying the disease has not been reported. Using a genome-wide association study (GWAS), we identified a 708 kb region on CFA 18 harboring 226 SNPs to be significantly associated with CLP. The significant SNPs explained 47.06% of the heritability of CLP. We narrowed the region to 431 kb through autozygosity mapping and found 18 of the 20 cases to be homozygous for the risk haplotype. Whole genome sequencing of two cases and a control ASD, and comparison with the genome of 657 dogs from various breeds, confirmed the homozygous status of the risk haplotype to be unique to the CLP cases. Most of the dogs that were homozygous for the risk allele had blue eyes. Gene annotation and a gene-based association study showed that the risk haplotype encompasses genes implicated in developmental and neurodegenerative disorders. Pathway analysis showed enrichment of glycoproteins and glycosaminoglycans biosynthesis, which play a key role in repairing damaged nerves. In conclusion, our results suggest an important role for the identified candidate region in CLP.

A CNTNAP1 Missense Variant Is Associated with Canine Laryngeal Paralysis and Polyneuropathy

Laryngeal paralysis associated with a generalized polyneuropathy (LPPN) most commonly exists in geriatric dogs from a variety of large and giant breeds. The purpose of this study was to discover the underlying genetic and molecular mechanisms in a younger-onset form of this neurodegenerative disease seen in two closely related giant dog breeds, the Leonberger and Saint Bernard. Neuropathology of an affected dog from each breed showed variable nerve fiber loss and scattered inappropriately thin myelinated fibers. Using across-breed genome-wide association, haplotype analysis, and whole-genome sequencing, we identified a missense variant in the CNTNAP1 gene (c.2810G>A; p.Gly937Glu) in which homozygotes in both studied breeds are affected. CNTNAP1 encodes a contactin-associated protein important for organization of myelinated axons. The herein described likely pathogenic CNTNAP1 variant occurs in unrelated breeds at variable frequencies. Individual homozygous mutant LPPN-affected Labrador retrievers that were on average four years younger than dogs affected by geriatric onset laryngeal paralysis polyneuropathy could be explained by this variant. Pathologic changes in a Labrador retriever nerve biopsy from a homozygous mutant dog were similar to those of the Leonberger and Saint Bernard. The impact of this variant on health in English bulldogs and Irish terriers, two breeds with higher CNTNAP1 variant allele frequencies, remains unclear. Pathogenic variants in CNTNAP1 have previously been reported in human patients with lethal congenital contracture syndrome and hypomyelinating neuropathy, including vocal cord palsy and severe respiratory distress. This is the first report of contactin-associated LPPN in dogs characterized by a deleterious variant that most likely predates modern breed establishment.

YARS2 Missense Variant in Belgian Shepherd Dogs with Cardiomyopathy and Juvenile Mortality

Dog puppy loss by the age of six to eight weeks after normal development is relatively uncommon. Necropsy findings in two spontaneously deceased Belgian Shepherd puppies indicated an abnormal accumulation of material in several organs. A third deceased puppy exhibited mild signs of an inflammation in the central nervous system and an enteritis. The puppies were closely related, raising the suspicion of a genetic cause. Pedigree analysis suggested a monogenic autosomal recessive inheritance. Combined linkage and homozygosity mapping assigned the most likely position of a potential genetic defect to 13 genome segments totaling 82 Mb. The genome of an affected puppy was sequenced and compared to 645 control genomes. Three private protein changing variants were found in the linked and homozygous regions. Targeted genotyping in 96 Belgian Shepherd dogs excluded two of these variants. The remaining variant, YARS2:1054G>A or p.Glu352Lys, was perfectly associated with the phenotype in a cohort of 474 Belgian Shepherd dogs. YARS2 encodes the mitochondrial tyrosyl-tRNA synthetase 2 and the predicted amino acid change replaces a negatively charged and evolutionary conserved glutamate at the surface of the tRNA binding domain of YARS2 with a positively charged lysine. Human patients with loss-of-function variants in YARS2 suffer from myopathy, lactic acidosis, and sideroblastic anemia 2, a disease with clinical similarities to the phenotype of the studied dogs. The carrier frequency was 27.2% in the tested Belgian Shepherd dogs. Our data suggest YARS2:1054G>A as the candidate causative variant for the observed juvenile mortality.