Reviewe: Genetics and genomics in equine exercise physiology: an overview of the new applications of molecular biology as positive and negative markers of performance and health

Gene doping in horse racing and equine sports: Current landscape and future perspectives

Gene doping, the use of gene therapy or genetic manipulation to enhance athletic performance, has emerged as a potential threat to the integrity and welfare of equine sports, such as horse racing and equestrian sports. This review aims to provide an overview of gene doping in horses, including the underlying technologies, potential applications, detection methods, ethical concerns and future perspectives. By understanding the current landscape of gene doping in horses, stakeholders can work together to develop strategies to safeguard the integrity of equine sports.

Equine respiratory viruses, airway inflammation and performance in thoroughbred racehorses

Equine asthma is a common cause of poor performance in racehorses but it is unclear if respiratory viruses contribute to its etiology. The objective of the study was to determine if respiratory viruses were associated with clinical signs and bronchoalveolar lavage fluid (BALF) cytology in Thoroughbred racehorses. Equine herpesviruses (EHV-1, 2, 4, 5) and equine rhinitis A and B viruses (ERBV, ERAV) genomes were quantified by qPCR in nasopharyngeal, tracheal, and BALF samples collected after racing. The relationships between virus detection and load and clinical signs, performance, BALF cytology, and environmental exposures were examined with generalized linear mixed models. Ninety-two samples were collected from 31 horses. EHV-1 and ERAV were not found; EHV-4 was detected in only one sample. EHV-2, EHV-5 and ERBV were more likely to be detected in upper airway samples than in BALF (P < 0.0001). Neither respiratory virus detection nor load was associated with clinical signs or performance. Nasopharyngeal detection and load of ERBV and tracheal detection and load of EHV-5 were associated with increased proportions of neutrophils in BALF (P < 0.003). However, nasopharyngeal detection and load of EHV-5 was not (P = 0.11). Nasopharyngeal detection and load of EHV-2 were associated with decreased BALF mast cell proportions. Respirable dust exposures were significantly higher in horses with detection of ERBV when compared to horses with no detectable ERBV (P < 0.001). Our results suggest that ERBV, EHV-2 and EHV-5 are commonly present in upper airways of healthy racehorses; however, the role they play in the etiology of equine asthma remains unclear.

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.

Molecular characterization of the apoptosis-related SH3RF1 and SH3RF2 genes and their association with exercise performance in Arabian horses

BACKGROUND

Apoptosis plays an important role in the regulation of healthy tissue growth and development as well as in controlling the maintenance of homeostasis in exercising muscles. During an intensive physical effort, the regulation of cell death by apoptosis results in the replacement of unaccustomed muscle cells by new cells that are better suited to exercise. The aim of this study was to determine the expression of two genes (SH3FR1 and SH3RF2) that control apoptosis in muscle tissues during training periods characterized by different intensities. The gene expression levels were estimated using real-time PCR method in skeletal muscle biopsies collected from 15 Arabian horses (untrained, after an intense gallop phase, and at the end of the racing season). An association study was performed on 250 Arabian horses to assess the effect of the SH3RF2:c.796 T > C (p.Ser266Pro) variant on race performance traits in flat gallop-racing.

RESULTS

A gene expression analysis confirmed a significant decrease (p < 0.01) in the anti-apoptotic SH3RF2 (POSHER) gene during training periods that differed in intensity. The highest SH3RF2 expression level was detected in the muscles of untrained horses, whereas the lowest expression was identified at the end of the racing season in horses that were fully adapted to the exercise. A non-significant decrease in SH3RF1 gene expression following the training periods was observed. Moreover, a serine substitution by proline at amino acid position 266 (CC genotype) was negatively associated with the probability of winning races, the number of races in which a horse occurred and the financial value of the prizes. Horses with the TT genotype achieved the highest financial benefits, both for total winnings and for winnings per race in which the horses participated.

CONCLUSIONS

The present study showed the supposed regulation mechanism of exercise-induced apoptosis in horses at the molecular level. The identified SH3RF2: c.796 T > C missense variant was associated with selected racing performance traits, which is important information during the evaluation of horses' exercise predisposition. The association results and frequencies of the CT and TT genotypes suggest the possibility of using SH3RF2 variant in selection to improve the racing performance of Arabian horses.

Integrated analysis of microRNA and mRNA expressions in peripheral blood leukocytes of Warmblood horses before and after exercise

