Leader of the pack: gene mapping in dogs and other model organisms

A variant in the 5'UTR of ERBB4 is associated with lifespan in Golden Retrievers

Genome-wide association studies (GWAS) in long-lived human populations have led to identification of variants associated with Alzheimer's disease and cardiovascular disease, the latter being the most common cause of mortality in people worldwide. In contrast, naturally occurring cancer represents the leading cause of death in pet dogs, and specific breeds like the Golden Retriever (GR) carry up to a 65% cancer-related death rate. We hypothesized that GWAS of long-lived GRs might lead to the identification of genetic variants capable of modifying longevity within this cancer-predisposed breed. A GWAS was performed comparing GR dogs ≥ 14 years to dogs dying prior to age 12 which revealed a significant association to ERBB4, the only member of the epidermal growth factor receptor family capable of serving as both a tumor suppressor gene and an oncogene. No coding variants were identified, however, distinct haplotypes in the 5'UTR were associated with reduced lifespan in two separate populations of GR dogs. When all GR dogs were analyzed together (n = 304), the presence of haplotype 3 was associated with shorter survival (11.8 years vs. 12.8 years, p = 0.024). GRs homozygous for haplotype 3 had the shortest survival, and GRs homozygous for haplotype 1 had the longest survival (11.6 years vs. 13.5 years, p = 0.0008). Sub-analyses revealed that the difference in lifespan for GRs carrying at least 1 copy of haplotype 3 was specific to female dogs (p = 0.009), whereas survival remained significantly different in both male and female GRs homozygous for haplotype 1 or haplotype 3 (p = 0.026 and p = 0.009, respectively). Taken together, these findings implicate a potential role for ERBB4 in GR longevity and provide evidence that within-breed canine lifespan studies could serve as a mechanism to identify favorable or disease-modifying variants important to the axis of aging and cancer.

Genome-wide association study suggests genetic candidate loci of insulin dysregulation in Finnhorses

Equine metabolic syndrome (EMS) is a common welfare problem in horses worldwide. It is characterized by insulin dysregulation (ID), predisposition to laminitis and often obesity. EMS is multifactorial by nature, with both the environment and genetics contributing to the phenotype. Environmental factors, such as feeding and exercise, can be controlled, thus forming the basis for treatment and prevention. Genetic factors, by contrast, are less well-known and not easily controllable. The aim of this study was to identify potential genetic loci influencing ID/EMS in Finnhorses. A single-breed (Finnhorse) case-control genome-wide association study (GWAS) of ID was conducted with controls that included age-appropriate non-ID horses. ID status was determined with an oral sugar test (OST) for fasted horses. Seventy-one Finnhorses participated (n = 34 ID, n = 37 control). DNA samples (hair roots) were genotyped for 65 157 single-nucleotide polymorphisms (SNPs) with the Illumina Equine SNP70 BeadChip, and these data were analysed for association and FST outliers with genomic tools. P-values that exceeded the suggestive threshold (P = 1.00 ×10-5) were found in SNP BIEC2_383954 (P = 3.45 ×10-6) in chromosome 17 and SNP BIEC2_312374 (P = 1.89 ×10-5) in chromosome 15. Hierarchical and Bayesian FST outlier tests also detected these SNPs. Potential candidate genes associated with the ID close to SNP BIEC2_383954, with functions in carbohydrate metabolism, were Arginine and Glutamate Rich 1 (ARGLU1) and Ephrin-B2 (EFNB2).

Genome-wide methylation patterns from canine nanopore assemblies

Recent advances in long-read sequencing have enabled the creation of reference-quality genome assemblies for multiple individuals within a species. In particular, 8 long-read genome assemblies have recently been published for the canine model (dogs and wolves). These assemblies were created using a range of sequencing and computational approaches, with only limited comparisons described among subsets of the assemblies. Here we present 3 high-quality de novo reference assemblies based upon Oxford Nanopore long-read sequencing: 2 Bernese Mountain Dogs (BD & OD) and a Cairn terrier (CA611). These breeds are of particular interest due to the enrichment of unresolved genetic disorders. Leveraging advancement in software technologies, we utilized published data of Labrador Retriever (Yella) to generate a new assembly, resulting in a ∼280-fold increase in continuity (N50 size of 91 kbp vs 25.75 Mbp). In conjunction with these 4 new assemblies, we uniformly assessed 8 existing assemblies for generalized quality metrics, sequence divergence, and a detailed BUSCO assessment. We identified a set of ∼400 conserved genes during the BUSCO analysis missing in all assemblies. Genome-wide methylation profiles were generated from the nanopore sequencing, resulting in broad concordance with existing whole-genome and reduced-representation bisulfite sequencing, while highlighting superior overage of mobile elements. These analyses demonstrate the ability of Nanopore sequencing to resolve the sequence and epigenetic profile of canine genomes.

Identification of genetic variants associated with anterior cruciate ligament rupture and AKC standard coat color in the Labrador Retriever

Canine anterior cruciate ligament (ACL) rupture is a common complex disease. Prevalence of ACL rupture is breed dependent. In an epidemiological study, yellow coat color was associated with increased risk of ACL rupture in the Labrador Retriever. ACL rupture risk variants may be linked to coat color through genetic selection or through linkage with coat color genes. To investigate these associations, Labrador Retrievers were phenotyped as ACL rupture case or controls and for coat color and were single nucleotide polymorphism (SNP) genotyped. After filtering, ~ 697 K SNPs were analyzed using GEMMA and mvBIMBAM for multivariate association. Functional annotation clustering analysis with DAVID was performed on candidate genes. A large 8 Mb region on chromosome 5 that included ACSF3, as well as 32 additional SNPs, met genome-wide significance at P < 6.07E-7 or Log10(BF) = 3.0 for GEMMA and mvBIMBAM, respectively. On chromosome 23, SNPs were located within or near PCCB and MSL2. On chromosome 30, a SNP was located within IGDCC3. SNPs associated with coat color were also located within ADAM9, FAM109B, SULT1C4, RTDR1, BCR, and RGS7. DZIP1L was associated with ACL rupture. Several significant SNPs on chromosomes 2, 3, 7, 24, and 26 were located within uncharacterized regions or long non-coding RNA sequences. This study validates associations with the previous ACL rupture candidate genes ACSF3 and DZIP1L and identifies novel candidate genes. These variants could act as targets for treatment or as factors in disease prediction modeling. The study highlighted the importance of regulatory SNPs in the disease, as several significant SNPs were located within non-coding regions.

Special Issue: "Canine Genetics 2"

Wolves were the first animal species to become domesticated by humans, approximately 30,000-50,000 years ago. Human-directed dog breeding over thousands of generations has generated more than 350 recognized breeds displaying surprisingly different phenotypes with respect to morphology, behavior and disease predispositions. The domestication of wolves and the subsequent breeding of dogs can be viewed as one of humankind's oldest and largest genetic experiments and provides us with unique opportunities for research. Dogs have not only become human's best friend but were also described as geneticists' best friend in a past issue of Science. In recognition of the importance of canine genetics, this Special Issue, entitled "Canine Genetics 2", was compiled. It represents a sequel to the former Special Issue "Canine Genetics", which was published in 2019. During the last 15 years, the canine community has heavily relied on a reference genome derived from the female Boxer Tasha. "Canine Genetics 2" includes an article describing a greatly improved version of this important community resource. This Special Issue further contains several reports related to monogenic or complex inherited diseases in dogs. Finally, important aspects of wild canid research, genetic diversity in different populations and canine morphology were investigated.

