The TNNT2:c.95-108G>A variant is common in Maine Coons and shows no association with hypertrophic cardiomyopathy.

Hypertrophic cardiomyopathy (HCM) is a common and potentially fatal heart disease in many cat breeds. An intronic variant in TNNT2, c.95-108G>A, was recently reported as the cause of HCM in the Maine Coon. The aim of this study was to determine this variant's allele frequency in different populations and its possible association with HCM. Based on 160 Maine Coon samples collected in Belgium, Italy, Sweden and the USA, the variant's allele frequency was estimated to be 0.32. Analysis of the 99 Lives feline whole genome sequencing database showed that the TNNT2 variant also occurs in other breeds, as well as mixed-breed cats. Comparison of 31 affected and 58 healthy cats did not reveal significantly increased odds for HCM in homozygotes. Based on the combined evidence and in agreement with the standards and guidelines for the interpretation of sequence variants, this variant is currently classified as a variant of unknown significance and should not be used for breeding decisions regarding HCM.

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.

Sequence dynamics of three influenza A virus strains grown in different MDCK cell lines, including those expressing different sialic acid receptors.

Viruses are often cultured in cell lines for research and vaccine development, and those often differ from the natural hosts or tissues. Cell lines can also differ in the presence of virus receptors, such as the sialic acid (Sia) receptors used by influenza A viruses (IAV), which can vary in linkage (α2,3- or α2,6-linkage) and form (N-glycolylneuraminic acid [Neu5Gc] or N-acetylneuraminic acid [Neu5Ac]). The selective pressures resulting from passaging viruses in cell types with host-specific variations in viral receptors are still only partially understood. IAV are commonly cultured in MDCK cells which are both derived from canine kidney tubule epithelium and inherently heterogeneous. MDCK cells naturally present Neu5Ac and α2,3-linked Sia forms. Here, we examine natural MDCK variant lineages, as well as engineered variants that synthesize Neu5Gc and/or α2,6-linkages. We determined how viral genetic variation occurred within human H3N2, H1N1 pandemic and canine H3N2 IAV populations when serially passaged in MDCK cell lines that vary in cell type (MDCK-Type I or MDCK-Type II clones) and in Sia display. Deep sequencing of viral genomes showed small numbers of consensus-level mutations, mostly within the hemagglutinin (HA) gene. Both human IAV showed variants in the HA stem and the HA receptor-binding site of populations passaged in cells displaying Neu5Gc. Canine H3N2 showed variants near the receptor-binding site when passaged in cells displaying Neu5Gc or α2,6-linkages. Viruses replicated to low titres in MDCK-Type II cells, suggesting that not all cell types in heterogeneous MDCK cell populations are equally permissive to infection.

Body size, inbreeding, and lifespan in domestic dogs.

Inbreeding poses a real or potential threat to nearly every species of conservation concern. Inbreeding leads to loss of diversity at the individual level, which can cause inbreeding depression, and at the population level, which can hinder ability to respond to a changing environment. In closed populations such as endangered species and ex situ breeding programs, some degree of inbreeding is inevitable. It is therefore vital to understand how different patterns of breeding and inbreeding can affect fitness in real animals. Domestic dogs provide an excellent model, showing dramatic variation in degree of inbreeding and in lifespan, an important aspect of fitness that is known to be impacted by inbreeding in other species. There is a strong negative correlation between body size and lifespan in dogs, but it is unknown whether the higher rate of aging in large dogs is due to body size per se or some other factor associated with large size. We used dense genome-wide SNP array data to calculate average inbreeding for over 100 dog breeds based on autozygous segment length and found that large breeds tend to have higher coefficients of inbreeding than small breeds. We then used data from the Veterinary medical Database and other published sources to estimate life expectancies for pure and mixed breed dogs. When controlling for size, variation in inbreeding was not associated with life expectancy across breeds. When comparing mixed versus purebred dogs, however, mixed breed dogs lived about 1.2 years longer on average than size-matched purebred dogs. Furthermore, individual pedigree coefficients of inbreeding and lifespans for over 9000 golden retrievers showed that inbreeding does negatively impact lifespan at the individual level. Registration data from the American Kennel Club suggest that the molecular inbreeding patterns observed in purebred dogs result from specific breeding practices and/or founder effects and not the current population size. Our results suggest that recent inbreeding, as reflected in variation within a breed, is more likely to affect fitness than historic inbreeding, as reflected in variation among breeds. Our results also indicate that occasional outcrosses, as in mixed breed dogs, can have a substantial positive effect on fitness.

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.

