Parallel Mapping and Simultaneous Sequencing Reveals Deletions in BCAN and FAM83H Associated with Discrete Inherited Disorders in a Domestic Dog Breed

The domestic dog (Canis familiaris) segregates more naturally-occurring diseases and phenotypic variation than any other species and has become established as an unparalled model with which to study the genetics of inherited traits. We used a genome-wide association study (GWAS) and targeted resequencing of DNA from just five dogs to simultaneously map and identify mutations for two distinct inherited disorders that both affect a single breed, the Cavalier King Charles Spaniel. We investigated episodic falling (EF), a paroxysmal exertion-induced dyskinesia, alongside the phenotypically distinct condition congenital keratoconjunctivitis sicca and ichthyosiform dermatosis (CKCSID), commonly known as dry eye curly coat syndrome. EF is characterised by episodes of exercise-induced muscular hypertonicity and abnormal posturing, usually occurring after exercise or periods of excitement. CKCSID is a congenital disorder that manifests as a rough coat present at birth, with keratoconjunctivitis sicca apparent on eyelid opening at 10–14 days, followed by hyperkeratinisation of footpads and distortion of nails that develops over the next few months. We undertook a GWAS with 31 EF cases, 23 CKCSID cases, and a common set of 38 controls and identified statistically associated signals for EF and CKCSID on chromosome 7 (Praw 1.9×10−14; Pgenome = 1.0×10−5) and chromosome 13 (Praw 1.2×10−17; Pgenome = 1.0×10−5), respectively. We resequenced both the EF and CKCSID disease-associated regions in just five dogs and identified a 15,724 bp deletion spanning three exons of BCAN associated with EF and a single base-pair exonic deletion in FAM83H associated with CKCSID. Neither BCAN or FAM83H have been associated with equivalent disease phenotypes in any other species, thus demonstrating the ability to use the domestic dog to study the genetic basis of more than one disease simultaneously in a single breed and to identify multiple novel candidate genes in parallel

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

Author summary Hereditary ataxias, which are a group of disorders characterised by incoordination of movement, are typically incurable and there are often no disease-modifying treatments available. Canine hereditary ataxias are a notable group of movement disorders in dogs, and represent well characterised naturally occurring disease models of ataxia that can help improve our understanding of the underlying biology of the disorder in both dogs and humans. We used the whole genome sequences of two affected siblings to investigate the genetic cause of a slowly progressive form of hereditary ataxia in the Norwegian Buhund dog breed, and identified a single base change within the KCNIP4 gene. We have characterised the expression of KCNIP4 in the dog, and investigated the effect of the identified mutation. This gene has not previously been implicated in inherited ataxia in any species, and our findings suggest that this and related genes represent potential candidates for ataxia in future studies in other species. Our findings will allow dog breeders to avoid producing affected dogs, reduce the disease allele frequency, and eventually eliminate the disease from the breed, through the use of a DNA test. 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.Peer reviewe

An intronic LINE-1 insertion in MERTK is strongly associated with retinopathy in Swedish Vallhund dogs

