Standards and guidelines for canine clinical genetic testing laboratories.

This publication represents a proposed approach to quality standards and guidelines for canine clinical genetic testing laboratories. Currently, there are no guidelines for laboratories performing clinical testing on dogs. Thus, there is no consensus set of protocols that set the minimal standards of quality among these laboratories, potentially causing variable results between laboratories, inconsistencies in reporting, and the inability to share information that could impact testing among organizations. A minimal standard for quality in testing is needed as breeders use the information from genetic testing to make breeding choices and irreversible decisions regarding spay, neuter or euthanasia. Incorrect results can have significant impact on the health of the dogs being tested and on their subsequent progeny. Because of the potentially serious consequences of an incorrect result or incorrect interpretation, results should be reviewed by and reported by individuals who meet a minimum standard of qualifications. Quality guidelines for canine genetic testing laboratories should include not only the analytical phase, but also the preanalytical and postanalytical phases, as this document attempts to address.

Heritability of Recurrent Exertional Rhabdomyolysis in Standardbred and Thoroughbred Racehorses Derived From SNP Genotyping Data.

Recurrent exertional rhabdomyolysis (RER) in Thoroughbred and Standardbred racehorses is characterized by episodes of muscle rigidity and cell damage that often recur upon strenuous exercise. The objective was to evaluate the importance of genetic factors in RER by obtaining an unbiased estimate of heritability in cohorts of unrelated Thoroughbred and Standardbred racehorses. Four hundred ninety-one Thoroughbred and 196 Standardbred racehorses were genotyped with the 54K or 74K SNP genotyping arrays. Heritability was calculated from genome-wide SNP data with a mixed linear and Bayesian model, utilizing the standard genetic relationship matrix (GRM). Both the mixed linear and Bayesian models estimated heritability of RER in Thoroughbreds to be approximately 0.34 and in Standardbred racehorses to be approximately 0.45 after adjusting for disease prevalence and sex. To account for potential differences in the genetic architecture of the underlying causal variants, heritability estimates were adjusted based on linkage disequilibrium weighted kinship matrix, minor allele frequency and variant effect size, yielding heritability estimates that ranged between 0.41-0.46 (Thoroughbreds) and 0.39-0.49 (Standardbreds). In conclusion, between 34-46% and 39-49% of the variance in RER susceptibility in Thoroughbred and Standardbred racehorses, respectively, can be explained by the SNPs present on these 2 genotyping arrays, indicating that RER is moderately heritable. These data provide further rationale for the investigation of genetic mutations associated with RER susceptibility.

A high density SNP array for the domestic horse and extant Perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies.

An equine SNP genotyping array was developed and evaluated on a panel of samples representing 14 domestic horse breeds and 18 evolutionarily related species. More than 54,000 polymorphic SNPs provided an average inter-SNP spacing of ∼43 kb. The mean minor allele frequency across domestic horse breeds was 0.23, and the number of polymorphic SNPs within breeds ranged from 43,287 to 52,085. Genome-wide linkage disequilibrium (LD) in most breeds declined rapidly over the first 50-100 kb and reached background levels within 1-2 Mb. The extent of LD and the level of inbreeding were highest in the Thoroughbred and lowest in the Mongolian and Quarter Horse. Multidimensional scaling (MDS) analyses demonstrated the tight grouping of individuals within most breeds, close proximity of related breeds, and less tight grouping in admixed breeds. The close relationship between the Przewalski's Horse and the domestic horse was demonstrated by pair-wise genetic distance and MDS. Genotyping of other Perissodactyla (zebras, asses, tapirs, and rhinoceros) was variably successful, with call rates and the number of polymorphic loci varying across taxa. Parsimony analysis placed the modern horse as sister taxa to Equus przewalski. The utility of the SNP array in genome-wide association was confirmed by mapping the known recessive chestnut coat color locus (MC1R) and defining a conserved haplotype of ∼750 kb across all breeds. These results demonstrate the high quality of this SNP genotyping resource, its usefulness in diverse genome analyses of the horse, and potential use in related species.

