ORIGINAL RESEARCH article
Sedentary and Trained Older Men Have Distinct Circulating Exosomal microRNA Profiles at Baseline and in Response to Acute Exercise
- 1Department of Neurology, Center for Advanced Research on Diagnostic Assays, Icahn School of Medicine at Mount Sinai, New York, NY, United States
- 2Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, United States
- 3Translational Research Institute, AdventHealth, Orlando, FL, United States
Exercise has multi-systemic benefits and attenuates the physiological impairments associated with aging. Emerging evidence suggests that circulating exosomes mediate some of the beneficial effects of exercise via the transfer of microRNAs between tissues. However, the impact of regular exercise and acute exercise on circulating exosomal microRNAs (exomiRs) in older populations remains unknown. In the present study, we analyzed circulating exomiR expression in endurance-trained elderly men (n = 5) and age-matched sedentary males (n = 5) at baseline (Pre), immediately after a forty minute bout of aerobic exercise on a cycle ergometer (Post), and three hours after this acute exercise (3hPost). Following the isolation and enrichment of exosomes from plasma, exosome-enriched preparations were characterized and exomiR levels were determined by sequencing. The effect of regular exercise on circulating exomiRs was assessed by comparing the baseline expression levels in the trained and sedentary groups. The effect of acute exercise was determined by comparing baseline and post-training expression levels in each group. Regular exercise resulted in significantly increased baseline expression of three exomiRs (miR-486-5p, miR-215-5p, miR-941) and decreased expression of one exomiR (miR-151b). Acute exercise altered circulating exomiR expression in both groups. However, exomiRs regulated by acute exercise in the trained group (7 miRNAs at Post and 8 at 3hPost) were distinct from those in the sedentary group (9 at Post and 4 at 3hPost). Pathway analysis prediction and reported target validation experiments revealed that the majority of exercise-regulated exomiRs are targeting genes that are related to IGF-1 signaling, a pathway involved in exercise-induced muscle and cardiac hypertrophy. The immediately post-acute exercise exomiR signature in the trained group correlates with activation of IGF-1 signaling, whereas in the sedentary group it is associated with inhibition of IGF-1 signaling. While further validation is needed, including measurements of IGF-1/IGF-1 signaling in blood or skeletal muscle, our results suggest that training status may counteract age-related anabolic resistance by modulating circulating exomiR profiles both at baseline and in response to acute exercise.
Aging is associated with loss of skeletal muscle mass and function, physiological alterations of the cardiorespiratory system, weakening of the immune system, and impaired cognitive function (Garatachea et al., 2015; Weyand and Goronzy, 2016). While a sedentary lifestyle may accelerate the hallmarks of aging, regular exercise can lessen those multisystem age-related declines, thereby decreasing the risk of cardiovascular diseases, cancers, and metabolic disorders (Febbraio, 2017; Harridge and Lazarus, 2017; Howden et al., 2018; Chowdhury et al., 2019). Major long-term adaptations to exercise training include improved cardiovascular and cardiorespiratory fitness, enhancement of skeletal muscle blood flow and oxidative capacities (via an increased density of capillaries and of mitochondria, respectively) (Rivera-Brown and Frontera, 2012), and muscle hypertrophy (Konopka and Harber, 2014). Even a single bout of exercise can induce several health benefits in skeletal muscle, such as enhanced insulin sensitivity, increased secretion of the anti-inflammatory cytokine interleukin-6, and stimulation of immune cells (Garatachea et al., 2015; Sellami et al., 2018).
While the molecular mechanisms underlying the health benefits of exercise remain unclear, accumulating evidence suggests the importance of tissue crosstalk via the release of molecules into the bloodstream (i.e., circulation). Many of these molecules are packaged into extracellular vesicles (EVs) such as microvesicles and exosomes (Whitham et al., 2018). Exosomes are nano-sized (30–150 nm) EVs that transport proteins and various RNA species, notably miRNAs (Kowal et al., 2014). Exosomes are found in most body fluids, including plasma, urine, and saliva. There is growing evidence that exercise leads to the release of exosomes by skeletal muscle and other tissues into circulation, suggesting that these EVs act as mediators of the systemic adaptations to exercise (Chaturvedi et al., 2015; Fruhbeis et al., 2015; Safdar and Tarnopolsky, 2018; Whitham et al., 2018).
Among the molecules transported by exosomes, microRNAs (miRNAs) are small non-coding RNAs that mediate post-transcriptional gene silencing. miRNAs are important regulators of gene expression (Selbach et al., 2008; O’Brien et al., 2018). Extracellular (or circulating) miRNAs consist of two populations: those associated with vesicles such as exosomes, and those bound to proteins like AGO2 (O’Brien et al., 2018). The majority of previous studies on circulating miRNA responses to exercise measured changes in vesicle-free extracellular miRNAs in either plasma or serum (Baggish et al., 2011, 2014; Aoi et al., 2013; Mooren et al., 2014; Nielsen et al., 2014; Russell and Lamon, 2015; Xu et al., 2015; D’Souza et al., 2017; Shah et al., 2017; Horak et al., 2018; Barber et al., 2019). Those studies have suggested that circulating miRNAs may mediate transient metabolic changes and adaptive responses induced by exercise. However, recent reports involving healthy young men suggest that the exosomal fraction is distinctively modulated by exercise and may play an important role in exercise-induced systemic adaptations (Guescini et al., 2015; D’Souza et al., 2018; Lovett et al., 2018; Hou et al., 2019).
Aging has been associated with alterations in intramuscular and circulating miRNA profiles (Drummond et al., 2011; Noren Hooten et al., 2013; Lai et al., 2014; Olivieri et al., 2014b; Zhang et al., 2015). However, the effect of exercise on the circulating exosomal miRNA (exomiR) expression profiles of aging individuals has not been investigated. Here, we studied the impact of long-term exercise training and acute endurance exercise on circulating exosomal miRNAs using a comprehensive small RNA sequencing-based approach. The results indicate that long-term exercise affects basal (resting) circulating exomiR profiles and demonstrate that acute exercise generates distinct patterns of exomiR profiles specific to trained and sedentary groups.
Materials and Methods
Older men (>65 years old) were recruited from the Orlando, FL, United States area. Participants were in good general health, defined as having no chronic medical conditions (i.e., type I or II diabetes), no contraindications for exercise, stable weight for the last 6 months, and controlled hypertension (<150 mmHg systolic, <90 mmHg diastolic). Participants (n = 10) were assigned to one of two groups based on their physical activity levels: 5 older trained and 5 older sedentary. Subject characteristics are summarized in Table 1. Trained subjects were engaged in endurance activities (running, cycling) for at least 5 years without an extensive layoff (>6 months) due to injury. Sedentary individuals completed ≤1 continuous exercise bout per week. Each participant provided written informed consent before completing any data collection procedures. The study was approved by the Florida Hospital, AdventHealth Orlando (Orlando, FL, United States) Institutional Review Board (IRBNet #554559) and performed in compliance with the Declaration of Helsinki.
All testing was completed at the Translational Research Institute at AdventHealth, Orlando. Visit #1 consisted of a fasting blood draw, physical measurements, medical history/physical activity questionnaires, and resting electrocardiography (ECG). On visit #2, participants completed a VO2 max test with ECG to determine cardiorespiratory fitness. On visit #3, participants completed a dual energy X-ray absorptiometry (DXA) scan and were given a triaxial accelerometer to determine free living physical activity over a 1-week period. On visit #4, fasting blood samples were drawn into K2-EDTA tubes (pre-exercise), and participants then consumed a small low glycemic index meal (200 kcal, 15% protein, 35% fat, and 50% carbohydrate). Participants then completed the acute exercise bout on a cycle ergometer, as described in the next sections.
Cardiorespiratory Testing and Assessment of Body Composition
VO2 max was determined as peak aerobic capacity measured during a graded exercise protocol on an electronically braked cycle ergometer, as previously described (Coen et al., 2013; Dube et al., 2014). Heart rate, blood pressure, and ECG (12-lead) were recorded throughout this test. The test was terminated according to the criteria outlined in the American College of Sports Medicine exercise testing guidelines (Thompson et al., 2010). DXA scans were performed using a GE Lunar iDXA whole-body scanner and analyzed with enCORE software. Coefficients of variation (CV) vary by tissue type and range from 0.5% to 1.1%. This CV data was provided by the manufacturer (for reference, see Faulkner et al., 2006).
Daily Physical Activity
The monitor used for this study was the SenseWear® Pro Armband (BodyMedia Inc., Pittsburgh, PA, United States). This activity monitor uses accelerometry, galvanic skin response, and skin temperature to estimate energy expenditure at a 1-min resolution (Berntsen et al., 2010). The SenseWear® Armband has been demonstrated to be a valid, accurate, and reliable method when compared with indirect calorimetry (Berntsen et al., 2010) and doubly labeled water (Mackey et al., 2011). SenseWear professional version 8.1 software was used for data analysis. Participants were instructed to wear the armband on their upper left arm at all times over the seven-day period, with the exception of showering/bathing. Only days with a wear time of at least 85% were considered for further analysis.
Acute Exercise Testing
Study participants performed an acute bout of aerobic exercise on a cycle ergometer as previously described (Distefano and Goodpaster, 2018). Following a 6-min warm-up, participants cycled for 40 min at ∼70% heart rate reserve, which was calculated during the maximum aerobic capacity test. Heart rate, perceived exertion, and blood pressure were measured every 5 min. Indirect calorimetry was examined with an open-circuit spirometry metabolic monitoring system (Parvo Medics, Sandy, UT), and used to verify exercise intensity. Blood was drawn before exercise (Pre), immediately after exercise (Post), and after 3 h of recovery (3hPost) via direct venipuncture using a 21-gauge butterfly needle, and blood samples were collected into PAXgene blood tubes with spray-coated K2-EDTA (Becton-Dickson, Catalog #366643).
Collected blood samples were processed according to the manufacturer’s instructions. Briefly, they were centrifuged at 3,000 ×g for 10 min at room temperature. Plasma was aliquoted and immediately stored at −80°C.
Isolation and Enrichment of Exosomes From Plasma
Plasma was thawed on ice and centrifuged at 3,000 × g for 15 min at 4°C to remove particulate matter. Pre-cleared plasma (250 μl) was then mixed with an equal volume of thromboplastin D and incubated for 15 min at 37°C. Supernatants were collected following centrifugation at 13,000 × g for 10 min at 4°C, mixed with 100 μl of ExoQuick (System Biosciences, Palo Alto, CA, United States) and 6.25 μl of RNase (1 mg/ml), and incubated overnight at 4°C. Exosome-enriched preparations were spun down at 1,500 × g for 25 min at 4°C. Supernatants were then gently removed, and exosome-enriched pellets were resuspended in 50 μl of PBS with 1.25 μl of RNAse inhibitor (ThermoFisher Scientific, Waltham, MA, United States) and stored at −20°C.
Nanoparticle Tracking Analysis
Exosome-enriched preparations were measured and quantified using a NanoSight NS300 instrument (Malvern Panalytical Ltd, Malvern, United Kingdom) equipped with a 532 nm laser and running NTA v3.2 software. These preparations were thawed on ice and gently resuspended immediately before sampling to prevent the settling of particles. Immediately before measurement, the exosome-enriched aggregates were broken up by gentle passage through a needle. Samples were diluted empirically with 0.22 mm filtered PBS to achieve 20–60 particles/frame (manufacturer’s optimum measurement range) using 1.5 ml low-adhesion tubes (United States Scientific, Orlando, FL, United States) to prevent exosome adsorption onto the walls of the tube. Diluted samples were loaded into the assembled sample chamber of a NanoSight NS300 instrument through a 1-ml syringe, applying each sample slowly to expel air without bubble formation. The chamber was then connected to the instrument and laser-engaged, and microparticles were brought into focus using the thumb print region as a reference; three 60-second video images/sample were acquired and analyzed using NanoSight NTA 2.3 software. Individual members of the prepared serial dilutions were analyzed until the detected raw concentration was within the recommended range for the instrument (1 × 108 – 1 × 109 particles/ml).
RNA Isolation From Exosome-Enriched Preparations
For RNA isolation from exosome-enriched preparations, 700 μl of QIAzol lysis reagent (Qiagen, Germantown, MD, United States) were added to 50 μl of exosome preparation, and the mixture was incubated for 5 min at room temperature. After the addition of 140 μl of chloroform, the mixture was vortexed for 15 s and incubated for 3 min at room temperature. Following incubation, the mixture was centrifuged at 12,000 × g for 15 min at 4°C. After the upper aqueous phase was transferred to a fresh tube, a two-fold volume of absolute ethanol (Sigma-Aldrich, St. Louis, MO, United States) was added, and the mixture was vortexed briefly. Samples were then transferred to RNeasy Mini spin columns (Qiagen) that were centrifuged at 10,000 × g for 30 s at 4°C, and flow-through was discarded. An 80-μl volume of DNase incubation mix was added directly onto the column membrane and samples were incubated for 15 min at room temperature. Following wash steps with RWT buffer, RPE buffer (Qiagen), and 85% ethanol, RNA pellets were dried by centrifugation at 10,000 × g for 2 min at room temperature and eluted in 15 μl of water. Electrophoretic analysis of the RNA preparations was conducted on an Agilent 2100 Bioanalyzer using the Agilent Small RNA kit (Agilent, Santa Clara, CA, United States).
Small RNA Library Preparation
Small RNA libraries were prepared using the NEBNext Small RNA Library Prep Set for Illumina (New England Biolabs, Ipswich, MA, United States), according to the manufacturer’s instructions. Briefly, RNA-Seq libraries were constructed using 2 ng of RNA isolated from exosome-enriched preparations and 15 cycles of PCR amplification. Libraries were run on a 6% TBE gel for 60 min at 120 V, and bands between 140 to 160 bp were excised from the gel using a razor blade. Gel pieces were crushed into smaller pieces and incubated in water overnight on a rotator. Libraries were isolated using gel filtration columns. Following ethanol precipitation, the resulting pellets were each resuspended in 11 μl of ultrapure water.
Small RNA Library Sequencing and Data Analysis
All libraries were analyzed on an Agilent 2100 Bioanalyzer using the Agilent High Sensitivity DNA kit. For sequencing, libraries were pooled at equimolar concentrations, mixed with 5% PhiX (Illumina, San Diego, CA, United States), and single-end sequencing was performed on an Illumina MiSeq platform following the manufacturer’s recommendations. Sequencing data have been deposited to GEO (GSE144627).
Small RNA-seq reads were processed using the exceRpt pipeline (Rozowsky et al., 2019), as recommended by the ERC Consortium1. The pipeline initially trims the reads using FASTX-Toolkit to remove remnants of the 3′-adapter sequence followed by removal of the reads of low-quality and homopolymer repeats, and subsequently aligns the trimmed reads to NCBI’s Univec database to filter out contaminants. The good reads were mapped to rRNA, and then to the human genome and transcriptome to obtain the counts of different gencode and RNA transcripts including miRNA and tRNA. The order of annotations is as follows: miRNAs (miRBase v21), tRNAs, piRNAs, gencode transcripts (hg38 Gencode v24; snRNAs, snoRNAs, miscRNAs, protein coding mRNAs, and lncRNAs), and circular RNAs. Mapping of the reads to the NCBI’s Univec database or rRNA or human genome/transcriptome was performed using STAR 2.4.2a (Dobin et al., 2013).
Subject characteristics, physiological and exomiR data were reported as the mean ± standard error of the mean (SEM). For subject characteristics and physiological data, paired variables were compared using a Student’s t-test or a Wilcoxon signed-rank test as appropriate for the data distribution. When comparing miRNA expression changes either between time points (within the same group) or between groups (baseline), the Bioconductor DESeq2 package was used (Love et al., 2014). A log2 fold change ≥ 1 in either direction together with a Wald test p-value ≤ 0.05 indicated a statistically significant difference. Differentially expressed miRNAs were displayed in volcano plots and a heatmap. The computation was done in R 3.6.1 (R Core Team, 2019) and Bioconductor 3.10 (Gentleman et al., 2004) packages.
Anthropometric and Physiological Characteristics of Participants
The study cohort (n = 10) comprised 5 older endurance-trained and 5 age-matched sedentary subjects, all of whom completed a 40 min bout of cycling exercise. Subject characteristics are provided in Table 1. The study groups were matched for age (all were male) and trained participants had significantly higher cardiorespiratory fitness (VO2 max), compared to the sedentary group, as per the study design. The sedentary group tended to have greater total fat mass and weight compared to the trained group. In contrast, trained and sedentary groups had comparable skeletal muscle index (SMI), and total leg lean mass.
Physiological characteristics measured during the acute exercise bout are presented in Table 2. While steady state VO2 and power output were significantly higher in the trained group, the participants performed the exercise bout at a similar relative intensity (% heart rate reserve (HRR) and % VO2 max). Additionally, the average respiratory quotient (RQ), an index of substrate metabolic utilization, was similar among the groups.
Characterization of Exosome-Enriched Preparations
We sought to validate our methodology for the isolation and enrichment of exosomes and for RNA extraction from the exosome-enriched preparations. Exosome-enriched preparations were isolated from blood plasma samples obtained from study participants, as described in the Materials and Methods. Analysis of exosome-enriched preparations using nanoparticle tracking analysis (NTA) allowed us to obtain the size distribution of EVs and an estimate of particle concentration with high resolution (Dragovic et al., 2011). Size distribution ranged from 68 to 105 nm, with a mean size of 85.1 nm (Figures 1A,B), which falls into the size range of exosomes (Jayaseelan, 2019), suggesting that the highly purified EVs contained mostly exosomes. Exosome concentrations in our preparations ranged from 5 × 108 to 1 × 109 particles/ml. We next examined RNA composition in the exosome-enriched preparations. Characterization of the RNA from exosome-enriched preparations by Bioanalyzer analysis revealed that the majority of RNA species (60–80%) were approximately 22 nt in size (Figures 1C,D). Those small RNA preparations were not degraded by DNase. However, they were completely degraded by RNase treatment (Figures 1C,D), thus demonstrating the integrity of small RNAs present in the exosome-enriched preparations.
Figure 1. Characterization of exosome-enriched preparations isolated from plasma and analysis of their small RNA contents. (A) Analysis of EV size distribution and concentration in exosome-enriched preparations using a NanoSight NS300 instrument. (B) Representative image from the NanoSight instrument showing optimal light scatter from plasma-derived nanovesicles. (C,D) Bioanalyzer traces of the small RNAs extracted from exosome-enriched preparations following incubation with DNase or RNase. Lane 1, Control (non-treated); lane 2, DNase-treated; lane 3, RNase-treated. (E) Sequencing analysis of the small RNAs extracted from exosome-enriched preparations. Blood was drawn from each study participant at 3 time points: before exercise (Pre; open bars), immediately after exercise (Post; diagonal line bars), and 3 h after exercise (3 hPost; vertical line bars). Participants comprised groups of older trained and sedentary adults, as indicated. Shown at the top are the total number of reads, the percentage of mappable reads after adapter clipping, and the percentage of reads mapping to the human genome/transcriptome. Stacked-bar plots represent the percentage of mapped reads assigned to miRNA, rRNA, tRNA, other gencode transcripts (which are not miRNA, rRNA or tRNA), and other mapped sequences (mapped to the genome, but not to the four transcript types mentioned above). The different RNA species are color-coded. Shown at the bottom is the percentage of transcriptome complexity, defined as the ratio of the number of unique transcript sequences to the total number of transcript sequences.
Characterization of Small RNAs Isolated From Exosome-Enriched Preparations
We performed sequencing of the small RNA libraries to characterize the circulating exosomal miRNAs from older trained and sedentary subjects and study their dynamic regulation following exercise. In our sequencing data, the percentage of mappable reads after adapter clipping ranged from 60 to over 90%. About 70–90% of total reads were mapped to the human genome and the transcriptome complexity was homogeneous across all samples. The QC of the sequencing data did not show bias towards either a group (trained vs. sedentary) or a time point. Overall, our small RNA sequencing data (Figure 1E) passed the raw read data QC when compared with the ERC Consortium QC standards for small RNA-seq data (see2).
Using mapped sequence read counts, we quantified the RNA contents of exosome-enriched preparations. As illustrated in Figure 1E, exosome-enriched preparations isolated from plasma contained a very diverse RNA ‘cargo’. As expected from previous research (Li et al., 2014; Perez-Boza et al., 2018), a large proportion of miRNAs, tRNAs, and rRNAs were detected, as well as a small fraction of other gencode transcripts including snRNAs, snoRNAs, miscRNAs, protein coding mRNAs, and lncRNAs. MicroRNAs were among the topmost abundant RNA species, accounting for 28.69 ± 8.76% of all mappable reads for all samples (all subjects and all time points). At the current read depth, 371 ± 71 miRNAs were detected. The five most common miRNAs (miR-99a-5p, miR-129-5p, miR-128-3p, let-7b-5p, miR-423-5p) collectively accounted for approximately 50% of all mappable miRNA sequences of all samples together. One of the most abundant miRNAs in red blood cells, miR-16 (Kirschner et al., 2011), was hardly detectable amongst all samples (about 2.63%), thus corroborating the lack of hemolysis during sample preparation. Altogether, our characterization of RNA content in exosome-enriched preparations is in agreement with previous reports (Li et al., 2014; Yeri et al., 2017) and supports the validity of our experimental approach for analyzing miRNAs in exosomes.
Older Trained and Sedentary Subjects Show Different Circulating ExomiR Profiles at Baseline
To assess similarities and differences in circulating exomiR profiles in the trained and sedentary groups, we employed a hierarchical clustering method comparing circulating exomiR profiles at baseline (Pre), immediately after acute exercise (Post), and 3 h after acute exercise (3hPost). We initially identified the miRNA species that were differentially expressed (i) at Pre between groups, (ii) between time points (Post vs. Pre, 3hPost vs. Pre, and 3hPost vs. Post) within each group. We then applied hierarchical clustering to these thirty-two differentially expressed miRNAs. As shown in Figure 2A, miRNA clusters in trained vs. sedentary subjects differed dramatically at baseline. For instance, miR-486-5p and miR-215-5p showed higher levels of basal expression in trained individuals than in the sedentary.
Figure 2. Differential expression of circulating exomiRs at baseline in the trained and sedentary groups. (A) Heat map with dendrogram illustrating clustering results for the 32 exomiRs identified as differentially expressed when comparing miRNA expression changes between groups (sedentary vs. trained) at baseline (Pre) and between time points (Post vs. Pre and 3hPost vs. Pre) within each group. Hierarchical clustering was applied using the average linkage method and Euclidean distance. Values signify the log2-transformed cpm (counts per million reads of the library size) normalized within each row such that the mean is zero and the standard deviation is one. Colors represent the level of miRNA expression; red: high expression; green: low expression. Groups numbered from 1 to 4 on the right highlight exomiRs that were differentially expressed at baseline (1 and 2) and at Post (4) between sedentary and trained subjects, and those differentially expressed at either Post or 3 hPost relative to Pre in trained subjects (3). (B) Differential expression of circulating exomiRs in the trained relative to sedentary group at baseline (Pre). The p-value in log10 scale of each miRNA is plotted against its fold change in log2 scale and each circle denotes a miRNA. The miRNAs with p-values ≤ 0.05 and log2 (fold change) ≥ 1 in either direction is represented by either red circles (upregulated) or green circles (downregulated). Three circulating exomiRs were upregulated and one was downregulated in the trained compared to sedentary group.
Comparison of the circulating exomiR profiles in trained and sedentary subjects at baseline revealed that four miRNA species were differentially expressed in the trained group (Figure 2B). Of those, three miRNAs (miR-486-5p, miR-215-5p, and miR-941) were upregulated, and one (miR-151b) was downregulated. Differential expression of muscle-enriched miR-486-5p was consistent with various studies reporting significant changes in miR-486-5p levels following exercise (Aoi et al., 2013; Denham and Prestes, 2016; D’Souza et al., 2017, 2018; Barber et al., 2019). Overall, these distinct patterns of miRNA expression indicate that regular exercise may be associated with significant and durable expression changes in resting plasma exomiRs.
Acute Exercise Regulates Distinct Circulating ExomiRs in Trained and Sedentary Subjects
In order to identify acute exercise-regulated miRNAs, we next compared post-exercise profiles (either Post or 3hPost) with the resting profile (Pre) within each group. In the trained group, we identified seven differentially regulated miRNAs at Post relative to Pre (Figures 2A, 3A). Among those, three miRNAs were upregulated (miR-383-5p, miR-339-5p, and miR-874-3p) and four were downregulated (miR-206, miR-486-5p, miR-148a-3p, and let-7b-5p). At 3hPost, eight exomiRs were differentially regulated (Figures 2A, 3B), five of which were upregulated (miR-34b-3p, miR-129-2-3p, miR-138-1-3p, miR-671-3p, and miR-885-5p), and three downregulated (miR-486-5p, miR-629-5p, and miR-16-2-3p). Remarkably, miR-486-5p was the only miRNA species to be shared by both post-exercise signatures. The levels of miR-486-5p in the trained group were significantly decreased at both Post and 3hPost, whereas at baseline miR-486-5p expression was higher than in the untrained (Figure 2B).
Figure 3. Differential regulation of circulating exomiRs following acute exercise in the trained and sedentary groups. (A) Differentially regulated circulating exomiRs immediately after acute exercise and 3 h post-exercise (B) compared to baseline (Pre) in trained group. (C) Differentially regulated circulating exomiRs immediately after exercise and 3 h post-exercise (D) compared to baseline in sedentary group. The p-value in log10 scale of each miRNA is plotted against its fold change in log2 scale and each circle symbolizes a miRNA. The miRNAs with p-values ≤ 0.05 and log2 (fold change) ≥ 1 in either direction are represented by either red circles (upregulated) or green circles (downregulated). In (A) three miRNAs were upregulated and four were downregulated. In (B) five miRNAs were upregulated and three were downregulated. In (C) eight miRNAs were upregulated and one was downregulated. In (D) all four miRNAs were downregulated.
Comparison of the Post and Pre exomiR profiles in the sedentary group identified nine differentially regulated exomiRs (Figures 2A, 3C). Except for miR-4433b-3p, all miRNAs (miR-505-3p, miR-29b-3p, miR-203a-3p, miR-384, miR-451a, miR-223-3p, miR-218-5p, and miR-495-3p) were upregulated in response to acute exercise. By contrast, all differentially regulated exomiRs at 3hPost (miR-4433b-3p, mir-378c, miR-151b, miR-151a-5p| hsa-miR-151b) were downregulated (Figures 2A, 3D). Except for miR-4433b-3p, there was no overlap between the two post-exercise signatures. Of note, miR-151b basal expression level was significantly greater in the sedentary compared to the trained group (Figures 2A, 3D). Overall, our data show that circulating exomiRs modulated by acute exercise in the trained and sedentary groups are totally distinct (Figure 4), and suggest that long-term exercise has a significant impact on the regulation of plasma exomiRs in response to acute exercise.
Figure 4. Acute exercise regulates distinct circulating exomiRs in the trained and sedentary groups. (A) Venn diagram depicting the exomiRs that were differentially regulated at Post and at 3 hPost in the trained group (7 and 8, respectively) and in the sedentary (9 and 4, respectively). There is no overlap between the trained and sedentary groups at either time point. (B) and (C) Shown are the expression levels of miR-383-5p (B) and miR-505-3p (C) in the circulating exosomes of sedentary and endurance-trained subjects before exercise (Pre), immediately after exercise (Post), and 3 h after exercise (3 hPost). Bar plots represent log2-transformed normalized read counts from small RNA sequencing.
Correlation of Exercise-Regulated ExomiRs With IGF1 Signaling
We employed pathway enrichment analysis (IPA), performed miRNA target prediction via miRDB (Wong and Wang, 2015), and mined the literature for experimentally validated miRNA targets to determine which molecular pathways would likely to be altered by the ExomiR signatures at baseline and in response to acute exercise. We found that the majority of differentially regulated miRNAs either target or are predicted to target components of the IGF1 signaling pathway. The IGF1 pathway, which signals through AKT-mediated activation of mTOR and repression of FOXO, and via the Ras/MAP kinase pathway, is a major regulator of muscle growth and glucose homeostasis (Schiaffino and Mammucari, 2011; O’Neill et al., 2015). IGF1 signaling is also implicated in exercise-induced muscle and cardiac hypertrophy (Velloso, 2008; Weeks and McMullen, 2011).
Three of the four miRNAs in the trained subjects’ baseline exomiR signature (miR-486-5p, miR-215-5p, and miR-941) were previously shown to target components of the IGF1 signaling pathway (Pichiorri et al., 2010; Hu et al., 2012; Peng et al., 2013; Figure 5A). Additionally, downregulated miR-151b is predicted to target RALGAPA1, an inhibitor of Ras signaling (Yan and Theodorescu, 2018). Thus, the data suggest the hypothesis that resting exomiRs differentially modulate IGF-1 signaling in trained vs. untrained elderly subjects.
Figure 5. Association of circulating exomiR signatures with IGF1 signaling. Shown are the circulating exomiRs that were differentially regulated at Post, along with their reported/predicted target genes linked to IGF1 signaling, when applicable. Red border-rectangles indicate miRNAs that are upregulated; blue border-rectangles indicate miRNAs that are downregulated. Colors represent the level of miRNA expression; red: high expression; green: low expression. Up/down arrows indicate whether a target gene is up- or downregulated, respectively. (A) Circulating exomiRs that are differentially expressed in exercise-trained individuals relative to the sedentary lead to inhibition of IGF1 signaling. (B) In trained individuals, acute exercise-regulated exomiRs increase tumor suppressor activity and activate IGF1 signaling. (C) In sedentary individuals, the majority of acute exercise-regulated exomiRs are upregulated and result in downregulation of IGF1 signaling. Note that downregulated miR-4433b-3p is not shown because its target is undetermined.
In the trained group’s Post vs. Pre signature, all three upregulated miRNAs (miR-383-5p, miR-339-5p, and miR-874-3p) have a tumor suppressive function (Zhang et al., 2014; Zhao et al., 2017; Jiang D. et al., 2019; Jiang J et al.,2019b; Figure 5B). In contrast, all four downregulated miRNA species target components of IGF1 signaling (Peng et al., 2013; Xu et al., 2013; Gao et al., 2014; Xu et al., 2014; Liu et al., 2018), suggesting that the exercise-induced exomiR signature in trained older adults activates IGF1 signaling. Except for miR-34b-3p (Hermeking, 2010), the majority of differentially regulated miRNAs at 3hPost vs. Pre target components of IGF1 signaling (Peng et al., 2013; Ye et al., 2015; Meier et al., 2016; Ntoumou et al., 2017; Yan et al., 2017; Zhu et al., 2018; Han et al., 2019), suggesting that they may also modulate the IGF1 pathway.
In contrast with the trained group, all but one miRNA species were upregulated in the sedentary’s Post vs. Pre signature. Downregulated miR-4433b-3p has no established or predicted gene target related to IGF1 signaling. The majority of increased miRNAs (seven out of eight) repress components of IGF1 signaling (Uesugi et al., 2011; Olivieri et al., 2014a; Yuan et al., 2015; Kim et al., 2016; Lai et al., 2016; Minna et al., 2016; Li et al., 2017; Chen et al., 2018; Wang Z. et al., 2018; Tang et al., 2019; Zeng et al., 2019), suggesting that the exercise-induced exomiR signature in the sedentary inhibits IGF1 signaling (Figure 5C). The majority of differentially regulated miRNAs at 3hPost vs. Pre either target or are predicted to target components of IGF1 signaling (Knezevic et al., 2012), implying a modulation of this pathway. Collectively, our findings suggest that exercise-regulated plasma exomiRs modulate IGF1 signaling. Particularly, the plasma exomiR signatures induced immediately after exercise may have opposing effects on IGF1 signaling in trained vs. sedentary aging adults.
In the present study, we examined for the first time the effects of chronic and acute endurance exercise on global circulating exomiR profiles in older adults, using a next-generation sequencing-based approach. We identified a subset of differentially expressed plasma-derived exomiRs that target genes the IGF-1 signaling pathway in a long-term exercise (trained) compared to an untrained (sedentary) older subjects. We also found that acute exercise-induced plasma exomiR signatures in trained and untrained groups are distinct and have opposing effects on IGF-1 signaling. IGF-1 signaling is an important regulator of muscle growth and glucose homeostasis (Schiaffino and Mammucari, 2011; O’Neill et al., 2015; Vassilakos et al., 2019). In light of the known relationship between exercise and IGF-1 signaling (Velloso, 2008) and the dysregulation of IGF-1 signaling associated with aging (Rebelo-Marques et al., 2018), exercise-regulated circulating exomiRs might contribute to a differential regulation of IGF-1 signaling in trained compared to sedentary elderly subjects, suggesting their role as mediators of exercise-induced adaptations in the elderly population. Thus, our data could provide novel insights into the molecular mechanisms underlying the physiology of exercise in the elderly population.
Among the species regulated with exercise, we identified two muscle-enriched miRNAs (Horak et al., 2016), miR-486-5p and miR-206, whose exosomal levels were downregulated in trained individuals immediately after a single bout of acute exercise. Additionally, miR-486 baseline level was greater in the older trained group than in the sedentary group (Figure 2A). Previous studies have shown that miR-486 was more abundant in resting athletes’ whole blood than in healthy controls, whereas it was decreased following acute exercise in healthy controls (Denham and Prestes, 2016). In skeletal muscle, miR-206 was unaltered following acute exercise, but significantly decreased following 12 weeks of endurance training (Nielsen et al., 2010). Validated targets of miR-486 include both positive (IGF-1, IGF-1R, PIK3R1) and negative (PTEN) components of IGF-1 signaling (Small et al., 2010; Peng et al., 2013), suggesting that miR-486 may have both activating and inhibiting effects on PI3K/AKT signaling. Because activation of IGF-1/PI3K/AKT signaling is recognized to promote skeletal and cardiac muscle growth (Antos et al., 2002; Stitt et al., 2004; Vassilakos and Barton, 2018), we anticipate miR-486 to be involved in the regulation of muscle growth. MiR-206, which inhibits the PI3K/AKT signaling pathway in rat (Liu et al., 2018), was previously shown to promote skeletal muscle differentiation (Kim et al., 2006). Overall, the presence of these muscle enriched miRNAs in post-exercise signature is consistent with the impact of exercise on skeletal muscle hypertrophy and regeneration in response to long-term regular exercise (Bazgir et al., 2017).
In addition to muscle-enriched miRNAs, we found various non-muscle-specific miRNAs in the plasma exomiR profiles of elderly subjects following exercise. Some of these miRNA species, e.g., miR-215-5p, miR-941, miR-671-3p, miR-203a-3p, miR-4433b, and miR-384, have not been linked to exercise thus far, while others, like miR-29b, miR-451, miR-223, miR-505, and miR-629 have. Modulation of miR-29b and miR-451 following exercise training was formerly reported in skeletal muscle (Keller et al., 2011; Russell et al., 2013), and miR-29b was shown to promote muscle atrophy via the targeting of IGF-1 and PI3K (Li et al., 2017). MiR-451a also targets the AKT/mTOR pathway (Minna et al., 2016; Olioso et al., 2019), although its inhibitory effect on myogenesis in mice is not AKT-mediated (Munk et al., 2019). Our observations of exercise-induced alterations in exosomal miR-29b and miR-451 are consistent with previous studies and suggest that they could be involved in the regulation of skeletal muscle adaptations to exercise. Remarkably, upregulation of miR-29b following exercise was previously demonstrated in heart tissue-derived exosomes, suggesting a potential role in cardiac remodeling (Chaturvedi et al., 2015). Exercise-induced alterations in miR-223, miR-505, miR-629, and miR-486 expression were reported in immune cells (blood neutrophils, natural killers, and PBMCs), supporting their role as potential mediators of the effects of acute exercise on the immune system and inflammatory mechanisms (Radom-Aizik et al., 2010, 2012, 2013; Simpson et al., 2012; Nieman and Wentz, 2019). Both miR-223 and miR-629 were shown to target FOXO3A, a key mediator of the PI3K/AKT signaling pathway (Olivieri et al., 2014a; Kim et al., 2016; Yan et al., 2017), while miR-505 was recently shown to inhibit AKT phosphorylation (Tang et al., 2019). Collectively, our data support the concept that circulating exomiRs may be released by various tissues besides skeletal muscle, thereby contributing to the multisystemic benefits of exercise (for review, see Safdar and Tarnopolsky, 2018).
There is a paucity of data on the regulation of circulating exomiR expression by exercise (Guescini et al., 2015; D’Souza et al., 2018; Lovett et al., 2018; Hou et al., 2019). Inversely, numerous studies have reported the modulation of cell-free circulating miRNAs following exercise (Baggish et al., 2011; Aoi et al., 2013; Mooren et al., 2014; Nielsen et al., 2014; D’Souza et al., 2017; Shah et al., 2017; Horak et al., 2018; Barber et al., 2019). In the present work, older trained individuals displayed significant changes in exosomal miR-383-5p, miR-874-3p, miR-206, miR-486, miR-148, and let-7b levels immediately after exercise. Previous studies on circulating miRNAs reported significant expression changes in these six miRNA species following acute or chronic aerobic exercise (Nielsen et al., 2010; Aoi et al., 2013; Mooren et al., 2014; Nielsen et al., 2014; Denham and Prestes, 2016; D’Souza et al., 2017, 2018; Barber et al., 2019; Faraldi et al., 2019; Tung et al., 2019), thus strengthening our findings. Comparison of our exercise-induced plasma exomiR alterations with those from previous studies on circulating exosomes revealed both similarities and differences. Similar to the D’Souza study (D’Souza et al., 2018), we found an increase in miR-451a, a miRNA expressed in various tissues (Ludwig et al., 2016), in sedentary individuals following exercise. However, while we observed a decrease in muscle-specific miR-486-5p in trained individuals post-exercise (at both Post and 3hPost), other groups reported either an increase (D’Souza et al., 2018) or no significant change in untrained subjects (Lovett et al., 2018). Contrary to a recent study (Hou et al., 2019), we saw no significant change in baseline miR-342-5p in trained individuals. Since all of those previous studies involved young men, differences in the age of the cohort and in exercise regimens may likely contribute to these variations.
The mechanisms involved in exosome formation and secretion are poorly understood (Stoorvogel, 2012). Exercise-induced downregulation of miR-486 and miR-206 in older trained individuals might reflect an intake by cardiac or skeletal muscle cells, as these species are enriched in cardiac and skeletal muscle (Small et al., 2010; Kirby and McCarthy, 2013). By contrast, the higher miR-486 levels observed in the trained group at baseline could reflect an increase in muscle miRNA export, and/or either an increase in miRNA biogenesis or a decrease in its degradation. Acute exercise was previously shown to cause a transient increase in the level of circulating exosomes (Fruhbeis et al., 2015; Lovett et al., 2018; Whitham et al., 2018). Although variation in particle number or size was not evaluated in this study, we anticipate that there be a similar increase in EV levels immediately after exercise. Additionally, it would be important to determine whether the increase in exomiRs is due to an increase in a subpopulation of exosomes containing these specific miRNAs, or/and to increased packaging of these miRNAs. New technologies are being developed to analyze individual vesicles and identify heterogeneous EV populations (in terms of size and composition) and EV subtypes (Hartjes et al., 2019). Thus, further investigation is needed to elucidate the molecular basis of miRNA transfer by exosomes and determine the kinetics of miRNA turnover and the functional roles of differentially regulated exomiRs.
The differential expression of resting plasma-derived exomiR in older trained relative to the sedentary subjects suggests that exercise training exerts long-term alterations on circulating exomiR profiles. This idea is supported by a recent study by Hou et al. demonstrating increased levels of a circulating exomiR (miR-342-5p) in young athletes compared with untrained students (Hou et al., 2019). Similarly, the levels of circulating miRNAs and intramuscular muscle-specific miRNAs are altered following exercise training (Nielsen et al., 2010, 2014; Baggish et al., 2011; Aoi et al., 2013; Horak et al., 2016; Barber et al., 2019), and pronounced intramuscular gene expression changes have been described following two months of endurance training (Schmutz et al., 2006; Popov et al., 2019). Hence, in older adults, circulating exomiRs that are induced by exercise training could mediate some of the long-term physiological adaptations in muscle, heart, and other tissues (Rivera-Brown and Frontera, 2012), conceivably via regulation of IGF-1 signaling.
Exercise stimulates IGF-1 synthesis directly in muscle, and indirectly in most tissues via the stimulation of growth hormone secretion (Velloso, 2008). In muscle, both exercise and IGF-1 activate the AKT/mTOR signaling pathway, which stimulates protein synthesis and results in muscle hypertrophy. In the heart, exercise-induced cardiac hypertrophy is dependent on IGF-1 signaling (Troncoso et al., 2014; Wang L. et al., 2018). The divergent effects of plasma exomiR signatures on IGF-1 signaling (both at pre and post) in trained compared to untrained older subjects echo previous reports showing: (i) a higher anabolic response (greater muscle mass) and higher insulin sensitivity in endurance-trained old individuals compared to the untrained (Mikkelsen et al., 2013), and (ii) an activation of IGF-1 signaling in senescent rats subjected to acceleration training (based on muscle protein quantification) (Launay et al., 2017). As aging is associated with reduced IGF-1 levels and decreased IGF-1 signaling (Rebelo-Marques et al., 2018), it is tempting to speculate that regular exercise induces plasma exomiR changes that may help to counteract age-related anabolic resistance. Further investigation will be needed to ascertain this hypothesis, including measurements of IGF-1 levels in blood and/or IGF-1 signaling components in skeletal muscle biopsies, and an evaluation of the effect of regular exercise in a group of untrained elderly subjects (e.g., before and after completing 20 weeks of endurance training).
Identification of several regulated miRNAs that have been reported in other studies provides support for the validity and reproducibility of our findings. However, we note that our study is limited in cohort size, and that replication and extension of our findings in a larger cohort would be valuable. We acknowledge that without a non-exercise control group the study design does not control for the effect of dietary status. Furthermore, our study is limited to males. Given the differences between men and women with regard to the physiological response to exercise (Sheel, 2016), it would be interesting to evaluate the impact of exercise training and acute exercise on plasma exomiRs in older females as well.
We also note methodological limitations in our study. Our EV characterization lacks the characterization of EV protein contents (e.g., by Western blot) as well as the visualization of single EVs using high resolution electron microscopy (e.g., by transmission electron microscopy) (for reference, see Thery et al., 2018). While the ExoQuick precipitation-based method used for exosome isolation we used elicits a lower protein purity compared to other methods (Tang et al., 2017; Macias et al., 2019), it shows a higher particle recovery, which is most appropriate for downstream analysis exosomal miRNAs. We used a low-bias protocol for small RNA library preparation. However, some sequence bias may still occur during adapter ligation (Baran-Gale et al., 2015; Dard-Dascot et al., 2018). We did not confirm our RNA-seq-based results with qPCR measurement of regulated transcripts, although others have reported good correspondence of miRNA levels measured by sequencing and by qPCR (Mitchell et al., 2018).
Beyond the potential utility of plasma-derived exomiRs as biomarkers of physical fitness and exercise physiology (e.g., muscle contractile activity), our work suggests their important role in the regulation of IGF-1 signaling in other tissues. Potential therapeutic applications include the assessment of exercise capacity in specific groups (e.g., overweight people), and predicting the risk for cardiovascular diseases (Gomes et al., 2015). Furthermore, this study sets the stage for a deeper inquiry into the dynamic biology of circulating exomiRs in response to exercise in the aging population.
Data Availability Statement
The data that support the findings of this study are available from the corresponding author upon request.
The studies involving human participants were reviewed and approved by Florida Hospital, AdventHealth Orlando (Orlando, FL) Institutional Review Board (IRBNet #554559). The patients/participants provided their written informed consent to participate in this study.
PC, VN, and SS designed the study. PC and BG recruited participants, collected samples, and performed experiments. VN, YG, SL, NJ, NS, MA, and MW performed experiments, collected and analyzed data, and revised the final version of the manuscript. VN and HP analyzed and interpreted data and drafted the manuscript. VN, HP, YG, PC, and SS critically revised the manuscript. All authors read and approved the final version of the manuscript.
This work was supported by funding from the National Institute on Aging (K01AG044437, to PC) and the Molecular Transducers of Physical Activity Consortium (MoTrPAC) program, which is supported by the National Institutes of Health (NIH) Common Fund (U24DK112331, to SS and MW).
Conflict of Interest
BG and PC were employed by the company AdventHealth.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Antos, C. L., McKinsey, T. A., Frey, N., Kutschke, W., McAnally, J., Shelton, J. M., et al. (2002). Activated glycogen synthase-3 beta suppresses cardiac hypertrophy in vivo. Proc. Natl. Acad. Sci. U.S.A. 99, 907–912. doi: 10.1073/pnas.231619298
Aoi, W., Ichikawa, H., Mune, K., Tanimura, Y., Mizushima, K., Naito, Y., et al. (2013). Muscle-enriched microRNA miR-486 decreases in circulation in response to exercise in young men. Front. Physiol. 4:80. doi: 10.3389/fphys.2013.00080
Baggish, A. L., Hale, A., Weiner, R. B., Lewis, G. D., Systrom, D., Wang, F., et al. (2011). Dynamic regulation of circulating microRNA during acute exhaustive exercise and sustained aerobic exercise training. J. Physiol. 589(Pt 16), 3983–3994. doi: 10.1113/jphysiol.2011.213363
Baggish, A. L., Park, J., Min, P. K., Isaacs, S., Parker, B. A., Thompson, P. D., et al. (2014). Rapid upregulation and clearance of distinct circulating microRNAs after prolonged aerobic exercise. J. Appl. Physiol. 116, 522–531. doi: 10.1152/japplphysiol.01141.2013
Baran-Gale, J., Kurtz, C. L., Erdos, M. R., Sison, C., Young, A., Fannin, E. E., et al. (2015). Addressing bias in small RNA library preparation for sequencing: a new protocol recovers microRNAs that evade capture by current methods. Front. Genet. 6:352. doi: 10.3389/fgene.2015.00352
Barber, J. L., Zellars, K. N., Barringhaus, K. G., Bouchard, C., Spinale, F. G., and Sarzynski, M. A. (2019). The effects of regular exercise on circulating cardiovascular-related microRNAs. Sci. Rep. 9:7527. doi: 10.1038/s41598-019-43978-x
Berntsen, S., Hageberg, R., Aandstad, A., Mowinckel, P., Anderssen, S. A., Carlsen, K. H., et al. (2010). Validity of physical activity monitors in adults participating in free-living activities. Br. J. Sports Med. 44, 657–664. doi: 10.1136/bjsm.2008.048868
Chaturvedi, P., Kalani, A., Medina, I., Familtseva, A., and Tyagi, S. C. (2015). Cardiosome mediated regulation of MMP9 in diabetic heart: role of mir29b and mir455 in exercise. J. Cell. Mol. Med. 19, 2153–2161. doi: 10.1111/jcmm.12589
Chen, S., Wu, J., Jiao, K., Wu, Q., Ma, J., Chen, D., et al. (2018). MicroRNA-495-3p inhibits multidrug resistance by modulating autophagy through GRP78/mTOR axis in gastric cancer. Cell Death Dis. 9:1070. doi: 10.1038/s41419-018-0950-x
Chowdhury, M. A., Sholl, H. K., Sharrett, M. S., Haller, S. T., Cooper, C. C., Gupta, R., et al. (2019). Exercise and cardioprotection: a natural defense against lethal myocardial ischemia-reperfusion injury and potential guide to cardiovascular prophylaxis. J. Cardiovasc. Pharmacol. Ther. 24, 18–30. doi: 10.1177/1074248418788575
Coen, P. M., Jubrias, S. A., Distefano, G., Amati, F., Mackey, D. C., Glynn, N. W., et al. (2013). Skeletal muscle mitochondrial energetics are associated with maximal aerobic capacity and walking speed in older adults. J. Gerontol. A Biol. Sci. Med. Sci. 68, 447–455. doi: 10.1093/gerona/gls196
Dard-Dascot, C., Naquin, D., d’Aubenton-Carafa, Y., Alix, K., Thermes, C., and van Dijk, E. (2018). Systematic comparison of small RNA library preparation protocols for next-generation sequencing. BMC Genomics 19:118. doi: 10.1186/s12864-018-4491-6
Denham, J., and Prestes, P. R. (2016). Muscle-enriched microRNAs isolated from whole blood are regulated by exercise and are potential biomarkers of cardiorespiratory fitness. Front. Genet. 7:196. doi: 10.3389/fgene.2016.00196
Dragovic, R. A., Gardiner, C., Brooks, A. S., Tannetta, D. S., Ferguson, D. J., Hole, P., et al. (2011). Sizing and phenotyping of cellular vesicles using Nanoparticle Tracking Analysis. Nanomedicine 7, 780–788. doi: 10.1016/j.nano.2011.04.003
Drummond, M. J., McCarthy, J. J., Sinha, M., Spratt, H. M., Volpi, E., Esser, K. A., et al. (2011). Aging and microRNA expression in human skeletal muscle: a microarray and bioinformatics analysis. Physiol. Genomics 43, 595–603. doi: 10.1152/physiolgenomics.00148.2010
D’Souza, R. F., Markworth, J. F., Aasen, K. M. M., Zeng, N., Cameron-Smith, D., and Mitchell, C. J. (2017). Acute resistance exercise modulates microRNA expression profiles: combined tissue and circulatory targeted analyses. PLoS One 12:e0181594. doi: 10.1371/journal.pone.0181594
D’Souza, R. F., Woodhead, J. S. T., Zeng, N., Blenkiron, C., Merry, T. L., Cameron-Smith, D., et al. (2018). Circulatory exosomal miRNA following intense exercise is unrelated to muscle and plasma miRNA abundances. Am. J. Physiol. Endocrinol. Metab. 315, E723–E733. doi: 10.1152/ajpendo.00138.2018
Dube, J. J., Coen, P. M., DiStefano, G., Chacon, A. C., Helbling, N. L., Desimone, M. E., et al. (2014). Effects of acute lipid overload on skeletal muscle insulin resistance, metabolic flexibility, and mitochondrial performance. Am. J. Physiol. Endocrinol. Metab. 307, E1117–E1124. doi: 10.1152/ajpendo.00257.2014
Faraldi, M., Gomarasca, M., Sansoni, V., Perego, S., Banfi, G., and Lombardi, G. (2019). Normalization strategies differently affect circulating miRNA profile associated with the training status. Sci. Rep. 9:1584. doi: 10.1038/s41598-019-38505-x
Fruhbeis, C., Helmig, S., Tug, S., Simon, P., and Kramer-Albers, E. M. (2015). Physical exercise induces rapid release of small extracellular vesicles into the circulation. J. Extracell. Vesicles 4:28239. doi: 10.3402/jev.v4.28239
Gao, L., Wang, X., Wang, X., Zhang, L., Qiang, C., Chang, S., et al. (2014). IGF-1R, a target of let-7b, mediates crosstalk between IRS-2/Akt and MAPK pathways to promote proliferation of oral squamous cell carcinoma. Oncotarget 5, 2562–2574. doi: 10.18632/oncotarget.1812
Garatachea, N., Pareja-Galeano, H., Sanchis-Gomar, F., Santos-Lozano, A., Fiuza-Luces, C., Moran, M., et al. (2015). Exercise attenuates the major hallmarks of aging. Rejuvenation Res. 18, 57–89. doi: 10.1089/rej.2014.1623
Gentleman, R. C., Carey, V. J., Bates, D. M., Bolstad, B., Dettling, M., Dudoit, S., et al. (2004). Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 5, R80. doi: 10.1186/gb-2004-5-10-r80
Gomes, C. P., Kim, T. K., Wang, K., and He, Y. (2015). The implications on clinical diagnostics of using microRNA-based biomarkers in exercise. Expert Rev. Mol. Diagn. 15, 761–772. doi: 10.1586/14737159.2015.1039517
Guescini, M., Canonico, B., Lucertini, F., Maggio, S., Annibalini, G., Barbieri, E., et al. (2015). Muscle releases alpha-sarcoglycan positive extracellular vesicles carrying miRNAs in the bloodstream. PLoS One 10:e0125094. doi: 10.1371/journal.pone.0125094
Han, T., Wu, N., Wang, Y., Shen, W., and Zou, J. (2019). miR1623p inhibits cell proliferation and migration and induces apoptosis by targeting PDPK1 in maxillary primordium mesenchymal cells. Int. J. Mol. Med. 43, 1441–1451. doi: 10.3892/ijmm.2019.4070
Hartjes, T. A., Mytnyk, S., Jenster, G. W., van Steijn, V., and van Royen, M. E. (2019). Extracellular vesicle quantification and characterization: common methods and emerging approaches. Bioengineering 6:7. doi: 10.3390/bioengineering6010007
Horak, M., Zlamal, F., Iliev, R., Kucera, J., Cacek, J., Svobodova, L., et al. (2018). Exercise-induced circulating microRNA changes in athletes in various training scenarios. PLoS One 13:e0191060. doi: 10.1371/journal.pone.0191060
Howden, E. J., Sarma, S., Lawley, J. S., Opondo, M., Cornwell, W., Stoller, D., et al. (2018). Reversing the cardiac effects of sedentary aging in middle age-A randomized controlled trial: implications for heart failure prevention. Circulation 137, 1549–1560. doi: 10.1161/CIRCULATIONAHA.117.030617
Jiang, D., Sun, X., Wang, S., and Man, H. (2019). Upregulation of miR-874-3p decreases cerebral ischemia/reperfusion injury by directly targeting BMF and BCL2L13. Biomed. Pharmacother. 117:108941. doi: 10.1016/j.biopha.2019.108941
Jiang, J., Xie, C., Liu, Y., Shi, Q., and Chen, Y. (2019). Up-regulation of miR-383-5p suppresses proliferation and enhances chemosensitivity in ovarian cancer cells by targeting TRIM27. Biomed. Pharmacother. 109, 595–601. doi: 10.1016/j.biopha.2018.10.148
Keller, P., Vollaard, N. B., Gustafsson, T., Gallagher, I. J., Sundberg, C. J., Rankinen, T., et al. (2011). A transcriptional map of the impact of endurance exercise training on skeletal muscle phenotype. J. Appl. Physiol. 110, 46–59. doi: 10.1152/japplphysiol.00634.2010
Kim, H. Y., Kwon, H. Y., Ha Thi, H. T., Lee, H. J., Kim, G. I., Hahm, K. B., et al. (2016). MicroRNA-132 and microRNA-223 control positive feedback circuit by regulating FOXO3a in inflammatory bowel disease. J. Gastroenterol. Hepatol. 31, 1727–1735. doi: 10.1111/jgh.13321
Kirschner, M. B., Kao, S. C., Edelman, J. J., Armstrong, N. J., Vallely, M. P., van Zandwijk, N., et al. (2011). Haemolysis during sample preparation alters microRNA content of plasma. PLoS One 6:e24145. doi: 10.1371/journal.pone.0024145
Knezevic, I., Patel, A., Sundaresan, N. R., Gupta, M. P., Solaro, R. J., Nagalingam, R. S., et al. (2012). A novel cardiomyocyte-enriched microRNA, miR-378, targets insulin-like growth factor 1 receptor: implications in postnatal cardiac remodeling and cell survival. J. Biol. Chem. 287, 12913–12926. doi: 10.1074/jbc.M111.331751
Lai, C. Y., Wu, Y. T., Yu, S. L., Yu, Y. H., Lee, S. Y., Liu, C. M., et al. (2014). Modulated expression of human peripheral blood microRNAs from infancy to adulthood and its role in aging. Aging Cell 13, 679–689. doi: 10.1111/acel.12225
Lai, Y. Y., Shen, F., Cai, W. S., Chen, J. W., Feng, J. H., Cao, J., et al. (2016). MiR-384 regulated IRS1 expression and suppressed cell proliferation of human hepatocellular carcinoma. Tumour Biol. 37, 14165–14171. doi: 10.1007/s13277-016-5233-5
Launay, T., Momken, I., Carreira, S., Mougenot, N., Zhou, X. L., De Koning, L., et al. (2017). Acceleration-based training: a new mode of training in senescent rats improving performance and left ventricular and muscle functions. Exp. Gerontol. 95, 71–76. doi: 10.1016/j.exger.2017.05.002
Li, M., Zeringer, E., Barta, T., Schageman, J., Cheng, A., and Vlassov, A. V. (2014). Analysis of the RNA content of the exosomes derived from blood serum and urine and its potential as biomarkers. Philos. Trans. R. Soc. Lond. B Biol. Sci. 369:20130502. doi: 10.1098/rstb.2013.0502
Liu, T. J., Wang, B., Li, Q. X., Dong, X. L., Han, X. L., and Zhang, S. B. (2018). Effects of microRNA-206 and its target gene IGF-1 on sevoflurane-induced activation of hippocampal astrocytes in aged rats through the PI3K/AKT/CREB signaling pathway. J. Cell. Physiol. 233, 4294–4306. doi: 10.1002/jcp.26248
Lovett, J. A. C., Durcan, P. J., and Myburgh, K. H. (2018). Investigation of circulating extracellular vesicle microRNA following two consecutive bouts of muscle-damaging exercise. Front. Physiol. 9:1149. doi: 10.3389/fphys.2018.01149
Ludwig, N., Leidinger, P., Becker, K., Backes, C., Fehlmann, T., Pallasch, C., et al. (2016). Distribution of miRNA expression across human tissues. Nucleic Acids Res. 44, 3865–3877. doi: 10.1093/nar/gkw116
Macias, M., Rebmann, V., Mateos, B., Varo, N., Perez-Gracia, J. L., Alegre, E., et al. (2019). Comparison of six commercial serum exosome isolation methods suitable for clinical laboratories. Effect in cytokine analysis. Clin. Chem. Lab. Med. 57, 1539–1545. doi: 10.1515/cclm-2018-1297
Mackey, D. C., Manini, T. M., Schoeller, D. A., Koster, A., Glynn, N. W., Goodpaster, B. H., et al. (2011). Validation of an armband to measure daily energy expenditure in older adults. J. Gerontol. A Biol. Sci. Med. Sci. 66, 1108–1113. doi: 10.1093/gerona/glr101
Meier, C., Hardtstock, P., Joost, S., Alla, V., and Putzer, B. M. (2016). p73 and IGF1R regulate emergence of aggressive cancer stem-like features via miR-885-5p COntrol. Cancer Res. 76, 197–205. doi: 10.1158/0008-5472.CAN-15-1228
Mikkelsen, U. R., Couppe, C., Karlsen, A., Grosset, J. F., Schjerling, P., Mackey, A. L., et al. (2013). Life-long endurance exercise in humans: circulating levels of inflammatory markers and leg muscle size. Mech. Ageing Dev. 134, 531–540. doi: 10.1016/j.mad.2013.11.004
Minna, E., Romeo, P., Dugo, M., De Cecco, L., Todoerti, K., Pilotti, S., et al. (2016). miR-451a is underexpressed and targets AKT/mTOR pathway in papillary thyroid carcinoma. Oncotarget 7, 12731–12747. doi: 10.18632/oncotarget.7262
Mitchell, C. J., D’Souza, R. F., Schierding, W., Zeng, N., Ramzan, F., O’Sullivan, J. M., et al. (2018). Identification of human skeletal muscle miRNA related to strength by high-throughput sequencing. Physiol. Genomics 50, 416–424. doi: 10.1152/physiolgenomics.00112.2017
Mooren, F. C., Viereck, J., Kruger, K., and Thum, T. (2014). Circulating microRNAs as potential biomarkers of aerobic exercise capacity. Am. J. Physiol. Heart Circ. Physiol. 306, H557–H563. doi: 10.1152/ajpheart.00711.2013
Munk, R., Martindale, J. L., Yang, X., Yang, J. H., Grammatikakis, I., Di Germanio, C., et al. (2019). Loss of miR-451a enhances SPARC production during myogenesis. PLoS One 14:e0214301. doi: 10.1371/journal.pone.0214301
Nielsen, S., Akerstrom, T., Rinnov, A., Yfanti, C., Scheele, C., Pedersen, B. K., et al. (2014). The miRNA plasma signature in response to acute aerobic exercise and endurance training. PLoS One 9:e87308. doi: 10.1371/journal.pone.0087308
Nielsen, S., Scheele, C., Yfanti, C., Akerstrom, T., Nielsen, A. R., Pedersen, B. K., et al. (2010). Muscle specific microRNAs are regulated by endurance exercise in human skeletal muscle. J. Physiol. 588(Pt 20), 4029–4037. doi: 10.1113/jphysiol.2010.189860
Ntoumou, E., Tzetis, M., Braoudaki, M., Lambrou, G., Poulou, M., Malizos, K., et al. (2017). Serum microRNA array analysis identifies miR-140-3p, miR-33b-3p and miR-671-3p as potential osteoarthritis biomarkers involved in metabolic processes. Clin. Epigenetics 9:127. doi: 10.1186/s13148-017-0428-1
Olioso, D., Dauriz, M., Bacchi, E., Negri, C., Santi, L., Bonora, E., et al. (2019). Effects of aerobic and resistance training on circulating micro-RNA expression profile in subjects with type 2 diabetes. J. Clin. Endocrinol. Metab. 104, 1119–1130. doi: 10.1210/jc.2018-01820
Olivieri, F., Ahtiainen, M., Lazzarini, R., Pollanen, E., Capri, M., Lorenzi, M., et al. (2014a). Hormone replacement therapy enhances IGF-1 signaling in skeletal muscle by diminishing miR-182 and miR-223 expressions: a study on postmenopausal monozygotic twin pairs. Aging Cell 13, 850–861. doi: 10.1111/acel.12245
Olivieri, F., Bonafe, M., Spazzafumo, L., Gobbi, M., Prattichizzo, F., Recchioni, R., et al. (2014b). Age- and glycemia-related miR-126-3p levels in plasma and endothelial cells. Aging 6, 771–787. doi: 10.18632/aging.100693
O’Neill, B. T., Lauritzen, H. P., Hirshman, M. F., Smyth, G., Goodyear, L. J., and Kahn, C. R. (2015). Differential role of insulin/IGF-1 receptor signaling in muscle growth and glucose homeostasis. Cell Rep. 11, 1220–1235. doi: 10.1016/j.celrep.2015.04.037
Peng, Y., Dai, Y., Hitchcock, C., Yang, X., Kassis, E. S., Liu, L., et al. (2013). Insulin growth factor signaling is regulated by microRNA-486, an underexpressed microRNA in lung cancer. Proc. Natl. Acad. Sci. U.S.A. 110, 15043–15048. doi: 10.1073/pnas.1307107110
Pichiorri, F., Suh, S. S., Rocci, A., De Luca, L., Taccioli, C., Santhanam, R., et al. (2010). Downregulation of p53-inducible microRNAs 192, 194, and 215 impairs the p53/MDM2 autoregulatory loop in multiple myeloma development. Cancer Cell 18, 367–381. doi: 10.1016/j.ccr.2010.09.005
Popov, D. V., Makhnovskii, P. A., Shagimardanova, E. I., Gazizova, G. R., Lysenko, E. A., Gusev, O. A., et al. (2019). Contractile activity-specific transcriptome response to acute endurance exercise and training in human skeletal muscle. Am. J. Physiol. Endocrinol. Metab. 316, E605–E614. doi: 10.1152/ajpendo.00449.2018
Radom-Aizik, S., Zaldivar, F. Jr., Oliver, S., Galassetti, P., and Cooper, D. M. (2010). Evidence for microRNA involvement in exercise-associated neutrophil gene expression changes. J. Appl. Physiol. 109, 252–261. doi: 10.1152/japplphysiol.01291.2009
Radom-Aizik, S., Zaldivar, F., Haddad, F., and Cooper, D. M. (2013). Impact of brief exercise on peripheral blood NK cell gene and microRNA expression in young adults. J. Appl. Physiol. 114, 628–636. doi: 10.1152/japplphysiol.01341.2012
Radom-Aizik, S., Zaldivar, F Jr., Leu, S. Y., Adams, G. R., Oliver, S., and Cooper, D. M. (2012). Effects of exercise on microRNA expression in young males peripheral blood mononuclear cells. Clin. Transl. Sci. 5, 32–38. doi: 10.1111/j.1752-8062.2011.00384.x
Rebelo-Marques, A., De Sousa Lages, A., Andrade, R., Ribeiro, C. F., Mota-Pinto, A., Carrilho, F., et al. (2018). Aging hallmarks: the benefits of physical exercise. Front. Endocrinol. 9:258. doi: 10.3389/fendo.2018.00258
Rozowsky, J., Kitchen, R. R., Park, J. J., Galeev, T. R., Diao, J., Warrell, J., et al. (2019). exceRpt: a comprehensive analytic platform for extracellular RNA profiling. Cell Syst. 8, 352–357.e3. doi: 10.1016/j.cels.2019.03.004
Russell, A. P., Lamon, S., Boon, H., Wada, S., Guller, I., Brown, E. L., et al. (2013). Regulation of miRNAs in human skeletal muscle following acute endurance exercise and short-term endurance training. J. Physiol. 591, 4637–4653. doi: 10.1113/jphysiol.2013.255695
Schmutz, S., Dapp, C., Wittwer, M., Vogt, M., Hoppeler, H., and Fluck, M. (2006). Endurance training modulates the muscular transcriptome response to acute exercise. Pflugers Arch. 451, 678–687. doi: 10.1007/s00424-005-1497-0
Selbach, M., Schwanhausser, B., Thierfelder, N., Fang, Z., Khanin, R., and Rajewsky, N. (2008). Widespread changes in protein synthesis induced by microRNAs. Nature 455, 58–63. doi: 10.1038/nature07228
Sellami, M., Gasmi, M., Denham, J., Hayes, L. D., Stratton, D., Padulo, J., et al. (2018). Effects of acute and chronic exercise on immunological parameters in the elderly aged: can physical activity counteract the effects of aging? Front. Immunol. 9:2187. doi: 10.3389/fimmu.2018.02187
Shah, R., Yeri, A., Das, A., Courtright-Lim, A., Ziegler, O., Gervino, E., et al. (2017). Small RNA-seq during acute maximal exercise reveal RNAs involved in vascular inflammation and cardiometabolic health: brief report. Am. J. Physiol. Heart Circ. Physiol. 313, H1162–H1167. doi: 10.1152/ajpheart.00500.2017
Small, E. M., O’Rourke, J. R., Moresi, V., Sutherland, L. B., McAnally, J., Gerard, R. D., et al. (2010). Regulation of PI3-kinase/Akt signaling by muscle-enriched microRNA-486. Proc. Natl. Acad. Sci. U.S.A. 107, 4218–4223. doi: 10.1073/pnas.1000300107
Stitt, T. N., Drujan, D., Clarke, B. A., Panaro, F., Timofeyva, Y., Kline, W. O., et al. (2004). The IGF-1/PI3K/Akt pathway prevents expression of muscle atrophy-induced ubiquitin ligases by inhibiting FOXO transcription factors. Mol. Cell 14, 395–403. doi: 10.1016/s1097-2765(04)00211-4
Tang, H., Lv, W., Sun, W., Bi, Q., and Hao, Y. (2019). miR505 inhibits cell growth and EMT by targeting MAP3K3 through the AKTNFkappaB pathway in NSCLC cells. Int. J. Mol. Med. 43, 1203–1216. doi: 10.3892/ijmm.2018.4041
Tang, Y. T., Huang, Y. Y., Zheng, L., Qin, S. H., Xu, X. P., An, T. X., et al. (2017). Comparison of isolation methods of exosomes and exosomal RNA from cell culture medium and serum. Int. J. Mol. Med. 40, 834–844. doi: 10.3892/ijmm.2017.3080
Thery, C., Witwer, K. W., Aikawa, E., Alcaraz, M. J., Anderson, J. D., Andriantsitohaina, R., et al. (2018). Minimal information for studies of extracellular vesicles 2018 (MISEV2018): a position statement of the International Society for Extracellular Vesicles and update of the MISEV2014 guidelines. J. Extracell. Vesicles 7:1535750. doi: 10.1080/20013078.2018.1535750
Troncoso, R., Ibarra, C., Vicencio, J. M., Jaimovich, E., and Lavandero, S. (2014). New insights into IGF-1 signaling in the heart. Trends Endocrinol. Metab. 25, 128–137. doi: 10.1016/j.tem.2013.12.002
Tung, Y. T., Hsu, Y. J., Liao, C. C., Ho, S. T., Huang, C. C., and Huang, W. C. (2019). Physiological and biochemical effects of intrinsically high and low exercise capacities through multiomics approaches. Front. Physiol. 10:1201. doi: 10.3389/fphys.2019.01201
Uesugi, A., Kozaki, K., Tsuruta, T., Furuta, M., Morita, K., Imoto, I., et al. (2011). The tumor suppressive microRNA miR-218 targets the mTOR component Rictor and inhibits AKT phosphorylation in oral cancer. Cancer Res. 71, 5765–5778. doi: 10.1158/0008-5472.CAN-11-0368
Vassilakos, G., Lei, H., Yang, Y., Puglise, J., Matheny, M., Durzynska, J., et al. (2019). Deletion of muscle IGF-I transiently impairs growth and progressively disrupts glucose homeostasis in male mice. FASEB J. 33, 181–194. doi: 10.1096/fj.201800459R
Wang, Z., Zhao, Z., Yang, Y., Luo, M., Zhang, M., Wang, X., et al. (2018). MiR-99b-5p and miR-203a-3p function as tumor suppressors by targeting IGF-1R in gastric cancer. Sci. Rep. 8:10119. doi: 10.1038/s41598-018-27583-y
Whitham, M., Parker, B. L., Friedrichsen, M., Hingst, J. R., Hjorth, M., Hughes, W. E., et al. (2018). Extracellular vesicles provide a means for tissue crosstalk during exercise. Cell Metab. 27, 237–251.e4. doi: 10.1016/j.cmet.2017.12.001
Xu, H., Liu, C., Zhang, Y., Guo, X., Liu, Z., Luo, Z., et al. (2014). Let-7b-5p regulates proliferation and apoptosis in multiple myeloma by targeting IGF1R. Acta Biochim. Biophys. Sin. 46, 965–972. doi: 10.1093/abbs/gmu089
Xu, Q., Jiang, Y., Yin, Y., Li, Q., He, J., Jing, Y., et al. (2013). A regulatory circuit of miR-148a/152 and DNMT1 in modulating cell transformation and tumor angiogenesis through IGF-IR and IRS1. J. Mol. Cell Biol. 5, 3–13. doi: 10.1093/jmcb/mjs049
Ye, X. W., Yu, H., Jin, Y. K., Jing, X. T., Xu, M., Wan, Z. F., et al. (2015). miR-138 inhibits proliferation by targeting 3-phosphoinositide-dependent protein kinase-1 in non-small cell lung cancer cells. Clin. Respir. J. 9, 27–33. doi: 10.1111/crj.12100
Yeri, A., Courtright, A., Reiman, R., Carlson, E., Beecroft, T., Janss, A., et al. (2017). Total extracellular small RNA Profiles from Plasma, Saliva, and urine of healthy subjects. Sci. Rep. 7:44061. doi: 10.1038/srep44061
Yuan, L., Zhou, C., Lu, Y., Hong, M., Zhang, Z., Zhang, Z., et al. (2015). IFN-gamma-mediated IRF1/miR-29b feedback loop suppresses colorectal cancer cell growth and metastasis by repressing IGF1. Cancer Lett. 359, 136–147. doi: 10.1016/j.canlet.2015.01.003
Zeng, Q., Wang, Y., Gao, J., Yan, Z., Li, Z., Zou, X., et al. (2019). miR-29b-3p regulated osteoblast differentiation via regulating IGF-1 secretion of mechanically stimulated osteocytes. Cell. Mol. Biol. Lett. 24:11. doi: 10.1186/s11658-019-0136-2
Zhang, C., Liu, J., Wang, X., Wu, R., Lin, M., Laddha, S. V., et al. (2014). MicroRNA-339-5p inhibits colorectal tumorigenesis through regulation of the MDM2/p53 signaling. Oncotarget 5, 9106–9117. doi: 10.18632/oncotarget.2379
Zhang, H., Yang, H., Zhang, C., Jing, Y., Wang, C., Liu, C., et al. (2015). Investigation of microRNA expression in human serum during the aging process. J. Gerontol. A Biol. Sci. Med. Sci. 70, 102–109. doi: 10.1093/gerona/glu145
Zhao, S., Gao, X., Zang, S., Li, Y., Feng, X., and Yuan, X. (2017). MicroRNA-383-5p acts as a prognostic marker and inhibitor of cell proliferation in lung adenocarcinoma by cancerous inhibitor of protein phosphatase 2A. Oncol. Lett. 14, 3573–3579. doi: 10.3892/ol.2017.6603
Keywords: Acute aerobic exercise, regular endurance training, aging, plasma exosomes, microRNA profiling, insulin growth factor-1 signaling
Citation: Nair VD, Ge Y, Li S, Pincas H, Jain N, Seenarine N, Amper MAS, Goodpaster BH, Walsh MJ, Coen PM and Sealfon SC (2020) Sedentary and Trained Older Men Have Distinct Circulating Exosomal microRNA Profiles at Baseline and in Response to Acute Exercise. Front. Physiol. 11:605. doi: 10.3389/fphys.2020.00605
Received: 13 March 2020; Accepted: 14 May 2020;
Published: 10 June 2020.
Edited by:Christopher S. Fry, University of Kentucky, United States
Reviewed by:Randall Felix D’Souza, The University of Auckland, New Zealand
Michael D. Roberts, Auburn University, United States
Copyright © 2020 Nair, Ge, Li, Pincas, Jain, Seenarine, Amper, Goodpaster, Walsh, Coen and Sealfon. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Venugopalan D. Nair, email@example.com