Elizabeth Corey

Genetic structure in village dogs reveals a Central Asian domestication origin

Significance Dogs were the first domesticated species, but the precise timing and location of domestication are hotly debated. Using genomic data from 5,392 dogs, including a global set of 549 village dogs, we find strong evidence that dogs were domesticated in Central Asia, perhaps near present-day Nepal and Mongolia. Dogs in nearby regions (e.g., East Asia, India, and Southwest Asia) contain high levels of genetic diversity due to their proximity to Central Asia and large population sizes. Indigenous dog populations in the Neotropics and South Pacific have been largely replaced by European dogs, whereas those in Africa show varying degrees of European vs. indigenous African ancestry. 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.

Quantitative trait loci for hip dysplasia in a crossbreed canine pedigree

Canine hip dysplasia is a common developmental inherited trait characterized by hip laxity, subluxation or incongruity of the femoral head and acetabulum in affected hips. The inheritance pattern is complex and the mutations contributing to trait expression are unknown. In the study reported here, 240 microsatellite markers distributed in 38 autosomes and the X chromosome were genotyped on 152 dogs from three generations of a crossbred pedigree based on trait-free Greyhound and dysplastic Labrador Retriever founders. Interval mapping was undertaken to map the QTL underlying the quantitative dysplastic traits of maximum passive hip laxity (the distraction index), the dorsolateral subluxation score, and the Norberg angle. Permutation testing was used to derive the chromosome-wide level of significance at p < 0.05 for each QTL. Chromosomes 4, 9, 10, 11 (p < 0.01), 16, 20, 22, 25, 29 (p < 0.01), 30, 35, and 37 harbor putative QTL for one or more traits. Successful detection of QTL was due to the crossbreed pedigree, multiple-trait measurements, control of environmental background, and marked advancement in canine mapping tools.

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.

Differential Genetic Regulation of Canine Hip Dysplasia and Osteoarthritis

Background Canine hip dysplasia (HD) is a common polygenic trait characterized by hip malformation that results in osteoarthritis (OA). The condition in dogs is very similar to developmental dysplasia of the human hip which also leads to OA. Methodology/Principal Findings A total of 721 dogs, including both an association and linkage population, were genotyped. The association population included 8 pure breeds (Labrador retriever, Greyhounds, German Shepherd, Newfoundland, Golden retriever, Rottweiler, Border Collie and Bernese Mountain Dog). The linkage population included Labrador retrievers, Greyhounds, and their crosses. Of these, 366 dogs were genotyped at ∼22,000 single nucleotide polymorphism (SNP) loci and a targeted screen across 8 chromosomes with ∼3,300 SNPs was performed on 551 dogs (196 dogs were common to both sets). A mixed linear model approach was used to perform an association study on this combined association and linkage population. The study identified 4 susceptibility SNPs associated with HD and 2 SNPs associated with hip OA. Conclusion/Significance The identified SNPs included those near known genes (PTPRD, PARD3B, and COL15A1) reported to be associated with, or expressed in, OA in humans. This suggested that the canine model could provide a unique opportunity to identify genes underlying natural HD and hip OA, which are common and debilitating conditions in both dogs and humans.

Evaluation of a fibrillin 2 gene haplotype associated with hip dysplasia and incipient osteoarthritis in dogs

OBJECTIVE To determine whether a mutation in the fibrillin 2 gene (FBN2) is associated with canine hip dysplasia (CHD) and osteoarthritis in dogs. ANIMALS 1,551 dogs. Procedures-Hip conformation was measured radiographically. The FBN2 was sequenced from genomic DNA of 21 Labrador Retrievers and 2 Greyhounds, and a haplotype in intron 30 of FBN2 was sequenced in 90 additional Labrador Retrievers and 143 dogs of 6 other breeds. Steady-state values of FBN2 mRNA and control genes were measured in hip joint tissues of fourteen 8-month-old Labrador Retriever-Greyhound crossbreeds. RESULTS The Labrador Retrievers homozygous for a 10-bp deletion haplotype in intron 30 of FBN2 had significantly worse CHD as measured via higher distraction index and extended-hip joint radiograph score and a lower Norberg angle and dorsolateral subluxation score. Among 143 dogs of 6 other breeds, those homozygous for the same deletion haplotype also had significantly worse radiographic CHD. Among the 14 crossbred dogs, as the dorsolateral subluxation score decreased, the capsular FBN2 mRNA increased significantly. Those dogs with incipient hip joint osteoarthritis had significantly increased capsular FBN2 mRNA, compared with those dogs without osteoarthritis. Dogs homozygous for the FBN2 deletion haplotype had significantly less FBN2 mRNA in their femoral head articular cartilage. CONCLUSIONS AND CLINICAL RELEVANCE The FBN2 deletion haplotype was associated with CHD. Capsular gene expression of FBN2 was confounded by incipient secondary osteoarthritis in dysplastic hip joints. Genes influencing complex traits in dogs can be identified by genome-wide screening, fine mapping, and candidate gene screening.

