- 1The August Krogh Section for Human Physiology, Department of Nutrition, Exercise and Sports, University of Copenhagen, Copenhagen, Denmark
- 2Institute for Health and Sport, Victoria University, Melbourne, VIC, Australia
- 3College of Physical Education, Southwest University, Chongqing, China
- 4College of Sport Science, Zhuhai College of Science and Technology, Zhuhai, China
- 5Department of Functional Biology, University of Oviedo, Oviedo, Spain
- 6Department of Physiology and Pharmacology, Karolinska Institute, Stockholm, Sweden
Introduction: The effect of muscle glycogen stores on performance during intense short-duration exercises in humans is unclear. We hypothesized that low initial muscle glycogen levels would impair constant-load intense one-legged knee extensor exercise lasting approximately 5 min and human muscle contractile function, as determined by maximal voluntary contraction (MVC), electrically induced single-twitch maximal force, rate of force development (RFD), and rate of relaxation. Furthermore, alter phosphorylation of the Na+/K+-ATPase (NKA) regulatory proteins AMPK and FXYD1 indicating attenuated NKA activity.
Methods: With one leg in a glycogen-depleted state and one leg in a glycogen-loaded state, ten healthy young males (age: 25 ± 2 years) performed three intense exercise trials including (i) two-legged cycling for ∼5 min and (ii) 2× one-legged knee extensor exercise to task failure. MVC determination, electrical muscle stimulation, blood sample testing, and vastus lateralis biopsies were performed to assess the muscle composition and function.
Results: Time to task failure during the one-legged knee-extensor exercise was reduced by approximately 40% (n = 10, P < 0.05) with exercise- and diet-induced glycogen depletion. At rest (n = 10), MVC, twitch force, RFD, and rate of relaxation were unaffected by glycogen content. After exercise to task failure, the single-twitch contractile characteristics were impaired to a greater extent (n = 10, P < 0.05) in the glycogen-loaded leg than in the glycogen-depleted leg, probably induced by longer exercise duration. Concomitantly, MVC (n = 10, except for 15 s: n = 5 and 8) was reduced to similar levels under both conditions. The exercise-induced increase in the nonspecific phosphorylation of FXYD1 (n = 10, P < 0.001), which regulates NKA activity, tended to be greater (n = 10, P = 0.06) in the glycogen-loaded legs, indicating augmented potassium handling.
Conclusion: Exercise- and diet-induced low muscle glycogen content impairs high-intensity constant-load exercise performance lasting approximately 5 min. This was observed even with concomitantly better single-twitch contractile characteristics and similar reduction in MVC after task failure compared to glycogen-loaded legs. At rest, glycogen levels did not affect MVC and contractile characteristics.
1 Introduction
Task failure is related to reaching critically low muscle glycogen levels during exercise at a broad range of intensities (Bergstrom et al., 1967; Coyle et al., 1986; Greenhaff et al., 1987; Rockwell et al., 2003). At rest, approximately 400 mmol·kgd.w.−1 of glycogen is stored in the human vastus lateralis muscle, and the absolute content is highly dependent on dietary interventions, training status, and preceding exercise. In highly trained endurance athletes, exhaustion due to very low muscle glycogen content can be postponed if carbohydrates are ingested during prolonged strenuous exercise (Coyle et al., 1986).
However, the importance of the glycogen content in relation to short-duration intense exercise with combined aerobic and anaerobic energy contribution remains unclear. A low initial glycogen content (<200 mmol·kgd.w.−1) has been associated with impaired repeated 6 s maximal sprint performance (Balsom et al., 1999), whereas increasing the muscle glycogen content to approximately 700 mmol·kgd.w−1 has no effect on intense exercise performance lasting approximately 3 min (Bangsbo et al., 1992). Currently, the importance of an initially low muscle glycogen content for intense constant-load exercise capability is unclear (Vigh-Larsen et al., 2022), and no studies have included direct measurements of the muscle glycogen content. Therefore, the first aim of this study was to examine the impact of exercise- and diet-induced low muscle glycogen content on performance during brief (∼5 min) constant-load intense exercise with a pronounced aerobic and anaerobic energy production and a notable glycogen breakdown (Chidnok et al., 2013). We hypothesized that low initial muscle glycogen levels would impair performance determined as time to task failure in constant-load intense one-legged knee extensor exercise lasting approximately 5 min.
The mechanism by which low glycogen content compromises muscle function has not been fully elucidated (Vigh-Larsen et al., 2021). However, a reduced sarcoplasmic reticulum (SR) Ca2+ release rate (Ortenblad et al., 2011; Gejl et al., 2014) and impaired Ca2+ uptake caused by reduced Ca2+ ATPase (SERCA) activity (Duhamel et al., 2006b; Nielsen et al., 2009) and compromised membrane excitability (Nielsen et al., 2009; Renaud et al., 2023) are possible explanations. Reduced muscle glycogen content has been associated with a compromised SR Ca2+ release rate and half-relaxation time, as measured in vitro in muscle homogenate preparations (Nielsen et al., 2009; Ortenblad et al., 2011). In humans, contractile measurements are affected by fatigue (Booth et al., 1997), and accelerated reductions in Ca2+ release and uptake rates have been observed when exercising with low compared to higher muscle glycogen levels, as determined using in vitro assays (Duhamel et al., 2006a). If similar in vivo impairments in SR function occur during intense contractions in humans, the consequences of reduced glycogen content during constant-load exercise may include exacerbated rates of force development (RFD) and relaxation. However, few studies have investigated the effects of experimentally manipulated low glycogen content on human muscle contractile features. Thus, the second aim of this study was to examine whether human muscle contractile function, determined by the maximal voluntary contraction (MVC), electrically induced single-twitch maximal force, RFD, and rate of relaxation, is affected by the glycogen content. We hypothesized that exercise- and diet-induced low initial muscle glycogen content would impair human muscle contractile function, as determined by the MVC, single-twitch maximal force, RFD, and rate of relaxation before and immediately after the one-legged exercises to task failure.
Regarding compromised membrane excitability, increased Na+/K+-ATPase (NKA) activity during contractile activity is a prerequisite for delaying extracellular K+ build-up and preserving excitability (Hostrup et al., 2021). When glycogen is utilized during contractions, protein kinase A (PKA) activation causes phosphorylation of several proteins, including phospholemman (FXYD1) at Ser68, which inhibits NKA in the unphosphorylated state (Crambert et al., 2002). This results in increased NKA activity and maximal Na+/K+ transport capacity (Bibert et al., 2008). Exercise with low glycogen content upregulates the 5′AMP-activated protein kinase (AMPK) activity (Bartlett et al., 2013), which in contrast stimulates NKA activity (Ingwersen et al., 2011; Benziane et al., 2012). With low glycogen levels, a reduction in ATP provision from glycogenolysis and glycolysis can be expected (Dutka and Lamb, 2007; Jensen et al., 2020; Gaitanos et al., 1993). Also, ATP concentrations per see must be expected to be severely reduced (Karatzaferi et al., 2001). Reduced ATP availability compromise NKA activity. Therefore, the net effect of low glycogen content on the regulation of NKA activity remains unclear. Thus, the third aim of this study was to determine the extent to which the muscle glycogen content alters phosphorylation of the NKA regulatory proteins AMPK and FXYD1. We hypothesized that low muscle glycogen content alter exercise-induced phosphorylation of the NKA regulatory proteins, AMPK and FXYD1 indicating lower NKA activity and attenuated ion handling.
We investigated these hypotheses in healthy young male participants who performed two different exercise tasks, with one leg in an exercise- and diet-induced glycogen-depleted state and one leg in a glycogen-loaded state. The first task was to perform 5 min of intense two-legged cycling, where muscle biopsies were taken at rest and after exercise to investigate the NKA protein abundance and phosphorylation of NKA regulatory proteins. The second task involved two sessions of one-legged knee extensor exercise continued until task failure, in order to access possible differences in exercise performance and muscle contractile characteristics at rest, task failure, and recovery between the glycogen-loaded and glycogen-depleted legs.
2 Materials and methods
2.1 Participants
Ten healthy males provided written informed consent to participate in this study. Studies assessing muscle characteristics in humans, especially those incorporating muscle biopsies, frequently operate with small sample sizes due to ethical constraints, invasive methodology, and the high level of detail captured per participant, with similar studies also using comparable participant numbers (n = 8–10) (Duhamel et al., 2006a; Benziane et al., 2011; Bartlett et al., 2013). Participants were recreationally active, with a mean ± SD age of 25 ± 2 years, height of 184 ± 7 cm, and body mass of 79.0 ± 10.1 kg. The study was approved by the Ethics Committee of the Copenhagen and Frederiksberg communities (KF 01-183/02) and conducted in accordance with the Code of Ethics of the World Medical Association (Declaration of Helsinki) at the Department of Nutrition, Exercise and Sports, University of Copenhagen, Denmark.
2.2 Peak oxygen uptake and familiarization trials
The cycling tests were conducted using a Monark 839E ergometer (Monark Exercise AB, Vansbro, Sweden). All participants completed incremental cycle exercises comprising 10 min at 90 W and 10 min at 160 W, followed by 25 W increments each minute until exhaustion, to determine the peak oxygen uptake (
2.3 Glycogen depletion trials and dietary manipulations
The purpose of the first part of the experiment was to achieve one leg in a glycogen-loaded state and the other leg in a glycogen-depleted state prior to the experimental trials involving muscle biopsies and determination of exercise performance. The legs chosen to be glycogen-loaded or glycogen-depleted were randomized and counterbalanced.
For the glycogen-loaded leg, the participants first completed a glycogen depletion exercise, followed by consuming a high-carbohydrate diet, based on a previous protocol (Pilegaard et al., 2002). Both continuous and intermittent one-legged cycling exercise were performed. The one-legged

Figure 1. Schematic overview of the experimental setup. Red arrows indicate muscle biopsy sampling times, large black arrows indicate maximal voluntary contractions (MVC) and electrical stimulations (ES), and small black arrows indicate arm venous blood sampling times. LL: glycogen-loaded leg. DL: glycogen-depleted leg.
On the evening of day 3, participants completed a similar glycogen depletion protocol by exercising one leg (the glycogen-depleted leg) with minor workload adjustments. Participants started with a 6 min warm-up at 50 W, followed by 20 min of one-legged cycling at 60%
2.4 Exercise trials, muscle function determination, and blood sampling
2.4.1 Two-legged cycling
To determine both contractile characteristics with MVC and electrical stimulation as well as protein phosphorylation as close to task failure as possible, we performed two different exercise trials on the experimental day. First, participants completed a two-legged cycling session leading to task failure (i.e., participants were unable to sustain the initial power output) in 367 ± 111 s (Figure 1). The initial power output (288 ± 61 W) was set at approximately 90%–105%
2.4.2 One-legged knee extensor exercise
After 2 hours of rest with ad libitum water intake but no food ingestion, participants sat in a seat with one foot strapped in a custom boot attached to a one-legged knee extensor ergometer (Andersen et al., 1985). Participants performed a one-legged knee extensor exercise (63 ± 7 W) to task failure, when the RPM decreased from 60 to 45, with one leg. After 45 min of rest, the protocol was repeated with the other leg, in a randomized and counterbalanced order. The quadriceps MVC was determined before exercise as the highest of three 3–5 s maximal contractions, each separated by 30 s, and in recovery as single maximal contractions at 15, 30, 60, 120, 300, and 600 s post-exercise. For MVC measurements, the knee was fixed at an angle of 90°. The force was recorded using a calibrated strain gauge (Model 615; Tedea-Huntleigh Electronics, UK) attached to a custom boot, sampled at 1 kHz, and reported in Newtons (N). Immediately after the MVC was determined, at approximately 1 s, double-pulse electrical stimulation (400 V, approximately 100 mA) was delivered (2 × 1 ms separated by 10 ms, corresponding to 100 Hz) using a Digitimer DS7AH and DG2A generator (Digitimer, Hertfordshire, UK) and 5 × 9 cm2 electrodes (PALS platinum, Axelgaard, Lystrup, Denmark) placed over the belly of the rectus femoris. Prior to the experiment, the optimal stimulation intensity was determined by increasing the stimulation current until the twitch response plateaued, and 110% of this intensity was used in all experiments. Blood samples were collected before, during (60 s), immediately after, and 300 s during recovery from the high-intensity exercise bouts and immediately analyzed for plasma pH, [K+], [Na+], [lactate], and [glucose] (ABL 700, Radiometer, Copenhagen, Denmark).
2.5 Muscle biopsies
Under local anesthesia induced by xylocaine injections (2 mL of lidocaine without epinephrine, xylocaine 20 mg/mL; Astra Zeneca, Cambridge, UK), two incisions were made in the vastus lateralis of each leg. Four biopsies were taken from each participant: one in the glycogen-depleted leg and one in the glycogen-loaded leg before the two-legged cycling exercise (before exercise) and one in each leg immediately (first biopsy: 18 ± 8 s, second biopsy: 43 ± 9 s, depleted leg: 33 ± 17 s, loaded leg: 27 ± 14 s) after exercise (after exercise) (Figure 1). Biopsy samples were rapidly frozen in liquid nitrogen. Most of the rapidly frozen muscle tissue was freeze-dried and dissected free from blood, fat, and connective tissue before further analysis; approximately 20 mg was stored as whole muscle tissue for NKA content determination. The analyses of the muscle tissue were prioritized in the following order [protein abundance and phosphorylation (n = 10), glycogen content (n = 8–9), NKA content (n = 7–8) and muscle pH (n = 5–7)] and due to limited amount of tissue in some biopsies n differs between the obtained data.
2.5.1 Glycogen and pH content
Muscle glycogen content was determined in muscle samples using a modified spectrophotometric method (Harris et al., 1974). First, 1–2 mgd.w. of the freeze-dried muscle tissue samples were extracted in 0.5 mL 1 M HCl and hydrolyzed at 100°C for 3 h. Then, the glucose content was determined by the hexokinase method using a glucose kit and a PentraC 400 analyzer (TrioLab). The amount of NADPH formed was directly proportional to the glucose content and was measured spectrophotometrically at 340 nm. In addition, the muscle pH was measured using a small glass electrode (XC 161, Radiometer-analytical, France) after homogenization of approximately 1 mgd.w. of the samples in 100 µL non-buffered solution containing 145 mM KCl, 10 mM NaCl, and 5 mM sodium fluoride (Mannion et al., 1993).
2.5.2 Western blotting
The protein abundance and phosphorylation status were determined, as previously described by our laboratory (Thomassen et al., 2010). Briefly, samples of freeze-dried muscle tissue were homogenized (Qiagen Tissuelyser II, Retsch GmbH, Haan Germany) in fresh buffer containing 10% glycerol, 20 mM Na-pyrophosphate, 150 mM NaCl, 50 mM HEPES (pH 7.5), 1% NP-40, 20 mM β-glycerophosphate, 2 mM Na3VO4, 10 mM NaF, 2 mM PMSF, 1 mM EDTA (pH 8), 1 mM EGTA (pH 8), 10 μg mL–1 aprotinin, 10 μg mL–1 leupeptin, and 3 mM benzamidine. Next, samples underwent end-over-end rotation for 1 h at 4°C and were centrifuged at 16,500 g for 30 min at 4°C; the supernatant (lysate) was used for further analyses. The total protein concentration in each sample was determined using a BCA standard kit (Pierce), and samples were mixed with 6× Laemmli buffer (7 mL 0.5 M Tris-base, 3 mL glycerol, 0.93 g DTT, 1 g SDS, and 1.2 mg bromophenol blue).
Equal amounts of total protein were loaded into each well of a precast gel (Bio-Rad Laboratories, USA). Two protein markers on each gel (all blue and dual colors; Bio-Rad Laboratories, USA) were used to determine the molecular weights of the proteins quantified. Proteins were separated according to their molecular weights by SDS-PAGE and semi-dry transferred to a PVDF membrane (Millipore A/S, Denmark). The membranes were blocked in either 2% skimmed milk or 3% BSA in Tris-buffered saline including 0.1% Tween-20 (TBST) before an overnight incubation with the primary antibody at 4°C. The membranes were washed in TBST and incubated for 1 h at room temperature with horseradish peroxidase (HRP)-conjugated secondary antibody. Membranes were then washed three times for 15 min each in TBST. The bands were visualized with ECL (Millipore) and recorded using a digital camera (Kodak Image Station 2000MM, Kodak, Denmark). The band intensities were quantified using 1D software (Kodak, Denmark) and were determined as the total band intensity minus the background intensity. Individual values were normalized to the mean values of the resting glycogen-loaded leg samples.
2.5.3 Antibodies
The primary antibodies used in the present experiment were optimized using mixed human standard lysates, to ensure that the amount of protein loaded would result in band signal intensities localized to the steep, linear part of the standard curve. To determine changes in abundance of proteins involved in Na+ and K+ muscle cell homeostasis, the following antibodies and corresponding migration of the quantified signal were used: NKAα1: 100 kDa, C464.6 (Millipore); NKAα2: 100 kDa, 07-647 (Upstate); NKAβ1: 40–50 kDa, MA3-930 (Affinity Bioreagents); N-terminal FXYD1: 12 kDa, FXYD1 (kindly donated by M. Shattock, (Fuller et al., 2004); and Na+/H+ exchanger (NHE1): 110 kDa, MAB3140 (Chemicon).
Alterations in the phosphorylation status were determined using the following antibodies: acetyl-CoA carboxylase (ACC) βser221 phosphorylation: 260 kDa, phospho-specific ACCαser79 antibody (#07-303, Millipore); α-AMPKthr172 phosphorylation: 62 kDa, #2531 (Cell Signaling Technology); unspecific FXYD1 phosphorylation: 12 kDa, AB_FXYD1 (kindly donated by Dr. J. Randall Moorman, University of Virginia); FXYD1ser68 phosphorylation: 12 kDa, AB_FXYD1ser68 (kindly donated by Dr. D. Bers, Loyola University); FXYD1ser63 phosphorylation: 12 kDa, #2261 (Cell Signaling Technology); and combined FXYD1ser68 and thr69 phosphorylation, 12 kDa, #9621 (Cell Signaling Technology). HRP-conjugated goat anti-mouse and goat anti-rabbit (P-0447 and P-0448, respectively; DAKO, Denmark) secondary antibodies were used.
2.5.4 Antibody specificity
The phospho-specific ACCαser79 antibody (#07-303, Millipore) recognizes the equivalent ser221 in human ACCβ (Wojtaszewski et al., 2003), and ACCβser221 phosphorylation is a downstream target of AMPK and well-known indicator of AMPK activity (Birk and Wojtaszewski, 2006).
FXYD1 has three known phosphorylation sites (ser63, ser68, and thr69) located at its C-terminus (Bibert et al., 2008; Fuller et al., 2009). The polyclonal AB_FXYD1 antibody recognizes total unphosphorylated FXYD1 proteins with a 16 amino acid epitope at the C-terminus. Thus, the more phosphorylated FXYD1 proteins and the more phosphorylated FXYD1 sites, the lower the binding of AB_FXYD1. Therefore, the level of nonspecific FXYD1 phosphorylation was determined as AB_FXYD1/N-terminus total FXYD1, where a decrease in the AB_FXYD1 signal intensity represented an increase in nonspecific phosphorylation (Thomassen et al., 2011). FXYD1 phosphorylation on ser68 was first determined as (AB_FXYD1ser68/FXYD1) and then normalized to the nonspecific FXYD1 phosphorylated content determined by AB_FXYD1 (AB_FXYD1ser68/AB_FXYD1), to consider the concomitant phosphorylation of ser63 and thr69 (Thomassen et al., 2016). Furthermore, FXYD1ser63 phosphorylation and combined FXYD1ser68 and FXYD1thr69 phosphorylation were determined using motif- and phospho-specific antibodies #2261 and #9621, respectively, which have shown specificity for FXYD1 at approximately 12 kDa (Thomassen et al., 2011). Representative western blots are shown in Figure 2.

Figure 2. Representative western blots shown for one subject in duplicates. Muscle biopsy tissue was divided into two tubes after freeze drying and dissection, thereafter, treated separately with the same protocol to reduce method variation. The band migration is indicated at the right side. R: Rest and TF: Task failure. NKA: Na+/K+-ATPase. FXYD1: phospholemman. AMPK: 5′AMP-activated protein kinase. ACC: acetyl-CoA carboxylase.
2.5.5 [3H]-ouabain binding site content
To determine the total NKA content, the [3H]-ouabain binding site content was measured in the muscles before and after exercise in both glycogen-loaded and glycogen-depleted legs. Approximately 20 mg of whole muscle was used to measure the [3H]-ouabain binding site content, as previously described (Norgaard et al., 1984; Petersen et al., 2012). Briefly, each sample was washed twice for 10 min at 37°C in vanadate buffer (250 mM sucrose, 10 mM Tris-HCl, 3 mM MgSO4, and 1 mM NaVO4; pH 7.3). Following washing, samples were incubated for 2 h at 37°C in vanadate buffer with the addition of [3H]-ouabain (2.0 μCi mL-1 and 10-6 M, PerkinElmer Inc., Boston, MA, USA). The muscle was then placed in ice-cold vanadate solution for 4 × 30 min to remove any unbound [3H]-ouabain. Muscle samples were then blotted on filter paper and weighed before soaking in 500 μL 5% trichloroacetic acid and 0.1 mM ouabain for approximately 20 h. Then, 2.5 mL of a scintillation cocktail (Ultima Gold, Packard, PerkinElmer Inc., Boston, MA, USA) was added before the liquid scintillation counting of [3H]-ouabain. The [3H]-ouabain binding site content was calculated using the sample wet weight and specific activity of the incubation buffer and samples and was expressed as pmol·gw.w.−1. The final [3H]-ouabain content was then calculated, accounting for nonspecific binding, using a correction factor for [3H]-ouabain impurities, loss of bound [3H]-ouabain during washout, and incomplete saturation, as previously described (Norgaard et al., 1984).
2.6 Statistics
All statistical analyses were performed using SPSS (version 27; IBM Corp., Armonk, NY, US). The data were inspected for normality using Q-Q plots. Differences in one-legged high-intensity exercise performance were tested using a paired t-test. Changes in arm venous plasma variables before and after the two-legged cycling exercise were tested using one-way ANOVA with SIDAK-adjusted post hoc analysis. To estimate the differences in all other variables between treatment (glycogen-loaded vs. glycogen-depleted) and time (rest vs. task failure or rest, exercise, and recovery), a two-factor linear mixed model was used, with the glycogen content and time as fixed effects and a random effect for subjects. If a main effect or interaction was observed, a SIDAK-adjusted post hoc analysis was performed to determine specific differences. Data are presented as mean ± SD. P-values < 0.05 indicated statistical significance. Exact P-values were only shown if P < 0.10 and P > 0.05.
3 Results
3.1 Two-legged cycling
3.1.1 Muscle glycogen content
Muscle glycogen content decreased (main effect of exercise, P < 0.001) during the two-legged cycling exercise and was higher (main effect of treatment, P < 0.001) in the glycogen-loaded leg than in the glycogen-depleted leg, both at rest (593 ± 100 vs. 177 ± 29 mmol·kgd.w.−1, respectively; P < 0.001, n = 9) and at task failure (438 ± 166 vs. 77 ± 23 mmol·kgd.w.−1, respectively; P < 0.001, n = 8) in the two-legged cycling exercise (Figure 3A). The absolute reduction in glycogen content was not significantly different between legs (152 ± 100 vs. 98 ± 31 mmolkgd.w.−1; P = 0.23, n = 8).

Figure 3. (A) Total muscle glycogen content (mmol kgd.w.−1) in the glycogen-loaded and glycogen-depleted legs at rest and task failure from two-legged intense exercise cycling (n = 9, pre-exercise; n = 8, post-exercise). (B) Muscle pH values in the glycogen-loaded and glycogen-depleted legs at rest (n = 6) and task failure (n = 5 and n = 7, respectively). Data are reported as means, with individual values marked as white (rest) and black circles (task failure). There was a main effect of exercise (P < 0.001) for both the muscle glycogen content and pH, a main effect of treatment (P < 0.001) for the muscle glycogen content, and a main tendency of treatment (P = 0.075) and interaction (P = 0.065) for the muscle pH. ###(P < 0.001) and ##(P < 0.01) indicate a significant difference at rest and task failure from the intense exercise session in the same leg. *Indicates a significant difference between the glycogen-loaded and glycogen-depleted legs at task failure (P < 0.05).
3.1.2 Muscle homogenate pH
The muscle homogenate pH decreased during the two-legged cycling exercise (main effect of exercise, P < 0.001) in the glycogen-loaded leg, from 7.43 ± 0.10 (n = 6) to 6.81 ± 0.08 (n = 5, P < 0.001), and in the glycogen-depleted leg, from 7.42 ± 0.11 (n = 6) to 7.00 ± 0.14 (n = 7, P < 0.001). The muscle homogenate pH was lower (main tendency of treatment, P = 0.075; interaction, P = 0.065) in the glycogen-loaded than in the glycogen-depleted leg at task failure (P < 0.05) (Figure 3B).
3.1.3 Arm venous plasma variables
The two-legged cycling exercise induced the following respective changes in arm venous plasma values from rest (n = 7) to task failure and 300 s after exercise (n = 5): [lactate] increased from 0.77 ± 0.14 to 7.86 ± 3.43 mM (P < 0.001) and further to 13.24 ± 0.72 mM (P < 0.01); plasma [K+] increased from 3.63 ± 0.19 to 4.12 ± 0.44 mM (P < 0.05) and decreased after exercise to 3.28 ± 0.2 mM (P < 0.01); [Na+] increased from 138.9 ± 1.1 to 145.8 ± 3.1 mM (P < 0.001) and decreased after exercise to 142.4 ± 1.1 mM (P < 0.05); and pH declined from 7.38 ± 0.01 to 7.23 ± 0.07 and 7.19 ± 0.02 (P < 0.001). The data from the one-legged exercise are shown below.
3.1.4 Muscle protein abundances
NKAα1- and α2-subunit isoform protein abundances decreased after exercise compared to those before the two-legged cycling exercise (main effect of exercise, P < 0.05), and the NKAα1 abundance was lower in the glycogen loaded leg at task failure than at rest (P < 0.05). Both NKAα1 and NKAα2 abundances in the depleted leg tended (P = 0.091 and P = 0.060, respectively) to be lower at task failure than at rest (n = 10; Figures 4A,B). The NKAβ1-subunit isoform protein abundance was unchanged with exercise and was not different in the glycogen-loaded and glycogen-depleted legs (n = 10, Figure 4C). The FXYD1 abundance was reduced at task failure from exercise compared to that at rest (main effect of exercise, P < 0.05) and was lower in the depleted leg (P < 0.05) after exercise than at rest; however, there was no difference between the legs (n = 10, Figure 4D). The NHE1 abundance was unaffected by exercise, while glycogen depletion induced a higher NHE1 abundance (main effect of treatment, P < 0.01), which was not significant at rest or task failure (n = 10; data not shown).

Figure 4. Protein abundance of the NKA α1 (A), α2 (B), β1 (C) isoforms and FXYD1 (D) in the glycogen-loaded and glycogen-depleted legs at rest and task failure from two-legged intense exercise cycling. Data are normalized to the mean rest value in the glycogen-loaded leg and reported in arbitrary units as means (n = 10), with individual values marked as white (rest) and black circles (task failure). There was a main effect of exercise (P < 0.05) for NKAα1, NKAα2, and FXYD1. #Indicates a significant difference at rest versus task failure from the intense exercise session in the same leg (P < 0.05).
The total NKA content, determined by the [3H]-ouabain binding site content, was not different at rest and task failure and was not different between the glycogen-loaded and glycogen-depleted legs at rest (310 ± 66 vs. 337 ± 47 pmol·gw.w.−1; n = 8) or task failure (337 ± 64 vs. 338 ± 56 pmol·gw.w.−1; n = 7 and n = 8, respectively).
3.1.5 Protein phosphorylation
Nonspecific FXYD1 phosphorylation increased with exercise (main effect of exercise, P < 0.001) and was increased at task failure compared to that at rest in both the glycogen-loaded (P < 0.001) and glycogen-depleted legs (P < 0.01; Figure 5A). Furthermore, nonspecific FXYD1 phosphorylation tended (main tendency of treatment, P = 0.064) to be higher in the depleted legs than in the loaded legs. At rest, the level of nonspecific FXYD1 phosphorylation was greater (P < 0.05) in the glycogen-depleted legs than in the glycogen-loaded legs, while no difference was observed at task failure between the legs (n = 10; Figure 5A).

Figure 5. Nonspecific and specific FXYD1 phosphorylation levels in the glycogen-loaded and glycogen-depleted legs at rest and task failure from two-legged intense exercise cycling. (A) AB_FXYD1/FXYD1 determines the nonspecific FXYD1 phosphorylation normalized to the total FXYD1 abundance (see methods). A decrease in the AB_FXYD1 signal intensity represents an increase in nonspecific phosphorylation (Thomassen et al., 2011). (B) FXYD1ser68 phosphorylation normalized to total FXYD1 abundance. (C) Normalized FXYD1ser68 phosphorylation determines the FXYD1ser68 phosphorylation by considering concomitant phosphorylation on FXYD1ser63 and thr69 by normalizing to the AB_FXYD1 signal intensity (Thomassen et al., 2016). (D) FXYD1ser68+thr69 phosphorylation normalized to total FXYD1 abundance. Data are normalized to mean values at rest in the glycogen-loaded leg and reported in arbitrary units as means (n = 10), with individual values marked as white (rest) and black circles (task failure). There was a main effect of exercise for nonspecific FXYD1 phosphorylation (P < 0.001) and normalized FXYD1ser68 phosphorylation (P < 0.05) and a main effect of treatment (P < 0.05) for combined FXYD1ser68+thr69 phosphorylation and a main tendency of treatment (P = 0.064) for nonspecific FXYD1 phosphorylation. ###(P < 0.001) and ##(P < 0.01) indicate a significant difference between rest and task failure from exercise in the same leg. *Indicates a significant difference between the glycogen-depleted and glycogen-loaded legs at rest (P < 0.05).
The FXYD1 phosphorylation at ser68 was unaltered by exercise or glycogen intervention (n = 10; Figure 5B). When considering concomitant phosphorylation at ser63 and thr69 (see Section 2.5.4), FXYD1ser68 phosphorylation increased (main effect of exercise, P < 0.05) during exercise and tended to be higher (P = 0.055) in the loaded leg at task failure but was unaffected by the glycogen level (n = 10; Figure 5C).
The combined FXYD1ser68 and FXYD1thr69 phosphorylation was unaffected by exercise and was higher (main effect of treatment, P < 0.05) in the glycogen-loaded leg than in the depleted leg but not at rest or task failure (n = 10; Figure 5D). The FXYD1ser63 phosphorylation remained constant during exercise and was unaffected by glycogen levels (data not shown).
AMPKthr172 phosphorylation tended to be increased by exercise (main tendency of exercise, P = 0.071) and tended to be higher (P = 0.084) at task failure in the depleted leg than that at rest (Figure 6A). The AMPKthr172 phosphorylation was higher (main effect of treatment, P < 0.05) in the glycogen-depleted leg than in the loaded leg but not before or after exercise (n = 10; Figure 6A). The ACCβser221 phosphorylation increased with exercise (main effect of exercise, P < 0.001) and was higher at task failure in both the glycogen-loaded (P < 0.001) and glycogen-depleted (P < 0.001) legs. The glycogen level had no effect on the ACCβser221 phosphorylation (n = 10; Figure 6B).

Figure 6. AMPK thr172 (A) and ACCβ ser221 (B) phosphorylation levels in the glycogen-loaded and glycogen-depleted legs at rest and task failure from two-legged intense exercise cycling. Data are normalized to mean values at rest in the glycogen-loaded leg and reported in arbitrary units as means (n = 10), with individual values marked as white (rest) and black circles (task failure). There was a main effect of exercise (P < 0.001) for ACCβ ser221 phosphorylation and a main effect of treatment (P < 0.05) and main tendency of exercise (P = 0.071) for AMPK thr172 phosphorylation. ###Indicates a significant difference between rest and task failure from intense exercise in the same leg (P < 0.001).
3.2 One-legged knee extensor exercise to exhaustion
3.2.1 Performance
Time to task failure in the one-legged knee extensor exercise was 65% longer (P < 0.05) in glycogen-loaded legs than in glycogen-depleted legs (258 ± 181 vs. 156 ± 75 s; n = 10; Figure 7).

Figure 7. Time to task failure in the one-legged knee extensor exercise performed with the glycogen-loaded legs (blue) and the glycogen-depleted legs (green). Data are shown as means with individual values marked as black circles. *Indicates a significant difference between the glycogen-loaded and glycogen-depleted legs (P < 0.05).
3.2.2 MVC
The quadriceps MVC was similar between the glycogen-loaded and glycogen-depleted legs, both at rest and task failure (Table 1). Specifically, the MVC was reduced (main effect of time, P < 0.001) to approximately 62% (P < 0.01), 75% (P < 0.01), and 80% (P < 0.05) at 15, 30, and 60 s of recovery, respectively. The MVC tended (main tendency for treatment, P = 0.06) to be higher (P = 0.08) in the glycogen-loaded leg than in the glycogen-depleted leg after 15 s of recovery.

Table 1. Quadriceps MVC pre and post dynamic one-legged knee extensor exercise to task failure in glycogen-depleted and glycogen-loaded conditions.
3.2.3 Electrically induced isometric twitch before and after dynamic one-legged knee extensor exercise
The quadriceps electrically induced maximal twitch force was similar at rest in the glycogen-loaded and glycogen-depleted legs (114 ± 40 vs. 112 ± 32 N, respectively; n = 10) and was reduced (main effect of time, P < 0.001) at 30 s after exercise (P < 0.001) in both legs to 60%–65% of the resting single-twitch maximal force. A main effect of treatment (P < 0.05) was evident for the single-twitch maximal force, which was greater in the depleted leg than in the loaded leg, but with no detectable time-specific differences between the glycogen-loaded and glycogen-depleted legs (Table 2).

Table 2. Muscle contraction characteristics in recovery from dynamic one-legged knee extensor exercise to task failure in glycogen-depleted and glycogen-loaded conditions.
The RFD, determined from the electrically induced twitch, was similar in the glycogen-loaded and glycogen-depleted legs at rest (1,448 ± 542 vs. 1,417 ± 455 N/s, respectively; n = 10). The RFD was reduced (main effect of time, P < 0.001) at 30 s after exercise (P < 0.001) to 59%–64% of resting values (Table 2). The main effect of the treatment (P < 0.01) was evident as a higher RFD in the glycogen-depleted leg than in the glycogen-loaded leg, with a significant difference between legs detected 60 s post-exercise (P < 0.05; Table 2).
The rate of relaxation, determined from the electrically induced twitch, was similar in the glycogen-loaded and glycogen-depleted legs at rest (−787 ± 296 vs. −771 ± 204 N/s; n = 10). The rate of relaxation was markedly reduced (main effect of time, P < 0.001) at 30 s after exercise (P < 0.001) to 45%–55% of resting values (Table 2). A main effect of treatment (P < 0.01) and an interaction between treatment and time (P < 0.05) were evident as a lower rate of relaxation in the glycogen-loaded leg than in the glycogen-depleted leg, with a significant difference between the legs at 60 s post-exercise (P < 0.05; Table 2).
3.2.4 Arm venous plasma concentrations and pH
The arm venous plasma [lactate] increased (main effect of time, P < 0.001) during the one-legged knee extensor exercise and was higher at task failure and recovery than before exercise in both legs. The arm venous plasma [lactate] was higher (main effect of treatment, P < 0.05) at task failure and recovery (P < 0.05) in the glycogen-loaded leg than in the glycogen-depleted leg (Table 3). The plasma [K+] was higher (main effect of time, P < 0.001) at task failure than before exercise (60 s) and during recovery from the one-legged knee extensor exercise in the glycogen-loaded (P < 0.05) and glycogen-depleted legs (P < 0.001), with no differences between treatments (Table 3). The arm-venous plasma [Na+] was higher (main effect of time, P < 0.001) at task failure than at rest and after 60 s (P < 0.05) of the one-legged knee extensor exercise, with no differences between treatments (Table 3). The arm venous pH decreased (main effect of time, P < 0.001) at 300 s of recovery from the one-legged knee extensor exercise compared to that at rest in both the glycogen-loaded (P < 0.05) and glycogen-depleted legs (P < 0.01) (Table 3). The arm venous [glucose] was similar before, during, and after the one-legged knee extensor exercise in both the glycogen-loaded and glycogen-depleted legs (Table 3).

Table 3. Arm venous plasma electrolyte, glucose concentrations and pH before, during and after one-legged knee extensor exercise to task failure in glycogen-depleted and glycogen-loaded conditions.
4 Discussion
The main finding of the present study was that an exercise- and diet-induced low muscle glycogen content, compared to a high content, impaired performance during high-intensity one-legged knee extensor exercise, with an approximately 40% reduction in the time to task failure. Furthermore, resting muscle contractile characteristics, determined by the MVC as well as the electrically-induced single-twitch maximal force, RFD, and rate of relaxation, were unaffected by the glycogen content. Despite the higher muscle glycogen content in the glycogen-loaded condition, the single-twitch maximal force, RFD, and rate of relaxation were impaired to a greater extent at task failure than in the glycogen-depleted condition, whereas the MVC was similar. This suggests that a longer exercise duration in the glycogen-loaded leg was a decisive factor for impaired electrically induced contractile characteristics, while the MVC was unaffected by differences in the exercise duration. Manipulation of the muscle glycogen content did not affect the NKA content, whereas the exercise-induced increase in nonspecific phosphorylation of the NKA activity regulating FXYD1 protein after 5–6 min of cycling tended to be augmented in the glycogen-loaded leg.
4.1 Low glycogen content impairs continuous brief intense exercise performance
In this study, we showed, for the first time, that exercise-induced low muscle glycogen content impairs performance during continuous brief intense exercise lasting approximately 5 min (2.6 vs. 4.3 min). This supports the findings of Balsom et al. (1999) that the capacity to perform 15 repeated 6 s intense cycle intervals is impaired with low (180 mmol·kgd.w.−1) compared to normal (397 mmol·kgd.w.−1) glycogen levels (Balsom et al., 1999), albeit with markedly different exercise protocols. In contrast, elevated glycogen content compared to normal (approximately 700 vs. 350 mmol·kgd.w.−1) does not affect the time to exhaustion during intense exercise lasting approximately 3 min (Bangsbo et al., 1992), indicating that low glycogen content has a detrimental effect, whereas high glycogen content does not improve performance during continuous intense exercise. The importance of the muscle glycogen content in relation to endurance exercise performance is well established (Allen et al., 2008) and is exemplified by a strong correlation between the cycling capacity (75% of
4.2 Glycogen content and muscle contractile characteristics
Reduced glycogen content may impair muscle function at major ATP-dependent sites, including myosin ATPases, SR Ca2+ ATPases, and NKAs (Nielsen et al., 2022). However, from a mechanistic perspective, how low glycogen content impairs intense exercise capacity remains unclear. A glycogen content of less than approximately 300 mmol·kgd.w.−1 can impair in vitro Ca2+ release rates in SR vesicles (Ortenblad et al., 2011) and reduce the tetanic cytosolic Ca2+ concentration and force output (Chin and Allen, 1997). Our findings that the MVC and electrically induced twitches in a non-fatigued human major muscle group are unaffected by severely depressed muscle glycogen content suggest that sufficient Ca2+ is released for maximal activation of the contractile components, indicating that in vivo rates of SR Ca2+ release and reuptake are not impaired, at least in non-fatigued muscles.
At task failure from the intense one-legged knee extensor exercise, the MVC was similar in the two trials and tended to be higher in the glycogen-loaded leg (Table 1) despite the longer exercise duration. Whether this was caused by a faster gradual reduction in MVC in the glycogen-depleted trial or a similarly fast initial MVC reduction in both trials remains unclear. However, we recently reported that the loss of twitch force (Rannou et al., 2019), RFD, and rate of relaxation (Rannou et al., 2021) during a similar exercise protocol occurs gradually, with a fast initial component followed by a slower component. Based on these observations, it is expected that the markedly shorter exercise time in the glycogen-depleted leg is associated with a less pronounced loss of evoked twitch contractile properties (Table 2). Thus, a similar MVC observed at task failure in the present study may be caused by a faster decay of muscle contractile ability in the glycogen-depleted trial, even though the electrically evoked twitch muscle contractile characteristics seem to be dependent on the exercise duration. Noteworthy, the glycogen concentrations were not determined concomitant with the contractile characteristics at task failure from the one-legged exercise, however since subjects did not receive carbohydrates during the experiment, we expected similar differences in glycogen content between the two legs, as observed after the two-legged cycling.
The recovery patterns of MVC, single-twitch maximal force, RFD, and rate of relaxation were similar in the glycogen-loaded and glycogen-depleted legs. However, a significant treatment effect was apparent and appeared to be related to the less compromised muscle function, as determined by electrical stimulations (Table 2) in the glycogen-depleted trial at 60 s of recovery and concomitant similar MVC with a tendency to be higher in the glycogen-loaded leg. The reason for this is unclear but may be related to an expected higher severity of metabolic disturbances with a longer exercise duration (Juel, 1997; Sugaya et al., 2011; Hostrup et al., 2021), as indicated by a slightly lower muscle homogenate pH in the glycogen-loaded trial (Figure 3B). Taken together, these results suggest that the recovery pattern from task failure is unrelated to the prevailing glycogen content.
4.3 Protein abundance and phosphorylation
Manipulating the glycogen content did not alter either the NKAα1, α2, and β1 isoform or FXYD1 protein abundances or the total NKA content determined by [3H]-ouabain binding. Thus, changes in abundance per se of the most important K+-transporting protein (Hostrup et al., 2021) did not play a role in the glycogen-induced difference in muscle performance. Importantly, nonspecific FXYD1 phosphorylation was lower at rest in the glycogen-loaded leg than in the glycogen-depleted leg (Figure 5A). Furthermore, from rest to the end of the two-legged cycling, the nonspecific FXYD1 phosphorylation increased in both legs, as seen previously (Benziane et al., 2011; Thomassen et al., 2011; Thomassen et al., 2013; Thomassen et al., 2016), and this increase tended (P = 0.06) to be augmented in the glycogen-loaded compared to in the glycogen-depleted leg, as indicated by a greater decrease in signal intensity (−0.33 ± 0.07 vs. −0.19 ± 0.07, respectively). In contrast, FXYD1ser68 phosphorylation did not increase (Figure 5B), most likely due to concomitant higher FXYD1ser63 and thr69 phosphorylation (Figures 5C,D), which influenced the binding of the AB_FXYD1ser68 antibody, blunting and underestimating the FXYD1ser68 phosphorylation measurements (Fuller et al., 2009; Thomassen et al., 2016). Thus, augmented FXYD1 phosphorylation may indicate a greater increase in NKA activity (Bibert et al., 2008) during intense exercise in the glycogen-loaded leg, which may have improved K+ and Na+ handling (Gunnarsson et al., 2013; Thomassen et al., 2016) and attenuated the increase in the rate of arm venous [K+]. However, further investigation is warranted, including NKA activity determination. AMPK activity, determined by ACCβser211 phosphorylation (Birk and Wojtaszewski, 2006), was increased by exercise but was unaffected by the glycogen levels (Figures 6A,B) and therefore probably did not differentially affect the NKA activity in the glycogen-loaded and -depleted legs.
The abundance of NHE1 was higher in the depleted leg than in the loaded leg, which may have counteracted the findings of impaired K+ handling. However, a higher abundance of NHE1 is not necessary for better K+ handling and improved performance (Bangsbo et al., 2009; Gunnarsson et al., 2013). Thus, changes in NHE1 abundance did probably not affect the overall effects of glycogen depletion in this study.
4.4 Exercise-induced metabolic disturbances and potential interaction between observations
The small changes in plasma [K+] and [Na+] after intense two-legged cycling and one-legged knee extensor exercise (Table 3) most likely reflect the large gradient between the antecubital venous sampling site and the active muscles as well as the post-exercise sampling time (McKenna et al., 2024). One-legged knee-extensor exercise to exhaustion with a similar duration showed an interstitial [K+] of 8.8–13.7 mM (Nordsborg et al., 2003), while the arm-venous plasma [K+] reached the same concentration in both trials (Table 3) with a shorter exercise duration in the glycogen-depleted trial. A potential consequence of increased interstitial [K+] is reduced muscle cell excitability (Cairns et al., 1995), which has been related to the development of fatigue, especially in type II fibers (Nielsen et al., 2017), due to the intermittent firing of action potentials during trains or completely inexcitable fibers (Lindinger and Cairns, 2021). Thus, far greater increases in interstitial [K+] are expected than the observed increases in venous [K+], and they may reach levels close to where a depressive effect may be seen, especially since intracellular metabolic stress would be expected (Renaud et al., 2023).
The in vivo MVC is assessed over a longer duration (3–4 s) than electrically induced isometric twitches (2–300 ms), involves voluntary tetanic activation, is conducted as trains, and is often initiated by a high-frequency doublet (i.e., two close spikes <5 ms apart) discharge (Van Cutsem et al., 1998). Therefore, the MVC during contractions would generally be more sensitive to disturbances in the ability to regulate excitability than isolated peripheral muscle function. During exercise, the same decrease in MVC was obtained over a longer duration in the glycogen-loaded leg, and in recovery, the MVC in the loaded leg tended (P = 0.06) to be higher, especially immediately after task failure (Table 1). This occurred concomitantly with greater impairments in peripheral contractile muscle function, probably due to more severe metabolic intracellular disturbances (e.g., lower pH; Figure 3B) and higher Pi and lactate, in part given the approximately 65% longer exercise duration at similar exercise intensities (Juel, 1997; Sugaya et al., 2011; Hostrup et al., 2021), which impaired Ca2+ handling and muscle twitch function (Cheng et al., 2018). The attenuated decrease in MVC in the loaded leg may have resulted from improved protection of excitability during the 3–4 s of isometric contractions, partly due to better ion handling in the recruited fibers. Notably, shortly after task failure from intense exercise, the K+ balance and membrane potential may be partially restored, as seen in both humans (Vollestad et al., 1994) and mice, where the resting membrane potential recovers with a time constant of 0.9 and 1.5 min after repetitive stimulation in the soleus and EDL muscles, respectively (Juel, 1986). This might explain why differences in the MVC between the glycogen-loaded and glycogen-depleted legs in this study were most pronounced immediately after exercise (Table 1) but probably also occurred during exercise.
Directly comparing stimulation motor output between humans and rodents is difficult because of size differences and because rodents have approximately 3–4 times higher stimulation frequencies than humans (Manuel et al., 2019). In mouse fast-twitch muscles, initial stimulation frequencies occur above 250 Hz (Gorassini et al., 2000). While in rats during train stimulations, fiber membranes become partially depolarized (Fraser et al., 2011), which increases the repriming period in fast-twitch muscles from 4.1 to 7.8 ms (Dutka and Lamb, 2007). In addition, glycogen depletion can increase the action potential repriming period from 4.5 to 5.3 ms (Watanabe and Wada, 2019) and from 4 (25th percentile:75th percentile, 4:5) to 8 (7:10) ms in rat EDL skinned fibers (Jensen et al., 2020). Thus, the combined effect of exercise-induced membrane depolarization and glycogen depletion may increase the repriming period sufficiently to influence the force outcome induced by stimulation frequencies above approximately 100 Hz (corresponding to 10 ms intervals). This is most crucial in type II fibers with faster glycogen depletion at high intensities (Gollnick et al., 1974) and less protected excitability (Pedersen et al., 2009), despite more pronounced FXYD1 phosphorylation (Thomassen et al., 2013). The electrically induced single twitch force and RFD determined post-exercise are unaffected because of the short time frame of 2–300 ms. In contrast, during the MVC (3–4 s), an increase in the action potential repriming period affects the force production, which is more pronounced during and immediately after intense exercise than during resting conditions. Due to the longer single-legged exercise duration compared to that of the glycogen-depleted leg, the augmented increase (P = 0.06) in nonspecific FXYD1 phosphorylation in the glycogen-loaded leg may induce an even greater nonspecific FXYD1 phosphorylation, which increases over time during continuous exercise (Thomassen et al., 2011) and thus suggests a possible greater NKA activity during exercise in the loaded leg. In support of this, differences in performance, determined using the pedaling frequencies, were only significant in the last part of the 6 s intense exercise sessions in glycogen-depleted participants (Balsom et al., 1999). Furthermore, a pronounced role of glycogen in excitability is observed in vitro in rat skeletal muscle when SR Ca2+ release is mediated by action potentials (Nielsen et al., 2009), while no effect is observed by direct stimulation of the voltage sensors by Na+ depolarization (Goodman et al., 2005), indicating that the coupling between glycogen and fatigue is at least partly due to t-tubular system depolarization. In mouse muscles, a synergistic interaction between increased extracellular [K+] and decreased muscle glycogen content has been demonstrated (Cairns and Renaud, 2023). Interestingly, a high abundance of NKA, especially the NKAα2 isoform, as well as the SR Ca2+ release channels are co-localized within the diffusion-restricted space between the SR and t-tubular system, and in rat skeletal muscles, the NKAs utilize glycogen only from the intramyofibrillar compartment (Nielsen et al., 2022). Furthermore, in humans, high-intensity exercise reveals a compartmentalized glycogen metabolism, including a higher glycogen breakdown rate in type II fibers than in type I fibers, especially in the intra-and intermyofibrillar glycogen fractions (Vigh-Larsen et al., 2022). Together, this indicates that preserved cell membrane excitability, possibly due to higher NKA activity suggested by increased FXYD1 phosphorylation, may play an important role in the improved high-intensity exercise performance in the glycogen-loaded leg; however, this warrants further exploration.
In conclusion, we showed that exercise- and diet-induced low muscle glycogen content impaired continuous high-intensity exercise performance lasting a few minutes. Resting muscle contractile characteristics were unaffected by the muscle glycogen content, whereas at task failure, the electrically-induced single-twitch maximal force, RFD, and rate of relaxation were impaired in the glycogen-loaded leg, probably because of the longer exercise duration. Concomitantly, the MVC was similar and tended to be augmented in the glycogen-loaded leg. Phosphorylation of the NKA activity-regulating FXYD1 tended to be more pronounced in the glycogen-loaded leg during exercise. However, how low glycogen content impairs brief and intense exercise capacity remains unclear.
4.5 Limitations
In this study, muscle glycogen status was manipulated using a combination of dietary control and exercise timing. We observed clear differences in the muscle glycogen content in the legs. However, the glycogen-depleted leg was exercised the evening before the experiment, which may have induced inflammation and increased production of reactive oxygen species, potentially affecting the measurements the following day. Regardless, the influence of exercise timing was not observed at rest on the experimental day, as determined by both the MVC and electrical stimulations. In addition, both legs completed the two-legged intense cycling task before the one-legged knee extensor exercise to task failure, which was used to determine performance.
Furthermore, because of difficulties in obtaining femoral blood samples during cycling, antecubital venous sampling was used. Although this sampling site was non-optimal when investigating leg exercises, any detected differences would have been relevant. As well, due to the many different measurements obtained, sample size and statistical power may be limited for some data, especially for muscle pH and plasma metabolites. Thus, some constraints may be taken in interpretation of blood values, muscle content and muscle contractile characteristics.
Finally, the absolute muscle pH values reported in this study were higher than expected. However, the pH was determined using homogenates at room temperature, which probably induced higher values than intracellular pH at 37°C due to evolution of CO2 in homogenates (Spriet et al., 1986). As well, previous studies have reported individual resting pH values near 7.4 (Mannion et al., 1993; Mannion et al., 1995). Although the reported muscle homogenate pH values were high, the patterns between the legs and from rest to task failure are still valuable.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Ethics Committee of Copenhagen and Frederiksberg communities (KF 01-183/02). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
MT: Data curation, Formal Analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing, Validation, Visualization. MM: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review and editing, Formal Analysis, Validation. HO: Data curation, Formal Analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing, Validation. VW: Data curation, Formal Analysis, Methodology, Writing – original draft, Writing – review and editing, Validation. JB: Conceptualization, Investigation, Writing – original draft, Writing – review and editing, Validation. NN: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing, Validation.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the Danish Ministry of Culture (KIF).
Acknowledgments
We are grateful to Freja Bach Kampmann (Drumm-Hakim), Jonas Salling Quist, and Jens Jung Nielsen who participated in the exercise trials and sample analyses. The results of this study are presented clearly, honestly, and without fabrication, falsification, or inappropriate data manipulation.
Conflict of interest
The 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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Allen D. G., Lamb G. D., Westerblad H. (2008). Skeletal muscle fatigue: cellular mechanisms. Physiol. Rev. 88, 287–332. doi:10.1152/physrev.00015.2007
Andersen P., Adams R. P., Sjogaard G., Thorboe A., Saltin B. (1985). Dynamic knee extension as model for study of isolated exercising muscle in humans. J. Appl. Physiol. 59, 1647–1653. doi:10.1152/jappl.1985.59.5.1647
Balsom P. D., Gaitanos G. C., Soderlund K., Ekblom B. (1999). High-intensity exercise and muscle glycogen availability in humans. Acta Physiol. Scand. 165, 337–345. doi:10.1046/j.1365-201x.1999.00517.x
Bangsbo J., Graham T. E., Kiens B., Saltin B. (1992). Elevated muscle glycogen and anaerobic energy production during exhaustive exercise in man. J. Physiol. (Lond) 451, 205–227. doi:10.1113/jphysiol.1992.sp019161
Bangsbo J., Gunnarsson T. P., Wendell J., Nybo L., Thomassen M. (2009). Reduced volume and increased training intensity elevate muscle Na+-K+ pump alpha2-subunit expression as well as short- and long-term work capacity in humans. J. Appl. Physiol. 107, 1771–1780. doi:10.1152/japplphysiol.00358.2009
Bartlett J. D., Louhelainen J., Iqbal Z., Cochran A. J., Gibala M. J., Gregson W., et al. (2013). Reduced carbohydrate availability enhances exercise-induced p53 signaling in human skeletal muscle: implications for mitochondrial biogenesis. Am. J. Physiol. Regul. Integr. Comp. Physiol. 304, R450–R458. doi:10.1152/ajpregu.00498.2012
Benziane B., Bjornholm M., Pirkmajer S., Austin R. L., Kotova O., Viollet B., et al. (2012). Activation of AMP-activated protein kinase stimulates Na+,K+-ATPase activity in skeletal muscle cells. J. Biol. Chem. 287, 23451–23463. doi:10.1074/jbc.M111.331926
Benziane B., Widegren U., Pirkmajer S., Henriksson J., Stepto N. K., Chibalin A. V. (2011). Effect of exercise and training on phospholemman phosphorylation in human skeletal muscle. Am. J. Physiol. Endocrinol. Metab. 301, E456–E466. doi:10.1152/ajpendo.00533.2010
Bergstrom J., Hermansen L., Hultman E., Saltin B. (1967). Diet, muscle glycogen and physical performance. Acta Physiol. Scand. 71, 140–150. doi:10.1111/j.1748-1716.1967.tb03720.x
Bibert S., Roy S., Schaer D., Horisberger J. D., Geering K. (2008). Phosphorylation of phospholemman (FXYD1) by protein kinases A and C modulates distinct Na,K-ATPase isozymes. J. Biol. Chem. 283, 476–486. doi:10.1074/jbc.M705830200
Birk J. B., Wojtaszewski J. F. (2006). Predominant alpha2/beta2/gamma3 AMPK activation during exercise in human skeletal muscle. J. Physiol. 577, 1021–1032. doi:10.1113/jphysiol.2006.120972
Booth J., McKenna M. J., Ruell P. A., Gwinn T. H., Davis G. M., Thompson M. W., et al. (1997). Impaired calcium pump function does not slow relaxation in human skeletal muscle after prolonged exercise. J. Appl. Physiol. 83, 511–521. doi:10.1152/jappl.1997.83.2.511
Cairns S. P., Flatman J. A., Clausen T. (1995). Relation between extracellular [K+], membrane potential and contraction in rat soleus muscle: modulation by the Na+-K+ pump. Pflugers Arch. 430, 909–915. doi:10.1007/BF01837404
Cairns S. P., Renaud J. M. (2023). The potassium-glycogen interaction on force and excitability in mouse skeletal muscle: implications for fatigue. J. Physiol. 601, 5669–5687. doi:10.1113/jp285129
Cheng A. J., Place N., Westerblad H. (2018). Molecular basis for exercise-induced fatigue: the importance of strictly controlled cellular Ca(2+) handling. Cold Spring Harb. Perspect. Med. 8, a029710. doi:10.1101/cshperspect.a029710
Chidnok W., DiMenna F. J., Fulford J., Bailey S. J., Skiba P. F., Vanhatalo A., et al. (2013). Muscle metabolic responses during high-intensity intermittent exercise measured by (31)P-MRS: relationship to the critical power concept. Am. J. Physiol. Regul. Integr. Comp. Physiol. 305, R1085–R1092. doi:10.1152/ajpregu.00406.2013
Chin E. R., Allen D. G. (1997). Effects of reduced muscle glycogen concentration on force, Ca2+ release and contractile protein function in intact mouse skeletal muscle. J. Physiol. 498 (Pt 1), 17–29. doi:10.1113/jphysiol.1997.sp021838
Coyle E. F., Coggan A. R., Hemmert M. K., Ivy J. L. (1986). Muscle glycogen utilization during prolonged strenuous exercise when fed carbohydrate. J. Appl. Physiol. 61, 165–172. doi:10.1152/jappl.1986.61.1.165
Crambert G., Fuzesi M., Garty H., Karlish S., Geering K. (2002). Phospholemman (FXYD1) associates with Na,K-ATPase and regulates its transport properties. Proc. Natl. Acad. Sci. U. S. A. 99, 11476–11481. doi:10.1073/pnas.182267299
Duhamel T. A., Green H. J., Perco J. G., Ouyang J. (2006a). Comparative effects of a low-carbohydrate diet and exercise plus a low-carbohydrate diet on muscle sarcoplasmic reticulum responses in males. Am. J. Physiol. Cell Physiol. 291, C607–C617. doi:10.1152/ajpcell.00643.2005
Duhamel T. A., Perco J. G., Green H. J. (2006b). Manipulation of dietary carbohydrates after prolonged effort modifies muscle sarcoplasmic reticulum responses in exercising males. Am. J. Physiol. Regul. Integr. Comp. Physiol. 291, R1100–R1110. doi:10.1152/ajpregu.00858.2005
Dutka T. L., Lamb G. D. (2007). Na+-K+ pumps in the transverse tubular system of skeletal muscle fibers preferentially use ATP from glycolysis. Am. J. Physiol. Cell Physiol. 293, C967–C977. doi:10.1152/ajpcell.00132.2007
Fraser J. A., Huang C. L., Pedersen T. H. (2011). Relationships between resting conductances, excitability, and t-system ionic homeostasis in skeletal muscle. J. Gen. Physiol. 138, 95–116. doi:10.1085/jgp.201110617
Fuller W., Eaton P., Bell J. R., Shattock M. J. (2004). Ischemia-induced phosphorylation of phospholemman directly activates rat cardiac Na/K-ATPase. FASEB J. 18 (1), 197–199. doi:10.1096/fj.03-0213fje
Fuller W., Howie J., McLatchie L. M., Weber R. J., Hastie C. J., Burness K., et al. (2009). FXYD1 phosphorylation in vitro and in adult rat cardiac myocytes: threonine 69 is a novel substrate for protein kinase C. Am. J. Physiol. Cell Physiol. 296, C1346–C1355. doi:10.1152/ajpcell.00523.2008
Gaitanos G. C., Williams C., Boobis L. H., Brooks S. (1993). Human muscle metabolism during intermittent maximal exercise. J. Appl. Physiol. 75, 712–719. doi:10.1152/jappl.1993.75.2.712
Gejl K. D., Hvid L. G., Frandsen U., Jensen K., Sahlin K., Ortenblad N. (2014). Muscle glycogen content modifies SR Ca2+ release rate in elite endurance athletes. Med. Sci. Sports Exerc 46, 496–505. doi:10.1249/MSS.0000000000000132
Gollnick P. D., Piehl K., Saltin B. (1974). Selective glycogen depletion pattern in human muscle fibres after exercise of varying intensity and at varying pedalling rates. J. Physiol. 241, 45–57. doi:10.1113/jphysiol.1974.sp010639
Goodman C., Blazev R., Stephenson G. (2005). Glycogen content and contractile responsiveness to T-system depolarization in skinned muscle fibres of the rat. Clin. Exp. Pharmacol. Physiol. 32, 749–756. doi:10.1111/j.1440-1681.2005.04260.x
Gorassini M., Eken T., Bennett D. J., Kiehn O., Hultborn H. (2000). Activity of hindlimb motor units during locomotion in the conscious rat. J. Neurophysiol. 83, 2002–2011. doi:10.1152/jn.2000.83.4.2002
Greenhaff P. L., Gleeson M., Maughan R. J. (1987). The effects of dietary manipulation on blood acid-base status and the performance of high intensity exercise. Eur. J. Appl. Physiol. Occup. Physiol. 56, 331–337. doi:10.1007/BF00690901
Gunnarsson T. P., Christensen P. M., Thomassen M., Nielsen L. R., Bangsbo J. (2013). Effect of intensified training on muscle ion kinetics, fatigue development, and repeated short-term performance in endurance-trained cyclists. Am. J. Physiol. Regul. Integr. Comp. Physiol. 305, R811–R821. doi:10.1152/ajpregu.00467.2012
Harris R. C., Hultman E., Nordesjo L. O. (1974). Glycogen, glycolytic intermediates and high-energy phosphates determined in biopsy samples of musculus quadriceps femoris of man at rest. Methods and variance of values. Scand. J. Clin. Lab. Invest. 33, 109–120. doi:10.1080/00365517409082477
Hawley J. A., Schabort E. J., Noakes T. D., Dennis S. C. (1997). Carbohydrate-loading and exercise performance. An update. Sports Med. 24, 73–81. doi:10.2165/00007256-199724020-00001
Hostrup M., Cairns S. P., Bangsbo J. (2021). Muscle ionic shifts during exercise: implications for fatigue and exercise performance. Compr. Physiol. 11, 1895–1959. doi:10.1002/j.2040-4603.2021.tb00169.x
Ingwersen M. S., Kristensen M., Pilegaard H., Wojtaszewski J. F., Richter E. A., Juel C. (2011). Na,K-ATPase activity in mouse muscle is regulated by AMPK and PGC-1α. J. Membr. Biol. 242, 1–10. doi:10.1007/s00232-011-9365-7
Jensen R., Nielsen J., Ortenblad N. (2020). Inhibition of glycogenolysis prolongs action potential repriming period and impairs muscle function in rat skeletal muscle. J. Physiol. 598, 789–803. doi:10.1113/JP278543
Juel C. (1986). Potassium and sodium shifts during in vitro isometric muscle contraction, and the time course of the ion-gradient recovery. Pflugers Arch. 406, 458–463. doi:10.1007/BF00583367
Juel C. (1997). Lactate-proton cotransport in skeletal muscle. Physiol. Rev. 77, 321–358. doi:10.1152/physrev.1997.77.2.321
Karatzaferi C., de H. A., Ferguson R. A., van M. W., Sargeant A. J. (2001). Phosphocreatine and ATP content in human single muscle fibres before and after maximum dynamic exercise. Pflugers Arch. 442, 467–474. doi:10.1007/s004240100552
Lindinger M. I., Cairns S. P. (2021). Regulation of muscle potassium: exercise performance, fatigue and health implications. Eur. J. Appl. Physiol. 121, 721–748. doi:10.1007/s00421-020-04546-8
Mannion A. F., Jakeman P. M., Willan P. L. (1993). Determination of human skeletal muscle buffer value by homogenate technique: methods of measurement. J. Appl. Physiol. 75, 1412–1418. doi:10.1152/jappl.1993.75.3.1412
Mannion A. F., Jakeman P. M., Willan P. L. (1995). Skeletal muscle buffer value, fibre type distribution and high intensity exercise performance in man. Exp. Physiol. 80, 89–101. doi:10.1113/expphysiol.1995.sp003837
Manuel M., Chardon M., Tysseling V., Heckman C. J. (2019). Scaling of motor output, from mouse to humans. Physiol. (Bethesda) 34, 5–13. doi:10.1152/physiol.00021.2018
McKenna M. J., Renaud J. M., Ortenblad N., Overgaard K. (2024). A century of exercise physiology: effects of muscle contraction and exercise on skeletal muscle Na(+),K(+)-ATPase, Na(+) and K(+) ions, and on plasma K(+) concentration-historical developments. Eur. J. Appl. Physiol. 124, 681–751. doi:10.1007/s00421-023-05335-9
Nielsen J., Dubillot P., Stausholm M. H., Ortenblad N. (2022). Specific ATPases drive compartmentalized glycogen utilization in rat skeletal muscle. J. Gen. Physiol. 154. doi:10.1085/jgp.202113071
Nielsen J., Schroder H. D., Rix C. G., Ortenblad N. (2009). Distinct effects of subcellular glycogen localization on tetanic relaxation time and endurance in mechanically skinned rat skeletal muscle fibres. J. Physiol. 587, 3679–3690. doi:10.1113/jphysiol.2009.174862
Nielsen O. B., de Paoli F. V., Riisager A., Pedersen T. H. (2017). Chloride channels take center stage in acute regulation of excitability in skeletal muscle: implications for fatigue. Physiol. (Bethesda) 32, 425–434. doi:10.1152/physiol.00006.2015
Nordsborg N., Mohr M., Pedersen L. D., Nielsen J. J., Langberg H., Bangsbo J. (2003). Muscle interstitial potassium kinetics during intense exhaustive exercise: effect of previous arm exercise. Am. J. Physiol. Regul. Integr. Comp. Physiol. 285, R143–R148. doi:10.1152/ajpregu.00029.2003
Norgaard A., Kjeldsen K., Clausen T. (1984). A method for the determination of the total number of 3H-ouabain binding sites in biopsies of human skeletal muscle. Scand. J. Clin. Lab. Invest. 44, 509–518. doi:10.3109/00365518409083604
Ortenblad N., Nielsen J., Saltin B., Holmberg H. C. (2011). Role of glycogen availability in sarcoplasmic reticulum Ca2+ kinetics in human skeletal muscle. J. Physiol. 589, 711–725. doi:10.1113/jphysiol.2010.195982
Pedersen T. H., Macdonald W. A., de Paoli F. V., Gurung I. S., Nielsen O. B. (2009). Comparison of regulated passive membrane conductance in action potential-firing fast- and slow-twitch muscle. J. General Physiology 134, 323–337. doi:10.1085/jgp.200910291
Pernow B., Saltin B. (1971). Availability of substrates and capacity for prolonged heavy exercise in man. J. Appl. Physiology 31, 416–422. doi:10.1152/jappl.1971.31.3.416
Petersen A. C., Leikis M. J., McMahon L. P., Kent A. B., Murphy K. T., Gong X., et al. (2012). Impaired exercise performance and muscle Na(+),K(+)-pump activity in renal transplantation and haemodialysis patients. Nephrol. Dial. Transpl. 27, 2036–2043. doi:10.1093/ndt/gfr586
Pilegaard H., Keller C., Steensberg A., Helge J. W., Pedersen B. K., Saltin B., et al. (2002). Influence of pre-exercise muscle glycogen content on exercise-induced transcriptional regulation of metabolic genes. J. Physiol. 541, 261–271. doi:10.1113/jphysiol.2002.016832
Rannou F., Nybo L., Andersen J. E., Nordsborg N. B. (2019). Monitoring muscle fatigue progression during dynamic exercise. Med. Sci. Sports Exerc 51, 1498–1505. doi:10.1249/MSS.0000000000001921
Rannou F., Nybo L., Andersen J. E., Nordsborg N. B. (2021). Muscle contractile characteristics during exhaustive dynamic exercise and recovery. Front. Physiol. 12, 660099. doi:10.3389/fphys.2021.660099
Renaud J. M., Ortenblad N., McKenna M. J., Overgaard K. (2023). Exercise and fatigue: integrating the role of K(+), Na(+) and Cl(-) in the regulation of sarcolemmal excitability of skeletal muscle. Eur. J. Appl. Physiol. 123, 2345–2378. doi:10.1007/s00421-023-05270-9
Rockwell M. S., Rankin J. W., Dixon H. (2003). Effects of muscle glycogen on performance of repeated sprints and mechanisms of fatigue. Int. J. Sport Nutr. Exerc Metab. 13, 1–14. doi:10.1123/ijsnem.13.1.1
Spriet L. L., Soderlund K., Thomson J. A., Hultman E. (1986). pH measurement in human skeletal muscle samples: effect of phosphagen hydrolysis. J. Appl. Physiol. 61, 1949–1954. doi:10.1152/jappl.1986.61.5.1949
Sugaya M., Yasuda T., Suga T., Okita K., Abe T. (2011). Change in intramuscular inorganic phosphate during multiple sets of blood flow-restricted low-intensity exercise. Clin. Physiol. Funct. Imaging 31, 411–413. doi:10.1111/j.1475-097X.2011.01033.x
Thomassen M., Christensen P. M., Gunnarsson T. P., Nybo L., Bangsbo J. (2010). Effect of 2-wk intensified training and inactivity on muscle Na+-K+ pump expression, phospholemman (FXYD1) phosphorylation, and performance in soccer players. J. Appl. Physiology 108, 898–905. doi:10.1152/japplphysiol.01015.2009
Thomassen M., Gunnarsson T. P., Christensen P. M., Pavlovic D., Shattock M. J., Bangsbo J. (2016). Intensive training and reduced volume increases muscle FXYD1 expression and phosphorylation at rest and during exercise in athletes. Am. J. Physiol. Regul. Integr. Comp. Physiol. 310, R659–R669. doi:10.1152/ajpregu.00081.2015
Thomassen M., Murphy R. M., Bangsbo J. (2013). Fibre type-specific change in FXYD1 phosphorylation during acute intense exercise in humans. J. Physiol. 591, 1523–1533. doi:10.1113/jphysiol.2012.247312
Thomassen M., Rose A. J., Jensen T. E., Maarbjerg S. J., Bune L., Leitges M., et al. (2011). Protein kinase Cα activity is important for contraction-induced FXYD1 phosphorylation in skeletal muscle. Am. J. Physiol. Regul. Integr. Comp. Physiol. 301, R1808–R1814. doi:10.1152/ajpregu.00066.2011
Van Cutsem M., Duchateau J., Hainaut K. (1998). Changes in single motor unit behaviour contribute to the increase in contraction speed after dynamic training in humans. J. Physiol. 513 (Pt 1), 295–305. doi:10.1111/j.1469-7793.1998.295by.x
Vigh-Larsen J. F., Ortenblad N., Emil Andersen O., Thorsteinsson H., Kristiansen T. H., Bilde S., et al. (2022). Fibre type- and localisation-specific muscle glycogen utilisation during repeated high-intensity intermittent exercise. J. Physiol. 600, 4713–4730. doi:10.1113/jp283225
Vigh-Larsen J. F., Ortenblad N., Spriet L. L., Overgaard K., Mohr M. (2021). Muscle glycogen metabolism and high-intensity exercise performance: a narrative review. Sports Med. 51, 1855–1874. doi:10.1007/s40279-021-01475-0
Vollestad N. K., Hallen J., Sejersted O. M. (1994). Effect of exercise intensity on potassium balance in muscle and blood of man. J. Physiol. (Lond) 475, 359–368. doi:10.1113/jphysiol.1994.sp020077
Watanabe D., Wada M. (2019). Effects of reduced muscle glycogen on excitation-contraction coupling in rat fast-twitch muscle: a glycogen removal study. J. Muscle Res. Cell Motil. 40, 353–364. doi:10.1007/s10974-019-09524-y
Wojtaszewski J. F., MacDonald C., Nielsen J. N., Hellsten Y., Hardie D. G., Kemp B. E., et al. (2003). Regulation of 5’AMP-activated protein kinase activity and substrate utilization in exercising human skeletal muscle. Am. J. Physiol. Endocrinol. Metab. 284, E813–E822. doi:10.1152/ajpendo.00436.2002
Keywords: MVC, potassium, Na+/K+-ATPase, FXYD1, excitability
Citation: Thomassen M, McKenna MJ, Olmedillas H, Wyckelsma V, Bangsbo J and Nordsborg NB (2025) Exercise- and diet-induced glycogen depletion impairs performance during one-legged constant-load, high-intensity exercise in humans. Front. Physiol. 16:1564523. doi: 10.3389/fphys.2025.1564523
Received: 21 January 2025; Accepted: 25 July 2025;
Published: 15 August 2025.
Edited by:
Gregory C. Bogdanis, National and Kapodistrian University of Athens, GreeceReviewed by:
Llion Arwyn Roberts, Griffith University, AustraliaJeffrey J. Brault, Indiana University Bloomington, United States
Copyright © 2025 Thomassen, McKenna, Olmedillas, Wyckelsma, Bangsbo and Nordsborg. 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: Martin Thomassen, bXRob21hc3NlbkBuZXhzLmt1LmRr