Identification of the myostatin locus (MSTN) as having a major effect on optimum racing distance in the Thoroughbred horse in the USA

Unveiling Genetic Potential for Equine Meat Production: A Bioinformatics Approach

In view of the predicted significant increase in global meat production, alternative sources such as horsemeat are becoming increasingly important due to their lower environmental impact and high nutritional value. This study aimed to identify SNP markers on the GeneSeek® Genomic Profiler™ Equine (Neogen, Lansing, MI, USA) that are important for horsemeat production traits. First, orthologous genes related to meat yield in cattle and common genes between horses and cattle within QTLs for body size and weight were identified. Markers for these genes were then evaluated based on predicted variant consequences, GERP scores, and positions within constrained elements and orthologous regulatory regions in pigs. A total of 268 markers in 57 genes related to meat production were analyzed. This resulted in 27 prioritized SNP markers in 22 genes, including notable markers in LCORL, LASP1, IGF1R, and MSTN. These results will benefit smallholder farmers by providing genetic insights for selective breeding that could improve meat yield. This study also supports future large-scale genetic analyses such as GWAS and Genomic Best Linear Unbiased Prediction (GBLUP). The results of this study may be helpful in improving the accuracy of genomic breeding values. However, limitations include reliance on bioinformatics without experimental validation. Future research can validate these markers and consider a wider range of traits to ensure accuracy in equine breeding.

Genomic tools for early selection among Thoroughbreds and Polo Argentino horses for practicing polo

The Polo Argentino (PA) horse is a recognized breed, developed originally by mixing crossbred and Thoroughbred (TB) horses to play polo. Early PA selection is difficult due to unreliable performance estimations. This study investigated the usefulness of genomic markers previously linked to morphological and functional traits as a tool for the early selection of PA. To this, we genotyped 520 PA and 30 TB horses using the Equine GGPArray (Illumina, n = 71,778 SNPs). Analyses included a genetic characterization of six genetic markers associated with behavioral (DRD4), muscular development (MSTN), and body size (LCORL, HMGA6, ZFAT, and LASP1) genes. Genetic differences in the DRD4, MSTN, and LCORL SNP were found between the two breeds, in the last two FST index between breeds was 0.13 and 0.6, respectively (p < 0.01). In DRD4, G allele was the more prevalent in PA (0.56 vs 0.45 in TB, p < 0.05), but no differences were observed between the genotypes associated with phenotypes. In MSTN, heterozygous genotypes were the most common in PA (48 %), with a significant decrease in AA (Hardy-Weinberg p < 0.05), suggesting a negative selection against it in polo horses. In body size, HMGA2 was monomorphic in all horses, while ZFAT and LASP1 SNP showed higher variability. Interestingly, 99 % of PA showed a TT genotype in LCORL (only 66 % in TB), demonstrating selection for smaller horses. Our results suggest that empirical selection in PA has generated an incipient genomic differentiation in discrete traits which could be used as a marker-assisted selection tool for early selection of polo horses.

Conditionally immortalised equine skeletal muscle cell lines for in vitro analysis

Background

Thoroughbred racehorse performance is largely influenced by a major quantitative trait locus at the myostatin (MSTN) gene which determines aptitude for certain race distances due to a promoter region insertion mutation influencing functional phenotypes in skeletal muscle. To develop an in vitro system for functional experiments we established three novel equine skeletal muscle cell lines reflecting the variation in phenotype associated with MSTN genotype (CC/II, CT/IN and TT/NN for SNP g.66493737C > T/SINE insertion 227 bp polymorphism). Primary equine skeletal muscle myoblasts, isolated from Thoroughbred horse gluteus medius, were conditionally immortalised and evaluated to determine whether cell phenotype and metabolic function were comparable to functional characteristics previously reported for ex vivo skeletal muscle isolated from Thoroughbred horses with each genotype.

Results

Primary myoblasts conditionally immortalised with the temperature sensitive SV40TtsA58 lentivirus vector successfully proliferated and could revert to their primary cell phenotype and differentiate into multinucleated myotubes. Skeletal muscle fibre type, MSTN gene expression, mitochondrial abundance, and mitochondrial function of the three MSTN genotype cell lines, were consistent with equivalent characterisation of ex vivo skeletal muscle samples with these genotypes. Furthermore, addition of coenzyme Q₁₀ (CoQ₁₀) to the cell lines improved mitochondrial function, an observation consistent with ex vivo skeletal muscle samples with these genotypes following supplementation with CoQ₁₀ in the diet.

Conclusions

The observation that the phenotypic characteristics and metabolic function of the cells lines are equivalent to ex vivo skeletal muscle indicates that this in vitro system will enable efficient and cost-effective analyses of equine skeletal muscle for a range of different applications including understanding metabolic function, testing of nutritional supplements, drug test development and gene doping test development. In the multi-billion-euro international Thoroughbred horse industry research advances in the biological function of skeletal muscle are likely to have considerable impact. Furthermore, this novel genotype-specific system may be adapted and applied to human biomedicine to improve understanding of the effects of myostatin in human physiology and medicine.

Locomotory Profiles in Thoroughbreds: Peak Stride Length and Frequency in Training and Association with Race Outcomes