Genomic analysis establishes correlation between growth and laryngeal neuropathy in Thoroughbreds

BackgroundEquine recurrent laryngeal neuropathy (RLN) is a bilateral mononeuropathy with an unknown pathogenesis that significantly affects performance in Thoroughbreds. A genetic contribution to the pathogenesis of RLN is suggested by the higher prevalence of the condition in offspring of RLN-affected than unaffected stallions. To better understand RLN pathogenesis and its genetic basis, we performed a genome-wide association (GWAS) of 282 RLN-affected and 268 control Thoroughbreds.ResultsWe found a significant association of RLN with the LCORL/NCAPG locus on ECA3 previously shown to affect body size in horses. Using height at the withers of 505 of these horses, we confirmed the strong association of this locus with body size, and demonstrated a significant phenotypic and genetic correlation between height and RLN grade in this cohort. Secondary genetic associations for RLN on ECA18 and X did not correlate with withers height in our cohort, but did contain candidate genes likely influencing muscle physiology and growth: myostatin (MSTN) and integral membrane protein 2A (ITM2A).ConclusionsThis linkage between body size and RLN suggests that selective breeding to reduce RLN prevalence would likely reduce adult size in this population. However, our results do not preclude the possibility of modifier loci that attenuate RLN risk without reducing size or performance, or that the RLN risk allele is distinct but tightly linked to the body size locus on ECA3. This study is both the largest body size GWAS and the largest RLN GWAS within Thoroughbred horses to date, and suggests that improved understanding of the relationship between genetics, equine growth rate, and RLN prevalence may significantly advance our understanding and management of this disease.

A novel iterative mixed model to remap three complex orthopedic traits in dogs

Hip dysplasia (HD), elbow dysplasia (ED), and rupture of the cranial (anterior) cruciate ligament (RCCL) are the most common complex orthopedic traits of dogs and all result in debilitating osteoarthritis. We reanalyzed previously reported data: the Norberg angle (a quantitative measure of HD) in 921 dogs, ED in 113 cases and 633 controls, and RCCL in 271 cases and 399 controls and their genotypes at ~185,000 single nucleotide polymorphisms. A novel fixed and random model with a circulating probability unification (FarmCPU) function, with marker-based principal components and a kinship matrix to correct for population stratification, was used. A Bonferroni correction at p<0.01 resulted in a P< 6.96 ×10−8. Six loci were identified; three for HD and three for RCCL. An associated locus at CFA28:34,369,342 for HD was described previously in the same dogs using a conventional mixed model. No loci were identified for RCCL in the previous report but the two loci for ED in the previous report did not reach genome-wide significance using the FarmCPU model. These results were supported by simulation which demonstrated that the FarmCPU held no power advantage over the linear mixed model for the ED sample but provided additional power for the HD and RCCL samples. Candidate genes for HD and RCCL are discussed. When using FarmCPU software, we recommend a resampling test, that a positive control be used to determine the optimum pseudo quantitative trait nucleotide-based covariate structure of the model, and a negative control be used consisting of permutation testing and the identical resampling test as for the non-permuted phenotypes.

Reply to Wang et al.: Sequencing datasets do not refute Central Asian domestication origin of dogs

We welcome the additional data and analyses of Wang et al. (1), but believe there are some misunderstandings regarding the methods and findings of Shannon et al. (2). First, although we merged Nepal and Mongolia when plotting linkage disequilibrium (LD) decay in figure 5B of ref. 2 for legibility, we did not assume Nepal and Mongolia represented a single, interbreeding population, and indeed computed separate LD scores for each population (figure 5A of ref. 2), matching Wang et al.’s (1) observation of slightly lower LD in Nepal than Mongolia. Although Nepal (along with India) is commonly considered part of South Asia, Nepal borders Central Asia. Dog populations in two Central Asian countries, Mongolia and Afghanistan, both have lower LD than India. Nepal does not border Southeast Asia. Because we cannot, given the resolution of current sampling … [↵][1]1To whom correspondence should be addressed. Email: arb359{at}cornell.edu. [1]: #xref-corresp-1-1