Characterisation of equine satellite cell transcriptomic profile response to β-hydroxy-β-methylbutyrate (HMB)

β-Hydroxy-β-methylbutyrate (HMB) is a popular ergogenic aid used by human athletes and as a supplement to sport horses, because of its ability to aid muscle recovery, improve performance and body composition. Recent findings suggest that HMB may stimulate satellite cells and affect expressions of genes regulating skeletal muscle cell growth. Despite the scientific data showing benefits of HMB supplementation in horses, no previous study has explained the mechanism of action of HMB in this species. The aim of this study was to reveal the molecular background of HMB action on equine skeletal muscle by investigating the transcriptomic profile changes induced by HMB in equine satellite cells in vitro. Upon isolation from the semitendinosus muscle, equine satellite cells were cultured until the 2nd day of differentiation. Differentiating cells were incubated with HMB for 24 h. Total cellular RNA was isolated, amplified, labelled and hybridised to microarray slides. Microarray data validation was performed with real-time quantitative PCR. HMB induced differential expressions of 361 genes. Functional analysis revealed that the main biological processes influenced by HMB in equine satellite cells were related to muscle organ development, protein metabolism, energy homoeostasis and lipid metabolism. In conclusion, this study demonstrated for the first time that HMB has the potential to influence equine satellite cells by controlling global gene expression. Genes and biological processes targeted by HMB in equine satellite cells may support HMB utility in improving growth and regeneration of equine skeletal muscle; however, the overall role of HMB in horses remains equivocal and requires further proteomic, biochemical and pharmacokinetic studies.

Despite the large amount of literature linked to HMB, only two reports have supported anecdotal data showing HMB's benefits in thoroughbred racing horses. In one of them, exercising thoroughbred race horses receiving daily 15 g Ca salt of HMB during a 16-week training season showed a significant decrease in post-exercise blood creatinine phosphokinase and lactate levels over both training and racing seasons (13) . Miller et al. (14) observed similar results when supplementing racing horses with 10 g of HMB daily, with a significantly improved win rate after the 1st month of racing. Taken together, the present experiment meets the demand for more detailed studies concerning HMB's effectiveness in horses.
In adult skeletal muscle, regeneration and hypertrophy depend on the activation of mononucleated, muscle precursor cells called satellite cells (SC) (15) , embedded between the sarcolemma and the basement membrane of muscle fibres. Previous in vitro and in vivo studies indicate that HMB may activate SC (8,10,16,17) , but the mechanism underlying this action remains unclear. Some evidence suggests that HMB regulates the expression of myogenesis-related genes (8) ; however, until now, no one has demonstrated any effect of HMB on global gene expression.
The horse is a valuable animal model for studying exercise physiology. Gene expression determines most of the phenotype; therefore, the present study focused on revealing the molecular background of HMB action in equine skeletal muscle by investigating the impact of HMB on global gene expression in differentiating equine satellite cells (ESC) in vitro. To our knowledge, this is the first study where HMB's trancriptomic profile was described. This in vitro model can help identify and better understand the potential therapeutic options to promote muscle regeneration and energy metabolism in horses and other mammals.