Racehorses competing in short (i.e., ‘sprinters’), middle- or longer-distance (i.e., ‘stayers’) flat races are assumed to have natural variation in locomotion; sprinters having an innately shorter stride than stayers. No study has objectively tested this theory. Here, racehorses (n = 421) were categorised as sprinters, milers or stayers based on known race distance (n = 3269 races). Stride parameters (peak length and frequency) of those racehorses were collected from prior race-pace training sessions on turf (n = 2689; ‘jumpout’, n = 1013), using a locomotion monitoring device. Pedigree information for all 421 racehorses was extracted to three-generations. In training, sprinters had a shorter stride of higher frequency and covered consecutive furlongs faster than stayers (p < 0.001). Relatively short or longer stride did not predict race success, but stayers had greater race success than sprinters (p < 0.001). Peak stride length and frequency were moderately heritable (h2 = 0.15 and 0.20, respectively). In conclusion, differences in stride were apparent between sprinters and stayers (e.g., shorter stride in sprinters) during routine training, even after accounting for their pedigree. Objective data on stride characteristics could supplement other less objectively obtained parameters to benefit trainers in the appropriate selection of races for each individual racehorse.

Genetic Parameters of Performance and Conformation Traits of 3-Year-Old Warmblood Sport Horses in the Czech Republic

The aim of this study was to estimate the genetic parameters of a one-day performance test together with the linear type traits of 3-year-old warmblood horses. The study of genetic parameters was based on 5958 tested horses in the period 1998-2021. A total of 22 traits of linear description, three quantitatively measured traits, and one summary mark from the performance test were tested. The model equation included the fixed effect of gender and combination effects of classifier-year of evaluation-place. A single-trait animal model was used for the estimation of heritability and genetic variance, while the two-trait animal model was applied for the estimation of variance and covariance between all traits. The heritability of the overall score of the performance test was 0.25. The range for heritability was between 0.04 and 0.33 for the linear type traits and between 0.46 and 0.57 for the quantitatively measured traits. Genetic correlations were between -0.47 and 0.92. The estimated genetic parameters suggest that the results from the performance test can be incorporated into genetic evaluation in the Czech Republic.

MSTN Regulatory Network in Mongolian Horse Muscle Satellite Cells Revealed with miRNA Interference Technologies

Myostatin (MSTN), a member of the transforming growth factor-β superfamily, inhibits the activation of muscle satellite cells. However, the role and regulatory network of MSTN in equine muscle cells are not well understood yet. We discovered that MSTN knockdown significantly reduces the proliferation rate of equine muscle satellite cells. In addition, after the RNA sequencing of equine satellite cells transfected with MSTN-interference plasmid and control plasmid, an analysis of the differentially expressed genes was carried out. It was revealed that MSTN regulatory networks mainly involve genes related to muscle function and cell-cycle regulation, and signaling pathways, such as Notch, MAPK, and WNT. Subsequent real-time PCR in equine satellite cells and immunohistochemistry on newborn and adult muscle also verified the MSTN regulatory network found in RNA sequencing analysis. The results of this study provide new insight into the regulatory mechanism of equine MSTN.

Longitudinal assessment of blood-borne musculoskeletal disease biomarkers in the DE50-MD dog model of Duchenne muscular dystrophy

Background: Duchenne muscular dystrophy (DMD) is a fatal muscle wasting disease caused by mutations in the dystrophin gene. Due to their phenotypic similarity to human patients, large animal models are invaluable tools for pre-clinical trials. The DE50-MD dog is a relatively new model of DMD, and carries a therapeutically-tractable mutation lying within the hotspot for human patients, making it especially valuable. Prior to conducting therapeutic trials using this novel animal model, it is essential to establish a panel of viable biomarkers. Methods: We evaluated a panel of blood-borne biomarkers of musculoskeletal disease in the DE50-MD dog. Venous blood samples were obtained monthly throughout an 18-month study period in DE50-MD (N=18) and wild-type (WT) control (N=14) dogs. A panel of potential plasma/serum biomarkers of DMD was measured and their theoretical utility in future clinical trials determined using sample size calculations. Results: Compared to WT dogs, DE50-MD dogs had substantially higher circulating creatine kinase (CK) activities, myomesin-3 (MYOM3), and the dystromiRs miR-1, miR-133a and miR-206, but significantly lower serum myostatin concentrations. An age-associated pattern, similar to that observed in DMD patients, was seen for CK and MYOM3. Sample size calculations suggested that low cohort sizes (N≤3) could be used to detect up to a 50% improvement in DE50-MD results towards WT levels for each biomarker or a combination thereof (via principal component analysis); as few as N=3 animals should enable detection of a 25% improvement using a combined biomarker approach (alpha 0.05, power 0.8). Conclusions: We have established a panel of blood-borne biomarkers that could be used to monitor musculoskeletal disease or response to a therapeutic intervention in the DE50-MD dog using low numbers of animals. The blood biomarker profile closely mimics that of DMD patients, supporting the hypothesis that this DMD model would be suitable for use in pre-clinical trials.

Longitudinal assessment of blood-borne musculoskeletal disease biomarkers in the DE50-MD dog model of Duchenne muscular dystrophy

