Insights into morphology and disease from the dog genome project.

Although most modern dog breeds are less than 200 years old, the symbiosis between man and dog is ancient. Since prehistoric times, repeated selection events have transformed the wolf into man's guardians, laborers, athletes, and companions. The rapid transformation from pack predator to loyal companion is a feat that is arguably unique among domesticated animals. How this transformation came to pass remained a biological mystery until recently: Within the past decade, the deployment of genomic approaches to study population structure, detect signatures of selection, and identify genetic variants that underlie canine phenotypes is ushering into focus novel biological mechanisms that make dogs remarkable. Ironically, the very practices responsible for breed formation also spurned morbidity; today, many diseases are correlated with breed identity. In this review, we discuss man's best friend in the context of a genetic model to understand paradigms of heritable phenotypes, both desirable and disadvantageous.

Insights into digit evolution from a fate map study of the forearm using Chameleon, a new transgenic chicken line.

The cellular and genetic networks that contribute to the development of the zeugopod (radius and ulna of the forearm, tibia and fibula of the leg) are not well understood, although these bones are susceptible to loss in congenital human syndromes and to the action of teratogens such as thalidomide. Using a new fate-mapping approach with the Chameleon transgenic chicken line, we show that there is a small contribution of SHH-expressing cells to the posterior ulna, posterior carpals and digit 3. We establish that although the majority of the ulna develops in response to paracrine SHH signalling in both the chicken and mouse, there are differences in the contribution of SHH-expressing cells between mouse and chicken as well as between the chicken ulna and fibula. This is evidence that, although zeugopod bones are clearly homologous according to the fossil record, the gene regulatory networks that contribute to their development and evolution are not fixed.

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).

The impact of the COVID-19 pandemic on a cohort of Labrador retrievers in England.

BACKGROUND: The COVID-19 pandemic is likely to have affected the welfare and health of dogs due to surges in adoptions and purchases, changes in the physical and mental health and financial status of dog owners, changes in dogs' lifestyle and routines and limited access to veterinary care. The aims of this study were to investigate whether COVID-19 restrictions were associated with differences in Labrador retrievers' lifestyle, routine care, insurance status, illness incidence or veterinary attendance with an illness, who were living in England and enrolled in Dogslife, an owner-based cohort study. Longitudinal questionnaire data from Dogslife that was relevant to the dates between the 23rd of March and the 4th of July 2020, during COVID-19 restrictions in England, were compared to data between the same dates in previous years from 2011 to 2019 using mixed regression models and adjusted chi-squared tests. RESULTS: Compared with previous years (March 23rd to July 4th, 2010 to 2019), the COVID-19 restrictions study period (March 23rd to July 4th 2020) was associated with owners reporting increases in their dogs' exercise and worming and decreases in insurance, titbit-feeding and vaccination. Odds of owners reporting that their dogs had an episode of coughing (0.20, 95% CI: 0.04-0.92) and that they took their dogs to a veterinarian with an episode of any illness (0.58, 95% CI: 0.45-0.76) were lower during the COVID-19 restrictions compared to before. During the restrictions period, owners were less likely to report that they took their dogs to a veterinarian with certain other illnesses, compared to before this period. CONCLUSIONS: Dogslife provided a unique opportunity to study prospective questionnaire data from owners already enrolled on a longitudinal cohort study. This approach minimised bias associated with recalling events prior to the pandemic and allowed a wider population of dogs to be studied than is available from primary care data. Distinctive insights into owners' decision making about their dogs' healthcare were offered. There are clear implications of the COVID-19 pandemic and associated restrictions for the lifestyle, care and health of dogs.

An ADAMTS3 missense variant is associated with Norwich Terrier upper airway syndrome.

In flat-faced dog breeds, air resistance caused by skull conformation is believed to be a major determinant of Brachycephalic Obstructive Airway Syndrome (BOAS). The clinical presentation of BOAS is heterogeneous, suggesting determinants independent of skull conformation contribute to airway disease. Norwich Terriers, a mesocephalic breed, are predisposed to Upper Airway Syndrome (UAS), a disease whose pathological features overlap with BOAS. Our health screening clinic examined and scored the airways of 401 Norwich terriers by laryngoscopy. Genome-wide association analyses of UAS-related pathologies revealed a genetic association on canine chromosome 13 (rs9043975, p = 7.79x10-16). Whole genome resequencing was used to identify causal variant(s) within a 414 kb critical interval. This approach highlighted an error in the CanFam3.1 dog assembly, which when resolved, led to the discovery of a c.2786G>A missense variant in exon 20 of the positional candidate gene, ADAM metallopeptidase with thrombospondin type 1 motif 3 (ADAMTS3). In addition to segregating with UAS amongst Norwich Terriers, the ADAMTS3 c.2786G>A risk allele frequency was enriched among the BOAS-susceptible French and (English) Bulldogs. Previous studies indicate that ADAMTS3 loss of function results in lymphoedema. Our results suggest a new paradigm in the understanding of canine upper airway disease aetiology: airway oedema caused by disruption of ADAMTS3 predisposes dogs to respiratory obstruction. These findings will enhance breeding practices and could refine the prognostics of surgical interventions that are often used to treat airway obstruction.

