Characterisation of canine KCNIP4: A novel gene for cerebellar ataxia identified by whole-genome sequencing two affected Norwegian Buhund dogs.

A form of hereditary cerebellar ataxia has recently been described in the Norwegian Buhund dog breed. This study aimed to identify the genetic cause of the disease. Whole-genome sequencing of two Norwegian Buhund siblings diagnosed with progressive cerebellar ataxia was carried out, and sequences compared with 405 whole genome sequences of dogs of other breeds to filter benign common variants. Nine variants predicted to be deleterious segregated among the genomes in concordance with an autosomal recessive mode of inheritance, only one of which segregated within the breed when genotyped in additional Norwegian Buhunds. In total this variant was assessed in 802 whole genome sequences, and genotyped in an additional 505 unaffected dogs (including 146 Buhunds), and only four affected Norwegian Buhunds were homozygous for the variant. The variant identified, a T to C single nucleotide polymorphism (SNP) (NC_006585.3:g.88890674T>C), is predicted to cause a tryptophan to arginine substitution in a highly conserved region of the potassium voltage-gated channel interacting protein KCNIP4. This gene has not been implicated previously in hereditary ataxia in any species. Evaluation of KCNIP4 protein expression through western blot and immunohistochemical analysis using cerebellum tissue of affected and control dogs demonstrated that the mutation causes a dramatic reduction of KCNIP4 protein expression. The expression of alternative KCNIP4 transcripts within the canine cerebellum, and regional differences in KCNIP4 protein expression, were characterised through RT-PCR and immunohistochemistry respectively. The voltage-gated potassium channel protein KCND3 has previously been implicated in spinocerebellar ataxia, and our findings suggest that the Kv4 channel complex KCNIP accessory subunits also have an essential role in voltage-gated potassium channel function in the cerebellum and should be investigated as potential candidate genes for cerebellar ataxia in future studies in other species.

Improving the resolution of canine genome-wide association studies using genotype imputation: A study of two breeds.

Genotype imputation using a reference panel that combines high-density array data and publicly available whole genome sequence consortium variant data is potentially a cost-effective method to increase the density of extant lower-density array datasets. In this study, three datasets (two Border Collie; one Italian Spinone) generated using a legacy array (Illumina CanineHD, 173 662 SNPs) were utilised to assess the feasibility and accuracy of this approach and to gather additional evidence for the efficacy of canine genotype imputation. The cosmopolitan reference panels used to impute genotypes comprised dogs of 158 breeds, mixed breed dogs, wolves and Chinese indigenous dogs, as well as breed-specific individuals genotyped using the Axiom Canine HD array. The two Border Collie reference panels comprised 808 individuals including 79 Border Collies and 426 326 or 426 332 SNPs; and the Italian Spinone reference panel comprised 807 individuals including 38 Italian Spinoni and 476 313 SNPs. A high accuracy for imputation was observed, with the lowest accuracy observed for one of the Border Collie datasets (mean R2  = 0.94) and the highest for the Italian Spinone dataset (mean R2  = 0.97). This study's findings demonstrate that imputation of a legacy array study set using a reference panel comprising both breed-specific array data and multi-breed variant data derived from whole genomes is effective and accurate. The process of canine genotype imputation, using the valuable growing resource of publicly available canine genome variant datasets alongside breed-specific data, is described in detail to facilitate and encourage use of this technique in canine genetics.

A LINE-1 insertion situated in the promoter of IMPG2 is associated with autosomal recessive progressive retinal atrophy in Lhasa Apso dogs.

