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

Athletic prognosis following arthroscopic removal of osteochondral fragments of the proximal dorsal sagittal ridge of the fetlock in Thoroughbreds intended for flat racing

OBJECTIVE

To establish racing prognosis in Thoroughbred yearlings with proximal sagittal ridge osteochondral lesions and compare them to dorsoproximal and palmar/plantar first phalanx osteochondral lesions.

ANIMALS

A total of 47 horses had proximal sagittal ridge lesions, 34 had palmar/plantar first phalanx lesions, and 115 had dorsoproximal first phalanx lesions.

STUDY DESIGN

Retrospective case series.

METHODS

Medical records of Thoroughbred yearling racing prospects treated arthroscopically for fetlock osteochondral lesions were reviewed. Data were collected from a public database, including the ability to train and race, earnings, starts, wins, and placed races. Racing prognosis was analyzed and compared between three lesion locations.

RESULTS

Of the proximal sagittal ridge group, 76.6% raced, 65.71% of the palmar/plantar first phalanx group, and 74.58% of the dorsoproximal first phalanx group. Career length was similar for all lesions. A lower number of "starts" in the group "five years old and older" was found for the dorsoproximal first phalanx group compared to the cohort with proximal sagittal ridge lesions. Mares had more "earnings" at two years old but fewer "starts" at "five years old" compared to colts and geldings. Forelimb lesions were predictive for fewer total career starts.

CONCLUSION

Proximal sagittal ridge lesions had similar prognoses to dorsoproximal and palmar/plantar first phalanx lesions. In older horses, there was an increase in the number of starts for the proximal sagittal ridge group compared to the dorsoproximal first phalanx group. Mares and forelimb lesions were associated with decreased racing starts.

CLINICAL SIGNIFICANCE

These findings aid in prognostication for Thoroughbreds with osteochondral lesions removed arthroscopically in the fetlock joint.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Equine Genotyping Arrays

High-quality genomic tools have been integral in understanding genomic architecture and function in the modern-day horse. The equine genetics community has a long tradition of pooling resources to develop genomic tools. Since the equine genome was sequenced in 2006, several iterations of high throughput genotyping arrays have been developed and released, enabling rapid and cost-effective genotyping. This review highlights the design considerations of each iteration, focusing on data available during development and outlining considerations in selecting the genetic variants included on each array. Additionally, we outline recent applications of equine genotyping arrays as well as future prospects and applications.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Differences between horse selection based on two forms of osteochondrosis in fetlock

Horses lose potential opportunities because of health problems. Available breeding strategies are not effective enough, probably also because of the different definition used and its genetic usefulness. The aim of the study was to compare the genetic background estimated by the genome-wide association study (GWAS) for osteochondrosis using two different scaling osteochondrosis (OC)/healthy and osteochondrosis dissecans (OCD)/healthy systems for evaluating the disease status of investigated fetlock joints. Two hundred one Warmblood horses trained for performance tests (87 stallions and 114 mares) were phenotyped and genotyped. Four fetlock x-ray images per horse were collected using the RTG Girth HF 80 and Vet Scan ray 3600. The DNA of each horse was genotyped using the BeadChip 70K. To identify SNPs that significantly affect the probability of osteochondrosis, two different methods were applied: the Cochran-Armitage test based on an additive mode of inheritance and logistic regression. The genetic background for osteochondrosis, expressed in the number of SNPs found with significant associations with osteochondrosis, was higher by evaluation in the scale of OCD/healthy horses (16 SNPs on several chromosomes mainly on the ECA1 and ECA10) than OC/healthy (2 SNPs on the ECA15 and one SNP on the ECA10). Detailed definition of osteochondrosis is needed in breeding and in veterinary practice. The genetic background for osteochondrosis and osteochondrosis dissecans seems not the same. Suggestive SNPs could be the candidate markers for osteochondrosis but should be checked on a larger population before usage.

Developing a 670k genotyping array to tag ~2M SNPs across 24 horse breeds

Background

To date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array.

Results

Using whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation.

Conclusions

Here, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3943-8) contains supplementary material, which is available to authorized users.

Single nucleotide polymorphisms associated with osteochondrosis dissecans in Warmblood horses at different stages of training