Host-specialized fibrinogen-binding by a bacterial surface protein promotes biofilm formation and innate immune evasion.

Fibrinogen is an essential part of the blood coagulation cascade and a major component of the extracellular matrix in mammals. The interface between fibrinogen and bacterial pathogens is an important determinant of the outcome of infection. Here, we demonstrate that a canine host-restricted skin pathogen, Staphylococcus pseudintermedius, produces a cell wall-associated protein (SpsL) that has evolved the capacity for high strength binding to canine fibrinogen, with reduced binding to fibrinogen of other mammalian species including humans. Binding occurs via the surface-expressed N2N3 subdomains, of the SpsL A-domain, to multiple sites in the fibrinogen α-chain C-domain by a mechanism analogous to the classical dock, lock, and latch binding model. Host-specific binding is dependent on a tandem repeat region of the fibrinogen α-chain, a region highly divergent between mammals. Of note, we discovered that the tandem repeat region is also polymorphic in different canine breeds suggesting a potential influence on canine host susceptibility to S. pseudintermedius infection. Importantly, the strong host-specific fibrinogen-binding interaction of SpsL to canine fibrinogen is essential for bacterial aggregation and biofilm formation, and promotes resistance to neutrophil phagocytosis, suggesting a key role for the interaction during pathogenesis. Taken together, we have dissected a bacterial surface protein-ligand interaction resulting from the co-evolution of host and pathogen that promotes host-specific innate immune evasion and may contribute to its host-restricted ecology.

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

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

Derived variants at six genes explain nearly half of size reduction in dog breeds.

Selective breeding of dogs by humans has generated extraordinary diversity in body size. A number of multibreed analyses have been undertaken to identify the genetic basis of this diversity. We analyzed four loci discovered in a previous genome-wide association study that used 60,968 SNPs to identify size-associated genomic intervals, which were too large to assign causative roles to genes. First, we performed fine-mapping to define critical intervals that included the candidate genes GHR, HMGA2, SMAD2, and STC2, identifying five highly associated markers at the four loci. We hypothesize that three of the variants are likely to be causative. We then genotyped each marker, together with previously reported size-associated variants in the IGF1 and IGF1R genes, on a panel of 500 domestic dogs from 93 breeds, and identified the ancestral allele by genotyping the same markers on 30 wild canids. We observed that the derived alleles at all markers correlated with reduced body size, and smaller dogs are more likely to carry derived alleles at multiple markers. However, breeds are not generally fixed at all markers; multiple combinations of genotypes are found within most breeds. Finally, we show that 46%-52.5% of the variance in body size of dog breeds can be explained by seven markers in proximity to exceptional candidate genes. Among breeds with standard weights <41 kg (90 lb), the genotypes accounted for 64.3% of variance in weight. This work advances our understanding of mammalian growth by describing genetic contributions to canine size determination in non-giant dog breeds.

Variation of BMP3 contributes to dog breed skull diversity.

Since the beginnings of domestication, the craniofacial architecture of the domestic dog has morphed and radiated to human whims. By beginning to define the genetic underpinnings of breed skull shapes, we can elucidate mechanisms of morphological diversification while presenting a framework for understanding human cephalic disorders. Using intrabreed association mapping with museum specimen measurements, we show that skull shape is regulated by at least five quantitative trait loci (QTLs). Our detailed analysis using whole-genome sequencing uncovers a missense mutation in BMP3. Validation studies in zebrafish show that Bmp3 function in cranial development is ancient. Our study reveals the causal variant for a canine QTL contributing to a major morphologic trait.

The DoGA consortium expression atlas of promoters and genes in 100 canine tissues.

The dog, Canis lupus familiaris, is an important model for studying human diseases. Unlike many model organisms, the dog genome has a comparatively poor functional annotation, which hampers gene discovery for development, morphology, disease, and behavior. To fill this gap, we established a comprehensive tissue biobank for both the dog and wolf samples. The biobank consists of 5485 samples representing 132 tissues from 13 dogs, 12 dog embryos, and 24 wolves. In a subset of 100 tissues from nine dogs and 12 embryos, we characterized gene expression activity for each promoter, including alternative and novel, i.e., previously not annotated, promoter regions, using the 5' targeting RNA sequencing technology STRT2-seq. We identified over 100,000 promoter region candidates in the recent canine genome assembly, CanFam4, including over 45,000 highly reproducible sites with gene expression and respective tissue enrichment levels. We provide a promoter and gene expression atlas with interactive, open data resources, including a data coordination center and genome browser track hubs. We demonstrated the applicability of Dog Genome Annotation (DoGA) data and resources using multiple examples spanning canine embryonic development, morphology and behavior, and diseases across species.