Genome-wide association study of mammary gland tumors in Maltese dogs

Background

A genome-wide association study (GWAS) is a valuable tool for investigating genetic and phenotypic variation in many diseases.

Objective

The objective of this study was to identify variations in the genomes of Maltese dogs that are associated with the mammary gland tumor (MGT) phenotype and to assess the association between each biological condition and MGT phenotype in Maltese dogs.

Methods

DNA was extracted from 22 tumor samples and 11 whole blood samples from dogs with MGTs. Genome-wide single-nucleotide polymorphism (SNP) genotyping was performed, and the top 20 SNPs associated with various conditions and genetic variations were mapped to their corresponding gene locations.

Results

The genotyping process successfully identified 173,662 loci, with an overall genotype completion rate of 99.92%. Through the quality control analysis, 46,912 of these SNPs were excluded. Allelic tests were conducted to generate Manhattan plots, which showed several significant SNPs associated with MGT phenotype in intergenic region. The most prominent SNP, located within a region associated with transcription and linked to the malignancy grade of MGT, was identified on chromosome 5 (p = 0.00001) though there may be lack of statistical significance. Other SNPs were also found in several genes associated with oncogenesis, including TNFSF18, WDR3, ASIC5, STAR, and IL1RAP.

Conclusion

To our knowledge, this is the first GWAS to analyze the genetic predisposition to MGT in Maltese dogs. Despite the limited number of cases, these analyzed data could provide the basis for further research on the genetic predisposition to MGTs in Maltese dogs.

Systematically identifying genetic signatures including novel SNP-clusters, nonsense variants, frame-shift INDELs, and long STR expansions that potentially link to unknown phenotypes existing in dog breeds

BACKGROUND

In light of previous studies that profiled breed-specific traits or used genome-wide association studies to refine loci associated with characteristic morphological features in dogs, the field has gained tremendous genetic insights for known dog traits observed among breeds. Here we aim to address the question from a reserve perspective: whether there are breed-specific genotypes that may underlie currently unknown phenotypes. This study provides a complete set of breed-specific genetic signatures (BSGS). Several novel BSGS with significant protein-altering effects were highlighted and validated.

RESULTS

Using the next generation whole-genome sequencing technology coupled with unsupervised machine learning for pattern recognitions, we constructed and analyzed a high-resolution sequence map for 76 breeds of 412 dogs. Genomic structures including novel single nucleotide polymorphisms (SNPs), SNP clusters, insertions, deletions (INDELs) and short tandem repeats (STRs) were uncovered mutually exclusively among breeds. We also partially validated some novel nonsense variants by Sanger sequencing with additional dogs. Four novel nonsense BSGS were found in the Bernese Mountain Dog, Samoyed, Bull Terrier, and Basset Hound, respectively. Four INDELs resulting in either frame-shift or codon disruptions were found in the Norwich Terrier, Airedale Terrier, Chow Chow and Bernese Mountain Dog, respectively. A total of 15 genomic regions containing three types of BSGS (SNP-clusters, INDELs and STRs) were identified in the Akita, Alaskan Malamute, Chow Chow, Field Spaniel, Keeshond, Shetland Sheepdog and Sussex Spaniel, in which Keeshond and Sussex Spaniel each carried one amino-acid changing BSGS in such regions.

CONCLUSION

Given the strong relationship between human and dog breed-specific traits, this study might be of considerable interest to researchers and all. Novel genetic signatures that can differentiate dog breeds were uncovered. Several functional genetic signatures might indicate potentially breed-specific unknown phenotypic traits or disease predispositions. These results open the door for further investigations. Importantly, the computational tools we developed can be applied to any dog breeds as well as other species. This study will stimulate new thinking, as the results of breed-specific genetic signatures may offer an overarching relevance of the animal models to human health and disease.

Social determinants of health and disease in companion dogs: a cohort study from the Dog Aging Project

Exposure to social environmental adversity is associated with health and survival across many social species, including humans. However, little is known about how these health and mortality effects vary across the lifespan and may be differentially impacted by various components of the environment. Here, we leveraged a relatively new and powerful model for human aging, the companion dog, to investigate which components of the social environment are associated with dog health and how these associations vary across the lifespan. We drew on comprehensive survey data collected on 21,410 dogs from the Dog Aging Project and identified five factors that together explained 33.7% of the variation in a dog's social environment. Factors capturing financial and household adversity were associated with poorer health and lower physical mobility in companion dogs, while factors that captured social support, such as living with other dogs, were associated with better health when controlling for dog age and weight. Notably, the effects of each environmental component were not equal: the effect of social support was 5× stronger than financial factors. The strength of these associations depended on the age of the dog, including a stronger relationship between the owner's age and the dog's health in younger as compared to older dogs. Taken together, these findings suggest the importance of income, stability and owner's age on owner-reported health outcomes in companion dogs and point to potential behavioral and/or environmental modifiers that can be used to promote healthy aging across species.

Translational veterinary epilepsy: A win-win situation for human and veterinary neurology

Epilepsy is a challenging multifactorial disorder with a complex genetic background. Our current understanding of the pathophysiology and treatment of epilepsy has substantially increased due to animal model studies, including canine studies, but additional basic and clinical research is required. Drug-resistant epilepsy is an important problem in both dogs and humans, since seizure freedom is not achieved with the available antiseizure medications. The evaluation and exploration of pharmacological and particularly non-pharmacological therapeutic options need to remain a priority in epilepsy research. Combined efforts and sharing knowledge and expertise between human medical and veterinary neurologists are important for improving the treatment outcomes or even curing epilepsy in dogs. Such interactions could offer an exciting approach to translate the knowledge gained from people and rodents to dogs and vice versa. In this article, a panel of experts discusses the similarities and knowledge gaps in human and animal epileptology, with the aim of establishing a common framework and the basis for future translational epilepsy research.

Identification of Genetic Risk Factors for Monogenic and Complex Canine Diseases

Advances in DNA sequencing and other technologies have greatly facilitated the identification of genetic risk factors for inherited diseases in dogs. We review recent technological developments based on selected examples from canine disease genetics. The identification of disease-causing variants in dogs with monogenic diseases may become a widely employed diagnostic approach in clinical veterinary medicine in the not-too-distant future. Diseases with complex modes of inheritance continue to pose challenges to researchers but have also become much more tangible than in the past. In addition to strategies for identifying genetic risk factors, we provide some thoughts on the interpretation of sequence variants that are largely inspired by developments in human clinical genetics.

Across-breed genetic investigation of canine hip dysplasia, elbow dysplasia, and anterior cruciate ligament rupture using whole-genome sequencing