Background: Duchenne muscular dystrophy (DMD) is a fatal muscle wasting disease caused by mutations in the dystrophin gene. Due to their phenotypic similarity to human patients, large animal models are invaluable tools for pre-clinical trials. The DE50-MD dog is a relatively new model of DMD, and carries a therapeutically-tractable mutation lying within the hotspot for human patients, making it especially valuable. Prior to conducting therapeutic trials using this novel animal model, it is essential to establish a panel of viable biomarkers. Methods: We evaluated a panel of blood-borne biomarkers of musculoskeletal disease in the DE50-MD dog. Venous blood samples were obtained monthly throughout an 18-month study period in DE50-MD (N=18) and wild-type (WT) control (N=14) dogs. A panel of potential plasma/serum biomarkers of DMD was measured and their theoretical utility in future clinical trials determined using sample size calculations. Results: Compared to WT dogs, DE50-MD dogs had substantially higher circulating creatine kinase (CK) activities, myomesin-3 (MYOM3), and the dystromiRs miR-1, miR-133a and miR-206, but significantly lower serum myostatin concentrations. An age-associated pattern, similar to that observed in DMD patients, was seen for CK and MYOM3. Sample size calculations suggested that low cohort sizes (N≤3) could be used to detect up to a 50% improvement in DE50-MD results towards WT levels for each biomarker or a combination thereof (via principal component analysis); as few as N=3 animals should enable detection of a 25% improvement using a combined biomarker approach (alpha 0.05, power 0.8). Conclusions: We have established a panel of blood-borne biomarkers that could be used to monitor musculoskeletal disease or response to a therapeutic intervention in the DE50-MD dog using low numbers of animals. The blood biomarker profile closely mimics that of DMD patients, supporting the hypothesis that this DMD model would be suitable for use in pre-clinical trials.

Genetics of Thoroughbred Racehorse Performance

Thoroughbred horses have been selected for racing performance for more than 400 years. Despite continued selection, race times have not improved significantly during the past 60 years, raising the question of whether genetic variation for racing performance still exists. Studies using phenotypes such as race time, money earned, and handicapping, however, demonstrate that there is extensive variation within these traits and that they are heritable. Even so, these are poor measures of racing success since Thoroughbreds race at different ages and distances and on different types of tracks, and some may not race at all. With the advent of genomic tools, DNA variants are being identified that contribute to racing success. Aside from strong associations for myostatin variants with best racing distance, weak to modest associations with racing phenotypes are reported for other genomic regions. These data suggest that diverse genetic strategies have contributed to producing a successful racehorse, and genetic variation contributing to athleticism remains important. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Role of genes related to performance and reproduction of Thoroughbreds in training and breeding - A review

Thoroughbreds have been selected for speed and stamina since the 1700s. This selection resulted in structural and functional system-wide adaptations that enhanced physiological characteristics for outstanding speed of 61-71 kph (38-44 mph) between 1,000 and 3,200 m (5 furlongs - 2 miles). At present, horseracing is still an economically important industrial sector, therefore intensive research is underway to explore genes that allow the utilisation of genetic abilities and are significant in breeding and training. This study aims to provide an overview of genetic research and its applicability related to Thoroughbreds.

Generation of myostatin edited horse embryos using CRISPR/Cas9 technology and somatic cell nuclear transfer

The application of new technologies for gene editing in horses may allow the generation of improved sportive individuals. Here, we aimed to knock out the myostatin gene (MSTN), a negative regulator of muscle mass development, using CRISPR/Cas9 and to generate edited embryos for the first time in horses. We nucleofected horse fetal fibroblasts with 1, 2 or 5 µg of 2 different gRNA/Cas9 plasmids targeting the first exon of MSTN. We observed that increasing plasmid concentrations improved mutation efficiency. The average efficiency was 63.6% for gRNA1 (14/22 edited clonal cell lines) and 96.2% for gRNA2 (25/26 edited clonal cell lines). Three clonal cell lines were chosen for embryo generation by somatic cell nuclear transfer: one with a monoallelic edition, one with biallelic heterozygous editions and one with a biallelic homozygous edition, which rendered edited blastocysts in each case. Both MSTN editions and off-targets were analyzed in the embryos. In conclusion, CRISPR/Cas9 proved an efficient method to edit the horse genome in a dose dependent manner with high specificity. Adapting this technology sport advantageous alleles could be generated, and a precision breeding program could be developed.

Genetic Testing in the Horse

Genetic testing in horses began in the 1960s, when parentage testing using blood group markers became the standard. In the 1990s, parentage testing shifted from evaluating blood groups to DNA testing. The development of genetics and genomics in both human and veterinarian medicine, along with continued technological advances in the last 2 decades, has helped unravel the causal variants for many horse traits. Genetic testing is also now possible for a variety of phenotypic and disease traits and is used to assist in breeding and clinical management decisions. This article describes the genetic tests that are currently available for horses.

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.

The use of the SLC16A1 gene as a potential marker to predict race performance in Arabian horses

BACKGROUND

Arabian horses are commonly believed to be one of the oldest and the most popular horse breeds in the world, characterized by favourable stamina traits and exercise phenotypes. During intensive training, the rates of lactate production and utilization are critical to avoid muscle fatigue and a decrease in exercise performance. The key factor determining transmembrane lactate transport is the monocarboxylate transporter 1 protein coded for by the SLC16A1 gene. The aim of the present research was to identify polymorphisms in the coding sequence and UTRs in the equine SLC16A1 gene and to evaluate their potential association with race performance traits in Arabian horses. Based on RNA-seq data, SNPs were identified and genotyped using PCR-RFLP or PCR-HRM methods in 254 Arabian horses that competed in flat races. An association analysis between polymorphisms and racing results was performed.

RESULTS