The domestic dog segregates a significant number of inherited progressive retinal diseases, several of which mirror human retinal diseases and which are collectively termed progressive retinal atrophy (PRA). In 2014, a novel form of PRA was reported in the Swedish Vallhund breed, and the disease was mapped to canine chromosome 17. The causal mutation was not identified, but expression analyses of the retinas of affected Vallhunds demonstrated a 6-fold increased expression of the MERTK gene compared to unaffected dogs. Using 24 retinopathy cases and 97 controls with no clinical signs of retinopathy, we replicated the chromosome 17 association in Swedish Vallhunds from the UK and aimed to elucidate the causal variant underlying this association using whole genome sequencing (WGS) of an affected dog. This revealed a 6-8 kb insertion in intron 1 of MERTK that was not present in WGS of 49 dogs of other breeds. Sequencing and BLASTN analysis of the inserted segment was consistent with the insertion comprising a full-length intact LINE-1 retroelement. Testing of the LINE-1 insertion for association with retinopathy in the UK set of 24 cases and 97 controls revealed a strong statistical association (P-value 6.0 x 10(-11)) that was subsequently replicated in the original Finnish study set (49 cases and 89 controls (P-value 4.3 x 10(-19)). In a pooled analysis of both studies (73 cases and 186 controls), the LINE-1 insertion was associated with a similar to 20-fold increased risk of retinopathy (odds ratio 23.41, 95% confidence intervals 10.99-49.86, P-value 1.3 x 10(-27)). Our study adds further support for regulatory disruption of MERTK in Swedish Vallhund retinopathy; however, further work is required to establish a functional overexpression model. Future work to characterise the mechanism by which this intronic mutation disrupts gene regulation will further improve the understanding of MERTK biology and its role in retinal function.Peer reviewe

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

Funder: Italian Spinone Club of Great BritainFunder: Kennel Club Charitable TrustGenotype 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

Both Paraoxonase-1 Genotype and Activity Do Not Predict the Risk of Future Coronary Artery Disease; the EPIC-Norfolk Prospective Population Study

Paraoxonase-1 (PON1) is an antioxidant enzyme, that resides on high-density lipoprotein (HDL). PON1-activity, is heavily influenced by the PON1-Q192R polymorphism. PON1 is considered to protect against atherosclerosis, but it is unclear whether this relation is independent of its carrier, HDL. In order to evaluate the atheroprotective potential of PON1, we assessed the relationships among PON1-genotype, PON1-activity and risk of future coronary artery disease (CAD), in a large prospective case-control study. Methodology/Principal Findings: Cases (n = 1138) were apparently healthy men and women aged 45-79 years who developed fatal or nonfatal CAD during a mean follow-up of 6 years. Controls (n = 2237) were matched by age, sex and enrollment time. PON1-activity was similar in cases and controls (60.7 +/- 645.3 versus 62.6 +/- 645.8 U/L, p = 0.3) and correlated with HDL-cholesterol levels (r = 0.16, p < 0.0001). The PON1-Q192R polymorphism had a profound impact on PON1-activity, but did not predict CAD risk (Odds Ratio [OR] per R allele 0.98[0.84-1.15], p = 0.8). Using conditional logistic regression, quartiles of PON1-activity showed a modest inverse relation with CAD risk (OR for the highest versus the lowest quartile 0.77[0.63-0.95], p = 0.01; p-trend = 0.06). PON1-activity adjusted for Q192R polymorphism correlated better with HDL-cholesterol (r = 0.26, p < 0.0001) and more linearly predicted CAD risk (0.79[0.64-0.98], p = 0.03; p-trend = 0.008). However, these relationships were abolished after adjustment for HDL (particles-cholesterol-size) and apolipoprotein A-l (0.94[0.74-1.18], p-trend = 0.3). Conclusions/Significance: This study, shows that PON1-activity inversely relates to CAD risk, but not independent of HDL, due to its close association with the HDL-particle. These data strongly suggest that a low PON1-activity is not a causal factor in atherogenesi

A synonymous germline variant in a gene encoding a cell adhesion molecule is associated with cutaneous mast cell tumour development in Labrador and Golden Retrievers

Mast cell tumours are the most common type of skin cancer in dogs, representing a significant concern in canine health. The molecular pathogenesis is largely unknown, but breed-predisposition for mast cell tumour development suggests the involvement of inherited genetic risk factors in some breeds. In this study, we aimed to identify germline risk factors associated with the development of mast cell tumours in Labrador Retrievers, a breed with an elevated risk of mast cell tumour development. Using a methodological approach that combined a genome-wide association study, targeted next generation sequencing, and TaqMan genotyping, we identified a synonymous variant in the DSCAM gene on canine chromosome 31 that is associated with mast cell tumours in Labrador Retrievers. DSCAM encodes a cell-adhesion molecule. We showed that the variant has no effect on the DSCAM mRNA level but is associated with a significant reduction in the level of the DSCAM protein, suggesting that the variant affects the dynamics of DSCAM mRNA translation. Furthermore, we showed that the variant is also associated with mast cell tumours in Golden Retrievers, a breed that is closely related to Labrador Retrievers and that also has a predilection for mast cell tumour development. The variant is common in both Labradors and Golden Retrievers and consequently is likely to be a significant genetic contributor to the increased susceptibility of both breeds to develop mast cell tumours. The results presented here not only represent an important contribution to the understanding of mast cell tumour development in dogs, as they highlight the role of cell adhesion in mast cell tumour tumourigenesis, but they also emphasise the potential importance of the effects of synonymous variants in complex diseases such as cancer. Author summary The combination of various genetic and environmental risk factors makes the understanding of the molecular circuitry behind complex diseases, like cancer, a major challenge. The homogeneous nature of pedigree dog breed genomes makes these dogs ideal for the identification of both simple disease-causing genetic variants and genetic risk factors for complex diseases. Mast cell tumours are the most common type of canine skin cancer, and one of the most common cancers affecting dogs of most breeds. Several breeds, including Labrador Retrievers (which represent one of the most popular dog breeds), have an elevated risk of mast cell tumour development. Here, by using a methodological approach that combined different techniques, we identified a common inherited synonymous variant, that predisposes Labrador Retrievers to mast cell tumour development. Interestingly, we showed that this variant, despite its synonymous nature, appears to have an effect on translation dynamics as it is associated with reduced levels of DSCAM, a cell adhesion molecule. The results presented here reveal dysregulation of cell adhesion to be an important factor in mast cell tumour pathogenesis, and also highlight the important role that synonymous variants can play in complex diseases

Common variants at 10 Genomic loci influence hemoglobin A1C levels via glycemic and nonglycemic pathways

OBJECTIVE: Glycated hemoglobin (HbA(1c)), used to monitor and diagnose diabetes, is influenced by average glycemia over a 2- to 3-month period. Genetic factors affecting expression, turnover, and abnormal glycation of hemoglobin could also be associated with increased levels of HbA(1c). We aimed to identify such genetic factors and investigate the extent to which they influence diabetes classification based on HbA(1c) levels. RESEARCH DESIGN AND METHODS: We studied associations with HbA(1c) in up to 46,368 nondiabetic adults of European descent from 23 genome-wide association studies (GWAS) and 8 cohorts with de novo genotyped single nucleotide polymorphisms (SNPs). We combined studies using inverse-variance meta-analysis and tested mediation by glycemia using conditional analyses. We estimated the global effect of HbA(1c) loci using a multilocus risk score, and used net reclassification to estimate genetic effects on diabetes screening. RESULTS: Ten loci reached genome-wide significant association with HbA(1c), including six new loci near FN3K (lead SNP/P value, rs1046896/P = 1.6 × 10(−26)), HFE (rs1800562/P = 2.6 × 10(−20)), TMPRSS6 (rs855791/P = 2.7 × 10(−14)), ANK1 (rs4737009/P = 6.1 × 10(−12)), SPTA1 (rs2779116/P = 2.8 × 10(−9)) and ATP11A/TUBGCP3 (rs7998202/P = 5.2 × 10(−9)), and four known HbA(1c) loci: HK1 (rs16926246/P = 3.1 × 10(−54)), MTNR1B (rs1387153/P = 4.0 × 10(−11)), GCK (rs1799884/P = 1.5 × 10(−20)) and G6PC2/ABCB11 (rs552976/P = 8.2 × 10(−18)). We show that associations with HbA(1c) are partly a function of hyperglycemia associated with 3 of the 10 loci (GCK, G6PC2 and MTNR1B). The seven nonglycemic loci accounted for a 0.19 (% HbA(1c)) difference between the extreme 10% tails of the risk score, and would reclassify ∼2% of a general white population screened for diabetes with HbA(1c). CONCLUSIONS: GWAS identified 10 genetic loci reproducibly associated with HbA(1c). Six are novel and seven map to loci where rarer variants cause hereditary anemias and iron storage disorders. Common variants at these loci likely influence HbA(1c) levels via erythrocyte biology, and confer a small but detectable reclassification of diabetes diagnosis by HbA(1c)

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

Funder: Lhasa Apso Breed CouncilFunder: Kennel Club Charitable Trust UKAbstract: 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

Genome-wide analyses implicate 33 loci in heritable dog osteosarcoma, including regulatory variants near CDKN2A/B

Background: Canine osteosarcoma is clinically nearly identical to the human disease, but is common and highly heritable, making genetic dissection feasible. Results: Through genome-wide association analyses in three breeds (greyhounds, Rottweilers, and Irish wolfhounds), we identify 33 inherited risk loci explaining 55% to 85% of phenotype variance in each breed. The greyhound locus exhibiting the strongest association, located 150 kilobases upstream of the genes CDKN2A/B, is also the most rearranged locus in canine osteosarcoma tumors. The top germline candidate variant is found at a >90% frequency in Rottweilers and Irish wolfhounds, and alters an evolutionarily constrained element that we show has strong enhancer activity in human osteosarcoma cells. In all three breeds, osteosarcoma-associated loci and regions of reduced heterozygosity are enriched for genes in pathways connected to bone differentiation and growth. Several pathways, including one of genes regulated by miR124, are also enriched for somatic copy-number changes in tumors. Conclusions: Mapping a complex cancer in multiple dog breeds reveals a polygenic spectrum of germline risk factors pointing to specific pathways as drivers of disease