- Research
- Open access
- Published:
Transcriptomic profile adaptations following exposure of equine satellite cells to nutriactive phytochemical gamma-oryzanol
Genes & Nutrition volume 11, Article number: 5 (2016)
Abstract
Background
Adult skeletal muscle myogenesis depends on the activation of satellite cells that have the potential to differentiate into new fibers. Gamma-oryzanol (GO), a commercially available nutriactive phytochemical, has gained global interest on account of its muscle-building and regenerating effects. Here, we investigated GO for its potential influence on myogenesis, using equine satellite cell culture model, since the horse is a unique animal, bred and exercised for competitive sport. To our knowledge, this is the first report where the global gene expression in cultured equine satellite cells has been described.
Methods
Equine satellite cells were isolated from semitendinosus muscle and cultured until the second day of differentiation. Differentiating cells were incubated with GO for the next 24 h. Subsequently, total RNA from GO-treated and control cells was isolated, amplified, labeled, and hybridized to two-color Horse Gene Expression Microarray slides. Quantitative PCR was used for the validation of microarray data.
Results
Our results revealed 58 genes with changed expression in GO-treated vs. control cells. Analysis of expression changes suggests that various processes are reinforced by GO in differentiating equine satellite cells, including inhibition of myoblast differentiation, increased proliferation and differentiation, stress response, and increased myogenic lineage commitment.
Conclusions
The present study may confirm putative muscle-enhancing abilities of GO; however, the collective role of GO in skeletal myogenesis remains equivocal. The diversity of these changes is likely due to heterogenous growth rate of cells in primary culture. Genes identified in our study, modulated by the presence of GO, may become potential targets of future research investigating impact of this supplement in skeletal muscle on proteomic and biochemical level.
Background
Since horses were domesticated on the Eurasian Steppe over 6000 years ago, they had been selected for strength, speed, and endurance exercise [23]. Nowadays in Europe, working horses are very seldom used for farming purposes; they became exclusively domestic animals exercised for competitive sport. In both equine and human athletes, nutritional ergogenic aids have been used to improve physical ability with an appropriate training regimen. Although training increases physical condition, the ease of taking a nutritional additive to improve training results explains the demand for supplementation, which may increase mechanical energy of work, delay onset of fatigue, and improve neuromuscular coordination [26].
Gamma-oryzanol (GO) is a unique mixture of triterpene alcohol and sterol ferulates present in rice bran oil, a by-product of rice processing. GO has a wide range of potential therapeutically useful activities including prevention of coronary atherosclerosis and anticancerogenic and antidiabetic action [28]; however, in our study, we focused on the anabolic properties of GO.
GO has been shown to be very safe with no major side effects being reported in either animal or human studies and is a popular ergogenic aid approved for commercial use in the USA [17]. Moreover, GO is often used by veterinary practitioners in order to improve muscle mass gain and recovery in horses. It was temporarily accounted as a prohibited substance by the Federation Equestre Internationale (FEI) unit governing rules of fair play in equestrian competition (FEI Veterinary Regulations 11th edition, 1st January 2009) [44].
Despite this, only a few studies of ergogenic use of GO were described in the peer-reviewed literature. In one study, authors looked on the resistance-weight-trained male athletes supplemented with 500 mg of GO daily [19]. However, this study failed to demonstrate any effect of GO on training performance. In a recent study, Eslami et al. [17] revealed that 600 mg/day of GO supplementation changed muscular strength in young healthy males in the 9-week resistance training without any significant alteration in anthropometric measurements. Ostaszewski et al. [49] showed that GO supplementation prevented exercise-induced muscle damage in thoroughbred race horses during 16 weeks of training. Animals received 3 g of GO a day. However, no paper describing GO’s influence on equine satellite cells was published.
Skeletal muscle satellite cells (SC) are mononucleated cells located under the basal lamina of myofiber. They play an indispensable role during adult skeletal muscle regeneration and hypertrophy. SC are quiescent during most of their normal adult life. Upon activation, they proliferate, differentiate, and finally fuse with adjacent muscle fiber or with other satellite cells to form new multinucleated fibers [40]. The procedure of harvesting satellite cells from equine muscles was first described in 1992 [22], but until now, except for that single report, studies concerning equine satellite cells in vitro are scarce. Previous transcriptomic analyses in horses were conducted on peripheral blood cells and muscle samples obtained by muscular biopsy; thus, our report is the first trial where the global gene expression in cultured equine satellite cells has been described [51]. Satellite cells require an activation signal to proliferate and differentiate and can be controlled by a variety of hormones, growth factors, and metabolites. In vitro reports suggest that satellite cells may also be regulated directly by specific dietary components including ergogenic aids [18, 50, 62]. Determining the factors that independently regulate satellite cells is important to understand overall mechanisms involved in muscle growth and regeneration.
This study focused on revealing possible molecular mechanism underlying GO’s influence on skeletal myoblasts. For this purpose, transcriptomic profile analysis of primary culture of equine satellite cells incubated with GO was performed. This in vitro model can help to identify and better understand potential therapeutic to promote muscle regeneration in mammals, among them also in sport horses.
Methods
Cell culture media and reagents
The following reagents/materials were used during cell culture: GO powder was purchased from TCI Chemicals, USA; Penicillinum crystalicum was purchased from Polfa Tarchomin, Poland; and Dulbecco’s Modified Eagle’s medium (DMEM) (1×) with Glutamax, fetal bovine serum (FBS), horse serum (HS), and antibiotics (AB), penicillin-streptomycin and amphotericin B, were purchased from Gibco, Life Technologies, USA. Phosphate-buffered saline (PBS), protease from Streptomyces griseus, and DMSO were purchased from Sigma-Aldrich, USA. Primaria tissue culture flasks (25 and 75 cm2) and Collagen I Cellware six-well plates were purchased from Becton Dickinson, USA.
Animals and muscle samples
Samples of semitendinosus muscle were collected from six horses (6-month-old, healthy colts), aseptically in a commercial abattoir during routine slaughter. Next, they were dissected free of connective and adipose tissues to minimize contamination with adipocytes and fibroblasts. The tissue was sliced to pieces, washed four times in PBS with gradually decreasing antibiotic concentration (40,000 units (U) Penicillinum crystalicum/100 ml PBS, 20,000 U Penicillinum crystalicum/100 ml PBS), suspended in sterile FBS with addition of 10 % DMSO, gradually cooled down to −80 °C, and then stored in liquid nitrogen until isolation.
Satellite cell isolation, cell culture, and experimental design
Since the purpose of the study was to evaluate the impact of GO on in vitro differentiation, primary satellite cell cultures from semitendinosus muscle of all horses were isolated and assessed for cell viability (MTT assay) and fusion index (data not shown). For further analysis, the cell line with the best scores was chosen.
To isolate equine satellite cells, six samples of semitendinosus muscle were thawed in water bath in 37 °C, centrifuged, and washed three times with PBS with Penicillinum crystalicum. Then, each sample was incubated (1.5 h) with DMEM/AB/protease from S. griseus (pH 7.3) and sieved in order to separate tissue debris. The filtrates were centrifuged three times, resuspended in proliferation medium (10 % FBS/10 % HS/DMEM/AB), and transferred to polypropylene Petri culture disks. To minimize possible fibroblast contamination, 1.5-h pre-plating was used. Subsequently, supernatant containing myoblasts was transferred to Primaria culture flasks. The growth medium was changed every 2 days. On the 10th day of proliferation, cells were trypsinized and counted by Scepter Cell Counter (Merck Millipore, Germany). Equal number of cells (30,000) from each isolation was transferred to separate the well of Collagen I Cellware six-well plates. After obtaining 80 % of confluence, proliferation medium was replaced by differentiation medium (2 % HS/DMEM/AB). Figure 1 presents equine SC at various stages of growth. At the 48th hour of differentiation, 0.125 μM GO was added and cells were incubated for 24 h. Because GO is insoluble in water, DMSO (0.04 μl/ml) was used as a vehicle. The concentration of GO has been chosen on the basis of assessing cell viability with 3-(4,5-dimethylthiazol-yl)-2-5-diphenyltetrazolium bromide (MTT) colorimetric assay (data not shown). Control medium contained DMSO in the same dose as in GO-treated cells. Following the GO treatment, medium from each plate was discarded and plates were stored at −80 °C until further analysis. Cell culture’s schema is illustrated in Fig. 2.
RNA isolation, validation, labeling, and hybridization
Total RNA (from GO-treated and control groups, n = 6) was isolated according to the protocol supplied with the miRNeasy Mini Kit (Qiagen, USA). RNA quantity was measured spectrophotometrically using NanoDrop (NanoDrop Technologies, USA). The analysis of the final RNA quality and integrity was performed with BioAnalyzer 2100 (Agilent Technologies, USA). To ensure optimal microarray data quality, four RNA samples with the highest RIN ≥ 9.2 were included into the analysis.
Analysis of gene expression profile was performed using Horse Gene Expression Microarray, 4 × 44 K (n = 4) (Agilent Technologies, USA). Low Input Quick Amp Labeling Kit (Agilent Technologies, USA) was used to amplify and label total RNA (100 ng) to generate complementary RNA (cRNA). On each two-color microarray, 825 ng of cRNA from GO-exposed cells (labeled by Cy5, n = 4) and 825 ng of cRNA from control cells (labeled by Cy3, n = 4) were hybridized to the arrays (Gene Expression Hybridization Kit, Agilent Technologies, USA) according to the manufacturer’s protocol.
RNA Spike-In Kit (Agilent Technologies, USA) was used as an internal control to efficiently monitor microarray workflow for linearity, sensitivity, and accuracy. Acquisition and analysis of hybridization intensities were performed using the Agilent DNA microarray scanner and Agilent Feature Extraction software 10.7.3.1 according to the standard manufacturer’s procedures and Linear Lowess normalization (Cy5/Cy3 dye bias compensation).
The statistical analysis was performed using Gene Spring 13 software (Agilent Technologies, USA) with default setting for two-color microarrays. The estimated significance level (p value) was corrected for multiple hypotheses testing using Benjamini and Hochberg false discovery rate (FDR) adjustment. Messenger RNAs (mRNAs) with FDR ≤ 0.05 and fold change (FC) ≥1.3 were selected as significantly differentially expressed.
The microarray experiment was performed according to MIAME guidelines [6]. The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus [16] and are accessible through GEO Series accession number GSE73730 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE73730).
cDNA synthesis and quantitative real-time PCR
To verify microarray results, the quantitative real-time polymerase chain reaction (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 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 the primers are listed in Table 1. The primers were purchased from Oligo IBB (Polish Academy of Science, Warsaw, Poland). 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.
One microgram of total RNA from GO-treated and control samples (n = 6) was reverse transcribed using Trancription First-Strand cDNA Synthesis Kit (Roche, USA). All analyses were performed on individual samples of total RNA using a SensiFAST SYBR lo-ROX Kit (Blirt, Bioline, Poland) following the manufacturer’s protocol. Assays for each gene were conducted in duplicate in a Stratagene Mx3005p thermal cycler according to the following protocol: pre-incubation for 2 min at 95 °C and amplification (40 cycles) containing denaturation at 95 °C for 5 s and annealing at a temperature specified in Table 1 for 15 s. Dissociation curve setting was as follows: denaturation at 95 °C for 0 s, annealing at a temperature 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. Gapdh was used as a reference gene. The relative expression of the target gene was calculated according to the following formula:
where ΔCT is the difference in CT between the targeted gene and reference control. Results were calculated as 2−ΔΔCT [38] using GenEx 6.0 (MultiD Analyses, Sweden). The amplification efficiency (E = 10(−1/slope) − 1) was determined by performing a comparative quantitation standard curve and was >0.9 for each target gene and the reference gene. Standard curves were generated using a four-point 1:10 dilution series starting with cDNA representing 10 ng of input total RNA. qPCR analysis has been conducted according to a standardized approach [10].
Functional analysis
The list of significantly modulated by GO genes was analyzed by Functional Analysis tool in the Database for Annotation, Visualization and Integrated Discovery (DAVID version 6.7) (da Huang et al. 2009) in order to identify their annotations in gene ontology divided into three functional classes: biological processes, cellular components, and molecular functions. Gene ontology enrichment was calculated by EASE score corrected for multiple hypotheses testing using Benjamini and Hochberg FDR.
Relationships between all differentially expressed genes (microarray and RT-PCR) were visualized with Pathway Studio’s Build Pathway functionality (Elsevier, USA) which is based on the wave-propagation algorithm developed for the navigation through complex networks. In this analysis, Find Direct Links/All Objects Directions Algorithm was used. Protein entity type and direct regulation, regulation, expression, and promoter-binding relationships were selected in the analysis. Such created graph was expanded with five muscle-related biological processes.
Results
The number of differentially expressed genes
Comparison of gene expression between GO-treated and control cells revealed statistically significant (FDR ≤ 0.05, FC ≥ 1.3) differences in the expression of 97 transcripts. Sequences of differentially expressed probes not annotated by microarray manufacturer were compared with the NCBI nucleotide databases (using blastn) to assign a gene name. Probes not matching in 100 % with any eukaryotic mRNA/gene were excluded from gene lists used for functional analysis. That resulted in 58 identified, unduplicated, transcript IDs including 17 up- and 41 down-regulated genes in GO vs. control group. The list of identified transcripts can be found in Table 2. All array data are plotted and shown in Additional file 1.
In this paper, symbols of genes and their protein products have been presented in their original form referring to specificity of the cited paper, as lowercase italic or uppercase letters, when they refer to transcriptomic or proteomic research, respectively. In the following discussion, arrows indicate up- (↑) or down-regulation (↓) of gene expression. Validated (from microarray) or tested (mstn, myog, igf1) by qPCR genes are marked with hash (#).
Real-time qPCR
According to the ontological classification and available literature, seven genes involved in the skeletal muscle development, sarcomere development, actin binding, and regulation of cell proliferation and differentiation were selected for qPCR validation. Except for gap junction protein alpha 1 (gja1), all expression changes from qPCR overlapped microarray results (Table 3, Fig. 3). In our study, we focused on the impact of GO on muscle development; thus, we also measured mRNA levels of three genes from muscle development profile myog, igf1, and mstn, which did not reach the FC ≥ 1.3 threshold (Fig. 4).
Functional analysis
Significantly enriched gene ontologies retrieved from DAVID output files are presented in Additional file 2. In Table 4, selection of most significant biological processes is shown (EASE score <0.001), providing a comprehensive overview of important processes possibly induced by GO in differentiating ESC.
Using Pathway Studio’s Build Pathway algorithm, we identified 116 relationships within selected entities, which are displayed in a graph view on Fig. 4. Seventy of identified relationships had the highest (third) confidence level which means that the number of references confirming this relationship was ≥3. Detected types of regulation were regulation (86), expression (28), direct regulation (1), and promoter binding (1). Fifty-five relations were positive, 22 negative, and 39 not defined. Chosen relationships are further discussed in the following chapter.
Discussion
According to Fernyhough et al. [18], oral supplements may affect myofibers in five potential ways: (1) by direct interaction, (2) by direct influence and producing paracrine agents that affect SC, (3) by affecting SC, (4) by producing hormones that exert indirect effects in muscle cells, and (5) by protecting all cell types against oxidative damage. Although, the exact mode of GO action on muscles was not revealed, GO has been previously postulated to act through mechanisms 4 [12, 44] and 5 [17, 49]. In further discussion, we suggest that the direct interaction of GO with SC (3) may be a part of its skeletal muscle-enhancing mechanism [63]. For this purpose, genes important for myogenesis have been selected and are discussed in four functional groups.
Ontology of differentially expressed genes
According to functional analysis, the five most significant biological processes induced by 24-h incubation of differentiating equine SC with 0.125 μmol GO are related to skeletal muscles; however, muscle organ development is the only record showing significance after FDR correction (Table 4). Moreover, GO influenced genes associated with contractile apparatus development (Additional file 2; CC and MF sheet). All this together indicates that GO, indeed, may directly influence SC activity and in vitro differentiation. Other significantly enriched processes referred to macromolecule metabolism and positive regulation of transcription (Table 4). The last has also a confirmation in enriched molecular functions: transcription regulator activity and transcription cofactor activity (Additional file 2; MF sheet). Functional analysis revealed also that GO influenced growth factor activity (Additional file 2; MF sheet). Axon and neuron projection-related genes (Additional file 2; CC sheet) could be connected with recent in vivo findings by Eslami et al. [17] who suggested that GO-induced increase in muscle strength may arise from neural adaptations without altering muscle mass.
Interaction of GO with myogenic regulatory factors
Pre- and postnatal skeletal muscle development depends on the expression of four basic helix-loop-helix transcription factors called myogenic regulatory factors (MRFs). This group shares the ability to induce myogenic differentiation when expressed in non-muscle cells [43]. Three MRFs, myf5, myogenin (myog), and myod, were investigated in our study, but only two of them possessed significant changes in expression.
In adult skeletal muscle, the highest expression of myf5 is present in committed satellite cells and decreases when progression to myocytes and myotubes occurs [4]. Likewise, in cultured equine SC, qPCR showed that myf5 mRNA levels were highest during intensive proliferation and decreased dramatically at the onset of differentiation [25]. Compared to control cells, 24-h incubation of equine SC with GO induced down-regulation of myf5 (↓#). According to this observation, we suggest that GO could potentially advance equine SC differentiation.
It was shown that up-regulation of MYF5 triggers the expression of myogenin (myog). That contributes to the withdrawal of myoblasts from the cell cycle [55]. Undetectable during proliferation, Myog expression level is highest when myocytes differentiate and fuse into myotubes [25] and remains at a slightly lower level in adult skeletal muscle cells [4]. Disregarded of this, additional qPCR showed down-regulation of myog (↓#) mRNA in experimental cells. Down-regulation of both MRFs may indicate that GO actually inhibits myogenic differentiation. Our experiment was performed, 48 h after changing proliferation medium to differentiation medium. At this stage, proliferating myocytes begin to differentiate and fuse. However, the exact cell cycle of experimental cells remains unequal and undefined. Thus, it is possible that GO has multiple ways of action and by down-regulation of myog, it delays terminal differentiation and simulates proliferation of equine SC. Moreover, transcription factors involved in myogenic lineage progression are not strictly acting in a linear manner but are organized in complex feedback and feed-forward networks with other muscle development genes [4]. These relationships will be discussed further in this paper.
Genes belonging to the transforming growth factor β superfamily
Transforming growth factor β2 (tgf-β2) (↓#) belongs to a subset of transforming growth factor β superfamily members, which exert strong control over proliferation, migration, and adhesion of satellite cells [35]. TGF-β2 delays myoblast differentiation while increasing cellular proliferation. Molecular mechanism of TGF-β2 includes an increased rate of degradation of myod (Schabort et al. 2009). However, it has not been examined whether a similar relationship exists in equine skeletal muscle. In a recent study, de Mello et al. [41] showed tgf-β2 mRNA levels were high only in proliferative C2C12 cells and began to decrease after differentiation induction. Tgf-β2 mRNA levels were 10-fold lower after 3 days of differentiation that strongly agrees with our observations, which indicate that GO down-regulated tgf-β2 (↓#) in equine myoblast during the third day of differentiation, when compared to control cells.
This result encouraged us to investigate another member of transforming growth factor superfamily, playing a key role in the negative regulation of muscle growth in mammals: myostatin (mstn). Although microarray analysis did not passed FDR correction test, additional qPCR demonstrated decreased mstn (↓#) mRNA level in GO vs. control cells. The deletion of myostatin in mice induces dramatic and widespread increase in skeletal muscle mass due to both muscle hypertrophy and hyperplasia [32], and it also cause double-muscling phenotype of some cattle and sheep breeds [57]. Wicik et al. [67] showed inhibitory effect of exogenous MSTN on differentiating C2C12 myoblasts.
We suggest that decreased expression of these two genes could contribute to the putative GO-mediated muscle growth enhancement, by increasing differentiation process. However, this needs to be reconciled with decreased expression of differentiation marker myog.
Other genes involved in cell proliferation and differentiation processes
Our interest was also focused on genes involved in cell proliferation and differentiation, described in the literature as important for myogenesis.
Mitogen-activated protein 14 (mapk14) (↑) also called p-38α belongs to one of the major regulators of gene transcription and metabolism in response to oxidative, energetic, and mechanical stress in skeletal muscle. Chronic activation of MAPK14 signaling pathway has been implicated in the development of adaptive and maladaptive response in skeletal muscle pathologies, such as diabetes [7] and protein catabolism [53], and physiological states, such as growth and differentiation [34]. Mapk14 expression coincides with myog expression during differentiation [8], which could potentially indicate the role of GO in muscle differentiation; however, in our study, expression of this genes was opposite.
According to the qPCR analysis, exposition of SC to GO induced significant down-regulation of insulin-like growth factor 1 (igf-1) (↓#). Igf-1 is, up to date, the only known growth factor positively regulating proliferation and differentiation of these cells [60]. IGF-1 can indirectly affect satellite cells, by increasing the expression of myf5, myog, and myod, and embryonic myosin heavy chain [70] and decreasing expression of mstn [56]. Analysis of the relationships between the abovementioned genes places igf-1 on the top of the signaling cascade hypothetically induced by GO (Fig. 5). Thus, it may be suggested that GO decreases expression of MRF through the inhibition of igf-1 signaling. These findings are in contrast with presumed abilities of GO.
Haddad and Adams [24] indicated that intact signaling through MAPK is necessary for the development of IGF-I-induced skeletal muscle hypertrophy in vivo. However, in our study, this relationship was not observed, as indicated by different changes in the expression of these two molecules. In contrast, Trendelenburg et al. [66] showed that pro-inflammatory cytokines may inhibit myoblast growth via p38 signaling activity. This observation corresponds with up-regulation of interleukin 1 receptor type II (il1r2) (↑) expression and MRF’s down-regulation in our culture cells, suggesting that GO can activate stress-related signaling in cultured cells.
Another gene deserving attention is gja1 (↓#) gene encoding gap junction protein alpha 1—component of intercellular channels that provides a route for the diffusion of low molecular weight particles from cell to cell. Except for that, GJA1 affects cell proliferation and differentiation in various tissues. Although not present in normal adult skeletal muscle, GJA1 is necessary during regenerative myogenesis in rodents [1, 21]. The relationship between gja1 and MRF expression was observed in the L6 myogenic cell line; pharmacological gap junction blockers were found to inhibit myogenin expression and myotube formation, reversibly [3]. Likewise, delayed expression of myog was found in myoblast cultures prepared from mice with an inducible deletion of gja1 [1]. Similar function in differentiating myoblasts was reported for B cell translocation gene 1 (btg1) (↓#). The btg1 overexpression inhibited proliferation of several cell lines including NIH3T3, PC12, and QM7 quail myoblast lines [52]. Therefore, it was initially believed that these proteins possess antiproliferative properties. However, other studies demonstrated that btg1 up-regulation is crucial at the onset of myoblast differentiation [54]. More recently, Busson et al. [9] described the molecular basis of btg1 action and proved that it stimulates the activity of not only MRFs but also nuclear receptors already known as positive myogenic regulators (T3 and all-trans retinoic acid receptors).
From all these studies, it becomes evident that the outcome of gja1 and btg1 expression in muscle cells is associated with the onset of differentiation while irreversible withdrawal from the cell cycle occurs. Thus, the context of the signal is essential for understanding their role in GO signaling. Taking into account that GO was added to differentiating myoblasts, gja1 and btg1 down-regulation indicate that GO may decrease equine myoblast fusion and could suggest pro-proliferative action of btg1 and gja1 in equine SC exposed to GO, which agrees with the reduced levels of MRF’s expression.
Herein, we demonstrate down-regulation of brain-derived neurotrophic factor (bdnf) (↓#). Bdnf is a member of the mammalian family of neurotrophins, studied for their roles in the nervous system function [13]. Further expression-profiling studies have shown that bdnf is differentially expressed in skeletal muscle under various physiological and pathological conditions [11] but bdnf expression pattern in certain physiological states remained elusive. Mousavi and Jasmin [45] showed that bdnf expression during myogenic differentiation has been significantly reduced. Moreover, siRNA knockdown of bdnf expression evoked precocious differentiation of myoblasts. This effect was inhibited by the administration of recombinant bdnf. However, more recently, Clow and Jasmin [13] indicated that instead of its inhibitory role, bdnf−/− mice demonstrated impaired induction of several molecular markers of regeneration, including pax7, myod, and myog when proliferating myoblasts began to differentiate and fuse to repair damaged myofibers. The role of bdnf remains questionable; however, attenuated MRF’s expression in our culture may confirm that expression of this gene is required for normal myogenic differentiation. This finding can also be applied to the next specified molecule down-regulated by GO: αβ-crystallin (cryab), (↓) a small heat shock protein involved in preventing protein aggregation. Expression of cryab was induced selectively in myoblasts during an early stage of their differentiation in vitro [33]. Moreover, mice lacking αβ-crystallin died prematurely with extensive muscle wastage [59] which indicates its essential role during myogenesis. Indeed, the αβ-crystallin promoter contains canonical skeletal muscle-specific E-box element that binds MRFs [20]. On the other hand, overexpression of cryab inhibits myogenic differentiation by negatively regulating differentiation-related apoptosis, inhibiting caspase-3 activation and altering myod levels [33, 59].
Next, gene affected by the presence of GO which deserves attention is nuclear receptor 2F2 (nr2f2) (↓#) alias chicken ovalbumin upstream promoter transcription factor II (coup-tfII). It encodes nuclear orphan receptor of the steroid-thyroid hormone receptor superfamily. Recently, growing evidence suggests the pivotal role of COUP-TFs in mesenchymal stem cell (MSC) regulation and development [69]. As illustrated both in vitro and in vivo, loss of COUP-TFII reprograms the precursors by shifting the cell identity to osteogenic and myogenic lineages at the expense of adipogenic and chondrogenic programs, supporting the notion that COUP-TFII is essential for fate choice and commitment [68]. The molecular mechanism by which COUP-TFII affects MSC development is, at least partially, through the down-regulation of Wnt cascade substrate: Wnt10b. This phenomenon resulted in reduced fat tissue and enhanced glucose tolerance and insulin sensitivity, as well as increased bone mineralization and muscle mass in COUP-TFII-deficient mice [68]. Moreover, previous reports indicated that COUP-TFII inhibits myogenesis by the transcriptional and posttranscriptional regulation of myod [2, 46]. In conclusion, down-regulation of this gene may indicate that GO stimulates in vitro myogenesis.
Finally, in our study, increased expression (↑) of activated leukocyte cell adhesion molecule (alcam), also known as cluster of differentiation 166 (CD166), has been detected. Alcam plays an important role in cell proliferation and the differentiation of mesenchymal tissues in multiple species [61]. Moreover, alcam may increase cell survival in several tumor types as well as monocytes [61]. These abilities may result from the antiapoptotic capacity of this protein [39]. Other authors speculate that ALCAM-ALCAM interactions, in cooperation with other adhesion molecules, facilitate stroma progenitor cell contact, resulting in proliferation of progenitor cell subsets [47]. It may be presumed that expression of alcam could mediate GO action on SC, but the functional roles in muscle physiology of this marker remains unknown and needs further investigation.
Other genes positively regulating muscle differentiation, whose expression was changed by GO, include neogenin (Neo1) (↓), actin-binding Rho-activating protein (Abra) (↓), actin γ-2 (actg2) (↓), integrin beta-1-binding protein 2 (melusin 2, Itgb1bp2) (↓), muscle-specific tyrosine kinase receptor (Musk) (↓), nexilin (nexn) (↓), histone H3-K4 methyltransferase (setd7) (↓), and calpain 6 (capn6) (↓)-myogenesis-suppressing gene [5, 14, 27, 36, 37, 42, 64, 65]. Down-regulation of these genes may suggest that 24-h incubation with 0.125 μM of GO decreased differentiation of equine satellite cells.
Other gamma-oryzanol treatment-related genes
GO has been reported to possess strong antioxidant activity; thus, it is postulated that GO might enhance endurance and muscle-building capacity by hindering the production of free radicals, which in theory could lessen muscle exhaustion and fatigue in reaction to anaerobic exercise [17]. Consistent with this, significantly lower post-exercise total antioxidant status and thiobarbituric acid-reactive substance level were observed in racing horses receiving GO compared to horses from other groups [49]. GO’s influence on oxidative stress-related genes has been shown in other tissues [30, 31]. However, in our study, gene ontology analysis did not show any significant involvement of differentially expressed genes in antioxidative processes. It will be important in future studies to determine if antioxidant activity is involved in the protective effects of GO under oxidative stress conditions.
Lowered expression of stearoyl-CoA desaturase and isoform 5 gene (scd5) (↓) as well as of coup-tfII (↓) may reflect positive influence of GO on energy homeostasis in skeletal muscle. As mentioned before, Coup-tfII-deficient mice demonstrated reduced fat tissue and enhanced glucose tolerance and insulin sensitivity [68]. Likewise, SCD deficiency is evidenced to activate metabolic pathways which promote β-oxidation and decrease lipogenesis both in the liver and skeletal muscles. SCD mutation results also in general changes in the expression of genes involved in lipid metabolism. SCD1-deficient mice have increased energy expenditure and reduced body adiposity and are resistant to diet-induced obesity (Dobrzyn and Ntambi 2005). These results may confirm GO utility in the treatment of diabetes (Ohara et al. 2009) as well as positive effect on lean body mass gain.
Conclusions
In conclusion, we demonstrated for the first time that GO can affect the transcriptomic profile of equine satellite cells in vitro. Many of differentially expressed genes are reported to be crucial for skeletal muscle development. Analysis of expression changes proves that various processes are reinforced by GO in equine SC (Fig. 6); however, inhibition of myoblast differentiation is prevailing observation. Incubation with GO was performed at the onset of differentiation; thus, decreased expression of differentiation markers may indicate that GO delays this process in cultured equine SC that may contribute to increased field of cells that will later differentiate into myotubes. This could be explained in further investigation of GO’s influence on satellite cells in proliferation stage. Beyond this, 24-h incubation with GO induced different expression of genes responsible for stress response, increased myogenic lineage commitment, cell proliferation and differentiation, and finally positive influence on energy homeostasis. Our study may confirm putative muscle-enhancing abilities of GO; however, analysis of changes in gene expression did not give equivocal results. Gene transcription is only one step in the regulatory pathway that leads to the functional protein synthesis. Our results encourage for investigation of GO-skeletal muscle relation in proteomic and biochemical level in the future. The collective role of GO in skeletal myogenesis remains unclear and needs further investigation.
Compliance with ethics guidelines
This article does not contain any studies with human or animal subjects performed by the any of the authors.
References
Araya R, Eckardt D, Maxeiner S, Krüger O, Theis M, Willecke K, et al. Expression of connexins during differentiation and regeneration of skeletal muscle: functional relevance of connexin43. J Cell Sci. 2005;118(1):27–37.
Bailey P, Muscat GE, Sartorelli V, Hamamori Y. The orphan nuclear receptor, COUP-TF II, inhibits myogenesis by post-transcriptional regulation of MyoD function: COUP-TF II directly interacts with p300 and myoD. Nucleic Acids Res. 1998;26(23):5501–10.
Balogh S, Naus CC, Merrifield PA. Expression of gap junctions in cultured rat L6 cells during myogenesis. Dev Biol. 1993;155(2):351–60.
Bentzinger CF, Wang YX, Rudnicki MA. Building muscle: molecular regulation of myogenesis. Cold Spring Harbor Perspect Biol. 2012;4(2):a008342.
Brancaccio M, Guazzone S, Menini N, Sibona E, Hirsch E, De Andrea M, et al. Melusin is a new muscle-specific interactor for β1 integrin cytoplasmic domain. J Biol Chem. 1999;274(41):29282–8.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, et al. Minimum information about a microarray experiment (MIAME)—toward standards for microarray data. Nat Genet. 2001;29(4):365–71.
Brown AE, Palsgaard J, Borup R, Avery P, Gunn DA, De Meyts P, et al. p38 MAPK activation upregulates proinflammatory pathways in skeletal muscle cells from insulin-resistant type 2 diabetic patients. Am J Physiol Endocrinol Metab. 2015;308(1):E63–70.
Brown D, Hikim APS, Kovacheva EL, Sinha-Hikim I. Mouse model of testosterone-induced muscle fiber hypertrophy: involvement of p38 mitogen-activated protein kinase-mediated Notch signaling. J Endocrinol. 2009;201(1):129–39.
Busson M, Carazo A, Seyer P, Grandemange S, Casas F, Pessemesse L, et al. Coactivation of nuclear receptors and myogenic factors induces the major BTG1 influence on muscle differentiation. Oncogene. 2005;24(10):1698–710.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22.
Chevrel G, Hohlfeld R, Sendtner M. The role of neurotrophins in muscle under physiological and pathological conditions. Muscle Nerve. 2006;33(4):462–76.
Cicero AFG, Gaddi A. Rice bran oil and γ‐oryzanol in the treatment of hyperlipoproteinaemias and other conditions. Phytother Res. 2001;15(4):277–89.
Clow C, Jasmin BJ. Brain-derived neurotrophic factor regulates satellite cell differentiation and skeletal muscle regeneration. Mol Biol Cell. 2010;21(13):2182–90.
DeChiara TM, Bowen DC, Valenzuela DM, Simmons MV, Poueymirou WT, Thomas S, et al. The receptor tyrosine kinase MuSK is required for neuromuscular junction formation in vivo. Cell. 1996;85(4):501–12.
Dobrzyn A, & Ntambi JM. (2005). The role of stearoyl-CoA desaturase in the control of metabolism. Prostaglandins, leukotrienes and essential fatty acids, 73(1), 35-41.
Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.
Eslami S, Esa NM, Marandi SM, Ghasemi G, Eslami S. Effects of gamma oryzanol supplementation on anthropometric measurements & muscular strength in healthy males following chronic resistance training. Indian J Med Res. 2014;139(6):857.
Fernyhough ME, Helterline DL, Vierck JL, Dodson MV, Bucci LR, Feliciano JEFF. Myogenic satellite cell proliferative and differentiative responses to components of common oral ergogenic supplements. Res Sports Med. 2004;12(3):161–90.
Fry AC, Bonner E, Lewis DL, Johnson RL, Stone MH, Kraemer WJ. The effects of gamma-oryzanol supplementation during resistance exercise training. Int J Sport Nutr. 1997;7(4):318–29.
Gopal-Srivastava R, Haynes J, Piatigorsky J. Regulation of the murine alpha B-crystallin/small heat shock protein gene in cardiac muscle. Mol Cell Biol. 1995;15(12):7081–90.
Gorbe A, Krenacs T, Cook JE, Becker DL. Myoblast proliferation and syncytial fusion both depend on connexin43 function in transfected skeletal muscle primary cultures. Exp Cell Res. 2007;313(6):1135–48.
Greene EA, Raub RH. Procedures for harvesting satellite cells from equine skeletal muscle. J Equine Vet Sci. 1992;12(1):33–5.
Gu J, Orr N, Park SD, Katz LM, Sulimova G, MacHugh DE, et al. A genome scan for positive selection in thoroughbred horses. PLoS One. 2009;4(6):e5767.
Haddad F, Adams GR. Inhibition of MAP/ERK kinase prevents IGF-I-induced hypertrophy in rat muscles. J Appl Physiol. 2004;96(1):203–10.
Hansen, Abigail. (2014). Myostatin mRNA expression in cultured equine satellite cells. Retrieved from the University of Minnesota Digital Conservancy, http://hdl.handle.net/11299/163114
Harris PA, Harris RC. Ergogenic potential of nutritional strategies and substances in the horse. Livest Prod Sci. 2005;92(2):147–65.
Hong M, Schachter K, Jiang G, Krauss R. Neogenin regulates Shh pathway activity during digit patterning. Dev Biol. 2011;356(1):147.
Huang CCJ. Potential functionality and digestibility of oryzanol as determined using in vitro cell culture models. Doctoral dissertation, Louisiana State University. 2003: http://etd.lsu.edu/docs/available/etd-0609103-135757/.
da Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nat Protoc. 2011;4(1):44–57.
Ismail M, Al-Naqeeb G, Mamat WA, Ahmad Z. Gamma-oryzanol rich fraction regulates the expression of antioxidant and oxidative stress related genes in stressed rat’s liver. Nutr Metab. 2010;7(23):1–13.
Ismail N, Ismail M, Imam MU, Azmi NH, Fathy SF, Foo JB. Mechanistic basis for protection of differentiated SH-SY5Y cells by oryzanol-rich fraction against hydrogen peroxide-induced neurotoxicity. BMC Complement Altern Med. 2014;14(1):467.
Jank M, Zwierzchowski L, Siadkowska E, Budasz-Świderska M, Sadkowski T, Oprządek J, et al. Polymorphism in the 5’flanking region of the myostatin gene affects myostatin and TGF-β. J Anim Feed Sci. 2006;15:381–91.
Kamradt MC, Chen F, Sam S, Cryns VL. The small heat shock protein αB-crystallin negatively regulates apoptosis during myogenic differentiation by inhibiting caspase-3 activation. J Biol Chem. 2002;277(41):38731–6.
Keren A, Tamir Y, Bengal E. The p38 MAPK signaling pathway: a major regulator of skeletal muscle development. Mol Cell Endocrinol. 2006;252(1):224–30.
Kollias HD, McDermott JC. Transforming growth factor-β and myostatin signaling in skeletal muscle. J Appl Physiol. 2008;104(3):579–87.
Lamon S, Wallace MA, Léger B, Russell AP. Regulation of STARS and its downstream targets suggest a novel pathway involved in human skeletal muscle hypertrophy and atrophy. J Physiol. 2009;587(8):1795–803.
Lee A, Hakuno F, Northcott P, Pessin JE, Adcock MR. Nexilin, a cardiomyopathy-associated F-actin binding protein, binds and regulates IRS1 signaling in skeletal muscle cells. PLoS One. 2013;8(1):e55634.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8.
Ma L, Lin J, Qiao Y, Weng W, Liu W, Wang J, et al. Serum CD166: a novel hepatocellular carcinoma tumor marker. Clin Chim Acta. 2015;441:156–62.
Mauro A. Satellite cell of skeletal muscle fibers. J Biophys Biochem Cytol. 1961;9(2):493–5.
de Mello F, Streit DP, Sabin N, Gabillard JC. Dynamic expression of tgf-β2, tgf-β3 and inhibin βA during muscle growth resumption and satellite cell differentiation in rainbow trout (Oncorhynchus mykiss). Gen Comp Endocrinol. 2015;210:23–9.
Miwa T, Manabe Y, Kurokawa K, Kamada S, Kanda N, Bruns G, et al. Structure, chromosome location, and expression of the human smooth muscle (enteric type) gamma-actin gene: evolution of six human actin genes. Mol Cell Biol. 1991;11(6):3296–306.
Mok GF, Sweetman D. Many routes to the same destination: lessons from skeletal muscle development. Reproduction. 2011;141(3):301–12.
Mösseler A, Licht S, Wilhelm L, Kamphues J, Ellis AD et al. Can oral intake of gamma-oryzanol (experimentally given orally as pure substance) result in doping relevant testosterone levels in the urine of mares and geldings?. The Impact of Nutrition on the Health and Welfare of Horses, Wageningen Academic Publishers. 2010; 128;293-296.
Mousavi K, Jasmin BJ. BDNF is expressed in skeletal muscle satellite cells and inhibits myogenic differentiation. J Neurosci. 2006;26(21):5739–49.
Muscat GE, Rea S, Downes M. Identification of a regulatory function for an orphan receptor in muscle: COUP-TF II affects the expression of the myoD gene family during myogenesis. Nucleic Acids Res. 1995;23(8):1311–8.
Nelissen JM, Torensma R, Pluyter M, Adema GJ, Raymakers RA, van Kooyk Y, et al. Molecular analysis of the hematopoiesis supporting osteoblastic cell line U2-OS. Exp Hematol. 2000;28(4):422–32.
Ohara K, Uchida A, Nagasaka R, Ushio H, Ohshima T. The effects of hydroxycinnamic acid derivatives on adiponectin secretion. Phytomedicine. 2009;16(2):130–7.
Ostaszewski P, Kowalska A, Szarska E, Szpotański P, Cywinska A, Bałasińska B, et al. Effects of β-hydroxy-β-methylbutyrate and γ-oryzanol on blood biochemical markers in exercising thoroughbred race horses. J Equine Vet Sci. 2012;32(9):542–51.
Powell DJ, McFarland DC, Cowieson AJ, Muir WI, Velleman SG. The effect of nutritional status and muscle fiber type on myogenic satellite cell fate and apoptosis. Poult Sci. 2014;93(1):163–73.
Ramery E, Closset R, Art T, Bureau F, Lekeux P. Expression microarrays in equine sciences. Vet Immunol Immunopathol. 2009;127(3):197–202.
Rodier A, Marchal-Victorion S, Rochard P, Casas F, Cassar-Malek I, Rouault MJP, et al. BTG1: a triiodothyronine target involved in the myogenic influence of the hormone. Exp Cell Res. 1999;249(2):337–48.
Rom O, Kaisari S, Reznick AZ, & Aizenbud D. Peroxynitrite induces degradation of myosin heavy chain via p38 MAPK and muscle-specific E3 ubiquitin ligases in C2 skeletal myotubes. Advs. Exp. Medicine, Biology. 2014;832:1-8.
Rouault JP, Rimokh R, Tessa C, Paranhos G, Ffrench M, Duret L, et al. BTG1, a member of a new family of antiproliferative genes. EMBO J. 1992;11(4):1663.
Rudnicki MA, Schnegelsberg PN, Stead RH, Braun T, Arnold HH, Jaenisch R. MyoD or Myf-5 is required for the formation of skeletal muscle. Cell. 1993;75(7):1351–9.
Sadkowski T, Jank M, Zwierzchowski L, Oprzadek J, Motyl T. Comparison of skeletal muscle transcriptional profiles in dairy and beef breeds bulls. J Appl Genet. 2009;50(2):109–23.
Sadkowski T, Jank M, Zwierzchowski L, Siadkowska E, Oprządek J, Motyl T. Gene expression profiling in skeletal muscle of Holstein-Friesian bulls with single-nucleotide polymorphism in the myostatin gene 5’-flanking region. J Appl Genet. 2008;49(3):237–50.
Schabort EJ, van der Merwe M, Loos B, Moore FP, Niesler CU. TGF-β's delay skeletal muscle progenitor cell differentiation in an isoform-independent manner. Exp Cell Res. 2009;315(3):373–84.
Singh BN, Rao KS, Rao CM. Ubiquitin–proteasome-mediated degradation and synthesis of MyoD is modulated by αB-crystallin, a small heat shock protein, during muscle differentiation. Biochim Biophys Acta. 2010;1803(2):288–99.
Suleiman MS, Singh RJR, Stewart CEH. Apoptosis and the cardiac action of insulin-like growth factor I. Pharmacol Ther. 2007;114(3):278–94.
Swart GW. Activated leukocyte cell adhesion molecule (CD166/ALCAM): developmental and mechanistic aspects of cell clustering and cell migration. Eur J Cell Biol. 2002;81(6):313–21.
Szcześniak KA, Ostaszewski P, Fuller JC, Ciecierska A, Sadkowski T. Dietary supplementation of β‐hydroxy‐β‐methylbutyrate in animals—a review. J Anim Physiol Anim Nutr. 2015;99(3):405–17.
Szcześniak KA, Ostaszewski P, Ciecierska A, & Sadkowski T. Investigation of nutriactive phytochemical–gamma‐oryzanol in experimental animal models. J Anim Physiol Animal Nutr. 2015: doi:10.1111/jpn.12428.
Tao Y, Neppl RL, Huang ZP, Chen J, Tang RH, Cao R, et al. The histone methyltransferase Set7/9 promotes myoblast differentiation and myofibril assembly. J Cell Biol. 2011;194(4):551–65.
Tonami K, Hata S, Ojima K, Ono Y, Kurihara Y, Amano T, et al. Calpain-6 deficiency promotes skeletal muscle development and regeneration. PLoS Genet. 2013;9:e1003668. doi:10.1371/journal.pgen.1003668.
Trendelenburg AU, Meyer A, Jacobi C, Feige JN, Glass DJ. TAK-1/p38/nNF B signaling inhibits myoblast differentiation by increasing levels of Activin A. Skelet Muscle. 2012;2(1):3.
Wicik Z, Sadkowski T, Jank M, Motyl T. The transcriptomic signature of myostatin inhibitory influence on the differentiation of mouse C2C12 myoblasts. Pol J Vet Sci. 2011;14(4):643–52.
Xie X, Qin J, Lin SH, Tsai SY, Tsai MJ. Nuclear receptor chicken ovalbumin upstream promoter-transcription factor II (COUP-TFII) modulates mesenchymal cell commitment and differentiation. Proc Natl Acad Sci. 2011;108(36):14843–8.
Xie X, Tang K, Yu CT, Tsai SY, Tsai MJ. Regulatory potential of COUP-TFs in development: stem/progenitor cells. Semin Cell Dev Biol. 2013;24(10):687–93. Academic Press.
Zhang L, Wang XH, Wang H, Du J, Mitch WE. Satellite cell dysfunction and impaired IGF-1 signaling cause CKD-induced muscle atrophy. J Am Soc Nephrol. 2010;21(3):419–27.
Acknowledgements
Not applicable.
Funding
This research was funded by National Science Centre (Poland), Grant No. 2011/03/B/NZ5/05697. Publication of this manuscript was supported by KNOW (Leading National Research Centre) Scientific Consortium “Healthy Animal - Safe Food,” decision of the Ministry of Science and Higher Education No. 05-1/KNOW2/2015.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
KS carried out the muscle sampling, QPCR validation of microarray results, ontological analysis, and interpretation of the obtained data and drafted the manuscript. AC carried out the equine satellite cell isolation and culture, RNA isolation, and microarray analysis. PO participated in the study design and helped in the manuscript preparation. TS designed the study, supervised the project, performed the muscle sampling and statistical analysis of the microarray and QPCR data, and assisted in the manuscript preparation and revision. All authors read and approved the final manuscript.
Additional files
Additional file 1:
Full list of differentially expressed transcripts, gamma-oryzanol vs. control. (XLS 55.0 kb)
Additional file 2:
Full list of significantly enriched (EASE score p value ≤0.05) gene ontologies (DAVID). (XLS 42.5 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Szcześniak, K.A., Ciecierska, A., Ostaszewski, P. et al. Transcriptomic profile adaptations following exposure of equine satellite cells to nutriactive phytochemical gamma-oryzanol. Genes Nutr 11, 5 (2016). https://doi.org/10.1186/s12263-016-0523-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12263-016-0523-5