A simple genetic architecture underlies morphological variation in dogs.

Domestic dogs exhibit tremendous phenotypic diversity, including a greater variation in body size than any other terrestrial mammal. Here, we generate a high density map of canine genetic variation by genotyping 915 dogs from 80 domestic dog breeds, 83 wild canids, and 10 outbred African shelter dogs across 60,968 single-nucleotide polymorphisms (SNPs). Coupling this genomic resource with external measurements from breed standards and individuals as well as skeletal measurements from museum specimens, we identify 51 regions of the dog genome associated with phenotypic variation among breeds in 57 traits. The complex traits include average breed body size and external body dimensions and cranial, dental, and long bone shape and size with and without allometric scaling. In contrast to the results from association mapping of quantitative traits in humans and domesticated plants, we find that across dog breeds, a small number of quantitative trait loci (< or = 3) explain the majority of phenotypic variation for most of the traits we studied. In addition, many genomic regions show signatures of recent selection, with most of the highly differentiated regions being associated with breed-defining traits such as body size, coat characteristics, and ear floppiness. Our results demonstrate the efficacy of mapping multiple traits in the domestic dog using a database of genotyped individuals and highlight the important role human-directed selection has played in altering the genetic architecture of key traits in this important species.

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.

Identification of Novel Coloboma Candidate Genes through Conserved Gene Expression Analyses across Four Vertebrate Species.

Ocular coloboma (OC) is a failure of complete optic fissure closure during embryonic development and presents as a tissue defect along the proximal-distal axis of the ventral eye. It is classed as part of the clinical spectrum of structural eye malformations with microphthalmia and anophthalmia, collectively abbreviated to MAC. Despite deliberate attempts to identify causative variants in MAC, many patients remain without a genetic diagnosis. To reveal potential candidate genes, we utilised transcriptomes experimentally generated from embryonic eye tissues derived from humans, mice, zebrafish, and chicken at stages coincident with optic fissure closure. Our in-silico analyses found 10 genes with optic fissure-specific enriched expression: ALDH1A3, BMPR1B, EMX2, EPHB3, NID1, NTN1, PAX2, SMOC1, TENM3, and VAX1. In situ hybridization revealed that all 10 genes were broadly expressed ventrally in the developing eye but that only PAX2 and NTN1 were expressed in cells at the edges of the optic fissure margin. Of these conserved optic fissure genes, EMX2, NID1, and EPHB3 have not previously been associated with human MAC cases. Targeted genetic manipulation in zebrafish embryos using CRISPR/Cas9 caused the developmental MAC phenotype for emx2 and ephb3. We analysed available whole genome sequencing datasets from MAC patients and identified a range of variants with plausible causality. In combination, our data suggest that expression of genes involved in ventral eye development is conserved across a range of vertebrate species and that EMX2, NID1, and EPHB3 are candidate loci that warrant further functional analysis in the context of MAC and should be considered for sequencing in cohorts of patients with structural eye malformations.

Vessel and blood specification override cardiac potential in anterior mesoderm.

Organ progenitors arise within organ fields, embryonic territories that are larger than the regions required for organ formation. Little is known about the regulatory pathways that define organ field boundaries and thereby limit organ size. Here we identify a mechanism for restricting heart size through confinement of the developmental potential of the heart field. Via fate mapping in zebrafish, we locate cardiac progenitors within hand2-expressing mesoderm and demonstrate that hand2 potentiates cardiac differentiation within this region. Beyond the rostral boundary of hand2 expression, we find progenitors of vessel and blood lineages. In embryos deficient in vessel and blood specification, rostral mesoderm undergoes a fate transformation and generates ectopic cardiomyocytes. Therefore, induction of vessel and blood specification represses cardiac specification and delimits the capacity of the heart field. This regulatory relationship between cardiovascular pathways suggests strategies for directing progenitor cell differentiation to facilitate cardiac regeneration.

Complex population structure in African village dogs and its implications for inferring dog domestication history.