Genetic recombination is targeted towards gene promoter regions in dogs.

The identification of the H3K4 trimethylase, PRDM9, as the gene responsible for recombination hotspot localization has provided considerable insight into the mechanisms by which recombination is initiated in mammals. However, uniquely amongst mammals, canids appear to lack a functional version of PRDM9 and may therefore provide a model for understanding recombination that occurs in the absence of PRDM9, and thus how PRDM9 functions to shape the recombination landscape. We have constructed a fine-scale genetic map from patterns of linkage disequilibrium assessed using high-throughput sequence data from 51 free-ranging dogs, Canis lupus familiaris. While broad-scale properties of recombination appear similar to other mammalian species, our fine-scale estimates indicate that canine highly elevated recombination rates are observed in the vicinity of CpG rich regions including gene promoter regions, but show little association with H3K4 trimethylation marks identified in spermatocytes. By comparison to genomic data from the Andean fox, Lycalopex culpaeus, we show that biased gene conversion is a plausible mechanism by which the high CpG content of the dog genome could have occurred.

Classification of feline hypertrophic cardiomyopathy-associated gene variants according to the American College of Medical Genetics and Genomics guidelines.

Introduction: The correct labeling of a genetic variant as pathogenic is important as breeding decisions based on incorrect DNA tests can lead to the unwarranted exclusion of animals, potentially compromising the long-term health of a population. In human medicine, the American college of Medical Genetics (ACMG) guidelines provide a framework for variant classification. This study aims to apply these guidelines to six genetic variants associated with hypertrophic cardiomyopathy (HCM) in certain cat breeds and to propose a modified criterion for variant classification. Methods: Genetic samples were sourced from five cat breeds: Maine Coon, Sphynx, Ragdoll, Devon Rex, and British Short- and Longhair. Allele frequencies were determined, and in the subset with phenotypes available, odds ratios to determine the association with HCM were calculated. In silico evaluation followed with joint evidence and data from other publications assisting in the classification of each variant. Results: Two variants, MYBPC3:c.91G > C [A31P] and MYBPC3:c.2453C > T [R818W], were designated as pathogenic. One variant, MYH7:c.5647G > A [E1883K], was found likely pathogenic, while the remaining three were labeled as variants of unknown significance. Discussion: Routine genetic testing is advised solely for the MYBPC3:c.91G > C [A31P] in the Maine Coon and MYBPC3:c.2453C > T [R818W] in the Ragdoll breed. The human ACMG guidelines serve as a suitable foundational tool to ascertain which variants to include; however, refining them for application in veterinary medicine might be beneficial.

Common Orthopedic Traits and Screening for Breeding Programs.

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

Genomic and Transcriptomic Characterization of Atypical Recurrent Flank Alopecia in the Cesky Fousek.

Non-inflammatory alopecia is a frequent skin problem in dogs, causing damaged coat integrity and compromised appearance of affected individuals. In this study, we examined the Cesky Fousek breed, which displays atypical recurrent flank alopecia (aRFA) at a high frequency. This type of alopecia can be quite severe and is characterized by seasonal episodes of well demarcated alopecic areas without hyperpigmentation. The genetic component responsible for aRFA remains unknown. Thus, here we aimed to identify variants involved in aRFA using a combination of histological, genomic, and transcriptomic data. We showed that aRFA is histologically similar to recurrent flank alopecia, characterized by a lack of anagen hair follicles and the presence of severely shortened telogen or kenogen hair follicles. We performed a genome-wide association study (GWAS) using 216 dogs phenotyped for aRFA and identified associations on chromosomes 19, 8, 30, 36, and 21, highlighting 144 candidate genes, which suggests a polygenic basis for aRFA. By comparing the skin cell transcription pattern of six aRFA and five control dogs, we identified 236 strongly differentially expressed genes (DEGs). We showed that the GWAS genes associated with aRFA are often predicted to interact with DEGs, suggesting their joint contribution to the development of the disease. Together, these genes affect four major metabolic pathways connected to aRFA: collagen formation, muscle structure/contraction, lipid metabolism, and the immune system.