Novel polymorphisms in the equine SLC16A1 locus have been identified (missense and 5'UTR variants: g.55601543C > T and g.55589063 T > G). Analysis showed a significant association between the 5'UTR polymorphism and several racing results as follows: the possibility of winning first or second place, the number of races in which horses started and total financial benefits. The analysis also showed differences in genotype distribution depending on race distance. In the studied population, the shorter distance races were only won by TT horses. The GG and TG horses took first and second places in middle- and long-distance races, and the percentage of winning heterozygotes increased from 19.5 to 27% at the middle and long distances, respectively. The p.Val432Ile (g.55601543C > T) polymorphism was not significantly related to the analysed racing results.

CONCLUSION

Our results showed that g.55589063 T > G polymorphism affected the possibility of winning first or second place and of competing in more races. The different distribution of genotypes depending on race distance indicated the possibility of using a SNP in the SLC16A1 gene as a marker to predict the best race distance for a horse. The presented results provide a basis for further research to validate the use of the SLC16A1 gene as a potential marker associated with racing performance.

Whole genome sequencing analysis of horse populations inhabiting the Korean Peninsula and Przewalski's horse

BACKGROUND

The Jeju horse is an indigenous horse breed in Korea. However, there is a severe lack of genomic studies on Korean horse breeds.

OBJECTIVE

The objective of this study was to report genomic characteristics of domestic horse populations that inhabit South Korea (Jeju, Jeju crossbred, and Thoroughbred) and a wild horse breed (Przewalski's horse).

RESULTS

Using the equine reference genome assembly (EquCab 2.0), more than ~ 6.5 billion sequence reads were successfully mapped, which generated an average of 40.87-fold coverage throughout the genome. Using these data, we detected a total of 12.88 million SNPs, of which 73.7% were found to be novel. All the detected SNPs were deeply annotated to retrieve SNPs in gene regions using the RefSeq and Ensemble gene sets. Approximately 27% of the total SNPs were located within genes, whereas the remaining 73% were found in intergenic regions. Using 129,776 coding SNPs, we retrieved a total of 49,171 nonsynonymous SNPs in 12,351 genes. Furthermore, we identified a total of 10,770 deleterious nonsynonymous SNPs which are predicted to affect protein structure or function.

CONCLUSION

We showed numerous genomic variants from domestic and wild horse breeds. These results provide a valuable resource for further studies on functions of SNP-containing genes, and can aid in determining the molecular basis underlying variation in economically important traits of horses.

The "speed gene" effect of myostatin arises in Thoroughbred horses due to a promoter proximal SINE insertion

Thoroughbred horses are finely-tuned athletes with a high aerobic capacity relative to skeletal muscle mass, attributable to centuries of genetic selection for speed and stamina. Polymorphisms in the myostatin gene (MSTN), a pronounced inhibitor of skeletal muscle growth, have been shown to almost singularly account for gene-based race distance aptitude in racehorses. In Thoroughbreds, two MSTN polymorphisms, a single nucleotide variation in the first intron (SNP g.66493737C>T) and a non-coding transposable element within the promoter region (a 227 bp SINE insertion) are of particular interest. Until now, it has not been clear which of these variants affect skeletal muscle phenotypes or whether both can impact racing performance. In a large cohort of Thoroughbreds, we observed a complete concordance between the SNP and the SINE insertion. By means of in vitro assays in C₂C₁₂ myoblasts, we isolated the SNP variant from the SINE polymorphism and showed the latter is exclusively responsible for adversely affecting transcription initiation and gene expression thereby limiting myostatin protein production. Mapping the MSTN transcription start site in horse skeletal muscle likewise revealed anomalous transcription initiation in the presence of the SINE insertion. Our data provides mechanistic evidence that the SINE insertion uniquely accounts for the MSTN “speed gene” effect on race distance aptitude in the Thoroughbred horse.

A genome-wide association study for harness racing success in the Norwegian-Swedish coldblooded trotter reveals genes for learning and energy metabolism

Background

Although harness racing is of high economic importance to the global equine industry, significant genomic resources have yet to be applied to mapping harness racing success. To identify genomic regions associated with harness racing success, the current study performs genome-wide association analyses with three racing performance traits in the Norwegian-Swedish Coldblooded Trotter using the 670 K Axiom Equine Genotyping Array.

Results

Following quality control, 613 horses and 359,635 SNPs were retained for further analysis. After strict Bonferroni correction, nine genome-wide significant SNPs were identified for career earnings. No genome-wide significant SNPs were identified for number of gallops or best km time. However, four suggestive genome-wide significant SNPs were identified for number of gallops, while 19 were identified for best km time. Multiple genes related to intelligence, energy metabolism, and immune function were identified as potential candidate genes for harness racing success.

Conclusions

Apart from the physiological requirements needed for a harness racing horse to be successful, the results of the current study also advocate learning ability and memory as important elements for harness racing success. Further exploration into the mental capacity required for a horse to achieve racing success is likely warranted.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0670-3) contains supplementary material, which is available to authorized users.

A potential regulatory region near the EDN3 gene may control both harness racing performance and coat color variation in horses

The Swedish-Norwegian Coldblooded trotter and the heavier North-Swedish draught horse both descend from the North-Swedish horse, but the Coldblooded trotters have been selected for racing performance while the North-Swedish draught horse is mainly used for agricultural and forestry work. By comparing the genomes of Coldblooded trotters, North-Swedish draught horses and Standardbreds for a large number of single-nucleotide polymorphisms (SNPs), the aim of the study was to identify genetic regions that may be under selection for racing performance. We hypothesized that the selection for racing performance, in combination with unauthorized crossbreeding of Coldblooded trotters and Standardbreds, has created regions in the genome where the Coldblooded trotters and Standardbreds are similar, but differ from the North-Swedish draught horse. A fixation index (Fst) analysis was performed and sliding window Delta Fst values were calculated across the three breeds. Five windows, where the average Fst between Coldblooded trotters and Standardbreds was low and the average Fst between Coldblooded trotters and North-Swedish draught horses was high, were selected for further investigation. Associations between the most highly ranked SNPs and harness racing performance were analyzed in 400 raced Coldblooded trotters with race records. One SNP showed a significant association with racing performance, with the CC genotype appearing to be negatively associated. The SNP identified was genotyped in 1915 horses of 18 different breeds. The frequency of the TT genotype was high in breeds typically used for racing and show jumping while the frequency of the CC genotype was high in most pony breeds and draught horses. The closest gene in this region was the Endothelin3 gene (EDN3), a gene mainly involved in melanocyte and enteric neuron development. Both functional genetic and physiological studies are needed to fully understand the possible impacts of the gene on racing performance.

Effect of selection for eventing on the MSTN gene in Brazilian sport horses

Polymorphisms in MSTN have previously been associated with equine performance. Therefore, the aim of this study was to identify variants in MSTN intron 1 in 16 Brazilian Sport Horses selected for competition in eventing and their possible effects of selection on performance. Among the nine variants identified, eight had already been reported in previous studies or genomic databases, although they showed differences in frequencies when compared with other horse breeds. Moreover, a new mutation was identified in two horses, both in heterozygous form. Considering the absence of molecular studies in this valuable Brazilian breed, these findings represent an important contribution to the characterization of its genetic profile and may possibly aid in further genotype-phenotype association studies.

Skeletal muscle mitochondrial bioenergetics and associations with myostatin genotypes in the Thoroughbred horse

Variation in the myostatin (MSTN) gene has been reported to be associated with race distance, body composition and skeletal muscle fibre composition in the horse. The aim of the present study was to test the hypothesis that MSTN variation influences mitochondrial phenotypes in equine skeletal muscle. Mitochondrial abundance and skeletal muscle fibre types were measured in whole muscle biopsies from the gluteus medius of n = 82 untrained (21 ± 3 months) Thoroughbred horses. Skeletal muscle fibre type proportions were significantly (p < 0.01) different among the three MSTN genotypes and mitochondrial content was significantly (p < 0.01) lower in the combined presence of the C-allele of SNP g.66493737C>T (C) and the SINE insertion 227 bp polymorphism (I). Evaluation of mitochondrial complex activities indicated higher combined mitochondrial complex I+III and II+III activities in the presence of the C-allele / I allele (p ≤ 0.05). The restoration of complex I+III and complex II+III activities following addition of exogenous coenzyme Q1 (ubiquinone1) (CoQ1) in vitro in the TT/NN (homozygous T allele/homozygous no insertion) cohort indicated decreased coenzyme Q in these animals. In addition, decreased gene expression in two coenzyme Q (CoQ) biosynthesis pathway genes (COQ4, p ≤ 0.05; ADCK3, p ≤ 0.01) in the TT/NN horses was observed. This study has identified several mitochondrial phenotypes associated with MSTN genotype in untrained Thoroughbred horses and in addition, our findings suggest that nutritional supplementation with CoQ may aid to restore coenzyme Q activity in TT/NN horses.

Equine skeletal muscle adaptations to exercise and training: evidence of differential regulation of autophagosomal and mitochondrial components

BACKGROUND

A single bout of exercise induces changes in gene expression in skeletal muscle. Regular exercise results in an adaptive response involving changes in muscle architecture and biochemistry, and is an effective way to manage and prevent common human diseases such as obesity, cardiovascular disorders and type II diabetes. However, the biomolecular mechanisms underlying such responses still need to be fully elucidated. Here we performed a transcriptome-wide analysis of skeletal muscle tissue in a large cohort of untrained Thoroughbred horses (n = 51) before and after a bout of high-intensity exercise and again after an extended period of training. We hypothesized that regular high-intensity exercise training primes the transcriptome for the demands of high-intensity exercise.

RESULTS

An extensive set of genes was observed to be significantly differentially regulated in response to a single bout of high-intensity exercise in the untrained cohort (3241 genes) and following multiple bouts of high-intensity exercise training over a six-month period (3405 genes). Approximately one-third of these genes (1025) and several biological processes related to energy metabolism were common to both the exercise and training responses. We then developed a novel network-based computational analysis pipeline to test the hypothesis that these transcriptional changes also influence the contextual molecular interactome and its dynamics in response to exercise and training. The contextual network analysis identified several important hub genes, including the autophagosomal-related gene GABARAPL1, and dynamic functional modules, including those enriched for mitochondrial respiratory chain complexes I and V, that were differentially regulated and had their putative interactions 're-wired' in the exercise and/or training responses.

CONCLUSION

Here we have generated for the first time, a comprehensive set of genes that are differentially expressed in Thoroughbred skeletal muscle in response to both exercise and training. These data indicate that consecutive bouts of high-intensity exercise result in a priming of the skeletal muscle transcriptome for the demands of the next exercise bout. Furthermore, this may also lead to an extensive 're-wiring' of the molecular interactome in both exercise and training and include key genes and functional modules related to autophagy and the mitochondrion.

The Evolutionary Origin and Genetic Makeup of Domestic Horses

The horse was domesticated only 5.5 KYA, thousands of years after dogs, cattle, pigs, sheep, and goats. The horse nonetheless represents the domestic animal that most impacted human history; providing us with rapid transportation, which has considerably changed the speed and magnitude of the circulation of goods and people, as well as their cultures and diseases. By revolutionizing warfare and agriculture, horses also deeply influenced the politico-economic trajectory of human societies. Reciprocally, human activities have circled back on the recent evolution of the horse, by creating hundreds of domestic breeds through selective programs, while leading all wild populations to near extinction. Despite being tightly associated with humans, several aspects in the evolution of the domestic horse remain controversial. Here, we review recent advances in comparative genomics and paleogenomics that helped advance our understanding of the genetic foundation of domestic horses.

Association of low race performance with mtDNA haplogroup L3b of Australian thoroughbred horses

Mitochondrial DNA (mtDNA) encodes the genes for respiratory chain sub-units that determine the efficiency of oxidative phosphorylation in mitochondria. The aim of this study was to determine if there were any haplogroups and variants in mtDNA that could be associated with athletic performance of Thoroughbred horses. The whole mitochondrial genomes of 53 maternally unrelated Australian Thoroughbred horses were sequenced and an association study was performed with the competition histories of 1123 horses within their maternal lineages. A horse mtDNA phylogenetic tree was constructed based on a total of 195 sequences (including 142 from previous reports). The association analysis showed that the sample groups with poor racing performance history were enriched in haplogroup L3b (p = .0003) and its sub-haplogroup L3b1a (p = .0007), while those that had elite performance appeared to be not significantly associated with haplogroups G2 and L3a1a1a (p > .05). Haplogroup L3b and L3b1a bear two and five specific variants of which variant T1458C (site 345 in 16s rRNA) is the only potential functional variant. Furthermore, secondary reconstruction of 16s RNA showed considerable differences between two types of 16s RNA molecules (with and without T1458C), indicating a potential functional effect. The results suggested that haplogroup L3b, could have a negative association with elite performance. The T1458C mutation harboured in haplogroup L3b could have a functional effect that is related to poor athletic performance.

Conformation Traits and Gaits in the Icelandic Horse are Associated with Genetic Variants in Myostatin (MSTN)

Many genes are known to have an influence on conformation and performance traits; however, the role of one gene, Myostatin (MSTN), has been highlighted in recent studies on horses. Myostatin acts as a repressor in the development and regulation of differentiation and proliferative growth of skeletal muscle. Several studies have examined the link between MSTN, conformation, and performance in racing breeds, but no studies have investigated the relationship in Icelandic horses. Icelandic horses, a highly unique breed, are known both for their robust and compact conformation as well as their additional gaits tölt and pace. Three SNPs (g.65868604G>T [PR8604], g.66493737C>T [PR3737], and g.66495826A>G [PR5826]) flanking or within equine MSTN were genotyped in 195 Icelandic horses. The SNPs and haplotypes were analyzed for association with official estimated breeding values (EBV) for conformation traits (n = 11) and gaits (n = 5). The EBV for neck, withers, and shoulders was significantly associated with both PR8604 and PR3737 (P < 0.05). PR8604 was also associated with EBV for total conformation (P = 0.05). These associations were all supported by the haplotype analysis. However, while SNP PR5826 showed a significant association with EBVs for leg stance and hooves (P < 0.05), haplotype analyses for these traits failed to fully support these associations. This study demonstrates the possible role of MSTN on both the form and function of horses from non-racing breeds. Further analysis of Icelandic horses as well as other non-racing breeds would be beneficial and likely help to completely understand the influence of MSTN on conformation and performance in horses.

A Genome-wide Scan for Selective Sweeps in Racing Horses

Using next-generation sequencing, we conducted a genome-wide scan of selective sweeps associated with selection toward genetic improvement in Thoroughbreds. We investigated potential phenotypic consequence of putative candidate loci by candidate gene association mapping for the finishing time in 240 Thoroughbred horses. We found a significant association with the trait for Ral GApase alpha 2 (RALGAP2) that regulates a variety of cellular processes of signal trafficking. Neighboring genes around RALGAP2 included insulinoma-associated 1 (INSM1), pallid (PLDN), and Ras and Rab interactor 2 (RIN2) genes have similar roles in signal trafficking, suggesting that a co-evolving gene cluster located on the chromosome 22 is under strong artificial selection in racehorses.

Genome-wide evolutionary and functional analysis of the Equine Repetitive Element 1: an insertion in the myostatin promoter affects gene expression

BACKGROUND

In mammals, an important source of genomic variation is insertion polymorphism of retrotransposons. These may acquire a functional role when inserted inside genes or in their proximity. The aim of this work was to carry out a genome wide analysis of ERE1 retrotransposons in the horse and to analyze insertion polymorphism in relation to evolution and function. The effect of an ERE1 insertion in the promoter of the myostatin gene, which is involved in muscle development, was also investigated.

RESULTS

In the horse population, the fraction of ERE1 polymorphic loci is related to the degree of similarity to their consensus sequence. Through the analysis of ERE1 conservation in seven equid species, we established that the level of identity to their consensus is indicative of evolutionary age of insertion. The position of ERE1s relative to genes suggests that some elements have acquired a functional role. Reporter gene assays showed that the ERE1 insertion within the horse myostatin promoter affects gene expression. The frequency of this variant promoter correlates with sport aptitude and racing performance.

CONCLUSIONS

Sequence conservation and insertion polymorphism of ERE1 elements are related to the time of their appearance in the horse lineage, therefore, ERE1s are a useful tool for evolutionary and population studies. Our results suggest that the ERE1 insertion at the myostatin locus has been unwittingly selected by breeders to obtain horses with specific racing abilities. Although a complex combination of environmental and genetic factors contributes to athletic performance, breeding schemes may take into account ERE1 insertion polymorphism at the myostatin promoter.

RNA sequencing as a powerful tool in searching for genes influencing health and performance traits of horses

RNA sequencing (RNA-seq) by next-generation technology is a powerful tool which creates new possibilities in whole-transcriptome analysis. In recent years, with the use of the RNA-seq method, several studies expanded transcriptional gene profiles to understand interactions between genotype and phenotype, supremely contributing to the field of equine biology. To date, in horses, massive parallel sequencing of cDNA has been successfully used to identify and quantify mRNA levels in several normal tissues, as well as to annotate genes. Moreover, the RNA-seq method has been applied to identify the genetic basis of several diseases or to investigate organism adaptation processes to the training conditions. The use of the RNA-seq approach has also confirmed that horses can be useful as a large animal model for human disease, especially in the field of immune response. The presented review summarizes the achievements of profiling gene expression in horses (Equus caballus).

Putative regulation mechanism for the MSTN gene by a CpG island generated by the SINE marker Ins227bp

Background

A single nucleotide polymorphism (SNP) in the first intron of the myostatin gene (MSTN) is associated with aptness of elite Thoroughbreds to race over sprint, middle or long distances. This intronic marker (g.66493737 T ≻ C), a short interspersed nuclear element (SINE) of 227 bp (Ins227bp) insertion polymorphism in the MSTN promoter, and the adjacent SNP BIEC2-417495 have not been studied for their association with racing aptness of the average Thoroughbreds raced in countries with lower status of the racing industry. This study investigated these markers regarding their prevalence and association with performance in common race horses. Markers were genotyped by amplification refractory mutation system-quantitative PCR (ARMS-qPCR) or amplicon melting. Furthermore, we asked whether the Ins227bp marker might theoretically regulate the expression of myostatin by generating a novel target for DNA methylation or by changing binding sites for transcription factors. Putative sites for DNA methylation or binding of transcription factors were predicted by MethPrimer and by the softwares JASPAR, MatInspector and UniPROBE, respectively.

Results

Pairwise linkage disequilibrium between g.66493737 T ≻ C and Ins227bp was high (r² = 0.93). A lower linkage was determined for g.66493737 T ≻ C and BIEC2-417495 (r² = 0.69) as well as for BIEC2-417495 and Ins227bp (r² = 0.76). The estimated frequencies for the presence of Ins227bp (I) indel and the C alleles at g.66493737 T ≻ C and BIEC2-417495 were 0.46, 0.47 and 0.43, respectively. Heterozygotes represented the most abundant genotype at each locus. The best racing distance (BRD) was significantly different between the homozygotes of each SNP (p = 0.01 to 0.03). C allele homozygotes at BIEC2-417495 or g.66493737 T ≻ C, as well as Ins227bp homozygotes earned most money on a mean distance ranging from 1211 to 1230 m. Heterozygotes earned most money on races over 1690 to 1709 m. The BRD for the T/T carriers at both SNP loci and for the SINE-free genotype was 1812 to 1854 m. Other performance parameters were not significantly different between the genotypes, except of the relative success score (RSS). The RSS was significantly slightly better on a distance of ≤1300 m for all carriers of the C allele and the Ins227bp compared to homozygous T genotypes and SINE-negative horses (p = 0.037 to 0.046). For distances of more than 1300 m the RSS was not significantly different between genotypes.

In silico assessment indicated that the Ins227bp promoter insertion might have generated a CpG island and a few novel putative binding sites for transcription factors.

Conclusions

All three target polymorphisms (Ins227bp, g.66493737 T ≻ C, BIEC2-417495) are suitable markers to assess the ability of non-elite Thoroughbreds to race at short or longer distances. The CpG island generated by Ins227bp may cause training-induced silencing of MSTN expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-015-0428-3) contains supplementary material, which is available to authorized users.

