Genome-wide analyses reveals an association between invasive urothelial carcinoma in the Shetland sheepdog and NIPAL1.

Naturally occurring canine invasive urinary carcinoma (iUC) closely resembles human muscle invasive bladder cancer in terms of histopathology, metastases, response to therapy, and low survival rate. The heterogeneous nature of the disease has led to the association of large numbers of risk loci in humans, however most are of small effect. There exists a need for new and accurate animal models of invasive bladder cancer. In dogs, distinct breeds show markedly different rates of iUC, thus presenting an opportunity to identify additional risk factors and overcome the locus heterogeneity encountered in human mapping studies. In the association study presented here, inclusive of 100 Shetland sheepdogs and 58 dogs of other breeds, we identify a homozygous protein altering point mutation within the NIPAL1 gene which increases risk by eight-fold (OR = 8.42, CI = 3.12-22.71), accounting for nearly 30% of iUC risk in the Shetland sheepdog. Inclusion of six additional loci accounts for most of the disease risk in the breed and explains nearly 75% of the phenotypes in this study. When combined with sequence data from tumors, we show that variation in the MAPK signaling pathway is an overarching cause of iUC susceptibility in dogs.

DNA methylation clocks for dogs and humans.

DNA methylation profiles have been used to develop biomarkers of aging known as epigenetic clocks, which predict chronological age with remarkable accuracy and show promise for inferring health status as an indicator of biological age. Epigenetic clocks were first built to monitor human aging, but their underlying principles appear to be evolutionarily conserved, as they have now been successfully developed for many mammalian species. Here, we describe reliable and highly accurate epigenetic clocks shown to apply to 93 domestic dog breeds. The methylation profiles were generated using the mammalian methylation array, which utilizes DNA sequences that are conserved across all mammalian species. Canine epigenetic clocks were constructed to estimate age and also average time to death. We also present two highly accurate human­dog dual species epigenetic clocks (R = 0.97), which may facilitate the ready translation from canine to human use (or vice versa) of antiaging treatments being developed for longevity and preventive medicine. Finally, epigenome-wide association studies here reveal individual methylation sites that may underlie the inverse relationship between breed weight and lifespan. Overall, we describe robust biomarkers to measure aging and, potentially, health status in canines.

Natural and human-driven selection of a single non-coding body size variant in ancient and modern canids.

Domestic dogs (Canis lupus familiaris) are the most variable-sized mammalian species on Earth, displaying a 40-fold size difference between breeds.1 Although dogs of variable size are found in the archeological record,2-4 the most dramatic shifts in body size are the result of selection over the last two centuries, as dog breeders selected and propagated phenotypic extremes within closed breeding populations.5 Analyses of over 200 domestic breeds have identified approximately 20 body size genes regulating insulin processing, fatty acid metabolism, TGFß signaling, and skeletal formation.6-10 Of these, insulin-like growth factor 1 (IGF1) predominates, controlling approximately 15% of body size variation between breeds.8 The identification of a functional mutation associated with IGF1 has thus far proven elusive.6,10,11 Here, to identify and elucidate the role of an ancestral IGF1 allele in the propagation of modern canids, we analyzed 1,431 genome sequences from 13 species, including both ancient and modern canids, thus allowing us to define the evolutionary history of both ancestral and derived alleles at this locus. We identified a single variant in an antisense long non-coding RNA (IGF1-AS) that interacts with the IGF1 gene, creating a duplex. While the derived mutation predominates in both modern gray wolves and large domestic breeds, the ancestral allele, which predisposes to small size, was common in small-sized breeds and smaller wild canids. Our analyses demonstrate that this major regulator of canid body size nearly vanished in Pleistocene wolves, before its recent resurgence resulting from human-imposed selection for small-sized breed dogs.

Multi-omics approach identifies germline regulatory variants associated with hematopoietic malignancies in retriever dog breeds.

Histiocytic sarcoma is an aggressive hematopoietic malignancy of mature tissue histiocytes with a poorly understood etiology in humans. A histologically and clinically similar counterpart affects flat-coated retrievers (FCRs) at unusually high frequency, with 20% developing the lethal disease. The similar clinical presentation combined with the closed population structure of dogs, leading to high genetic homogeneity, makes dogs an excellent model for genetic studies of cancer susceptibility. To determine the genetic risk factors underlying histiocytic sarcoma in FCRs, we conducted multiple genome-wide association studies (GWASs), identifying two loci that confer significant risk on canine chromosomes (CFA) 5 (Pwald = 4.83x10-9) and 19 (Pwald = 2.25x10-7). We subsequently undertook a multi-omics approach that has been largely unexplored in the canine model to interrogate these regions, generating whole genome, transcriptome, and chromatin immunoprecipitation sequencing. These data highlight the PI3K pathway gene PIK3R6 on CFA5, and proximal candidate regulatory variants that are strongly associated with histiocytic sarcoma and predicted to impact transcription factor binding. The CFA5 association colocalizes with susceptibility loci for two hematopoietic malignancies, hemangiosarcoma and B-cell lymphoma, in the closely related golden retriever breed, revealing the risk contribution this single locus makes to multiple hematological cancers. By comparison, the CFA19 locus is unique to the FCR and harbors risk alleles associated with upregulation of TNFAIP6, which itself affects cell migration and metastasis. Together, these loci explain ~35% of disease risk, an exceptionally high value that demonstrates the advantages of domestic dogs for complex trait mapping and genetic studies of cancer susceptibility.