Exercise capacity is a valuable trait in horses, and it has been used as a horse selection criterion. Although exercise affects molecular homeostasis and adaptation in horses, the mechanisms underlying these effects are not fully described. This study was carried out to identify changes in the blood profiles of microRNAs (miRNAs) and mRNAs induced by exercise in horse leukocytes. Total RNAs isolated from the peripheral blood leukocytes of four Warmblood horses before and after exercise were subjected to next-generation sequencing (NGS) and microarray analyses to determine the miRNA and mRNA expression profiles, respectively. The expressions of 6 miRNAs, including 4 known and 2 novel miRNAs, were altered by exercise. The predicted target genes of the differentially expressed miRNAs identified by NGS were matched to the exercise-induced mRNAs determined by microarray analysis. Five genes (LOC100050849, LOC100054517, KHDRBS3, LOC100053996, and LOC100062720) from the microarray analysis were matched to the predicted target genes of the 6 miRNAs. The subset of mRNAs and miRNAs affected by exercise in peripheral blood leukocytes may be useful in elucidating the molecular mechanisms of exercise-associated physiology in horses.

Endurance Exercise Ability in the Horse: A Trait with Complex Polygenic Determinism

Endurance horses are able to run at more than 20 km/h for 160 km (in bouts of 30-40 km). This level of performance is based on intense aerobic metabolism, effective body heat dissipation and the ability to endure painful exercise. The known heritabilities of endurance performance and exercise-related physiological traits in Arabian horses suggest that adaptation to extreme endurance exercise is influenced by genetic factors. The objective of the present genome-wide association study (GWAS) was to identify single nucleotide polymorphisms (SNPs) related to endurance racing performance in 597 Arabian horses. The performance traits studied were the total race distance, average race speed and finishing status (qualified, eliminated or retired). We used three mixed models that included a fixed allele or genotype effect and a random, polygenic effect. Quantile-quantile plots were acceptable, and the regression coefficients for actual vs. expected log10p-values ranged from 0.865 to 1.055. The GWAS revealed five significant quantitative trait loci (QTL) corresponding to 6 SNPs on chromosomes 6, 1, 7, 16, and 29 (two SNPs) with corrected p-values from 1.7 × 10-6 to 1.8 × 10-5. Annotation of these 5 QTL revealed two genes: sortilin-related VPS10-domain-containing receptor 3 (SORCS3) on chromosome 1 is involved in protein trafficking, and solute carrier family 39 member 12 (SLC39A12) on chromosome 29 is active in zinc transport and cell homeostasis. These two coding genes could be involved in neuronal tissues (CNS). The other QTL on chromosomes 6, 7, and 16 may be involved in the regulation of the gene expression through non-coding RNAs, CpG islands and transcription factor binding sites. On chromosome 6, a new candidate equine long non-coding RNA (KCNQ1OT1 ortholog: opposite antisense transcript 1 of potassium voltage-gated channel subfamily Q member 1 gene) was predicted in silico and validated by RT-qPCR in primary cultures of equine myoblasts and fibroblasts. This lncRNA could be one element of the cardiac rhythm regulation. Our GWAS revealed that equine performance during endurance races is a complex polygenic trait, and is partially governed by at least 5 QTL: two coding genes involved in neuronal tissues and three other loci with many regulatory functions such as slowing down heart rate.

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.

Prospection of genomic regions divergently selected in racing line of Quarter Horses in relation to cutting line

Selection of Quarter Horses for different purposes has led to the formation of lines, including racing and cutting horses. The objective of this study was to identify genomic regions divergently selected in racing line of Quarter Horses in relation to cutting line applying relative extended haplotype homozygosity (REHH) analysis, an extension of extended haplotype homozygosity (EHH) analysis, and the fixation index (F ST) statistic. A total of 188 horses of both sexes, born between 1985 and 2009 and registered at the Brazilian Association of Quarter Horse Breeders, including 120 of the racing line and 68 of the cutting line, were genotyped using single nucleotide polymorphism arrays. On the basis of 27 genomic regions identified as selection signatures by REHH and F ST statistics, functional annotations of genes were made in order to identify those that could have been important during formation of the racing line and that could be used subsequently for the development of selection tools. Genes involved in muscle growth (n=8), skeletal growth (n=10), muscle energy metabolism (n=15), cardiovascular system (n=14) and nervous system (n=23) were identified, including the FKTN, INSR, GYS1, CLCN1, MYLK, SYK, ANG, CNTFR and HTR2B.