A Multi-Breed Genome-Wide Association Analysis for Canine Hypothyroidism Identifies a Shared Major Risk Locus on CFA12

Hypothyroidism is a complex clinical condition found in both humans and dogs, thought to be caused by a combination of genetic and environmental factors. In this study we present a multi-breed analysis of predisposing genetic risk factors for hypothyroidism in dogs using three high-risk breeds—the Gordon Setter, Hovawart and the Rhodesian Ridgeback. Using a genome-wide association approach and meta-analysis, we identified a major hypothyroidism risk locus shared by these breeds on chromosome 12 (p = 2.1x10-11). Further characterisation of the candidate region revealed a shared ~167 kb risk haplotype (4,915,018–5,081,823 bp), tagged by two SNPs in almost complete linkage disequilibrium. This breed-shared risk haplotype includes three genes (LHFPL5, SRPK1 and SLC26A8) and does not extend to the dog leukocyte antigen (DLA) class II gene cluster located in the vicinity. These three genes have not been identified as candidate genes for hypothyroid disease previously, but have functions that could potentially contribute to the development of the disease. Our results implicate the potential involvement of novel genes and pathways for the development of canine hypothyroidism, raising new possibilities for screening, breeding programmes and treatments in dogs. This study may also contribute to our understanding of the genetic etiology of human hypothyroid disease, which is one of the most common endocrine disorders in humans.

[1]  A. Toker,et al.  Akt-ing up on SRPK1: oncogene or tumor suppressor? , 2014, Molecular cell.

[2]  Leif Andersson,et al.  Domestic-animal genomics: deciphering the genetics of complex traits , 2004, Nature Reviews Genetics.

[3]  G. Happ.  Thyroiditis--a model canine autoimmune disease. , 1995, Advances in veterinary science and comparative medicine.

[4]  G. Szinnai.  Genetics of normal and abnormal thyroid development in humans. , 2014, Best practice & research. Clinical endocrinology & metabolism.

[5]  James A. Cuff,et al.  Genome sequence, comparative analysis and haplotype structure of the domestic dog , 2005, Nature.

[6]  R. Fisher.  On the Interpretation of χ2 from Contingency Tables, and the Calculation of P , 2018, Journal of the Royal Statistical Society Series A (Statistics in Society).

[7]  E. Ostrander,et al.  Dog star rising: the canine genetic system , 2004, Nature Reviews Genetics.

[8]  S. Riazuddin,et al.  Mutations of human TMHS cause recessively inherited non-syndromic hearing loss , 2006, Journal of Medical Genetics.

[9]  F. Trimarchi,et al.  Thyroid function in children with cystic fibrosis , 1982, European Journal of Pediatrics.

[10]  V. Lucke,et al.  Thyroid pathology in canine hypothyroidism. , 1983, Journal of comparative pathology.

[11]  R. Nachreiner,et al.  Prevalence of serum thyroid hormone autoantibodies in dogs with clinical signs of hypothyroidism. , 2002, Journal of the American Veterinary Medical Association.

[12]  Martyn Plummer,et al.  The R Journal , 2012 .

[13]  C. Mallinckrodt,et al.  Associations between Lymphocytic Thyroiditis, Hypothyroidism, and Thyroid Neoplasia in Beagles , 1996, Veterinary pathology.

[14]  D. Ferguson.  Testing for hypothyroidism in dogs. , 2007, The Veterinary clinics of North America. Small animal practice.

[15]  K. Lindblad-Toh,et al.  Genome-wide analyses implicate 33 loci in heritable dog osteosarcoma, including regulatory variants near CDKN2A/B , 2013, Genome Biology.

[16]  CT Mooney,et al.  Canine hypothyroidism: A review of aetiology and diagnosis , 2011, New Zealand veterinary journal.

[17]  Å. Hedhammar,et al.  Increased genetic risk or protection for canine autoimmune lymphocytic thyroiditis in Giant Schnauzers depends on DLA class II genotype. , 2010, Tissue antigens.

[18]  H. Shwachman,et al.  Abnormal thyroid function and response to iodides in patients with cystic fibrosis. , 1974, Transactions of the Association of American Physicians.

[19]  J. Massey.  Mapping the genes for complex canine autoimmune diseases , 2012 .

[20]  K. Johnson,et al.  A missense mutation in the previously undescribed gene Tmhs underlies deafness in hurry-scurry (hscy) mice , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[21]  J. Khan,et al.  A Compendium of Canine Normal Tissue Gene Expression , 2011, PloS one.

