- Research article
- Open Access
Evaluation of microRNA expression in plasma and skeletal muscle of thoroughbred racehorses in training
BMC Veterinary Researchvolume 13, Article number: 347 (2017)
Circulating miRNAs (ci-miRNAs) are endogenous, non-coding RNAs emerging as potential diagnostic biomarkers. Equine miRNAs have been previously identified including subsets of tissue-specific miRNAs. In order to investigate ci-miRNAs as diagnostic tools, normal patterns of expression for different scenarios including responses to exercise need to be identified. Human studies have demonstrated that many ci-miRNAs are up-regulated following exercise with changes in expression patterns in skeletal muscle. However, technical challenges such as haemolysis impact on accurate plasma ci-miRNA quantification, with haemolysis often occurring naturally in horses following moderate-to-intense exercise. The objectives of this study were to identify plasma ci-miRNA profiles and skeletal muscle miRNAs before and after exercise in Thoroughbreds (Tb), and to evaluate for the presence and effect of haemolysis on plasma ci-miRNA determination. Resting and post-exercise plasma ci-miRNA profiles and haemolysis were evaluated in twenty 3 year-old Tbs in sprint training. Resting and post-exercise skeletal muscle miRNA abundance was evaluated in a second cohort of eleven 2 year-old Tbs just entering sprint training. Haemolysis was further quantified in resting blood samples from twelve Tbs in sprint training. A human plasma panel containing 179 miRNAs was used for profiling, with haemolysis assessed spectrophotometrically. Data was analysed using a paired Student’s t-test and Pearson’s rank correlation.
Plasma ci-miRNA data for 13/20 horses and all skeletal muscle miRNA data passed quality control. From plasma, 52/179 miRNAs were detected at both time-points. Haemolysis levels were greater than the threshold for accurate quantification of ci-miRNAs in 18/25 resting and all post-exercise plasma samples. Positive correlations (P < 0.05) between haemolysis and miRNA abundance were detected for all but 4 miRNAs, so exercise-induced changes in plasma ci-miRNA expression could not be quantified. In skeletal muscle samples, 97/179 miRNAs were detected with 5 miRNAs (miR-21-5p, let-7d-3p, let-7d-5p, miR-30b-5p, miR-30e-5p) differentially expressed (DE, P < 0.05) between time-points.
The degree of haemolysis needs to be determined prior to quantifying plasma ci-miRNA expression from horses in high-intensity exercise training. Identification of DE miRNAs in skeletal muscle indicates modification of miRNA expression may contribute to adaptive training responses in Tbs. Using a human plasma panel likely limited detection of equine-specific miRNAs.
MicroRNAs (miRNAs) are endogenous, non-coding RNAs with complex roles causing translational repression/degradation of bound mRNA , regulating many cellular processes . They have intracellular effects, are extremely stable when secreted into the bloodstream  and have been identified in multiple body fluids . MiRNAs are often tissue-specific and can be secreted from their tissue of origin into the bloodstream bound to lipoproteins  or packaged into extracellular microparticles . When in blood they are termed circulating miRNAs (ci-miRNAs)  and are emerging as potential diagnostic biomarkers due to their stability  and ease of access for sample collection.
Human studies have identified a multitude of ci-miRNAs up-regulated following intense exercise training [9,10,11,12]. Furthermore, exercise training has been found to alter the expression pattern of miRNAs in skeletal muscle [11, 13, 14]. It has been hypothesised that skeletal muscle may contribute to changes in differentially expressed (DE) ci-miRNAs in response to exercise training by releasing exosomes into plasma [15, 16], supporting the possibility that DE plasma ci-miRNAs could be used to assess adaptive training responses of skeletal muscle.
MiRNA analysis in horses is an emerging area of research. To-date, >700 equine miRNAs have been identified with subsets of tissue-specific DE miRNAs isolated  including skeletal muscle and blood [18, 19]. Recently, 197 miRNAs were identified in equine skeletal muscle, with 76 found to be muscle-specific . However, to-date there has only been one study evaluating ci-MRNA in exercising horses . In this report, 167 DE miRNAs were identified from pre- and post-endurance exercise blood samples; however, whole-blood rather than plasma/serum was evaluated making it difficult to determine if ci-miRNAs were expressed from red blood cells (RBC) or other tissue sites such as skeletal muscle.
There are technical challenges involved in plasma/serum ci-miRNA measurement such as inaccurately quantifying plasma RNA due to low yields, haemolysis impacting on ci-miRNA quantification and a lack of validated stable reference miRNAs [22,23,24,25]. Since haemolysis may occur naturally in horses following moderate-to-intense exercise [26, 27], this needs to be considered when attempting to evaluate equine plasma/serum ci-miRNA expression.
The aim was to identify DE plasma ci-miRNAs and skeletal muscle miRNAs before and after exercise in Thoroughbreds (Tb), and to evaluate for the presence and effect of haemolysis on plasma ci-miRNA determination.
A Department of Health License (B100/3525), ethical approval from the University College Dublin Animal Research Ethics Committee and owner consent were obtained for all horses and procedures in this study.
All horses in the study were from the same stable and were managed and trained in a similar way . Horses trained 6 days per week with gradual introduction of 800–1000 m sprint training (work day [WD]) alternating with submaximal training, following which horses entered competitive racing.
Cohort A: Plasma ci-miRNA
Plasma lactate concentrations ([LA]) and ci-miRNA were evaluated in n = 20 sex- and fitness-matched 3 year-old Tbs (n = 10 males, n = 10 females) in active training and racing before (T0) and 5 mins after (T5min) a WD training session. The degree of haemolysis was evaluated in all samples.
Cohort B: Skeletal muscle miRNA
Resting (T0) and T5min plasma [LA] along with T0 and 4 h (T4hr) post-exercise skeletal muscle miRNAs were evaluated in n = 11 2 year-old Tbs (n = 7 males, n = 4 females) before and after their first WD training session.
Cohort C: Degree of plasma haemolysis
The degree of haemolysis was further evaluated in T0 plasma samples from an additional n = 12 Tbs in active training and racing (n = 5 males, n = 7 females, 2–4 years old).
Resting jugular venous blood samples were collected before feeding and exercise. Samples for plasma [LA] and ci-miRNA measurement were placed into fluoride oxalate and ethylenediamine tetraacetic acid (EDTA) tubes, respectively. Once collected all blood tubes were immediately centrifuged with plasma separated. [LA] was measured on-site with an autoanalyser.Footnote 3 Plasma for haemoglobin concentrations [Hb] and ci-miRNA analysis was stored at 4 °C until transported to be stored at -80 °C (within 6 h of collection) and analysed as a batch at a later date.
Blood samples for miRNA
MiRNAs were extracted from plasma samplesFootnote 4 with 4 μL amounts then reversely transcribed in 20 μL reactions.Footnote 5 Exogenous RNA spike-ins (UniSp2, UniSp4, UniSp5, cel-miR-39-3p) provided by the manufacturer of the human serum/plasma focus miRNA PCR panel kit were added during RNA extraction and analysed according to the manufacturer’s instructions5.
The degree of haemolysis was assessed by quantifying plasma [Hb] using a UV/VIS spectrophotometer,Footnote 6 measuring oxyhaemoglobin absorbance at k = 414 nm . A value of ≤0.2 absorption units (AU) was used as the threshold for acceptable haemolysis and accurate quantification of ci-miRNA detection .
Skeletal muscle biopsy samples for miRNA
Skeletal muscle biopsies were obtained from the gluteus medius using a 6 mm-diameter, modified Bergstrom biopsy needleFootnote 7 according to previously described methods . Samples were preserved in RNAlaterFootnote 8 until further processing.
Total RNA was extracted from 50 to 100 mg of tissue using a protocol combining TRIzol reagent,Footnote 9 DNase treatmentFootnote 10 and RNeasy.Footnote 11 RNA was quantified using a Qubit fluorimeter with the Qubit RNA HS assay. RNA quality was assessed using the 18S/28S ratio and RNA integrity number on an Agilent Bioanalyzer with the RNA 6000 Nano Lab-Chip kit.Footnote 12 RNA (2 μL@5 ng/ul) was reversely transcribed in 10 μL reactions5. UniSp6 (exogenous RNA spike-in) was added during cDNA conversion.
Real-time quantitative polymerase chain reaction measurements
cDNA was assayed in 10 μL PCR reactions5. All miRNAs were assayed once by RT-qPCR on the Exiqon Human Serum/Plasma Focus miRNA PCR panel I and II. Each panel consisted of a 96 well plate with 179 miRNAs of interest and 15 controls (inter-plate calibrators, negative and spike-in controls). Amplification was performed in an Applied Biosystems StepOnePlus™ Real-Time PCR System. Amplification curves were analysed using Applied Biosystems software to obtain raw PCR cycle threshold (C t) values.
Quality control measures (inter-plate calibration, ommission of outliers, miRNAs with C t values >34 from analysis) were carried out using GeneEx software. MiRNAs with >80% missing data were removed from the analysis, with missing values replaced with the mean value for that miRNA. Results from muscle biopsies were normalised to the global mean.
A paired Student’s t-test was used to compare T0 and T5min plasma haemolysis and ci-miRNAs and T0 and T4h muscle miRNA abundance. A Shapiro-Wilk test was performed to confirm that the C t data was normally distributed and the Benjamini and Hochberg method used to correct for multiple testing. Pearson’s rank correlation was used to evaluate correlation between T0 and T5min ci-miRNA expression level and absorbance at 414 nm, and a paired Student’s t-test used to assess the significance of that correlation. Significance was set at P ≤ 0.05.
Prediction and annotation of miRNA targets
Prediction of potential target genes and pathway analysis was completed using the DIANA-mirPath software  in the Diana Tools suite. The DIANA-microT-CDS algorithm was used to predict potential target genes for DE miRNAs; these were then utilised by miRPath for KEGG pathway analysis, enriched through analysis of all predicted targets. A 0.8 default setting was used as the MicroT prediction score threshold with species set to human since there was no equine option. Since there is a high conservation of miRNAs accross species suggesting the sharing of biological functions, it was believed that inferences could be made on the function of equine miRNAs based on their human homologues. Significance for over-represented KEGG pathways was set as P ≤ 0.05 following Benjamini and Hochberg correction for multiple testing.
Quality control for miRNA quantification
The miRNA data for 7/20 horses from cohort A failed quality control and were omitted from analysis. All samples from cohort B passed quality control.
Physiological exercise parameters
For all horses undertaking an exercise test in the study, the average Vpeak was 16.9 ± 0.98 m/s (range 15.2–18.5), the average HRpeak 221.4 ± 7.1 bpm (range 209–232) and the average total exercise distance 1130.19 ± 145.1 m/s (range 876–1471). The average pre- and 5 mins post-exercise plasma [LA] were 0.46 ± 0.09 mmol/l (range 0.33–0.63) and 26.3 ± 3.68 mmol/l (range 19.2–32.1), respectively. There were no significant differences between males and females regarding any of the exercise and physiological measurements.
Cohort A: Plasma ci-miRNA
In the equine plasma samples, 29% (52/179) of the miRNAs on the human plasma panel were detected at both timepoints, all of which had increased expression (3.23 ± 1.5 fold mean increase, 1.1–10.3 range increase) at T5min when compared to T0. However, absorbance >0.2 AU at 414 nm were detected in 11/13 T0 and all T5min plasma samples for horses in cohort A, with absorbance peaks significantly higher in the T5min samples (P = 1.2−05; Fig. 1). The T0 absorbance ranged from 0.15–0.39 AU with an average of 0.25 ± 0.06 AU while the T5min values ranged from 0.26–0.76 AU with an average of 0.46 ± 0.13 AU. The expression levels of 35 miRNAs were significantly correlated (r > 0.55, P < 0.05) with absorbance and thus haemolysis levels in both T0 and T5min samples while 13 miRNA expression levels were only significantly correlated with haemolysis in T0 samples; 4 miRNAs were not correlated with haemolysis at either time point. Due to the high levels of haemolysis in the plasma samples and the strong positive correlation between haemolysis and miRNA abundance (Fig. 2), it was not possible to accurately quantify exercise-induced changes in plasma ci-miRNA expression in the present study. There were no significant differences between males and females regarding any of the haemolysis and ci-miRNA measurements.
Cohort B: Skeletal muscle miRNA
In the equine skeletal muscle biopsy samples, 55% (98/179) of the miRNAs on the human plasma panel were detected at both T0 and T4h. Following the Benjamini and Hochberg correction for multiple testing, 5 of these miRNAs were significantly DE at T4h (P < 0.05; Fig. 3). The greatest expression difference was detected in miRNA miR-21-5p which increased 0.92 ± 0.2 fold. The expression of the miR-30b-5p and miR-30e-5p family members increased 0.52 ± 0.1 and 0.64 ± 0.2 fold, respectively. Expression of let-7d-3p and let-7d-5p decreased 0.68 ± 0.1 and 0.44 ± 0.1 fold, respectively. Using the DIANA-microT-CDS algorithm a total of 1687 gene targets for these 5 miRNAs were identified. The largest number of gene targets were identified for the miR-30 family members. In total 1077 of the predicted targets were unique to these 2 miRNAs while only 3 gene targets were predicted for let-7d-3p. The most over-represented KEGG pathways for each miRNA based on a predicted target gene list are summarised in Table 1 [31, 32]. Pathways previously implicated in the equine exercise response were significantly enriched among predicted targets of the DE miRNAs (Table 2). These included the Jak-STAT signaling, MAPK signaling, insulin signaling and mTOR signaling pathways and long-term potentiation.
Cohort C: Degree of plasma haemolysis
The average absorbance at 414 nm for this cohort was 0.21 ± 0.09 AU, ranging from 0.10–0.41 AU (Fig. 1); 7/12 samples had absorbtion peaks >0.2 AU. There was no significant difference in resting haemolysis levels between cohorts A and C, confirming that haemolysis at rest is prevelant in Tb racehorses in high-intensity exercise training.
The stability and ease of access for collecting samples for ci-miRNA measurement make them attractive as potential diagnostic biomarkers. However, in order to investigate ci-miRNAs as diagnostic tools it is vital to identify normal patterns of expression for different scenarios including responses to exercise. Furthermore, rigourous quality control measures must be in place to ensure reliable results. One adverse situation identified to cause variation in plasma ci-miRNA measurement is RBC haemolysis resulting in leakage of miRNAs into the plasma [22, 23, 33]. This is important to consider since there appears to be distinct miRNA profiles for the various cell types in blood . Keeping this in mind, the use of stringent quality controls and careful sample handling has resulted in the accurate quantification of resting and post-exercise plasma ci-miRNAs from human subjects .
In the present study, haemolysis above the threshold for accurate quantification of ci-miRNAs was identified in the majority of the resting plasma samples as well as in all of the T5min plasma samples. Intravascular haemolysis has been identified to occur in horses following racing , exercise on the treadmill [34,35,36] and after work on the gallops . Exercise-induced haemolysis in horses primarily results from increased fragility of erythrocytes [37, 38], hypothesized to be due to the frequent accumulation of RBCs in the spleen [35, 37]. Changes in pH due to anaerobic exercise has also been associated with increased erythrocyte fragility , which also occurs in human althletes following maximal exercise . It was interesting to note that in the present study even resting plasma samples had a degree of haemolysis above the threshold for accurate quantification of ci-miRNA detection . We hypothesised that since the horses were in active training for 6 days of the week, which is typical for a flat racehorse training establishment, this resulted in almost daily and consistent exercise-induced insult to the RBCs resulting in persistent levels of haemolysis. Thus, the degree of plasma haemolysis and subsequent effect on miRNA abundance is important to consider when investigating plasma ci-miRNA levels in horses that are in active exercise programmes, even when evaluating pre-exercise samples.
In the present study plasma samples were analysed spectrophotometrically to determine levels of free Hb, as this is an extremely simple and accurate method to determine the extent of plasma haemolysis . For the 48/52 miRNAs detected in plasma in the current study, miRNA abundance was significantly correlated with absorbance at one or both sampling time points supporting correlation between haemolysis and plasma ci-miRNA detection. Since measurement of ≥90% of the detected plasma ci-miRNAs in the current study was correlated with haemolysis, an accurate assessment of the changes in ci-miRNA abundance in response to exercise could not be made. It is believed, however, that this problem could be overcome through the quantification of exosomal ci-miRNAs (plasma miRNAs packaged in vesicles) rather than measurement of all plasma ci-miRNAs. This technique would exclude all unpackaged miRNAs which may have been leaked into the plasma from haemolysed erythrocytes.
Our research group previously demonstrated that mRNA in equine skeletal muscle responds to a single bout of exercise as well as exercise training [31, 32]. These results underpinned our interest in evaluating the effect that exercise may have on equine skeletal muscle miRNA expression. The initial profiling for equine skeletal muscle miRNA expression in the present study was carried out on horses just entering sprint training. This was based on results from previous studies in our laboratory in which Tb skeletal muscle was found to have the greatest transcriptional response to intense exercise at this training stage . The 4 h post-exercise sampling time was also chosen based on previous work performed by our research group where it was identified to be the time-point exhibiting the most transcriptional activity . These observations are in agreement with several human exercise studies [40,41,42].
In this study we observed changes in equine skeletal muscle miRNA abundance following exercise. All five of the DE miRNAs from skeletal muscle expressed at T4 in the present study had been identified amongst the 48 miRNAs detected from haemolysed plasma samples from horses in cohort A. Furthermore, all five are known equine miRNAs with miR-30e-5p previously reported to be expresesed in muscle in a tissue-specific manner . Of great interest is that these five DE miRNAs are amongst 167 DE miRNAs reported from pre- and 30 min post-endurance exercise whole-blood samples in a group of n = 14 pure-breed/half-breed Arabian horses .
In the present study, two members of the let-7 family of miRNAs, let-7d-3p and let-7d-5p, and miR-21-5p had increased expression following exercise. The Let-7 family of miRNAs are involved in the regulation of glucose homeostasis and insulin sensitivity , both key processes in energy metabolism during exercise. Mir-21 has been shown to be involved in the fibrogenic pathway and the progression of Duchenne muscular dystrophy , so may play a role in muscle remodelling during exercise training.
Two members of the miR-30 family, miR-30b-5p and miR-30e-5p, had decreased expression post-exercise. The miR-30 family has been shown to be down-regulated in mdx4cv mice (models for Duchenne muscular dystrophy), with in vitro analysis indicating that miR-30 miRNAs are decreased following injury and are increased during myoblast differentiation . These findings suggest that these miRNAs may play an important role in skeletal muscle growth and repair . The decreased expression in response to exercise may be related to repair of exercise-induced muscle damage. The over-representation of miRNA gene targets in exercise-related KEGG pathways certainly suggests that miRNAs play a key role in the modulation of gene expression in response to exercise.
In conclusion, our results indicate even resting plasma from horses in high-intensity training may not be suitable for ci-miRNA quantification due to haemolysis, although it is possible this problem could be overcome through the quantification of exosomal ci-miRNAs. However, identification of DE miRNAs in skeletal muscle indicates that modification of miRNA expression may contribute to adaptive training responses in Tbs. It is likely that additional miRNAs are differentially regulated in response to exercise but were undetectable due to the use of a human plasma panel. Future studies using a global platform such as miRNA sequencing will be required to generate a comprehensive profile of miRNAs in equine skeletal muscle at rest and following exercise. The most powerful approach to such an investigation would be an analysis of the miRNA-Seq transcriptome.
Polar Equine S810i, Polar Electro Ltd., Warwick, UK
GPSports Systems SPI10, Canberra, Australia
YSI 2300 stat plus, YSI UK Ltd., Hampshire, UK
biofluids miRCURYTM RNA Isolation Kit, Exiqon, Denmark
miRCURY LNA™ Universal RT microRNA PCR, polyadenylation and cDNA synthesis kit, Exiqon, Denmark
NanoDrop™ Bioanalyzer ND2000, Labtech, Uckfield, UK
Jørgen KRUUSE A/S, Langeskov, Denmark
Ambion, ThermoFisher Scientific Inc., Darmstadt, Germany
Invitrogen, ThermoFisher Scientific Inc., Darmstadt, Germany
RNase-free DNase, Qiagen Ltd., Manchester, UK
Mini-kit, Qiagen Ltd., Manchester, UK.
Agilent Technologies Ireland, Ltd.
- C t :
Ethylenediamine tetraacetic acid
Messenger ribonucleic acid
Red blood cells
Reverse transcription quantitative real-time polymerase chain reaction
- T0 :
before exercise (rest)
- T4hr :
4 h after exercise
- T5min :
5 mins after exercise
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Cai Y, Xiaomin Y, Songnian H, Yu JA. Brief review on the mechanisms of miRNA regulation. Genomics Proteomics Bioinformatics. 2009;7:147–54.
Cortez MA, Calin GA. MicroRNA identification in plasma and serum: a new tool to diagnose and monitor diseases. Expert Opin Biol Ther. 2009;9:703–11.
Chen X, Ba Y, Ma L, Cai X, Yin Y, Wang K, Guo J, Zhang Y, Chen J, Guo X, Li Q, Li X, Wang W, Zhang Y, Wang J, Jiang X, Xiang Y, Xu C, Zheng P, Zhang J, Li R, Zhang H, Shang X, Gong T, Ning G, Wang J, Zen K, Zhang J, Zhang C-Y. Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008;18:997–1006.
Vickers KC, Palmisano BT, Shoucri BM, Shamburek RD, Remaley AT. MicroRNAs are transported in plasma and delivered to recipient cells by high-density lipoproteins. Nat Cell Biol. 2011;13:423–33.
Turchinovich A, Weiz L, Langheinz A, Burwinkel B. Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011;39:7223–33.
Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, Peterson A, Noteboom J, O'Briant KC, Allen A, Lin DW. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci. 2008;105:10513–8.
Reid G, Kirschner MB, van Zandwijk N. Circulating microRNAs: association with disease and potential use as biomarkers. Crit Rev Oncol Hematol. 2011;80:193–208.
Baggish AL, Hale A, Weiner RB, Lewis GD, Systrom D, Wang F, Wang TJ, Chan SY. Dynamic regulation of circulating microRNAs during acute exhaustive exercise and sustained aerobic exercise training. J Physiol. 2011;589:3983–94.
Bye A, Rosjo H, Aspenes ST, Condorelli G, Omland T, Wisloff U. Circulating microRNAs and aerobic fitness–the HUNT-study. PLoS One. 2013;8:e57496.
Russell AP, Lamon S, Boon H, Wada S, Guller I, Brown EL, Chibalin AV, Zierath JR, Snow RJ, Stepto N, Wadley GD, Akimoto T. Regulation of miRNAs in human skeletal muscle following acute endurance exercise and short-term endurance training. J Physiol. 2013;591:4637–53.
Nielsen S, Akerstrom T, Rinnov A, Yfanti C, Scheele C, Pedersen BK, Laye MJ. The miRNA plasma signature in response to acute aerobic exercise and endurance training. PLoS One. 2014;9:e87308. doi:10.1371/journal.pone.0087308.
Davidsen PK, Gallagher IJ, Hartman JW, Tarnopolsky MA, Dela F, Helge JW, Timmons JA, Phillips SM. High responders to resistance exercise training demonstrate differential regulation of skeletal muscle microRNA expression. J Appl Physiol. 2011;110:309–17.
Keller P, Vollaard NB, Gustafsson T, Gallagher IJ, Sundberg CJ, Rankinen T, Britton SL, Bouchard C, Koch LG, Timmons JA. A transcriptional map of the impact of endurance exercise training on skeletal muscle phenotype. J Appl Physiol. 2011;110:46–59.
Guescini M, Guidolin D, Vallorani L, Casadei L, Gioacchini AM, Tibollo P, Battistelli M, Falcieri E, Battistin L, Agnati LF, Stocchi V. C2C12 myoblasts release micro-vesicles containing mtDNA and proteins involved in signal transduction. Exp Cell Res. 2010;316:1977–2984.
Pedersen BK, Febbraio MA. Muscles, exercise and obesity: skeletal muscle as a secretory organ. Nat Rev Endocrinol. 2012;8:457–65.
Van der Kolk JH, Pacholewska A, Gerber V. The role of microRNAs in equine medicine: a review. Vet Q. 2015;35:88–96.
Barrey E, Bonnamy B, Barrey EJ, Mata X, Chaffaux S, Guerin G. Muscular microRNA expressions in healthy and myopathic horses suffering from polysaccharide storage myopathy or recurrent exertional rhabdomyolysis. Equine Vet J. 2010;42:303–10.
Gim JA, Ayarpadikannan S, Eo J, Kwon YJ, Choi Y, Lee HK, Park KD, Yang YM, Cho BW, Kim HS. Transcriptional expression changes of glucose metabolism genes after exercise in thoroughbred horses. Gene. 2014;547:152–8.
Kim MC, Lee SW, Ryu DY, Cui FJ, Bhak J, Kim Y. Identification and characterization of microRNAs in normal equine tissues by next generation sequencing. PLoS One. 2014;9:e93662.
Mach N, Plancade S, Pacholewska A, Lecardonnel J, Riviere J, Moroldo M, Vaiman A, Morgenthaler C, Beinat M, Nevot A, Robert C, Barrey E. Integrated mRNA and miRNA expression profiling in blood reveals candiate biomarkers associated with endurance exercise in the horse. Sci Rep. 2016;6:22932.
Kirschner MB, Kao SC, Edelman JJ, Armstrong NJ, Vallely MP, van Zandwijk N, Reid G. Haemolysis during sample preparation alters microRNA content of plasma. PLoS One. 2011;6:e24145.
McDonald JS, Milosevic D, Reddi HV, Grebe SK, Algeciras-Schimnich A. Analysis of circulating microRNA: preanalytical and analytical challenges. Clin Chem. 2011;57:833–40.
Blondal T, Jensby Nielsen S, Baker A, Andreasen D, Mouritzen P, Wrang Teilum M, Dahlsveen IK. Assessing sample and miRNA profile quality in serum and plasma or other biofluids. Methods. 2013;59:S1–6.
Kirschner MB, Edelman JJ, Kao SC, Vallely MP, van Zandwijk N, Reid G. The impact of hemolysis on cell-free microRNA biomarkers. Front Genet. 2013;4:94.
Masini AP, Tedeschi D, Badagli P, Sighieri C, Lubas C. Exercise-induced intravascular haemolysis in standardbred horses. Comp Clin Path. 2003;12:45–8.
Cywinska A, Szarska E, Kowalska A, Ostaszewski P, Schollenberger A. Gender differences in exercise--induced intravascular haemolysis during race training in thoroughbred horses. Res Vet Sci. 2011;90:133–7.
Fonseca RG, Kenny DA, Hill EW, Katz LM. The association of various speed indices to training responses in thoroughbred flat racehorses measured with a global positioning and heart rate monitoring system. Equine Vet J. 2010;42:51–7.
Dingboom EG, Dijkstra G, Enzerink E, van Oudheusden HC, Weijs WA. Postnatal muscle fibre composition of the gluteus medius muscle of Dutch Warmblood foals; maturation and the influence of exercise. Equine Vet J. 1999;31:95–100.
Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M, Paraskevopoulou MD, Prionidis K, Dalamagas T, Hatzigeorgiou AGDIANA. miRPath v.2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. 2012;40:W498–504.
McGivney BA, Eivers SS, MacHugh DE, MacLeod JN, O'Gorman GM, Park SD, Katz LM, Hill EW. Transcriptional adaptations following exercise in thoroughbred horse skeletal muscle highlights molecular mechanisms that lead to muscle hypertrophy. BMC Genomics. 2009;10:638.
McGivney BA, McGettigan PA, Browne JA, Evans AC, Fonseca RG, Loftus BJ, Lohan A, MacHugh DE, Murphy BA, Katz LM, Hill EW. Characterization of the equine skeletal muscle transcriptome identifies novel functional responses to exercise training. BMC Genomics. 2010;11:398.
Pritchard CC, Kroh E, Wood B, Arroyo JD, Dougherty KJ, Miyaji MM, Tait JF, Tewari M. Blood cell origin of circulating microRNAs: a cautionary note for cancer biomarker studies. Cancer Prev Res. 2012;5:492–7.
Schott HC, Hodgson DR, Bayly WM. Haematuria, pigmenturia and proteinuria in exercising horses. Equine Vet J. 1995;27:67–72.
Hanzawa K, Higara A, Yoshida Y, Hara H, Kai M, Kubo K, Watanabe S. Effects of exercise on plasma haptoglobin composition in control and splenectomized thoroughbred horses. J Equine Sci. 2002;13:89–92.
Inoue Y, Matsui A, Asai Y, Aoki F, Matsui T, Yano H. Effect of exercise on iron metabolism in horses. Biol Trace Elem Res. 2005;107:33–42.
Hanzawa K, Kai M, Hirag A, Watanabe S. Fragility of red cells during exercise is affected by blood pH and temperature. Equine Vet J. 1999;31:610–1.
Hanzawa K, Watanabe S. Changes in osmotic fragility of erythrocytes during exercise in athletic horses. J Equine Sci. 2000;11:51–61.
Smith JA, Telford RD, Kolbuch-Braddon M, Weidermann MJ. Lactate/H+ uptake by red blood cells during exercise alters their physical properties. Eur J Appl Physiol Occup Physiol. 1995;75:54–61.
Yang Y, Creer A, Jemiolo B, Trappe S. Time course of myogenic and metabolic gene expression in response to acute exercise in human skeletal muscle. J Appl Physiol. 2005;98:1745–52.
Mahoney DJ, Safdar A, Parise G, Melov S, Fu M, MacNeil L, Kaczor J, Payne ET, Tarnopolsky MA. Gene expression profiling in human skeletal muscle during recovery from eccentric exercise. Am J Physiol Regul Integr Comp Physiol. 2008;294:R1901–10.
Nader GA, von Walden F, Liu C, Lindvall J, Gutmann L, Pistilli EE, Gordon PM. Resistance exercise training modulates acute gene expression during human skeletal muscle hypertrophy. J Appl Physiol. 2014;116:693–702.
Frost RJ, Olson EN. Control of glucose homeostasis and insulin sensitivity by the Let-7 family of microRNAs. Proc Natl Acad Sci. 2011;108:21075–80.
Ardite E, Perdiguero E, Vidal B, Gutarra S, Serrano AL, Munoz-Canoves P. PAI-1-regulated miR-21 defines a novel age-associated fibrogenic pathway in muscular dystrophy. J Cell Biol. 2012;196:163–75.
Guess MG, Barthel KK, Harrison BC, Leinwand LA. miR-30 family microRNAs regulate myogenic differentiation and provide negative feedback on the microRNA pathway. PLoS One. 2015;10:e0118229.
The authors wish to express their appreciation to the trainer Mr. Jim Bolger and his staff.
This study has emanated from research conducted with the financial support of Science Foundation Ireland under Grant Number 11/PI/1166.
Availability of data and materials
All data generated or analysed during this study are included in this published article. Accession numbers for the miRNA data generated in this study are: MIMAT0004484 (hsa-let-7d-3p), MIMAT0000076 (hsa-miR-21-5p), MIMAT0000065 (hsa-let-7d-5p), MIMAT0000420 (hsa-miR-30b-5p) and MIMAT0000692 (hsa-miR-30e-5p).
University College Dublin Animal Research Ethics Committee approval, a Department of Health License (B100/3525) and explicit owner/trainer informed consent were obtained for the use of the horses and procedures in this study.
Consent for publication
None of the authors has any financial or personal relationships that could inappropriately influence or bias the content of the paper. E.W.H. is a shareholder in Plusvital Ltd., an equine nutrition and genetic testing company. Plusvital Ltd. has been granted a license for commercial use of data contained within patent applications: United States Provisional Serial Number 61/136553 and Irish patent application number 2008/0735, Patent Cooperation Treaty filing: A method for predicting athletic performance potential, September 7, 2009. E.W.H. and L.M.K. are named on the applications. The patent contents are not related to this manuscript. Plusvital Ltd. had no part in the research in the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.