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

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

Genetic structure in village dogs reveals a Central Asian domestication origin.

Dogs were the first domesticated species, originating at least 15,000 y ago from Eurasian gray wolves. Dogs today consist primarily of two specialized groups--a diverse set of nearly 400 pure breeds and a far more populous group of free-ranging animals adapted to a human commensal lifestyle (village dogs). Village dogs are more genetically diverse and geographically widespread than purebred dogs making them vital for unraveling dog population history. Using a semicustom 185,805-marker genotyping array, we conducted a large-scale survey of autosomal, mitochondrial, and Y chromosome diversity in 4,676 purebred dogs from 161 breeds and 549 village dogs from 38 countries. Geographic structure shows both isolation and gene flow have shaped genetic diversity in village dog populations. Some populations (notably those in the Neotropics and the South Pacific) are almost completely derived from European stock, whereas others are clearly admixed between indigenous and European dogs. Importantly, many populations--including those of Vietnam, India, and Egypt-show minimal evidence of European admixture. These populations exhibit a clear gradient of short--range linkage disequilibrium consistent with a Central Asian domestication origin.

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.

Disease severity in a mouse model of ataxia telangiectasia is modulated by the DNA damage checkpoint gene Hus1.

The human genomic instability syndrome ataxia telangiectasia (A-T), caused by mutations in the gene encoding the DNA damage checkpoint kinase ATM, is characterized by multisystem defects including neurodegeneration, immunodeficiency and increased cancer predisposition. ATM is central to a pathway that responds to double-strand DNA breaks, whereas the related kinase ATR leads a parallel signaling cascade that is activated by replication stress. To dissect the physiological relationship between the ATM and ATR pathways, we generated mice defective for both. Because complete ATR pathway inactivation causes embryonic lethality, we weakened the ATR mechanism to different degrees by impairing HUS1, a member of the 911 complex that is required for efficient ATR signaling. Notably, simultaneous ATM and HUS1 defects caused synthetic lethality. Atm/Hus1 double-mutant embryos showed widespread apoptosis and died mid-gestationally. Despite the underlying DNA damage checkpoint defects, increased DNA damage signaling was observed, as evidenced by H2AX phosphorylation and p53 accumulation. A less severe Hus1 defect together with Atm loss resulted in partial embryonic lethality, with the surviving double-mutant mice showing synergistic increases in genomic instability and specific developmental defects, including dwarfism, craniofacial abnormalities and brachymesophalangy, phenotypes that are observed in several human genomic instability disorders. In addition to identifying tissue-specific consequences of checkpoint dysfunction, these data highlight a robust, cooperative configuration for the mammalian DNA damage response network and further suggest HUS1 and related genes in the ATR pathway as candidate modifiers of disease severity in A-T patients.

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.

The insulin-like growth factor 1 receptor (IGF1R) contributes to reduced size in dogs.

Domestic dog breeds have undergone intense selection for a variety of morphologic features, including size. Among small-dog breeds, defined as those averaging less than ~15 in. at the withers, there remains still considerable variation in body size. Yet essentially all such dogs are fixed for the same allele at the insulin-like growth factor 1 gene, which we and others previously found to be a size locus of large effect. In this study we sought to identify additional genes that contribute to tiny size in dogs using an association scan with the single nucleotide polymorphism (SNP) dataset CanMap, in which 915 purebred dogs were genotyped at 60,968 SNP markers. Our strongest association for tiny size (defined as breed-average height not more than 10 in. at the withers) was on canine chromosome 3 (p = 1.9 × 10(-70)). Fine mapping revealed a nonsynonymous SNP at chr3:44,706,389 that changes a highly conserved arginine at amino acid 204 to histidine in the insulin-like growth factor 1 receptor (IGF1R). This mutation is predicted to prevent formation of several hydrogen bonds within the cysteine-rich domain of the receptor's ligand-binding extracellular subunit. Nine of 13 tiny dog breeds carry the mutation and many dogs are homozygous for it. This work underscores the central importance of the IGF1 pathway in controlling the tremendous size diversity of dogs.

Genome sequence, comparative analysis and haplotype structure of the domestic dog.