[22]  W. Beierwaltes,et al.  Dog thyroiditis: occurrence and similarity to Hashimoto's struma. , 1968, Endocrinology.

[23]  K. Lindblad-Toh,et al.  Genome-wide association mapping identifies multiple loci for a canine SLE-related disease complex , 2010, Nature Genetics.

[24]  Jon W. Huss,et al.  BioGPS: an extensible and customizable portal for querying and organizing gene annotation resources , 2009, Genome Biology.

[25]  H. Huson,et al.  Association of hypothyroid disease in Doberman Pinscher dogs with a rare major histocompatibility complex DLA class II haplotype. , 2006, Tissue antigens.

[26]  DY Gaitonde,et al.  Hypothyroidism: an update , 2012, American family physician.

[27]  Karl Pearson F.R.S..  X. On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling , 2009 .

[28]  S. Reid,et al.  Epidemiological, clinical, haematological and biochemical characteristics of canine hypothyroidism , 1999, Veterinary Record.

[29]  R. Nachreiner,et al.  Etiopathologic findings of canine hypothyroidism. , 2007, The Veterinary clinics of North America. Small animal practice.

[30]  Å. Hedhammar,et al.  Gender, age, breed and distribution of morbidity and mortality in insured dogs in Sweden during 1995 and 1996 , 2000, Veterinary Record.

[31]  J. Fyfe,et al.  Congenital hypothyroidism with goiter in Tenterfield terriers. , 2012, Journal of veterinary internal medicine.

[32]  Jussi Taipale,et al.  Human kinome analysis reveals novel kinases contributing to virus infection and retinoic-acid inducible gene I-induced type I and type III IFN gene expression , 2013, Innate immunity.

[33]  R. Fisher.  On the Interpretation of χ2 from Contingency Tables, and the Calculation of P , 2010 .

[34]  D. Bannasch,et al.  Expanded dog leukocyte antigen (DLA) single nucleotide polymorphism (SNP) genotyping reveals spurious class II associations. , 2011, Veterinary journal.

[35]  Y. Tomer,et al.  The influence of human leucocyte antigen (HLA) genes on autoimmune thyroid disease (AITD): results of studies in HLA‐DR3 positive AITD families , 2002, Clinical endocrinology.

[36]  Moudud Alam,et al.  hglm: A Package for Fitting Hierarchical Generalized Linear Models , 2010, R J..

[37]  John P A Ioannidis,et al.  Meta-analysis in genome-wide association studies. , 2009, Pharmacogenomics.

[38]  Yurii S. Aulchenko,et al.  Region-Based Association Analysis of Human Quantitative Traits in Related Individuals , 2013, PloS one.

[39]  J. Seddon,et al.  Evolutionary history of DLA class II haplotypes in canine diabetes mellitus through single nucleotide polymorphism genotyping. , 2010, Tissue antigens.

[40]  K. Lindblad-Toh,et al.  Genome-Wide Analysis in German Shepherd Dogs Reveals Association of a Locus on CFA 27 with Atopic Dermatitis , 2013, PLoS genetics.

[41]  C. Cremers,et al.  Mutations in the lipoma HMGIC fusion partner‐like 5 (LHFPL5) gene cause autosomal recessive nonsyndromic hearing loss , 2006, Human mutation.

[42]  C. Pineau,et al.  Tat1, a Novel Sulfate Transporter Specifically Expressed in Human Male Germ Cells and Potentially Linked to RhoGTPase Signaling* , 2001, The Journal of Biological Chemistry.

[43]  P. Laurberg,et al.  Iodine intake as a determinant of thyroid disorders in populations. , 2010, Best practice & research. Clinical endocrinology & metabolism.

[44]  Ross Ihaka,et al.  Gentleman R: R: A language for data analysis and graphics , 1996 .

[45]  B. Everitt,et al.  Statistical methods for rates and proportions , 1973 .

[46]  Albert J. Vilella,et al.  A high-resolution map of human evolutionary constraint using 29 mammals , 2011, Nature.

[47]  P. Anderson,et al.  Human Autoimmune Sera as Molecular Probes for the Identification of an Autoantigen Kinase Signaling Pathway , 2002, The Journal of experimental medicine.

[48]  E. Mauceli,et al.  A Novel Unstable Duplication Upstream of HAS2 Predisposes to a Breed-Defining Skin Phenotype and a Periodic Fever Syndrome in Chinese Shar-Pei Dogs , 2011, PLoS genetics.