Multiple Genes Related to Muscle Identified through a Joint Analysis of a Two-stage Genome-wide Association Study for Racing Performance of 1,156 Thoroughbreds

Thoroughbred, a relatively recent horse breed, is best known for its use in horse racing. Although myostatin (MSTN) variants have been reported to be highly associated with horse racing performance, the trait is more likely to be polygenic in nature. The purpose of this study was to identify genetic variants strongly associated with racing performance by using estimated breeding value (EBV) for race time as a phenotype. We conducted a two-stage genome-wide association study to search for genetic variants associated with the EBV. In the first stage of genome-wide association study, a relatively large number of markers (~54,000 single-nucleotide polymorphisms, SNPs) were evaluated in a small number of samples (240 horses). In the second stage, a relatively small number of markers identified to have large effects (170 SNPs) were evaluated in a much larger number of samples (1,156 horses). We also validated the SNPs related to MSTN known to have large effects on racing performance and found significant associations in the stage two analysis, but not in stage one. We identified 28 significant SNPs related to 17 genes. Among these, six genes have a function related to myogenesis and five genes are involved in muscle maintenance. To our knowledge, these genes are newly reported for the genetic association with racing performance of Thoroughbreds. It complements a recent horse genome-wide association studies of racing performance that identified other SNPs and genes as the most significant variants. These results will help to expand our knowledge of the polygenic nature of racing performance in Thoroughbreds.