The genetic background of osteochondrosis dissecans (OCD) has been studied for years, but the compatibility of the position of markers has not been reached between results, probably because of unknown additional effects that may influence the results, such as definition of the trait, gene–environmental interactions and the dynamics of trait development. The aim of the study was to identify single nucleotide polymorphisms (SNP) associated with the occurrence of OCD in Polish Warmblood sport breed horses in two different stages of training. Warmblood horses (87 stallions and 114 mares) were phenotyped and genotyped. Horses were X-rayed twice, at the beginning and at the end of the tests (100 days for stallions and 60 days for mares). Ten images per horse were collected using digital equipment for the fetlocks, stifles and hocks. The DNA was genotyped using the Illumina Neogen Equine Array. Statistical analysis included the Cochran–Armitage test and logistic regression assuming an additive model of inheritance. The Monte Carlo Markov Chain method was also applied to determine heritability coefficients. Nineteen and twenty SNP were identified that were significantly associated with OCD using logistic regression at the first and second stage of training, respectively. Four SNP were significant for both stages of training. The estimation of the heritability of a horse’s OCD status does not achieve the same level at different stages of training. The study on the genetic background of horse OCD should include as much detailed information on their training as possible.

Polymorphisms in ten candidate genes are associated with conformational and locomotive traits in Spanish Purebred horses

The Spanish Purebred horses, also known as Andalusian horses, compete to the highest standards in international dressage events. Gait and conformation could be used as early selection criteria to detect young horses with promising dressage ability. Although the genetic background of equine size variation has been recently uncovered, the genetic basis of horse conformational and locomotive traits is not known, hampered by the complex genetic architecture underlying quantitative traits and the lack of phenotypic data. The aim of this study was to validate the loci associated with size in 144 Spanish Purebred horses, and to seek novel associations between loci previously associated with the development of osteochondrosis (OC) lesions and 20 conformational and locomotive traits. Ten loci were associated with different conformational and locomotive traits (LCORL/NCAPG, HMGA2, USP31, MECR, COL24A1, MGP, FAM184B, PTH1R, KLF3 and SGK1), and the LCORL/NCAPG association with size in the Spanish Purebred horse was validated. Except for HMGA2, all polymorphisms seem to influence both the prevalence of OC lesions and morphological characters, supporting the link between conformation and OC. Also, the implication of most genes in either immune and inflammatory responses and cellular growth, or ossification processes, reinforces the role that these mechanisms have in the aetiology of OC, as well as their reflection on the general conformation of the individual. These polymorphisms could be used in marker-assisted selection (MAS) programmes to improve desirable conformational traits, but taking into account their possible detrimental effect on OC prevalence.

Heritability and prevalence of selected osteochondrosis lesions in yearling Thoroughbred horses

BACKGROUND

Osteochondrosis is considered multifactorial in origin, with factors such as nutrition, conformation, body size, trauma and genetics thought to contribute to its pathogenesis. Few studies have investigated the effects of genetic variability of osteochondrosis in Thoroughbreds.

OBJECTIVES

To describe the prevalence and genetic variability of a subset of osteochondrosis lesions in a group of Thoroughbred yearlings.

STUDY DESIGN

Retrospective cohort study.

METHODS

Radiographs of 1962 Thoroughbred yearlings were retrieved from clinical records obtained between 2005 and 2013. Pedigree information was obtained from the Australian Stud Book. Osteochondrosis lesions were documented in selected joints and estimates of heritability were obtained by fitting linear mixed models in ASREML software.

RESULTS

The overall prevalence of osteochondrosis was 23%. Osteochondrosis was identified in 10% of stifle joints, 6% of hock joints and 8% of fetlock joints. The heritability estimates ranged from 0 to 0.21. The largest estimates were 0.10, 0.14, 0.16 and 0.21 for lesions of the distal intermediate ridge of the tibia, dorso-proximal proximal phalanx (P1), any stifle osteochondrosis, and lesions of the lateral trochlear ridge of the distal femur, respectively. Although calculated heritability estimates had high standard errors, meta-analyses combining the present results with published estimates were significant at 0.10, 0.17, 0.15 and 0.20 for stifle, tarsal, fetlock and these joints combined, respectively. In addition, there was a permanent environment attributable to the dam effect.

