Direct-to-consumer DNA testing of 6,000 dogs reveals 98.6-kb duplication associated with blue eyes and heterochromia in Siberian Huskies

Common variants in the CPT1A gene are associated with cataracts in Northern breeds of domestic dog

Primary hereditary cataract affects many purebred domestic dog breeds and is a major cause of visual impairment in dogs. Cataracts are common in Northern breeds such as the Siberian Husky, Alaskan Malamute and Samoyed, but their aetiology is currently unknown. Only two genetic loci are known to be causally related to primary hereditary cataracts in the dog. To search for genetic loci associated with cataracts in Northern breeds, we used a genome-wide association study approach in three breeds-Siberian Husky, Alaskan Malamute and Samoyed. Cases were defined as dogs with bilateral posterior polar subcapsular cataracts and controls were at least four years of age with no evidence of cataracts or other ocular abnormality. We found a genome-wide statistical association for cataracts in the Siberian Husky on canine chromosome 18 (P-value: 1.1 x 10 - 7), which was independently replicated in a second larger case-control set (P-value 9.8 x 10 - 29). The Samoyed breed also showed evidence for association in the same genomic region (P-value: 2.4 x 10 - 5). We subsequently used targeted resequencing of the associated region (6.5 Mb) in ten Siberian Huskies and whole genome sequencing of a Husky, Malamute, Samoyed and Norwegian Buhund case to conduct fine-mapping and screen for candidate causal variants. These analyses identified a region of linkage disequilibrium in the four breeds containing common variants in the carnitine palmitoyltransferase 1A (CPT1A) gene that are strongly associated with bilateral posterior polar subcapsular cataracts in the Siberian Husky, Samoyed, Icelandic Sheepdog and Norwegian Buhund and we demonstrate that CPT1A is expressed in the dog lens and retina through RNAseq. Our findings represent a novel locus for cataracts in dogs. However, further work is needed to elucidate the pathophysiology underlying the association between CPT1A and cataracts in Northern breeds.

Is There Hybrid Vigor in Dogs? Oxidative Stress and Cytokine Concentrations in Low- To Mid-Content Wolf-Dog Hybrids

AbstractChanges of deleterious traits in mitochondria within hybrids of two different species are likely to be reflected across whole-animal phenotypes. Two processes linked to mitochondrial fitness are oxidative stress and inflammation. Here, plasma oxidative stress (lipid oxidative damage, total antioxidant capacity, and activities of catalase, glutathione peroxidase, and superoxide dismutase concentration) and cytokine concentrations (IL-1β, IL-6, and TNF-α) from wolf-dog hybrids of low to mid wolf content are presented and then compared with previously published values for similar-sized dogs and gray wolves. Results indicate that lipid oxidative damage and catalase activity were not significantly different across species and hybrids. Total antioxidant capacity and glutathione peroxidase were significantly lower in wolf-dog hybrids than in domestic dogs and wolves (although total antioxidant capacity significantly increased with wolf content), but superoxide dismutase concentration was higher in hybrids. Thus, it seems that a low percentage of wolf content decreases several aspects of antioxidants but without any accumulating lipid oxidative damage. Additionally, wolf-dog hybrids had higher IL-1β and IL-6 concentrations but lower TNF-α concentrations than domestic dogs. And there was a significantly positive correlation between percentage of wolf content and IL-1β. These data imply either a mitonuclear incompatibility or a nuclear-nuclear incompatibility within wolf-dog hybrids, a case that does not support heterosis.

The LASSIE MPS panel: Predicting externally visible traits in dogs for forensic purposes

Predicting the outward appearance of dogs via their DNA, also known as Canine DNA Phenotyping, is a young, emerging field of research in forensic genetics. The few previous studies published in this respect were restricted to the consecutive analysis of single DNA markers, a process that is time- and sample-consuming and therefore not a viable option for limited forensic specimens. Here, we report on the development and evaluation of a Massively Parallel Sequencing (MPS) based molecular genetic assay, the LASSIE MPS Panel. This panel aims to predict externally visible as well as skeletal traits, which include coat color, coat pattern, coat structure, tail morphology, skull shape, ear shape, eye color and body size from DNA using 44 genetic markers in a single molecular genetic assay. A biostatistical naïve Bayes classification approach was applied to identify the most informative marker combinations for predicting phenotypes. Overall, the predictive performance was characterized by a very high classification success for some of the trait categories, and high to moderate success for others. The performance of the developed predictive framework was further evaluated using blind samples from three randomly selected dog individuals, whose appearance was well predicted.