Here, we report the use of genome-wide association study (GWAS) for the analysis of canine whole-genome sequencing (WGS) repository data using breed phenotypes. Single-nucleotide polymorphisms (SNPs) were called from WGS data from 648 dogs that included 119 breeds from the Dog10K Genomes Project. Next, we assigned breed phenotypes for hip dysplasia (Orthopedic Foundation for Animals (OFA) HD, n = 230 dogs from 27 breeds; hospital HD, n = 279 dogs from 38 breeds), elbow dysplasia (ED, n = 230 dogs from 27 breeds), and anterior cruciate ligament rupture (ACL rupture, n = 279 dogs from 38 breeds), the three most important canine spontaneous complex orthopedic diseases. Substantial morbidity is common with these diseases. Previous within- and between-breed GWAS for HD, ED, and ACL rupture using array SNPs have identified disease-associated loci. Individual disease phenotypes are lacking in repository data. There is a critical knowledge gap regarding the optimal approach to undertake categorical GWAS without individual phenotypes. We considered four GWAS approaches: a classical linear mixed model, a haplotype-based model, a binary case-control model, and a weighted least squares model using SNP average allelic frequency. We found that categorical GWAS was able to validate HD candidate loci. Additionally, we discovered novel candidate loci and genes for all three diseases, including FBX025, IL1A, IL1B, COL27A1, SPRED2 (HD), UGDH, FAF1 (ED), TGIF2 (ED & ACL rupture), and IL22, IL26, CSMD1, LDHA, and TNS1 (ACL rupture). Therefore, categorical GWAS of ancestral dog populations may contribute to the understanding of any disease for which breed epidemiological risk data are available, including diseases for which GWAS has not been performed and candidate loci remain elusive.

How size and genetic diversity shape lifespan across breeds of purebred dogs

While the lifespan advantage of small body size and mixed breed status has been documented repeatedly, evidence for an effect of genetic diversity across dog breeds is equivocal. We hypothesized that this might be due to a strong right-censoring bias in available breed-specific lifespan estimates where early-dying dogs from birth cohorts that have not died off completely at the time of data collection are sampled disproportionately, especially in breeds with rapidly growing populations. We took advantage of data on owner reported lifespan and cause of death from a large public database to quantify the effect of size and genetic diversity (heterozygosity) on mortality patterns across 118 breeds based on more than 40,000 dogs. After documenting and removing the right-censoring bias from the breed-specific lifespan estimates by including only completed birth cohorts in our analyses, we show that small size and genetic diversity are both linked to a significant increase in mean lifespan across breeds. To better understand the proximate mechanisms underlying these patterns, we then investigated two major mortality causes in dogs - the cumulative pathophysiologies of old age and cancer. Old age lifespan, as well as the percentage of old age mortality, decreased with increasing body size and increased with increasing genetic diversity. The lifespan of dogs dying of cancer followed the same patterns, but while large size significantly increased proportional cancer mortality, we could not detect a significant signal for lowered cancer mortality with increasing diversity. Our findings suggest that outcross programs will be beneficial for breed health and longevity. They also emphasize the need for high-quality mortality data for veterinary epidemiology as well as for developing the dog as a translational model for human geroscience.

Four novel genes associated with longevity found in Cane corso purebred dogs

Background

Longevity-related genes have been found in several animal species as well as in humans. The goal of this study was to perform genetic analysis of long-lived Cane corso dogs with the aim to find genes that are associated with longevity.

Results

SNPs with particular nucleotides were significantly overrepresented in long-lived dogs in four genes, TDRP, MC2R, FBXO25 and FBXL21. In FBXL21, the longevity-associated SNP localises to the exon. In the FBXL21 protein, tryptophan in long-lived dogs replaced arginine present in reference dogs.

Conclusions

Four SNPs associated with longevity in dogs were identified using GWAS and validated by DNA sequencing. We conclude that genes TDRP, MC2R, FBXO25 and FBXL21 are associated with longevity in Cane corso dogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03290-9.

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.

Canine Cancer Cells Activate Platelets via the Platelet P2Y12 Receptor

In addition to their well-known functions in haemostasis, anucleated platelets have a critical role in cancer biology. Many human and non-human cancer types can directly interact with and activate platelets, promoting cancer malignancy and progression. Although naturally occurring canine neoplastic diseases mimic the biologically complex conditions of human cancers more closely than laboratory-bred mice, studies evaluating the relationship between cancer cells and platelets in dogs are scarce, and the effects of tumour cells on platelets in these animals are unknown. To evaluate whether cancer cells could activate canine platelets, we assessed the response of platelet-rich plasma to cultured canine cancer cells using light transmittance aggregometry. Similar to human and murine cancer cell research, we demonstrated that both canine osteosarcoma and mammary carcinoma cells activated canine platelets in vitro, resulting in platelet aggregation. The degree of aggregation was most pronounced at a cancer cell to platelet ratio of 1:200 for most cell lines. Mechanistic studies revealed that the platelet adenosine diphosphate (ADP) receptor P2Y12 is essential for canine platelet aggregation induced by canine cancer. ADP receptor blockage on platelets inhibited >50% of cancer cell-induced maximum platelet aggregation in all cell lines evaluated. As in other species, our results suggest that canine cancers may activate canine platelets in vivo. This mechanism is likely relevant for the biology and progression of cancer in the dog.

Genetics of inherited skin disorders in dogs

Canine genodermatoses represent a broad spectrum of diseases with diverse phenotypes. Modern genetic technology including whole genome sequencing has expedited the identification of novel genes and greatly simplified the establishment of genetic diagnoses in such heterogeneous disorders. The precise genetic diagnosis of a heritable skin disorder is essential for the appropriate counselling of owners regarding the course of the disease, prognosis and implications for breeding. Understanding the underlying pathophysiology is a prerequisite to developing specific, targeted or individualized therapeutic approaches. This review aims to create a comprehensive overview of canine genodermatoses and their respective genetic aetiology known to date. Raising awareness of genodermatoses in dogs is important and this review may help clinicians to apply modern genetics in daily clinical practice, so that a precise diagnoses can be established in suspected genodermatoses.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

The Rat Genome Database (RGD) facilitates genomic and phenotypic data integration across multiple species for biomedical research

