MiRNAs are essential mediators of many biological processes. The aim of this study was to investigate the dynamics of miRNA-mRNA regulatory networks during exercise and the subsequent recovery period.
Results
Here we monitored the transcriptome changes using microarray analysis of the whole blood of eight highly trained athletes before and after 30 min of moderate exercise followed by 30 min and 60 min of recovery period. We combined expression profiling and bioinformatics and analysed metabolic pathways enriched with differentially expressed mRNAs and mRNAs which are known to be validated targets of differentially expressed miRNAs. Finally we revealed four dynamically regulated networks comprising differentially expressed miRNAs and their known target mRNAs with anti-correlated expression profiles over time. The data suggest that hsa-miR-21-5p regulated TGFBR3, PDGFD and PPM1L mRNAs. Hsa-miR-24-2-5p was likely to be responsible for MYC and KCNJ2 genes and hsa-miR-27a-5p for ST3GAL6. The targets of hsa-miR-181a-5p included ROPN1L and SLC37A3. All these mRNAs are involved in processes highly relevant to exercise response, including immune function, apoptosis, membrane traffic of proteins and transcription regulation.
Conclusions
We have identified metabolic pathways involved in response to exercise and revealed four miRNA-mRNA networks dynamically regulated following exercise. This work is the first study to monitor miRNAs and mRNAs in parallel into the recovery period. The results provide a novel insight into the regulatory role of miRNAs in stress adaptation.
Keywords
ExerciseRegulationmiRNA-mRNA networks
Background
MiRNAs are one family of small (20–22 nucleotides) non-coding RNAs. They regulate gene expression post-transcriptionally through binding to the complementary sites of target mRNAs in the 3′-UTR, and play an important role in regulating diverse biological processes [1].
Recently, miRNA have been demonstrated as regulators of processes involved in physiological stress adaptation, including inflammation [2], angiogenesis [3], mitochondrial metabolism [4], muscle force generation [5]. However, just a few studies were published to date describing the changes in miRNA expression during exercise of different intensity [6–12]. They did not include the analysis of post-exercise recovery period and thus provided no information concerning dynamics of the predicted miRNA-mRNA regulatory pairs. Detailed investigation of miRNA-mRNA networks specifically regulated by exercise could reveal important biomarkers of exercise physiology and would provide for deep insight into the molecular control of the stress response.
MiRNAs regulate target gene expression in different ways including mRNA degradation and translation inhibition [1]. The target genes which were regulated by miRNAs through mRNA degradation are anti-correlated with the miRNA regulators. In this study, for the first time whole transcriptome changes were monitored during exercise followed by 30 min and 60 min of recovery period and differentially expressed mRNAs and miRNAs were analysed resulting in identification of four dynamically regulated miRNA-mRNA networks.
Results and discussion
Anthropometric and physiological data
To exclude possible effects of gender, only male subjects were recruited for this study. Anthropometric and physiological parameters of athletes are presented in Table 1. Before exercise the serum lactate level was 1.7 ± 0.4 mM. After exercise, it was mildly elevated, but still below 4.0 mM, confirming that the exercise performed was moderate without transgression of the anaerobic threshold.
Table 1
Anthropometric and physiological data
Parameter
Value
Age (year)
21.7 ± 2.6
Body mass (kg)
74.9 ± 2.3
Height (cm)
185.3 ± 3.5
Body mass index (kg/m2)
21.8 ± 1.2
Heart rate before exercise(bpm)
57.3 ± 1.9
Heart rate after exercise (bpm)
179.4 ± 3.2
Blood pressure before exercise ( mmHg)
120/66 ± 3/2
VO2max (ml min-1 kg-1)
74.8 ± 3.3
Values are mean ± SD.
Branched-chain amino acids (BCAA) include three structurally related amino acids Leucine (Leu), Isoleucine (Ile), and Valine (Val). The initial steps of their degradation are catalyzed by the same set of mitochondrial enzymes, and therefore, the BCAA behave as a very homogenous group. Their regulation is performed by short-term metabolic control reflecting consuming of energy sources. It has been shown previously that an acute bout of prolonged exercise increases the rate of BCAA oxidation by skeletal muscle [13]. We observed a slight increase in the BCAA level immediately after exercise followed by a decline below base level during recovery (Table 2). The ratio of citrulline (Cit) to ornithine (Orn) is indicative of the ornithine carbamoylphosphate transferase activity and characterizes the regulation of the urea cycle pathway [14]. This ratio had a tendency to increase (Table 2). The ratios of methionine sulfoxide (Met-SO) to methionine (Met) and tyrosine (Tyr) to phenylalanine (Phe) indicate oxidative stress [15, 16]. We found a decrease in Met-SO/Met and a slight increase in Tyr/Phe (Table 2).
The data summarized in Table 2 confirmed that the exercise was moderate and athletes reacted as normal healthy subjects [17].
Flow cytometry analysis
Changes in white blood cell subpopulations in response to exercise are presented in Table 3. Total white blood cell counts revealed the expected exercise-induced leukocytosis. NK lymphocytes (defined as CD3–, CD16/56+) substantially contributed to the observed changes which was consistent with the published data [18]. However NK-specific mRNAs (e.g., coding for KIR receptors) in the whole blood did not demonstrate similar nearly 3-fold increase (see Additional file 1) thus confirming that our subsequent microarray analysis showed true changes in RNA expression.
Table 3
Changes in white blood cell subpopulations in response to an exercise
Pre
Post (E)
Fold change relative
Absolute, cells/ul
Relative, %
Absolute, cells/ul
Relative, %
WBC
5710 ± 1170
7705 ± 2000
LY
2387 ± 674
41.8 ± 11.8
3989 ± 932
69.9 ± 12.1
1.67
СD3+
1616 ± 549
28.3 ± 9.6
2349 ± 901
41.1 ± 11.7
1.45
СD3 + CD4+
1020 ± 413
17.9 ± 7.2
1189 ± 419
20.8 ± 5.4
1.17
СD3 + CD8+
596 ± 265
10.4 ± 4.6
1124 ± 532
19.7 ± 6.9
1.89
СD3-CD16+
494 ± 203
8.7 ± 3.6
1380 ± 523
24.2 ± 6.8
2.79
СD3-CD56+
393 ± 217
6.9 ± 3.8
1111 ± 468
19.5 ± 6.1
2.83
CD19+
227 ± 124
4.0 ± 2.2
308 ± 221
5.4 ± 2.9
1.36
Values are mean ± SD.
miRNA and mRNA differential expression profiles
The miRNA and mRNA expression profiles in the whole blood for each time point and each athlete were determined using microarray analysis. PAXgene RNA tubes enable the isolation of intracellular RNA of circulating leukocytes, including B cells, T cells, neutrophils, monocytes, and other less abundant cell types. Furthermore, a large proportion of reticulocyte-derived globin mRNA is prepared from PAXgene blood RNA tubes, as it has been demonstrated previously [19]. Affymetrix GeneChip Human Gene 1.0 ST arrays contain both miRNA and mRNA probes and are capable of measuring about 20,000 mRNAs and 200 miRNAs.
The data are presented in Additional file 1. 298 mRNAs and 5 miRNAs were changed significantly (at least 40% change with the adjusted P-value threshold of 0.05), including hsa-miR-21-5p, hsa-miR-24-2-5p, hsa-miR-27a-5p, hsa-miR-181a-5p and hsa-miR-181b-5p. Remarkably, hsa-miR-24-2-5p is clustered with hsa-miR-27a-5p, and hsa-miR-181a-5p is clustered with hsa-miR-181b-5p. Consistently, the clustered miRNAs exhibited similar expression profiles over time (see Additional file 1 and Figures 1, 2, 3 and 4).
Figure 1
Expression profiles of hsa-miR-21-5p and its target mRNAs. (A-H) the tested athletes. E, after 30 min of exercise; R1, after 30 min relaxation; R2, after 60 min relaxation. Bold blue, hsa-miR-21-5p; brown, TGFBR3; green, PDGFD; violet, PPM1L; light blue, RHOBTB3. Solid lines indicate validated mRNA targets and dashed line indicates predicted potential mRNA target.
Figure 2
Expression profiles of hsa-miR-24-2-5p and its target mRNAs. (A-H) the tested athletes. E, after 30 min of exercise; R1, after 30 min relaxation; R2, after 60 min relaxation. Bold blue, hsa-miR-24-2-5p; brown, MYC; green, KCNJ2. Solid line indicates validated mRNA target and dashed line indicates predicted potential mRNA target.
Figure 3
Expression profiles of hsa-miR-27a-5p and its target mRNA. (A-H) the tested athletes. E, after 30 min of exercise; R1, after 30 min relaxation; R2, after 60 min relaxation. Bold blue, hsa-miR-27a-5p; dashed brown, predicted potential target mRNA ST3GAL6.
Figure 4
Expression profiles of hsa-miR-181a-5p and its validated target mRNAs. (A-H) the tested athletes. E, after 30 min of exercise; R1, after 30 min relaxation; R2, after 60 min relaxation. Bold blue, hsa-miR-181a-5p; brown, ROPN1L; green, SCL37A3.
Pathway analysis of differentially expressed mRNAs
All 298 differentially expressed mRNAs were analysed for enriched metabolic pathways. Table 4 indicates the revealed pathways including immune response and glycoproteins. As expected, they were previously reported to be associated with exercise [6, 8] thus confirming the relevance of our experimental model.
Table 4
Pathway analysis of differentially expressed mRNAs
Pathway
Number of genes
Adjusted P-value
Glycoprotein
76
3.00E-09
Natural killer cell mediated cytotoxicity
22
1.80E-15
Immune response
13
3.10E-03
Response to wounding
7
8.70E-01
Inflammatory response
7
1.90E-01
Regulation of lymphocyte mediated immunity
5
1.40E-01
Leukocyte activation
5
8.90E-01
Response to mechanical stimulus
3
3.80E-02
Stress response
3
5.90E-02
Cytolysis
3
3.20E-02
Pathway analysis of mRNA targets of differentially expressed miRNAs
All 5 differentially expressed miRNAs have 1136 validated target mRNAs in total. We performed pathway enrichment analysis for these mRNAs (Table 5). Again, the revealed pathways are highly relevant to exercise, e.g. transcription regulation, apoptosis, response to stress etc.
Table 5
Pathway analysis of all validated mRNA targets of differentially expressed miRNAs
Pathway
Number of genes
Adjusted P-value
Transcription regulation
283
2.03E-13
Apoptotic process
157
2.59E-11
Cell cycle
148
5.95E-12
Regulation of kinase activity
76
1.08E-08
Regulation of cellular response to stress
34
1.63E-04
p53 signaling pathway
19
4.26E-08
E2F transcription factor network
19
1.51E-07
Validated targets of C-MYC transcriptional activation
14
1.08E-03
It has been demonstrated that the same mRNAs can be targeted by more than one miRNA which provides for more efficient and specific regulation [20, 21]. We found 49 mRNAs which are known to be validated targets for 2 or even 3 differentially expressed miRNAs. They have higher potential to be involved in exercise-induced regulation. Table 6 shows the pathways enriched with some of these mRNAs. Notably, these exercise-relevant pathways (cell death, stress response, proliferation) comprise significant number of intersecting mRNAs. Figure 5 presents the identified regulatory miRNA-mRNA network for all 3 revealed pathways. Interestingly, some of these genes are known to interact with each other. Namely transcription factor MYC was reported to be functionally associated with RNA helicase DDX3X [22], apoptosis regulator BCL2 [23] and tumor suppressor BRCA1 [24]. BRCA1 in turn interacts itself with BCL2 [25] and transcription factor E2F1 [26]. The presented data support the regulatory role of identified miRNAs in response to exercise.
Figure 5
Regulatory network of mRNAs validated to be the targets for 2 or 3 differentially expressed miRNAs. All mRNAs belong to three enriched pathways listed in Table 6. Red dashed lines connect mRNAs of proteins reported to functionally interact with each other. MYC is highlighted with red oval indicating that this mRNA was further identified to be a component of dynamically regulated miRNA-mRNA network.
Table 6
Pathway analysis of mRNAs validated to be targets for 2 or 3 differentially expressed miRNAs
Pathway
Number of genes
Adjusted P-value
Cell death
17
7.64E-07
Regulation of cell proliferation
15
2.10E-07
Cellular response to stress
12
5.35E-05
Dynamically regulated miRNA-mRNA networks
We monitored the transcriptome expression level before and following exercise and this allowed us to reveal dynamically regulated miRNA-mRNA networks. We used a two-step approach to identify the mRNA targets for the differentially expressed miRNAs. First, we looked for anti-correlating groups of miRNAs and mRNAs expression of which over time tended to have opposite profiles. Target mRNA degradation is one of the mechanisms of miRNA action when their perfect complementarity occurs [1, 27]. Thus, the second step of our strategy was either analysis of published data and selection of experimentally validated target mRNAs for a given miRNA or theoretical prediction of miRNA targets based on their complementarity, using one of the most popular web resources TargetScan [28]. The expression profiles and biological function of selected miRNAs and mRNAs were analysed in more detail. Based on this analysis, our final identified miRNA-mRNA pairs have a high probability of being involved in the regulation of exercise-related physiological processes.
hsa-miR-21-5p
MiRNA hsa-miR-21-5p demonstrated different expression profiles over time (adjusted P-value 0.0039) however remarkably anti-correlating with experimentally validated target mRNAs TGFBR3 (adjusted P-value 1.079E-06), PDGFD (adjusted P-value 1.28E-06) and PPM1L (adjusted P-value 6.78E-07) (Figure 1). RHOBTB3 mRNA (adjusted P-value 0.002) predicted by TargetScan to be a potential target for hsa-miR-21-2-5p behaved similarly. We observed the up-regulation of hsa-miR-21-5p one hour after exercise. The differences in kinetics can be potentially explained by individualities of each athlete.
The expression level of hsa-miR-21-5p itself is known to be stress-responsive and play an important role in heart failure [29] and renal ischemia reperfusion injury [30]. Notably the up-regulation of circulating hsa-miR-21-5p was recently reported to occur in plasma upon exercise [6]. The overall action of hsa-miR-21-5p has been described by several authors to be strongly anti-inflammatory [6, 31]. Note, that 60 min into relaxation, there was an up-regulation of hsa-miR-21-5p in all 8 subjects. This may reflect the self-protective anti-inflammatory reaction to exercise.
We identified four pairing targets for this miRNA namely TGFBR3, PDGFD, PPM1L, and RHOBTB3. TGFBR3 is a transforming growth factor (TGF)-beta type III receptor, mRNA of which is known to be up-regulated in the peripheral blood leukocytes in allograft rejection-prone recipients after intestinal transplantation thus mediating innate and adaptive inflammatory functions of leukocytes [32]. PDGFD encodes for platelet derived growth factor D [33], a member of the platelet-derived growth factor family which can regulate many cellular processes, including cell proliferation, apoptosis, transformation, migration, invasion, angiogenesis and metastasis [34]. PPM1L encodes a protein phosphatase gene [35] responsible for the regulation of stress-activated protein kinase signaling cascade and apoptosis [36]. Finally, the fourth identified target for hsa-miR-21-5p is RHOBTB3 mRNA [33] which encodes for the Rho GTPase regulating the membrane traffic of proteins [37]. Interestingly, this mRNA proved to be a blood biomarker of psychosis and shows a decreased expression level in high hallucinations states [38].
hsa-miR-24-2-5p
MiRNA hsa-miR-24-2-5p (adjusted P-value 0.00017) was up-regulated immediately after exercise, then tended to decrease during the recovery period except athlete D (Figure 2). In our study for the first time we report the reaction to exercise of this miRNA which is known to have protective effects on myocytes against myocardial ischaemia/reperfusion-induced apoptosis [39].
MYC mRNA (adjusted P-value 0.00013) which is known to be the target for this miRNA from the literature, and KCNJ2 mRNA (adjusted P-value 0.00023) predicted by TargetScan to be a potential target essentially followed anti-profile of hsa-miR-24-2-5p. The protein encoded by the MYC gene is a multifunctional, nuclear phosphoprotein that plays a role in cell cycle progression, apoptosis and cellular transformation. It functions as a transcription factor that regulates transcription of specific target genes. Interestingly, hsa-miR-24-2-5p is known to be up-regulated during hematopoietic cell terminal differentiation suppressing MYC expression [40]. Thus the increase of this miRNA we observed (Figure 2) may reflect the regulation of hematopoiesis upon exercise. Remarkably, MYC mRNA is among those which are regulated by 2 differentially expressed miRNAs, namely hsa-miR-21-5p and hsa-miR-24-2-5p (Figure 5), however its expression profile is anti-correlated with only the profile of hsa-miR-24-2-5p (Figures 1 and 2).
KCNJ2 protein is an integral membrane protein and inward-rectifier type potassium channel participating in establishing the action potential waveform and excitability of neuronal and muscle tissues [41]. This mRNA expressed in peripheral blood lymphocytes is a biomarker for Parkinson’s disease [42].
hsa-miR-27a-5p
MiRNA hsa-miR-27a-5p is clustered with hsa-miR-24-2-5p and behaved similarly to it increasing after exercise and decreasing during the recovery period except athlete D (Figure 3) with the adjusted P-value 0.00012. This miRNA was reported to promote myoblast proliferation by reducing the expression of myostatin [43].
The only mRNA target identified is predicted by TargetScan ST3GAL6. The encoded protein belongs to the sialyltransferase family and is responsible for the synthesis of selectin ligands [44].
hsa-miR-181a-5p
hsa-miR-181a-5p tended to increase after exercise and then to down-regulate during the first as well as the second period of the relaxation time (Figure 4) with an adjusted P-value of 5.83E-05. The observed differential expression of hsa-miR-181a-5p in our athletes is consistent with previously published results [7, 8]. This miRNA is characterized as a regulator of hematopoietic lineage differentiation [45] and a modulator of T cell sensitivity and selection [46]. Radom-Aizik showed up-regulation of this miRNA after 30 min interval exercise. They related it to increased T cell responsiveness and reduced susceptibility to infection due to physical activity. In our study 7 subjects showed up-regulation immediately after exercise (Figure 4).
The mRNAs ROPN1L (adjusted P-value 0.00024) and SLC37A3 (adjusted P-value 0.0019) were previously validated to be targets for hsa-miR-181a-5p and demonstrated pronounced anti-correlation with the miRNA expression profile. The ROPN1L gene encodes a member of the ropporin family. The encoded protein is involved in the targeting towards specific physiological substrates of Protein Kinase A, regulating glycogen, sugar, and lipid metabolism [47]. The SLC37A3 protein belongs to transmembrane sugar transporters and is responsible for sugar metabolism [48].
Conclusion
We have identified metabolic pathways enriched with differentially expressed mRNAs and with mRNA targets of differentially expressed miRNAs, including mRNAs known to be regulated by 2 or 3 miRNAs described here. The result supports previously published data. Moreover, we revealed four miRNA-mRNA networks dynamically regulated following exercise. These observations provide a novel insight into the potential regulatory role of miRNAs in the numerous physiological processes involved in stress adaptation.
Methods
Ethical approval and study participants
Eight national level ski athletes took part in this study. None of them suffered from acute or chronic diseases or reported intake of medication. Participants were informed about the nature, purpose, and potential risks of the experiments and signed an informed consent statement approved by the ethics committee of Scientific Research Center Bioclinicum (Moscow, Russia).
Anthropometric measurements
Height, weight, medical historical data and resting vital signs were recorded at the time of enrolment.
Exercise test protocol
In order to determine the VO2max values, each subject performed a treadmill test with an incremental step protocol until exhaustion as described previously [47]. VO2max was calculated as described [49]. The anaerobic threshold was calculated using the standard V-slope method [50].
Two weeks later, athletes participated in the main exercise, consisting of running at 80% VO2 peak for 30 min on a treadmill. The exercise was performed during the morning hours (between 8 and 11 a.m.) keeping the exact test time for each participant constant.
Lactate concentration
Lactate concentration in capillary blood was measured electrochemically using the automated analyzer Biosen C_Line (EFK Diagnostic, Germany).
Analysis of aminoacids
Sera derived from the venous blood samples were analyzed for aminoacids in Genome Analysis Centre (Helmholtz Zentrum, Munich, Germany) using the BIOCRATES AbsoluteIDQ p150 Kit (Biocrates Life Sciences AG, Innsbruck, Austria) in a combined FIA-MS/MS and LC-MS/MS assay as recommended by the manufacturer and described previously [51].
Blood sampling
Venous blood was collected at four time points during the exercise. A 20-gauge intravenous catheter was placed antiseptically into a dorsal hand vein or a vein in the distal forearm as dictated by favourable anatomy using the Seldinger technique and then secured with tape. During ME 2.5 ml of blood was collected in PAXgene blood RNA tube, 7 ml in a Serum separator tube (BD, USA) and 4.5 ml in a tube containing buffered tri-sodium citrate (BD, USA) for flow cytometry analysis at baseline (prior to exercise testing) and immediately post-exercise (within 1 min of completion of exercise testing). After 30 min of rest and after 60 min of rest following exercise testing 2.5 ml of blood was collected in PAXgene blood RNA tubes.
RNA extraction
According to the Affymetrix Manual P/N 701880 Rev. 4 total RNA was extracted using the PAXgene Blood RNA kit as recommended by the manufacturer. RNA concentrations were determined by the Nanodrop photometer (NanoDrop, USA). RNA quality was checked using the Agilent Bioanalyser 2100 System (Agilent Technologies, USA). For all samples RNA integrity number (RIN) was greater than 7.
Flow cytometry analysis
Flow cytometry analysis of blood samples was performed using fluorescently labeled antibodies against B and T cell receptors and natural killer (NK) cell markers. Cells were labeled with antibodies against CD3 (FITC), CD4 (PE), CD8 (PE), CD16 (FITC), CD56 (PE), CD19 (FITC) (Sorbent, Russia), where NK cells were distinguished from the rest of the lymphocytes via positive expression of CD56 and negative expression of CD3.
The samples were analyzed on a FACScan Calibur flow cytometer (BD Biosciences, USA) and leukocytes were gated based on forward and side scatter properties. Events in the range of 40,000–200,000 were collected depending on the occurrence of the investigated leukocyte population, and analyzed with CELLQuest Pro analysis software (BD Biosciences, USA). To ensure flow cytometric standardization, the voltage settings were updated daily using ‘Calibrate Beads’ (BD Biosciences, USA).
Microarray analysis
RNA samples were prepared according to manufacturer’s instructions (Affymetrix Manual P/N 701880 Rev. 4) as described elsewhere [52]. The samples were hybridized on GeneChip Human Gene 1.0 ST Arrays containing both miRNA and mRNA probes (Affymetrix, USA) for 16 h at 45°C. Arrays were washed to remove non-specifically bound nucleic acids and stained on Fluidics Station 450 (Affymetrix) using FS450_0007 protocol followed by scanning on a GeneChip Scanner 3000 7G system (Affymetrix). The microarray CEL files have been deposited in the GEO database (accession GSE46075).
Microarray data processing
Microarray data was processed using bioconductor [53] xps package implementation of RMA [54]. At the first step background correction was performed based on a global model for the distribution of probe intensities [54]. Then a quantile normalization algorithm [55] (so-called probe-level normalization) was applied to the preprocessed data. Finally, fitting a robust linear model using Tukey’s median polish procedure [56] was done to convert probe intensities to the expression levels of probesets.
The statistical analysis of microarray data was performed using bioconductor [53] package limma [57]. The analysis was based on a generalized linear model [57, 58] approach. In this approach one constructs a linear data model with a structure determined by the experiment layout, and then fits this model to the actual data. The linear model is defined in terms of a so-called design matrix. The number of rows in this matrix coincides with the number of experiment samples, and the number of columns coincides with the number of factors that have an essential influence on the measurable values. The value at the i-th row and j-th column of a design matrix specifies an effect of a factor j on a sample i. Each measurable value (i.e., each probeset) in this approach is analysed independently. For each probeset a vector of its expression values E is represented in the form E = Dβ + ϵ, where D is a design matrix, β is a vector of coefficients indicating values of each factor’s actual influence on the analyzed probeset, ϵ is a vector of error, and model fitting consists in the minimization of “error term” ϵ by finding optimal coefficients β. After the coefficients β are computed for each probeset, one can test various hypotheses on the structure of considered factors. For example, in order to find probesets that are affected by a factor i, one should search for the probesets with βi statistically different from zero.
The general linear model was applied in the analysis of the studied transcriptome changes as follows. The model that takes into account both experiment time points and athletes individual features was used. Thus, the total number of factors was eleven: c1, …, c8 correspond to an expression level of each athlete in a normal state (for all samples the effect of these factors is set to 1 if a sample and a factor correspond to the same athlete, and set to 0 otherwise), and d1, d2, d3 correspond to the changes induces by the exercises, exercises and 30-minutes relaxation, and exercises and 60-minutes relaxation respectively (for all samples the effect of these factors is set to 1 if a sample and a factor correspond to the same experiment time point, and set to 0 otherwise). The total number of samples was 32: 4 samples for each athlete.
For each pair of experiment time points the detection of probesets with reliable difference between time points was performed. The probeset was considered to have a reliable difference between time points k, m if an adjusted p-value of an equality βk = βm (where βl was an actual influence of exercises at experiment time point for l = 1,2,3, and β0 = 0) was less than 0.05, and a log-fold change was greater than 0.484 (this threshold corresponds to an intensity change by more than 40%). The Benjamini-Hochberg [59] algorithm was used for multiple testing adjustment. The minimum adjusted p-value for all pairs of time points is indicated for each differentially expressed probeset in the Table of differentially expressed mRNAs and miRNAs (Additional file 1) as a statistical significance value.
Pathway analysis
Bioinformatic analysis was performed using DAVID online tool (http://david.abcc.ncifcrf.gov) as described elsewhere [60]. So all the analyzed genes were classified into several functional groups, and the groups that may be potentially associated with physiological stress were considered and listed in the tabs of excel document. P-values on the tabs are modified Fisher Exact P-Values. When members of two independent groups can fall into one of two mutually exclusive categories, Fisher Exact test is used to determine whether the proportions of those falling into each category differs by group. In DAVID annotation system, Fisher Exact test is adopted to measure the gene-enrichment in annotation terms.
Abbreviations
miRNA:
microRNA
mRNA:
messenger RNA.
Declarations
Acknowledgements
We thank Svetlana Vinogradova for the help with bioinformatic data analysis. This work was supported by BMBF grant RUS 10/040 and by Russian Ministry of Science Grants No. 16.522.12.2015 and 16.522.11.2004.
Authors’ Affiliations
(1)
Moscow State University
(2)
The Institute of General Pathology and Pathophysiology, Russian Academy of Medical Sciences
(3)
SRC Bioclinicum
(4)
Institute of Clinical and Experimental Transfusion Medicine (IKET), University of Tübingen
(5)
Institute for Biomedical Problems, Russian Academy of Sciences
Davidson-Moncada J, Papavasiliou FN, Tam W: MicroRNAs of the immune system: roles in inflammation and cancer.Ann N Y Acad Sci 2010, 1183: 183-194. 10.1111/j.1749-6632.2009.05121.xView ArticlePubMedPubMed CentralGoogle Scholar
Dang CV: Rethinking the Warburg effect with Myc micromanaging glutamine metabolism.Cancer Res 2010, 70: 859-862. 10.1158/0008-5472.CAN-09-3556View ArticlePubMedPubMed CentralGoogle Scholar
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-317. 10.1152/japplphysiol.00901.2010View ArticlePubMedGoogle Scholar
Baggish AL, Hale A, Weiner RB, Lewis GD, Systrom D, Wang F, Wang TJ, Chan SY: Dynamic regulation of circulating microRNA during acute exhaustive exercise and sustained aerobic exercise training.J Physiol 2011, 589: 3983-3994. 10.1113/jphysiol.2011.213363View ArticlePubMedPubMed CentralGoogle Scholar
Radom-Aizik S, Zaldivar F Jr, Oliver S, Galassetti P, Cooper DM: Evidence for microRNA involvement in exercise-associated neutrophil gene expression changes.J Appl Physiol 2010, 109: 252-261. 10.1152/japplphysiol.01291.2009View ArticlePubMedPubMed CentralGoogle Scholar
Radom-Aizik S, Zaldivar F Jr, Leu SY, Adams GR, Oliver S, Cooper DM: Effects of exercise on microRNA expression in young males peripheral blood mononuclear cells.Clin Transl Sci 2012, 5: 32-38. 10.1111/j.1752-8062.2011.00384.xView ArticlePubMedPubMed CentralGoogle Scholar
Fernandes T, Magalhães FC, Roque FR, Phillips MI, Oliveira EM: Exercise training prevents the microvascular rarefaction in hypertension balancing angiogenic and apoptotic factors: role of microRNAs-16, -21, and −126.Hypertension 2012, 59: 513-520. 10.1161/HYPERTENSIONAHA.111.185801View ArticlePubMedGoogle Scholar
Radom-Aizik S, Zaldivar FP, Haddad F, Cooper DM: Impact of brief exercise on peripheral blood NK cell gene and microRNA expression in young adults.J Appl Physiol 2013, 114: 628-636. 10.1152/japplphysiol.01341.2012View ArticlePubMedPubMed CentralGoogle Scholar
Uhlemann M, Möbius-Winkler S, Fikenzer S, Adam J, Redlich M, Möhlenkamp S, Hilberg T, Schuler GC, Adams V: Circulating microRNA-126 increases after different forms of endurance exercise in healthy adults.Eur J Prev Cardiol 2012. 10.1177/2047487312467902Google Scholar
Bye A, Røsjø H, Aspenes ST, Condorelli G, Omland T, Wisløff U: Circulating MicroRNAs and Aerobic Fitness - The HUNT-Study.PLoS One 2013, 8: e57496. 10.1371/journal.pone.0057496View ArticlePubMedPubMed CentralGoogle Scholar
McKenzie S, Phillips SM, Carter SL, Lowther S, Gibala MJ, Tarnopolsky MA: Endurance exercise training attenuates leucine oxidation and BCOAD activation during exercise in humans.Am J Physiol Endocrinol Metab 2000, 278: E580-587.PubMedGoogle Scholar
Huang CC, Lin WT, Hsu FL, Tsai PW, Hou CC: Metabolomics investigation of exercise-modulated changes in metabolism in rat liver after exhaustive and endurance exercises.Eur J Appl Physiol 2010, 108: 557-566. 10.1007/s00421-009-1247-7View ArticlePubMedGoogle Scholar
Ribas GS, Sitta A, Wajner M, Vargas CR: Oxidative stress in phenylketonuria: what is the evidence?Cell Mol Neurobiol 2011, 31: 653-662. 10.1007/s10571-011-9693-2View ArticlePubMedGoogle Scholar
Mashima R, Nakanishi-Ueda T, Yamamoto Y: Simultaneous determination of methionine sulfoxide and methionine in blood plasma using gas chromatography–mass spectrometry.Anal Biochem 2003, 313: 28-33. 10.1016/S0003-2697(02)00537-7View ArticlePubMedGoogle Scholar
Gibala MJ: Protein metabolism and endurance exercise.Sports Med 2007, 37: 337-340. 10.2165/00007256-200737040-00016View ArticlePubMedGoogle Scholar
Horn P, Kalz A, Lim CL, Pyne D, Saunders P, Mackinnon L, Peake J, Suzuki K: Exercise-recruited NK cells display exercise-associated eHSP-70.Exerc Immunol Rev 2007, 13: 100-111.PubMedGoogle Scholar
Feezor RJ, Baker HV, Mindrinos M, Hayden D, Tannahill CL, Brownstein BH, Fay A, MacMillan S, Laramie J, Xiao W, Moldawer LL, Cobb JP, Laudanski K, Miller-Graziano CL, Maier RV, Schoenfeld D, Davis RW, Tompkins RG, Inflammation and Host Response to Injury, Large-Scale Collaborative Research Program: Whole blood and leukocyte RNA isolation for gene expression analyses.Physiol Genomics 2004, 19: 247-254. 10.1152/physiolgenomics.00020.2004View ArticlePubMedGoogle Scholar
Broderick JA, Salomon WE, Ryder SP, Aronin N, Zamore PD: Argonaute protein identity and pairing geometry determine cooperativity in mammalian RNA silencing.RNA 2011, 17: 1858-1869. 10.1261/rna.2778911View ArticlePubMedPubMed CentralGoogle Scholar
Takimoto K, Wakiyama M, Yokoyama S: Mammalian GW182 contains multiple Argonaute-binding sites and functions in microRNA-mediated translational repression.RNA 2009, 15: 1078-1089. 10.1261/rna.1363109View ArticlePubMedPubMed CentralGoogle Scholar
Davidson ME, Kerepesi LA, Soto A, Chan VT: D-Serine exposure resulted in gene expression changes implicated in neurodegenerative disorders and neuronal dysfunction in male Fischer 344 rats.Arch Toxicol 2009, 83: 747-762. 10.1007/s00204-009-0405-3View ArticlePubMedGoogle Scholar
Jin Z, May WS, Gao F, Flagg T, Deng X: Bcl2 suppresses DNA repair by enhancing c-Myc transcriptional activity.J Biol Chem 2006, 281: 14446-14456. 10.1074/jbc.M511914200View ArticlePubMedGoogle Scholar
Tonic I, Yu WN, Park Y, Chen CC, Hay N: Akt activation emulates Chk1 inhibition and Bcl2 overexpression and abrogates G2 cell cycle checkpoint by inhibiting BRCA1 foci.J Biol Chem 2010, 285: 23790-23798. 10.1074/jbc.M110.104372View ArticlePubMedPubMed CentralGoogle Scholar
De Siervi A, De Luca P, Byun JS, Di LJ, Fufa T, Haggerty CM, Vazquez E, Moiola C, Longo DL, Gardner K: Transcriptional autoregulation by BRCA1.Cancer Res 2010, 70: 532-542. 10.1158/0008-5472.CAN-09-1477View ArticlePubMedPubMed CentralGoogle Scholar
Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets.Cell 2005, 120: 15-20. 10.1016/j.cell.2004.12.035View ArticlePubMedGoogle Scholar
Thum T, Gross C, Fiedler J, Fischer T, Kissler S, Bussen M, Galuppo P, Just S, Rottbauer W, Frantz S, Castoldi M, Soutschek J, Koteliansky V, Rosenwald A, Basson MA, Licht JD, Pena JT, Rouhanifard SH, Muckenthaler MU, Tuschl T, Martin GR, Bauersachs J, Engelhardt S: MicroRNA-21 contributes to myocardial disease by stimulating MAP kinase signalling in fibroblasts.Nature 2008, 456: 980-984. 10.1038/nature07511View ArticlePubMedGoogle Scholar
Godwin JG, Ge X, Stephan K, Jurisch A, Tullius SG, Iacomini J: Identification of a microRNA signature of renal ischemia reperfusion injury.Proc Natl Acad Sci USA 2010, 107: 14339-14344. 10.1073/pnas.0912701107View ArticlePubMedPubMed CentralGoogle Scholar
Wang Z, Ahmad A, Li Y, Kong D, Azmi AS, Banerjee S, Sarkar FH: Emerging roles of PDGF-D signaling pathway in tumor development and progression.Biochim Biophys Acta 1806, 2010: 122-130.Google Scholar
Hafner M, Landthaler M, Burger L, Khorshid M, Hausser J, Berninger P, Rothballer A, Ascano M Jr, Jungkamp AC, Munschauer M, Ulrich A, Wardle GS, Dewell S, Zavolan M, Tuschl T: Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP.Cell 2010, 141: 129-141. 10.1016/j.cell.2010.03.009View ArticlePubMedPubMed CentralGoogle Scholar
Espinosa EJ, Calero M, Sridevi K, Pfeffer SR: RhoBTB3: a Rho GTPase-family ATPase required for endosome to Golgi transport.Cell 2009, 137: 938-948. 10.1016/j.cell.2009.03.043View ArticlePubMedPubMed CentralGoogle Scholar
Kurian SM, Le-Niculescu H, Patel SD, Bertram D, Davis J, Dike C, Yehyawi N, Lysaker P, Dustin J, Caligiuri M, Lohr J, Lahiri DK, Nurnberger JI Jr, Faraone SV, Geyer MA, Tsuang MT, Schork NJ, Salomon DR, Niculescu AB: Identification of blood biomarkers for psychosis using convergent functional genomics.Mol Psychiatry 2011, 16: 37-58. 10.1038/mp.2009.117View ArticlePubMedGoogle Scholar
Zhu H, Fan GC: Role of microRNAs in the reperfused myocardium towards post-infarct remodelling.Cardiovasc Res 2012, 94: 284-292. 10.1093/cvr/cvr291View ArticlePubMedPubMed CentralGoogle Scholar
Lal A, Navarro F, Maher CA, Maliszewski LE, Yan N, O’Day E, Chowdhury D, Dykxhoorn DM, Tsai P, Hofmann O, Becker KG, Gorospe M, Hide W, Lieberman J: miR-24 Inhibits cell proliferation by targeting E2F2, MYC, and other cell-cycle genes via binding to “seedless” 3‘UTR microRNA recognition elements.Mol Cell 2009, 35: 610-625. 10.1016/j.molcel.2009.08.020View ArticlePubMedPubMed CentralGoogle Scholar
Burge JA, Hanna MG: Novel insights into the pathomechanisms of skeletal muscle channelopathies.Curr Neurol Neurosci Rep 2012, 12: 62-69. 10.1007/s11910-011-0238-3View ArticlePubMedGoogle Scholar
Gui YX, Wan Y, Xiao Q, Wang Y, Wang G, Chen SD: Verification of expressions of Kir2 as potential peripheral biomarkers in lymphocytes from patients with Parkinson’s disease.Neurosci Lett 2011, 505: 104-108. 10.1016/j.neulet.2011.09.070View ArticlePubMedGoogle Scholar
Huang Z, Chen X, Yu B, He J, Chen D: MicroRNA-27a promotes myoblast proliferation by targeting myostatin.Biochem Biophys Res Commun 2012, 423: 265-269. 10.1016/j.bbrc.2012.05.106View ArticlePubMedGoogle Scholar
Yang WH, Nussbaum C, Grewal PK, Marth JD, Sperandio M: Coordinated roles of ST3Gal-VI and ST3Gal-IV sialyltransferases in the synthesis of selectin ligands.Blood 2012, 120: 1015-1026. 10.1182/blood-2012-04-424366View ArticlePubMedPubMed CentralGoogle Scholar
Li QJ, Chau J, Ebert PJ, Sylvester G, Min H, Liu G, Braich R, Manoharan M, Soutschek J, Skare P, Klein LO, Davis MM, Chen CZ: miR-181a is an intrinsic modulator of T cell sensitivity and selection.Cell 2007, 129: 147-161. 10.1016/j.cell.2007.03.008View ArticlePubMedGoogle Scholar
Chen L, Kass RS: A-kinase anchoring proteins: different partners, different dance.Nat Cell Biol 2005, 7: 1050-1051.View ArticlePubMedGoogle Scholar
Bartoloni L, Antonarakis SE: The human sugar-phosphate/phosphate exchanger family SLC37.Pflugers Arch 2004, 447: 780-783. 10.1007/s00424-003-1105-0View ArticlePubMedGoogle Scholar
Cooper DM, Weiler-Ravell D, Whipp BJ, Wasserman K: Aerobic parameters of exercise as a function of body size during growth in children.J Appl Physiol 1984, 56: 628-634.PubMedGoogle Scholar
Beaver WL, Wasserman K, Whipp BJ: A new method for detecting anaerobic threshold by gas exchange.J Appl Physiol 1986, 60: 2020-2027.PubMedGoogle Scholar
Floegel A, Stefan N, Yu Z, Mühlenbruch K, Drogan D, Joost HG, Fritsche A, Häring HU, Hrabe de Angelis M, Peters A, Roden M, Prehn C, Wang-Sattler R, Illig T, Schulze MB, Adamski J, Boeing H, Pischon T: Identification of serum metabolites associated with risk of type 2 diabetes using a targeted metabolomic approach.Diabetes 2013, 62: 639-648. 10.2337/db12-0495View ArticlePubMedPubMed CentralGoogle Scholar
Sakharov DA, Maltseva DV, Riabenko EA, Shkurnikov MU, Northoff H, Tonevitsky AG, Grigoriev AI: Passing the anaerobic threshold is associated with substantial changes in the gene expression profile in white blood cells.Eur J Appl Physiol 2012, 112: 963-972. 10.1007/s00421-011-2048-3View ArticlePubMedGoogle Scholar
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics.Genome Biol 2004, 5: R80. 10.1186/gb-2004-5-10-r80View ArticlePubMedPubMed CentralGoogle Scholar
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data.Biostatistics 2003, 4: 249-264. 10.1093/biostatistics/4.2.249View ArticlePubMedGoogle Scholar
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias.Bioinformatics 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185View ArticlePubMedGoogle Scholar
Tukey JW: Exploratory data analysis. Reading: Addison-Wesley; 1977.Google Scholar
Smyth GK: Limma: linear models for microarray data. In Bioinformatics and computational biology solutions using R and bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. New York: Springer; 2005:397-420.View ArticleGoogle Scholar
Yang YH, Speed TP: Design and analysis of comparative microarray experiments. In: Statistical analysis of gene expression microarray data. Edited by Speed TP. Chapman and Hall/CRC Press; 2003:35-93.Google Scholar
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing.J R Statist Soc B (Methodological) 1995, 57: 289-300.Google Scholar
da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.Nat Protoc 2009, 4: 44-57.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
