Genetic Variants Affecting Skeletal Morphology in Domestic Dogs

Large-scale genomic analysis of the domestic dog informs biological discovery

Recent advances in genomics, coupled with a unique population structure and remarkable levels of variation, have propelled the domestic dog to new levels as a system for understanding fundamental principles in mammalian biology. Central to this advance are more than 350 recognized breeds, each a closed population that has undergone selection for unique features. Genetic variation in the domestic dog is particularly well characterized compared with other domestic mammals, with almost 3000 high-coverage genomes publicly available. Importantly, as the number of sequenced genomes increases, new avenues for analysis are becoming available. Herein, we discuss recent discoveries in canine genomics regarding behavior, morphology, and disease susceptibility. We explore the limitations of current data sets for variant interpretation, tradeoffs between sequencing strategies, and the burgeoning role of long-read genomes for capturing structural variants. In addition, we consider how large-scale collections of whole-genome sequence data drive rare variant discovery and assess the geographic distribution of canine diversity, which identifies Asia as a major source of missing variation. Finally, we review recent comparative genomic analyses that will facilitate annotation of the noncoding genome in dogs.

Identification of consensus homozygous regions and their associations with growth and feed efficiency traits in American mink

The recent chromosome-based genome assembly and the newly developed 70K single nucleotide polymorphism (SNP) array for American mink (Neogale vison) facilitate the identification of genetic variants underlying complex traits in this species. The objective of this study was to evaluate the association between consensus runs of homozygosity (ROH) with growth and feed efficiency traits in American mink. A subsample of two mink populations (n = 2,986) were genotyped using the Affymetrix Mink 70K SNP array. The identified ROH segments were included simultaneously, concatenated into consensus regions, and the ROH-based association studies were carried out with linear mixed models considering a genomic relationship matrix for 11 growth and feed efficiency traits implemented in ASReml-R version 4. In total, 298,313 ROH were identified across all individuals, with an average length and coverage of 4.16 Mb and 414.8 Mb, respectively. After merging ROH segments, 196 consensus ROH regions were detected and used for genome-wide ROH-based association analysis. Thirteen consensus ROH regions were significantly (P < 0.01) associated with growth and feed efficiency traits. Several candidate genes within the significant regions are known for their involvement in growth and body size development, including MEF2A, ADAMTS17, POU3F2, and TYRO3. In addition, we found ten consensus ROH regions, defined as ROH islands, with frequencies over 80% of the population. These islands harbored 12 annotated genes, some of which were related to immune system processes such as DTX3L, PARP9, PARP14, CD86, and HCLS1. This is the first study to explore the associations between homozygous regions with growth and feed efficiency traits in American mink. Our findings shed the light on the effects of homozygosity in the mink genome on growth and feed efficiency traits, that can be utilized in developing a sustainable breeding program for mink.

Sleep-physiological correlates of brachycephaly in dogs

The shape of the cranium is one of the most notable physical changes induced in domestic dogs through selective breeding and is measured using the cephalic index (CI). High CI (a ratio of skull width to skull length > 60) is characterized by a short muzzle and flat face and is referred to as brachycephaly. Brachycephalic dogs display some potentially harmful changes in neuroanatomy, and there are implications for differences in behavior, as well. The path from anatomy to cognition, however, has not been charted in its entirety. Here, we report that sleep-physiological markers of white-matter loss (high delta power, low frontal spindle frequency, i.e., spindle waves/s), along with a spectral profile for REM (low beta, high delta) associated with low intelligence in humans, are each linked to higher CI values in the dog. Additionally, brachycephalic subjects spent more time sleeping, suggesting that the sleep apnea these breeds usually suffer from increases daytime sleepiness. Within sleep, more time was spent in the REM sleep stage than in non-REM, while REM duration was correlated positively with the number of REM episodes across dogs. It is currently not clear if the patterns of sleep and sleep-stage duration are mainly caused by sleep-impairing troubles in breathing and thermoregulation, present a juvenile-like sleeping profile, or are caused by neuro-psychological conditions secondary to the effects of brachycephaly, e.g., frequent REM episodes are known to appear in human patients with depression. While future studies should more directly address the interplay of anatomy, physiology, and behavior within a single experiment, this represents the first description of how the dynamics of the canine brain covary with CI, as measured in resting companion dogs using a non-invasive sleep EEG methodology. The observations suggest that the neuroanatomical changes accompanying brachycephaly alter neural systems in a way that can be captured in the sleep EEG, thus supporting the utility of the latter in the study of canine brain health and function.

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.

