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

Imputation of ancient canid genomes reveals inbreeding history over the past 10,000 years

The multi-millenia long history between dogs and humans has placed them at the forefront of archeological and genomic research. Despite ongoing efforts including the analysis of ancient dog and wolf genomes, many questions remain regarding their geographic and temporal origins, and the microevolutionary processes that led to the diversity of breeds today. Although ancient genomes provide valuable information, their use is hindered by low depth of coverage and post-mortem damage, which inhibits confident genotype calling. In the present study, we assess how genotype imputation of ancient dog and wolf genomes, utilising a large reference panel, can improve the resolution provided by ancient datasets. Imputation accuracy was evaluated by down-sampling high coverage dog and wolf genomes to 0.05-2x coverage and comparing concordance between imputed and high coverage genotypes. We measured the impact of imputation on principal component analyses and runs of homozygosity. Our findings show high (R2>0.9) imputation accuracy for dogs with coverage as low as 0.5x and for wolves as low as 1.0x. We then imputed a dataset of 90 ancient dog and wolf genomes, to assess changes in inbreeding during the last 10,000 years of dog evolution. Ancient dog and wolf populations generally exhibited lower inbreeding levels than present-day individuals. Interestingly, regions with low ROH density maintained across ancient and present-day samples were significantly associated with genes related to olfaction and immune response. Our study indicates that imputing ancient canine genomes is a viable strategy that allows for the use of analytical methods previously limited to high-quality genetic data.

A cautionary tale of low-pass sequencing and imputation with respect to haplotype accuracy

BACKGROUND

Low-pass whole-genome sequencing and imputation offer significant cost savings, enabling substantial increases in sample size and statistical power. This approach is particularly promising in livestock breeding, providing an affordable means of screening individuals for deleterious alleles or calculating genomic breeding values. Consequently, it may also be of value in companion animal genomics to support pedigree breeding. We sought to evaluate in dogs the impact of low coverage sequencing and reference-guided imputation on genotype concordance and association analyses.

RESULTS

DNA isolated from saliva of 30 Labrador retrievers was sequenced at low (0.9X and 3.8X) and high (43.5X) coverage, and down-sampled from 43.5X to 9.6X and 17.4X. Genotype imputation was performed using a diverse reference panel (1021 dogs), and two subsets of the former panel (256 dogs each) where one had an excess of Labrador retrievers relative to other breeds. We observed little difference in imputed genotype concordance between reference panels. Association analyses for a locus acting as a disease proxy were performed using single-marker (GEMMA) and haplotype-based (XP-EHH) tests. GEMMA results were highly correlated (r ≥ 0.97) between 43.5X and ≥ 3.8X depths of coverage, while for 0.9X the correlation was lower (r ≤ 0.8). XP-EHH results were less well correlated, with r ranging from 0.58 (0.9X) to 0.88 (17.4X). Across a random sample of 10,000 genomic regions averaging 17 kb in size, we observed a median of three haplotypes per dog across the sequencing depths, with 5% of the regions returning more than eight haplotypes. Inspection of one such region revealed genotype and phasing inconsistencies across sequencing depths.

CONCLUSIONS

We demonstrate that saliva-derived canine DNA is suitable for whole-genome sequencing, highlighting the feasibility of client-based sampling. Low-pass sequencing and imputation require caution as incorrect allele assignments result when the subject possesses alleles that are absent in the reference panel. Larger panels have the capacity for greater allelic diversity, which should reduce the potential for imputation error. Although low-pass sequencing can accurately impute allele dosage, we highlight issues with phasing accuracy that impact haplotype-based analyses. Consequently, if accurately phased genotypes are required for analyses, we advocate sequencing at high depth (> 20X).

GWAS using low-pass whole genome sequence reveals a novel locus in canine congenital idiopathic megaesophagus