MAIN LIMITATIONS

Inclusion criteria were based on radiographic findings in specific joints at a specific age range in Thoroughbreds.

CONCLUSIONS

The present results indicate that only a proportion of osteochondrosis in Thoroughbreds is heritable. The permanent environment effects of the dam were observed to have effects on some categories of osteochondrosis.

Identification and validation of risk loci for osteochondrosis in standardbreds

BACKGROUND

Osteochondrosis (OC), simply defined as a failure of endochondral ossification, is a complex disease with both genetic and environmental risk factors that is commonly diagnosed in young horses, as well as other domestic species. Although up to 50 % of the risk for developing OC is reportedly inherited, specific genes and alleles underlying risk are thus far completely unknown. Regions of the genome identified as associated with OC vary across studies in different populations of horses. In this study, we used a cohort of Standardbred horses from the U.S. (n = 182) specifically selected for a shared early environment (to reduce confounding factors) to identify regions of the genome associated with tarsal OC. Subsequently, putative risk variants within these regions were evaluated in both the discovery population and an independently sampled validation population of Norwegian Standardbreds (n = 139) with tarsal OC.

RESULTS

After genome-wide association analysis of imputed data with information from >200,000 single nucleotide polymorphisms, two regions on equine chromosome 14 were associated with OC in the discovery cohort. Variant discovery in these and 30 additional regions of interest (including 11 from other published studies) was performed via whole-genome sequencing. 240 putative risk variants from 10 chromosomes were subsequently genotyped in both the discovery and validation cohorts. After correction for population structure, gait (trot or pace) and sex, the variants most highly associated with OC status in both populations were located within the chromosome 14 regions of association.

CONCLUSIONS

The association of putative risk alleles from within the same regions with disease status in two independent populations of Standardbreds suggest that these are true risk loci in this breed, although population-specific risk factors may still exist. Evaluation of these loci in other populations will help determine if they are specific to the Standardbred breed, or to tarsal OC or are universal risk loci for OC. Further work is needed to identify the specific variants underlying OC risk within these loci. This is the first step towards the long-term goal of constructing a genetic risk model for OC that allows for genetic testing and quantification of risk in individuals.

Genome-wide association study for jumping performances in French sport horses

A genome-wide association study was performed to identify single nucleotide polymorphisms (SNPs) associated with jumping performances of warmbloods in France. The 999 horses included in the study for jumping performances were sport horses [mostly Selle Français (68%), Anglo-Arabians (13%) and horses from the other European studbooks]. Horses were genotyped using the Illumina EquineSNP50 BeadChip. Of the 54,602 SNPs available on this chip, 44,424 were retained after quality testing. Phenotypes were obtained by deregressing official breeding values for jumping competitions to use all available information, that is, the performances of each horse as well as those of its relatives. Two models were used to test the effects of the genotypes on deregressed phenotypes: a single-marker mixed model and a haplotype-based mixed model (significant: P < 1E-05; suggestive: P < 1E-04). Both models included a polygenic effect to take into account familial structures. For jumping performances, one suggestive quantitative trait locus (QTL) located on chromosome 1 (BIEC2_31196 and BIEC2_31198) was detected with both models. This QTL explains 0.7% of the phenotypic variance. RYR2, a gene encoding a major calcium channel in cardiac muscle in humans and mice, is located 0.55 Mb from this potential QTL.

Risk of false positive genetic associations in complex traits with underlying population structure: a case study

Genome-wide association (GWA) studies are widely used to investigate the genetic etiology of diseases in domestic animals. In the horse, GWA studies using 40-50,000 single nucleotide polymorphisms (SNPs) in sample sizes of 30-40 individuals, consisting of only 6-14 affected horses, have led to the discovery of genetic mutations for simple monogenic traits. Equine neuroaxonal dystrophy is a common inherited neurological disorder characterized by symmetric ataxia. A case-control GWA study was performed using genotypes from 42,819 SNP marker loci distributed across the genome in 99 clinically phenotyped Quarter horses (37 affected, 62 unaffected). A significant GWA was not achieved although a suggestive association was uncovered when only the most stringently phenotyped NAD-affected horses (n = 10) were included (chromosome 8:62130605 and 62134644 [log(1/P) = 5.56]). Candidate genes (PIK3C3, RIT2, and SYT4) within the associated region were excluded through sequencing, association testing of uncovered variants and quantitative RT-PCR. It was concluded that variants in PIK3C3, RIT2, and SYT4 are not responsible for equine neuroaxonal dystrophy. This study demonstrates the risk of false positive associations when performing GWA studies on complex traits and underlying population structure when using 40-50,000 SNP markers and small sample size.