Cell culture
Media and reagents. The following materials were used during cell culture: the Ca salt (monohydrate) of HMB (Ca-HMB) was purchased from Metabolic Technologies; Dulbecco's Modified Eagle Medium (DMEM) (1×) with glutamax, fetal bovine serum (FBS), horse serum (HS) and antibiotics (AB)penicillinstreptomycin and fungizonewere purchased from Gibco, Life Technologies; penicillium crystalicum (AB) was purchased from Polfa Tarchomin; PBS, protease from Streptomyces griseus and DMSO were purchased from Sigma Aldrich. Tissue culture flasks Primaria (25, 75 cm 2 ) and Collagen I Cellware six-well plates were purchased from Becton Dickinson. Ca-HMB was transformed to the acid form by acidification with 1 N-HCl. HMB was then extracted four times with diethyl ether. The pooled organic layer was dried under vacuum for 24 h at 38°C. The resulting free acid was 99 % HMB as assessed by HPLC.
Muscle sampling and satellite cells isolation. Semitendinosus muscle samples were collected ex vivo from six horses (6-month-old, healthy colts). Muscle sampling and ESC isolation are described in detail by Szcześniak et al. (18) . In brief, semitendinosus muscle samples were dissected free of surrounding tissues, sliced, washed in PBS with decreasing antibiotics concentration, suspended in FBS with 10 % DMSO, cooled to −80°C and stored in liquid N 2 . Before isolation, the samples were thawed, centrifuged and washed three times with PBS along with antibiotics. Samples were incubated with DMEM/ AB/protease from S. griseus and sieved in order to separate tissue debris. The filtrates were centrifuged three times, re-suspended in proliferation medium (10 %FBS/10 %HS/DMEM/AB) and transferred to polypropylene Petri culture disks. One-and-a-half hours of preplating was performed to minimise possible fibroblast contamination. Subsequently, the supernatant containing ESC was transferred to Primaria culture flasks.
Cell culture and experimental design. The experimental design is presented in Fig. 1. Upon isolation, samples of ESC (n 6) were incubated for 10 d in Primaria culture flasks. The proliferation medium was changed every 2 d. On the 10th day, cells were trypsinised, and 30 000 cells (counted by Scepter Cell Counter; Merck Millipore) from each flask were transferred to the respective wells of two six-well plates. One plate was dedicated to HMB treatment and one served as the control. After obtaining 80 % of confluence, the proliferation medium was replaced with a differentiation medium (2 % HS/DMEM/ AB). Immediately after 48 h of differentiation, the medium from one plate was replaced by a differentiation medium containing 50 µM of HMB, whereas in the second plate the standard differentiation medium was used as a control. After 24 h, the medium from each plate was discarded, plates were washed with PBS and stored at −80°C until further analysis. The concentration of HMB was based on the available literature values and cell viability colourimetric assay test with 3-(4,5-dimethylthiazol-2-yl)-2-5-diphenyltetrazolium bromide (data not shown).
Microarray analysis and real-time quantitative PCR validation RNA isolation, validation, labelling hybridisation and microarray analysis. Total RNA from HMB and control cells was isolated according to the protocol supplied with the miR-Neasy Mini Kit (Qiagen). RNA quantity was measured spectrophotometrically using NanoDrop (NanoDrop Technologies). The analysis of final RNA quality and integrity was performed with BioAnalyzer 2100 (Agilent Technologies). To ensure optimal microarray data quality, only samples with the highest RNA integrity number (RIN) ≥ 9·2 were included in the analysis.
Analysis of gene expression profiles was performed using Horse Gene Expression Microarray, 4× 44K (Agilent Technologies). Low Input Quick Amp Labeling Kit (Agilent Technologies) was used to amplify and label total RNA (100 ng) to generate complementary RNA (cRNA). On each two-colour microarray, 825 ng of cRNA from HMB-exposed cells (labelled by Cy5, n 4) and 825 ng of cRNA from control cells (labelled by Cy3, n 4) were hybridised to the arrays (Gene Expression Hybridization Kit; Agilent Technologies) according to the manufacturer's protocol.
RNA Spike-In Kit (Agilent Technologies) was used as an internal control to efficiently monitor microarray workflow for