The Revolution of Animal Genomics in Forensic Sciences

Nowadays, the coexistence between humans and domestic animals (especially dogs and cats) has become a common scenario of daily life. Consequently, during a forensic investigation in civil or criminal cases, the biological material from a domestic animal could be considered "evidence" by law enforcement agencies. Animal genomics offers an important contribution in attacks and episodes of property destruction or in a crime scene where the non-human biological material is linked to the victim or perpetrator. However, only a few animal genetics laboratories in the world are able to carry out a valid forensic analysis, adhering to standards and guidelines that ensure the admissibility of data before a court of law. Today, forensic sciences focus on animal genetics considering all domestic species through the analysis of STRs (short tandem repeats) and autosomal and mitochondrial DNA SNPs (single nucleotide polymorphisms). However, the application of these molecular markers to wildlife seems to have gradually gained a strong relevance, aiming to tackle illegal traffic, avoid the loss of biodiversity, and protect endangered species. The development of third-generation sequencing technologies has glimmered new possibilities by bringing "the laboratory into the field", with a reduction of both the enormous cost management of samples and the degradation of the biological material.

Transcriptomics and Selection Pressure Analysis Reveals the Influence Mechanism of PLIN1 Protein on the Development of Small Size in Min Pigs

Body size is an important biological phenotypic trait that has attracted substantial attention. Small domestic pigs can serve as excellent animal models for biomedicine and also help meet sacrificial culture needs in human societies. Although the mechanisms underlying vertebral development regulating body size variation in domestic pigs during the embryonic period have been well described, few studies have examined the genetic basis of body size variation in post embryonic developmental stages. In this study, seven candidate genes-PLIN1, LIPE, PNPLA1, SCD, FABP5, KRT10 and IVL-significantly associated with body size were identified in Min pigs, on the basis of weighted gene co-expression network analysis (WGCNA), and most of their functions were found to be associated with lipid deposition. Six candidate genes except for IVL were found to have been subjected to purifying selection. PLIN1 had the lowest ω value (0.139) and showed heterogeneous selective pressure among domestic pig lineages with different body sizes (p < 0.05). These results suggested that PLIN1 is an important genetic factor regulating lipid deposition and consequently affecting body size variation in pigs. The culture of whole pig sacrifice in Manchu during the Qing Dynasty in China might have contributed to the strong artificial domestication and selection of Hebao pigs.

That brachycephalic look: Infant-like facial appearance in short-muzzled dog breeds

Brachycephalic dog breeds are highly popular, yet their conformation-related disorders represent a major welfare concern. It has been suggested that the current popularity of such breeds can be explained by their cute, infant-like facial appearances. The concept of 'kindchenschema' refers to the observation that certain physical features of infant humans and other animals can automatically stimulate positive and nurturant feelings in adult observers. But the proposal that brachycephalic dogs possess heightened 'kindchenschema' facial features, even into adulthood, has never been formally investigated. Here, we hypothesised that relative muzzle shortening across a range of breeds would be associated with known 'kindchenschema' facial features, including a relatively larger forehead, larger eyes and smaller nose. Relative fronto-facial feature sizes in exemplar photographs of adult dogs from 42 popular breeds were measured and associated with existing data on the relative muzzle length and height-at-withers of the same breeds. Our results show that, in adulthood, shorter-muzzled breeds have relatively larger (taller) foreheads and relatively larger eyes (i.e. area of exposed eyeball relative to overall face area) than longer-muzzled breeds, and that this effect is independent of breed size. In sum, brachycephalic dog breeds do show exaggeration of some, but not all, known fronto-facial 'kindchenschema' features, and this may well contribute to their apparently cute appearance and to their current popularity as companion animals. We conclude that the challenge of addressing conformation-related disorders in companion dogs needs to take account of the cute, 'kindchenschema' looks that many owners are likely to be attracted to.

PCYT1A Missense Variant in Vizslas with Disproportionate Dwarfism