A genome-wide association study demonstrates significant genetic variation for fracture risk in Thoroughbred racehorses.

BACKGROUND: Thoroughbred racehorses are subject to non-traumatic distal limb bone fractures that occur during racing and exercise. Susceptibility to fracture may be due to underlying disturbances in bone metabolism which have a genetic cause. Fracture risk has been shown to be heritable in several species but this study is the first genetic analysis of fracture risk in the horse. RESULTS: Fracture cases (n = 269) were horses that sustained catastrophic distal limb fractures while racing on UK racecourses, necessitating euthanasia. Control horses (n = 253) were over 4 years of age, were racing during the same time period as the cases, and had no history of fracture at the time the study was carried out. The horses sampled were bred for both flat and National Hunt (NH) jump racing. 43,417 SNPs were employed to perform a genome-wide association analysis and to estimate the proportion of genetic variance attributable to the SNPs on each chromosome using restricted maximum likelihood (REML). Significant genetic variation associated with fracture risk was found on chromosomes 9, 18, 22 and 31. Three SNPs on chromosome 18 (62.05 Mb - 62.15 Mb) and one SNP on chromosome 1 (14.17 Mb) reached genome-wide significance (p < 0.05) in a genome-wide association study (GWAS). Two of the SNPs on ECA 18 were located in a haplotype block containing the gene zinc finger protein 804A (ZNF804A). One haplotype within this block has a protective effect (controls at 1.95 times less risk of fracture than cases, p = 1 × 10(-4)), while a second haplotype increases fracture risk (cases at 3.39 times higher risk of fracture than controls, p = 0.042). CONCLUSIONS: Fracture risk in the Thoroughbred horse is a complex condition with an underlying genetic basis. Multiple genomic regions contribute to susceptibility to fracture risk. This suggests there is the potential to develop SNP-based estimators for genetic risk of fracture in the Thoroughbred racehorse, using methods pioneered in livestock genetics such as genomic selection. This information would be useful to racehorse breeders and owners, enabling them to reduce the risk of injury in their horses.

Genomic analysis establishes correlation between growth and laryngeal neuropathy in Thoroughbreds.

BACKGROUND: Equine 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. RESULTS: We 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). CONCLUSIONS: This 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.

The utility of low-density genotyping for imputation in the Thoroughbred horse.

BACKGROUND: Despite the dramatic reduction in the cost of high-density genotyping that has occurred over the last decade, it remains one of the limiting factors for obtaining the large datasets required for genomic studies of disease in the horse. In this study, we investigated the potential for low-density genotyping and subsequent imputation to address this problem. RESULTS: Using the haplotype phasing and imputation program, BEAGLE, it is possible to impute genotypes from low- to high-density (50K) in the Thoroughbred horse with reasonable to high accuracy. Analysis of the sources of variation in imputation accuracy revealed dependence both on the minor allele frequency of the single nucleotide polymorphisms (SNPs) being imputed and on the underlying linkage disequilibrium structure. Whereas equidistant spacing of the SNPs on the low-density panel worked well, optimising SNP selection to increase their minor allele frequency was advantageous, even when the panel was subsequently used in a population of different geographical origin. Replacing base pair position with linkage disequilibrium map distance reduced the variation in imputation accuracy across SNPs. Whereas a 1K SNP panel was generally sufficient to ensure that more than 80% of genotypes were correctly imputed, other studies suggest that a 2K to 3K panel is more efficient to minimize the subsequent loss of accuracy in genomic prediction analyses. The relationship between accuracy and genotyping costs for the different low-density panels, suggests that a 2K SNP panel would represent good value for money. CONCLUSIONS: Low-density genotyping with a 2K SNP panel followed by imputation provides a compromise between cost and accuracy that could promote more widespread genotyping, and hence the use of genomic information in horses. In addition to offering a low cost alternative to high-density genotyping, imputation provides a means to combine datasets from different genotyping platforms, which is becoming necessary since researchers are starting to use the recently developed equine 70K SNP chip. However, more work is needed to evaluate the impact of between-breed differences on imputation accuracy.