High genetic diversity of East Asian village dogs has recently been used to argue for an East Asian origin of the domestic dog. However, global village dog genetic diversity and the extent to which semiferal village dogs represent distinct, indigenous populations instead of admixtures of various dog breeds has not been quantified. Understanding these issues is critical to properly reconstructing the timing, number, and locations of dog domestication. To address these questions, we sampled 318 village dogs from 7 regions in Egypt, Uganda, and Namibia, measuring genetic diversity >680 bp of the mitochondrial D-loop, 300 SNPs, and 89 microsatellite markers. We also analyzed breed dogs, including putatively African breeds (Afghan hounds, Basenjis, Pharaoh hounds, Rhodesian ridgebacks, and Salukis), Puerto Rican street dogs, and mixed breed dogs from the United States. Village dogs from most African regions appear genetically distinct from non-native breed and mixed-breed dogs, although some individuals cluster genetically with Puerto Rican dogs or United States breed mixes instead of with neighboring village dogs. Thus, African village dogs are a mosaic of indigenous dogs descended from early migrants to Africa, and non-native, breed-admixed individuals. Among putatively African breeds, Pharaoh hounds, and Rhodesian ridgebacks clustered with non-native rather than indigenous African dogs, suggesting they have predominantly non-African origins. Surprisingly, we find similar mtDNA haplotype diversity in African and East Asian village dogs, potentially calling into question the hypothesis of an East Asian origin for dog domestication.

From head to hind: Elucidating function through contrasting morphometrics of ancient and modern pedigree dogs.

Used together, caliper- and geometric-based morphometric analyses provide complimentary approaches to classifying form and function of archaeozoological remains. Here we apply these analytical tools to the skeletal remains of an ancient male dog unearthed from a rural farm settlement of Roman date near present day Warmington, United Kingdom. Our comparisons of the Warmington Roman dog against the morphological characteristics of modern dog breeds enabled us to establish the former's size and shape. It was of medium stature. Analysis of viscerocrania and neurocrania indicate it falls within the meso- to dolichocephalic rankings of modern dogs. The neurocranium shape and the dimensions of its long bones strongly suggest that the Warmington dog shares similarities to modern sight hounds. Historically sight hounds were bred for speed, as necessitated of a hunter that runs down small prey. Our analysis suggests that the Warmington dog was likely bred for, or derived from, Roman hunting stock. By revealing the Warmington Roman dog's form from cranial and postcranial analyses, we shed light on Roman life in one of the furthest outposts of the Roman Empire.

Surveillance of a vomiting outbreak in dogs in the UK using owner-derived and internet search data.

BACKGROUND: In early 2020, the Small Animal Veterinary Surveillance Network reported evidence of an outbreak of acute prolific vomiting in dogs in the UK. The aims of this study were to investigate whether there was evidence for a vomiting outbreak in Dogslife and Google Trends data and to describe its characteristics. METHODS: Incidence of Dogslife vomiting reports and the Google search index for 'dog vomiting' and 'puppy vomiting' between December 2019 and March 2020 was compared to the respective data from the same months in previous years. Risks for dogs vomiting and factors influencing veterinary attendance in Dogslife were identified using multivariable logistic regression. RESULTS: This study confirmed a vomiting outbreak was evident in UK dogs between December 2019 and March 2020 using data from Dogslife and Google Trends. The odds of a vomiting incident being reported to Dogslife was 1.51 (95% CI: 1.24-1.84) in comparison to previous years. Dogslife data identified differences in owner-decision making when seeking veterinary attention and identified factors associated with dogs at higher odds of experiencing a vomiting episode. CONCLUSION: Owner-derived data including questionnaires and internet search queries should be considered a valid, valuable source of information for veterinary population health surveillance.

A Highly Conserved Shh Enhancer Coordinates Hypothalamic and Craniofacial Development.

Enhancers that are conserved deep in evolutionary time regulate characteristics held in common across taxonomic classes. Here, deletion of the highly conserved Shh enhancer SBE2 (Shh brain enhancer 2) in mouse markedly reduced Shh expression within the embryonic brain specifically in the rostral diencephalon; however, no abnormal anatomical phenotype was observed. Secondary enhancer activity was subsequently identified which likely mediates low levels of expression. In contrast, when crossing the SBE2 deletion with the Shh null allele, brain and craniofacial development were disrupted; thus, linking SBE2 regulated Shh expression to multiple defects and further enabling the study of the effects of differing levels of Shh on embryogenesis. Development of the hypothalamus, derived from the rostral diencephalon, was disrupted along both the anterior-posterior (AP) and the dorsal-ventral (DV) axes. Expression of DV patterning genes and subsequent neuronal population induction were particularly sensitive to Shh expression levels, demonstrating a novel morphogenic context for Shh. The role of SBE2, which is highlighted by DV gene expression, is to step-up expression of Shh above the minimal activity of the second enhancer, ensuring the necessary levels of Shh in a regional-specific manner. We also show that low Shh levels in the diencephalon disrupted neighbouring craniofacial development, including mediolateral patterning of the bones along the cranial floor and viscerocranium. Thus, SBE2 contributes to hypothalamic morphogenesis and ensures there is coordination with the formation of the adjacent midline cranial bones that subsequently protect the neural tissue.