Disproportionate dwarfism phenotypes represent a heterogeneous subset of skeletal dysplasias and have been described in many species including humans and dogs. In this study, we investigated Vizsla dogs that were affected by disproportionate dwarfism that we propose to designate as skeletal dysplasia 3 (SD3). The most striking skeletal changes comprised a marked shortening and deformation of the humerus and femur. An extended pedigree with six affected dogs suggested autosomal recessive inheritance. Combined linkage and homozygosity mapping localized a potential genetic defect to a ~4 Mb interval on chromosome 33. We sequenced the genome of an affected dog, and comparison with 926 control genomes revealed a single, private protein-changing variant in the critical interval, PCYT1A:XM_038583131.1:c.673T>C, predicted to cause an exchange of a highly conserved amino acid, XP_038439059.1:p.(Y225H). We observed perfect co-segregation of the genotypes with the phenotype in the studied family. When genotyping additional Vizslas, we encountered a single dog with disproportionate dwarfism that did not carry the mutant PCYT1A allele, which we hypothesize was due to heterogeneity. In the remaining 130 dogs, we observed perfect genotype-phenotype association, and none of the unaffected dogs were homozygous for the mutant PCYT1A allele. PCYT1A loss-of-function variants cause spondylometaphyseal dysplasia with cone-rod dystrophy (SMD-CRD) in humans. The skeletal changes in Vizslas were comparable to human patients. So far, no ocular phenotype has been recognized in dwarf Vizslas. We propose the PCYT1A missense variant as a candidate causative variant for SD3. Our data facilitate genetic testing of Vizslas to prevent the unintentional breeding of further affected puppies.

Domestic dog lineages reveal genetic drivers of behavioral diversification

Selective breeding of domestic dogs has generated diverse breeds often optimized for performing specialized tasks. Despite the heritability of breed-typical behavioral traits, identification of causal loci has proven challenging due to the complexity of canine population structure. We overcome longstanding difficulties in identifying genetic drivers of canine behavior by developing a framework for understanding relationships between breeds and the behaviors that define them, utilizing genetic data for over 4,000 domestic, semi-feral, and wild canids and behavioral survey data for over 46,000 dogs. We identify ten major canine genetic lineages and their behavioral correlates and show that breed diversification is predominantly driven by non-coding regulatory variation. We determine that lineage-associated genes converge in neurodevelopmental co-expression networks, identifying a sheepdog-associated enrichment for interrelated axon guidance functions. This work presents a scaffold for canine diversification that positions the domestic dog as an unparalleled system for revealing the genetic origins of behavioral diversity.

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.

Forensic Tools for Species Identification of Skeletal Remains: Metrics, Statistics, and OsteoID

Although nonhuman remains constitute a significant portion of forensic anthropological casework, the potential use of bone metrics to assess the human origin and to classify species of skeletal remains has not been thoroughly investigated. This study aimed to assess the utility of quantitative methods in distinguishing human from nonhuman remains and present additional resources for species identification. Over 50,000 measurements were compiled from humans and 27 nonhuman (mostly North American) species. Decision trees developed from the long bone data can differentiate human from nonhuman remains with over 90% accuracy (>98% accuracy for the human sample), even if all long bones are pooled. Stepwise discriminant function results were slightly lower (>87.4% overall accuracy). The quantitative models can be used to support visual identifications or preliminarily assess forensic significance at scenes. For species classification, bone-specific discriminant functions returned accuracies between 77.7% and 89.1%, but classification results varied highly across species. From the study data, we developed a web tool, OsteoID, for users who can input measurements and be shown photographs of potential bones/species to aid in visual identification. OsteoID also includes supplementary images (e.g., 3D scans), creating an additional resource for forensic anthropologists and others involved in skeletal species identification and comparative osteology.

LTBP3 Frameshift Variant in British Shorthair Cats with Complex Skeletal Dysplasia

We investigated a highly inbred family of British Shorthair cats in which two offspring were affected by deteriorating paraparesis due to complex skeletal malformations. Radiographs of both affected kittens revealed vertebral deformations with marked stenosis of the vertebral canal from T11 to L3. Additionally, compression of the spinal cord, cerebellar herniation, coprostasis and hypogangliosis were found. The pedigree suggested monogenic autosomal recessive inheritance of the trait. We sequenced the genome of an affected kitten and compared the data to 62 control genomes. This search yielded 55 private protein-changing variants of which only one was located in a likely functional candidate gene, LTBP3, encoding latent transforming growth factor β binding protein 3. This variant, c.158delG or p.(Gly53Alafs*16), represents a 1 bp frameshift deletion predicted to truncate 95% of the open reading frame. LTBP3 is a known key regulator of transforming growth factor β (TGF-β) and is involved in bone morphogenesis and remodeling. Genotypes at the LTBP3:c.158delG variant perfectly co-segregated with the phenotype in the investigated family. The available experimental data together with current knowledge on LTBP3 variants and their functional impact in human patients and mice suggest LTBP3:c.158delG as a candidate causative variant for the observed skeletal malformations in British Shorthair cats. To the best of our knowledge, this study represents the first report of LTBP3-related complex skeletal dysplasia in domestic animals.

