Participants’ characteristics at baseline
Out of 49 persons, 20 T2D and 23 individuals without diabetes (NDM), completed the study (Fig. S1). Defined by a cutoff-value of 4.9 mg*kg-1*min-1 for the M-value18, which reflects insulin-stimulated skeletal muscle glucose uptake16, 11 healthy persons were insulin-resistant (IR NDM) and 12 insulin-sensitive (IS NDM). Cardiorespiratory fitness (VO2max) was comparable between the three groups (Fig. 1a). M-value was lower in T2D and IR NDM than IS NDM (p<0.001, Fig. 1b). Insulin-mediated suppression of endogenous glucose production (iEGP), reflecting hepatic insulin sensitivity was lower (p<0.001, Fig. 1c), whereas liver fat content was higher in T2D than in both NDM groups (p<0.001, Fig. 1d). As expected, T2D also had higher HbA1c and lower HDL-cholesterol than both NDM groups. Compared to IS NDM, T2D had higher body mass index (BMI), partly due to due to higher visceral fat mass (Table 1), whereas the higher BMI of IR NDM resulted from both higher subcutaneous and visceral fat mass.
HIIT uniformly improves whole body oxidative and skeletal muscle mitochondrial capacity
After 12-week of HIIT, all participants uniformly increased their VO2max (Fig. 1a) despite unchanged body weight and whole-body fat content (Table 1). Maximal uncoupled respiration, as assessed ex vivo from permeabilized skeletal muscle fibers tended to be lower in T2D than in IS NDM at baseline (p=0.06), but rose by 28-45% in all groups after HIIT (Fig. 1e). HIIT neither affected leak control ratio (LCR), reflecting uncoupling at constant electron transport capacity, nor mitochondrial efficiency, assessed from respiratory control ratio (RCR), or reactive oxygen species (ROS) emission in any group (Fig. 1f,g). Muscle citrate synthase activity (CSA), as a surrogate marker of mitochondrial content19, at least doubled after HIIT in all groups without any difference between groups (Fig. 1h).
Mainly insulin-resistant persons (T2D, IR NDM) exhibit increased peripheral insulin sensitivity upon HIIT
Despite consistent increase in whole-body and skeletal muscle oxidative capacity, the response of M-value to HIIT substantially varied, separating groups of responders and non-responders. Logistic regression analysis showed an approximately doubled probability of the responder status among T2D (16/20; 80%) and IR NDM (9/11; 82%) as compared to IS NDM (5/12; 42%) (β=1.94, p=0.02 and β=-2.2, p=0.03, respectively). Likewise, only the insulin-resistant groups, T2D and IR NDM, improved their mean peripheral (Fig. 1b) and hepatic insulin sensitivity (Fig. 1c) as well as liver fat content (Fig. 1d). Only in T2D, HIIT also decreased visceral fat content, serum free fatty acids (FFA), triglycerides (TG) and showed a trend towards decreased glutamic-pyruvic transaminase (SGPT, p=0.06) (Table 1).
HIIT differently affects myocellular pathways of insulin sensitivity in the insulin-resistant responders (T2D-R, IR-R)
To further examine the heterogeneous responses of the peripheral insulin sensitivity (M-value) within the groups, we employed skeletal muscle biopsies from subgroups of insulin-resistant responders, (T2D-R, N=16; IR-R, N=9) and insulin-sensitive non-responders (IS-NR, N=7) (Fig. 2a). We primarily addressed key mechanisms known to underlie human myocellular insulin resistance, such as abnormal mitochondrial function and inhibitory lipid signaling via activation of the diacylglycerol/novel protein kinase C isoform (DAG/nPKC) pathway and/or inhibitory inflammatory pathways16,20.
Similar to the whole groups, HIIT increased mitochondrial function (Fig. 2b), but did not affect ROS in all subgroups, although ROS was lower in IR-R than in T2D-R (Fig. 2c).
Activation of the nPKCε and θ isoforms was assessed from translocation of the respective proteins to the myocellular membrane (Fig. S2a). PKCε activity negatively associated with M-value (log transformed; β=-0.55; p=0.002) at baseline, but not after HIIT. At baseline, PKCε activity was doubled in T2D-R compared to IR-R and IS-NR (Fig. 2d). After HIIT, T2D-R exhibited 50% and 40% decreases in PKCε and θ activities, respectively (Fig. 2d,e). Additional measurement of the nPKC activities in the whole groups of T2D, IR NDM and IS NDM, showed an identical pattern (Fig. S2b,c), indicating that the improvement in insulin resistance by HIIT is mediated by the nPKC pathway, at least in overt T2D.
Inflammatory pathways, specifically nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), not only relate to insulin resistance and oxidative stress, but may be also downregulated by chronic exercising21,22. Thus, we also analyzed muscle NF-κB protein levels in all groups, which did not differ between responders and non-responders at baseline, but were reduced by 24% after HIIT exclusively in IR-R (Fig. 2f). This exercise effect was confirmed in the whole groups of T2D, IR NDM and IS NDM, although baseline NF-κB levels were 40% higher in T2D (Fig. S2d). However, circulating pro- and anti-inflammatory cytokines (interleukin-1β, -6, -15, -1ra, tumor necrosis factor (TNF)α) were not significantly different between groups and after HIIT (data not shown). Collectively, these data suggest that the metabolic response to HIIT is at least partly mediated by reduction of nPKC activity in responders with diabetes (T2D-R) and by lower activation of NF-κB-dependent pathways in responders without diabetes (IR-R).
HIIT induces an increase in circulating SEV in the insulin-resistant responders (T2D-R, IR-R)
To clarify the role of SEV in exercise response, we measured the effect of HIIT on SEV release in representative subgroups of T2D-R (N=8), IR-R (N=8) and IS-NR (N=6). These subgroups showed overall similar changes as the whole cohort, i. e. increased VO2max in all, but increased peripheral and hepatic insulin sensitivity only in T2D-R and IR-R (Fig. S3a,b,c). Mitochondrial function increased in responders and was nominally higher compared to baseline in IS-NR (Fig. S3d), whereas ROS production remained unchanged in all three subgroups. Also, only T2D-R showed reduced PKCε activation (Fig. S3e), while only IR-R had lower NF-κB levels (Fig. S3f).
We measured the size and the number of circulating SEV in serum of responders and non-responders, collected at baseline and 72 h after the last bout of the 12-week HIIT. The purity of the SEV preparations was validated by nanoparticle tracking analysis (NTA) (Fig. 3a). The mean diameter of the isolated SEV was about 100 nm with a peak at about 80 nm and comparable between groups before and after HIIT (Fig. 3b). Interestingly, the estimated SEV concentration (number of circulating SEV per protein) was higher in IS-NR compared to T2D-R and numerically higher than in IR-R at baseline (Fig. 3). After HIIT, the SEV concentration rose only in the insulin-resistant groups, T2D-R and IR-R (Fig. 3d), with an increase of 52% (expressed as log2 fold change, FC) for both groups (Fig. S4a).
As HIIT increased the number of circulating SEV only in the insulin-resistant responders, we hypothesized that SEV and their cargo could contribute to the different cellular metabolic adaptations to this mode of chronic exercise training. We therefore characterized the proteome of the SEV isolated from T2D-R (N=5), IR-R (N=5) and IS-NR (N=5) individuals at baseline and after HIIT and identified a total of 1707 proteins including the 24 exosomal-enriched proteins (Table S1). Furthermore, 1589 of the identified proteins (98%) overlapped with the human exosome-associated proteins, previously identified in the Vesiclepedia database (13550 unique protein entries), whereas 39 were newly discovered as SEV-carried proteins (Fig. 4a, Table S2). Among the total SEV-associated proteins, we then selected the candidates with low variability (q value <0.05) and high degree of regulation (absolute log FC between groups or baseline vs HIIT >0.5) and we found that the proteomic profile differed between groups at baseline (Supplementary information) and after HIIT.
HIIT affects the proteomic profile of SEV
Quantitative proteomic analysis revealed that HIIT regulates the expression of 262 SEV proteins (n=122 in T2D-R, n=130 in IR-R, n=89 in IS-NR), of which 102 were downregulated and 160 upregulated (Fig. 4b, Table S3). Among the regulated SEV proteins, we identified proteins typically associated with exosomes23, such as biogenesis markers (ALIX), signaling proteins (GTPase, Ras-related protein), proteins associated with membrane trafficking and fusion (Rab proteins, annexins), lipid rafts (flotilin), cytoskeleton components (moesin, tubulin) and cell adhesion molecules (integrins). Moreover, the cellular component (CC) enrichment analysis of regulated SEV proteins confirmed a significant enrichment of proteins associated with extracellular exosomes (Fig. 4c). Of note, 13 of the 262 SEV proteins, differentially expressed after HIIT, were shared between all groups (Fig. 4d, Table S3), including antithrombin III, kininogen I, histidine-rich glycoprotein and α1-antitrypsin, which relate to inflammatory and immune responses24. Moreover, 29 SEV proteins were up- or down-regulated after HIIT exclusively in responders, but not in IS-NR (Fig. 4d, Table S3). Among these proteins, fibrinogen α, b and g chains (FGA, FGB, FGG) were similarly upregulated in SEV isolated from T2D-R and IR-R after HIIT. These acute-phase proteins are not only associated with insulin resistance and acutely increased by insulin in T2D25, but have been also described as myokine candidates carried by EV and released from the exercising limb after recovery14, suggesting that exercise might activate insulin sensitizing pathways.
Proteomic profiling of SEV suggests a new mode of myokine release for the metabolic adaptation to exercising
Since 12-week HIIT triggered the release of SEV, we assumed that SEV might represent an alternative to the release of biologically active proteins by classical secretory pathway. Indeed, we found that only 34% (89/262) of the proteins differentially expressed after HIIT had a predicted secretory signal peptide (SP) and only 12% (32/262) were predicted to follow a non-classical secretion pathway, based on the bioinformatics tools SignalP and SecretomeP (Fig. 4e, Table S3). These findings suggest that a large number of proteins can enter the circulation within SEV and that SEV can thereby contribute to inter-organ communication affecting cellular metabolism. In order to examine whether HIIT-induced increase in circulating SEV - at least partially - originated from skeletal muscle, we performed another quantitative proteomic analysis of SEV released from primary human skeletal muscle cells (hSkMC) using electrical pulse stimulation (EPS) to simulate the in vivo exercise intervention under in vitro conditions (Fig. 4f, Table S4). Interestingly, 44 SEV proteins regulated after HIIT were also present in SEV collected from the media of EPS-trained hSkMC, suggesting that these SEV proteins could represent novel myokines. Some of the SEV proteins upregulated, both after HIIT in vivo and after EPS in vitro, are indeed involved in signal transduction (GO:0007165; Ras-related protein Rap-2b and 1b), in oxidation-reduction processes (GO:0055114; aldehyde dehydrogenase family 16 member A1) and in cellular response to oxidative stress (GO:0034599; protein/nucleic acid deglycase DJ1,peroxiredoxin-2).
HIIT enriches SEV with proteins related to insulin sensitivity in T2D-R and inflammation and oxidative metabolism in IR-R
The 262 SEV proteins differentially regulated by HIIT were subsequently subjected to gene ontology (GO) analysis of biological processes (BP) and molecular function (MF) and to Ingenuity Pathway Analysis (IPA), as described in the supplementary information. In order to assess whether SEV protein cargo is responsible for the different metabolic adaptations induced by exercise in responders and non-responders, we performed a functional analysis of the exercise-regulated SEV proteins in each group (T2D-R, IR-R and IS-NR) and found different enriched GO-terms and pathways (Table S5).
In T2D-R, functional analysis revealed a significant enrichment of proteins linked to glycolytic process and the inhibition of glycolysis (z-score -2), PLC (z-score -2), mitogen-activated protein kinase (ERK/MAPK) (z-score -2.236) and protein kinase A (PKA) signaling (z-score -2.249), since the proteins associated to these pathways were downregulated after HIIT (glycolysis: ENO1, PGK1, PKM, TPI1; ERK/MAPK: HSPB1, PRKACB, PRKAR1A, RAB1B, PRKCB; PLC: BTK, GNAQ, PRKCB, RAB1B; PKA: GNAQ, MYLK, PRKACB, PKACB, PRKAR1A, PRKCB, RAB1B) (Fig. 5a). In line with the pathways analysis, IPA identified MAPK1 and interleukin-15 (IL15) as upstream molecules inhibited after HIIT (z-score -2.000 for both).
In IR-R, we found an overrepresentation of the non-canonical NF-κB and TNF-mediated pathways with upregulation of the 20S core proteasome complex (PSMA1, PSMA3, PSMA5, PSMA6, PSMA7, PSMB1, PSMB5, PSMB8)26,27 as well as the process “cellular response to oxidative stress” with an upregulation of proteins belonging to the cellular antioxidant system (CAT, CCS, G6PD, NME2, PRDX1, PRDX2, SOD2) (Fig. 5b). In line, IPA predicted the overrepresentation of pentose phosphate pathway, with an upregulation of glucose-6-phosphate dehydrogenase (G6PD), the rate-controlling enzyme of this pathway28, and the activation of the upstream molecule NFE2L2 (z-score 2.607) and the NRF2-mediated oxidative stress response (z-score 2) since the downstream targets of NRF2 (CAT, PRDX1, PRKCA, RAB1B, SOD2, AKR7A2) were upregulated after HIIT. IPA revealed also IL15 as activated upstream regulator (z-score 2.224), probably stimulated by the antioxidant system29. Finally, in the IR groups we found an enrichment of proteins associated to “response to calcium ion”, extending the concept of the involvement of Ca2+ in EV release30(Fig. 5a,b).
In IS-NR, we found a “lipid metabolism and transport” signature (Fig. 5c) and IPA revealed “autophagy of cells” as a cellular function associated only with SEV derived from IS-NR after 12-week HIIT (B-H p-value <0.05; z-score 1.432).
Taken together, these results indicate that HIIT (i) stimulates the insulin downstream pathway by inhibiting the MAPK, PLC and PKA signaling in the SEV from T2D-R and (ii) enhances the antioxidant system in the SEV released by IR-R, which eventually leads to improved peripheral insulin sensitivity, while (iii) it does not affect any pathway in SEV from IS-NR.
SEV cargoed proteins released after HIIT may directly affect cellular pathways in skeletal muscle
Finally, for investigating whether the specific SEV proteins also affect molecular pathways in skeletal muscle, we measured selected candidates related to insulin signaling for validation in muscle biopsies of T2D-R, IR-R and IS-NR individuals. Activating Thr172 phosphorylation of AMP activated protein kinase (AMPK) was higher in T2D-R than in IR-R at baseline and after HIIT (Fig. 6a). Inhibitory Ser307-phosphorylation of insulin receptor substrate 1 (IRS1) was numerically lower in IR-R after HIIT (p=0.13), whereas Ser1101-phosphorylation of IRS1 was not different between the groups (Fig. 6b,c). In addition, NRF2 and its downstream target NADPH quinone dehydrogenase 1 (NQO1) levels were lower at baseline in IR-R compared to T2D-R and tended to rise only in IR-R (Fig. 6d,e). After HIIT, also expression of both p38 mitogen-activated protein kinase (MAPK) and p44/p42 MAPK, as inflammation mediators target of TNFα31, were lower in IR-R and in IS-NR (Fig. 6f,g). Expression levels of microtubule-associated proteins 1A/1B light chain 3B (LC3) and ubiquitin-binding protein p62 (p62) was increased after HIIT only in IS-NR (Fig. 6h,i).