​Comparative genomics of Balto, a famous historic dog, captures lost diversity of 1920s sled dogs

We reconstruct the phenotype of Balto, the heroic sled dog renowned for transporting diphtheria antitoxin to Nome, Alaska, in 1925, using evolutionary constraint estimates from the Zoonomia alignment of 240 mammals and 682 genomes from dogs and wolves of the 21st century. Balto shares just part of his diverse ancestry with the eponymous Siberian husky breed. Balto's genotype predicts a combination of coat features atypical for modern sled dog breeds, and a slightly smaller stature. He had enhanced starch digestion compared with Greenland sled dogs and a compendium of derived homozygous coding variants at constrained positions in genes connected to bone and skin development. We propose that Balto's population of origin, which was less inbred and genetically healthier than that of modern breeds, was adapted to the extreme environment of 1920s Alaska.

A Study Evaluating Consumer Motivations, Perceptions, and Responses to Direct-to-Consumer Canine Genetic Test Results

Direct-to-consumer genetic services allow companion animal guardians to purchase a DNA test and receive detailed results about their pet's ancestry, health, and traits results. In collaboration with Wisdom Panel, we present novel findings about consumer motivations, perceptions, and responses to their use of canine genomic services. Wisdom Panel customers were invited to complete an online survey anonymously in which they were asked about their reasons for using a genetic test for their dog, how they perceived the test's results, and how they responded to the results they received. Participant data revealed most utilized a test that provided more ancestry/breed results (75.9%) as compared to health-related results. The majority of participants perceived the breed test results as accurate (52.0% strongly agree, 27.6% somewhat agree) and the genetic services provided as having great value (49.6% strongly agree, 32.7% somewhat agreed). In responding to their dog's results, participants indicated they shared the information with family (88.1%) and friends (84.2%). Collectively, our study indicates consumers are more focused on their dog's ancestry than other test results. Using these findings and previous literature on human direct-to-consumer genetic testing, human-animal dyads, and identity construction, we consider the possibility of "breed options theory" and future areas of research.

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.

Grey wolf genomic history reveals a dual ancestry of dogs

The grey wolf (Canis lupus) was the first species to give rise to a domestic population, and they remained widespread throughout the last Ice Age when many other large mammal species went extinct. Little is known, however, about the history and possible extinction of past wolf populations or when and where the wolf progenitors of the present-day dog lineage (Canis familiaris) lived^1–8. Here we analysed 72 ancient wolf genomes spanning the last 100,000 years from Europe, Siberia and North America. We found that wolf populations were highly connected throughout the Late Pleistocene, with levels of differentiation an order of magnitude lower than they are today. This population connectivity allowed us to detect natural selection across the time series, including rapid fixation of mutations in the gene IFT88 40,000–30,000 years ago. We show that dogs are overall more closely related to ancient wolves from eastern Eurasia than to those from western Eurasia, suggesting a domestication process in the east. However, we also found that dogs in the Near East and Africa derive up to half of their ancestry from a distinct population related to modern southwest Eurasian wolves, reflecting either an independent domestication process or admixture from local wolves. None of the analysed ancient wolf genomes is a direct match for either of these dog ancestries, meaning that the exact progenitor populations remain to be located.

DNA from ancient wolves spanning 100,000 years sheds light on wolves’ evolutionary history and the genomic origin of dogs.

Ancestry-inclusive dog genomics challenges popular breed stereotypes

Behavioral genetics in dogs has focused on modern breeds, which are isolated subgroups with distinctive physical and, purportedly, behavioral characteristics. We interrogated breed stereotypes by surveying owners of 18,385 purebred and mixed-breed dogs and genotyping 2155 dogs. Most behavioral traits are heritable [heritability (h²) > 25%], and admixture patterns in mixed-breed dogs reveal breed propensities. Breed explains just 9% of behavioral variation in individuals. Genome-wide association analyses identify 11 loci that are significantly associated with behavior, and characteristic breed behaviors exhibit genetic complexity. Behavioral loci are not unusually differentiated in breeds, but breed propensities align, albeit weakly, with ancestral function. We propose that behaviors perceived as characteristic of modern breeds derive from thousands of years of polygenic adaptation that predates breed formation, with modern breeds distinguished primarily by aesthetic traits.