Identification of a mutation associated with fatal Foal Immunodeficiency Syndrome in the Fell and Dales pony.

The Fell and Dales are rare native UK pony breeds at risk due to falling numbers, in-breeding, and inherited disease. Specifically, the lethal Mendelian recessive disease Foal Immunodeficiency Syndrome (FIS), which manifests as B-lymphocyte immunodeficiency and progressive anemia, is a substantial threat. A significant percentage (∼10%) of the Fell ponies born each year dies from FIS, compromising the long-term survival of this breed. Moreover, the likely spread of FIS into other breeds is of major concern. Indeed, FIS was identified in the Dales pony, a related breed, during the course of this work. Using a stepwise approach comprising linkage and homozygosity mapping followed by haplotype analysis, we mapped the mutation using 14 FIS-affected, 17 obligate carriers, and 10 adults of unknown carrier status to a ∼1 Mb region (29.8 - 30.8 Mb) on chromosome (ECA) 26. A subsequent genome-wide association study identified two SNPs on ECA26 that showed genome-wide significance after Bonferroni correction for multiple testing: BIEC2-692674 at 29.804 Mb and BIEC2-693138 at 32.19 Mb. The associated region spanned 2.6 Mb from ∼29.6 Mb to 32.2 Mb on ECA26. Re-sequencing of this region identified a mutation in the sodium/myo-inositol cotransporter gene (SLC5A3); this causes a P446L substitution in the protein. This gene plays a crucial role in the regulatory response to osmotic stress that is essential in many tissues including lymphoid tissues and during early embryonic development. We propose that the amino acid substitution we identify here alters the function of SLC5A3, leading to erythropoiesis failure and compromise of the immune system. FIS is of significant biological interest as it is unique and is caused by a gene not previously associated with a mammalian disease. Having identified the associated gene, we are now able to eradicate FIS from equine populations by informed selective breeding.

The same ELA class II risk factors confer equine insect bite hypersensitivity in two distinct populations.

Insect bite hypersensitivity (IBH) is a chronic allergic dermatitis common in horses. Affected horses mainly react against antigens present in the saliva from the biting midges, Culicoides ssp, and occasionally black flies, Simulium ssp. Because of this insect dependency, the disease is clearly seasonal and prevalence varies between geographical locations. For two distinct horse breeds, we genotyped four microsatellite markers positioned within the MHC class II region and sequenced the highly polymorphic exons two from DRA and DRB3, respectively. Initially, 94 IBH-affected and 93 unaffected Swedish born Icelandic horses were tested for genetic association. These horses had previously been genotyped on the Illumina Equine SNP50 BeadChip, which made it possible to ensure that our study did not suffer from the effects of stratification. The second population consisted of 106 unaffected and 80 IBH-affected Exmoor ponies. We show that variants in the MHC class II region are associated with disease susceptibility (p (raw) = 2.34 × 10(-5)), with the same allele (COR112:274) associated in two separate populations. In addition, we combined microsatellite and sequencing data in order to investigate the pattern of homozygosity and show that homozygosity across the entire MHC class II region is associated with a higher risk of developing IBH (p = 0.0013). To our knowledge this is the first time in any atopic dermatitis suffering species, including man, where the same risk allele has been identified in two distinct populations.

A genome-wide association study of osteochondritis dissecans in the Thoroughbred.

