Dynamically regulated miRNA-mRNA networks revealed by exercise
© Tonevitsky et al.; licensee BioMed Central Ltd. 2013
Received: 11 February 2013
Accepted: 28 May 2013
Published: 7 June 2013
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.
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.
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.
KeywordsExercise Regulation miRNA-mRNA networks
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 .
Recently, miRNA have been demonstrated as regulators of processes involved in physiological stress adaptation, including inflammation , angiogenesis , mitochondrial metabolism , muscle force generation . 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 . 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
Anthropometric and physiological data
21.7 ± 2.6
Body mass (kg)
74.9 ± 2.3
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
Aminoacids before and in response to exercise
567.63 ± 172.16
0.35 ± 0.09
2.28 ± 1.73
0.83 ± 0.15
612.83 ± 99.68
0.34 ± 0.11
0.92 ± 0.21
0.91 ± 0.18
Rest 30 min (R1)
496.89 ± 128.31
0.42 ± 0.09
1.71 ± 0.92
0.91 ± 0.19
Rest 60 min (R2)
417.75 ± 58.2
0.5 ± 0.15
1.21 ± 0.47
0.9 ± 0.16
Flow cytometry analysis
Changes in white blood cell subpopulations in response to an exercise
Fold change relative
5710 ± 1170
7705 ± 2000
2387 ± 674
41.8 ± 11.8
3989 ± 932
69.9 ± 12.1
1616 ± 549
28.3 ± 9.6
2349 ± 901
41.1 ± 11.7
СD3 + CD4+
1020 ± 413
17.9 ± 7.2
1189 ± 419
20.8 ± 5.4
СD3 + CD8+
596 ± 265
10.4 ± 4.6
1124 ± 532
19.7 ± 6.9
494 ± 203
8.7 ± 3.6
1380 ± 523
24.2 ± 6.8
393 ± 217
6.9 ± 3.8
1111 ± 468
19.5 ± 6.1
227 ± 124
4.0 ± 2.2
308 ± 221
5.4 ± 2.9
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 . 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.
Pathway analysis of differentially expressed mRNAs
Pathway analysis of differentially expressed mRNAs
Number of genes
Natural killer cell mediated cytotoxicity
Response to wounding
Regulation of lymphocyte mediated immunity
Response to mechanical stimulus
Pathway analysis of mRNA targets of differentially expressed miRNAs
Pathway analysis of all validated mRNA targets of differentially expressed miRNAs
Number of genes
Regulation of kinase activity
Regulation of cellular response to stress
p53 signaling pathway
E2F transcription factor network
Validated targets of C-MYC transcriptional activation
Pathway analysis of mRNAs validated to be targets for 2 or 3 differentially expressed miRNAs
Number of genes
Regulation of cell proliferation
Cellular response to stress
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 . 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.
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  and renal ischemia reperfusion injury . Notably the up-regulation of circulating hsa-miR-21-5p was recently reported to occur in plasma upon exercise . 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 . PDGFD encodes for platelet derived growth factor D , 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 . PPM1L encodes a protein phosphatase gene  responsible for the regulation of stress-activated protein kinase signaling cascade and apoptosis . Finally, the fourth identified target for hsa-miR-21-5p is RHOBTB3 mRNA  which encodes for the Rho GTPase regulating the membrane traffic of proteins . Interestingly, this mRNA proved to be a blood biomarker of psychosis and shows a decreased expression level in high hallucinations states .
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 .
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 . 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 . This mRNA expressed in peripheral blood lymphocytes is a biomarker for Parkinson’s disease .
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 .
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 .
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  and a modulator of T cell sensitivity and selection . 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 . The SLC37A3 protein belongs to transmembrane sugar transporters and is responsible for sugar metabolism .
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.
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).
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 . VO2max was calculated as described . The anaerobic threshold was calculated using the standard V-slope method .
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 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 .
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.
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).
RNA samples were prepared according to manufacturer’s instructions (Affymetrix Manual P/N 701880 Rev. 4) as described elsewhere . 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  xps package implementation of RMA . At the first step background correction was performed based on a global model for the distribution of probe intensities . Then a quantile normalization algorithm  (so-called probe-level normalization) was applied to the preprocessed data. Finally, fitting a robust linear model using Tukey’s median polish procedure  was done to convert probe intensities to the expression levels of probesets.
The statistical analysis of microarray data was performed using bioconductor  package limma . 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  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.
Bioinformatic analysis was performed using DAVID online tool (http://david.abcc.ncifcrf.gov) as described elsewhere . 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.
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.
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5View ArticlePubMedGoogle Scholar
- 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
- Zhang C: MicroRNAs in vascular biology and vascular disease. J Cardiovasc Transl Res 2010, 3: 235-240. 10.1007/s12265-010-9164-zView 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
- Chen Y, Olopade OI: MYC in breast tumor progression. Expert Rev Anticancer Ther 2008, 8: 1689-1698. 10.1586/14737126.96.36.1999View ArticlePubMedPubMed CentralGoogle 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
- Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant microRNA targets. Cell 2002, 110: 513-520. 10.1016/S0092-8674(02)00863-2View ArticlePubMedGoogle 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
- O’Neill LA, Sheedy FJ, McCoy CE: MicroRNAs: the fine-tuners of Toll-like receptor signalling. Nat Rev Immunol 2011, 11: 163-175. 10.1038/nri2957View ArticlePubMedGoogle Scholar
- Ashokkumar C, Ningappa M, Ranganathan S, Higgs BW, Sun Q, Schmitt L, Snyder S, Dobberstein J, Branca M, Jaffe R, Zeevi A, Squires R, Alissa F, Shneider B, Soltys K, Bond G, Abu-Elmagd K, Humar A, Mazariegos G, Hakonarson H, Sindhi R: Increased expression of peripheral blood leukocyte genes implicate CD14+ tissue macrophages in cellular intestine allograft rejection. Am J Pathol 2011, 179: 1929-1938. 10.1016/j.ajpath.2011.06.040View ArticlePubMedPubMed CentralGoogle Scholar
- Gabriely G, Wurdinger T, Kesari S, Esau CC, Burchard J, Linsley PS, Krichevsky AM: MicroRNA 21 promotes glioma invasion by targeting matrix metalloproteinase regulators. Mol Cell Biol 2008, 28: 5369-5380. 10.1128/MCB.00479-08View 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
- Saito J, Toriumi S, Awano K, Ichijo H, Sasaki K, Kobayashi T, Tamura S: Regulation of apoptosis signal-regulating kinase 1 by protein phosphatase 2Cepsilon. Biochem J 2007, 405: 591-596. 10.1042/BJ20070231View 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
- Chen CZ, Li L, Lodish HF, Bartel DP: MicroRNAs modulate hematopoietic lineage differentiation. Science 2004, 303: 83-86. 10.1126/science.1091903View ArticlePubMedGoogle 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