Haplotype diversity in the equine myostatin gene with focus on variants associated with race distance propensity and muscle fiber type proportions

Two variants in the equine myostatin gene (MSTN), including a T/C SNP in the first intron and a 227-bp SINE insertion in the promoter, are associated with muscle fiber type proportions in the Quarter Horse (QH) and with the prediction of race distance propensity in the Thoroughbred (TB). Genotypes from these loci, along with 18 additional variants surrounding MSTN, were examined in 301 horses of 14 breeds to evaluate haplotype relationships and diversity. The C allele of intron 1 was found in 12 of 14 breeds at a frequency of 0.27; the SINE was observed in five breeds, but common in only the TB and QH (0.73 and 0.48 respectively). Haplotype data suggest the SINE insertion is contemporary to and arose upon a haplotype containing the intron 1 C allele. Gluteal muscle biopsies of TBs showed a significant association of the intron 1 C allele and SINE with a higher proportion of Type 2B and lower proportion of Type 1 fibers. However, in the Belgian horse, in which the SINE is not present, the intron 1 SNP was not associated with fiber type proportions, and evaluation of fiber type proportions across the Belgian, TB and QH breeds shows the significant effect of breed on fiber type proportions is negated when evaluating horses without the SINE variant. These data suggest the SINE, rather than the intron 1 SNP, is driving the observed muscle fiber type characteristics and is the variant targeted by selection for short-distance racing.