Osteochondrosis is a developmental orthopaedic disease that occurs in horses, other livestock species, companion animal species, and humans. The principal aim of this study was to identify quantitative trait loci (QTL) associated with osteochondritis dissecans (OCD) in the Thoroughbred using a genome-wide association study. A secondary objective was to test the effect of previously identified QTL in the current population. Over 300 horses, classified as cases or controls according to clinical findings, were genotyped for the Illumina Equine SNP50 BeadChip. An animal model was first implemented in order to adjust each horse's phenotypic status for average relatedness among horses and other potentially confounding factors which were present in the data. The genome-wide association test was then conducted on the residuals from the animal model. A single SNP on chromosome 3 was found to be associated with OCD at a genome-wide level of significance, as determined by permutation. According to the current sequence annotation, the SNP is located in an intergenic region of the genome. The effects of 24 SNPs, representing QTL previously identified in a sample of Hanoverian Warmblood horses, were tested directly in the animal model. When fitted alongside the significant SNP on ECA3, two of these SNPs were found to be associated with OCD. Confirmation of the putative QTL identified on ECA3 requires validation in an independent sample. The results of this study suggest that a significant challenge faced by equine researchers is the generation of sufficiently large data sets to effectively study complex diseases such as osteochondrosis.

Comparison of genomic and proteomic data in recurrent airway obstruction affected horses using Ingenuity Pathway Analysis®.

BACKGROUND: Recurrent airway obstruction (RAO) is a severe chronic respiratory disease affecting horses worldwide, though mostly in the Northern hemisphere. Environmental as well as genetic factors strongly influence the course and prognosis of the disease. Research has been focused on characterization of immunologic factors contributing to inflammatory responses, on genetic linkage analysis, and, more recently, on proteomic analysis of airway secretions from affected horses. The goal of this study was to investigate the interactions between eight candidate genes previously identified in a genetic linkage study and proteins expressed in bronchoalveolar lavage fluid (BALF) collected from healthy and RAO-affected horses. The analysis was carried out with Ingenuity Pathway Analysis® bioinformatics software. RESULTS: The gene with the greatest number of indirect interactions with the set of proteins identified is Interleukin 4 Receptor (IL-4R), whose protein has also been detected in BALF. Interleukin 21 receptor and chemokine (C-C motif) ligand 24 also showed a large number of interactions with the group of detected proteins. Protein products of other genes like that of SOCS5, revealed direct interactions with the IL-4R protein. The interacting proteins NOD2, RPS6KA5 and FOXP3 found in several pathways are reported regulators of the NFκB pathway. CONCLUSIONS: The pathways generated with IL-4R highlight possible important intracellular signaling cascades implicating, for instance, NFκB. Furthermore, the proposed interaction between SOCS5 and IL-4R could explain how different genes can lead to identical clinical RAO phenotypes, as observed in two Swiss Warmblood half sibling families because these proteins interact upstream of an important cascade where they may act as a functional unit.

A whole-genome scan for recurrent airway obstruction in Warmblood sport horses indicates two positional candidate regions.

Recurrent airway obstruction (RAO), or heaves, is a naturally occurring asthma-like disease that is related to sensitisation and exposure to mouldy hay and has a familial basis with a complex mode of inheritance. A genome-wide scanning approach using two half-sibling families was taken in order to locate the chromosome regions that contribute to the inherited component of this condition in these families. Initially, a panel of 250 microsatellite markers, which were chosen as a well-spaced, polymorphic selection covering the 31 equine autosomes, was used to genotype the two half-sibling families, which comprised in total 239 Warmblood horses. Subsequently, supplementary markers were added for a total of 315 genotyped markers. Each half-sibling family is focused around a severely RAO-affected stallion, and the phenotype of each individual was assessed for RAO and related signs, namely, breathing effort at rest, breathing effort at work, coughing, and nasal discharge, using an owner-based questionnaire. Analysis using a regression method for half-sibling family structures was performed using RAO and each of the composite clinical signs separately; two chromosome regions (on ECA13 and ECA15) showed a genome-wide significant association with RAO at P < 0.05. An additional 11 chromosome regions showed a more modest association. This is the first publication that describes the mapping of genetic loci involved in RAO. Several candidate genes are located in these regions, a number of which are interleukins. These are important signalling molecules that are intricately involved in the control of the immune response and are therefore good positional candidates.

Exclusion of linkage of the RYR1, CACNA1S, and ATP2A1 genes to recurrent exertional rhabdomyolysis in Thoroughbreds.

