Variants of the MSTN gene encoding myostatin are associated with muscle hypertrophy phenotypes in a range of mammalian species, most notably cattle, dogs, mice, and humans. Using a sample of registered Thoroughbred horses (n = 148), we have identified a novel MSTN sequence polymorphism that is strongly associated (g.66493737C>T, P = 4.85×10−8) with best race distance among elite racehorses (n = 79). This observation was independently validated (P = 1.91×10−6) in a resampled group of Thoroughbreds (n = 62) and in a cohort of Thoroughbreds (n = 37, P = 0.0047) produced by the same trainer. We observed that C/C horses are suited to fast, short-distance races; C/T horses compete favorably in middle-distance races; and T/T horses have greater stamina. Evaluation of retrospective racecourse performance (n = 142) and stallion progeny performance predict that C/C and C/T horses are more likely to be successful two-year-old racehorses than T/T animals. Here we describe for the first time the identification of a gene variant in Thoroughbred racehorses that is predictive of genetic potential for an athletic phenotype.

BackgroundThoroughbred horses have been selected for traits contributing to speed and stamina for centuries. It is widely recognized that inherited variation in physical and physiological characteristics is responsible for variation in individual aptitude for race distance, and that muscle phenotypes in particular are important.ResultsA genome-wide SNP-association study for optimum racing distance was performed using the EquineSNP50 Bead Chip genotyping array in a cohort of n = 118 elite Thoroughbred racehorses divergent for race distance aptitude. In a cohort-based association test we evaluated genotypic variation at 40,977 SNPs between horses suited to short distance (≤ 8 f) and middle-long distance (> 8 f) races. The most significant SNP was located on chromosome 18: BIEC2-417495 ~690 kb from the gene encoding myostatin (MSTN) [Punadj. = 6.96 × 10-6]. Considering best race distance as a quantitative phenotype, a peak of association on chromosome 18 (chr18:65809482-67545806) comprising eight SNPs encompassing a 1.7 Mb region was observed. Again, similar to the cohort-based analysis, the most significant SNP was BIEC2-417495 (Punadj. = 1.61 × 10-9; PBonf. = 6.58 × 10-5). In a candidate gene study we have previously reported a SNP (g.66493737C>T) in MSTN associated with best race distance in Thoroughbreds; however, its functional and genome-wide relevance were uncertain. Additional re-sequencing in the flanking regions of the MSTN gene revealed four novel 3' UTR SNPs and a 227 bp SINE insertion polymorphism in the 5' UTR promoter sequence. Linkage disequilibrium was highest between g.66493737C>T and BIEC2-417495 (r2 = 0.86).ConclusionsComparative association tests consistently demonstrated the g.66493737C>T SNP as the superior variant in the prediction of distance aptitude in racehorses (g.66493737C>T, P = 1.02 × 10-10; BIEC2-417495, Punadj. = 1.61 × 10-9). Functional investigations will be required to determine whether this polymorphism affects putative transcription-factor binding and gives rise to variation in gene and protein expression. Nonetheless, this study demonstrates that the g.66493737C>T SNP provides the most powerful genetic marker for prediction of race distance aptitude in Thoroughbreds.

Thoroughbred horses have been selected for exceptional racing performance resulting in system-wide structural and functional adaptations contributing to elite athletic phenotypes. Because selection has been recent and intense in a closed population that stems from a small number of founder animals Thoroughbreds represent a unique population within which to identify genomic contributions to exercise-related traits. Employing a population genetics-based hitchhiking mapping approach we performed a genome scan using 394 autosomal and X chromosome microsatellite loci and identified positively selected loci in the extreme tail-ends of the empirical distributions for (1) deviations from expected heterozygosity (Ewens-Watterson test) in Thoroughbred (n = 112) and (2) global differentiation among four geographically diverse horse populations (FST). We found positively selected genomic regions in Thoroughbred enriched for phosphoinositide-mediated signalling (3.2-fold enrichment; P<0.01), insulin receptor signalling (5.0-fold enrichment; P<0.01) and lipid transport (2.2-fold enrichment; P<0.05) genes. We found a significant overrepresentation of sarcoglycan complex (11.1-fold enrichment; P<0.05) and focal adhesion pathway (1.9-fold enrichment; P<0.01) genes highlighting the role for muscle strength and integrity in the Thoroughbred athletic phenotype. We report for the first time candidate athletic-performance genes within regions targeted by selection in Thoroughbred horses that are principally responsible for fatty acid oxidation, increased insulin sensitivity and muscle strength: ACSS1 (acyl-CoA synthetase short-chain family member 1), ACTA1 (actin, alpha 1, skeletal muscle), ACTN2 (actinin, alpha 2), ADHFE1 (alcohol dehydrogenase, iron containing, 1), MTFR1 (mitochondrial fission regulator 1), PDK4 (pyruvate dehydrogenase kinase, isozyme 4) and TNC (tenascin C). Understanding the genetic basis for exercise adaptation will be crucial for the identification of genes within the complex molecular networks underlying obesity and its consequential pathologies, such as type 2 diabetes. Therefore, we propose Thoroughbred as a novel in vivo large animal model for understanding molecular protection against metabolic disease.