The DMRT3 'Gait keeper' mutation affects performance of Nordic and Standardbred trotters

In a previous study it was shown that a nonsense mutation in the DMRT3 gene alters the pattern of locomotion in horses and that this mutation has a strong positive impact on trotting performance of Standardbreds. One aim of this study was to test if racing performance and trotting technique in the Nordic (Coldblood) trotters are also influenced by the DMRT3 genotype. Another aim was to further investigate the effect of the mutation on performance in Standardbreds, by using a within-family analysis and genotype-phenotype correlations in a larger horse material than in the previous study. We genotyped 427 Nordic trotters and 621 Standardbreds for the DMRT3 nonsense mutation and a SNP in strong linkage disequilibrium with it. In Nordic trotters, we show that horses homozygous for the DMRT3 mutation (A) had significantly higher EBV for trotting performance traits than heterozygous (CA) or homozygous wild-type (CC) horses (P = 0.001). Furthermore, AA homozygotes had a higher proportion of victories and top 3 placings than horses heterozygous or homozygous wild-type, when analyzing performance data for the period 3 to 6 yr of age (P = 0.06 and P = 0.05, respectively). Another finding in the Nordic trotters was that the DMRT3 mutation influenced trotting technique (P = 2.1 × 10(-8)). Standardbred horses homozygous AA had significantly higher EBV for all traits than horses with at least 1 wild-type allele (CA and CC; P = 1.6 × 10(-16)). In a within-family analysis of Standardbreds, we found significant differences in several traits (e.g., earnings, P = 0.002; number of entered races, P = 0.004; and fraction of offspring that entered races, P = 0.002) among paternal half-sibs with genotype AA or CA sired by a CA stallion. For most traits, we found significant differences at young ages. For Nordic trotters, most of the results were significant at 3 yr of age but not for the older ages, and for the Standardbreds most of the results for the ages 3 to 5 were significant. For Nordic trotters, the proportion of victories and placings were the only traits that were significant for other ages than 3 yr.

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.

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.