Model organism research is essential for discovering the mechanisms of human diseases by defining biologically meaningful gene to disease relationships. The Rat Genome Database (RGD, ( https://rgd.mcw.edu )) is a cross-species knowledgebase and the premier online resource for rat genetic and physiologic data. This rich resource is enhanced by the inclusion and integration of comparative data for human and mouse, as well as other human disease models including chinchilla, dog, bonobo, pig, 13-lined ground squirrel, green monkey, and naked mole-rat. Functional information has been added to records via the assignment of annotations based on sequence similarity to human, rat, and mouse genes. RGD has also imported well-supported cross-species data from external resources. To enable use of these data, RGD has developed a robust infrastructure of standardized ontologies, data formats, and disease- and species-centric portals, complemented with a suite of innovative tools for discovery and analysis. Using examples of single-gene and polygenic human diseases, we illustrate how data from multiple species can help to identify or confirm a gene as involved in a disease and to identify model organisms that can be studied to understand the pathophysiology of a gene or pathway. The ultimate aim of this report is to demonstrate the utility of RGD not only as the core resource for the rat research community but also as a source of bioinformatic tools to support a wider audience, empowering the search for appropriate models for human afflictions.

The genomic profiling and MAMLD1 expression in human and canines with Cushing's disease

BACKGROUND

Cushing's disease (CD) is defined as hypercortisolemia caused by adrenocorticotropic hormone (ACTH)-secreting pituitary adenomas (corticotroph PA) that afflicts humans and dogs. In order to map common aberrant genomic features of CD between humans and dogs, we performed genomic sequencing and immunostaining on corticotroph PA.

METHODS

For inclusion, humans and dog were diagnosed with CD. Whole exome sequencing (WES) was conducted on 6 human corticotroph PA. Transcriptome RNA-Seq was performed on 6 human and 7 dog corticotroph PA. Immunohistochemistry (IHC) was complete on 31 human corticotroph PA. Corticotroph PA were compared with normal tissue and between species analysis were also performed.

RESULTS

Eight genes (MAMLD1, MNX1, RASEF, TBX19, BIRC5, TK1, GLDC, FAM131B) were significantly (P < 0.05) overexpressed across human and canine corticotroph PA. IHC revealed MAMLD1 to be positively (3+) expressed in the nucleus of ACTH-secreting tumor cells of human corticotroph PA (22/31, 70.9%), but absent in healthy human pituitary glands.

CONCLUSIONS

In this small exploratory cohort, we provide the first preliminary insights into profiling the genomic characterizations of human and dog corticotroph PA with respect to MAMLD1 overexpression, a finding of potential direct impact to CD microadenoma diagnosis. Our study also offers a rationale for potential use of the canine model in development of precision therapeutics.

The genetic consequences of dog breed formation-Accumulation of deleterious genetic variation and fixation of mutations associated with myxomatous mitral valve disease in cavalier King Charles spaniels

Selective breeding for desirable traits in strictly controlled populations has generated an extraordinary diversity in canine morphology and behaviour, but has also led to loss of genetic variation and random entrapment of disease alleles. As a consequence, specific diseases are now prevalent in certain breeds, but whether the recent breeding practice led to an overall increase in genetic load remains unclear. Here we generate whole genome sequencing (WGS) data from 20 dogs per breed from eight breeds and document a ~10% rise in the number of derived alleles per genome at evolutionarily conserved sites in the heavily bottlenecked cavalier King Charles spaniel breed (cKCs) relative to in most breeds studied here. Our finding represents the first clear indication of a relative increase in levels of deleterious genetic variation in a specific breed, arguing that recent breeding practices probably were associated with an accumulation of genetic load in dogs. We then use the WGS data to identify candidate risk alleles for the most common cause for veterinary care in cKCs–the heart disease myxomatous mitral valve disease (MMVD). We verify a potential link to MMVD for candidate variants near the heart specific NEBL gene in a dachshund population and show that two of the NEBL candidate variants have regulatory potential in heart-derived cell lines and are associated with reduced NEBL isoform nebulette expression in papillary muscle (but not in mitral valve, nor in left ventricular wall). Alleles linked to reduced nebulette expression may hence predispose cKCs and other breeds to MMVD via loss of papillary muscle integrity.

As a consequence of selective breeding, specific disease-causing mutations have become more frequent in certain dog breeds. Whether the breeding practice also resulted in a general increase in the overall number of disease-causing mutations per dog genome is however not clear. To address this question, we compare the amount of harmful, potentially disease-causing, mutations in dogs from eight common breeds that have experienced varying degrees of intense selective breeding. We find that individuals belonging to the breed affected by the most intense breeding—cavalier King Charles spaniel (cKCs)—carry more harmful variants than other breeds, indicating that past breeding practices may have increased the overall levels of harmful genetic variation in dogs. The most common disease in cKCs is myxomatous mitral valve disease (MMVD). To identify variants linked to this disease we next characterize mutations that are common in cKCs, but rare in other breeds, and then investigate if these mutations can predict MMVD in dachshunds. We find that variants that regulate the expression of the gene NEBL in papillary muscles may increase the risk of the disease, indicating that loss of papillary muscle integrity could contribute to the development of MMVD.

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 (P_wald = 4.83x10^-9) and 19 (P_wald = 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.

We have identified two regions of the canine genome that explain a striking 35% of risk for developing histiocytic sarcoma in FCRs. The disease is uniformly lethal, affects 20% of FCRs, and parallels a cancer of the same name in humans. Both regions harbor genes involved in cell migration and cancer-related pathways. The first includes variants in regulatory regions at the tumor suppressor PIK3R6 locus that are strongly associated with histiocytic sarcoma and likely confer risk for other hematopoietic cancers. FCRs with risk alleles at the second locus demonstrate increased expression of TNFAIP6, which correlates with poor prognosis in multiple human cancers. In identifying genomic differences between affected and unaffected dogs, we advance our understanding of both canine and human health biology and set the stage for the development of diagnostic and therapeutic strategies.

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.

Long-read assembly of a Great Dane genome highlights the contribution of GC-rich sequence and mobile elements to canine genomes

Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3' end of LINE-1_Cfs (i.e., LINE-1_Cf 3'-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

Biologically Enhanced Genome-Wide Association Study Provides Further Evidence for Candidate Loci and Discovers Novel Loci That Influence Risk of Anterior Cruciate Ligament Rupture in a Dog Model

Anterior cruciate ligament (ACL) rupture is a common condition that disproportionately affects young people, 50% of whom will develop knee osteoarthritis (OA) within 10 years of rupture. ACL rupture exhibits both hereditary and environmental risk factors, but the genetic basis of the disease remains unexplained. Spontaneous ACL rupture in the dog has a similar disease presentation and progression, making it a valuable genomic model for ACL rupture. We leveraged the dog model with Bayesian mixture model (BMM) analysis (BayesRC) to identify novel and relevant genetic variants associated with ACL rupture. We performed RNA sequencing of ACL and synovial tissue and assigned single nucleotide polymorphisms (SNPs) within differentially expressed genes to biological prior classes. SNPs with the largest effects were on chromosomes 3, 5, 7, 9, and 24. Selection signature analysis identified several regions under selection in ACL rupture cases compared to controls. These selection signatures overlapped with genome-wide associations with ACL rupture as well as morphological traits. Notable findings include differentially expressed ACSF3 with MC1R (coat color) and an association on chromosome 7 that overlaps the boundaries of SMAD2 (weight and body size). Smaller effect associations were within or near genes associated with regulation of the actin cytoskeleton and the extracellular matrix, including several collagen genes. The results of the current analysis are consistent with previous work published by our laboratory and others, and also highlight new genes in biological pathways that have not previously been associated with ACL rupture. The genetic associations identified in this study mirror those found in human beings, which lays the groundwork for development of disease-modifying therapies for both species.

Genetics of canine diabetes mellitus part 1: Phenotypes of disease

This two-part article discusses the mechanisms by which genetic variation can influence the risk of complex diseases, with a focus on canine diabetes mellitus. In Part 1, presented here, the importance of accurate methods for classifying different types of diabetes will be discussed, since this underpins the selection of cases and controls for genetic studies. Part 2 will focus on our current understanding of the genes involved in diabetes risk, and the way in which new genome sequencing technologies are poised to reveal new diabetes genes in veterinary species.

Genetics of canine diabetes mellitus part 2: Current understanding and future directions

Part 1 of this 2-part review outlined the importance of disease classification in diabetes genetic studies, as well as the ways in which genetic variants may contribute to risk of a complex disease within an individual, or within a particular group of individuals. Part 2, presented here, describes in more detail our current understanding of the genetics of canine diabetes mellitus compared to our knowledge of the human disease. Ongoing work to improve our knowledge, using new technologies, is also introduced.

Analysis of copy number variation in dogs implicates genomic structural variation in the development of anterior cruciate ligament rupture

Anterior cruciate ligament (ACL) rupture is an important condition of the human knee. Second ruptures are common and societal costs are substantial. Canine cranial cruciate ligament (CCL) rupture closely models the human disease. CCL rupture is common in the Labrador Retriever (5.79% prevalence), ~100-fold more prevalent than in humans. Labrador Retriever CCL rupture is a polygenic complex disease, based on genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers. Dissection of genetic variation in complex traits can be enhanced by studying structural variation, including copy number variants (CNVs). Dogs are an ideal model for CNV research because of reduced genetic variability within breeds and extensive phenotypic diversity across breeds. We studied the genetic etiology of CCL rupture by association analysis of CNV regions (CNVRs) using 110 case and 164 control Labrador Retrievers. CNVs were called from SNPs using three different programs (PennCNV, CNVPartition, and QuantiSNP). After quality control, CNV calls were combined to create CNVRs using ParseCNV and an association analysis was performed. We found no strong effect CNVRs but found 46 small effect (max(T) permutation P<0.05) CCL rupture associated CNVRs in 22 autosomes; 25 were deletions and 21 were duplications. Of the 46 CCL rupture associated CNVRs, we identified 39 unique regions. Thirty four were identified by a single calling algorithm, 3 were identified by two calling algorithms, and 2 were identified by all three algorithms. For 42 of the associated CNVRs, frequency in the population was <10% while 4 occurred at a frequency in the population ranging from 10–25%. Average CNVR length was 198,872bp and CNVRs covered 0.11 to 0.15% of the genome. All CNVRs were associated with case status. CNVRs did not overlap previous canine CCL rupture risk loci identified by GWAS. Associated CNVRs contained 152 annotated genes; 12 CNVRs did not have genes mapped to CanFam3.1. Using pathway analysis, a cluster of 19 homeobox domain transcript regulator genes was associated with CCL rupture (P = 6.6E-13). This gene cluster influences cranial-caudal body pattern formation during embryonic limb development. Clustered genes were found in 3 CNVRs on chromosome 14 (HoxA), 28 (NKX6-2), and 36 (HoxD). When analysis was limited to deletion CNVRs, the association was strengthened (P = 8.7E-16). This study suggests a component of the polygenic risk of CCL rupture in Labrador Retrievers is associated with small effect CNVs and may include aspects of stifle morphology regulated by homeobox domain transcript regulator genes.

Comparison of inherited neural tube defects in companion animals and livestock

Neural tube defects (NTDs) are congenital malformations resulting from the improper or incomplete closure of the neural tube during embryonic development. A number of similar malformations of the protective coverings surrounding the central nervous system are also often included under this umbrella term, which may not strictly fit this definition. A range of NTD phenotypes exist and have been reported in humans and a wide range of domestic and livestock species. In the veterinary literature, these include cases of anencephaly, encephalocele, dermoid sinus, spina bifida, and craniorachischisis. While environmental factors have a role, genetic predisposition may account for a significant part of the risk of NTDs in these animal cases. Studies of laboratory model species (fish, birds, amphibians, and rodents) have been instrumental in improving our understanding of the neurulation process. In mice, over 200 genes that may be involved in this process have been identified and variant phenotypes investigated. Like laboratory mouse models, domestic animals and livestock species display a wide range of NTD phenotypes. They remain, however, a largely underutilized population and could complement already established laboratory models. Here we review reports of NTDs in companion animals and livestock, and compare these to other animal species and human cases. We aim to highlight the potential of nonlaboratory animal models for mutation discovery as well as general insights into the mechanisms of neurulation and the development of NTDs.

The Genetic Basis of Obesity and Related Metabolic Diseases in Humans and Companion Animals

Obesity is one of the most prevalent health conditions in humans and companion animals globally. It is associated with premature mortality, metabolic dysfunction, and multiple health conditions across species. Obesity is, therefore, of importance in the fields of medicine and veterinary medicine. The regulation of adiposity is a homeostatic process vulnerable to disruption by a multitude of genetic and environmental factors. It is well established that the heritability of obesity is high in humans and laboratory animals, with ample evidence that the same is true in companion animals. In this review, we provide an overview of how genes link to obesity in humans, drawing on a wealth of information from laboratory animal models, and summarise the mechanisms by which obesity causes related disease. Throughout, we focus on how large-scale human studies and niche investigations of rare mutations in severely affected patients have improved our understanding of obesity biology and can inform our ability to interpret results of animal studies. For dogs, cats, and horses, we compare the similarities in obesity pathophysiology to humans and review the genetic studies that have been previously reported in those species. Finally, we discuss how veterinary genetics may learn from humans about studying precise, nuanced phenotypes and implementing large-scale studies, but also how veterinary studies may be able to look past clinical findings to mechanistic ones and demonstrate translational benefits to human research.

A 50-kb deletion disrupting the RSPO2 gene is associated with tetradysmelia in Holstein Friesian cattle

Background

Tetradysmelia is a rare genetic disorder that is characterized by an extremely severe reduction of all limb parts distal of the scapula and pelvic girdle. We studied a Holstein Friesian backcross family with 24 offspring, among which six calves displayed autosomal recessive tetradysmelia. In order to identify the genetic basis of the disorder, we genotyped three affected calves, five dams and nine unaffected siblings using a Bovine Illumina 50 k BeadChip and sequenced the whole genome of the sire.

Results

Pathological examination of four tetradysmelia cases revealed a uniform and severe dysmelia of all limbs. Applying a homozygosity mapping approach, we identified a homozygous region of 10.54 Mb on chromosome 14 (Bos taurus BTA14). Only calves that were diagnosed with tetradysmelia shared a distinct homozygous haplotype for this region. We sequenced the whole genome of the cases’ sire and searched for heterozygous single nucleotide polymorphisms (SNPs) and small variants on BTA14 that were uniquely present in the sire and absent from 3102 control whole-genome sequences of the 1000 Bull Genomes Project, but none were identified in the 10.54-Mb candidate region on BTA14. Therefore, we subsequently performed a more comprehensive analysis by also considering structural variants and detected a 50-kb deletion in the targeted chromosomal region that was in the heterozygous state in the cases’ sire. Using PCR, we confirmed that this detected deletion segregated perfectly within the family with tetradysmelia. The deletion spanned three exons of the bovine R-spondin 2 (RSPO2) gene, which encode three domains of the respective protein. R-spondin 2 is a secreted ligand of leucine-rich repeats containing G protein-coupled receptors that enhance Wnt signalling and is involved in a broad range of developmental processes during embryogenesis.

Conclusions

We identified a 50-kb deletion on BTA14 that disrupts the coding sequence of the RSPO2 gene and is associated with bovine tetradysmelia. To our knowledge, this is the first reported candidate causal mutation for tetradysmelia in a large animal model. Since signalling pathways involved in limb development are conserved across species, the observed inherited defect may serve as a model to further elucidate fundamental pathways of limb development.

Animal models of naturally occurring stone disease

The prevalence of urolithiasis in humans is increasing worldwide; however, non-surgical treatment and prevention options remain limited despite decades of investigation. Most existing laboratory animal models for urolithiasis rely on highly artificial methods of stone induction and, as a result, might not be fully applicable to the study of natural stone initiation and growth. Animal models that naturally and spontaneously form uroliths are an underused resource in the study of human stone disease and offer many potential opportunities for improving insight into stone pathogenesis. These models include domestic dogs and cats, as well as a variety of other captive and wild species, such as otters, dolphins and ferrets, that form calcium oxalate, struvite, uric acid, cystine and other stone types. Improved collaboration between urologists, basic scientists and veterinarians is warranted to further our understanding of how stones form and to consider possible new preventive and therapeutic treatment options.

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.

Bayesian and Machine Learning Models for Genomic Prediction of Anterior Cruciate Ligament Rupture in the Canine Model

Anterior cruciate ligament (ACL) rupture is a common, debilitating condition that leads to early-onset osteoarthritis and reduced quality of human life. ACL rupture is a complex disease with both genetic and environmental risk factors. Characterizing the genetic basis of ACL rupture would provide the ability to identify individuals that have high genetic risk and allow the opportunity for preventative management. Spontaneous ACL rupture is also common in dogs and shows a similar clinical presentation and progression. Thus, the dog has emerged as an excellent genomic model for human ACL rupture. Genome-wide association studies (GWAS) in the dog have identified a number of candidate genetic variants, but research in genomic prediction has been limited. In this analysis, we explore several Bayesian and machine learning models for genomic prediction of ACL rupture in the Labrador Retriever dog. Our work demonstrates the feasibility of predicting ACL rupture from SNPs in the Labrador Retriever model with and without consideration of non-genetic risk factors. Genomic prediction including non-genetic risk factors approached clinical relevance using multiple linear Bayesian and non-linear models. This analysis represents the first steps toward development of a predictive algorithm for ACL rupture in the Labrador Retriever model. Future work may extend this algorithm to other high-risk breeds of dog. The ability to accurately predict individual dogs at high risk for ACL rupture would identify candidates for clinical trials that would benefit both veterinary and human medicine.

Breed Differences in Dog Cognition Associated with Brain-Expressed Genes and Neurological Functions

Given their remarkable phenotypic diversity, dogs present a unique opportunity for investigating the genetic bases of cognitive and behavioral traits. Our previous work demonstrated that genetic relatedness among breeds accounts for a substantial portion of variation in dog cognition. Here, we investigated the genetic architecture of breed differences in cognition, seeking to identify genes that contribute to variation in cognitive phenotypes. To do so, we combined cognitive data from the citizen science project Dognition.com with published breed-average genetic polymorphism data, resulting in a dataset of 1654 individuals with cognitive phenotypes representing 49 breeds. We conducted a breed-average genome-wide association study to identify specific polymorphisms associated with breed differences in inhibitory control, communication, memory, and physical reasoning. We found five single nucleotide polymorphisms (SNPs) that reached genome-wide significance after Bonferroni correction, located in EML1, OR52E2, HS3ST5, a U6 spliceosomal RNA, and a long noncoding RNA. When we combined results across multiple SNPs within the same gene, we identified 188 genes implicated in breed differences in cognition. This gene set included more genes than expected by chance that were (1) differentially expressed in brain tissue and (2) involved in nervous system functions including peripheral nervous system development, Wnt signaling, presynapse assembly, and synaptic vesicle exocytosis. These results advance our understanding of the genetic underpinnings of complex cognitive phenotypes and identify specific genetic variants for further research.

A scan for genes associated with cancer mortality and longevity in pedigree dog breeds

Selective breeding of the domestic dog (Canis lupus familiaris) rigidly retains desirable features, and could inadvertently fix disease-causing variants within a breed. We combine phenotypic data from > 72,000 dogs with a large genotypic dataset to search for genes associated with cancer mortality and longevity in pedigree dog breeds. We validated previous findings that breeds with higher average body weight have higher cancer mortality rates and lower life expectancy. We identified a significant positive correlation between life span and cancer mortality residuals corrected for body weight, implying that long-lived breeds die more frequently from cancer compared to short-lived breeds. We replicated a number of known genetic associations with body weight (IGF1, GHR, CD36, SMAD2 and IGF2BP2). Subsequently, we identified five genetic variants in known cancer-related genes (located within SIPA1, ADCY7 and ARNT2) that could be associated with cancer mortality residuals corrected for confounding factors. One putative genetic variant was marginally significantly associated with longevity residuals that had been corrected for the effects of body weight; this genetic variant is located within PRDX1, a peroxiredoxin that belongs to an emerging class of pro-longevity associated genes. This research should be considered as an exploratory analysis to uncover associations between genes and longevity/cancer mortality.

Genome-wide association analysis of canine T zone lymphoma identifies link to hypothyroidism and a shared association with mast-cell tumors

Background

T zone lymphoma (TZL), a histologic variant of peripheral T cell lymphoma, represents about 12% of all canine lymphomas. Golden Retrievers appear predisposed, representing over 40% of TZL cases. Prior research found that asymptomatic aged Golden Retrievers frequently have populations of T zone-like cells (phenotypically identical to TZL) of undetermined significance (TZUS), potentially representing a pre-clinical state. These findings suggest a genetic risk factor for this disease and caused us to investigate potential genes of interest using a genome-wide association study of privately-owned U.S. Golden Retrievers.

Results

Dogs were categorized as TZL (n = 95), TZUS (n = 142), or control (n = 101) using flow cytometry and genotyped using the Illumina CanineHD BeadChip. Using a mixed linear model adjusting for population stratification, we found association with genome-wide significance in regions on chromosomes 8 and 14. The chromosome 14 peak included four SNPs (Odds Ratio = 1.18–1.19, p = .3 × 10^− 5–5.1 × 10^− 5) near three hyaluronidase genes (SPAM1, HYAL4, and HYALP1). Targeted resequencing of this region using a custom sequence capture array identified missense mutations in all three genes; the variant in SPAM1 was predicted to be damaging. These mutations were also associated with risk for mast cell tumors among Golden Retrievers in an unrelated study. The chromosome 8 peak contained 7 SNPs (Odds Ratio = 1.24–1.42, p = 2.7 × 10^− 7–7.5 × 10^− 5) near genes involved in thyroid hormone regulation (DIO2 and TSHR). A prior study from our laboratory found hypothyroidism is inversely associated with TZL risk. No coding mutations were found with targeted resequencing but identified variants may play a regulatory role for all or some of the genes.

Conclusions

The pathogenesis of canine TZL may be related to hyaluronan breakdown and subsequent production of pro-inflammatory and pro-oncogenic byproducts. The association on chromosome 8 may indicate thyroid hormone is involved in TZL development, consistent with findings from a previous study evaluating epidemiologic risk factors for TZL. Future work is needed to elucidate these mechanisms.

Whole-genome genotyping and resequencing reveal the association of a deletion in the complex interferon alpha gene cluster with hypothyroidism in dogs

Background

Hypothyroidism is a common complex endocrinopathy that typically has an autoimmune etiology, and it affects both humans and dogs. Genetic and environmental factors are both known to play important roles in the disease development. In this study, we sought to identify the genetic risk factors potentially involved in the susceptibility to the disease in the high-risk Giant Schnauzer dog breed.

Results

By employing genome-wide association followed by fine-mapping (top variant p-value = 5.7 × 10^− 6), integrated with whole-genome resequencing and copy number variation analysis, we detected a ~ 8.9 kbp deletion strongly associated (p-value = 0.0001) with protection against development of hypothyroidism. The deletion is located between two predicted Interferon alpha (IFNA) genes and it may eliminate functional elements potentially involved in the transcriptional regulation of these genes. Remarkably, type I IFNs have been extensively associated to human autoimmune hypothyroidism and general autoimmunity. Nonetheless, the extreme genomic complexity of the associated region on CFA11 warrants further long-read sequencing and annotation efforts in order to ascribe functions to the identified deletion and to characterize the canine IFNA gene cluster in more detail.

Conclusions

Our results expand the current knowledge on genetic determinants of canine hypothyroidism by revealing a significant link with the human counterpart disease, potentially translating into better diagnostic tools across species, and may contribute to improved canine breeding strategies.

Differential stromal reprogramming in benign and malignant naturally occurring canine mammary tumours identifies disease-modulating stromal components

While cancer-associated stroma (CAS) in malignant tumours is well described, stromal changes in benign forms of naturally occurring tumours remain poorly characterized. Spontaneous canine mammary carcinomas (mCA) are viewed as excellent models of human mCA. We have recently reported highly conserved stromal reprogramming between canine and human mCA based on transcriptome analysis of laser-capture-microdissected FFPE specimen. To identify stromal changes between benign and malignant mammary tumours, we have analysed matched normal and adenoma-associated stroma (AAS) from 13 canine mammary adenomas and compared them to previous data from 15 canine mCA. Our analyses reveal distinct stromal reprogramming even in small benign tumours. While similarities between AAS and CAS exist, the stromal signature clearly distinguished adenomas from mCA. The distinction between AAS and CAS is further substantiated by differential enrichment in several hallmark signalling pathways as well as differential abundance in cellular composition. Finally, we identify COL11A1, VIT, CD74, HLA-DRA, STRA6, IGFBP4, PIGR, and TNIP1 as strongly discriminatory stromal genes between adenoma and mCA, and demonstrate their prognostic value for human breast cancer. Given the relevance of canine CAS as a model for the human disease, our approach identifies disease-modulating stromal components with implications for both human and canine breast cancer.

Humanity's Best Friend: A Dog-Centric Approach to Addressing Global Challenges

No other animal has a closer mutualistic relationship with humans than the dog (Canis familiaris). Domesticated from the Eurasian grey wolf (Canis lupus), dogs have evolved alongside humans over millennia in a relationship that has transformed dogs and the environments in which humans and dogs have co-inhabited. The story of the dog is the story of recent humanity, in all its biological and cultural complexity. By exploring human-dog-environment interactions throughout time and space, it is possible not only to understand vital elements of global history, but also to critically assess our present-day relationship with the natural world, and to begin to mitigate future global challenges. In this paper, co-authored by researchers from across the natural and social sciences, arts and humanities, we argue that a dog-centric approach provides a new model for future academic enquiry and engagement with both the public and the global environmental agenda.

Whole Genome Analysis of a Single Scottish Deerhound Dog Family Provides Independent Corroboration That a SGK3 Coding Variant Leads to Hairlessness

The breeds of domestic dog, Canis lupus familiaris, display a range of coat types with variation in color, texture, length, curl, and growth pattern. One trait of interest is that of partial or full hairlessness, which is found in a small number of breeds. While the standard for some breeds, such as the Xoloitzcuintli, requires sparse hair on their extremities, others are entirely bald, including the American Hairless Terrier. We identified a small, rare family of Scottish Deerhounds in which coated parents produced a mixed litter of coated and hairless offspring. To identify the underlying variant, we performed whole genome sequencing of the dam and five offspring, comparing single nucleotide polymorphisms and small insertions/deletions against an established catalog of 91 million canine variants. Of 325 homozygous alternative alleles found in both hairless dogs, 56 displayed the expected pattern of segregation and only a single, high impact variant within a coding region was observed: a single base pair insertion in exon two of SGK3 leading to a potential frameshift, thus verifying recently published findings. In addition, we observed that gene expression levels between coated and hairless dogs are similar, suggesting a mechanism other than non-sense mediated decay is responsible for the phenotype.

Know Thy Model: Charting Molecular Homology in Stromal Reprogramming Between Canine and Human Mammary Tumors

Spontaneous canine simple mammary tumors (CMTs) are often viewed as models of human breast cancer. Cancer-associated stroma (CAS) is central for initiation and progression of human cancer, and is likely to play a key role in canine tumors as well. Until recently, however, canine CAS in general, and in CMT in particular, lacked detailed characterization and it remained unclear how canine and human CAS compare. This void in knowledge regarding canine CAS and the resulting lack of unbiased cross-species analysis of molecular homologies and differences undermined the validity of the canine model for human disease. To assess stromal reprogramming in canine breast tumors, we have recently established a protocol to specifically isolate and analyze CAS and matched normal stroma from archival, formalin-fixed paraffin embedded (FFPE) clinical tumor samples using laser-capture microdissection followed by next-generation RNA-sequencing. Using this approach, we have analyzed stromal reprogramming in both malignant canine mammary carcinomas (mCAs) as well as benign canine mammary adenomas in a series of studies. Our results demonstrate strong stromal reprogramming in CMTs and identify high-grade molecular homology between human and canine CAS. Here, I aim to give a short background on the value of comparative oncology in general, and spontaneous CMT in particular. This will be followed by a concise review of the current knowledge of stromal reprogramming in both malignant canine mCA and benign adenoma. Finally, I will conclude with insights on highly conserved aspects of stromal reprogramming between CMT and human breast cancer that accentuate the relevance of CAS in CMT as a model for the human disease.

Imputation of canine genotype array data using 365 whole-genome sequences improves power of genome-wide association studies

Genomic resources for the domestic dog have improved with the widespread adoption of a 173k SNP array platform and updated reference genome. SNP arrays of this density are sufficient for detecting genetic associations within breeds but are underpowered for finding associations across multiple breeds or in mixed-breed dogs, where linkage disequilibrium rapidly decays between markers, even though such studies would hold particular promise for mapping complex diseases and traits. Here we introduce an imputation reference panel, consisting of 365 diverse, whole-genome sequenced dogs and wolves, which increases the number of markers that can be queried in genome-wide association studies approximately 130-fold. Using previously genotyped dogs, we show the utility of this reference panel in identifying potentially novel associations, including a locus on CFA20 significantly associated with cranial cruciate ligament disease, and fine-mapping for canine body size and blood phenotypes, even when causal loci are not in strong linkage disequilibrium with any single array marker. This reference panel resource will improve future genome-wide association studies for canine complex diseases and other phenotypes.

Complex traits are controlled by more than one gene and as such are difficult to map. For complex trait mapping in the domestic dog, researchers use the current array of 173,000 variants, with only minimal success. Here, we use a method called imputation to increase the number of variants–from 173,000 to 24 million–that can be queried in canine association studies. We use sequence data from the whole genomes of 365 dogs and wolves to accurately predict variants, in a separate cohort of dogs, that are not present on the array. Using dog body size, blood phenotypes, and a common orthopedic disease that involves rupture of the cranial cruciate ligament, we show that the increase in variants results in an increase in mapping power, through the identification of new associations and the narrowing of regions of interest. This imputation panel is particularly important because of its usefulness in improving complex trait mapping in the dog, which has significant implications for discovery of variants in humans with similar diseases.

Next-generation RNA sequencing of FFPE subsections reveals highly conserved stromal reprogramming between canine and human mammary carcinoma

Spontaneous canine simple mammary carcinomas (mCA) are often viewed as models of human mCA. Cancer-associated stroma (CAS) is central for initiation and progression of human cancer, and is likely to play a key role in canine tumours as well. However, canine CAS lacks characterisation and it remains unclear how canine and human CAS compare. Formalin-fixed paraffin embedded (FFPE) tissue constitutes a valuable resource of patient material, but chemical crosslinking has largely precluded its analysis by next-generation RNA sequencing (RNAseq). We have recently established a protocol to isolate CAS and normal stroma from archival FFPE tumours using laser-capture microdissection followed by RNAseq. Using this approach, we have analysed stroma from 15 canine mCA. Our data reveal strong reprogramming of canine CAS. We demonstrate a high-grade molecular homology between canine and human CAS, and show that enrichment of upregulated canine CAS genes strongly correlates with the enrichment of an independently derived human stromal signature in the TCGA breast tumour dataset. Relationships between different gene signatures observed in human breast cancer are largely maintained in the canine model, suggesting a close interspecies similarity in the network of cancer signalling circuitries. Finally, we establish the prognostic potential of the canine CAS signature in human samples, emphasising the relevance of studying canine CAS as a model of the human disease. In conclusion, we provide a proof-of-principle to analyse specific subsections of FFPE tissue by RNAseq, and compare stromal gene expression between human and canine mCA to reveal molecular drivers in CAS supporting tumour growth and malignancy.

Summary: This study offers proof-of-principle for a novel protocol to analyse gene expression in subsections of FFPE patient tissue, supporting the use of spontaneous canine mammary tumours as models for the human disease.

Dog10K: an international sequencing effort to advance studies of canine domestication, phenotypes and health

Dogs are the most phenotypically diverse mammalian species, and they possess more known heritable disorders than any other non-human mammal. Efforts to catalog and characterize genetic variation across well-chosen populations of canines are necessary to advance our understanding of their evolutionary history and genetic architecture. To date, no organized effort has been undertaken to sequence the world's canid populations. The Dog10K Consortium (http://www.dog10kgenomes.org) is an international collaboration of researchers from across the globe who will generate 20× whole genomes from 10 000 canids in 5 years. This effort will capture the genetic diversity that underlies the phenotypic and geographical variability of modern canids worldwide. Breeds, village dogs, niche populations and extended pedigrees are currently being sequenced, and de novo assemblies of multiple canids are being constructed. This unprecedented dataset will address the genetic underpinnings of domestication, breed formation, aging, behavior and morphological variation. More generally, this effort will advance our understanding of human and canine health.

BarkBase: Epigenomic Annotation of Canine Genomes

Dogs are an unparalleled natural model for investigating the genetics of health and disease, particularly for complex diseases like cancer. Comprehensive genomic annotation of regulatory elements active in healthy canine tissues is crucial both for identifying candidate causal variants and for designing functional studies needed to translate genetic associations into disease insight. Currently, canine geneticists rely primarily on annotations of the human or mouse genome that have been remapped to dog, an approach that misses dog-specific features. Here, we describe BarkBase, a canine epigenomic resource available at barkbase.org. BarkBase hosts data for 27 adult tissue types, with biological replicates, and for one sample of up to five tissues sampled at each of four carefully staged embryonic time points. RNA sequencing is complemented with whole genome sequencing and with assay for transposase-accessible chromatin using sequencing (ATAC-seq), which identifies open chromatin regions. By including replicates, we can more confidently discern tissue-specific transcripts and assess differential gene expression between tissues and timepoints. By offering data in easy-to-use file formats, through a visual browser modeled on similar genomic resources for human, BarkBase introduces a powerful new resource to support comparative studies in dogs and humans.

Genome-wide association analysis uncovers variants for reproductive variation across dog breeds and links to domestication

BACKGROUND AND OBJECTIVES

The diversity of eutherian reproductive strategies has led to variation in many traits, such as number of offspring, age of reproductive maturity and gestation length. While reproductive trait variation has been extensively investigated and is well established in mammals, the genetic loci contributing to this variation remain largely unknown. The domestic dog, Canis lupus familiaris is a powerful model for studies of the genetics of inherited disease due to its unique history of domestication. To gain insight into the genetic basis of reproductive traits across domestic dog breeds, we collected phenotypic data for four traits, cesarean section rate, litter size, stillbirth rate and gestation length, from primary literature and breeders' handbooks.

METHODOLOGY

By matching our phenotypic data to genomic data from the Cornell Veterinary Biobank, we performed genome-wide association analyses for these four reproductive traits, using body mass and kinship among breeds as covariates.

RESULTS

We identified 12 genome-wide significant associations between these traits and genetic loci, including variants near CACNA2D3 with gestation length, MSRB3 and MSANTD1 with litter size, SMOC2 with cesarean section rate and UFM1 with stillbirth rate. A few of these loci, such as CACNA2D3 and MSRB3, have been previously implicated in human reproductive pathologies, whereas others have been associated with domestication-related traits, including brachycephaly (SMOC2) and coat curl (KRT71).

CONCLUSIONS AND IMPLICATIONS

We hypothesize that the artificial selection that gave rise to dog breeds also influenced the observed variation in their reproductive traits. Overall, our work establishes the domestic dog as a system for studying the genetics of reproductive biology and disease.

LAY SUMMARY

The genetic contributors to variation in mammalian reproductive traits remain largely unknown. We took advantage of the domestic dog, a powerful model system, to test for associations between genome-wide variants and four reproductive traits (cesarean section rate, litter size, stillbirth rate and gestation length) that vary extensively across breeds. We identified associations at a dozen loci, including ones previously associated with domestication-related traits, suggesting that selection on dog breeds also influenced their reproductive traits.

Cdrom Archive: A Gateway to Study Camel Phenotypes

Camels are livestock that exhibit unique morphological, biochemical, and behavioral traits, which arose by natural and artificial selection. Investigating the molecular basis of camel traits has been limited by: (1) the absence of a comprehensive record of morphological trait variation (e.g., diseases) and the associated mode of inheritance, (2) the lack of extended pedigrees of specific trait(s), and (3) the long reproductive cycle of the camel, which makes the cost of establishing and maintaining a breeding colony (i.e., monitoring crosses) prohibitively high. Overcoming these challenges requires (1) detailed documentation of phenotypes/genetic diseases and their likely mode of inheritance (and collection of related DNA samples), (2) conducting association studies to identify phenotypes/genetic diseases causing genetic variants (instead of classical linkage analysis, which requires extended pedigrees), and (3) validating likely causative variants by screening a large number of camel samples from different populations. We attempt to address these issues by establishing a systematic way of collecting camel DNA samples, and associated phenotypic information, which we call the "Cdrom Archive." Here, we outline the process of building this archive to introduce it to other camel researchers (as an example). Additionally, we discuss the use of this archive to study the phenotypic traits of Arabian Peninsula camel breeds (the "Mezayen" camels). Using the Cdrom Archive, we report variable phenotypic traits related to the coat (color, length, and texture), ear and tail lengths, along with other morphological measurements.

Genome-wide association study of degenerative mitral valve disease in Maltese dogs

Genome-wide association study (GWAS) is a powerful tool for identifying the genetic causes of various diseases. This study was conducted to identify genomic variation in Maltese dog genomes associated with degenerative mitral valve disease (DMVD) development and to evaluate the association of each biological condition with DMVD in Maltese dogs. DNA was extracted from blood samples obtained from 48 Maltese dogs (32 with DMVD and 16 controls). Genome-wide single nucleotide polymorphism (SNP) genotyping was performed. The top 30 SNPs from each association of various conditions and genetic variations were mapped to their gene locations. A total of 173,662 loci were successfully genotyped, with an overall genotype completion rate of 99.41%. Quality control analysis excluded 46,610 of these SNPs. Manhattan plots were produced using allelic tests with various candidate clinical conditions. A significant peak of association was observed between mitral valve prolapse (MVP) and SNPs on chromosome 17. The present study revealed significant SNPs in several genes associated with cardiac function, including PDZ2, Armadillo repeat protein detected in velo-cardio-facial syndrome, catenin (cadherin-associated protein) alpha 3, low-density lipoprotein receptor class A domain containing protein 4, and sterile alpha motif domain containing protein 3. To our knowledge, this is the first study of a genetic predisposition to DMVD in Maltese dogs. Although only a limited number of cases were analyzed, these data could be the basis for further research on the genetic predisposition to MVP and DMVD in Maltese dogs.