Copy number variation underlies complex phenotypes in domestic dog breeds and other canids.

Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3, which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds.

Genetic analysis of the modern Australian labradoodle dog breed reveals an excess of the poodle genome.

The genomic diversity of the domestic dog is an invaluable resource for advancing understanding of mammalian biology, evolutionary biology, morphologic variation, and behavior. There are approximately 350 recognized breeds in the world today, many established through hybridization and selection followed by intense breeding programs aimed at retaining or enhancing specific traits. As a result, many breeds suffer from an excess of particular diseases, one of many factors leading to the recent trend of "designer breed" development, i.e. crossing purebred dogs from existing breeds in the hope that offspring will be enriched for desired traits and characteristics of the parental breeds. We used a dense panel of 150,106 SNPs to analyze the population structure of the Australian labradoodle (ALBD), to understand how such breeds are developed. Haplotype and admixture analyses show that breeds other than the poodle (POOD) and Labrador retriever (LAB) contributed to ALBD formation, but that the breed is, at the genetic level, predominantly POOD, with all small and large varieties contributing to its construction. Allele frequency analysis reveals that the breed is enhanced for variants associated with a poodle-like coat, which is perceived by breeders to have an association with hypoallergenicity. We observed little enhancement for LAB-specific alleles. This study provides a blueprint for understanding how dog breeds are formed, highlighting the limited scope of desired traits in defining new breeds.

Canine neuropathies: powerful spontaneous models for human hereditary sensory neuropathies.

In humans, hereditary sensory neuropathies (HSN), also known as hereditary sensory and autonomic neuropathies (HSAN), constitute a clinically and genetically heterogeneous group of disorders characterized by progressive sensory loss, often accompanied by chronic skin ulcerations and nail dystrophic changes. To date, although around 20 genes have already been discovered, they do not explain the genetic causes of all patients. In dogs, similar neuropathies are also diagnosed, several breeds being predisposed to specific forms of the disease. Indeed, the breed specificity of most canine genetic diseases is due to the small numbers of founders and high levels of inbreeding. Recent knowledge and tools developed to study the canine genome efficiently allows deciphering the genetic bases of such diseases. To date, a dozen breeds are recognized to develop specific HSN. For the Border collie and hunting dog breeds, the genes involved have recently been discovered. Other affected breeds thus constitute potential genetic models, with new genes to be found in dogs that can be considered as candidate genes for human HSAN/HSN. Here, we review the different forms of human and canine HSAN/HSN and we present a novel form in Fox terrier cases, highlighting the advantages of the dog model for such rare human diseases.

Whole genome sequencing of canids reveals genomic regions under selection and variants influencing morphology.

Domestic dog breeds are characterized by an unrivaled diversity of morphologic traits and breed-associated behaviors resulting from human selective pressures. To identify the genetic underpinnings of such traits, we analyze 722 canine whole genome sequences (WGS), documenting over 91 million single nucleotide and small indels, creating a large catalog of genomic variation for a companion animal species. We undertake both selective sweep analyses and genome wide association studies (GWAS) inclusive of over 144 modern breeds, 54 wild canids and a hundred village dogs. Our results identify variants of strong impact associated with 16 phenotypes, including body weight variation which, when combined with existing data, explain greater than 90% of body size variation in dogs. We thus demonstrate that GWAS and selection scans performed with WGS are powerful complementary methods for expanding the utility of companion animal systems for the study of mammalian growth and biology.

Genetic selection of athletic success in sport-hunting dogs.

Modern dogs are distinguished among domesticated species by the vast breadth of phenotypic variation produced by strong and consistent human-driven selective pressure. The resulting breeds reflect the development of closed populations with well-defined physical and behavioral attributes. The sport-hunting dog group has long been employed in assistance to hunters, reflecting strong behavioral pressures to locate and pursue quarry over great distances and variable terrain. Comparison of whole-genome sequence data between sport-hunting and terrier breeds, groups at the ends of a continuum in both form and function, reveals that genes underlying cardiovascular, muscular, and neuronal functions are under strong selection in sport-hunting breeds, including ADRB1, TRPM3, RYR3, UTRN, ASIC3, and ROBO1 We also identified an allele of TRPM3 that was significantly associated with increased racing speed in Whippets, accounting for 11.6% of the total variance in racing performance. Finally, we observed a significant association of ROBO1 with breed-specific accomplishments in competitive obstacle course events. These results provide strong evidence that sport-hunting breeds have been adapted to their occupations by improved endurance, cardiac function, blood flow, and cognitive performance, demonstrating how strong behavioral selection alters physiology to create breeds with distinct capabilities.

Analysis of large versus small dogs reveals three genes on the canine X chromosome associated with body weight, muscling and back fat thickness.