A chromosomal region on ECA13 is associated with maxillary prognathism in horses

Hereditary variations in head morphology and head malformations are known in many species. The most common variation encountered in horses is maxillary prognathism. Prognathism and brachygnathism are syndromes of the upper and lower jaw, respectively. The resulting malocclusion can negatively affect teeth wear, and is considered a non-desirable trait in breeding programs. We performed a case-control analysis for maxillary prognathism in horses using 96 cases and 763 controls. All horses had been previously genotyped with a commercially available 50 k SNP array. We analyzed the data with a mixed-model considering the genomic relationships in order to account for population stratification. Two SNPs within a region on the distal end of chromosome ECA 13 reached the Bonferroni corrected genome-wide significance level. There is no known prognathism candidate gene located within this region. Therefore, our findings in the horse offer the possibility of identifying a novel gene involved in the complex genetics of prognathism that might also be relevant for humans and other livestock species.

Single nucleotide polymorphisms of myostatin gene in Chinese domestic horses

The myostatin gene (MSTN) is a genetic determinant of skeletal muscle growth. Single nucleotide polymorphisms (SNP) in MSTN are of importance due to their strong associations with horse racing performances. In this study, we screened the SNPs in MSTN gene in 514 horses from 15 Chinese horse breeds. Six SNPs (g.26T>C, g.156T>C, g.587A>G, g.598C>T, g.1485C>T, g.2115A>G) in MSTN gene were detected by sequencing and genotyped using PCR-RFLP method. The g.587A>G and g.598C>T residing in the 5'UTR region were novel SNPs identified by this study. The g.2115A>G which have previously been associated with racing performances were present in Chinese horse breeds, providing valuable genetic information for evaluating the potential racing performances in Chinese domestic breeds. The six SNPs together defined thirteen haplotypes, demonstrating abundant haplotype diversities in Chinese horses. Most of the haplotypes were shared among different breeds with no haplotype restricted to a specific region or a single horse breed. AMOVA analysis indicated that most of the genetic variance was attributable to differences among individuals without any significant contribution by the four geographical groups. This study will provide fundamental and instrumental genetic information for evaluating the potential racing performances of Chinese horse breeds.