SummaryAthletic performance is influenced by a complex interplay among the environment and a suite of genes, which contributes to system-wide structure and function. In a panel of elite and non-elite Thoroughbred horses (n = 148), we genotyped 68 SNPs in 17 putative exercise-relevant genes chosen from a genome scan for selection. We performed a series of case-control and quantitative association tests for relationships with racecourse performance. Thirteen SNPs in nine genes were significantly (P < 0.05) associated with a performance phenotype. We selected five SNPs in four genes (ACSS1, ACN9, COX4I1, PDK4) for validation in an independent sample set of elite and non-elite Thoroughbreds (n = 130). Two SNPs in the PDK4 gene were validated (P < 0.01) for associations with elite racing performance. When all samples were considered together (n = 278), the PDK4_ 38973231 SNP was strongly associated (P < 0.0005) with elite racing performance. Individuals with the A:A and A:G genotypes had a 16.2-16.6 lb advantage over G:G individuals in terms of handicap rating. Re-sequencing of the PDK4 gene and further genotyping will be required to identify the causative variant that is likely influencing exercise-induced variation in expression of the gene. Notwithstanding, this information may be employed as a marker for the selection of racehorses with the genetic potential for superior racing ability.

No abstract

SummaryReasons for performing study: The wild progenitors of the domestic horse were subject to natural selection for speed and stamina for millennia. Uniquely, this process has been augmented in Thoroughbreds, which have undergone at least 3 centuries of intense artificial selection for athletic phenotypes. While the phenotypic adaptations to exercise are well described, only a small number of the underlying genetic variants contributing to these phenotypes have been reported. Objectives: A panel of candidate performance-related genes was examined for DNA sequence variation in Thoroughbreds and the association with racecourse performance investigated. Materials and methods: Eighteen candidate genes were chosen for their putative roles in exercise. Re-sequencing in Thoroughbred samples was successful for primer sets in 13 of these genes. SNPs identified in this study and from the EquCab2.0 SNP database were genotyped in 2 sets of Thoroughbred samples (n = 150 and 148) and a series of population-based case-control investigations were performed by separating the samples into discrete cohorts on the basis of retrospective racecourse performance. Results: Twenty novel SNPs were detected in 3 genes: ACTN3, CKM and COX4I2. Genotype frequency distributions for 3 SNPs in CKM and COX4I2 were significantly (P<0.05) different between elite Thoroughbreds and racehorses that had never won a race. These associations were not validated when an additional (n = 130) independent set of samples was genotyped, but when analyses included all samples (n = 278) the significance of association at COX4I2 g.22684390C>T was confirmed (P<0.02). Conclusions: While molecular genetic information has the potential to become a powerful tool to make improved decisions in horse industries, it is vital that rigour is applied to studies generating these data and that adequate and appropriate sample sets, particularly for independent replication, are used.

Sequence variation at the equine myostatin gene (MSTN) locus has previously been shown to have a singular genomic influence on optimum race distance in Thoroughbred racehorses. Myostatin, encoded by the MSTN gene, is a member of the TGF-β superfamily that regulates skeletal muscle development in a range of mammalian species including the horse. In the Thoroughbred, the C-allele at the g.66493737C/T SNP has been found at significantly higher frequency in subgroups of the population that are suited to fast, short distance, sprint races and also influences body composition phenotypes. We investigated the influence of the g.66493737C/T SNP on speed indexes measured in a cohort of n = 85 Thoroughbred horses-in-training. We found significant associations between genotypes at the g.66493737C/T SNP and all measured speed variables: Dist(6) [distance travelled during 6 s before and after maximal velocity (V(max)); P = 0.0040], V(maxt) (duration at V(max); P = 0.0249), V(max) (P = 0.0265), Dist(6b) (distance travelled during 6 s before V(max); P = 0.0032), and Dist(6a)(distance travelled during 6 s after V(max); P = 0.0317). For each measure, horses with the C/C and C/T genotypes outperformed T/T horses, indicating the requirement for at least one C-allele to improve speed. For the most significantly associated variables (Dist(6) and Dist(6b)) the C/C cohort performed better than the T/T cohort with the heterozygotes intermediate, indicating a dose-dependent manifestation. These findings clearly indicate that variation at the MSTN gene influences speed in Thoroughbred horses.

SummaryThe recent completion of the horse genome and commercial availability of an equine SNP genotyping array has facilitated the mapping of disease genes. We report putative localization of the gene responsible for dwarfism, a trait in Friesian horses that is thought to have a recessive mode of inheritance, to a 2-MB region of chromosome 14 using just 10 affected animals and 10 controls. We successfully genotyped 34 429 SNPs that were tested for association with dwarfism using chi-square tests. The most significant SNP in our study, BIEC2-239376 (P 2df = 4.54 · 10 )5 , P rec = 7.74 · 10), is located close to a gene implicated in human dwarfism. Fine-mapping and resequencing analyses did not aid in further localization of the causative variant, and replication of our findings in independent sample sets will be necessary to confirm these results.