Insight into the Candidate Genes and Enriched Pathways Associated with Height, Length, Length to Height Ratio and Body-Weight of Korean Indigenous Breed, Jindo Dog Using Gene Set Enrichment-Based GWAS Analysis

As a companion and hunting dog, height, length, length to height ratio (LHR) and body-weight are the vital economic traits for Jindo dog. Human selection and targeted breeding have produced an extraordinary diversity in these traits. Therefore, the identification of causative markers, genes and pathways that help us to understand the genetic basis of this variability is essential for their selection purposes. Here, we performed a genome-wide association study (GWAS) combined with enrichment analysis on 757 dogs using 118,879 SNPs. The genomic heritability (h²) was 0.33 for height and 0.28 for weight trait in Jindo. At p-value < 5 × 10^-5, ten, six, thirteen and eleven SNPs on different chromosomes were significantly associated with height, length, LHR and body-weight traits, respectively. Based on our results, HHIP, LCORL and NCAPG for height, IGFI and FGFR3 for length, DLK1 and EFEMP1 for LHR and PTPN2, IGFI and RASAL2 for weight can be the potential candidate genes because of the significant SNPs located in their intronic or upstream regions. The gene-set enrichment analysis highlighted here nine and seven overlapping significant (p < 0.05) gene ontology (GO) terms and pathways among traits. Interestingly, the highlighted pathways were related to hormone synthesis, secretion and signalling were generally involved in the metabolism, growth and development process. Our data provide an insight into the significant genes and pathways if verified further, which will have a significant effect on the breeding of the Jindo dog's population.

PRKG2 Splice Site Variant in Dogo Argentino Dogs with Disproportionate Dwarfism

Dwarfism phenotypes occur in many species and may be caused by genetic or environmental factors. In this study, we investigated a family of nine Dogo Argentino dogs, in which two dogs were affected by disproportionate dwarfism. Radiographs of an affected dog revealed a decreased level of endochondral ossification in its growth plates, and a premature closure of the distal ulnar physes. The pedigree of the dogs presented evidence of monogenic autosomal recessive inheritance; combined linkage and homozygosity mapping assigned the most likely position of a potential genetic defect to 34 genome segments, totaling 125 Mb. The genome of an affected dog was sequenced and compared to 795 control genomes. The prioritization of private variants revealed a clear top candidate variant for the observed dwarfism. This variant, PRKG2:XM_022413533.1:c.1634+1G>T, affects the splice donor site and is therefore predicted to disrupt the function of the PKRG2 gene encoding protein, kinase cGMP-dependent type 2, a known regulator of chondrocyte differentiation. The genotypes of the PRKG2 variant were perfectly associated with the phenotype in the studied family of dogs. PRKG2 loss-of-function variants were previously reported to cause disproportionate dwarfism in humans, cattle, mice, and rats. Together with the comparative data from other species, our data strongly suggest PRKG2:c.1634+1G>T to be a candidate causative variant for the observed dwarfism phenotype in Dogo Argentino dogs.

Genomic insights into body size evolution in Carnivora support Peto's paradox

Background

The range of body sizes in Carnivora is unparalleled in any other mammalian order—the heaviest species is 130,000 times heavier than the lightest and the longest species is 50 times longer than the shortest. However, the molecular mechanisms underlying these huge differences in body size have not been explored.

Results

Herein, we performed a comparative genomics analysis of 20 carnivores to explore the evolutionary basis of the order’s great variations in body size. Phylogenetic generalized least squares (PGLS) revealed that 337 genes were significantly related to both head body length and body mass; these genes were defined as body size associated genes (BSAGs). Fourteen positively-related BSAGs were found to be associated with obesity, and three of these were under rapid evolution in the extremely large carnivores, suggesting that these obesity-related BSAGs might have driven the body size expansion in carnivores. Interestingly, 100 BSAGs were statistically significantly enriched in cancer control in carnivores, and 15 of which were found to be under rapid evolution in extremely large carnivores. These results suggested that large carnivores might have evolved an effective mechanism to resist cancer, which could be regarded as molecular evidence to support Peto’s paradox. For small carnivores, we identified 15 rapidly evolving genes and found six genes with fixed amino acid changes that were reported to reduce body size.

Conclusions

This study brings new insights into the molecular mechanisms that drove the diversifying evolution of body size in carnivores, and provides new target genes for exploring the mysteries of body size evolution in mammals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07732-w.

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.