Genome-wide analysis of DNA methylation patterns in horse

Background

DNA methylation is an epigenetic regulatory mechanism that plays an essential role in mediating biological processes and determining phenotypic plasticity in organisms. Although the horse reference genome and whole transcriptome data are publically available the global DNA methylation data are yet to be known.

Results

We report the first genome-wide DNA methylation characteristics data from skeletal muscle, heart, lung, and cerebrum tissues of thoroughbred (TH) and Jeju (JH) horses, an indigenous Korea breed, respectively by methyl-DNA immunoprecipitation sequencing. The analysis of the DNA methylation patterns indicated that the average methylation density was the lowest in the promoter region, while the density in the coding DNA sequence region was the highest. Among repeat elements, a relatively high density of methylation was observed in long interspersed nuclear elements compared to short interspersed nuclear elements or long terminal repeat elements. We also successfully identified differential methylated regions through a comparative analysis of corresponding tissues from TH and JH, indicating that the gene body regions showed a high methylation density.

Conclusions

We provide report the first DNA methylation landscape and differentially methylated genomic regions (DMRs) of thoroughbred and Jeju horses, providing comprehensive DMRs maps of the DNA methylome. These data are invaluable resource to better understanding of epigenetics in the horse providing information for the further biological function analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-598) contains supplementary material, which is available to authorized users.

Applied equine genetics

Genome sequencing of the domestic horse and subsequent advancements in the field of equine genomics have led to an explosion in the development of tools for mapping traits and diseases and evaluating gene expression. The objective of this review is to discuss the current progress in the field of equine genomics, with specific emphasis on assembly and analysis of the reference sequence and subsequent sequencing of a Quarter Horse mare; the genomic tools currently available to researchers and their implications in genomic investigations in the horse; the genomics of Mendelian and non-Mendelian traits; the genomics of performance traits and considerations regarding genetic testing in the horse. The whole-genome sequencing of a Quarter Horse mare has provided additional variants within the equine genome that extend past single nucleotide polymorphisms to include insertions/deletions and copy number variants. Equine single nucleotide polymorphism arrays have allowed for the investigation of both simple and complex genetic traits while DNA microarrays have provided a tool for examining gene expression across various tissues and with certain disease conditions. Recently, next-generation sequencing has become more affordable and both whole-genome DNA sequencing and transcriptome-wide RNA sequencing are methodologies that are being applied to equine genomic research. Research in the field of equine genomics continues to expand rapidly as the cost of genotyping and sequencing decreases, resulting in a need for quality bioinformatics software and expertise to appropriately handle both the size and complexity of these data.

Identification of genomic loci associated with Rhodococcus equi susceptibility in foals

Pneumonia caused by Rhodococcus equi is a common cause of disease and death in foals. Although agent and environmental factors contribute to the incidence of this disease, the genetic factors influencing the clinical outcomes of R. equi pneumonia are ill-defined. Here, we performed independent single nucleotide polymorphism (SNP)- and copy number variant (CNV)-based genome-wide association studies to identify genomic loci associated with R. equi pneumonia in foals. Foals at a large Quarter Horse breeding farm were categorized into 3 groups: 1) foals with R. equi pneumonia (clinical group [N = 43]); 2) foals with ultrasonographic evidence of pulmonary lesions that never developed clinical signs of pneumonia (subclinical group [N = 156]); and, 3) foals without clinical signs or ultrasonographic evidence of pneumonia (unaffected group [N = 49]). From each group, 24 foals were randomly selected and used for independent SNP- and CNV-based genome-wide association studies (GWAS). The SNP-based GWAS identified a region on chromosome 26 that had moderate evidence of association with R. equi pneumonia when comparing clinical and subclinical foals. A joint analysis including all study foals revealed a 3- to 4-fold increase in odds of disease for a homozygous SNP within the associated region when comparing the clinical group with either of the other 2 groups of foals or their combination. The region contains the transient receptor potential cation channel, subfamily M, member 2 (TRPM2) gene, which is involved in neutrophil function. No associations were identified in the CNV-based GWAS. Collectively, these data identify a region on chromosome 26 associated with R. equi pneumonia in foals, providing evidence that genetic factors may indeed contribute to this important disease of foals.

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.