BACKGROUND: Canine progressive retinal atrophies are a group of hereditary retinal degenerations in dogs characterised by depletion of photoreceptor cells in the retina, which ultimately leads to blindness. PRA in the Lhasa Apso (LA) dog has not previously been clinically characterised or described in the literature, but owners in the UK are advised to have their dog examined through the British Veterinary Association/ Kennel Club/ International Sheep Dog Society (BVA/KC/ISDS) eye scheme annually, and similar schemes that are in operation in other countries. After the exclusion of 25 previously reported canine retinal mutations in LA PRA-affected dogs, we sought to identify the genetic cause of PRA in this breed. RESULTS: Analysis of whole-exome sequencing data of three PRA-affected LA and three LA without signs of PRA did not identify any exonic or splice site variants, suggesting the causal variant was non-exonic. We subsequently undertook a genome-wide association study (GWAS), which identified a 1.3 Mb disease-associated region on canine chromosome 33, followed by whole-genome sequencing analysis that revealed a long interspersed element-1 (LINE-1) insertion upstream of the IMPG2 gene. IMPG2 has previously been implicated in human retinal disease; however, until now no canine PRAs have been associated with this gene. The identification of this PRA-associated variant has enabled the development of a DNA test for this form of PRA in the breed, here termed PRA4 to distinguish it from other forms of PRA described in other breeds. This test has been used to determine the genotypes of over 900 LA dogs. A large cohort of genotyped dogs was used to estimate the allele frequency as between 0.07-0.1 in the UK LA population. CONCLUSIONS: Through the use of GWAS and subsequent sequencing of a PRA case, we have identified a LINE-1 insertion in the retinal candidate gene IMPG2 that is associated with a form of PRA in the LA dog. Validation of this variant in 447 dogs of 123 breeds determined it was private to LA dogs. We envisage that, over time, the developed DNA test will offer breeders the opportunity to avoid producing dogs affected with this form of PRA.

Integration of disease association and eQTL data using a Bayesian colocalisation approach highlights six candidate causal genes in immune-mediated diseases.

The genes and cells that mediate genetic associations identified through genome-wide association studies (GWAS) are only partially understood. Several studies that have investigated the genetic regulation of gene expression have shown that disease-associated variants are over-represented amongst expression quantitative trait loci (eQTL) variants. Evidence for colocalisation of eQTL and disease causal variants can suggest causal genes and cells for these genetic associations. Here, we used colocalisation analysis to investigate whether 595 genetic associations to ten immune-mediated diseases are consistent with a causal variant that regulates, in cis, gene expression in resting B cells, and in resting and stimulated monocytes. Previously published candidate causal genes were over-represented amongst genes exhibiting colocalisation (odds ratio > 1.5), and we identified evidence for colocalisation (posterior odds > 5) between cis eQTLs in at least one cell type and at least one disease for six genes: ADAM15, RGS1, CARD9, LTBR, CTSH and SYNGR1. We identified cell-specific effects, such as for CTSH, the expression of which in monocytes, but not in B cells, may mediate type 1 diabetes and narcolepsy associations in the chromosome 15q25.1 region. Our results demonstrate the utility of integrating genetic studies of disease and gene expression for highlighting causal genes and cell types.

Copper toxicosis in Bedlington terriers is associated with multiple independent genetic variants.

BACKGROUND: Bedlington terrier copper toxicosis (CT) is due to a homozygous exon deletion in COMMD1. CT also occurs in Bedlingtons lacking this deletion. An association with two ABCA12 single nuceotide polymorphism (SNP) splice variants was reported. Labrador retriever CT is associated with a missense mutation in ATP7B, and with a protective mutation in ATP7A. METHODS: Liver and DNA samples from 24 affected and 10 unaffected Bedlingtons were assessed for copper and genetic variants. Allelic frequencies were compared. The ATP7B mutation frequency was investigated in 144 dogs of other breeds. RESULTS: The ABCA12 SNPs showed no differences between groups. The COMMD1 deletion was less frequent in unaffected than in affected dogs and in affected dogs post-2001 than pre-2001. The ATP7B mutation was more frequent in affected than unaffected Bedlingtons. Thirty-five of 144 dogs of other breeds were homo- or heterozygous for the ATP7B mutation. The ATP7A mutation was absent from Bedlingtons. LIMITATIONS: Clinical information and qualitative copper measurements were unavailable for most dogs. CONCLUSION: The COMMD1 deletion remains present in Bedlington terriers but is no longer the primary cause of CT. ABCA12 SNPs were not associated with CT. The ATP7B:c.4358G>A mutation was significantly associated with Bedlington CT and was more common in dogs of this breed than in the 144 dogs of other breeds.

Multiocular defect in the Old English Sheepdog: A canine form of Stickler syndrome type II associated with a missense variant in the collagen-type gene COL11A1.