A genome-wide association study to investigate genetic loci associated with primary glaucoma in American Cocker Spaniels.

OBJECTIVE: To identify genetic associations with primary glaucoma (PG) in American Cocker Spaniels using a genome-wide association study (GWAS). ANIMALS: A nationwide ambidirectional case-control cohort study was performed in American Cocker Spaniels that had an ophthalmic examination performed by a veterinarian. Ninety-four dogs with PG (cases) and 111 dogs without glaucoma (controls) met phenotypic criteria and had a blood sample collected after receiving informed owner consent. PROCEDURES: Genomic DNA was extracted from whole blood samples and genotyped (CanineHD BeadChip, Illumina Inc). A case-control GWAS using a linear mixed model was performed, and 3 significance thresholds were calculated (1) using a Bonferroni correction on all single nucleotide polymorphisms (SNPs) included in the GWAS, (2) using a Bonferroni correction on only the unlinked SNPs from a pruned data set, and (3) using 10,000 random phenotype permutations. RESULTS: Following genotype data quality control, 89 cases and 93 controls were included in the GWAS. We identified an association on canine chromosome (CFA10); however, it did not reach statistical significance. Potential candidate genes within the surrounding linkage disequilibrium interval include coiled-coil domain containing 85A (CCDC85A) and extracellular growth factor containing fibulin extracellular matrix protein 1 (EFEMP1). CLINICAL RELEVANCE: Primary glaucoma in the American Cocker Spaniel is a complex heterogeneous disease that may be influenced by a locus on CFA10. The candidate genes CCDC85A and EFEMP1 within the identified linkage disequilibrium interval have been shown to be involved in human open-angle glaucoma.

Complex Feline Disease Mapping Using a Dense Genotyping Array.

The current feline genotyping array of 63 k single nucleotide polymorphisms has proven its utility for mapping within breeds, and its use has led to the identification of variants associated with Mendelian traits in purebred cats. However, compared to single gene disorders, association studies of complex diseases, especially with the inclusion of random bred cats with relatively low linkage disequilibrium, require a denser genotyping array and an increased sample size to provide statistically significant associations. Here, we undertook a multi-breed study of 1,122 cats, most of which were admitted and phenotyped for nine common complex feline diseases at the Cornell University Hospital for Animals. Using a proprietary 340 k single nucleotide polymorphism mapping array, we identified significant genome-wide associations with hyperthyroidism, diabetes mellitus, and eosinophilic keratoconjunctivitis. These results provide genomic locations for variant discovery and candidate gene screening for these important complex feline diseases, which are relevant not only to feline health, but also to the development of disease models for comparative studies.

Microevolution of canine influenza virus in shelters and its molecular epidemiology in the United States.

Canine influenza virus (CIV) emerged around 2000 when an equine influenza virus (EIV) was transmitted to dogs in Florida. After 2003, the canine virus was carried by infected greyhounds to various parts of the United States and then became established in several large animal shelters, where it has continued to circulate. To better understand the evolution of CIV since its emergence, and particularly its microevolution in spatially restricted populations, we examined multiple gene segments of CIV from dogs resident in two large animal shelters in New York City during the period 2006 to 2009. In particular, we focused on viruses circulating in the two shelters in 2008 and 2009, which we found shared a common ancestor. While viruses in each shelter were generally monophyletic, we observed some gene flow between them. These shelter sequences were compared to earlier CIV isolates. The shelter viruses differed in 1 to 6 amino acids in each gene segment compared to viruses isolated in Florida between 2003 and 2005 and in Colorado in 2006 and 2008. A comparison of the sequences of equine and canine viruses revealed amino acid replacements that distinguished the viruses from the two hosts, but no clear evidence of positive selection indicative of host adaptation was detected, suggesting that any host range adaptation in CIV occurred early in the emergence of this virus or even before it transferred to dogs.

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

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

A genome-wide association study of deafness in three canine breeds.