Five genetic variants explain over 70% of hair coat pheomelanin intensity variation in purebred and mixed breed domestic dogs

In mammals, the pigment molecule pheomelanin confers red and yellow color to hair, and the intensity of this coloration is caused by variation in the amount of pheomelanin. Domestic dogs exhibit a wide range of pheomelanin intensity, ranging from the white coat of the Samoyed to the deep red coat of the Irish Setter. While several genetic variants have been associated with specific coat intensity phenotypes in certain dog breeds, they do not explain the majority of phenotypic variation across breeds. In order to gain further insight into the extent of multigenicity and epistatic interactions underlying coat pheomelanin intensity in dogs, we leveraged a large dataset obtained via a direct-to-consumer canine genetic testing service. This consisted of genome-wide single nucleotide polymorphism (SNP) genotype data and owner-provided photos for 3,057 pheomelanic mixed breed and purebred dogs from 63 breeds and varieties spanning the full range of canine coat pheomelanin intensity. We first performed a genome-wide association study (GWAS) on 2,149 of these dogs to search for additional genetic variants that underlie intensity variation. GWAS identified five loci significantly associated with intensity, of which two (CFA15 29.8 Mb and CFA20 55.8 Mb) replicate previous findings and three (CFA2 74.7 Mb, CFA18 12.9 Mb, CFA21 10.9 Mb) have not previously been reported. In order to assess the combined predictive power of these loci across dog breeds, we used our GWAS data set to fit a linear model, which explained over 70% of variation in coat pheomelanin intensity in an independent validation dataset of 908 dogs. These results introduce three novel pheomelanin intensity loci, and further demonstrate the multigenic nature of coat pheomelanin intensity determination in domestic dogs.

Man's best friend in life and death: scientific perspectives and challenges of dog brain banking

Biobanking refers to the systematic collection, storage, and distribution of pre- or post-mortem biological samples derived from volunteer donors. The demand for high-quality human specimens is clearly demonstrated by the number of newly emerging biobanking facilities and large international collaborative networks. Several animal species are relevant today in medical research; therefore, similar initiatives in comparative physiology could be fruitful. Dogs, in particular, are gaining increasing attention in translational research on complex phenomena, like aging, cancer, and neurodegenerative diseases. Therefore, biobanks gathering and storing dog biological materials together with related data could play a vital role in translational and veterinary research projects. To achieve these aims, a canine biobank should meet the same standards in sample quality and data management as human biobanks and should rely on well-designed collaborative networks between different professionals and dog owners. While efforts to create dog biobanks could face similar financial and technical challenges as their human counterparts, they can widen the spectrum of successful collaborative initiatives towards a better picture of dogs’ physiology, disease, evolution, and translational potential. In this review, we provide an overview about the current state of dog biobanking and introduce the “Canine Brain and Tissue Bank” (CBTB)—a new, large-scale collaborative endeavor in the field.

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.

Identification of High-Confidence Structural Variants in Domesticated Rainbow Trout Using Whole-Genome Sequencing

Genomic structural variants (SVs) are a major source of genetic and phenotypic variation but have not been investigated systematically in rainbow trout (Oncorhynchus mykiss), an important aquaculture species of cold freshwater. The objectives of this study were 1) to identify and validate high-confidence SVs in rainbow trout using whole-genome re-sequencing; and 2) to examine the contribution of transposable elements (TEs) to SVs in rainbow trout. A total of 96 rainbow trout, including 11 homozygous lines and 85 outbred fish from three breeding populations, were whole-genome sequenced with an average genome coverage of 17.2×. Putative SVs were identified using the program Smoove which integrates LUMPY and other associated tools into one package. After rigorous filtering, 13,863 high-confidence SVs were identified. Pacific Biosciences long-reads of Arlee, one of the homozygous lines used for SV detection, validated 98% (3,948 of 4,030) of the high-confidence SVs identified in the Arlee homozygous line. Based on principal component analysis, the 85 outbred fish clustered into three groups consistent with their populations of origin, further indicating that the high-confidence SVs identified in this study are robust. The repetitive DNA content of the high-confidence SV sequences was 86.5%, which is much higher than the 57.1% repetitive DNA content of the reference genome, and is also higher than the repetitive DNA content of Atlantic salmon SVs reported previously. TEs thus contribute substantially to SVs in rainbow trout as TEs make up the majority of repetitive sequences. Hundreds of the high-confidence SVs were annotated as exon-loss or gene-fusion variants, and may have phenotypic effects. The high-confidence SVs reported in this study provide a foundation for further rainbow trout SV studies.