Multiocular defect has been described in different canine breeds, including the Old English Sheepdog. Affected dogs typically present with multiple and various ocular abnormalities. We carried out whole genome sequencing on an Old English Sheepdog that had been diagnosed with hereditary cataracts at the age of five and then referred to a board-certified veterinary ophthalmologist due to owner-reported visual deterioration. An ophthalmic assessment revealed that there was bilateral vitreal degeneration, macrophthalmos, and spherophakia in addition to cataracts. Follow-up consultations revealed cataract progression, retinal detachment, uveitis and secondary glaucoma. Whole genome sequence filtered variants private to the case, shared with another Old English Sheepdog genome and predicted to be deleterious were genotyped in an initial cohort of six Old English Sheepdogs (three affected by multiocular defect and three control dogs without evidence of inherited eye disease). Only one of the twenty-two variants segregated correctly with multiocular defect. The variant is a single nucleotide substitution, located in the collagen-type gene COL11A1, c.1775T>C, that causes an amino acid change, p.Phe1592Ser. Genotyping of an additional 14 Old English Sheepdogs affected by multiocular defect revealed a dominant mode of inheritance with four cases heterozygous for the variant. Further genotyping of hereditary cataract-affected Old English Sheepdogs revealed segregation of the variant in eight out of nine dogs. In humans, variants in the COL11A1 gene are associated with Stickler syndrome type II, also dominantly inherited.

Globetrotting strangles: the unbridled national and international transmission of Streptococcus equi between horses.

The equine disease strangles, which is characterized by the formation of abscesses in the lymph nodes of the head and neck, is one of the most frequently diagnosed infectious diseases of horses around the world. The causal agent, Streptococcus equi subspecies equi, establishes a persistent infection in approximately 10 % of animals that recover from the acute disease. Such 'carrier' animals appear healthy and are rarely identified during routine veterinary examinations pre-purchase or transit, but can transmit S. equi to naïve animals initiating new episodes of disease. Here, we report the analysis and visualization of phylogenomic and epidemiological data for 670 isolates of S. equi recovered from 19 different countries using a new core-genome multilocus sequence typing (cgMLST) web bioresource. Genetic relationships among all 670 S. equi isolates were determined at high resolution, revealing national and international transmission events that drive this endemic disease in horse populations throughout the world. Our data argue for the recognition of the international importance of strangles by the Office International des Épizooties to highlight the health, welfare and economic cost of this disease. The Pathogenwatch cgMLST web bioresource described herein is available for tailored genomic analysis of populations of S. equi and its close relative S. equi subspecies zooepidemicus that are recovered from horses and other animals, including humans, throughout the world. This article contains data hosted by Microreact.

Clinical, histopathological and genetic characterisation of oculoskeletal dysplasia in the Northern Inuit Dog.

Seven Northern Inuit Dogs (NID) were diagnosed by pedigree analysis with an autosomal recessive inherited oculoskeletal dysplasia (OSD). Short-limbed dwarfism, angular limb deformities and a variable combination of macroglobus, cataracts, lens coloboma, microphakia and vitreopathy were present in all seven dogs, while retinal detachment was diagnosed in five dogs. Autosomal recessive OSD caused by COL9A3 and COL9A2 mutations have previously been identified in the Labrador Retriever (dwarfism with retinal dysplasia 1-drd1) and Samoyed dog (dwarfism with retinal dysplasia 2-drd2) respectively; both of those mutations were excluded in all affected NID. Nine candidate genes were screened in whole genome sequence data; only one variant was identified that was homozygous in two affected NID but absent in controls. This variant was a nonsense single nucleotide polymorphism in COL9A3 predicted to result in a premature termination codon and a truncated protein product. This variant was genotyped in a total of 1,232 dogs. All seven affected NID were homozygous for the variant allele (T/T), while 31/116 OSD-unaffected NID were heterozygous for the variant (C/T) and 85/116 were homozygous for the wildtype allele (C/C); indicating a significant association with OSD (p = 1.41x10-11). A subset of 56 NID unrelated at the parent level were analysed to determine an allele frequency of 0.08, estimating carrier and affected rates to be 15% and 0.6% respectively in NID. All 1,109 non-NID were C/C, suggesting the variant is rare or absent in other breeds. Expression of retinal mRNA was similar between an OSD-affected NID and OSD-unaffected non-NID. In conclusion, a nonsense variant in COL9A3 is strongly associated with OSD in NID, and appears to be widespread in this breed.

Whole Genome Sequencing of Giant Schnauzer Dogs with Progressive Retinal Atrophy Establishes NECAP1 as a Novel Candidate Gene for Retinal Degeneration.