Congenital deafness in the domestic dog is usually related to the presence of white pigmentation, which is controlled primarily by the piebald locus on chromosome 20 and also by merle on chromosome 10. Pigment-associated deafness is also seen in other species, including cats, mice, sheep, alpacas, horses, cows, pigs, and humans, but the genetic factors determining why some piebald or merle dogs develop deafness while others do not have yet to be determined. Here we perform a genome-wide association study (GWAS) to identify regions of the canine genome significantly associated with deafness in three dog breeds carrying piebald: Dalmatian, Australian cattle dog, and English setter. We include bilaterally deaf, unilaterally deaf, and matched control dogs from the same litter, phenotyped using the brainstem auditory evoked response (BAER) hearing test. Principal component analysis showed that we have different distributions of cases and controls in genetically distinct Dalmatian populations, therefore GWAS was performed separately for North American and UK samples. We identified one genome-wide significant association and 14 suggestive (chromosome-wide) associations using the GWAS design of bilaterally deaf vs. control Australian cattle dogs. However, these associations were not located on the same chromosome as the piebald locus, indicating the complexity of the genetics underlying this disease in the domestic dog. Because of this apparent complex genetic architecture, larger sample sizes may be needed to detect the genetic loci modulating risk in piebald dogs.

Genetic mapping of distal femoral, stifle, and tibial radiographic morphology in dogs with cranial cruciate ligament disease.

Cranial cruciate ligament disease (CCLD) is a complex trait. Ten measurements were made on orthogonal distal pelvic limb radiographs of 161 pure and mixed breed dogs with, and 55 without, cranial cruciate partial or complete ligament rupture. Dogs with CCLD had significantly smaller infrapatellar fat pad width, higher average tibial plateau angle, and were heavier than control dogs. The first PC weightings captured the overall size of the dog's stifle and PC2 weightings reflected an increasing tibial plateau angle coupled with a smaller fat pad width. Of these dogs, 175 were genotyped, and 144,509 polymorphisms were used in a genome-wide association study with both a mixed linear and a multi-locus model. For both models, significant (pgenome <3.46×10-7 for the mixed and< 6.9x10-8 for the multilocus model) associations were found for PC1, tibial diaphyseal length and width, fat pad base length, and femoral and tibial condyle width at LCORL, a known body size-regulating locus. Other body size loci with significant associations were growth hormone 1 (GH1), which was associated with the length of the fat pad base and the width of the tibial diaphysis, and a region on CFAX near IRS4 and ACSL4 in the multilocus model. The tibial plateau angle was associated significantly with a locus on CFA10 in the linear mixed model with nearest candidate genes BET1 and MYH9 and on CFA08 near candidate genes WDHD1 and GCH1. MYH9 has a major role in osteoclastogenesis. Our study indicated that tibial plateau slope is associated with CCLD and a compressed infrapatellar fat pad, a surrogate for stifle osteoarthritis. Because of the association between tibial plateau slope and CCLD, and pending independent validation, these candidate genes for tibial plateau slope may be tested in breeds susceptible to CCLD before they develop disease or are bred.

Gene expression in hip soft tissues in incipient canine hip dysplasia and osteoarthritis.

Canine hip dysplasia and developmental dysplasia of the human hip share demographic, phenotypic, and clinical features including the predisposition to develop osteoarthritis in affected joints. To support the results of genetic mapping studies for CHD and its concomitant osteoarthritis with functional information, we performed RNA-seq on hip capsule and teres ligament of affected and unaffected dogs. RNA seq showed that expressed genes segregated according age, capsule or ligament, and hip phenotype. Expression of HHIP, DACT2, and WIF1 was significantly higher in capsule from control hips than dysplastic hips indicating a disruption of the hedgehog signaling pathway. Expression of SPON 1, a key component of the WNT pathway, was increased significantly in both dysplastic capsule and ligament while FBN2 and EMILIN3 were significantly increased in dysplastic capsule. Of genes associated with human hip osteoarthritis, expression of ACAN, IGF1, CILP2, COL11A1, COL8A1, and HAPLN was increased significantly in dysplastic capsule. The significant increase in expression of PLA2F, TNFRSF, TMEM, and IGFBP in dysplastic capsule indicated an injury response. Gene set enrichment analysis revealed that genes involved in extracellular matrix structure, epithelial to mesenchymal transition, myogenesis, growth factor signaling, cancer and immune pathways were enriched in dysplastic capsule. For teres ligament from dysplastic joints, genes in retinoic signaling pathways and those encoding extracellular matrix molecules, but not proteoglycans, were enriched. Hip tissues respond to abnormal mechanics early in dysplastic hip development and these pathways present targets for intervention in the early synovitis and capsulitis secondary to canine and human hip dysplasia. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 37:313-324, 2019.