Genome-wide association studies for iris pigmentation and heterochromia patterns in Large White pigs

Eye colour genetics have been extensively studied in humans since the rediscovery of Mendel's laws. This trait was first interpreted using simplistic genetic models but soon it was realised that it is more complex. In this study, we analysed eye colour variability in a Large White pig population (n = 897) and report the results of GWASs based on several comparisons including pigs having four main eye colour categories (three with both pigmented eyes of different brown grades: pale, 17.9%; medium, 14.8%; and dark, 54.3%; another one with both eyes completely depigmented, 3.8%) and heterochromia patterns (heterochromia iridis - depigmented iris sectors in pigmented irises, 3.2%; heterochromia iridum - one whole eye iris of depigmented phenotype and the other eye with the iris completely pigmented, 5.9%). Pigs were genotyped with the Illumina PorcineSNP60 BeadChip and GEMMA was used for the association analyses. The results indicated that SLC45A2 (on chromosome 16, SSC16), EDNRB (SSC11) and KITLG (SSC5) affect the different grades of brown pigmentation of the eyes, the bilateral eye depigmentation defect and the heterochromia iridis defect recorded in this white pig population respectively. These genes are involved in several mechanisms affecting pigmentation. Significant associations for the eye depigmented patterns were also identified for SNPs on two SSC4 regions (including two candidate genes: NOTCH2 and PREX2) and on SSC6, SSC8 and SSC14 (including COL17A1 as candidate gene). This study provided useful information to understand eye pigmentation mechanisms, further valuing the pig as animal model to study complex phenotypes in humans.

Merle allele variations in the Mudi dog breed and their effects on phenotypes

A retrotransposon insertion in the SILV gene is associated with a peculiar phenotype of dog, known as a merle. It is characterised by various areas of their coat colour becoming diluted due to a malfunction in the eumelanin-producing pigment cells. Recent studies have shown that the exact size of the short interspersed element (SINE) insertion is in correlation with specific phenotypic attributes, but was not able to absolutely confine dogs to a certain colour pattern. Our study focused on the merle variations occurring in the Mudi breed. Altogether, 123 dog samples from 11 countries were tested and genotyped. The exact length of the merle alleles were determined by automated fluorescent capillary fragment analysis. The most frequent merle genotype in this Mudi sample collection was the 'classic' merle (m/M: 61.8%), whereas other variants, such as atypical (m/Ma and m/Ma+: 5.7%), harlequin (m/Mh: 13.8%), double merle (M/M: 0.8%) and mosaic profiles (17.9%) were also observed. The practical significance of testing this mutation is that, phenotypically, not only merle dogs are carriers of this insertion, but also the so-called hidden merle individuals (where the merle phenotype is fully covered by the pheomelanin-dominated colouration) are potentially capable of producing unintentionally homozygous 'double merle' progeny with ophthalmologic, viability and auditory impairments.

Inbreeding depression causes reduced fecundity in Golden Retrievers

Inbreeding depression has been demonstrated to impact vital rates, productivity, and performance in human populations, wild and endangered species, and in recent years, the domestic species. In all cases, standardized, high-quality phenotype data on all individuals are invaluable for longitudinal analyses such as those required to evaluate vital rates of a study cohort. Further, many investigators agree upon the preference for and utility of genomic measures of inbreeding in lieu of pedigree-based estimates of inbreeding. We evaluated the association of measures of reproductive fitness in 93 Golden Retrievers enrolled in the Golden Retriever Lifetime Study with a genomic measurement of inbreeding, F_ROH. We demonstrate a statistically significant negative correlation between fecundity and F_ROH. This work sets the stage for larger scale analyses to investigate genomic regions associated with fecundity and other measures of fitness.

Copy Number Variation in Domestication

Domesticated plants have long served as excellent models for studying evolution. Many genes and mutations underlying important domestication traits have been identified, and most causal mutations appear to be SNPs. Copy number variation (CNV) is an important source of genetic variation that has been largely neglected in studies of domestication. Ongoing work demonstrates the importance of CNVs as a source of genetic variation during domestication, and during the diversification of domesticated taxa. Here, we review how CNVs contribute to evolutionary processes underlying domestication, and review examples of domestication traits caused by CNVs. We draw from examples in plant species, but also highlight cases in animal systems that could illuminate the roles of CNVs in the domestication process.