Canine progressive retinal atrophies (PRA) are genetically heterogeneous diseases characterized by retinal degeneration and subsequent blindness. PRAs are untreatable and affect multiple dog breeds, significantly impacting welfare. Three out of seven Giant Schnauzer (GS) littermates presented with PRA around four years of age. We sought to identify the causal variant to improve our understanding of the aetiology of this form of PRA and to enable development of a DNA test. Whole genome sequencing of two PRA-affected full-siblings and both unaffected parents was performed. Variants were filtered based on those segregating appropriately for an autosomal recessive disorder and predicted to be deleterious. Successive filtering against 568 canine genomes identified a single nucleotide variant in the gene encoding NECAP endocytosis associated 1 (NECAP1): c.544G>A (p.Gly182Arg). Five thousand one hundred and thirty canids of 175 breeds, 10 cross-breeds and 3 wolves were genotyped for c.544G>A. Only the three PRA-affected GS were homozygous (allele frequency in GS, excluding proband family = 0.015). In addition, we identified heterozygotes belonging to Spitz and Dachshund varieties, demonstrating c.544G>A segregates in other breeds of German origin. This study, in parallel with the known retinal expression and role of NECAP1 in clathrin mediated endocytosis (CME) in synapses, presents NECAP1 as a novel candidate gene for retinal degeneration in dogs and other species.

Chromosome contacts in activated T cells identify autoimmune disease candidate genes.

BACKGROUND: Autoimmune disease-associated variants are preferentially found in regulatory regions in immune cells, particularly CD4+ T cells. Linking such regulatory regions to gene promoters in disease-relevant cell contexts facilitates identification of candidate disease genes. RESULTS: Within 4 h, activation of CD4+ T cells invokes changes in histone modifications and enhancer RNA transcription that correspond to altered expression of the interacting genes identified by promoter capture Hi-C. By integrating promoter capture Hi-C data with genetic associations for five autoimmune diseases, we prioritised 245 candidate genes with a median distance from peak signal to prioritised gene of 153 kb. Just under half (108/245) prioritised genes related to activation-sensitive interactions. This included IL2RA, where allele-specific expression analyses were consistent with its interaction-mediated regulation, illustrating the utility of the approach. CONCLUSIONS: Our systematic experimental framework offers an alternative approach to candidate causal gene identification for variants with cell state-specific functional effects, with achievable sample sizes.

Statistical colocalization of genetic risk variants for related autoimmune diseases in the context of common controls.

Determining whether potential causal variants for related diseases are shared can identify overlapping etiologies of multifactorial disorders. Colocalization methods disentangle shared and distinct causal variants. However, existing approaches require independent data sets. Here we extend two colocalization methods to allow for the shared-control design commonly used in comparison of genome-wide association study results across diseases. Our analysis of four autoimmune diseases--type 1 diabetes (T1D), rheumatoid arthritis, celiac disease and multiple sclerosis--identified 90 regions that were associated with at least one disease, 33 (37%) of which were associated with 2 or more disorders. Nevertheless, for 14 of these 33 shared regions, there was evidence that the causal variants differed. We identified new disease associations in 11 regions previously associated with one or more of the other 3 disorders. Four of eight T1D-specific regions contained known type 2 diabetes (T2D) candidate genes (COBL, GLIS3, RNLS and BCAR1), suggesting a shared cellular etiology.

Fine mapping of type 1 diabetes susceptibility loci and evidence for colocalization of causal variants with lymphoid gene enhancers.

Genetic studies of type 1 diabetes (T1D) have identified 50 susceptibility regions, finding major pathways contributing to risk, with some loci shared across immune disorders. To make genetic comparisons across autoimmune disorders as informative as possible, a dense genotyping array, the Immunochip, was developed, from which we identified four new T1D-associated regions (P < 5 × 10(-8)). A comparative analysis with 15 immune diseases showed that T1D is more similar genetically to other autoantibody-positive diseases, significantly most similar to juvenile idiopathic arthritis and significantly least similar to ulcerative colitis, and provided support for three additional new T1D risk loci. Using a Bayesian approach, we defined credible sets for the T1D-associated SNPs. The associated SNPs localized to enhancer sequences active in thymus, T and B cells, and CD34(+) stem cells. Enhancer-promoter interactions can now be analyzed in these cell types to identify which particular genes and regulatory sequences are causal.