Performance selection for Thoroughbreds racing in Hong Kong

REASONS FOR PERFORMING STUDY

Different indicators of racing performance are commonly used in the racing industry to assess the genetic superiority of racing Thoroughbreds. However, how well these indicators predict the performance of future progeny or siblings varies depending on the population and circumstances in which the indicators were recorded or achieved.

OBJECTIVES

To identify heritable indicators of racing performance for horses racing in Hong Kong.

STUDY DESIGN

Heritability analysis of racing performance traits.

METHODS

Performance data on the population of Thoroughbreds racing in Hong Kong between 3 September 2000 and 12 March 2011 (n = 4947) were acquired and used to estimate the heritabilities and probability values of fixed effects and covariates for a range of racing performance traits. Heritabilities for all performance traits were estimated using a single trait animal model. Each model included, as a minimum, the effects of sex, region of origin and trainer.

RESULTS

Heritability estimates for traits relating to finish position ranged from 0.01 to 0.06. Average handicap weight had a heritability of 0.07 ± 0.03. The effects of sex (fixed) and trainer (random) were significant (P<0.05) for all performance traits relating to earnings measures, handicap weights and finish positions. The heritability of win time at 1600 m was 0.52 ± 0.06 and was the only significant estimate of heritability for win time in the current study.

CONCLUSIONS

Although significantly affected by multiple environmental factors, certain indicators of Hong Kong racing performance can be reliably used to predict the performance of the individual's progeny or siblings. However, despite Hong Kong's controlled racing environment, these indicators appear to be no more heritable than in other less controlled racing environments.

Emerging genetic basis of osteochondritis dissecans

Genome-wide association studies (GWAS) provide an unbiased approach in the identification of genes that increase the risk for osteochondritis dissecans (OCD). Recent GWAS in humans, horses, and pigs are reviewed and genes identified. The identified genes tended to cluster with respect to function and biologic processes. GWAS in humans are a critical next step in the effort to provide a better understanding of the causes of OCD, which will, in turn, allow preventive strategies for treatment of adolescents and young adults who are at risk for the development of degenerative joint disease due to the effects of OCD.

Articular osteochondrosis: a comparison of naturally-occurring human and animal disease

BACKGROUND

Osteochondrosis (OC) is a common developmental orthopedic disease affecting both humans and animals. Despite increasing recognition of this disease among children and adolescents, its pathogenesis is incompletely understood because clinical signs are often not apparent until lesions have progressed to end-stage, and examination of cadaveric early lesions is not feasible. In contrast, both naturally-occurring and surgically-induced animal models of disease have been extensively studied, most notably in horses and swine, species in which OC is recognized to have profound health and economic implications. The potential for a translational model of human OC has not been recognized in the existing human literature.

OBJECTIVE

The purpose of this review is to highlight the similarities in signalment, predilection sites and clinical presentation of naturally-occurring OC in humans and animals and to propose a common pathogenesis for this condition across species.

STUDY DESIGN

Review.

METHODS

The published human and veterinary literature for the various manifestations of OC was reviewed. Peer-reviewed original scientific articles and species-specific review articles accessible in PubMed (US National Library of Medicine) were eligible for inclusion.

RESULTS

A broad range of similarities exists between OC affecting humans and animals, including predilection sites, clinical presentation, radiographic/MRI changes, and histological appearance of the end-stage lesion, suggesting a shared pathogenesis across species.

CONCLUSION

This proposed shared pathogenesis for OC between species implies that naturally-occurring and surgically-induced models of OC in animals may be useful in determining risk factors and for testing new diagnostic and therapeutic interventions that can be used in humans.

The genetics of equine osteochondrosis