Statistical analysis
Statistical analysis was performed using Gene Spring 13 software (Agilent Technologies) with the default setting for two-colour microarrays. The estimated significance level (P value) was corrected for multiple hypotheses testing using the Benjamini and Hochberg false discovery rate (FDR) adjustment. mRNA with FDR ≤ 0·05 were selected as significantly differentially expressed genes (DEG). The microarray experiment was performed according to Minimum information about a microarray experiment (MIAME) guidelines (19) . The data discussed in this publication have been deposited in National Center for Biotechnology Information's (NCBI's) Gene Expression Omnibus (GEO) (20) and are accessible through GEO Series accession number GSE74495 (http:// www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE74495).
Complementary DNA synthesis and real-time quantitative PCR. To independently assess expression changes for a selected group of genes obtained from the microarray data, the real-time quantitative PCR (RT-qPCR) method was applied. The sequences of verified genes, complementary to those on microarrays, were obtained from Ensembl database. Primers were designed using Primer-Blast software (NCBI database) and then checked for secondary structures using the Oligo Calculator (http://www.basic.northwestern.edu/biotools/oligocalc.html). The secondary structures of the amplicon were examined using m-fold Web Server (http://mfold.rna.albany. edu/?q=mfold). The sequences of primers are listed in Table 1. The primers were purchased from Oligo IBB (Polish Academy of Science). Each primer pair was quality tested to ensure that a single product was amplified (dissociation curve analysis) and that there was no primer-dimer coupling.
A quantity of 1 µg of total RNA from HMB-treated and control cells (n 6) was reverse transcribed using a Transcription First Strand cDNA Synthesis Kit (Agilent Technologies). All analyses were performed on individual samples of total RNA using a SensiFAST SYBR lo-ROX Kit (Blirt, Bioline) following the manufacturer's protocol. Assays for each gene were conducted in duplicate in a Stratagene Mx3005p thermal cycler (Agilent Technologies) according to the following protocol: pre-incubation for 2 min at 95°C and amplification (forty cycles), with denaturation at 95°C for 5 s and annealing at the temperatures specified in Table 1 for 15 s. The dissociation curve setting was as follows: denaturation at 95°C for 0 s, annealing (at the temperatures specified in Table 1), continuous melting up to 95°C for 0 s (slope = 0·1°C/s) and cooling at 40°C for 30 s. Glyceraldehyde 3-phosphate dehydrogenase (Gapdh) was used as a reference gene. The relative expression of the target gene was calculated according to the following formula: where ΔC T is the difference in C T between the targeted gene and the reference control. Results were calculated as 2 ÀΔΔCT using GenEx 6.0 (MultiD Analyses) (21) . The amplification efficiency (E = 10 (−1/slope) -1) was determined using a comparative quantitation standard curve and was >0·9 for each target gene and the reference gene. Standard curves were generated using a fourpoint 1:10 dilution series starting with cDNA representing 10 ng of input total RNA. RT-qPCR analysis was conducted according to a standardised approach (22) .

Functional analysis
The list of DEG was examined by the Functional Analysis tool in the Database for Annotation, Visualization and Integrated Discovery (DAVID version 6.7) to assign them to gene ontology (GO) terms and KEGG pathways (Kyoto Encyclopedia of Genes and Genomes) (23) . Human background was used for this analysis, because far more human genes are annotated and more information in databases is available for humans than for horses. Enrichment of DEG was calculated by EASE score (modified Fisher exact test). For further analysis and visualisation of data, the Pathway Studio Web Mammalian was used. This database of functional relationships between mammalian proteins is compiled using Med Scan technology from over twenty-four million PubMed abstracts and over 3·5 million Elsevier full-text papers. All identified relations were filtered by reference count (≥2) to ensure maximal confidence levels, which means that the number of publications confirming each relationship was ≥2.

Microarray analysis
Analysis of gene expression between HMB-treated and control cells revealed statistically significant (FDR ≤ 0·05) differences in the case of 627 records. Within them were 361 unduplicated, identified transcript ID including 159 up-and 202 down-  Fig. 1. Experiment design. Equine satellite cells (ESC) were cultured until they reached 80 % confluence; next, the proliferation medium was replaced with a differentiation medium. After the 2nd day of differentiation, cells were incubated for 24 h with β-hydroxy-β-methylbutyrate (HMB). Following the HMB treatment, differentiating cells were scraped and stored at −80°C until further analysis.
Satellite cell response to β-hydroxy-β-methylbutyrate regulated DEG, in the HMB v. the control group. All array data are plotted and shown in the online Supplementary Material S1. Table 2 presents genes selected for discussion, presumably involved in HMB action on ESC.

Functional analysis
DAVID functional analysis assigned DEG to seventy-five biological processes (BP), eleven cellular components and ten molecular functions as well as four KEGG pathways (EASE score P < 0·05). All GO considered significant are shown in the online Supplementary Material S2. KEGG pathways and the most significantly enriched (EASE score <0·01) GO retrieved from DAVID are presented in Table 3, providing a comprehensive overview of important processes, most likely induced by HMB in differentiating ESC.
Using Pathway Studio Web Mammalin Build Pathway Wizard Find Direct Links, we depicted all genes discussed in the present study that can directly or indirectly affect skeletal muscle cell functions (Fig. 3). Moreover, Pathway Studio Web Mammalian Build Pathway Wizard Find Common Targets algorithm allowed us to identify cell processes regulated by at least two of the DEG according to literature data. This resulted in fifty-six identified targets; among these, the twelve regulated by the highest number of genes were considered to be the most important for the HMB effect on ESC. A chart presenting these processes is presented in Fig. 4. From all targeted cell processes, we selected the most important relationships and are presented in Fig. 5. The online Supplementary Material S3  CCCGCAGAGTTGACACAATA  TGTGGCATCGTACAAAGCAT  60  282  2  Myf5  GGAGACGCCTGAAGAAAGTC  CCGGCAGGCTGTAGTAATTC  60  171  3  Rbfox  GAACCAGGAGGGATCTTCCA  TTGCCATACACAGGCTCTTG  60  213  4  S1pp1  CCCAAGTCAGTCCAACGAAA  GGCACAGCTGGTGTAAAAAC  60  143  5  Tgfb2  AGTACTACGCCAAGGAGGTT  TAGGCGGGATGGCATTTTCC  60  72  6  Trim63  AAGGAGGCAGCCAGGTAGAG  CACGGACACTGAGCCACTTC  62  220  7 Gapdh GTTTGTGATGGGCGTGAACC GTCTTCTGGGTGGCAGTGAT 60 198 Cfl2, coffilin 2; Myf5, myogenic factor 5; Rbfox, RNA binding protein, fox-1 homolog C. elegans; S1pp1, secreted phosphoprotein 1; Tgfb2, transforming growth factor, β2; Trim63, muscle-specific RING finger protein 1; Gapdh, glyceraldehyde 3-phosphate dehydrogenase. Table 2. List of selected differentially expressed genes in β-hydroxy-β-methylbutyrate-treated v. control equine satellite cells (false discovery rate ≤0·05, n 4) No. Gene symbol Fold change Description False discovery rate (corrected p-value) contains details of all identified relationships between DEG and cell processes.

Discussion
The objective of the present study was to identify the molecular background of HMB action on equine skeletal muscle. In order to cover all the salient points of functional analysis, only relations significant in DAVID and possessing the highest reference number in Pathway Studio analysis were considered to be important. To date, no official genome nomenclature has been established for the horse. According to the guidelines published by The International Society for Animal Genetics, for all genes with human orthologues, official human gene symbols (Human Genome Organisation (HUGO) Gene Nomenclature Committee) are applied. We decided to use a primary SC model because of its stem cell potential. SC are able to differentiate into multiple mesenchymal lineages (24) and to self-renew (25) , because of which they maintain extraordinary regenerative properties of skeletal muscles. However, the capacity of SC to proliferate and differentiate may vary depending on the origin of the muscle (26) , cell surface markers expression (27) , myogenic regulatory factors (MRF) expression (28) and muscle fibre type (29) . In our study, all samples of ESC were isolated from semitendinosus muscle, which in horses is composed mainly of type II fast-twitch fibre muscle (30) . SC originating from this type of muscle may have less adipogenic properties compared with SC from type I fibres (29) . Heterogeneity of the SC could limit in vivo significance of the data obtained in the present study.
In general, the present analysis underlined the role of HMB as a global regulator, which is shown by the strong over-representation of genes linked to the BP: 'regulation of developmental process' and 'positive regulation of BP'. Moreover, functional analysis revealed significant enrichment in ontology terms associated with cellular responses ( Table 3). The three main cellular processes include cell proliferation, apoptosis and differentiation, which suggest that HMB is an important cell growth regulator (Fig. 4 and 5).
In adult skeletal muscle, extracellular matrix proteins anchor SC between the basal lamina and the apical sarcolemma, which create a specialised micro-environment called a stem cell niche. It is able to produce factors controlling stem cell behaviour (31) . Impaired adhesion of SC to their niche can stimulate proliferation (32) . Thereby, enrichment of the terms 'regulation of cell adhesion' and 'cellular localisation' may suggest HMB's ability to indirectly control ESC proliferation by affecting their localisation in the niche.

Muscle development
The term 'muscle organ development' is the most significantly enriched annotation among genes regulated in ESC exposed to HMB (Table 3). This indicates that at least at the mRNA level HMB may affect muscle development (summarised on Fig. 3). A total of fourteen DEG were annotated to this term; however, among them, Mapk14 (mitogen-activated protein kinase 14) possessed the highest potential to regulate other genes and cell processes ( Fig. 3 and 5). Mapk14 is activated by extracellular stimuli such as pro-inflammatory cytokines or physical stress, leading to direct activation of multiple cellular processes such as proliferation, differentiation, apoptosis and transcription regulation (33) . In SC, phosphorylation of MAPK14 may induce initiation (34,35) or withdrawal (36) from the cell cycle. The second can lead either to terminal differentiation or to programmed cell death (37) depending on the nature of the stimulant and cell type. In vitro studies suggest that the two isoforms of Mapk14, p38α and p38β, appear to have different effects on cardiomyocyte hypertrophy: p38β seems to be more potent in inducing hypertrophy, whereas p38α appears to be more important in apoptosis (38) . The contribution of Mapk14 in cellular responses to HMB has already been reported by Kornasio et al. (8) , who suggested that the MAPK/ERK pathways mediate HMB's effects on myoblast proliferation. HMB-related increase in phosphorylation of MAPK14 was also observed in dexamethasone-induced muscle atrophy in rats (39) .
Another gene of particular importance to the 'muscle organ development' term is Myf5, belonging to the MRF family of transcription regulators (46) . The high expression of Myf5 in adult skeletal muscle features committed SC and decreases when differentiation to myotubes occurs (46,47) . Accordingly, decreased expression levels of Myf5 in ESC at the beginning of  differentiation may indicate that HMB enhanced withdrawal of equine myoblasts from the cell cycle, compared with control cells. This finding is accompanied by previous reports presenting an HMB-dependent increase in mRNA and protein levels of muscle differentiation markers such as MyoD and myogenin (8,16) . However, at the time of our analysis, none of the differentiation markers reached significance criteria in ESC, which may emphasise the need for time-course studies in the future. Another down-regulated gene in HMB-treated cells was Tgf-β2. Activity of Tgf-β2 has been recently linked with increased proliferation and delayed differentiation in C2C12 (48) ; thus, its down-regulation may confirm HMB-mediated enhancement of differentiation in ESC.

Muscle protein metabolism
One of the first described mechanisms of HMB action was the effect on muscle protein metabolism. Preliminary studies suggest that HMB protects the skeletal muscle by inhibiting protein degradation (5) and by stimulating protein synthesis (6) ; however, this issue is subjected to constant research (17) . Functional analyses have demonstrated significant DEG enrichment of terms associated with cellular protein maintenance (Table 3, Fig. 4). The three most important genes of this group are Cul3 (cullin 3), Trim63 and Mapk14 (Fig. 5). Cul3 is a scaffold protein of E3 ubiquitin-protein ligase complexes, which mediate the ubiquitination and subsequent proteasomal degradation of target proteins. Cul3 also interacts with Kelch family proteins, and disturbances in functioning of this complex are implicated in muscle myopathies (51) . E3 Ubiquitin ligase produced by Trim63 regulates the proteasomal degradation of muscle proteins and inhibits de novo skeletal muscle protein synthesis under amino acid starvation, consequently leading to muscle atrophy (52) . As observed in the present study, down-expression of Trim63 mediated by HMB confirms the results obtained by Aversa et al. (39) in a dexamethasone-induced muscle atrophy  (17,53) . This indicates that the effect of HMB on this gene expression could be species and/or condition related. Multiple studies suggest that Mapk14 signalling may be involved in HMB-mediated stimulation of protein synthesis in catabolic conditions (8,39,54) , which may be confirmed by the up-regulation of this gene in HMB-treated ESC.

Lipid metabolism and energy homoeostasis
Recent studies have revealed that HMB supplementation may alter metabolism, as evidenced by improved aerobic performance and increased fat loss during exercise (11,12) . This is confirmed in our study, which showed influence of DEG on cell processes such as 'energy homoeostasis', 'lipid metabolism', 'glucose import', 'fatty acid oxidation' and 'gluconeogenesis' (Fig. 4 and 5). An extensive amount of research describing the positive role of Mapk14 on glucose uptake (55) and gluconeogenesis (56) has been published. Thereby, we postulate that apart from the established role of Mapk14 in HMB-dependent influence on protein metabolism and cell growth it can mediate HMB influence on energy homoeostasis as well. The rate of post-exercise muscle glycogen synthesis is 2-3-fold slower in horses compared with other mammals (1) ; therefore, the positive impact of HMB on glucose uptake could enhance this process in equine skeletal muscles. This is an interesting aspect of our study, which deserves more attention in future investigations. Another salient point of HMB influence on metabolism may be the transcription factor Esrra (oestrogen-related receptor α), controlling vast gene networks involved in all aspects of energy homoeostasis, including lipid and glucose metabolism as well as mitochondrial biogenesis and function (57) . Common targets algorithm showed its strong association with 'fatty acid oxidation' and 'lipid metabolism' (Fig. 5). Essra is targeted by Ppargc1b (peroxisome proliferator-activated receptor γ, coactivator 1 β) (PPAR-γ coactivator), a well-established regulator of β-oxidation of fatty acids and oxidative phosphorylation in mitochondria, which is highly induced during myogenic differentiation (58) . Prkab2 (protein kinase, AMP-activated, β2 noncatalytic subunit) is essential for the regulation of a multitude of metabolic processes maintaining energy homoeostasis, especially in tissues with high metabolic rates, such as skeletal muscle (59) . Bruckbauer et al. (12) reported that HMB increases the activity of Prkab2 in adipocytes and muscle cells; however, our results showed that HMB slightly decreased its expression in ESC at the time of the analysis. Prkab2 senses cellular energy levels. In response to low cellular ATP levels, Prkab2 switches off ATP-consuming anabolic pathways (mechanistic target of rapamycin (mTOR) kinase pathway), which results in inhibition of cell growth, proliferation and macromolecules synthesis, and at the same time Prkab2 switches on catabolic pathways that generate ATP (e.g. glucose uptake, glycolysis, fatty acid oxidation) (59) . In regulation of the cellular process 'lipid metabolism', two genes appear to take the lead -Abca1 (ATP-binding cassette, sub-family A, member 1), encoding a membrane-associated protein belonging to the ATP-binding cassette transporters superfamily and Abhd5 (abhydrolase domain-containing protein 5). The analysis indicated up-regulation of both in ESC. The latter encodes a co-activator of adipose triglyceride lipase, thereby enhancing adipocyte and muscle lipolysis (60) . Abca1 is a key regulator of the reverse cholesterol transport process and HDL biogenesis. Increased Abca1 expression was demonstrated in skeletal and cardiac muscles in response to training (61) , which indicates the role of Abca1 in the reduction of CVD risk by physical exercise.
Several reports have established HMB's role in supporting muscle cell membrane integrity during exercise (13,14) . However, as already mentioned, our analysis showed that at least at mRNA levels HMB decreased the expressions of genes encoding sarcolemmal scaffold proteins (Dmd, Lama2, Lama5). Alternatively, functional analysis enrichment of terms associated with lipid maintenance, as well as KEGG pathways 'biosynthesis of unsaturated fatty acids' and 'glicosphingolipids biosynthesis', may indicate HMB's ability to support cell membrane integrity by decreasing its rigidity (62) . Moreover, this may have an indirect impact on the inflammatory processes, signal transduction and myoblast differentiation (62,63) (Fig. 3).

Conclusions
The results presented in this study suggest the capability of HMB to influence ESC proliferation, differentiation and apoptosis as well as inflammatory response, protein anabolism, sarcolemma integrity, and cell energy utilisation and storage. As we have summarised in Fig. 5, most of the above-mentioned processes could be controlled by the Mapk14 gene, which suggests that at least at the mRNA level HMB triggers its cellular responses by stress signalling pathways. It should be noted that in vivo response of ESC to HMB may differ from the presented results because of the heterogeneity of the SC population and undefined postprandial HMB concentrations in equine skeletal muscle. Moreover, transcription is only one step in the regulatory pathway that leads to functional protein synthesis, therefore, further research on the proteomic, biochemical and pharmacodynamic level is highly recommended.
In conclusion, this study demonstrated for the first time that HMB has the potential to influence ESC by controlling its global gene expression. Transcriptomic profile analysis identified valuable gene targets of HMB in ESC, which may support the role of HMB in improving skeletal muscle growth and regeneration in horses; however, the overall role of HMB in equine skeletal muscle remains equivocal and requires further research.
Food', a decision of the Ministry of Science and Higher Education, no. 05-1/KNOW2/2015. K. A. S. carried out muscle sampling, RT-qPCR validation of microarray results, ontological analysis, interpretation of the obtained data and wrote the manuscript. A. C. carried out equine satellite cell isolation and culture analysis, RNA isolation and microarray analysis. P. O. participated in the study design and helped in manuscript revision. T. S. participated in the study design, supervised the project, performed muscle sampling and statistical analysis of microarray and RT-qPCR data, as well as assisted in the manuscript revision. All the authors read and approved the final manuscript.
The authors declare that they have no conflicts of interest.