OBJECTIVE: To determine whether there was genetic linkage between the recurrent exertional rhabdomyolysis (RER) trait in Thoroughbred horse pedigrees and DNA markers in genes (the sarcoplasmic reticulum calcium release channel [RYR1] gene, the sarcoplasmic reticulum calcium ATPase [ATP2A1] gene, and the transverse tubule dihydropyridine receptor-voltage sensor [CACNA1S] gene) that are important in myoplasmic calcium regulation. ANIMALS: 34 horses in the University of Minnesota RER resource herd and 62 Thoroughbreds from 3 families of Thoroughbreds outside of the university in which RER-affected status was assigned after 2 or more episodes of ER had been observed. PROCEDURES: Microsatellite DNA markers from the RYR1, ATP2A1, and CACNA1S gene loci on equine chromosomes 10, 13, and 30 were identified. Genotypes were obtained for all horses in the 4 families affected by RER, and data were used to test for linkage of these 3 loci to the RER phenotype. RESULTS: Analysis of the RYR1, CACNA1S, and ATP2A1 microsatellites excluded a link between those markers and the RER trait. CONCLUSIONS AND CLINICAL RELEVANCE: It is likely that the heritable alterations in muscle contractility that are characteristic of RER are caused by a gene that is not yet known to cause related muscle disease in other species.

Using an Inbred Horse Breed in a High Density Genome-Wide Scan for Genetic Risk Factors of Insect Bite Hypersensitivity (IBH).

While susceptibility to hypersensitive reactions is a common problem amongst humans and animals alike, the population structure of certain animal species and breeds provides a more advantageous route to better understanding the biology underpinning these conditions. The current study uses Exmoor ponies, a highly inbred breed of horse known to frequently suffer from insect bite hypersensitivity, to identify genomic regions associated with a type I and type IV hypersensitive reaction. A total of 110 cases and 170 controls were genotyped on the 670K Axiom Equine Genotyping Array. Quality control resulted in 452,457 SNPs and 268 individuals being tested for association. Genome-wide association analyses were performed using the GenABEL package in R and resulted in the identification of two regions of interest on Chromosome 8. The first region contained the most significant SNP identified, which was located in an intron of the DCC netrin 1 receptor gene. The second region identified contained multiple top SNPs and encompassed the PIGN, KIAA1468, TNFRSF11A, ZCCHC2, and PHLPP1 genes. Although additional studies will be needed to validate the importance of these regions in horses and the relevance of these regions in other species, the knowledge gained from the current study has the potential to be a step forward in unraveling the complex nature of hypersensitive reactions.

Single linkage group per chromosome genetic linkage map for the horse, based on two three-generation, full-sibling, crossbred horse reference families.

A genetic linkage map of the horse consisting of 742 markers, which comprises a single linkage group for each of the autosomes and the X chromosome, is presented. The map has been generated from two three-generation full-sibling reference families, sired by the same stallion, in which there are 61 individuals in the F2 generation. Each linkage group has been assigned to a chromosome and oriented with reference to markers mapped by fluorescence in situ hybridization. The average interval between markers is 3.7 cM and the linkage groups collectively span 2772 cM. The 742 markers comprise 734 microsatellite and 8 gene-based markers. The utility of the microsatellite markers for comparative mapping has been significantly enhanced by comparing their flanking sequences with the human genome sequence; this enabled conserved segments between human and horse to be identified. The new map provides a valuable resource for genetically mapping traits of interest in the horse.

Assignment of the horse grey coat colour gene to ECA25 using whole genome scanning.

The dominant grey coat colour gene of horses has been mapped using a whole genome scanning approach. Samples from a large half-sibling pedigree of Thoroughbred horses were utilized in order to map the grey coat colour locus, G. Multiplex groups of microsatellite markers were developed and used to efficiently screen the horse genome at a resolution of approximately 22 cM, based on an estimated map length for the horse genome of 2720 cM. The grey gene was assigned to chromosome 25 (ECA25), one of the smaller acrocentric horse chromosomes. Based on the current state of knowledge of conserved synteny and coat colour genetics in other mammalian species, there are no obvious candidate genes for the grey gene in the region.