Diet modulates the development of insulin resistance during aging. This includes tissue-specific alterations in insulin signaling and mitochondrial function, which ultimately affect glucose homeostasis. Exercise stimulates glucose clearance and mitochondrial lipid oxidation and also enhances insulin sensitivity (IS). It is not well known how exercise interacts with age and diet in the development of insulin resistance. To investigate this, oral glucose tolerance tests with tracers were conducted in mice ranging from 4 to 21 months of age, fed a low-fat diet (LFD) or high-fat diet (HFD) with or without life-long voluntary access to a running wheel (RW). We developed a computational model to derive glucose fluxes, which were commensurate with independent values from steady-state tracer infusions. Values for an IS index derived for peripheral tissues (IS-P) and one for the liver (IS-L) were steeply decreased by aging and an HFD. This preceded the age-dependent decline in the mitochondrial capacity to oxidize lipids. In young animals fed an LFD, RW access enhanced the IS-P concomitantly with the muscle β-oxidation capacity. Surprisingly, RW access completely prevented the age-dependent IS-L decrease; however this only occurred in animals fed an LFD. Therefore, this study indicates that endurance exercise can improve the age-dependent decline in organ-specific IS if paired with a healthy diet.
Exercise is a known strategy to improve insulin sensitivity (IS), whereas aging and a lipid-rich diet decrease IS.
Using a tracer-based oral glucose tolerance test, we investigated how exercise, age, and diet interact in the development of tissue-specific insulin resistance.
Exercise (voluntary access to a running wheel) mainly improved IS in animals fed a low-fat diet. In these animals, exercise improved peripheral IS only at young age but fully prevented the age-dependent decline of hepatic IS.
The prevention of age-dependent decline in IS by exercise is tissue-specific and blunted by a lipid-rich diet.
Introduction
Aging, diet, and exercise have a profound effect on insulin sensitivity (IS) and glucose homeostasis (1–3). Consequently, they affect the emergence of insulin resistance, a pathophysiological process that precedes the development of type 2 diabetes. Insulin resistance is characterized by decreased insulin-stimulated glucose uptake and decreased insulin-dependent suppression of the endogenous glucose production (EGP) (2). In the postprandial state, together with the central nervous system, the skeletal muscle represents one of the major sites for glucose clearance, whereas the main site of EGP is the liver (4–6). Insulin resistance in the skeletal muscle and the liver can be induced by the accumulation of lipids, which involves the inhibition of the insulin signaling cascade by complex lipids (2). This lipid-induced insulin resistance can be elicited in mice fed a high-fat diet (HFD) (7). It is questionable, however, whether insulin resistance development is caused by an overload of substrate (8), by a loss of mitochondrial flexibility and/or function (9,10), or a combination of both.
Aging exacerbates the susceptibility to lipid-induced insulin resistance (11,12), which, among other factors, involves changes in lipid handling, both in the skeletal muscle (11,13) and the liver (14). The main route for lipid oxidation is the mitochondrial fatty-acid β-oxidation, with its flux being mostly controlled at the level of carnitine palmitoyltransferase 1 (CPT1) (13,15). We previously reported that an inability of old animals to upregulate muscle CPT1B (the muscle isoform) in response to an HFDfunctioned as a basis for insulin resistance development (13). Aging is also a risk factor for insulin resistance on its own (1), and its association with sarcopenia (loss of muscle mass) (16) may further exacerbate insulin resistance (17).
In humans, both resistance (strength) and endurance (aerobic) exercise have been associated with improved IS (18,19). Whereas resistance training rescues age-related muscle loss (20), endurance training is generally more effective in improving mitochondrial function in aging (16,21,22). The high mitochondrial capacity of endurance-trained athletes protects against lipid-induced insulin resistance (23,24). With aging, however, the link between mitochondrial capacity and IS as a result of exercise may disappear: It was found that aerobic exercise improved oxidative capacity irrespective of age, whereas it only improved IS in young people (24). Access to a running wheel (RW) is a widely used protocol to study the effects of endurance exercise in mice (25). In line with the above, in mice ranging from 6 to 24 months old, RW access increased mitochondrial content and oxidative capacity in the quadriceps, irrespective of whether a low-fat diet (LFD) or HFD was given. This was dissociated from the age-related decline in muscle weight (26). Finally, exercise affects IS and glucose handling in the muscle and the liver (27,28). After long-term endurance training, the ability of the liver to take up glucose increases and so does its sensitivity to insulin (27).
Although the effects of endurance exercise have been studied at different ages and on different diets, it has not been systematically addressed how these three factors interact in the development of muscle and liver insulin resistance. Moreover, a tissue-specific analysis of insulin resistance and glucose handling is crucial to dissect the role of mitochondrial function in insulin resistance. For this purpose, indices that distinguish muscle (peripheral) and liver IS are necessary. The frequently used HOMA-IR index is a measure of whole-body IS based on fasting glucose levels. An oral glucose tolerance test (OGTT) is used to derive the Muscle Insulin Sensitivity Index (MISI; focused on peripheral IS) or a Matsuda index (whole-body IS) during a glucose load (5,29,30). Including a stable-isotope-labeled glucose tracer to either oral or intravenous (IV) glucose tolerance tests enables one to dissect the contribution of the EGP from that of the glucose elimination by the muscle and other peripheral tissues and enables one to quantify IS in a tissue-specific manner (31–34). In the present study, we combine this advantage of a tracer with the OGTT.
Here, we investigate the combined effects of age, diet, and voluntary wheel running on glucose homeostasis and tissue specific IS in a cohort of aging mice about which we reported previously (26,35). To this end, we applied an existing OGTT protocol with [6,6-2H2]-glucose in mice (36) and adapted a computational model from earlier work in humans (37).
Research Design and Methods
Mouse Experiments and Ethics Approval
Mouse experiments were performed as previously described (26,35). Briefly, male C57BL/6JOlaHsd mice (Jackson Laboratory, Bar Harbor, ME) were maintained at 22°C on a 12-h light/12-h dark cycle and ad libitum access to an LFD (6% calories from fat, 1% w/w sugar; AMII 2141, HopeFarms BV, Woerden, the Netherlands). On postnatal day 28, mice were placed in individual cages (Makrolon type II; Bayer) and given either the same LFD or an HFD high in sucrose during their entire lifetime (45% calories from fat, 20% w/w sucrose; 4031.09, HopeFarms BV). A total of 128 mice were randomly subdivided into voluntary RW groups or sedentary controls (Ctrl). At the age of 4, 9, 15, or 21 months (used throughout the text to define the age groups), animals underwent an OGTT with a tracer (Fig. 1A). The same animals underwent steady-state IV infusion (SS-IV) experiments (38) at the age of 6, 12, 18, or 24 months, respectively, after which animals were terminated and the quadriceps and liver were collected, as previously described (26).
Schematic representation of experimental design (A) and compartment model for data analysis (B). GI, gastrointestinal; absorp., absorption flux into the plasma compartment, represented by the rate constant k1; loss, loss flux from the GI compartment, represented by the rate constant kL. Both tracer and unlabeled glucose in the plasma compartment can be cleared from the circulation, represented by the rate constant k2. EGP by the liver feeds into the unlabeled glucose pool in the plasma. Insulin enhances the clearance of both tracer and glucose and inhibits the EGP. The apparent absorption constant ka obtained from the tracer curves equals k1 + kL (See Research Design and Methods for details).
Schematic representation of experimental design (A) and compartment model for data analysis (B). GI, gastrointestinal; absorp., absorption flux into the plasma compartment, represented by the rate constant k1; loss, loss flux from the GI compartment, represented by the rate constant kL. Both tracer and unlabeled glucose in the plasma compartment can be cleared from the circulation, represented by the rate constant k2. EGP by the liver feeds into the unlabeled glucose pool in the plasma. Insulin enhances the clearance of both tracer and glucose and inhibits the EGP. The apparent absorption constant ka obtained from the tracer curves equals k1 + kL (See Research Design and Methods for details).
To determine the bioavailability of the tracer, independent experiments were conducted in 3-week-old male C57BL/6J mice (Jackson Laboratory), fed an LFD (Teklab Custom Diet TD.10098; Envigo Teklad Diets, Madison, WI) for 2 weeks prior to tracer measurements. Animal experiments were approved by the Dutch Central Authority for Scientific Procedures on Animals and by the University of Groningen Ethical Committee for Animal Experiments.
OGTT With Glucose Tracer
On the day of the experiment, food was removed in the morning at 7:00 (winter) or 8:00 (summertime). The OGTT was conducted after a fasting period of 6 h (i.e., at 13:00 or 14:00). Body weight was measured before the start of the experiment. Access to an RW was maintained during the OGTT time frame. In the RW groups, the last bout of intense exercise occurred 6–7 h before the ingestion of the glucose bolus, according to data on average RW revolutions at different ages (Supplementary Fig. 1). At time point zero, a glucose bolus of 1 g · kg−1 (5.5 · 103 µmol · kg−1) was administered orally. It consisted of 0.7 g · kg−1 (3.9 · 103 µmol · kg−1) unlabeled glucose and 0.3 g · kg−1 (1.6 · 103 µmol · kg−1) [6,6-2H2]-glucose (a tracer). Total blood glucose (Lifescan Euroflash; Lifescan Benelux, Beerse, Belgium) was measured from the tail vein, from which blood spots were also collected. Insulin was measured according to manufacturer instructions (Ultrasensitive Mouse Insulin kit; Mercodia, Uppsala, Sweden). The fractional contribution of the administered tracer in blood glucose (described below) was measured in air-dried blood spots. Glucose and tracer enrichment were measured at 0, 15, 30, 45, 60, 75, 90, 105, and 120 min and insulin was measured at 0, 30, and 60 min.
Bioavailability
Five-week-old animals from an independent cohort were fasted for 4 h and received an IV bolus of [6,6-2H2]-glucose (1.0 · 103 µmol · kg−1) injected retro-orbitally. This was directly followed by an oral bolus containing glucose (2.0 · 103 µmol · kg−1) and a tracer, [U-13C]-glucose (70 µmol · kg−1). Blood spots were collected at 0, 10, 20, 30, 40, 50, 60, and 90 min, and glucose was measured as described above. Air-dried blood spots were used for posterior analysis of fractional abundance.
Gas Chromatography Coupled to Mass Spectrometry and Isotope Correction
Fractional distributions of [6,6-2H2]-glucose and [U-13C]-glucose in blood spots were measured according to van Dijk et al. (38). Glucose was converted into its pentaacetate derivative by adding a pyridine–acetic anhydride mixture. All samples were analyzed by gas chromatography coupled to mass spectrometry (Agilent 9575C inert MSD; Agilent Technologies, Amstelveen, the Netherlands). Derivatives were separated on an AT-1701 30 m × 0.25 mm ID (0.25-µm film thickness) capillary column (Alltech, Breda, the Netherlands).
The monitored isotopologue distribution (m/z 408–414) was corrected for the natural abundance of isotopes based on the approach using the measured isotopologue distribution in the baseline samples, as described in the Supplementary Material. Tracer concentrations were obtained by multiplying the corrected fractional abundance of the m/z 410 isotopologue for [6,6-2H2]-glucose and m/z 414 for [U-13C]-glucose by the corresponding total glucose concentration. Unlabeled glucose concentrations were obtained by subtracting tracer concentrations from the total measured glucose concentrations.
Modeling Strategy
To derive the kinetic constants from the tracer data, we considered two compartments: the gastrointestinal tract compartment (compartment 1) and the blood plasma (compartment 2) (Fig. 1B). The pool sizes (q) and concentrations (c) of the tracer in these compartments are denoted by q1, q2, c1, and c2 and those of unlabeled glucose analogously by Q1, Q2, C1, and C2, with the subscript specifying the compartment. Pool sizes Q and q were expressed in µmol · kg−1 (referring to kg of body weight), whereas concentrations C and c were expressed in mmol/L. Inspired by the experimental setup, the model was confined to the pool size of the glucose that was administered in the stomach by oral gavage and dosed in µmol · kg−1 (Q1, q1) and by the plasma concentrations (C2, c2) that were measured in mmol/L. All modeling details can be found in the Supplementary Appendix.
CPT1 Quantification and Oxidative Capacity
Mitochondria were isolated from quadriceps and liver tissue, as previously described (26,39). CPT1A and CPT1B protein levels were quantified by targeted proteomics in isolated mitochondrial suspensions (40). Oxygen consumption rates were measured in MiR05 buffer (41) in a two-channel, high-resolution Oroboros oxygraph-2k (Innsbruck, Austria). Palmitoyl-CoA (25 μmol/L), l-carnitine (2 mmol/L) and malate (2 mmol/L) were used as substrates. Maximal ADP-stimulated respiration was measured after the addition of 1.5 U/mL hexokinase, 12.5 mmol/L glucose, and 1 mmol/L ATP. Protein content was measured with the Pierce BCA Protein Assay Kit (Thermo Fisher; catalog 23225). The ratio of cytosolic to mitochondrial protein was determined as previously described (26) and used to express CPT1 content and oxidative capacity per total tissue protein.
Software and Analysis
Correction for natural abundance was performed in Excel 2019. All other calculations and simulations were conducted in Python (version 3.7) using Jupyter Notebook (version 5.7.4). Parameters from the tracer kinetics and curves for the unlabeled glucose were estimated with the minimize function from the lmfit Python package (Levenberg-Marquardt method, nonlinear least squares). Data visualization and statistical analysis were conducted either in Python (matplotlib, scipy) or GraphPad Prism (version 9.0; GraphPad Software, San Diego, CA). Statistical analyses comparing the effects of age, diet, and exercise, as well as their interaction, were conducted by three-way ANOVA.
Data and Resource Availability
All data generated or analyzed during this study are included in the published article and its online supplementary files. The resources generated during this study can be found at github.com/mvieiralara/OGTT_modelling.
Results
HFD Reduces Glucose Elimination in an Age-Dependent Manner
Mice fed either an LFD or HFD were divided into a sedentary group (Ctrl) and a group with voluntary access to an RW. Physiological characteristics of these mice have been previously described (35). To study glucose kinetics and IS, a tracer OGTT was conducted in mice that were 4, 9, 15, or 21 months old (Fig. 1A). A computational model consisting of glucose absorption, loss, clearance, and EGP was built (Fig. 1B and Supplementary Material). The tracer time-course data were fitted to the equation in Fig. 1B (Fig. 2 and Supplementary Material).
Curve fits for time courses of the [6,6-2H2]-glucose tracer during OGTT. Each column represents a different diet and activity group, whereas each row represents a different age. Mean ± SD for each time point are shown in blue. The average curve fit from all animals per experimental group (c2 in the model) is shown in orange, with the SD of the fitted curves represented by the shaded area. n = 6–8 per group.
Curve fits for time courses of the [6,6-2H2]-glucose tracer during OGTT. Each column represents a different diet and activity group, whereas each row represents a different age. Mean ± SD for each time point are shown in blue. The average curve fit from all animals per experimental group (c2 in the model) is shown in orange, with the SD of the fitted curves represented by the shaded area. n = 6–8 per group.
The apparent rate constant of absorbance ka (Fig. 1B) was higher in the HFD than in the LFD groups (Pdiet < 0.01; Fig. 3A) and decreased with advanced age (Page < 0.001; Fig. 3A). The apparent rate constant k2 represents the fractional clearance rate of glucose from the plasma compartment (“glucose effectiveness”). The HFD decreased k2 (vs. LFD; Pdiet < 0.001), whereas k2 responded differently to aging in each diet group (Page × diet < 0.001), with no effect of RW access. In an independent tracer experiment in which IV and oral administration of the tracer were compared, the bioavailability F was determined to be 0.8 (Fig. 3C–F), implying that 80% of the tracer is absorbed and 20% is lost from our observation. Thus, k1 = 0.8 · ka and kL = 0.2 · ka (compare with Fig. 1B). An 80% bioavailability is in agreement with previous studies in humans and dogs (42,43). The apparent volume of distribution (see Vol, Eq. 14 in Supplementary Material) was lower in the HFD group than in the LFD group (Pdiet < 0.001) and responded in a diet-dependent manner to RW access (Pdiet × RW < 0.01) and age (Pdiet × age < 0.05) (Fig. 3G). In summary, all kinetic constants in the model could be obtained from the tracer OGTT if the bioavailability is independently determined. Moreover, the rate constant of glucose clearance was reduced by the HFD in an age-dependent manner.
The HFD reduced both the rate constant and the rate of glucose elimination. Kinetic constants ka (A) and k2 (B) calculated from curve fits. Equation for calculating the bioavailability F based on IV and oral experiments, in which the areas under the curve (AUC) are corrected for the dose and then divided (C). To calculate F, curves were fitted to tracer time courses after either IV (D) or oral administration (E) of glucose to 5-week-old mice. Data points are shown in blue as mean ± SD (n = 8), and the average curve fit from all animals are in orange; the SD of the fitted curve is represented by the shaded area. Calculated bioavailability F (F). Apparent volume of distribution (Vol) assuming a constant F for all groups (G). Basal rate of disappearance of glucose (Rd) calculated from the OGTT data (mean ± SEM) (H). Comparison between Rd calculated from the tracer OGTT in 4-month-old animals and SS-IV experiments conducted in the same animals at 6 months. Data are shown as mean ± SD; n = 6–8 (I).
The HFD reduced both the rate constant and the rate of glucose elimination. Kinetic constants ka (A) and k2 (B) calculated from curve fits. Equation for calculating the bioavailability F based on IV and oral experiments, in which the areas under the curve (AUC) are corrected for the dose and then divided (C). To calculate F, curves were fitted to tracer time courses after either IV (D) or oral administration (E) of glucose to 5-week-old mice. Data points are shown in blue as mean ± SD (n = 8), and the average curve fit from all animals are in orange; the SD of the fitted curve is represented by the shaded area. Calculated bioavailability F (F). Apparent volume of distribution (Vol) assuming a constant F for all groups (G). Basal rate of disappearance of glucose (Rd) calculated from the OGTT data (mean ± SEM) (H). Comparison between Rd calculated from the tracer OGTT in 4-month-old animals and SS-IV experiments conducted in the same animals at 6 months. Data are shown as mean ± SD; n = 6–8 (I).
The HFD treatment reduced the basal glucose disappearance rate (Rd) (Pdiet < 0.05), whereas Rd was affected by age in a diet-dependent manner (Page × diet < 0.001) (Fig. 3H). Subsequently, these outcomes were compared with those of classical SS-IV experiments (38), which allow for a model-independent measurement of basal glucose fluxes. In the 4-month-old group, basal values (Eq. 26 in Supplementary Material) were in the same range as those obtained from the tracer OGTT (Fig. 3I). Although this was true for all ages, the biological trends were not the same (Supplementary Fig. 2). This reflects the nature of the experiment: in the OGTT experiment, the k2 was estimated over the entire time course, during which the glucose bolus stimulated the release of hormones such as insulin. In the steady-state experiment, only tracer was infused, which should have barely affected hormone production.
Advanced Age Leads to Decreased EGP
The time courses of unlabeled glucose (Fig. 4) depend on the same kinetic parameters as that of the tracer andthe EGP (Fig. 1B). The variation within the groups was larger for unlabeled than for labeled glucose (compare with Figs. 2 and 4), indicating more variation in the EGP component, which may reflect stress due to animal handling (44). Individual curves that did not show the typical absorption followed by a clearance phase occurred in all experimental groups and were excluded (Supplementary Fig. 3). Fasting glucose values (Supplementary Fig. 4) were elevated in the HFD group (Pdiet < 0.01) and reduced by advanced aging (Page < 0.05), which, in turn, was modulated by the RW (PRW × age < 0.05).
Curve fits of unlabeled glucose time courses during OGTT. Each column represents a different diet and activity group, whereas each row represents a different age. Mean ± SD for each time point are shown in blue. The average curve fit (C2 in the model) is shown in orange, with the SD represented by the shaded area. n = 6–8 per group.
Curve fits of unlabeled glucose time courses during OGTT. Each column represents a different diet and activity group, whereas each row represents a different age. Mean ± SD for each time point are shown in blue. The average curve fit (C2 in the model) is shown in orange, with the SD represented by the shaded area. n = 6–8 per group.
After administration of the glucose bolus, the EGP decreased initially in most groups, reaching a minimum between 5 and 20 min, before reaching a new stationary level (Fig. 5A). The initial decline can be attributed to the inhibitory effect of increased plasma insulin and glucose (45). The steady-state EGP was calculated as the average of the last 30 min (Fig. 5B), and the overall average EGP was calculated over the entire period from 5 to 120 min (Fig. 5C). Higher age decreased both steady-state and average EGP (Page < 0.01) (Fig. 5B and C), which is in line with the usually lower glucose levels in older animals (Fig. 4). The RW generally decreased the steady-state and average EGP (vs. Ctrl, PRW < 0.05), mostly in the LFD group (PRW × diet < 0.05). The calculation of EGP in µmol · kg−1 · min−1 (Fig. 5) depends on the assumption of constant bioavailability. When EGP is expressed in mmol/L · min−1 (denoted EGP*), it was independent of the bioavailability (Supplementary Material) and followed the same trend (Supplementary Fig. 5A) as in Fig. 5A. Similar patterns were also observed for steady-state and average EGP* (Supplementary Fig. 5B and C). This highlights that the differences between groups are independent of the assumption of constant bioavailability.
EGP is reduced by aging in the LFD group. Time courses for EGP (A). Each column represents a different diet and activity group, whereas each row represents a different age. Mean EGP (line) ± SEM (shaded area) data are shown per experimental group. Steady-state EGP values (B) calculated from the curves (mean ± SEM). Time-averaged EGP values (C) obtained from OGTT time frame (5–120 min; mean ± SEM). Comparison between steady-state EGP values (mean ± SD) obtained from the tracer OGTT in 4-month-old animals and SS-IV experiments conducted in the same animals at 6 months (D). n = 2–8 per group.
EGP is reduced by aging in the LFD group. Time courses for EGP (A). Each column represents a different diet and activity group, whereas each row represents a different age. Mean EGP (line) ± SEM (shaded area) data are shown per experimental group. Steady-state EGP values (B) calculated from the curves (mean ± SEM). Time-averaged EGP values (C) obtained from OGTT time frame (5–120 min; mean ± SEM). Comparison between steady-state EGP values (mean ± SD) obtained from the tracer OGTT in 4-month-old animals and SS-IV experiments conducted in the same animals at 6 months (D). n = 2–8 per group.
The steady-state EGP derived from the tracer OGTT (EGPOGTT) was compared with independent SS-IV experiments conducted in the same cohort of animals. This yielded a model-independent basal EGP (EGPSS-IV). EGPOGTT and EGPSS-IV were similar to each other in 4-month-old animals (Fig. 5D). Although the same trend was observed in all age groups, the tracer OGTT usually yielded higher EGP values than the steady-state experiment (Fig. 5D and Supplementary Fig. 5D).
RW Access Mostly Affects Liver, Not Muscle, Insulin Sensitivity
Both age and an HFD strongly increased the insulin levels during the OGTT (for age and diet each, P < 0.001), in line with previous results (13,35), whereas RW access decreased insulin levels by ∼20% (PRW < 0.01) (Fig. 6A). HOMA-IR, a classic surrogate index for whole-body insulin resistance, followed a pattern similar to that of average insulin levels (Fig. 6B). Age and an HFD (for age and diet each, P < 0.001) increased HOMA-IR, whereas the RW access reduced it (PRW < 0.05). The effect of RW access was particularly pronounced in the LFD group.
Peripheral and central IS are regulated by age, diet, and RW access. Average insulin levels during the OGTT (0, 30, and 60 min) (A). HOMA-IR calculated from fasting glucose and fasting insulin levels (B). IS-P (C) and IS-L (D) for tracer OGTT, calculated according to the featured equations. INSOGTT is the time-average insulin levels per animal; EGP is the time-average EGP per animal (5–120 min). and are the average of EGP and INSOGTT values, respectively, for all 16 experimental groups. Data are shown as mean ± SEM; n = 2–8 per group.
Peripheral and central IS are regulated by age, diet, and RW access. Average insulin levels during the OGTT (0, 30, and 60 min) (A). HOMA-IR calculated from fasting glucose and fasting insulin levels (B). IS-P (C) and IS-L (D) for tracer OGTT, calculated according to the featured equations. INSOGTT is the time-average insulin levels per animal; EGP is the time-average EGP per animal (5–120 min). and are the average of EGP and INSOGTT values, respectively, for all 16 experimental groups. Data are shown as mean ± SEM; n = 2–8 per group.
Despite its wide application, the HOMA-IR does not capture organ-specific insulin resistance (5,46). We defined an IS index for peripheral tissues (IS-P) as the fractional clearance (k2; min−1) divided by the insulin concentration (mU · L−1) (Eq. 29 in Supplementary Material; Fig. 6). In the LFD groups, the IS-P decreased strongly between 4 and 9 months in sedentary and RW groups. In the HFD groups, the IS-P was already low when mice were young (4 months) and did not decline further during aging (Fig. 6C). Together, age, diet, and their interaction explained more than 50% of the variation in the data set (for age, diet, and age × diet each, P < 0.001). RW access significantly increased the IS-P (PRW < 0.05), particularly in the 4-month-old LFD group (∼50% increase vs. Ctrl) (Fig. 6C). RW access alone, however, explained <2% of the variation in the data. The interaction terms RW × diet and RW × age had P values of 0.06 and 0.07, respectively, which point to a trend of differential effects of RW access in older HFD groups.
Insulin also regulates the EGP, although this regulation is more complex to analyze (47). We adapted a hepatic IS index that was originally based on fasting EGP and insulin levels (29). Instead, we used the time-averaged EGP and insulin levels to apply the IS index for the liver (IS-L) to the postprandial state mimicked by the tracer OGTT (Eq. 30 in Supplementary Material and Fig. 6). Similar to the IS-P, the IS-L was strongly reduced by aging, and again, to a smaller extent, in the HFD group in which the IS-L was already low from 4 months of age (Fig. 6D). Surprisingly, in the LFD group, RW access strongly enhanced the IS-L, more than it affected the IS-P, thus delaying the age-related loss of the hepatic IS (Fig. 6D). Using the EGP* (Supplementary Fig. 5C) yielded the same pattern for the IS-L (Supplementary Fig. 6A), which reiterates that the effects of age, diet, and RW access do not depend on the assumption of constant tracer bioavailability.
Both IS-P and IS-L correlated with the HOMA-IR (Supplementary Fig. 6B and C). The differences between IS-P and IS-L curves warrant a strict dissection of whole-body IS by its peripheral and liver components.
RW Access Increases the β-Oxidation Capacity Primarily in the Muscle of LFD-Fed Animals
The quantitative dissection of peripheral and liver IS allowed us to relate these to mitochondrial function in specific tissues. Alterations in mitochondrial fatty-acid β-oxidation can modulate lipid clearance and impact insulin resistance development (48,49). The content of CPT1B in the skeletal muscle declined with age (Page < 0.001) in all analyzed groups (Fig. 7A), which aligns with results of an independent study (13). The HFD strongly upregulated muscle CPT1B content (Pdiet < 0.001), whereas RW access increased it only in LFD animals (PRW < 0.01 and PRW × diet < 0.001). In the skeletal muscle, the measured oxidative capacity with palmitoyl-CoA as substrate roughly correlated with the CPT1B levels (Fig. 7C, Supplementary Fig. 7A). However, LFD Ctrl mice did not show a decline in oxidative capacity with age, in contrast to the CPT1B level (compare with Fig. 7A and C). Of note, in the youngest and oldest groups, both CPT1B content and oxidative capacity followed the same pattern of age-dependent loss of flexibility to the HFD, as we previously noted (Supplementary Fig. 7C and D).
CPT1 content and oxidative capacity decline during aging. CPT1B protein content expressed per total skeletal (Sk.) muscle protein (A). CPT1A protein content expressed per total liver protein (B). Oxidative capacity (maximal ADP-stimulated O2 consumption) with the use of palmitoyl-CoA, carnitine, and malate as substrates for skeletal muscle (C) and liver (D) expressed per total tissue protein. CPT1B protein content was available from Stolle et al. (26). Skeletal muscle oxidative capacity data were reproduced from the same study (26). Data are shown as mean ± SEM; n = 4–8 per group.
CPT1 content and oxidative capacity decline during aging. CPT1B protein content expressed per total skeletal (Sk.) muscle protein (A). CPT1A protein content expressed per total liver protein (B). Oxidative capacity (maximal ADP-stimulated O2 consumption) with the use of palmitoyl-CoA, carnitine, and malate as substrates for skeletal muscle (C) and liver (D) expressed per total tissue protein. CPT1B protein content was available from Stolle et al. (26). Skeletal muscle oxidative capacity data were reproduced from the same study (26). Data are shown as mean ± SEM; n = 4–8 per group.
In the liver, CPT1A (liver isoform) content was also increased by the HFD (Pdiet < 0.05), and declined with age (Fig. 7B), even though there was no effect of RW access. A correlation between CPT1A and oxidative capacity was also observed in the liver, albeit weaker than in the muscle (Fig. 7D, Supplementary Fig. 7B). In the liver, the oxidative capacity was decreased by age and enhanced by the HFD (for age and diet each, P < 0.001), with no clear RW access effect (Fig. 7D).
The steep decline of the IS-P between 4 and 9 months in the LFD groups (Fig. 6C) preceded the decline of CPT1B content of the muscle. Nevertheless, RW access increased CPT1B content and oxidative capacity predominantly in the 4-month-old LFD group, the same group in which it rescued the IS-P the most. The correlation between CPT1B content and IS-P in the LFD group had an R2 of 0.55 (Supplementary Fig. 7E). which indicates that other factors contributed 45% to the modulation of peripheral IS by age and RW access. Altogether, the upregulation in β-oxidation capacity in the skeletal muscle by RW access only improved IS when there was no lipid overload (LFD groups).
Discussion
Here, we dissected in detail how age, diet, physical activity, and their interaction affect peripheral and hepatic IS and how these are associated with the mitochondrial fatty-acid oxidation capacity. To our knowledge, such a comprehensive analysis has never been conducted before. The most surprising result was that physical activity rescued the hepatic IS in aging mice, with only a minor age-dependent effect on the muscle. In addition, these beneficial effects of RW access were only observable in the LFD group. In agreement with earlier results in the muscle (13), the HFD-dependent induction of CPT1 and mitochondrial fatty-acid oxidation were lost in aged mice, both in the muscle and in the liver. This loss of mitochondrial flexibility was not directly reflected in IS, suggesting that the latter also depends on other factors. Finally, RW access, which results in a chronic exercise regimen, rescued this age-dependent loss of CPT1B and oxidative capacity only in the skeletal muscle of LFD animals, the same diet group in which some improvement of peripheral IS was observed. An overview of the main findings is given in Table 1.
Overview of main findings and strength of associations*
Parameter . | Age (old vs. young) . | Diet (HFD vs. LFD) . | Exercise (RW vs. Ctrl) . |
---|---|---|---|
k2 (clearance rate) | Diet-dependent effect | ↓ | — |
EGP | ↓↓ | — | ↓ (LFD group) |
IS-P | ↓↓ (LFD group) | ↓↓↓ | ↑ |
IS-L | ↓ | ↓↓ | ↑ (LFD group) |
Oxidative capacity quadriceps | ↓↓ | ↑↑↑ | ↑ |
Oxidative capacity liver | ↓↓↓↓ | ↑ | Age-dependent effect |
Parameter . | Age (old vs. young) . | Diet (HFD vs. LFD) . | Exercise (RW vs. Ctrl) . |
---|---|---|---|
k2 (clearance rate) | Diet-dependent effect | ↓ | — |
EGP | ↓↓ | — | ↓ (LFD group) |
IS-P | ↓↓ (LFD group) | ↓↓↓ | ↑ |
IS-L | ↓ | ↓↓ | ↑ (LFD group) |
Oxidative capacity quadriceps | ↓↓ | ↑↑↑ | ↑ |
Oxidative capacity liver | ↓↓↓↓ | ↑ | Age-dependent effect |
Arrows represent the direction of change (increase vs. decrease). The number of arrows represents to what extent a result is explained by each variable. 1 arrow: 0–10% variance, 2 arrows:10–20% variance, 3 arrows: 20–30% variance, 4 arrows: >30% variance, dash: no significant changes. Groups between parentheses indicate that the variance is further explained by interaction with the specified group.
The computational model to analyze the tracer OGTT was inspired by the Oral Minimal Model for humans (37). A difference between mouse and human studies is that less frequent and smaller blood volumes can be obtained from mice (38). Recently, a more complex computational model was fitted to mouse tracer-OGTT data that had been obtained with a similar experimental setup. The latter model included insulin-dependent and -independent glucose elimination (36). Although biologically more realistic, it is plausible that this higher complexity came at the expense of parameter identifiability (50), but this was not explicitly addressed (36). A limitation of the present study is that glucose bioavailability was considered constant, irrespective of age, diet, and physical exercise, even though this did not affect the relative differences between the groups. Moreover, the model does not accurately capture the initial dynamics of glucose absorption, given that transit compartments were not included. Despite these simplifications, the calculated steady-state fluxes were comparable to those obtained from independent steady-state infusion experiments. A final limitation concerns the use of only male mice because sex is known to play a role in the response to diet, age, and exercise (51,52).
Previously, peripheral IS has been defined as the ability of insulin to enhance glucose effectiveness (31). Because the latter is expressed by the apparent elimination rate constant k2, we defined the peripheral IS as IS-P = k2/INSOGTT. This parameter has the same units as the previously described insulin action parameter (31). IS-P is reminiscent of the MISI (5,30). The MISI is defined as the slope of the glucose decay during the clearance phase of the OGTT divided by the average insulin concentration during the same period. The slope depends on glucose absorption and elimination, which, in turn, depend on the glucose concentration. In contrast to the glucose decay slope used to calculate the MISI, k2 is specific for the elimination process and corrected for the glucose concentration. Similarly, the IS-L was adapted from an existing index for hepatic IS (29), the difference being that the IS-L analysis relies on the average response during the OGTT, whereas the classic analysis was based on fasting data. To capture the dynamics of hepatic and peripheral IS during the OGTT would require more data points and additional assumptions, or multiple tracer experiments (45).
HFDs have long been used to induce obesity and insulin resistance in animal models (53). In agreement with this, the HFD decreased both IS-P and IS-L (Fig. 6C and D). In previous studies, aging exacerbated HFD-induced insulin resistance (11,12). In this study, this was confirmed by the HOMA-IR, and it could be specifically attributed to the hepatic insulin response, as quantified by the IS-L. The peripheral insulin response (as indicated by the IS-P) was, however, equally low in all HFD groups, in contrast to our earlier findings that HFD-induced peripheral insulin resistance was aggravated by age (13). The studies differ, however, in terms of 1) diet duration, timing, and composition; and 2) mouse strain.
Despite the just-mentioned differences in experimental setup, the CPT1B content and the mitochondrial capacity to oxidize fatty acids in the skeletal muscle recapitulate our previous findings (Supplementary Fig. 7C and D). In both studies, older mice showed less flexibility to upregulate the fatty-acid oxidation capacity in the muscle when fed an HFD (13). A similar pattern was observed in the liver (Fig. 7D). This mechanism was proposed to underlie the high muscle-lipid accumulation in older HFDanimals, consequently aggravating insulin resistance. In the present study, however, the IS-P decline with age preceded the drop in CPT1B content and fatty-acid oxidation capacity (Fig. 7A and C). Consequently, additional insulin resistance–inducing mechanisms likely underlie the observed phenotype.
In this study, RW access was present throughout the life of the mice and was not removed during the OGTT time frame. However, because the mice did not exercise during the day when the OGTT was conducted (Supplementary Fig. 1), we address here the effects of chronic, and not acute, exercise on IS. The peripheral IS was barely enhanced by exercise in older groups, which corroborates the findings of Short et al. in humans (24). However, several human studies have pointed to resistance or high-intensity aerobic training as strategies to improve IS in the elderly (3,54,55). Moreover, training regimens that result in weight loss could be a more effective strategy (56). In the present cohort, weight was maintained (26). Different explanations might underlie the heterogeneity of outcomes. First, the type and duration of exercise differed among studies. Second, the quantification of insulin resistance ranged from simple whole-body HOMA-IR quantification to more complex tracer and clamp studies. Third, mice run less as they grow older (57,58), which was confirmed for the mice used here, especially for the HFD groups, and described in earlier reports about the same cohort (26,35). Nevertheless, voluntary running promoted beneficial health adaptations, including mitochondrial remodeling (25,26). We observed, however, that the insulin-desensitizing effects of the HFD outweighed the insulin-sensitizing response to endurance exercise. The net effect of diet and exercise on the insulin-signaling cascade dictates glucose uptake rates and, consequently, insulin resistance development or prevention.
Lifelong RW access decreased the EGP in the LFD group (Fig. 5B and C). This corroborates previous findings after a short-term RW intervention in mice (32). Moreover, the RW-access intervention remarkably sustained a prolonged, high hepatic IS during aging, albeit only in the LFD group. This was in stark contrast with the muscle, in which RW access did not prevent the age-dependent decline of the IS-P. The increase of the IS-L by RW access was not associated with CPT1A content or oxidative capacity. Rather, it may be attributed to the aforementioned increased capacity of hepatic glucose uptake, irrespective of insulin. Alternatively, it may derive from the action of insulin-sensitizing myokines, which can communicate the effect of exercise from the muscle to the liver (59).
In summary, based on a simple two-compartment model, EGP and muscle- and liver-specific IS could be reliably estimated from a noninvasive tracer-based OGTT. Hepatic and peripheral IS were both similarly decreased by aging and an HFD. In LFD-fed mice, exercise enhanced peripheral IS mostly in young animals and completely prevented the age-dependent loss of hepatic IS. These results, together with the applied methodology, may contribute to future research on the role of mitochondria in insulin resistance.
This article contains supplementary material online at https://doi.org/10.2337/figshare.22665826.
M.A.V.-L. and A.C.R. contributed equally to this work. G.v.D. and B.M.B. contributed equally to this work.
Article Information
Acknowledgments. We thank Maaike Oosterveer and Angela Tol of the Laboratory of Pediatrics of the University Medical Center Groningen for valuable discussions.
Funding. This study was funded by a grant from the Netherlands Organization for Scientific Research (NWO; grant 853.00.110), by DSM Nutritional products in the framework of the Complexity program (grant 645.001.001/3501), a University Medical Center Groningen (UMCG) MD/PhD fellowship to C.J.V., and a UMCG–Graduate School of Medical Sciences PhD fellowship to M.A.V.-L.
Duality of Interest. No potential conflicts of interest relevant to this article were reported.
Author Contributions. M.A.V.-L., S.K., and B.M.B. developed the computational model, advised by A.C.R., J.C., F.A., A.T., and T.H.v.D. conducted animal experiments. M.A.V.-L., A.C.R., and S.K. conducted data analysis. T.H.v.D. analyzed steady-state labeling data. J.C.W. conducted proteomics experiments. C.J.V. and R.H.J.B. contributed with the experiments on bioavailability. A.K.G., D.-J.R., G.v.D., and B.M.B. conceptualized and supervised the project. M.A.V.-L. and B.M.B. wrote the manuscript. All authors read and approved the final manuscript. B.M.B. is the guarantor of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Prior Presentation. Parts of this study were presented orally at the Bioinformatics & Systems Biology Conference 2021, virtual, 15–16 June 2021; and at the International Study Group for Systems Biology 2022, Innsbruck, Austria, 19–23 September 2022.