Here we report a high-quality draft genome sequence of the domestic dog (Canis familiaris), together with a dense map of single nucleotide polymorphisms (SNPs) across breeds. The dog is of particular interest because it provides important evolutionary information and because existing breeds show great phenotypic diversity for morphological, physiological and behavioural traits. We use sequence comparison with the primate and rodent lineages to shed light on the structure and evolution of genomes and genes. Notably, the majority of the most highly conserved non-coding sequences in mammalian genomes are clustered near a small subset of genes with important roles in development. Analysis of SNPs reveals long-range haplotypes across the entire dog genome, and defines the nature of genetic diversity within and across breeds. The current SNP map now makes it possible for genome-wide association studies to identify genes responsible for diseases and traits, with important consequences for human and companion animal health.

The IGF1 small dog haplotype is derived from Middle Eastern grey wolves.

BACKGROUND: A selective sweep containing the insulin-like growth factor 1 (IGF1) gene is associated with size variation in domestic dogs. Intron 2 of IGF1 contains a SINE element and single nucleotide polymorphism (SNP) found in all small dog breeds that is almost entirely absent from large breeds. In this study, we surveyed a large sample of grey wolf populations to better understand the ancestral pattern of variation at IGF1 with a particular focus on the distribution of the small dog haplotype and its relationship to the origin of the dog. RESULTS: We present DNA sequence data that confirms the absence of the derived small SNP allele in the intron 2 region of IGF1 in a large sample of grey wolves and further establishes the absence of a small dog associated SINE element in all wild canids and most large dog breeds. Grey wolf haplotypes from the Middle East have higher nucleotide diversity suggesting an origin there. Additionally, PCA and phylogenetic analyses suggests a closer kinship of the small domestic dog IGF1 haplotype with those from Middle Eastern grey wolves. CONCLUSIONS: The absence of both the SINE element and SNP allele in grey wolves suggests that the mutation for small body size post-dates the domestication of dogs. However, because all small dogs possess these diagnostic mutations, the mutations likely arose early in the history of domestic dogs. Our results show that the small dog haplotype is closely related to those in Middle Eastern wolves and is consistent with an ancient origin of the small dog haplotype there. Thus, in concordance with past archeological studies, our molecular analysis is consistent with the early evolution of small size in dogs from the Middle East.See associated opinion by Driscoll and Macdonald: http://jbiol.com/content/9/2/10.

Linkage disequilibrium and demographic history of wild and domestic canids.

Assessing the extent of linkage disequilibrium (LD) in natural populations of a nonmodel species has been difficult due to the lack of available genomic markers. However, with advances in genotyping and genome sequencing, genomic characterization of natural populations has become feasible. Using sequence data and SNP genotypes, we measured LD and modeled the demographic history of wild canid populations and domestic dog breeds. In 11 gray wolf populations and one coyote population, we find that the extent of LD as measured by the distance at which r2=0.2 extends <10 kb in outbred populations to >1.7 Mb in populations that have experienced significant founder events and bottlenecks. This large range in the extent of LD parallels that observed in 18 dog breeds where the r2 value varies from approximately 20 kb to >5 Mb. Furthermore, in modeling demographic history under a composite-likelihood framework, we find that two of five wild canid populations exhibit evidence of a historical population contraction. Five domestic dog breeds display evidence for a minor population contraction during domestication and a more severe contraction during breed formation. Only a 5% reduction in nucleotide diversity was observed as a result of domestication, whereas the loss of nucleotide diversity with breed formation averaged 35%.

A mutation in the myostatin gene increases muscle mass and enhances racing performance in heterozygote dogs.