Congenital idiopathic megaesophagus (CIM) is a gastrointestinal disorder of dogs wherein the esophagus is dilated and swallowing activity is reduced, causing regurgitation of ingesta. Affected individuals experience weight loss and malnourishment and are at risk for aspiration pneumonia, intussusception, and euthanasia. Great Danes have among the highest incidences of CIM across dog breeds, suggesting a genetic predisposition. We generated low-pass sequencing data for 83 Great Danes and used variant calls to impute missing whole genome single-nucleotide variants (SNVs) for each individual based on haplotypes phased from 624 high-coverage dog genomes, including 21 Great Danes. We validated the utility of our imputed data set for genome-wide association studies (GWASs) by mapping loci known to underlie coat phenotypes with simple and complex inheritance patterns. We conducted a GWAS for CIM with 2,010,300 SNVs, identifying a novel locus on canine chromosome 1 (P-val = 2.76 × 10-10). Associated SNVs are intergenic or intronic and are found in two clusters across a 1.7-Mb region. Inspection of coding regions in high-coverage genomes from affected Great Danes did not reveal candidate causal variants, suggesting that regulatory variants underlie CIM. Further studies are necessary to assess the role of these non-coding variants.

Best practices for genotype imputation from low-coverage sequencing data in natural populations

Monitoring genetic diversity in wild populations is a central goal of ecological and evolutionary genetics and is critical for conservation biology. However, genetic studies of nonmodel organisms generally lack access to species-specific genotyping methods (e.g. array-based genotyping) and must instead use sequencing-based approaches. Although costs are decreasing, high-coverage whole-genome sequencing (WGS), which produces the highest confidence genotypes, remains expensive. More economical reduced representation sequencing approaches fail to capture much of the genome, which can hinder downstream inference. Low-coverage WGS combined with imputation using a high-confidence reference panel is a cost-effective alternative, but the accuracy of genotyping using low-coverage WGS and imputation in nonmodel populations is still largely uncharacterized. Here, we empirically tested the accuracy of low-coverage sequencing (0.1-10×) and imputation in two natural populations, one with a large (n = 741) reference panel, rhesus macaques (Macaca mulatta), and one with a smaller (n = 68) reference panel, gelada monkeys (Theropithecus gelada). Using samples sequenced to coverage as low as 0.5×, we could impute genotypes at >95% of the sites in the reference panel with high accuracy (median r2  ≥ 0.92). We show that low-coverage imputed genotypes can reliably calculate genetic relatedness and population structure. Based on these data, we also provide best practices and recommendations for researchers who wish to deploy this approach in other populations, with all code available on GitHub (https://github.com/mwatowich/LoCSI-for-non-model-species). Our results endorse accurate and effective genotype imputation from low-coverage sequencing, enabling the cost-effective generation of population-scale genetic datasets necessary for tackling many pressing challenges of wildlife conservation.

Genome sequencing of 2000 canids by the Dog10K consortium advances the understanding of demography, genome function and architecture

BACKGROUND

The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function.

RESULTS

We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection.

CONCLUSIONS

We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available.

A new cheese population in Penicillium roqueforti and adaptation of the five populations to their ecological niche

Domestication is an excellent case study for understanding adaptation and multiple fungal lineages have been domesticated for fermenting food products. Studying domestication in fungi has thus both fundamental and applied interest. Genomic studies have revealed the existence of four populations within the blue-cheese-making fungus Penicillium roqueforti. The two cheese populations show footprints of domestication, but the adaptation of the two non-cheese populations to their ecological niches (i.e., silage/spoiled food and lumber/spoiled food) has not been investigated yet. Here, we reveal the existence of a new P. roqueforti population, specific to French Termignon cheeses, produced using small-scale traditional practices, with spontaneous blue mould colonisation. This Termignon population is genetically differentiated from the four previously identified populations, providing a novel source of genetic diversity for cheese making. The Termignon population indeed displayed substantial genetic diversity, both mating types, horizontally transferred regions previously detected in the non-Roquefort population, and intermediate phenotypes between cheese and non-cheese populations. Phenotypically, the non-Roquefort cheese population was the most differentiated, with specific traits beneficial for cheese making, in particular higher tolerance to salt, to acidic pH and to lactic acid. Our results support the view that this clonal population, used for many cheese types in multiple countries, is a domesticated lineage on which humans exerted strong selection. The lumber/spoiled food and silage/spoiled food populations were not more tolerant to crop fungicides but showed faster growth in various carbon sources (e.g., dextrose, pectin, sucrose, xylose and/or lactose), which can be beneficial in their ecological niches. Such contrasted phenotypes between P. roqueforti populations, with beneficial traits for cheese-making in the cheese populations and enhanced ability to metabolise sugars in the lumber/spoiled food population, support the inference of domestication in cheese fungi and more generally of adaptation to anthropized environments.

Common Orthopedic Traits and Screening for Breeding Programs

Orthopedic diseases are complex traits, meaning genetics and environmental factors affect risk, making identification of genetic associations difficult. In the United States, hip and elbow scores, patellar luxation scores, Legg-Calvé-Perthes disease, and shoulder osteochondrosis affectedness are available in the Orthopedic Foundation for Animals registry. Distraction indices and extended, ventrodorsal hip conformation scores are recorded by PennHIP. Application of estimated breeding values for hip and elbow dysplasia in breeder selection reduces the severity and prevalence of these traits. Genomic prediction and whole-genome sequence technologies and methods should improve knowledge of genetics underlying orthopedic diseases, leading to improved canine orthopedic genetic quality.

AmelHap: Leveraging drone whole-genome sequence data to create a honey bee HapMap

Honey bee, Apis mellifera, drones are typically haploid, developing from an unfertilized egg, inheriting only their queen's alleles and none from the many drones she mated with. Thus the ordered combination or 'phase' of alleles is known, making drones a valuable haplotype resource. We collated whole-genome sequence data for 1,407 drones, including 45 newly sequenced Scottish drones, collectively representing 19 countries, 8 subspecies and various hybrids. Following alignment to Amel_HAv3.1, variant calling and quality filtering, we retained 17.4 M high quality variants across 1,328 samples with a genotyping rate of 98.7%. We demonstrate the utility of this haplotype resource, AmelHap, for genotype imputation, returning >95% concordance when up to 61% of data is missing in haploids and up to 12% of data is missing in diploids. AmelHap will serve as a useful resource for the community for imputation from low-depth sequencing or SNP chip data, accurate phasing of diploids for association studies, and as a comprehensive reference panel for population genetic and evolutionary analyses.

Increased risk of cancer in dogs and humans: a consequence of recent extension of lifespan beyond evolutionarily-determined limitations?

Cancer is among the most common causes of death for dogs (and cats) and humans in the developed world, even though it is uncommon in wildlife and other domestic animals. We provide a rationale for this observation based on recent advances in our understanding of the evolutionary basis of cancer. Over the course of evolutionary time, species have acquired and fine-tuned adaptive cancer protective mechanisms that are intrinsically related to their energy demands, reproductive strategies, and expected lifespan. These cancer protective mechanisms are general across species and/or specific to each species and their niche, and they do not seem to be limited in diversity. The evolutionarily acquired cancer-free longevity that defines a species' life history can explain why the relative cancer risk, rate, and incidence are largely similar across most species in the animal kingdom despite differences in body size and life expectancy. The molecular, cellular, and metabolic events that promote malignant transformation and cancerous growth can overcome these adaptive, species-specific protective mechanisms in a small proportion of individuals, while independently, some individuals in the population might achieve exceptional longevity. In dogs and humans, recent dramatic alterations in healthcare and social structures have allowed increasing numbers of individuals in both species to far exceed their species-adapted longevities (by 2-4 times) without allowing the time necessary for compensatory natural selection. In other words, the cancer protective mechanisms that restrain risk at comparable levels to other species for their adapted lifespan are incapable of providing cancer protection over this recent, drastic and widespread increase in longevity.

The Effects of FGF4 Retrogenes on Canine Morphology

Two FGF4 retrogenes (FGF4L1 on chromosome 18 and FGF4L2 on chromosome 12) have been identified to cause dwarfism across many dog breeds. Some breeds are nearly homozygous for both retrogenes (e.g., Dachshunds) and others are homozygous for just one (e.g., Beagles and Scottish Terriers). Since most breeds do not segregate both of these retrogenes, it is challenging to evaluate their individual effects on long bone length and body size. We identified two dog breeds selected for hunting ability, the Alpine Dachsbracke and the Schweizer Niederlaufhund, that segregate both of these retrogenes. Using individual measurements of height at the shoulder, back length, head width, thorax depth and width, and thoracic limb measurements, we evaluated the combined effects of FGF4 retrogenes within these breeds. We applied multivariable linear regression analysis to determine the effects of retrogene copy numbers on the measurements. Copy numbers of both retrogenes had significant effects reducing height at the shoulders and antebrachial length, with FGF4L1 having a much greater effect than FGF4L2. FGF4L1 alone influenced the degree of carpal valgus and FGF4L2 alone increased head width. Neither retrogene had an effect on thorax width or depth. Selectively breeding dogs with FGF4L1 and without FGF4L2 would likely lead to a reduction in the FGF4L2-related risk of intervertebral disc herniation while maintaining the reduction in leg length resulting from FGF4L1.

Genetic testing of dogs predicts problem behaviors in clinical and nonclinical samples

BACKGROUND

Very little is known about the etiology of personality and psychiatric disorders. Because the core neurobiology of many such traits is evolutionarily conserved, dogs present a powerful model. We previously reported genome scans of breed averages of ten traits related to fear, anxiety, aggression and social behavior in multiple cohorts of pedigree dogs. As a second phase of that discovery, here we tested the ability of markers at 13 of those loci to predict canine behavior in a community sample of 397 pedigree and mixed-breed dogs with individual-level genotype and phenotype data.

RESULTS

We found support for all markers and loci. By including 122 dogs with veterinary behavioral diagnoses in our cohort, we were able to identify eight loci associated with those diagnoses. Logistic regression models showed subsets of those loci could predict behavioral diagnoses. We corroborated our previous findings that small body size is associated with many problem behaviors and large body size is associated with increased trainability. Children in the home were associated with anxiety traits; illness and other animals in the home with coprophagia; working-dog status with increased energy and separation-related problems; and competitive dogs with increased aggression directed at familiar dogs, but reduced fear directed at humans and unfamiliar dogs. Compared to other dogs, Pit Bull-type dogs were not defined by a set of our markers and were not more aggressive; but they were strongly associated with pulling on the leash. Using severity-threshold models, Pit Bull-type dogs showed reduced risk of owner-directed aggression (75th quantile) and increased risk of dog-directed fear (95th quantile).

CONCLUSIONS

Our association analysis in a community sample of pedigree and mixed-breed dogs supports the interbreed mapping. The modeling shows some markers are predictive of behavioral diagnoses. Our findings have broad utility, including for clinical and breeding purposes, but we caution that thorough understanding is necessary for their interpretation and use.

The construction of a haplotype reference panel using extremely low coverage whole genome sequences and its application in genome-wide association studies and genomic prediction in Duroc pigs

Extremely low coverage whole genome sequencing (lcWGS) is an economical technique to obtain high-density single nucleotide polymorphisms (SNPs). Here, we explored the feasibility of constructing a haplotype reference panel (lcHRP) using lcWGS and evaluated the effects of lcHRP through a genome-wide association study (GWAS) and genomic prediction in pigs. A total of 297 and 974 Duroc pigs were genotyped using lcWGS and a 50 K SNP array, respectively. We obtained 19,306,498 SNPs using lcWGS with an accuracy of 0.984. With the help of lcHRP, the accuracy of imputation from the SNP array to lcWGS was 0.922. Compared to the SNP array findings, those from the imputation-based GWAS identified more signals across four traits. With the integration of the top 1% imputation-based GWAS findings as genomic features, the accuracies of genomic prediction was improved by 6.0% to 13.2%. This study showed the great potential of lcWGS in pigs' molecular breeding.

Genomic Prediction of Two Complex Orthopedic Traits Across Multiple Pure and Mixed Breed Dogs

Canine hip dysplasia (CHD) and rupture of the cranial cruciate ligament (RCCL) are two complex inherited orthopedic traits of dogs. These two traits may occur concurrently in the same dog. Genomic prediction of these two diseases would benefit veterinary medicine, the dog’s owner, and dog breeders because of their high prevalence, and because both traits result in painful debilitating osteoarthritis in affected joints. In this study, 842 unique dogs from 6 breeds with hip and stifle phenotypes were genotyped on a customized Illumina high density 183 k single nucleotide polymorphism (SNP) array and also analyzed using an imputed dataset of 20,487,155 SNPs. To implement genomic prediction, two different statistical methods were employed: Genomic Best Linear Unbiased Prediction (GBLUP) and a Bayesian method called BayesC. The cross-validation results showed that the two methods gave similar prediction accuracy (r = 0.3–0.4) for CHD (measured as Norberg angle) and RCCL in the multi-breed population. For CHD, the average correlation of the AUC was 0.71 (BayesC) and 0.70 (GBLUP), which is a medium level of prediction accuracy and consistent with Pearson correlation results. For RCCL, the correlation of the AUC was slightly higher. The prediction accuracy of GBLUP from the imputed genotype data was similar to the accuracy from DNA array data. We demonstrated that the genomic prediction of CHD and RCCL with DNA array genotype data is feasible in a multiple breed population if there is a genetic connection, such as breed, between the reference population and the validation population. Albeit these traits have heritability of about one-third, higher accuracy is needed to implement in a natural population and predicting a complex phenotype will require much larger number of dogs within a breed and across breeds. It is possible that with higher accuracy, genomic prediction of these orthopedic traits could be implemented in a clinical setting for early diagnosis and treatment, and the selection of dogs for breeding. These results need continuous improvement in model prediction through ongoing genotyping and data sharing. When genomic prediction indicates that a dog is susceptible to one of these orthopedic traits, it should be accompanied by clinical and radiographic screening at an acceptable age with appropriate follow-up.

Best practices for analyzing imputed genotypes from low-pass sequencing in dogs

Although DNA array-based approaches for genome-wide association studies (GWAS) permit the collection of thousands of low-cost genotypes, it is often at the expense of resolution and completeness, as SNP chip technologies are ultimately limited by SNPs chosen during array development. An alternative low-cost approach is low-pass whole genome sequencing (WGS) followed by imputation. Rather than relying on high levels of genotype confidence at a set of select loci, low-pass WGS and imputation rely on the combined information from millions of randomly sampled low-confidence genotypes. To investigate low-pass WGS and imputation in the dog, we assessed accuracy and performance by downsampling 97 high-coverage (> 15×) WGS datasets from 51 different breeds to approximately 1× coverage, simulating low-pass WGS. Using a reference panel of 676 dogs from 91 breeds, genotypes were imputed from the downsampled data and compared to a truth set of genotypes generated from high-coverage WGS. Using our truth set, we optimized a variant quality filtering strategy that retained approximately 80% of 14 M imputed sites and lowered the imputation error rate from 3.0% to 1.5%. Seven million sites remained with a MAF > 5% and an average imputation quality score of 0.95. Finally, we simulated the impact of imputation errors on outcomes for case-control GWAS, where small effect sizes were most impacted and medium-to-large effect sizes were minorly impacted. These analyses provide best practice guidelines for study design and data post-processing of low-pass WGS-imputed genotypes in dogs.

Whole Genome Sequencing Reveals Multiple Linked Genetic Variants on Canine Chromosome 12 Associated with Risk for Symmetrical Lupoid Onychodystrophy (SLO) in the Bearded Collie

In dogs, symmetrical lupoid onychodystrophy (SLO) results in nail loss and an abnormal regrowth of the claws. In Bearded Collies, an autoimmune nature has been suggested because certain dog leukocyte antigen (DLA) class II haplotypes are associated with the condition. A genome-wide association study of the Bearded Collie revealed two regions of association that conferred risk for disease: one on canine chromosome (CFA) 12 that encompasses the DLA genes, and one on CFA17. Case-control association was employed on whole genome sequencing data to uncover putative causative variants in SLO within the CFA12 and CFA17 associated regions. Genotype imputation was then employed to refine variants of interest. Although no SLO-associated protein-coding variants were identified on CFA17, multiple variants, many with predicted damaging effects, were identified within potential candidate genes on CFA12. Furthermore, many potentially damaging alleles were fully correlated with the presence of DLA class II risk haplotypes for SLO, suggesting that the variants may reflect DLA class II haplotype association with disease or vice versa. Strong linkage disequilibrium in the region precluded the ability to isolate and assess the individual or combined effect of variants on disease development. Nonetheless, all were predictive of risk for SLO and, with judicious assessment, their application in selective breeding may prove useful to reduce the incidence of SLO in the breed.

Integration of selection signatures and multi-trait GWAS reveals polygenic genetic architecture of carcass traits in beef cattle

Carcass merits are widely considered as economically important traits affecting beef production in the beef cattle industry. However, the genetic basis of carcass traits remains to be well understood. Here, we applied multiple methods, including the Composite of Likelihood Ratio (CLR) and Genome-wide Association Study (GWAS), to explore the selection signatures and candidate variants affecting carcass traits. We identified 11,600 selected regions overlapping with 2214 candidate genes, and most of those were enriched in binding and gene regulation. Notably, we identified 66 and 110 potential variants significantly associated with carcass traits using single-trait and multi-traits analyses, respectively. By integrating selection signatures with single and multi-traits associations, we identified 12 and 27 putative genes, respectively. Several highly conserved missense variants were identified in OR5M13D, NCAPG, and TEX2. Our study supported polygenic genetic architecture of carcass traits and provided novel insights into the genetic basis of complex traits in beef cattle.

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

Genotype imputation using a reference panel that combines high-density array data and publicly available whole genome sequence consortium variant data is potentially a cost-effective method to increase the density of extant lower-density array datasets. In this study, three datasets (two Border Collie; one Italian Spinone) generated using a legacy array (Illumina CanineHD, 173 662 SNPs) were utilised to assess the feasibility and accuracy of this approach and to gather additional evidence for the efficacy of canine genotype imputation. The cosmopolitan reference panels used to impute genotypes comprised dogs of 158 breeds, mixed breed dogs, wolves and Chinese indigenous dogs, as well as breed-specific individuals genotyped using the Axiom Canine HD array. The two Border Collie reference panels comprised 808 individuals including 79 Border Collies and 426 326 or 426 332 SNPs; and the Italian Spinone reference panel comprised 807 individuals including 38 Italian Spinoni and 476 313 SNPs. A high accuracy for imputation was observed, with the lowest accuracy observed for one of the Border Collie datasets (mean R² = 0.94) and the highest for the Italian Spinone dataset (mean R² = 0.97). This study’s findings demonstrate that imputation of a legacy array study set using a reference panel comprising both breed-specific array data and multi-breed variant data derived from whole genomes is effective and accurate. The process of canine genotype imputation, using the valuable growing resource of publicly available canine genome variant datasets alongside breed-specific data, is described in detail to facilitate and encourage use of this technique in canine genetics.

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.

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.

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.

Identification of Genetic Susceptibility Factors Associated with Canine Gastric Dilatation-Volvulus

Canine gastric dilatation-volvulus (GDV) is a common life-threatening condition occurring primarily in large and giant breeds with a 3.9% to 36.7% lifetime risk. The genetic correlates of GDV have not previously been systematically explored. We undertook an inter-breed genome-wide association analysis (GWAS) of 253 dogs from ten breeds including 106 healthy dogs and 147 dogs with at least one GDV episode. SNP array genotyping followed by imputation was conducted on 241 samples to identify GDV-associated single-nucleotide polymorphisms (SNPs) and copy number variations (CNVs). A subset of 33 dogs (15 healthy dogs and 18 GDV patients from the three most represented breeds) was characterized by whole genome sequencing (WGS). After genome-wide Bonferroni correction, we identified a significant putatively protective intergenic SNP (rs851737064) across all breeds. The signal was most significant in Collies, German Shorthaired Pointers, and Great Danes. Subsequent focused analysis across these three breeds identified 12 significant additional putatively protective or deleterious SNPs. Notable significant SNPs included those occurring in genes involved in gastric tone and motility including VHL, NALCN, and PRKCZ. These data provide important new clues to canine GDV risk factors and facilitate generation of hypotheses regarding the genetic and molecular underpinnings this syndrome.

Variants in FtsJ RNA 2'-O-Methyltransferase 3 and Growth Hormone 1 are associated with small body size and a dental anomaly in dogs

Domesticated dogs show unparalleled diversity in body size across breeds, but within breeds variation is limited by selective breeding. Many heritable diseases of dogs are found among breeds of similar sizes, suggesting that as in humans, alleles governing growth have pleiotropic effects. Here, we conducted independent genome-wide association studies in the small Shetland Sheepdog breed and discovered a locus on chromosome 9 that is associated with a dental abnormality called maxillary canine-tooth mesioversion (MCM) (P = 1.53 × 10^-7) as well as two body size traits: height (P = 1.67 × 10^-5) and weight (P = 1.16 × 10^-7). Using whole-genome resequencing data, we identified variants in two proximal genes: FTSJ3, encoding an RNA methyltransferase, and GH1, encoding growth hormone. A substitution in FTSJ3 and a splice donor insertion in GH1 are strongly associated with MCM and reduced body size in Shetland Sheepdogs. We demonstrated in vitro that the GH1 variant leads to exon 3 skipping, predicting a mutant protein known to cause human pituitary dwarfism. Statistical modeling, however, indicates that the FTSJ3 variant is the stronger predictor of MCM and that each derived allele reduces body size by about 1 inch and 5 pounds. In a survey of 224 breeds, both FTSJ3 and GH1 variants are frequent among very small "toy" breeds and absent from larger breeds. Our findings indicate that a chromosome 9 locus harboring tightly linked variants in FTSJ3 and GH1 reduces growth in the Shetland Sheepdog and toy breed dogs and confers risk for MCM through vertical pleiotropy.

Current Understanding of the Genetics of Intervertebral Disc Degeneration

Premature degeneration of the intervertebral disc and its association with specific chondrodystrophic dog breeds has been recognized for over a century. Several lines of evidence including disease breed predisposition, studies suggesting heritability of premature intervertebral disc degeneration (IVDD) and association of a dog chromosome 12 (CFA 12) locus with intervertebral disc calcification have strongly supported a genetic component in IVDD in dogs. Recent studies documenting association of IVDD with an overexpressing FGF4 retrogene on CFA 12 have opened up new areas of investigation to further define the pathophysiology of premature IVDD. While preliminary data from studies investigating FGF4 retrogenes in IVDD implicate FGF4 overexpression as a major disease factor, they have also highlighted knowledge gaps in our understanding of intervertebral disc herniation which is a complex and multifactorial disease process.

SNP-based heritability and genetic architecture of cranial cruciate ligament rupture in Labrador Retrievers

Cranial cruciate ligament rupture (CCLR) is one of the leading causes of pelvic limb lameness in dogs. About 6% of Labrador Retrievers suffer from this orthopedic problem. The aim of this study was to determine the heritability of CCLR in this breed using SNP array genotyping data. DNA samples were collected from CCLR-affected dogs (n = 190) and unaffected dogs over the age of 8 years (n = 143). All 333 dogs were genotyped directly or imputed up to approximately 710k SNPs on the Affymetrix Axiom CanineHD SNP array. Heritability of CCLR was calculated using multiple methodologies, including linear mixed models, Bayesian models and a model that incorporates LD. The covariates of sex and sterilization status were added to each analysis to assess their impact. Across the algorithms of these models, heritability ranged from 0.550 to 0.886, depending on covariate inclusion. The relatively high heritability for this disease indicates that a substantial genetic component contributes to CCLR in the Labrador Retriever.

Association of Common Genetic Variants in the CPSF7 and SDHAF2 Genes with Canine Idiopathic Pulmonary Fibrosis in the West Highland White Terrier

Canine idiopathic pulmonary fibrosis (CIPF) is a chronic fibrotic lung disease that is observed at a higher frequency in the West Highland White Terrier dog breed (WHWT) and may have molecular pathological overlap with human lung fibrotic disease. We conducted a genome-wide association study (GWAS) in the WHWT using whole genome sequencing (WGS) to discover genetic variants associated with CIPF. Saliva-derived DNA samples were sequenced using the Riptide DNA library prep kit. After quality controls, 28 affected, 44 unaffected, and 1,843,695 informative single nucleotide polymorphisms (SNPs) were included in the GWAS. Data were analyzed both at the single SNP and gene levels using the GEMMA and GATES methods, respectively. We detected significant signals at the gene level in both the cleavage and polyadenylation specific factor 7 (CPSF7) and succinate dehydrogenase complex assembly factor 2 (SDHAF2) genes (adjusted p = 0.016 and 0.024, respectively), two overlapping genes located on chromosome 18. The top SNP for both genes was rs22669389; however, it did not reach genome-wide significance in the GWAS (adjusted p = 0.078). Our studies provide, for the first time, candidate loci for CIPF in the WHWT. CPSF7 was recently associated with lung adenocarcinoma, further highlighting the potential relevance of our results because IPF and lung cancer share several pathological mechanisms.

Genetic Variants Affecting Skeletal Morphology in Domestic Dogs

Purebred dog breeds provide a powerful resource for the discovery of genetic variants affecting skeletal morphology. Domesticated and subsequently purebred dogs have undergone strong artificial selection for a broad range of skeletal variation, which include both the size and shapes of their bones. While the phenotypic variation between breeds is high, within-breed morphological variation is typically low. Approaches for defining genetic variants associated with canine morphology include quantitative within-breed analyses, as well as across-breed analyses, using breed standards as proxies for individual measurements. The ability to identify variants across the genomes of individual dogs can now be paired with precise measures of morphological variation to define the genetic interactions and the phenotypic effect of variants on skeletal morphology.