[49]  R. Nachreiner,et al.  Lymphocytic thyroiditis. , 2020, The Veterinary clinics of North America. Small animal practice.

[50]  S. Wallenstein,et al.  Genetics of autoimmune thyroid disease: lack of evidence for linkage to HLA within families. , 1992, The Journal of clinical endocrinology and metabolism.

[51]  Gabriel E. Hoffman,et al.  Correcting for Population Structure and Kinship Using the Linear Mixed Model: Theory and Extensions , 2013, PloS one.

[52]  K. Lindblad-Toh,et al.  Efficient mapping of mendelian traits in dogs through genome-wide association , 2007, Nature Genetics.

[53]  R. Osman,et al.  Molecular amino acid signatures in the MHC class II peptide-binding pocket predispose to autoimmune thyroiditis in humans and in mice , 2008, Proceedings of the National Academy of Sciences.

[54]  A. Karlsson,et al.  Prevalence of diagnostic characteristics indicating canine autoimmune lymphocytic thyroiditis in giant schnauzer and hovawart dogs. , 2009, The Journal of small animal practice.

[55]  Y. Tomer,et al.  Analysis of HLA genes in families with autoimmune diabetes and thyroiditis. , 2004, Human immunology.

[56]  N. Farid,et al.  Hashimoto's thyroiditis is associated with HLA-DRw3. , 1978, The New England journal of medicine.

[57]  K. Lindblad-Toh,et al.  Identification of Genomic Regions Associated with Phenotypic Variation between Dog Breeds using Selection Mapping , 2011, PLoS genetics.

[58]  S J Gange,et al.  Epidemiology and estimated population burden of selected autoimmune diseases in the United States. , 1997, Clinical immunology and immunopathology.

[59]  E. Prescott,et al.  Human Papillomavirus Type 1 E1^E4 Protein Is a Potent Inhibitor of the Serine-Arginine (SR) Protein Kinase SRPK1 and Inhibits Phosphorylation of Host SR Proteins and of the Viral Transcription and Replication Regulator E2 , 2014, Journal of Virology.

[60]  W. Grove.  Statistical Methods for Rates and Proportions, 2nd ed , 1981 .

[61]  K. Lindblad-Toh,et al.  The genomic signature of dog domestication reveals adaptation to a starch-rich diet , 2013, Nature.

[62]  K. Lindblad-Toh,et al.  Localization of Canine Brachycephaly Using an Across Breed Mapping Approach , 2010, PloS one.

[63]  Amiel A. Dror,et al.  Atrophic thyroid follicles and inner ear defects reminiscent of cochlear hypothyroidism in Slc26a4-related deafness , 2014, Mammalian Genome.

[64]  P. Tufféry,et al.  Missense mutations in SLC26A8, encoding a sperm-specific activator of CFTR, are associated with human asthenozoospermia. , 2013, American journal of human genetics.

[65]  Stephen D. Turner.  qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots , 2014 .

[66]  K. Pearson.  On the Criterion that a Given System of Deviations from the Probable in the Case of a Correlated System of Variables is Such that it Can be Reasonably Supposed to have Arisen from Random Sampling , 1900 .

[67]  J. Ioannidis,et al.  Meta-analysis methods for genome-wide association studies and beyond , 2013, Nature Reviews Genetics.

[68]  W. Ollier,et al.  Association of canine hypothyroidism with a common major histocompatibility complex DLA class II allele. , 2006, Tissue antigens.

[69]  D. Scott,et al.  A survey of canine and feline skin disorders seen in a university practice: Small Animal Clinic, University of Montréal, Saint-Hyacinthe, Québec (1987-1988). , 1990, The Canadian veterinary journal = La revue veterinaire canadienne.

[70]  Paul Scheet,et al.  A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. , 2006, American journal of human genetics.

[71]  E. Lundberg,et al.  A Genecentric Human Protein Atlas for Expression Profiles Based on Antibodies* , 2008, Molecular & Cellular Proteomics.

[72]  Eric S. Lander,et al.  An Improved Canine Genome and a Comprehensive Catalogue of Coding Genes and Non-Coding Transcripts , 2014, PloS one.

[73]  Yurii S. Aulchenko,et al.  BIOINFORMATICS APPLICATIONS NOTE doi:10.1093/bioinformatics/btm108 Genetics and population analysis GenABEL: an R library for genome-wide association analysis , 2022 .