Double muscling is a trait previously described in several mammalian species including cattle and sheep and is caused by mutations in the myostatin (MSTN) gene (previously referred to as GDF8). Here we describe a new mutation in MSTN found in the whippet dog breed that results in a double-muscled phenotype known as the "bully" whippet. Individuals with this phenotype carry two copies of a two-base-pair deletion in the third exon of MSTN leading to a premature stop codon at amino acid 313. Individuals carrying only one copy of the mutation are, on average, more muscular than wild-type individuals (p = 7.43 x 10(-6); Kruskal-Wallis Test) and are significantly faster than individuals carrying the wild-type genotype in competitive racing events (Kendall's nonparametric measure, tau = 0.3619; p approximately 0.00028). These results highlight the utility of performance-enhancing polymorphisms, marking the first time a mutation in MSTN has been quantitatively linked to increased athletic performance.

Whole genome detection of sequence and structural polymorphism in six diverse horses.

The domesticated horse has played a unique role in human history, serving not just as a source of animal protein, but also as a catalyst for long-distance migration and military conquest. As a result, the horse developed unique physiological adaptations to meet the demands of both their climatic environment and their relationship with man. Completed in 2009, the first domesticated horse reference genome assembly (EquCab 2.0) produced most of the publicly available genetic variations annotations in this species. Yet, there are around 400 geographically and physiologically diverse breeds of horse. To enrich the current collection of genetic variants in the horse, we sequenced whole genomes from six horses of six different breeds: an American Miniature, a Percheron, an Arabian, a Mangalarga Marchador, a Native Mongolian Chakouyi, and a Tennessee Walking Horse, and mapped them to EquCab3.0 genome. Aside from extreme contrasts in body size, these breeds originate from diverse global locations and each possess unique adaptive physiology. A total of 1.3 billion reads were generated for the six horses with coverage between 15x to 24x per horse. After applying rigorous filtration, we identified and functionally annotated 17,514,723 Single Nucleotide Polymorphisms (SNPs), and 1,923,693 Insertions/Deletions (INDELs), as well as an average of 1,540 Copy Number Variations (CNVs) and 3,321 Structural Variations (SVs) per horse. Our results revealed putative functional variants including genes associated with size variation like LCORL gene (found in all horses), ZFAT in the Arabian, American Miniature and Percheron horses and ANKRD1 in the Native Mongolian Chakouyi horse. We detected a copy number variation in the Latherin gene that may be the result of evolutionary selection impacting thermoregulation by sweating, an important component of athleticism and heat tolerance. The newly discovered variants were formatted into user-friendly browser tracks and will provide a foundational database for future studies of the genetic underpinnings of diverse phenotypes within the horse.

Genetic architecture of the dog: sexual size dimorphism and functional morphology.

Purebred dogs are a valuable resource for genetic analysis of quantitative traits. Quantitative traits are complex, controlled by many genes that are contained within regions of the genome known as quantitative trait loci (QTL). The genetic architecture of quantitative traits is defined by the characteristics of these genes: their number, the magnitude of their effects, their positions in the genome and their interactions with each other. QTL analysis is a valuable tool for exploring genetic architecture, and highlighting regions of the genome that contribute to the variation of a trait within a population.

Mutation responsible for congenital photosensitivity and hyperbilirubinemia in Southdown sheep.

OBJECTIVE To identify the genetic cause for congenital photosensitivity and hyperbilirubinemia (CPH) in Southdown sheep. ANIMALS 73 Southdown sheep from a CPH research flock and 48 sheep of various breeds from commercial flocks without CPH. PROCEDURES Whole-genome sequencing was performed for a phenotypically normal Southdown sheep heterozygous for CPH. Heterozygous variants within Slco1b3 coding exons were identified, and exons that contained candidate mutations were amplified by PCR assay methods for Sanger sequencing. Blood samples from the other 72 Southdown sheep of the CPH research flock were used to determine plasma direct and indirect bilirubin concentrations. Southdown sheep with a plasma total bilirubin concentration < 0.3 mg/dL were classified as controls, and those with a total bilirubin concentration ≥ 0.3 mg/dL and signs of photosensitivity were classified as mutants. Sanger sequencing was used to determine the Slco1b3 genotype for all sheep. Genotypes were compared between mutants and controls of the CPH research flock and among all sheep. Protein homology was measured across 8 species to detect evolutionary conservation of Slco1b. RESULTS A nonsynonymous mutation at ovine Chr3:193,691,195, which generated a glycine-to-arginine amino acid change within the predicted Slco1b3 protein, was significantly associated with hyperbilirubinemia and predicted to be deleterious. That amino acid was conserved across 7 other mammalian species. CONCLUSIONS AND CLINICAL RELEVANCE Results suggested a nonsynonymous mutation in Slco1b3 causes CPH in Southdown sheep. This disease appears to be similar to Rotor syndrome in humans. Sheep with CPH might be useful animals for Rotor syndrome research.

Complex disease and phenotype mapping in the domestic dog.

The domestic dog is becoming an increasingly valuable model species in medical genetics, showing particular promise to advance our understanding of cancer and orthopaedic disease. Here we undertake the largest canine genome-wide association study to date, with a panel of over 4,200 dogs genotyped at 180,000 markers, to accelerate mapping efforts. For complex diseases, we identify loci significantly associated with hip dysplasia, elbow dysplasia, idiopathic epilepsy, lymphoma, mast cell tumour and granulomatous colitis; for morphological traits, we report three novel quantitative trait loci that influence body size and one that influences fur length and shedding. Using simulation studies, we show that modestly larger sample sizes and denser marker sets will be sufficient to identify most moderate- to large-effect complex disease loci. This proposed design will enable efficient mapping of canine complex diseases, most of which have human homologues, using far fewer samples than required in human studies.

Four loci explain 83% of size variation in the horse.

Horse body size varies greatly due to intense selection within each breed. American Miniatures are less than one meter tall at the withers while Shires and Percherons can exceed two meters. The genetic basis for this variation is not known. We hypothesize that the breed population structure of the horse should simplify efforts to identify genes controlling size. In support of this, here we show with genome-wide association scans (GWAS) that genetic variation at just four loci can explain the great majority of horse size variation. Unlike humans, which are naturally reproducing and possess many genetic variants with weak effects on size, we show that horses, like other domestic mammals, carry just a small number of size loci with alleles of large effect. Furthermore, three of our horse size loci contain the LCORL, HMGA2 and ZFAT genes that have previously been found to control human height. The LCORL/NCAPG locus is also implicated in cattle growth and HMGA2 is associated with dog size. Extreme size diversification is a hallmark of domestication. Our results in the horse, complemented by the prior work in cattle and dog, serve to pinpoint those very few genes that have played major roles in the rapid evolution of size during domestication.

DOG-SPOT database for comprehensive management of dog genetic research data.

Research laboratories studying the genetics of companion animals have no database tools specifically designed to aid in the management of the many kinds of data that are generated, stored and analyzed. We have developed a relational database, "DOG-SPOT," to provide such a tool. Implemented in MS-Access, the database is easy to extend or customize to suit a lab's particular needs. With DOG-SPOT a lab can manage data relating to dogs, breeds, samples, biomaterials, phenotypes, owners, communications, amplicons, sequences, markers, genotypes and personnel. Such an integrated data structure helps ensure high quality data entry and makes it easy to track physical stocks of biomaterials and oligonucleotides.

An expressed fgf4 retrogene is associated with breed-defining chondrodysplasia in domestic dogs.

Retrotransposition of processed mRNAs is a common source of novel sequence acquired during the evolution of genomes. Although the vast majority of retroposed gene copies, or retrogenes, rapidly accumulate debilitating mutations that disrupt the reading frame, a small percentage become new genes that encode functional proteins. By using a multibreed association analysis in the domestic dog, we demonstrate that expression of a recently acquired retrogene encoding fibroblast growth factor 4 (fgf4) is strongly associated with chondrodysplasia, a short-legged phenotype that defines at least 19 dog breeds including dachshund, corgi, and basset hound. These results illustrate the important role of a single evolutionary event in constraining and directing phenotypic diversity in the domestic dog.

Morphometrics within dog breeds are highly reproducible and dispute Rensch's rule.

Using 27 body measurements, we have identified 13 breed-defining metrics for 109 of 159 domestic dog breeds, most of which are recognized by the American Kennel Club (AKC). The data set included 1,155 dogs at least 1 year old (average 5.4 years), and for 53 breed populations, complete measurement data were collected from at least three males and three females. We demonstrate, first, that AKC breed standards are rigorously adhered to for most domestic breeds with little variation observed within breeds. Second, Rensch's rule, which describes a scaling among taxa such that sexual dimorphism is greater among larger species if males are the larger sex, with less pronounced differences in male versus female body size in smaller species, is not maintained in domestic dog breeds because the proportional size difference between males and females of small and large breeds is essentially the same. Finally, principal components (PCs) analysis describes both the overall body size (PC1) and the shape (length versus width) of the skeleton (PC2). That the integrity of the data set is sufficiently rich to discern PCs has strong implications for mapping studies, suggesting that individual measurements may not be needed for genetic studies of morphologic traits, particularly in the case of breed-defining traits that are typically under strong selection. Rather, phenotypes derived from data sets such as these, collected at a fraction of the effort and cost, may be used to direct whole-genome association studies aimed at understanding the genetic basis of fixed morphologic phenotypes defining distinct dog breeds.

A single IGF1 allele is a major determinant of small size in dogs.

The domestic dog exhibits greater diversity in body size than any other terrestrial vertebrate. We used a strategy that exploits the breed structure of dogs to investigate the genetic basis of size. First, through a genome-wide scan, we identified a major quantitative trait locus (QTL) on chromosome 15 influencing size variation within a single breed. Second, we examined genetic variation in the 15-megabase interval surrounding the QTL in small and giant breeds and found marked evidence for a selective sweep spanning a single gene (IGF1), encoding insulin-like growth factor 1. A single IGF1 single-nucleotide polymorphism haplotype is common to all small breeds and nearly absent from giant breeds, suggesting that the same causal sequence variant is a major contributor to body size in all small dogs.