Domestic dog breeds display significant diversity in both body mass and skeletal size, resulting from intensive selective pressure during the formation and maintenance of modern breeds. While previous studies focused on the identification of alleles that contribute to small skeletal size, little is known about the underlying genetics controlling large size. We first performed a genome-wide association study (GWAS) using the Illumina Canine HD 170,000 single nucleotide polymorphism (SNP) array which compared 165 large-breed dogs from 19 breeds (defined as having a Standard Breed Weight (SBW) >41 kg [90 lb]) to 690 dogs from 69 small breeds (SBW ≤41 kg). We identified two loci on the canine X chromosome that were strongly associated with large body size at 82-84 megabases (Mb) and 101-104 Mb. Analyses of whole genome sequencing (WGS) data from 163 dogs revealed two indels in the Insulin Receptor Substrate 4 (IRS4) gene at 82.2 Mb and two additional mutations, one SNP and one deletion of a single codon, in Immunoglobulin Superfamily member 1 gene (IGSF1) at 102.3 Mb. IRS4 and IGSF1 are members of the GH/IGF1 and thyroid pathways whose roles include determination of body size. We also found one highly associated SNP in the 5'UTR of Acyl-CoA Synthetase Long-chain family member 4 (ACSL4) at 82.9 Mb, a gene which controls the traits of muscling and back fat thickness. We show by analysis of sequencing data from 26 wolves and 959 dogs representing 102 domestic dog breeds that skeletal size and body mass in large dog breeds are strongly associated with variants within IRS4, ACSL4 and IGSF1.

A Point Mutation in a lincRNA Upstream of GDNF Is Associated to a Canine Insensitivity to Pain: A Spontaneous Model for Human Sensory Neuropathies.

Human Hereditary Sensory Autonomic Neuropathies (HSANs) are characterized by insensitivity to pain, sometimes combined with self-mutilation. Strikingly, several sporting dog breeds are particularly affected by such neuropathies. Clinical signs appear in young puppies and consist of acral analgesia, with or without sudden intense licking, biting and severe self-mutilation of the feet, whereas proprioception, motor abilities and spinal reflexes remain intact. Through a Genome Wide Association Study (GWAS) with 24 affected and 30 unaffected sporting dogs using the Canine HD 170K SNP array (Illumina), we identified a 1.8 Mb homozygous locus on canine chromosome 4 (adj. p-val = 2.5x10-6). Targeted high-throughput sequencing of this locus in 4 affected and 4 unaffected dogs identified 478 variants. Only one variant perfectly segregated with the expected recessive inheritance in 300 sporting dogs of known clinical status, while it was never present in 900 unaffected dogs from 130 other breeds. This variant, located 90 kb upstream of the GDNF gene, a highly relevant neurotrophic factor candidate gene, lies in a long intergenic non-coding RNAs (lincRNA), GDNF-AS. Using human comparative genomic analysis, we observed that the canine variant maps onto an enhancer element. Quantitative RT-PCR of dorsal root ganglia RNAs of affected dogs showed a significant decrease of both GDNF mRNA and GDNF-AS expression levels (respectively 60% and 80%), as compared to unaffected dogs. We thus performed gel shift assays (EMSA) that reveal that the canine variant significantly alters the binding of regulatory elements. Altogether, these results allowed the identification in dogs of GDNF as a relevant candidate for human HSAN and insensitivity to pain, but also shed light on the regulation of GDNF transcription. Finally, such results allow proposing these sporting dog breeds as natural models for clinical trials with a double benefit for human and veterinary medicine.

A mutation in the FAM83G gene in dogs with hereditary footpad hyperkeratosis (HFH).

Hereditary footpad hyperkeratosis (HFH) represents a palmoplantar hyperkeratosis, which is inherited as a monogenic autosomal recessive trait in several dog breeds, such as e.g. Kromfohrländer and Irish Terriers. We performed genome-wide association studies (GWAS) in both breeds. In Kromfohrländer we obtained a single strong association signal on chromosome 5 (p(raw) = 1.0×10(-13)) using 13 HFH cases and 29 controls. The association signal replicated in an independent cohort of Irish Terriers with 10 cases and 21 controls (p(raw) = 6.9×10(-10)). The analysis of shared haplotypes among the combined Kromfohrländer and Irish Terrier cases defined a critical interval of 611 kb with 13 predicted genes. We re-sequenced the genome of one affected Kromfohrländer at 23.5× coverage. The comparison of the sequence data with 46 genomes of non-affected dogs from other breeds revealed a single private non-synonymous variant in the critical interval with respect to the reference genome assembly. The variant is a missense variant (c.155G>C) in the FAM83G gene encoding a protein with largely unknown function. It is predicted to change an evolutionary conserved arginine into a proline residue (p.R52P). We genotyped this variant in a larger cohort of dogs and found perfect association with the HFH phenotype. We further studied the clinical and histopathological alterations in the epidermis in vivo. Affected dogs show a moderate to severe orthokeratotic hyperplasia of the palmoplantar epidermis. Thus, our data provide the first evidence that FAM83G has an essential role for maintaining the integrity of the palmoplantar epidermis.