Recombination in feline immunodeficiency virus from feral and companion domestic cats.

BACKGROUND: Recombination is a relatively common phenomenon in retroviruses. We investigated recombination in Feline Immunodeficiency Virus from naturally-infected New Zealand domestic cats (Felis catus) by sequencing regions of the gag, pol and env genes. RESULTS: The occurrence of intragenic recombination was highest in env, with evidence of recombination in 6.4% (n = 156) of all cats. A further recombinant was identified in each of the gag (n = 48) and pol (n = 91) genes. Comparisons of phylogenetic trees across genes identified cases of incongruence, indicating intergenic recombination. Three (7.7%, n = 39) of these incongruencies were found to be significantly different using the Shimodaira-Hasegawa test.Surprisingly, our phylogenies from the gag and pol genes showed that no New Zealand sequences group with reference subtype C sequences within intrasubtype pairwise distances. Indeed, we find one and two distinct unknown subtype groups in gag and pol, respectively. These observations cause us to speculate that these New Zealand FIV strains have undergone several recombination events between subtype A parent strains and undefined unknown subtype strains, similar to the evolutionary history hypothesised for HIV-1 "subtype E".Endpoint dilution sequencing was used to confirm the consensus sequences of the putative recombinants and unknown subtype groups, providing evidence for the authenticity of these sequences. Endpoint dilution sequencing also resulted in the identification of a dual infection event in the env gene. In addition, an intrahost recombination event between variants of the same subtype in the pol gene was established. This is the first known example of naturally-occurring recombination in a cat with infection of the parent strains. CONCLUSION: Evidence of intragenic recombination in the gag, pol and env regions, and complex intergenic recombination, of FIV from naturally-infected domestic cats in New Zealand was found. Strains of unknown subtype were identified in all three gene regions. These results have implications for the use of the current FIV vaccine in New Zealand.

Genetic mapping of principal components of canine pelvic morphology.

BACKGROUND: Concentrated breeding effort to produce various body structures and behaviors of dogs to suit human demand has inadvertently produced unwanted traits and diseases that accompany the morphological and behavioral phenotypes. We explored the relationship between pelvic conformation and canine hip dysplasia (HD) because purebred dogs which are predisposed, or not, to HD share common morphologic features, respectively. Thirteen unique bilateral anatomical features of the pelvis were measured on 392 dogs of 51 breeds and 95 mixed breed dogs. Principal components (PCs) were derived to describe pelvic morphology. Dogs were genotyped at ~183,000 single nucleotide polymorphisms and their hip conformation was measured by the Norberg angle and angle of inclination between the femoral neck and diaphysis. RESULTS: No associations reached genome wide significance for the Norberg angle when averaged over both hips. PC1 was negatively correlated with the Norberg angle (r = -0.31; P < 0.05) but not the angle of inclination (r = -0.08; P > 0.05). PC1, 2, 4, and 5 differed significantly between male and female dogs confirming pelvic sexual dimorphism. With sex as a covariate, the eigenvector contribution to PC1 reflected the overall size of the pelvis and was significantly associated with the IGF-1 locus, a known contributor to canine body size. PC3, which represented a tradeoff between ilial length and ischial length in which a longer ischium is associated with a shorter ilium, was significantly associated with a marker on canine chromosome 16:5181388 bp. The closest candidate gene is TPK1, a thiamine-dependent enzyme and part of the PKA complex. Associations with the remaining PCs did not reach genome wide significance. CONCLUSION: IGF-1 was associated with the overall size of the pelvis and sex is related to pelvic size. Ilial/ischial proportion is genetically controlled and the closest candidate gene is thiamine-dependent and affects birth weight and development of the nervous system. Dogs with larger pelves tend to have smaller NAs consistent with increased tendency toward HD in large breed dogs. Based on the current study, pelvic shape alone was not strongly associated with canine hip dysplasia.

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.