Osteochondrosis (OC) develops in growing horses due to disturbed differentiation and maturation of cartilage, particularly at the predilection sites of the fetlock, hock and stifle joints. Horses with osteochondrotic lesions are at a high risk of developing orthopaedic problems later in life. This article briefly reviews the published heritability estimates for OC and offers perspectives for selection in the horse industry. Heritabilities for OC in Warmblood and Standardbred horses have been estimated at 0.1-0.4 in animal threshold models. Whole genome scans using microsatellites have identified 14 quantitative trait loci (QTL) and the eight most important QTL have been refined using dense marker maps. Genome-wide association studies with single nucleotide polymorphisms revealed further QTL in Thoroughbred, Standardbred and Hanoverian horses. Only a few QTL have corresponding locations among the different breeds. Comparative genomics using positional candidate genes and next-generation-sequencing may lead to new insights into the genetic determination of equine OC and might help in understanding the molecular mechanisms of its pathogenesis. Implementation of selection schemes based on breeding values, or even genomic selection against OC, should be considered as an option for improving equine musculoskeletal health.

Equine developmental orthopaedic diseases--a genome-wide association study of first phalanx plantar osteochondral fragments in Standardbred trotters

Palmar/plantar osteochondral fragments (POF) in fetlock joints commonly affect and influence the athletic performance of horses. In this study, we used the Equine SNP50 BeadChip® to perform a genome-wide association study of metatarsophalangeal POF in 176 Norwegian Standardbred trotter yearlings. Putative quantitative trait loci (QTL) for medial and/or lateral POF, and medial POF only were identified on ECA1, 2, 7, 9 and 31, whereas for lateral POF, only on ECA7, 11, 27 and X. The moderate number of QTL evidences a complex inheritance and suggests various genes controlling POF development in medial and lateral locations.

Genome-wide association study implicates testis-sperm specific FKBP6 as a susceptibility locus for impaired acrosome reaction in stallions

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals. Despite compelling evidence about the genetic control over acrosome biogenesis and function, the genomics of IAR is as yet poorly understood, providing no molecular tools for diagnostics. Here we conducted Equine SNP50 Beadchip genotyping and GWAS using 7 IAR–affected and 37 control Thoroughbred stallions. A significant (P<6.75E-08) genotype–phenotype association was found in horse chromosome 13 in FK506 binding protein 6 (FKBP6). The gene belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. Direct sequencing of FKBP6 exons in cases and controls identified SNPs g.11040315G>A and g.11040379C>A (p.166H>N) in exon 4 that were significantly associated with the IAR phenotype both in the GWAS cohort (n = 44) and in a large multi-breed cohort of 265 horses. All IAR stallions were homozygous for the A-alleles, while this genotype was found only in 2% of controls. The equine FKBP6 was exclusively expressed in testis and sperm and had 5 different transcripts, of which 4 were novel. The expression of this gene in AC/AG heterozygous controls was monoallelic, and we observed a tendency for FKBP6 up-regulation in IAR stallions compared to controls. Because exon 4 SNPs had no effect on the protein structure, it is likely that FKBP6 relates to the IAR phenotype via regulatory or modifying functions. In conclusion, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in mammals. FKBP6 genotyping is recommended for the detection of IAR–susceptible individuals among potential breeding stallions. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.

Impaired acrosomal reaction (IAR) of sperm causes male subfertility in humans and animals, and currently the molecular causes of the condition are not known. Here we report the mapping, identification, and functional analysis of a susceptibility locus for IAR in stallions. The candidate region was mapped to horse chromosome 13 by SNP genotyping and GWAS of 7 IAR affected and 44 control Thoroughbred stallions. Re-sequencing and case-control analysis of functionally relevant candidate genes in the region identified FKBP6 gene as a significantly associated locus. The association was confirmed by genotyping 265 male horses of multiple breeds. FKBP6 belongs to the immunophilins FKBP family known to be involved in meiosis, calcium homeostasis, clathrin-coated vesicles, and membrane fusions. We showed that the equine FKBP6 is exclusively and monoallelically expressed in testis and sperm and has 5 different transcripts, of which 4 were novel. Overall, FKBP6 was considered a susceptibility gene of incomplete penetrance for IAR in stallions and a candidate gene for male subfertility in other mammals. Successful use of sperm as a source of DNA and RNA propagates non-invasive sample procurement for fertility genomics in animals and humans.