Advertisement
Clinical Study| Volume 22, ISSUE 4, P616-628, April 2022

Download started.

Ok

Intramuscular lipid concentration increased in localized regions of the lumbar muscles following 60 day bedrest

Open AccessPublished:November 20, 2021DOI:https://doi.org/10.1016/j.spinee.2021.11.007

      Highlights

      • Lipid accumulated in lumbar multifidus and erector spinae after bedrest
      • Lipid accumulation is distributed inhomogeneously in lumbar muscles after bedrest
      • Lipid accumulated more in medial and/or lateral than central region of lumbar multifidus
      • Inhomogeneous spatial lipid accumulation may influence lumbar spine function

      Abstract

      BACKGROUND CONTEXT

      Prolonged bedrest induces accumulation of intramuscular lipid concentration (ILC) in the lumbar musculature; however, spatial distribution of ILC has not been determined. Artificial gravity (AG) mitigates some adaptations induced by 60 day bedrest by creating a head-to-feet force while participants are in a supine position.

      PURPOSE

      To quantify the spatial distribution of accumulation of ILC in the lumbar musculature after 60 day bedrest, and whether this can be mitigated by AG exposure.

      STUDY DESIGN

      Prospective longitudinal study.

      PATIENT SAMPLE

      Twenty-four healthy individuals (8 females) participated in the study: Eight received 30 min continuous AG (cAG); Eight received 6 × 5 min AG (iAG), interspersed with rests; Eight were not exposed to AG (CRTL).

      OUTCOME MEASURES

      From 3T magnetic resonance imaging (MRI), axial images were selected to assess lumbar multifidus (LM), lumbar erector spinae (LES), quadratus lumborum (QL), and psoas major (PM) muscles from L1/L2 to L5/S1 intervertebral disc levels. Chemical shift-based 2‐echo lipid and/or water Dixon sequence was used to measure tissue composition. Each lumbar muscle was segmented into four equal quartiles (from medial to lateral).

      METHODS

      Participants arrived at the facility for the baseline data collection before undergoing a 60 day strict 6° head-down tilt (HDT) bedrest period. MRI of the lumbopelvic region was conducted at baseline and Day-59 of bedrest. Participants performed all activities, including hygiene, in 6° HDT and were discouraged from moving excessively or unnecessarily.

      RESULTS

      At the L4/L5 and L5/S1 intervertebral disc levels, 60-day bedrest induced a greater increase in ILC in medial and lateral regions (∼+4%) of the LM than central regions (∼+2%; p<.05). A smaller increase in ILC was induced in the lateral region of LES (∼+1%) at L1/L2 and L2/L3 than at the centro-medial region (∼+2%; p<.05). There was no difference between CRTL and intervention groups.

      CONCLUSIONS

      Inhomogeneous spatial distribution of accumulation of ILC was found in the lumbar musculature after 60 day bedrest. These findings might reflect pathophysiological mechanisms related to muscle disuse and contribute to localized lumbar spine dysfunction. Altered spatial distribution of ILC may impair lumbar spine function after prolonged body unloading, which could increase injury risk to vulnerable soft tissues, such as the lumbar intervertebral discs. These novel results may represent a new biomarker of lumbar deconditioning for astronauts, bedridden, sedentary individuals, or those with chronic back pain. Changes are potentially modifiable but not by the AG protocols tested here.

      Keywords

      Introduction

      Movement and gravitational loading appear necessary to preserve spinal health. When vertical loading of the human body is reduced, such as experienced during prolonged bedrest or spaceflight, the lumbar spine undergoes rapid remodelling [
      • Bailey JF
      • Miller SL
      • Khieu K
      • Neill CWO
      • Healey RM
      • Coughlin DG
      • et al.
      From the international space station to the clinic : how prolonged unloading may disrupt lumbar spine stability.
      ,
      • Chang DG
      • Healey RM
      • Snyder AJ
      • Sayson JV
      • Macias BR
      • Coughlin DG
      • et al.
      Lumbar spine paraspinal muscle and intervertebral disc height changes in astronauts after long-duration spaceflight on the International Space Station.
      ], which may increase the risk of low back pain or injury [
      • Belavý DL
      • Adams M
      • Brisby H
      • Cagnie B
      • Danneels L
      • Fairbank J
      • et al.
      Disc herniations in astronauts: What causes them, and what does it tell us about herniation on earth?.
      ].
      Although morphological changes of the lumbar spine (eg, reduced muscle size, reduced lumbar lordosis) after prolonged vertical unloading have been extensively investigated [
      • Hides JA
      • Lambrecht G
      • Sexton CT
      • Pruett C
      • Petersen N
      • Jaekel P
      • et al.
      The effects of exposure to microgravity and reconditioning of the lumbar multifidus and anterolateral abdominal muscles: implications for people with LBP.
      ,
      • Belavý DL
      • Armbrecht G
      • Richardson CA
      • Felsenberg D
      • Hides JA
      Muscle atrophy and changes in spinal morphology: Is the lumbar spine vulnerable after prolonged bed-rest?.
      ,
      • Burkhart K
      • Allaire B
      • Bouxsein ML
      Negative Effects of Long-duration Spaceflight on Paraspinal Muscle Morphology.
      ], changes in muscle composition, such as intramuscular lipid concentration (ILC), have received less attention. Increased ILC has been reported in the lumbar multifidus (LM) muscle of individuals with persistent low back pain (LBP) [
      • Hodges PW
      • Danneels L
      Changes in Structure and Function of the Back Muscles in Low Back Pain: Different Time Points, Observations, and Mechanisms.
      ,
      • Kjaer P
      • Bendix T
      • Sorensen JS
      • Korsholm L
      • Leboeuf-Yde C
      Are MRI-defined fat infiltrations in the multifidus muscles associated with low back pain?.
      ,
      • Hildebrandt M
      • Fankhauser G
      • Meichtry A
      • Luomajoki H
      Correlation between lumbar dysfunction and fat infiltration in lumbar multifidus muscles in patients with low back pain.
      ], in adults with sway-back posture [
      • Pezolato A
      • De Vasconcelos EE
      • Defino HLA
      • Defino A
      • Nogueira-Barbosa MH
      Fat infiltration in the lumbar multifidus and erector spinae muscles in subjects with sway-back posture.
      ], elderly individuals [
      • Dahlqvist JR
      • Vissing CR
      • Hedermann G
      • Thomsen C
      • Vissing J
      Fat Replacement of Paraspinal Muscles with Aging in Healthy Adults.
      ,
      • Crawford RJ
      • Filli L
      • Elliott JM
      • Nanz D
      • Fischer MA
      • Marcon M
      • et al.
      Age- and level-dependence of fatty infiltration in lumbar paravertebral muscles of healthy volunteers.
      ,
      • Crawford RJ
      • Volken ÃT
      • Elliott J
      • Hoggarth MA
      • Dpt MS
      • Samartzis D
      Geography of Lumbar Paravertebral Muscle.
      ], and astronauts after spaceflights [
      • McNamara KP
      • Greene KA
      • Moore AM
      • Lenchik L
      • Weaver AA
      Lumbopelvic muscle changes following long-duration spaceflight.
      ]. Animal studies have shown time-dependent lipid accumulation after intervertebral disc injury [
      • Hodges PW
      • James G
      • Blomster L
      • Hall L
      • Schmid A
      • Shu C
      • et al.
      Multifidus Muscle Changes after Back Injury Are Characterized by Structural Remodeling of Muscle, Adipose and Connective Tissue, but Not Muscle Atrophy: Molecular and Morphological Evidence.
      ,
      • Hodges PW
      • James G
      • Blomster L
      • Hall L
      • Schmid AB
      • Shu C
      • et al.
      Can proinflammatory cytokine gene expression explain multifidus muscle fiber changes after an intervertebral disc lesion?.
      ], also supported by human data [
      • Teichtahl AJ
      • Urquhart DM
      • Wang Y
      • Wluka AE
      • Wijethilake P
      • O'Sullivan R
      • et al.
      Fat infiltration of paraspinal muscles is associated with low back pain, disability, and structural abnormalities in community-based adults.
      ]. It has been proposed that the accumulation of intramuscular lipid may impact the capacity to meet functional demands to control the spine [
      • Hodges PW
      • James G
      • Blomster L
      • Hall L
      • Schmid A
      • Shu C
      • et al.
      Multifidus Muscle Changes after Back Injury Are Characterized by Structural Remodeling of Muscle, Adipose and Connective Tissue, but Not Muscle Atrophy: Molecular and Morphological Evidence.
      ,
      • Teichtahl AJ
      • Urquhart DM
      • Wang Y
      • Wluka AE
      • Wijethilake P
      • O'Sullivan R
      • et al.
      Fat infiltration of paraspinal muscles is associated with low back pain, disability, and structural abnormalities in community-based adults.
      ].
      The spatial distribution of ILC in the LM and lumbar erector spinae (LES) muscles appears dependent on the vertebral level (greater accumulation at the lower lumbar vertebral levels [
      • Crawford RJ
      • Volken ÃT
      • Elliott J
      • Hoggarth MA
      • Dpt MS
      • Samartzis D
      Geography of Lumbar Paravertebral Muscle.
      ,
      • Fortin M
      • Lazáry À
      • Varga PP
      • Battié MC
      Association between paraspinal muscle morphology, clinical symptoms and functional status in patients with lumbar spinal stenosis.
      ,
      • Lee SH
      • Park SW
      • Kim YB
      • Nam TK
      • Lee YS
      The fatty degeneration of lumbar paraspinal muscles on computed tomography scan according to age and disc level.
      ]) and also on the muscle location, where the most medial regions, adjacent to the spinous processes, have proportionally greater ILC [
      • Crawford RJ
      • Volken ÃT
      • Elliott J
      • Hoggarth MA
      • Dpt MS
      • Samartzis D
      Geography of Lumbar Paravertebral Muscle.
      ,
      • Smith AC
      • Albin SR
      • Abbott R
      • Crawford RJ
      • Hoggarth MA
      • Wasielewski M
      • et al.
      Confirming the geography of fatty infiltration in the deep cervical extensor muscles in whiplash recovery.
      ,
      • Abbott R
      • Peolsson A
      • West J
      • Elliott JM
      • Åslund U
      • Karlsson A
      • et al.
      The qualitative grading of muscle fat infiltration in whiplash using fat and water magnetic resonance imaging.
      ]. The spatial distribution of accumulation of ILC in lumbar musculature induced by prolonged disuse remains unclear. Investigation of such changes is necessary to guide the design of prevention and rehabilitation protocols. Importantly, as paraspinal muscles have a complex architectural anatomy [
      • Hodges PW
      • Danneels L
      Changes in Structure and Function of the Back Muscles in Low Back Pain: Different Time Points, Observations, and Mechanisms.
      ,
      • Macintosh JE
      • Valencia F
      • Bogduk N
      • Munro R
      The morphology multifidus of the human lumbar.
      ,
      • Macintosh JE
      • Bogduk N
      The morphology of the lumbar erector spinae.
      ] and because ILC has been shown to accumulate in the muscles of the lower lumbar spine [
      • Crawford RJ
      • Filli L
      • Elliott JM
      • Nanz D
      • Fischer MA
      • Marcon M
      • et al.
      Age- and level-dependence of fatty infiltration in lumbar paravertebral muscles of healthy volunteers.
      ,
      • Crawford RJ
      • Volken ÃT
      • Elliott J
      • Hoggarth MA
      • Dpt MS
      • Samartzis D
      Geography of Lumbar Paravertebral Muscle.
      ,
      • Lee SH
      • Park SW
      • Kim YB
      • Nam TK
      • Lee YS
      The fatty degeneration of lumbar paraspinal muscles on computed tomography scan according to age and disc level.
      ], a detailed analysis of the muscles at all lumbar vertebral levels is necessary to understand any potential changes due to prolonged disuse.
      Strict head-down tilt bedrest (HDTBR) have been used to investigate the effect of prolonged vertical unloading upon the human body [
      • Hargens AR
      • Vico L
      Long-duration bed rest as an analog to microgravity.
      ]. This position unloads the body's upright weight [
      • Hargens AR
      • Vico L
      Long-duration bed rest as an analog to microgravity.
      ], reduces energy requirements [
      • Pavy-Le Traon A
      • Heer M
      • Narici MV
      • Rittweger J
      • Vernikos J
      From space to Earth: Advances in human physiology from 20 years of bed rest studies.
      ] and overall sensory stimulation [
      • Pavy-Le Traon A
      • Heer M
      • Narici MV
      • Rittweger J
      • Vernikos J
      From space to Earth: Advances in human physiology from 20 years of bed rest studies.
      ]. In contrast, artificial gravity (AG), via centrifugation in a short-arm centrifuge, has been suggested to mitigate many of these effects by stimulating the proprioceptive, vestibular, and neuroendocrine systems [
      • Clément GR
      • Bukley AP
      • Paloski WH
      Artificial gravity as a countermeasure for mitigating physiological deconditioning during long-duration space missions.
      ,
      • Kotovskaya AR
      The problem of artificial gravity: The current state and prospects.
      ]. As AG is associated with a large acceleration gradient along the body axes (∼1Gz at the lumbar vertebral level), it has also been proposed that a mechanical compressive force applied to the lumbar spine could stimulate the lumbar muscles [
      • De Martino E
      • Hides J
      • Elliott JM
      • Hoggarth M
      • Zange J
      • Lindsay K
      • et al.
      Lumbar muscle atrophy and increased relative intramuscular lipid concentration are not mitigated by daily artificial gravity following 60-day head-down tilt bed rest.
      ]. It is unknown whether the increased ILC in the lumbar musculature in response to exposure to prolonged HDTBR is homogeneous within the muscles and spinal levels or with spatial variation of accumulation, or whether this can be mitigated by exposure to AG. This study aimed to investigate the spatial distribution of lipid accumulation patterns in the lumbar musculature in response to HDTBR. Furthermore, we aimed to examine whether AG, either continuously or intermittent, could mitigate the localized increases in ILC.

      Methods

      Study Design

      The study was conducted at the ":envihab" facility in Cologne (Germany) as part of the Artificial Gravity Bed Rest—European Space Agency (AGBRESA) study. Twenty-four healthy participants (16 males) were recruited and allocated to one of three groups: (1) Control (CTRL) group (N=8); (2) 30 minutes continuous AG (cAG) daily (N=8); (3) intermittent 6 × 5 minutes AG (iAG) daily (N=8). The sex, age, height, and weight of the participant groups were comparable (CRTL – 2 females, 34±8 years, 177±7 cm, 79±13 kg; cAG – 3 females, 32±10 years, 173±8 cm, 72±10 kg; iAG – 3 females, 34±11 years, 174±11 cm, 71±5 kg). The study size of 24 was selected based on previous HDTBR studies showing the protective effects of AG on bone resorption [
      • Rittweger J
      • Bareille MP
      • Clément G
      • Linnarsson D
      • Paloski WH
      • Wuyts F
      • et al.
      Short-arm centrifugation as a partially effective musculoskeletal countermeasure during 5-day head-down tilt bed rest—results from the BRAG1 study.
      ] and orthostatic tolerance [
      • Linnarsson D
      • Hughson RL
      • Fraser KS
      • Clément G
      • Karlsson LL
      • Mulder E
      • et al.
      Effects of an artificial gravity countermeasure on orthostatic tolerance, blood volumes and aerobic power after short-term bed rest (BR-AG1).
      ]. Other measures from this study have been published elsewhere and showed that HDTBR induced increases in spinal length and area of lumbar intervertebral discs, a reduction in the lumbar lordosis, and atrophy of the LM, LES, and QL muscles [
      • De Martino E
      • Hides J
      • Elliott JM
      • Hoggarth M
      • Zange J
      • Lindsay K
      • et al.
      Lumbar muscle atrophy and increased relative intramuscular lipid concentration are not mitigated by daily artificial gravity following 60-day head-down tilt bed rest.
      ].
      The current study was performed in accordance with the International Guidelines for Standardization of bedrest in the spaceflight context [
      • Sundblad P
      • Orlov O
      • Angerer O
      • Larina I
      • Cromwell R
      Guidelines for standardization of bed rest studies in the spaceflight context.
      ]. The study consisted of a baseline data collection (BDC) period, followed by 60 days of strict 6° HDTBR period and 14-day rehabilitation [
      • Ganse B
      • Bosutti A
      • Drey M
      • Degens H
      Sixty days of head-down tilt bed rest with or without artificial gravity do not affect the neuromuscular secretome.
      ]. Participants carried out all tasks, including hygiene, in a supine posture and were discouraged from making unnecessary movements [
      • Ganse B
      • Bosutti A
      • Drey M
      • Degens H
      Sixty days of head-down tilt bed rest with or without artificial gravity do not affect the neuromuscular secretome.
      ]. They were allowed to lie supine or on their side but were advised to have at least one shoulder touching the bed at all times [
      • Frett T
      • Green DA
      • Mulder E
      • Noppe A
      • Arz M
      • Pustowalow W
      • et al.
      Tolerability of daily intermittent or continuous short-arm centrifugation during 60-day 6o head down bed rest (AGBRESA study).
      ]. Participants received a controlled diet with a daily caloric intake of 1.3 times the metabolic rate at rest and were expected to eat all the food they were served [
      • Kramer A
      • Venegas-Carro M
      • Mulder E
      • Lee JK
      • Moreno-Villanueva M
      • Bürkle A
      • et al.
      Cardiorespiratory and Neuromuscular Demand of Daily Centrifugation: Results From the 60-Day AGBRESA Bed Rest Study.
      ]. The Ethics Committee of the Northern Rhine Medical Association approved this study (Düsseldorf, Germany, Application No. 2018143), and participants provided written informed consent to participate in the study. The study was registered at the German Clinical Trial Register under No. DKRS00015677.

      Artificial gravity

      Transfer to the centrifuge was accomplished with a specific tilt gurney, and participants were placed on the 3 meter centrifuge arm in supine position (6° head-down tilt). During centrifugation, participants were exposed to 1Gz at their estimated center of mass. They could perform anti-orthostatic maneuvers, such as heel raises and shallow knee bends, to avoid calf pain and maintain circulation while spinning but were otherwise instructed to remain still [
      • Frett T
      • Green DA
      • Mulder E
      • Noppe A
      • Arz M
      • Pustowalow W
      • et al.
      Tolerability of daily intermittent or continuous short-arm centrifugation during 60-day 6o head down bed rest (AGBRESA study).
      ]. Continuous medical monitoring to ensure participant safety was implemented during all centrifuge runs [
      • Frett T
      • Green DA
      • Mulder E
      • Noppe A
      • Arz M
      • Pustowalow W
      • et al.
      Tolerability of daily intermittent or continuous short-arm centrifugation during 60-day 6o head down bed rest (AGBRESA study).
      ].

      MRI Measurements

      MRIs were collected using a 3 Tesla Magnetom Vision system (Siemens, Erlangen, Germany). Participants were positioned in supine lying on the scanning table with their knees and hips supported in slight flexion by a pillow. MRIs were acquired two days before HDTBR (BDC) and on the 59th day of HDTBR (HDT59) and stored for offline analysis. A set of 64 transverse images were acquired from the T12 vertebra to the sacrum (T1 weighted Dixon sequence, total acquisition time=5 minutes; slice thickness=4 mm; distance factor=20%, TR=7.02 ms, TE1=2.46 ms, TE2=3.69 ms, flip angle=5 deg; field of view=400 mm x 400 mm at 1.0 mm x 1.0 mm pixel size). Images were obtained with the fat and water in-phase and out of phase; then, fat (F) and water images (W) were reconstructed. Regions of interest (ROI) were manually traced over the lumbar paravertebral muscles using a semi-automated Matlab-based program [
      • Crawford RJ
      • Volken ÃT
      • Elliott J
      • Hoggarth MA
      • Dpt MS
      • Samartzis D
      Geography of Lumbar Paravertebral Muscle.
      ,
      • Mhuiris ÁN
      • Volken T
      • Elliott JM
      • Hoggarth M
      • Samartzis D
      • Crawford RJ
      Reliability of quantifying the spatial distribution of fatty infiltration in lumbar paravertebral muscles using a new segmentation method for T1-weighted MRI.
      ]. The custom built Matlab (Natick, MA, USA) program automatically divided the ROI into quarters of equal area from medial (Q1) to lateral (Q4) based on the pixel number (Fig. 1) [
      • Crawford RJ
      • Volken ÃT
      • Elliott J
      • Hoggarth MA
      • Dpt MS
      • Samartzis D
      Geography of Lumbar Paravertebral Muscle.
      ,
      • Abbott R
      • Peolsson A
      • West J
      • Elliott JM
      • Åslund U
      • Karlsson A
      • et al.
      The qualitative grading of muscle fat infiltration in whiplash using fat and water magnetic resonance imaging.
      ]. ILC was calculated as the ratio of pixel intensities from the F and W images:
      ILC=F(W+F)*100


      Fig 1
      Fig. 1Sagittal projection of the lumbar spine. B: Characteristic location of lumbar paraspinal muscles identified for the measurement on axial images at BDC (average values of 24 participants). ILC was automatically quartiled based on equal pixel numbers. The colour scale represents the percentage of fat infiltration (0%-52%).
      Bilateral ILC measurements of the LM, LES, QL, and PM muscles were extracted from each transverse MRI [
      • Crawford RJ
      • Volken ÃT
      • Elliott J
      • Hoggarth MA
      • Dpt MS
      • Samartzis D
      Geography of Lumbar Paravertebral Muscle.
      ,
      • Abbott R
      • Peolsson A
      • West J
      • Elliott JM
      • Åslund U
      • Karlsson A
      • et al.
      The qualitative grading of muscle fat infiltration in whiplash using fat and water magnetic resonance imaging.
      ]. Four slices were identified for each of five lumbar intervertebral discs (L1/L2, L2/L3, L3/L4, L4/L5, L5/S1; Fig. 1). The ILC measurements for the MRIs were averaged for the 4 slices at each lumbar region and the left- and right-side. The changes were calculated as a difference to the BDC (HDT59 value – BDC value) for each muscle quartile for the statistical analysis. Reliability of quantification of fat distribution in the lumbar paravertebral muscles in the transverse plane has been demonstrated in humans [
      • Mhuiris ÁN
      • Volken T
      • Elliott JM
      • Hoggarth M
      • Samartzis D
      • Crawford RJ
      Reliability of quantifying the spatial distribution of fatty infiltration in lumbar paravertebral muscles using a new segmentation method for T1-weighted MRI.
      ]. This ILC evaluation technique has been validated in pig and rabbit models using the reference standard biopsy and/or histology [
      • Smith AC
      • Parrish TB
      • Abbott R
      • Hoggarth MA
      • Mendoza K
      • Chen YF
      • et al.
      Muscle-fat MRI: 1.5 Tesla and 3.0 Tesla versus histology.
      ].

      Statistical analysis

      Statistical analysis was undertaken using SPSS (Version 25, IBM, Chicago, USA). All results are presented as mean (standard deviation, SD). Statistical significance was set at the (2-sided) 0.05 level. Outcomes were assessed for normality using visual inspection (histograms and Q–Q plots) and Shapiro–Wilk tests. First, a two-way Analysis of Variance (ANOVA) was used to examine the distribution of ILC throughout the lumbar paraspinal muscles at BDC (before HDTBR) using Quartile (Q1, Q2, Q3, and Q4) and Level (L1/L2, L2/L3, L3/L4, L4/L5, and L5/S1) as within-group factors. Changes in response to HDTBR were assessed with three-way ANOVA using time (BDC and HDT59) and Quartile (Q1, Q2, Q3, and Q4) as repeated measures and Groups (iAG, cAG, CTRL) as between-group factors. Since the percentage of ILC in the quartiles was different in all muscles and all vertebral levels, a mixed-model ANOVA was used for change in ILC after HDTBR (HDT59 value – BDC value) between Groups (CRTL, cAG, iAG; between-group factor) and Quartile (Q1, Q2, Q3, and Q4; within-subject factor) to specifically investigate the spatial variation of lipid accumulation. The interaction effect of Group and Quartile was included in all models. The Greenhouse–Geisser approach was used to correct against violations of sphericity. Effect sizes (partial eta-squared: η2partial) were reported. Post hoc pairwise comparisons were performed using Bonferroni corrected multiple comparisons when significant main effect or interaction and corresponding 95% confidence intervals were generated.

      Results

      Intramuscular lipid concentration at BDC

      Fig. 1 shows the ILC patterns at BDC. For all paraspinal muscles, the 2-way ANOVA showed a main effect of Quartile (all – F3,69 > 35; p<.001) and Level*Quartile interaction (all – F12,276 > 3.5; p<.05), but the effect of Level was only found for LM and LES (both – F4,92 > 40; p<.001).
      For the LM muscle, pairwise comparisons showed higher values in ILC at the L3/L4, L4/L5, and L5/S1 vertebral levels compared with L1/L2 and L2/L3 (all – p<.01), with the highest value at L5/S1 (all – p<.01). Except for L5/S1, the ILC progressively decreased from Q1 (medial) to Q3 (all – p<.001) and increased again from Q3 to Q4 (all – p<.001). This pattern slightly differed at the L5/S1 vertebral level, where the highest values were found in Q1, Q2 and Q4 compared with Q3 (all – p<.01).
      For the LES muscle, pairwise comparisons showed a progressive increase in ILC from the L2/L3 to L5/S1 vertebral levels (all – p<.001), with the highest value at L5/S1 (all – p<.01). Except for L5/S1, the ILC progressively decreased from Q1 (medial) to Q4 (lateral) (all – p<.05). This pattern differed at the L5/S1vertebral level, where the highest values were found in Q1 and Q4 compared with Q2 and Q3 (all – p<.001).
      For the PM muscle, the ILC was higher in Q1 (medial) compared with Q2, Q3, and Q4 (p<.05) for all vertebral levels except for L5/S1, where no difference in quartiles was found. For the QL muscle, ILC was higher in Q1 (medial) and Q4 (lateral) compared with Q2 and Q3 (all - p<.05) for L1/L2 and L3/L4 vertebral levels. This pattern slightly differed at L2/L3 vertebral level, where the highest values were found in Q1 compared with Q2, Q3, and Q4 (all - p<.05).

      Intramuscular lipid concentration in LM muscle after HDTBR

      The changes in ILC at each vertebral level are shown in Fig. 2. At L4/L5 and L5/S1 vertebral levels, the three-way ANOVA showed a main effect of Time (F > 100; p<.001), Quartile (F > 25; p<.001), and Time*Quartile interaction (F > 5; p<.01). No significant main effect of Groups or other interactions were found (data and statistical analysis: Supplementary Table 1 and 2). Analysis of the change in ILC after HDTBR revealed a significant effect of Quartile at the L4/L5 (F2.4,51.2 = 7.1; p=.001) and L5/S1 (F2.0,42.2 = 5.7; p=.006) vertebral levels (Table 1 and Fig. 3). At the level of the L4/L5, pairwise comparisons showed a greater increase in ILC in Q1 than Q3 (p=.044; CI 95% [0.40, 3.99]), Q2 than Q3 (p=.001; CI 95% [0.77, 3.10]) and Q4 than Q3 (p=.001; CI 95% [1.10, 4.88]). At the L5/S1 vertebral level, pairwise comparisons showed a greater increase in ILC in Q4 than Q2 (p=.007; CI 95% [0.80, 7.10]) and Q4 than Q3 (p=.017; CI 95% [0.07, 6.71]). There was no difference between Groups (main effect – F2,21 < 2.5; p>.1) or interaction bewteen Quartile*Groups (all – F3.6,38.0 < 1.5; p>.2).
      Fig 2
      Fig. 2Absolute changes in ILC (HDT59-BDC). The colour scale represents delta ILC (from -2% to 8%). Note: Medial and lateral regions of LM accumulated more adipose tissue at the lower vertebral levels. Less adipose tissue accumulated in the lateral region of LES at the upper vertebral levels.
      Table 1Mean (± standard deviation) of lumbar multifidus change in ILC after HDTBR (HDT59 value – BDC value) from CTRL (N=8), cAG (N=8), and iAG (N=8). Q1 = medial; Q4 = lateral
      Lumbar multifidus
      VariableGroupQuartileMixed model repeated measures ANOVA
      Q1Q2Q3Q4QuartilesGroupsQuartiles*Groups
      L1/L2 intervertebral discCRTL3.3±2.91.2±2.20.8±1.60.5±3.3F1.8,38.0=2.01;

      p=.122

      η2partial=0.09
      F2,21=0.85;

      p=.441

      η2partial=0.75
      F3.6, 38.0=1.41;

      p=.224

      η2partial=0.12
      cAG2.5±2.31.9±2.82.1±2.32.0±2.8
      iAG1.0±2.80.2±2.10.5±2.41.9±4.4
      L2/L3 intervertebral discCRTL2.2±2.92.2±2.21.1±1.71.1±3.9F1.8,39.4=1.26;

      p=.297

      η2partial=0.06
      F2,21=;

      p=

      η2partial=
      F3.7, 39.4=0.95;

      p=.444

      η2partial=0.08
      cAG2.7±1.51.5±1.01.2±2.41.4±3.5
      iAG2.6±1.32.4±1.92.5±2.02.1±3.5
      L3/L4 intervertebral discCRTL2.9±2.12.3±2.11.2±1.31.8±1.3F1.8,36.9=2.53;

      p=.099

      η2partial=0.11
      F2,21=2.42;

      p=.113

      η2partial=0.18
      F3.5,36.9=1.17;

      p=.34

      η2partial=0.10
      cAG5.0±2.72.8±1.31.9±1.45.3±6.1
      iAG2.5±1.62.6±1.72.3±1.62.9±1.7
      L4/L5 intervertebral discCRTL3.6±2.63.2±3.11.7±2.44.2±3.3F3,63=7.1;

      p<.001

      η2partial=0.25
      F2,21=2.30;

      p=.125

      η2partial=0.18
      F6,63=1.08;

      p=.38

      η2partial=0.09
      cAG5.6±1.74.4±2.62.2±2.76.5±3.7
      iAG3.7±1.54.9±2.72.9±2.25.0±2.8
      L5/S1 intervertebral discCRTL4.5±1.84.2±3.33.9±3.66.4±4.9F2.0,42.2=5.74;

      p=.006

      η2partial=0.22
      F2,21=0.32;

      p=.729

      η2partial=0.03
      F4.0,42.2=0.77;

      p=.547

      η2partial=0.07
      cAG6.3±2.54.2±3.53.1±4.47.7±5.6
      iAG4.8±1.32.6±1.73.1±3.27.0±5.8
      Fig 3
      Fig. 3(A) ILC of the lumbar multifidus muscle at BDC and HDT59 (average values of 24 participants) over quartiles (Q1 = medial; Q4 = lateral). (B) Change in ILC after 60 days (HDT59 value – BDC value). * = Significant post hoc test within the Quartile (p<.05).
      At L1/L2, L2/L3 and L3/L4 vertebral levels, the three-way ANOVA showed a main effect of Time (F > 10; p<.005), Quartile (F > 25; p<.001). No significant main effect of Groups or other interactions were found (data and statistical analysis: Supplementary Table 1 and 2).

      Intramuscular lipid concentration in LES muscle after HDTBR

      At L1/L2 and L2/L3 vertebral levels, the three-way ANOVA showed a main effect of Time (F > 30; p<.001), Quartile (F > 8; p<.001), and Time*Quartile interaction (F > 3; p<.05). No significant main effect of Groups or other interactions were found (data and statistical analysis: Supplementary Table 3 and 4). Analysis of the change in ILC after HDTBR revealed a significant effect of Quartile at the levels of the L1/L2 (F2.2,46.3 = 3.3; p=.042) and L2/L3 (F2.0,42.2 = 5.7; p=.006) vertebral levels (Table 2 and Fig. 4). At the L1/L2, there was a greater increased in ILC in Q4 than Q2 (p=.04; CI 95% [0.01, 2.88]). At the L2/L3 level, pairwise comparisons showed a greater increase in ILC in Q1 than Q4 (p=.03; CI 95% [0.06, 1.42]), Q2 than Q4 (p=.01; CI 95% [0.65 2.73]) and Q3 than Q4 (p=.04; CI 95% [0.02, 1.39]). There was no difference between Group (main effect – F2,21 < 3; p>.05) or interaction between Quartile*Groups (all – F3.4,36.1 < 2.5; p>.05).
      Table 2Mean (± standard deviation) of lumbar erector spinae change in ILC after HDTBR (HDT59 value – BDC value) from CTRL (N=8), cAG (N=8), and iAG (N=8). Q1 = medial; Q4 = lateral
      Lumbar erector spinae
      VariableGroupQuartileMixed model repeated measures ANOVA
      Q1Q2Q3Q4QuartilesGroupsQuartiles*Groups
      L1/L2 intervertebral discCRTL1.9±2.92.2±2.61.6±2.10.5±0.6F2.7,46.4=3.28;

      p=.032

      η2partial=0.14
      F2,21=0.49;

      p=.62

      η2partial=0.04
      F4.4,46.4=0.18;

      p=.982

      η2partial=0.02
      cAG2.4±1.62.9±2.91.5±1.91.3±1.1
      iAG1.5±2.42.0±2.11.2±1.30.8±1.2
      L2/L3 intervertebral discCRTL2.1±1.33.2±2.61.8±1.51.4±1.5F1.9,41.2=9.82;

      p<.001

      η2partial=0.32
      F2,21=0.54;

      p=.594

      η2partial=0.05
      F3.9,41.2=0.79;

      p=.582

      η2partial=0.07
      cAG2.3±1.63.2±3.61.8±1.41.6±1.9
      iAG1.3±1.22.2±0.72.0±0.60.5±1.1
      L3/L4 intervertebral discCRTL3.3±3.22.7±1.52.5±1.92.1±1.9F2.0,36.1=1.89;

      p=.140

      η2partial=0.08
      F2,21=0.62;

      p=.548

      η2partial=0.06
      F3.4,36.1=2.54;

      p=.064

      η2partial=0.19
      cAG3.2±3.73.1±2.11.2±1.41.5±1.5
      iAG0.7±1.42.3±1.22.0±0.92.0±1.5
      L4/L5 intervertebral discCRTL4.6±2.92.9±2.32.2±1.52.5±3.5F2.2,45.3=0.80;

      p=.462

      η2partial=0.04
      F2,21=2.84;

      p=.081

      η2partial=0.21
      F4.32,45.3=2.04;

      p=.100

      η2partial=0.16
      cAG2.5±3.73.5±2.92.6±3.53.5±3.2
      iAG-0.6±3.11.5±2.60.8±2.70.1±3.5
      L5/S1 intervertebral discCRTL4.4±3.72.2±2.12.0±1.63.5±2.9F2.1,43.1=1.80;

      p=.156

      η2partial=0.08
      F2,21=1.12;

      p=.344

      η2partial=0.09
      F4.1,43.1=1.79;

      p=.951

      η2partial=0.02
      cAG3.8±4.52.2±3.73.1±4.22.8±4.7
      iAG3.1±3.90.9±3.00.4±2.32.1±4.9
      Fig 4
      Fig. 4(A) ILC of the lumbar erector spinae muscle at BDC and HDT59 (average values of 24 participants) over quartiles (Q1 = medial; Q4 = lateral). (B) Change in ILC after 60 days (HDT59 value – BDC value). * = Significant post hoc test within the Quartile (p<.05).
      At L3/L4, L4/L5 and L5/S1 vertebral levels, the three-way ANOVA showed a main effect of Time (F > 15; p<.001), Quartile (F > 40; p<.001), without any main effect of Groups or other interactions (data and statistical analysis: Supplementary Table 3 and 4).

      Intramuscular lipid concentration in PM and QL muscles after HDTBR

      At all vertebral levels, the three-way ANOVA showed a main effect of Quartile (F > 8; p<.001). No significant main effect of Time, Groups or other interactions were found (data and statistical analysis: Supplementary Table 5, 6, 7, and 8). Analysis of the change in ILC did not differ between Quartiles of PM or QL (Main effects – F2.2,46.7 < 2.5; p>.05). There was no difference between Groups (all – F2,21 < 1; p>.05) or interactions between Quartile*Groups (all – F3.9,41.1 < 2; p>.05) at any intervertebral disc level for these muscles (Table 3 and 4; Fig. 5, Fig. 6).
      Table 3Mean (± standard deviation) of psoas major change in ILC after HDTBR (HDT59 value – BDC value) from CTRL (N=8), cAG (N=8), and iAG (N=8). Q1 = medial; Q4 = lateral
      Psoas Major
      VariableGroupQuartileMixed model repeated measures ANOVA
      Q1Q2Q3Q4QuartilesGroupsQuartiles*Groups
      L1/L2 intervertebral discCRTL-1.5±0.6-0.8±1.9-1.5±1.9-1.0±2.5F1.6,33.6=0.03;

      p=.994

      η2partial=0.01
      F2,21=0.76;

      p=.481

      η2partial=0.07
      F3.2,33.6=0.55;

      p=.663

      η2partial=0.05
      cAG-0.4±4.3-0.3±1.80.6±1.90.8±3.9
      iAG0.3±2.90.2±1.3-0.4±1.1-0.8±1.6
      L2/L3 intervertebral discCRTL-1.0±3.1-0.7±1.10.7±1.2-0.2±1.5F1.7,36.0=0.14;

      p=.838

      η2partial=0.01
      F2,21=1.03;

      p=.373

      η2partial=0.09
      F3.4,36.0=0.25;

      p=.881

      η2partial=0.02
      cAG0.9±2.1-0.7±0.7-0.3±0.9-0.6±2.1
      iAG0.1±2.70.2±1.0-0.0±0.8-0.1±1.4
      L3/L4 intervertebral discCRTL0.1±1.7-0.3±1.5-0.2±1.30.08±1.1F1.9,41.1=0.52;

      p=.670

      η2partial=0.03
      F2,21=0.09;

      p=.911

      η2partial=0.01
      F3.9,41.1=1.17;

      p=.337

      η2partial=0.10
      cAG-0.2±1.60.4±0.8-0.3±0.80.6±0.7
      iAG0.6±2.2-0.2±0.7-0.1±0.9-0.4±1.2
      L4/L5 intervertebral discCRTL0.7±2.3-0.1±1.4-0.1±0.90.0±2.2F1.7,35.8=0.12;

      p=.854

      η2partial=0.01
      F2,21=1.12;

      p=.344

      η2partial=0.09
      F3.8, 35.8=0.41;

      p=.770

      η2partial=0.04
      cAG0.0±2.20.8±1.00.2±1.20.4±1.6
      iAG0.5±2.20.3±1.50.4±0.30.3±1.2
      L5/S1 intervertebral discCRTL-0.2±1.9-0.7±1.3-0.2±1.1-0.3±2.0F1.8,39.4=1.02;

      p=.390

      η2partial=0.05
      F2,21=1.14;

      p=.340

      η2partial=0.09
      F3.8,39.4=0.47;

      p=.745

      η2partial=0.04
      cAG0.8±2.70.8±1.81.0±2.42.0±2.9
      iAG1.3±3.30.8±3.11.0±3.01.2±3.4
      Table 4Mean (± standard deviation) of quadratus lumborum change in ILC after HDTBR (HDT59 value – BDC value) from CTRL (N=8), cAG (N=8), and iAG (N=8). Q1 = medial; Q4 = lateral
      Quadratus lumborum (%)
      VariableGroupQuartileMixed model repeated measures ANOVA
      Q1Q2Q3Q4QuartilesGroupsQuartiles*Groups
      L3/L4 intervertebral discCRTL-0.9±2.10.2±2.80.5±2.9-0.4±4.0F1.7,37.3=1.85;

      p=.148

      η2partial=0.08
      F2,21=0.36;

      p=.698

      η2partial=0.03
      F3.6,37.3=1.02;

      p=.403

      η2partial=0.09
      cAG0.6±2.60.1±1.0-0.7±2.9-2.3±2.9
      iAG0.5±3.40.7±1.6-0.1±1.1-0.5±3.1
      L4/L5 intervertebral discCRTL-0.9±3.30.0±2.90.6±1.80.0±3.0F2.1,43.5=0.67;

      p=.518

      η2partial=0.03
      F2,21=0.18;

      p=.836

      η2partial=0.02
      F4.1,43.5=2.27;

      p=.075

      η2partial=0.18
      cAG0.6±2.1-0.3±1.3-0.4±1.2-1.2±2.1
      iAG1.2±1.4-0.5±1.3-0.4±1.00.0±3.0
      L5/S1 intervertebral discCRTL-1.2±2.0-0.3±0.9-0.2±1.1-1.2±1.9F2.2,46.7=2.29;

      p=.107

      η2partial=0.10
      F2,21=0.64;

      p=.539

      η2partial=0.06
      F4.4,46.7=0.64;

      p=.655

      η2partial=0.06
      cAG-0.4±1.6-0.0±0.70.2±1.2-0.1±1.4
      iAG-0.3±1.8-0.1±0.90.1±0.8-1.6±3.2
      Fig 5
      Fig. 5(A) ILC of the psoas major muscle at BDC and HDT59 (average values of 24 participants) over quartiles (Q1 = medial; Q4 = lateral). (B) Change in ILC after 60 days (HDT59 value – BDC value).
      Fig 6
      Fig. 6(A) Intramuscular lipid concentration of the quadratus lumborum at BDC and HDT59 (average values of 24 participants) over quartiles (Q1 = medial; Q4 = lateral). (B) Change in ILC after 60 days (HDT59 value – BDC value).

      Discussion

      This study identified and quantified the spatial pattern of accumulation of ILC in the lumbar musculature after 60 day HDTBR. Our results demonstrated a greater increase in lipid proportion in medial and lateral regions of the LM muscle at the L4/L5 and S1/L5 vertebral levels than in other regions and levels of the muscle. In contrast, a smaller increase in lipid proportions was observed in the lateral region of the LES muscle at the L1/L2 and L2/L3 vertebral levels than in other regions and levels of the muscle. No changes in lipid proportion were detected in the PM and QL muscles after the HDTBR. Finally, exposure to AG did not change the accumulation of ILC.

      Intramuscular lipid concentration in LM muscle

      The current results showed that ILC increased in all LM muscle regions at all vertebral levels at the end of HDTBR compared with before. However, the lateral and medial regions and lower vertebral levels were most affected: At the L4/L5 vertebral level, more than 4% lipid accumulation was found in the centro-medial (Q1 and Q2) lateral (Q4) regions. At the L5/S1 vertebral level, ILC had increased by more than 5% in the lateral (Q4) and medial region (Q1) after HDTBR. This inhomogeneous accumulation of ILC may be related to the complex architectural structure and functional differentiation between regions of the LM muscle [
      • Macintosh JE
      • Bogduk N
      The biomechanics of the lumbar multifidus.
      ,
      • Moseley GL
      • Hons BP
      • Hodges PW
      • Gandevia SC
      Deep and Superficial Fibers of the Lumbar Multifidus Muscle Are Differentially Active During Voluntary Arm Movements.
      ,
      • Moseley GL
      • Hodges PW
      • Gandevia SC
      External perturbation of the trunk in standing humans differentially activates components of the medial back muscles.
      ]. The LM fascicles originate from the spinous process and adjacent lamina of each lumbar vertebra, descend caudolaterally, over 2-4 vertebral levels, and attach to mamillary processes, the iliac crest, and the dorsal surface of the sacrum [
      • Macintosh JE
      • Valencia F
      • Bogduk N
      • Munro R
      The morphology multifidus of the human lumbar.
      ,
      • Kalimo H
      • Rantanen J
      • Viljanen T
      • Einola S
      Lumbar muscles: Structure and function.
      ]. In the current study, the lateral region (Q4) at the L4/L5 and L5/S1 vertebral levels is likely to represent the distal portion of the longer fascicles originating from upper lumbar vertebrae (laminae of L1 and L2) and attaching to the L4 and L5 mamillary processes and the dorsal sacral surface [
      • Macintosh JE
      • Valencia F
      • Bogduk N
      • Munro R
      The morphology multifidus of the human lumbar.
      ]. These longer fascicles have a moment arm that can increase the lumbar lordosis and extend the lumbar spine, as their line of action falls behind the lumbar lordotic curve [
      • Bogduk N
      • Macintosh JE
      • Pearcy MJ
      A universal model of the lumbar back muscles in the upright position.
      ], and this is supported by electromyography (EMG) data [
      • Moseley GL
      • Hons BP
      • Hodges PW
      • Gandevia SC
      Deep and Superficial Fibers of the Lumbar Multifidus Muscle Are Differentially Active During Voluntary Arm Movements.
      ,
      • Moseley GL
      • Hodges PW
      • Gandevia SC
      External perturbation of the trunk in standing humans differentially activates components of the medial back muscles.
      ]. As other EMG recordings have shown that the paraspinal extensor muscles are not active in reclined/lying positions [
      • Andersson GBJ
      • Ortengren R
      • Nachemson A
      Intradiskal pressure, intra-abdominal pressure and myoelectric back muscle activity related to posture and loading.
      ], this may explain the ILC increase seen in the lateral, longer torque producing fibers in response to HDTBR.
      The deeper portion of the centro-medial regions (Q1-Q2) at L4/L5 vertebral level would represent some of the shorter fascicles originating from the laminae of L3 and L4 and attaching to the mamillary processes of L5, the capsule of the zygapophyseal joints, and the dorsal surface of the sacrum [
      • Macintosh JE
      • Valencia F
      • Bogduk N
      • Munro R
      The morphology multifidus of the human lumbar.
      ]. The deeper, shorter fascicles are less able to generate torque than the more superficial, longer fascicles. Some studies suggest that the shorter fascicles have a higher percentage of Type I muscle fibers [
      • Sirca A
      • Kostevc V
      The fibre type composition of thoracic and lumbar paravertebral muscles in man.
      ], indicating a propensity for low-level tonic activation. EMG investigations have demonstrated that the short deep fibers are less affected by the direction of applied force than the long fascicles in functional loading tasks [
      • Moseley GL
      • Hons BP
      • Hodges PW
      • Gandevia SC
      Deep and Superficial Fibers of the Lumbar Multifidus Muscle Are Differentially Active During Voluntary Arm Movements.
      ]. Although the mechanisms are unclear, our results may indicate that this deep fiber region is differently impacted by unloading than the adjacent longer fascicles. However, it is also important to consider that ILC accumulates adjacent to the bone (spinous process/lamina) and adjacent to the thick inter-muscular fascia (lateral aspect of muscle), as shown by our measurements before HDTBR, suggesting that these portions of the muscles may have greater potential to accumulate lipid.
      The level-specific accumulation of fat might be explained by regional differentiation in loading. The lower vertebral levels have close proximity to the sacropelvic complex [
      • Macintosh JE
      • Bogduk N
      The biomechanics of the lumbar multifidus.
      ,
      • Aspden RM
      Intra-abdominal pressure and its role in spinal mechanics.
      ,
      • Bartelink DL
      The role of abdominal pressure in relieving the pressure on the lumbar intervertebral discs.
      ] and are exposed to greater mechanical stress than the upper lumbar levels [
      • Shirazi-Adl A
      • Parnianpour M
      Effect of changes in lordosis on mechanics of the lumbar spine-lumbar curvature in lifting.
      ]. Animal studies have demonstrated that high mechanical stress regulates the expression of transcriptional factors influencing cell differentiation [
      • Akimoto T
      • Ushida T
      • Miyaki S
      • Tateishi T
      • Fukubayashi T
      Mechanical stretch is a down-regulatory signal for differentation of C2C12 myogenic cells.
      ]. Decreased cyclic mechanical stress leads to upregulation of adipogenic transcription factors [
      • Kook SH
      • Lee HJ
      • Chung WT
      • Hwang IH
      • Lee SA
      • Kim BS
      • et al.
      Cyclic mechanical stretch stimulates the proliferation of C2C12 myoblasts and inhibits their differentiation via prolonged activation of p38 MAPK.
      ] and reduced expression of factors that inhibit myoblast transdifferentiation to adipocytes [
      • Akimoto T
      • Ushida T
      • Miyaki S
      • Akaogi H
      • Tsuchiya K
      • Yan Z
      • et al.
      Mechanical stretch inhibits myoblast-to-adipocyte differentiation through Wnt signaling.
      ]. Taken together with the results of the current study, we hypothesize that the muscle fibers at the lower vertebral levels, which are subjected to higher mechanical stress [
      • Shirazi-Adl A
      • Parnianpour M
      Effect of changes in lordosis on mechanics of the lumbar spine-lumbar curvature in lifting.
      ], might have a higher propensity to accumulate lipids when the lumbar spine is unloaded.

      Intramuscular lipid concentration in LES muscle

      Similar to the LM muscle, ILC in the LES muscle also tended to accumulate at the lowest lumbar vertebral levels, and a high percentage change in ILC with HDTBR was observed in medial and lateral regions at the L4/L5 and L5/S1 vertebral levels. Although still unknown, structural muscle changes close to these regions may have important biomechanical implications, as internal forces and moments produced by the muscles across those joints can be altered.
      After HDTBR, lower intramuscular lipid accumulation was also found in the lateral region of LES muscle at the L1/L2 and L2/L3 vertebral levels (∼ 1%) than in the centro-medial regions. The LES muscle consists of two large muscles: longissimus lumborum, medially, and iliocostalis lumborum, laterally [
      • Macintosh JE
      • Bogduk N
      The morphology of the lumbar erector spinae.
      ,
      • Bogduk N
      A reappraisal of the anatomy of the human lumbar erector spinae.
      ]. At the L1 and L2 vertebral levels, the lateral division of the LES muscle contains many fibres that originate from the lumbar aponeurosis and insert on the 12-7th ribs (thoracic parts of iliocostalis lumborum muscle) [
      • Kalimo H
      • Rantanen J
      • Viljanen T
      • Einola S
      Lumbar muscles: Structure and function.
      ]. Other functions of this lateral region of LES, such as a role in respiration, might lead to a lesser effect of HDTBR on this region. Although quiet breathing in supine does not involve activation of the abdominal and paraspinal muscles [
      • Campbell EJM
      • Green JH
      The variations in intra-abdominal pressure and the activity of the abdominal muscles during breathing; a study in man.
      ,
      • Shirley D
      • Hodges PW
      • Eriksson AEM
      • Gandevia SC
      Spinal stiffness changes throughout the respiratory cycle.
      ], inspiratory demand is increased in HDTBR as the diaphragm has to overcome the weight of the abdominal organs [
      • Montmerle S
      • Spaak J
      • Linnarsson D
      Lung function during and after prolonged head-down bed rest.
      ], and this has been suggested to explain the absence of major deconditioning of respiratory muscles after HDTBR [
      • Montmerle S
      • Spaak J
      • Linnarsson D
      Lung function during and after prolonged head-down bed rest.
      ]. Other functions such as a role in thoracolumbar rotation to move around the bed might also be involved.

      Intramuscular lipid concentration in PM and QL muscles

      In contrast to the LM and LES muscles, ILC in the PM and QL muscles did not accumulate at any vertebral levels and did not show any inhomogeneity of the spatial distribution of accumulation of ILC after HDTBR. Although the absence of ILC changes in the PM muscle is consistent with previous studies in astronauts [
      • Burkhart K
      • Allaire B
      • Bouxsein ML
      Negative Effects of Long-duration Spaceflight on Paraspinal Muscle Morphology.
      ,
      • McNamara KP
      • Greene KA
      • Moore AM
      • Lenchik L
      • Weaver AA
      Lumbopelvic muscle changes following long-duration spaceflight.
      ], our results contrast with previous results showing an increase in the QL muscle after a 6 month spaceflight [
      • Burkhart K
      • Allaire B
      • Bouxsein ML
      Negative Effects of Long-duration Spaceflight on Paraspinal Muscle Morphology.
      ,
      • McNamara KP
      • Greene KA
      • Moore AM
      • Lenchik L
      • Weaver AA
      Lumbopelvic muscle changes following long-duration spaceflight.
      ]. The absence of ILC accumulation in the current investigation may be explained by some essential distinctions between HDTBR and spaceflight. Microgravity is accompanied by the complete unloading of postural muscles, whereas these muscles still need to work against gravity during HDTBR, particularly when changing postures. For instance, our participants are likely to have recruited the QL muscle to laterally flex their spine during hygiene activities and when rolling onto their sides, dressing, and eating. Daily side-lying activities during the current HDTBR may have been sufficient to prevent the accumulation of ILC in the QL muscle. Alternatively, ILC accumulation in QL muscle may require long-lasting body unloading.

      Clinical implications for patients and operational relevance for spaceflight

      The present results expand previous findings indicating that higher adipose tissue accumulation in low back muscles is associated with physical inactivity [
      • Sions JM
      • Coyle PC
      • Velasco TO
      • Elliott JM
      • Hicks GE
      Multifidi Muscle Characteristics and Physical Function Among Older Adults With and Without Chronic Low Back Pain.
      ]. Greater amounts of ILC in the LM muscle, rather than changes in muscle size, have been associated with an increased risk of high-intensity pain among adults with LBP [
      • Teichtahl AJ
      • Urquhart DM
      • Wang Y
      • Wluka AE
      • Wijethilake P
      • O'Sullivan R
      • et al.
      Fat infiltration of paraspinal muscles is associated with low back pain, disability, and structural abnormalities in community-based adults.
      ]. Future studies should investigate whether reconditioning interventions following a period of prolonged unloading, such as HDTBR or spaceflight, can reverse the localized accumulation of ILC in the LM and LES muscles.
      Space Agencies have renewed their focus on the Moon missions [
      • Landgraf M
      • et al.
      International Space Exploration Coordination Group
      Lunar surface concept of operations for the global exploration roadmap lunar surface exploration scenario. IAC-21,A5,1,5,x66702.
      ], and more astronauts are expected to go into space and work on the Moon in this decade. However, prolonged spaceflight may increase the risk of spinal injury, and optimal countermeasures to mitigate the deconditioning effect of mechanical unloading to the lumbar spine may be crucial. Our results showed that daily AG protocols did not significantly protect the muscles from these effects. Future studies should consider whether other paradigms are effective, such as prolonged AG exposure, higher compressive Gz, or the combination of AG with specific lumbar muscle exercises.

      Limitations

      The results of the current study should be interpreted in consideration of its limitations. Small sample sizes are a common limitation of HDTBR studies due to the intrinsically complex nature and expense of these studies. Nevertheless, similar HDT bed rest studies have included similar numbers of groups and participants and have found differences using the same statistical approach as the current study [
      • Belavý DL
      • Armbrecht G
      • Gast U
      • Richardson CA
      • Hides JA
      • Felsenberg D
      Countermeasures against lumbar spine deconditioning in prolonged bed rest : resistive exercise with and without whole body vibration.
      ,
      • Holguin N
      • Muir J
      • Rubin C
      • Judex S
      Short applications of very low-magnitude vibrations attenuate expansion of the intervertebral disc during extended bed rest.
      ,
      • De Martino E
      • Salomoni SE
      • Hodges PW
      • Hides J
      • Lindsay K
      • Debuse D
      • et al.
      Intermittent short-arm centrifugation is a partially effective countermeasure against upright balance deterioration following 60-day head-down tilt bed rest.
      ]. Because of the small sample size in the groups, only large effect sizes from countermeasures can be observed, and more subtle effects are likely to go unnoticed.
      Although quartile segmentation has frequently been used to study ILC in cervical and lumbar muscles, this segmentation may be insufficient to detect more localized changes, and, due to the complex architecture of the muscles studied, future studies could consider a higher spatial resolution.

      Conclusion

      The current investigation suggests that ILC increased more in the medial and lateral regions than central regions of the lower LM muscle after 60 day HDTBR. In contrast, the LES muscle accumulated less ILC in the lateral than medio-central region of the upper lumbar spine. The present findings may represent a new target for lumbar muscle reconditioning for those exposed to prolonged extreme physical inactivity, astronauts, elderly, or individuals with chronic LBP.

      Declaration of competing interests

      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.

      Acknowledgments

      The AGBRESA study was funded by the German Aerospace Center, the European Space Agency (contract number: 4000113871/15/NL/PG), and the National Aeronautics and Space Administration (contract number: 80JSC018P0078). The study was performed at the “:envihab” research facility of the DLR Institute of Aerospace Medicine. Funding for this ESA-selected project (ESA-HSO-U-LE-0629) was received from the STFC/UK Space Agency (ST/R005753/1).

      Appendix. Supplementary materials

      References

        • Bailey JF
        • Miller SL
        • Khieu K
        • Neill CWO
        • Healey RM
        • Coughlin DG
        • et al.
        From the international space station to the clinic : how prolonged unloading may disrupt lumbar spine stability.
        Spine J. 2018; 18: 7-14https://doi.org/10.1016/j.spinee.2017.08.261
        • Chang DG
        • Healey RM
        • Snyder AJ
        • Sayson JV
        • Macias BR
        • Coughlin DG
        • et al.
        Lumbar spine paraspinal muscle and intervertebral disc height changes in astronauts after long-duration spaceflight on the International Space Station.
        Spine (Phila Pa 1976). 2016; 41: 1917-1924https://doi.org/10.1097/BRS.0000000000001873
        • Belavý DL
        • Adams M
        • Brisby H
        • Cagnie B
        • Danneels L
        • Fairbank J
        • et al.
        Disc herniations in astronauts: What causes them, and what does it tell us about herniation on earth?.
        Eur Spine J. 2016; 25: 144-154https://doi.org/10.1007/s00586-015-3917-y
        • Hides JA
        • Lambrecht G
        • Sexton CT
        • Pruett C
        • Petersen N
        • Jaekel P
        • et al.
        The effects of exposure to microgravity and reconditioning of the lumbar multifidus and anterolateral abdominal muscles: implications for people with LBP.
        Spine J. 2020; : 1-14https://doi.org/10.1016/j.spinee.2020.09.006
        • Belavý DL
        • Armbrecht G
        • Richardson CA
        • Felsenberg D
        • Hides JA
        Muscle atrophy and changes in spinal morphology: Is the lumbar spine vulnerable after prolonged bed-rest?.
        Spine (Phila Pa 1976). 2011; 36: 137-145https://doi.org/10.1097/BRS.0b013e3181cc93e8
        • Burkhart K
        • Allaire B
        • Bouxsein ML
        Negative Effects of Long-duration Spaceflight on Paraspinal Muscle Morphology.
        Spine (Phila Pa 1976). 2019; 44: 879-886https://doi.org/10.1097/BRS.0000000000002959
        • Hodges PW
        • Danneels L
        Changes in Structure and Function of the Back Muscles in Low Back Pain: Different Time Points, Observations, and Mechanisms.
        J Orthop Sport Phys Ther. 2019; 49: 464-476https://doi.org/10.2519/jospt.2019.8827
        • Kjaer P
        • Bendix T
        • Sorensen JS
        • Korsholm L
        • Leboeuf-Yde C
        Are MRI-defined fat infiltrations in the multifidus muscles associated with low back pain?.
        BMC Med. 2007; 5: 1-10https://doi.org/10.1186/1741-7015-5-2
        • Hildebrandt M
        • Fankhauser G
        • Meichtry A
        • Luomajoki H
        Correlation between lumbar dysfunction and fat infiltration in lumbar multifidus muscles in patients with low back pain.
        BMC Musculoskelet Disord. 2017; : 1-9https://doi.org/10.1186/s12891-016-1376-1
        • Pezolato A
        • De Vasconcelos EE
        • Defino HLA
        • Defino A
        • Nogueira-Barbosa MH
        Fat infiltration in the lumbar multifidus and erector spinae muscles in subjects with sway-back posture.
        Eur Spine J. 2012; 21: 2158-2164https://doi.org/10.1007/s00586-012-2286-z
        • Dahlqvist JR
        • Vissing CR
        • Hedermann G
        • Thomsen C
        • Vissing J
        Fat Replacement of Paraspinal Muscles with Aging in Healthy Adults.
        Med Sci Sports Exerc. 2017; 49: 595-601https://doi.org/10.1249/MSS.0000000000001119
        • Crawford RJ
        • Filli L
        • Elliott JM
        • Nanz D
        • Fischer MA
        • Marcon M
        • et al.
        Age- and level-dependence of fatty infiltration in lumbar paravertebral muscles of healthy volunteers.
        Am J Neuroradiol. 2015; 37: 742-748https://doi.org/10.3174/ajnr.A4596
        • Crawford RJ
        • Volken ÃT
        • Elliott J
        • Hoggarth MA
        • Dpt MS
        • Samartzis D
        Geography of Lumbar Paravertebral Muscle.
        Spine (Phila Pa 1976). 2019; 44: 1294-1302https://doi.org/10.1097/BRS.000000 0000003060
        • McNamara KP
        • Greene KA
        • Moore AM
        • Lenchik L
        • Weaver AA
        Lumbopelvic muscle changes following long-duration spaceflight.
        Front Physiol. 2019; 10https://doi.org/10.3389/fphys.2019.00627
        • Hodges PW
        • James G
        • Blomster L
        • Hall L
        • Schmid A
        • Shu C
        • et al.
        Multifidus Muscle Changes after Back Injury Are Characterized by Structural Remodeling of Muscle, Adipose and Connective Tissue, but Not Muscle Atrophy: Molecular and Morphological Evidence.
        Spine (Phila Pa 1976). 2015; 40: 1057-1071https://doi.org/10.1097/BRS.0000000000000972
        • Hodges PW
        • James G
        • Blomster L
        • Hall L
        • Schmid AB
        • Shu C
        • et al.
        Can proinflammatory cytokine gene expression explain multifidus muscle fiber changes after an intervertebral disc lesion?.
        Spine (Phila Pa 1976). 2014; 39: 1010-1017https://doi.org/10.1097/BRS.0000000000000318
        • Teichtahl AJ
        • Urquhart DM
        • Wang Y
        • Wluka AE
        • Wijethilake P
        • O'Sullivan R
        • et al.
        Fat infiltration of paraspinal muscles is associated with low back pain, disability, and structural abnormalities in community-based adults.
        Spine J. 2015; 15: 1593-1601https://doi.org/10.1016/j.spinee.2015.03.039
        • Fortin M
        • Lazáry À
        • Varga PP
        • Battié MC
        Association between paraspinal muscle morphology, clinical symptoms and functional status in patients with lumbar spinal stenosis.
        Eur Spine J. 2017; 26: 2543-2551https://doi.org/10.1007/s00586-017-5228-y
        • Lee SH
        • Park SW
        • Kim YB
        • Nam TK
        • Lee YS
        The fatty degeneration of lumbar paraspinal muscles on computed tomography scan according to age and disc level.
        Spine J. 2017; 17: 81-87https://doi.org/10.1016/j.spinee.2016.08.001
        • Smith AC
        • Albin SR
        • Abbott R
        • Crawford RJ
        • Hoggarth MA
        • Wasielewski M
        • et al.
        Confirming the geography of fatty infiltration in the deep cervical extensor muscles in whiplash recovery.
        Sci Rep. 2020; 10: 1-8https://doi.org/10.1038/s41598-020-68452-x
        • Abbott R
        • Peolsson A
        • West J
        • Elliott JM
        • Åslund U
        • Karlsson A
        • et al.
        The qualitative grading of muscle fat infiltration in whiplash using fat and water magnetic resonance imaging.
        Spine J. 2018; 18: 717-725https://doi.org/10.1016/j.spinee.2017.08.233
        • Macintosh JE
        • Valencia F
        • Bogduk N
        • Munro R
        The morphology multifidus of the human lumbar.
        Clin Biomech. 1986; 1: 196-204
        • Macintosh JE
        • Bogduk N
        The morphology of the lumbar erector spinae.
        Spine (Phila Pa 1976). 1987; 12: 658-668
        • Hargens AR
        • Vico L
        Long-duration bed rest as an analog to microgravity.
        J Appl Physiol. 2016; 120: 891-903https://doi.org/10.1152/japplphysiol.00935.2015
        • Pavy-Le Traon A
        • Heer M
        • Narici MV
        • Rittweger J
        • Vernikos J
        From space to Earth: Advances in human physiology from 20 years of bed rest studies.
        Eur J Appl Physiol. 2007; 101: 143-194https://doi.org/10.1007/s00421-007-0474-z
        • Clément GR
        • Bukley AP
        • Paloski WH
        Artificial gravity as a countermeasure for mitigating physiological deconditioning during long-duration space missions.
        Front Syst Neurosci. 2015; 9: 1-11https://doi.org/10.3389/fnsys.2015.00092
        • Kotovskaya AR
        The problem of artificial gravity: The current state and prospects.
        Hum Physiol. 2010; 36: 780-787https://doi.org/10.1134/s0362119710070078
        • De Martino E
        • Hides J
        • Elliott JM
        • Hoggarth M
        • Zange J
        • Lindsay K
        • et al.
        Lumbar muscle atrophy and increased relative intramuscular lipid concentration are not mitigated by daily artificial gravity following 60-day head-down tilt bed rest.
        J Appl Physiol. 2021; https://doi.org/10.1152/japplphysiol.00990.2020
        • Rittweger J
        • Bareille MP
        • Clément G
        • Linnarsson D
        • Paloski WH
        • Wuyts F
        • et al.
        Short-arm centrifugation as a partially effective musculoskeletal countermeasure during 5-day head-down tilt bed rest—results from the BRAG1 study.
        Eur J Appl Physiol. 2015; 115: 1233-1244https://doi.org/10.1007/s00421-015-3120-1
        • Linnarsson D
        • Hughson RL
        • Fraser KS
        • Clément G
        • Karlsson LL
        • Mulder E
        • et al.
        Effects of an artificial gravity countermeasure on orthostatic tolerance, blood volumes and aerobic power after short-term bed rest (BR-AG1).
        J Appl Physiol. 2015; 118: 29-35https://doi.org/10.1152/japplphysiol.00061.2014
        • Sundblad P
        • Orlov O
        • Angerer O
        • Larina I
        • Cromwell R
        Guidelines for standardization of bed rest studies in the spaceflight context.
        International Academy of Astronautics (IAA). 2014; 1 (Internatio): 1-70
        • Ganse B
        • Bosutti A
        • Drey M
        • Degens H
        Sixty days of head-down tilt bed rest with or without artificial gravity do not affect the neuromuscular secretome.
        Exp Cell Res. 2021; 399: 112463https://doi.org/10.1016/j.yexcr.2020.112463
        • Frett T
        • Green DA
        • Mulder E
        • Noppe A
        • Arz M
        • Pustowalow W
        • et al.
        Tolerability of daily intermittent or continuous short-arm centrifugation during 60-day 6o head down bed rest (AGBRESA study).
        PLoS One. 2020; 15: 1-11https://doi.org/10.1371/journal.pone.0239228
        • Kramer A
        • Venegas-Carro M
        • Mulder E
        • Lee JK
        • Moreno-Villanueva M
        • Bürkle A
        • et al.
        Cardiorespiratory and Neuromuscular Demand of Daily Centrifugation: Results From the 60-Day AGBRESA Bed Rest Study.
        Front Physiol. 2020; 11: 1-9https://doi.org/10.3389/fphys.2020.562377
        • Mhuiris ÁN
        • Volken T
        • Elliott JM
        • Hoggarth M
        • Samartzis D
        • Crawford RJ
        Reliability of quantifying the spatial distribution of fatty infiltration in lumbar paravertebral muscles using a new segmentation method for T1-weighted MRI.
        BMC Musculoskelet Disord. 2016; 17: 1-7https://doi.org/10.1186/s12891-016-1090-z
        • Smith AC
        • Parrish TB
        • Abbott R
        • Hoggarth MA
        • Mendoza K
        • Chen YF
        • et al.
        Muscle-fat MRI: 1.5 Tesla and 3.0 Tesla versus histology.
        Muscle and Nerve. 2014; 50: 170-176https://doi.org/10.1002/mus.24255
        • Macintosh JE
        • Bogduk N
        The biomechanics of the lumbar multifidus.
        Clin Biomech. 1986; 1: 205-213https://doi.org/10.1016/0268-0033(86)90147-6
        • Moseley GL
        • Hons BP
        • Hodges PW
        • Gandevia SC
        Deep and Superficial Fibers of the Lumbar Multifidus Muscle Are Differentially Active During Voluntary Arm Movements.
        Spine (Phila Pa 1976). 2002; 27: 29-36
        • Moseley GL
        • Hodges PW
        • Gandevia SC
        External perturbation of the trunk in standing humans differentially activates components of the medial back muscles.
        J Physiol. 2003; 547: 581-587https://doi.org/10.1113/jphysiol.2002.024950
        • Kalimo H
        • Rantanen J
        • Viljanen T
        • Einola S
        Lumbar muscles: Structure and function.
        Ann Med. 1989; 21: 353-359https://doi.org/10.3109/07853898909149220
        • Bogduk N
        • Macintosh JE
        • Pearcy MJ
        A universal model of the lumbar back muscles in the upright position.
        Spine (Phila Pa 1976). 1992; 17: 897-913
        • Andersson GBJ
        • Ortengren R
        • Nachemson A
        Intradiskal pressure, intra-abdominal pressure and myoelectric back muscle activity related to posture and loading.
        Clin Orthop Relat Res. 1977; 129 (No): 156-164https://doi.org/10.1097/00003086-197711000-00018
        • Sirca A
        • Kostevc V
        The fibre type composition of thoracic and lumbar paravertebral muscles in man.
        J Anat. 1985; 141: 131-137
        • Aspden RM
        Intra-abdominal pressure and its role in spinal mechanics.
        Clin Biomech. 1987; 2: 168-174https://doi.org/10.1016/0268-0033(87)90010-6
        • Bartelink DL
        The role of abdominal pressure in relieving the pressure on the lumbar intervertebral discs.
        J Bone Joint Surg Br. 1957; 39: 718-725https://doi.org/10.1302/0301-620x.39b4.718
        • Shirazi-Adl A
        • Parnianpour M
        Effect of changes in lordosis on mechanics of the lumbar spine-lumbar curvature in lifting.
        J Spinal Disord. 1999; 12: 426-447
        • Akimoto T
        • Ushida T
        • Miyaki S
        • Tateishi T
        • Fukubayashi T
        Mechanical stretch is a down-regulatory signal for differentation of C2C12 myogenic cells.
        Mater Sci Eng C. 2001; 17: 75-78https://doi.org/10.1016/S0928-4931(01)00340-X
        • Kook SH
        • Lee HJ
        • Chung WT
        • Hwang IH
        • Lee SA
        • Kim BS
        • et al.
        Cyclic mechanical stretch stimulates the proliferation of C2C12 myoblasts and inhibits their differentiation via prolonged activation of p38 MAPK.
        Mol Cells. 2008; 25: 479-486
        • Akimoto T
        • Ushida T
        • Miyaki S
        • Akaogi H
        • Tsuchiya K
        • Yan Z
        • et al.
        Mechanical stretch inhibits myoblast-to-adipocyte differentiation through Wnt signaling.
        Biochem Biophys Res Commun. 2005; 329: 381-385https://doi.org/10.1016/j.bbrc.2005.01.136
        • Bogduk N
        A reappraisal of the anatomy of the human lumbar erector spinae.
        J Anat. 1980; 131: 525-540
        • Campbell EJM
        • Green JH
        The variations in intra-abdominal pressure and the activity of the abdominal muscles during breathing; a study in man.
        J Physiol. 1953; 122: 282-290
        • Shirley D
        • Hodges PW
        • Eriksson AEM
        • Gandevia SC
        Spinal stiffness changes throughout the respiratory cycle.
        J Appl Physiol. 2003; 95: 1467-1475https://doi.org/10.1152/japplphysiol.00939.2002
        • Montmerle S
        • Spaak J
        • Linnarsson D
        Lung function during and after prolonged head-down bed rest.
        J Appl Physiol. 2002; 92: 75-83https://doi.org/10.1152/jappl.2002.92.1.75
        • Sions JM
        • Coyle PC
        • Velasco TO
        • Elliott JM
        • Hicks GE
        Multifidi Muscle Characteristics and Physical Function Among Older Adults With and Without Chronic Low Back Pain.
        Arch Phys Med Rehabil. 2017; 98: 51-57https://doi.org/10.1016/j.apmr.2016.07.027
        • Landgraf M
        • et al.
        • International Space Exploration Coordination Group
        Lunar surface concept of operations for the global exploration roadmap lunar surface exploration scenario. IAC-21,A5,1,5,x66702.
        IAC-21,A5,1,5,x66702. IAA symposium on human exploration of the solar system, 2021
        • Belavý DL
        • Armbrecht G
        • Gast U
        • Richardson CA
        • Hides JA
        • Felsenberg D
        Countermeasures against lumbar spine deconditioning in prolonged bed rest : resistive exercise with and without whole body vibration.
        J Appl Physiol. 2010; 109: 1801-1811https://doi.org/10.1152/japplphysiol.00707.2010
        • Holguin N
        • Muir J
        • Rubin C
        • Judex S
        Short applications of very low-magnitude vibrations attenuate expansion of the intervertebral disc during extended bed rest.
        Spine J. 2009; 9: 470-477https://doi.org/10.1016/j.spinee.2009.02.009
        • De Martino E
        • Salomoni SE
        • Hodges PW
        • Hides J
        • Lindsay K
        • Debuse D
        • et al.
        Intermittent short-arm centrifugation is a partially effective countermeasure against upright balance deterioration following 60-day head-down tilt bed rest.
        J Appl Physiol. 2021; https://doi.org/10.1152/japplphysiol.00180.2021