The relationship between acceleration during running at optimal speed and changes in the T2 times of the lumbar intervertebral disc
===================================================================================================================================

* Takayoshi Hakkaku
* Yoshiaki Kubo
* Koji Koyama
* Koichi Nakazato
* Takashi Okada
* Kenji Hiranuma

## Abstract

T2 times by magnetic resonance imaging techniques has been used to examine the beneficial effects of exercise on the human intervertebral disc (IVDs). The magnitude, frequency, and duration of running can contribute to both beneficial and detrimental responses. Slow running (7-9 km/h) is considered the optimal range. By revealing the detailed loading direction and acute changes in T2 times that occur in the IVDs, the positive effects of running on the IVDs can be explored. This study aimed to investigate the relationship between load during slow running and changes in the T2 times of lumbar IVDs before and after running. Sixteen healthy male students were fitted with a triaxial accelerator and ran on a treadmill at 8 km/h for 1 min. Three lumbar T2 times from the L3/L4 to L5/S1 levels were measured before and after exercise using magnetic resonance imaging, and the analysis divided into five regions of interest. Acceleration was 0.23 ± 0.06 root mean square in the X-axis (mediolateral), 1.37 ± 0.08 in the Y-axis (vertical), and 0.30 ± 0.06 in the Z-axis (anteroposterior). An analysis of the correlation between the change in T2 time and acceleration showed a strong correlation, particularly in the Z-axis (anteroposterior direction) acceleration. At L3/L4, a positive correlation was observed for the posterior nucleus (r = 0.72, p = 0.002, *R2* = 0.59), at L4/L5, a positive correlation was observed for the central nucleus (r = 0.73, p = 0.003, R2 = 0.49); in L5/S1, a negative correlation was observed for the anterior annulus fibrosus (r = -0.73, p = 0.01, *R2* = 0.48). These results suggest that anteroposterior loading may play a significant role in the response of the IVDs.

## Introduction

Human intervertebral discs (IVDs) are primarily avascular tissues, and the maintenance of homeostasis in disc cells depends on nutrient exchange between the disc and the adjacent vertebrae. Nutrients such as glucose and oxygen are transported into the vertebral column by two mechanisms: diffusion of nutrients from the vertebral vascular supply terminating at the vertebral endplate and fluid convection due to loading [1–3]. Oxygen plays an important role in disc cell metabolism, and the diminished transport of these nutrients is often proposed as a factor in age-related disc degeneration (IDD) [4].

Quantitative magnetic resonance imaging (MRI) techniques using T2 values can assess IVDs, and T2 times are strongly correlated with its biochemical composition [5]. T2 times are based on two parameters: water content and collagen arrangement in the IVDs, and as the IDD, the water content decreases and the collagen arrangement becomes irregular, resulting in a decrease in the T2 times of the nucleus pulposus [6, 7]. Investigations using T2 time to elucidate the mechanism of IDD suggest that in the nucleus pulposus, which undergoes fewer degenerative changes, T2 time is temporarily decreased by exercise stress and recovered by rest. However, progressive degeneration does not cause changes in T2, suggesting irreversible progression [8].

Recently, T2 has been used to examine the beneficial effects of exercise on the IVDs. Belavy et al. [9] reported that athletic individuals had longer T2 times of lumbar IVDs, whereas long-distance runners exhibited longer T2 times of lumbar IVDs, suggesting greater vertebral height and IVDs hypertrophy. Another study found longer T2 times of the IVDs annulus fibrosus and better IDD status in men with more than 10 years of running experience, suggesting that long-term running may delay IDD [10]. In contrast, prolonged dynamic loading during running has been reported to cause signs of degeneration. Takatalo et al. [11] investigated the incidence of IDD in adolescents and reported that running at least twice a week among endurance runners was associated with IDD. This result is inconsistent with those of earlier studies that proposed that running benefits IVDs.

Thus, even during running, the magnitude, frequency, and duration of force may contribute to both beneficial and deleterious responses. Belavy et al. [9] used triaxial accelerometer to determine the mean absolute deviation (MAD) during running and found that 0.46 to 0.78 G is the optimal range for lumbar IVDs. In other words, by clarifying the detailed loading direction and T2 times of the acute changes that occur in IVDs during slow running at 7–9 km/h (the optimal range), the positive effects of running on IVDs can be explored. This also clarifies the inconsistency in the impact of running on the IVDs. Thus, the relationship between the direction of loading and the change in T2 times of the IVDs due to exercise loading in the optimal range of MAD has yet to be elucidated. Therefore, this study aimed to investigate the relationship between load during slow running and changes in the T2 times of lumbar IVDs before and after running. We hypothesized that investigating the load applied by the MAD in the optimal range would provide clues for improving T2 times. Our primary hypothesis was that the T2 times in the nucleus pulposus of the IVDs would improve after slow running (MAD: 0.4–0.7 G; running speed: 7–9 km/h), and that changes would be categorized into positive and negative effects depending on the direction of any of the triaxial accelerations.

## Materials and methods

### Participants

Our local ethics committee approved this study (no. 022-H173; January 31, 2023), which was performed in accordance with the guidelines for experimental studies involving human participants and met the ethical standards of the journal. The study was conducted in accordance with the principles of the Declaration of Helsinki. All participants were informed of the study’s purpose, experimental procedures, potential benefits, and possible risks, and were enrolled after providing written informed consent. Minor subjects were enrolled after obtaining parental consent. This study was conducted from February to May 2023. Sixteen healthy male university students were enrolled in this study (age, 20.8 ± 1.2 years; height, 172.2 ± 4.4 cm; weight, 72.3 ± 9.2 kg). The participants did not run daily and had no history of low back pain.

### Procedures

After arriving at the laboratory, the participants were briefed on the study, placed on their backs in bed, and allowed to rest for 10 min before pre-MRI testing was performed. Thereafter, patients were fitted with a triaxial accelerometer and performed a 1-minute run, followed by post-MRI testing. After running, the participants rested on their backs at all times, and MRI measurements were obtained 30 min after running to determine IVDs changes over time.

### Running exercise and triaxial accelerometer

Participants ran on a treadmill (Elite 5000; Johnson, Tokyo, Japan) for 1 min. The exercise duration in this study was selected based on the fact that the T2 values changed even with very low exercise loads (15 repetitions of trunk flexion, extension, and rotation exercises) [8]. The running speed was set at 8 km/h, within the range of MAD that had a positive effect on the IVDs based on the report by Bellaby et al. [9] No warm-up was performed before the run, and the experiment began as soon as the participants were ready. To measure the load caused by running, a triaxial accelerometer (wGT3X-BT; ActiGraph, Chiba, Japan) was fixed to the L4-L5 lumbar intervertebral region using a belt, and acceleration data were collected at 100 Hz. The x axis was oriented in the mediolateral direction. Acceleration was indicated on the X-axis as mediolateral acceleration, on the Y-axis as vertical acceleration, and on the Z-axis as anteroposterior acceleration.

### Testing and MRI scanning protocol

The participants were instructed not to perform any exercise on the day of the scan. T2 times before and after exercise were measured using MRI (ECHELON OVAL, 1.5 T; Hitachi Medical Systems, Tokyo, Japan). To avoid the effects of normal diurnal variations on the spine [12], all the tests were performed at approximately noon. The resulting images were evaluated by an orthopedic surgeon specializing in spinal diseases to assess IDD.

Sagittal T2 mapping using spin-echo multi-echo sequences (nine echo times: 12, 24, 36, 48, 60, 72, 84, 96, and 108 ms; repetition time, 1600 ms; number of slices, 12; slice thickness, 4 mm; gap, 5 mm; interslice distance, 5.0 mm; field of view, 200 mm; resolution, 1.04 × 1.25 × 4.00 mm per pixel; and acquisition time, 6 min 55 s).

### Acceleration and image data analysis

The recorded raw acceleration data collected during treadmill running were transferred to analysis software (OT BioLab, OT Bioelettronica, Italy) and filtered at a 10-Hz cut-off frequency. The rms of the X-, Y-axis-, and Z-axes, resultant acceleration ![Graphic][1]</img>, and mean amplitude deviation ![Graphic][2]</img> *x* = resultant acceleration) were calculated every 30 s on the analysis screen, excluding the 15 s before and after the 1-min running session.

T2 values were calculated using the MRI software. After the images were obtained in ECHELON OVAL format, they were loaded and viewed using the T2*RelaxMap feature of the MRI system (Hitachi Medical Systems). After the images of the nine echo times were superimposed using the software, the outlines of the lumbar IVDs were traced to determine the region of interest (ROI) in each image. The T2 times for each pixel were calculated using the software, whereas the T2 times were measured for each slice. The three IVDs analyzed were L3/L4 to L5/S1, which had a high incidence of IDD. The IVDs was segmented into five subregions from the anterior to posterior aspect (Fig 1) [13]: from front to back, the anterior annulus fibrosus (AF), anterior nucleus (AN), central nucleus (CN), posterior nucleus (PN), and posterior annulus fibrosus (PF). Image calculations were performed thrice for each ROI, and the average value was used as the T2 time.

![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/05/31/2024.05.30.24308204/F1.medium.gif)

[Fig 1.](http://medrxiv.org/content/early/2024/05/31/2024.05.30.24308204/F1)

Fig 1. 
Examples of IVDs and ROIs on MR images.

T2* RelaxMap of lumbar intervertebral discs in the central slice (left) and from L3/L4 to L5/S1 divided into five subregions (right): 1. anterior annulus fibrosus; 2. anterior nucleus; 3. central nucleus; 4. posterior nucleus; and 5. posterior annulus fibrosus.

### Statistical analysis

Data are expressed as mean ± standard deviation. Changes in T2 times and differences in IVDs levels were compared using one-way analysis of variance. Pearson’s correlation coefficient was used to examine the correlation between changes in T2 times and triaxial acceleration, synthetic acceleration, and MAD. The rate of change in T2 was determined by normalizing pre-MRI values. Where a correlation was found, the adjusted R-squared was determined using simple linear regression analysis. All tests were performed using statistical analysis software SPSS Statistics for Macintosh (ver. 29.0.2.0; IBM, Armonk, NY, USA). The significance level was set at p < .05.

## Results

### T2 times

Disc degeneration was characterized by the Pfirrmann grade [14]. Of the 16 subjects in the study, two discs were Pfirrmann grade III on L4/L5 and four discs were Pfirrmann grade III on L5/S1, and the corresponding IDD was excluded from this study. The T2 values for each IVDs level are shown in Table1. There were no significant differences in the changes over time in any of the ROI.

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/05/31/2024.05.30.24308204/T1)

Table 1. Changes in T2 relaxation times within the intervertebral disc pre-versus post-running

### Triaxial acceleration during running

The triaxial accelerations were 0.23 ± 0.06 rms in the X-axis (mediolateral), 1.37 ± 0.08 rms in the Y-axis (vertical), and 0.30 ± 0.06 rms in the Z-axis (anteroposterior). The resultant acceleration and MAD were 1.43 ± 0.07 rms and 0.71 ± 0.09 G, respectively.

### Correlation between T2 times and triaxial acceleration

The results of the correlation analysis between changes in T2 and acceleration are shown in Table 2. In the T2 times variation of L3/L4, negative correlations were observed for AN (r = -0.58, p = 0.02), CN (r = -0.54, p = 0.03) and PN (r = -0.59, p = 0.02) on the X-axis and positive correlations were observed for AN (r = 0.53, p = 0.04), CN (r = 0.62, p = 0.01) and PN (r = 0.72, p = 0.002) on the Z-axis. In addition, the Y-axis (r = 0.57, p = 0.02), resultant acceleration (r = 0.55, p = 0.03), and MAD (r = 0.68, p = 0.004) were positively correlated with PF. In the change in T2 time for L4/L5, a positive correlation was observed for CN (r = 0.73, p = 0.003), and a negative correlation was observed for PF at the post (r = -0.59, p = 0.03) on the Z-axis. In MAD, a negative correlation is observed for CN (r = -0.55, p = 0.04). The change in T2 times for L5/S1 showed a negative correlation with AF (r = -0.73, p = 0.01) and PN at the post-(r = -0.59, p = 0.05) and post30 (r = -0.69, p = 0.01) time points on the Z-axis. The strongest correlations were post30 PN (*R2* = 0.59) for L3-L4, CN (*R2* = 0.49) for L4-L5, and AF (*R2* = 0.48) for L5-S1 on the Z-axis. The regression equations were *y*=144.5 x +55.2 for the PF of L3/4, *y*=115.4 x +63.2 for the CN of L4/5, and *y*=-219.9 x +174.1 for the AF of L5/S1, respectively.

View this table:
[Table 2.](http://medrxiv.org/content/early/2024/05/31/2024.05.30.24308204/T2)

Table 2. Correlation between T2 relaxation time and triaxial acceleration

## Discussion

The positive effects of running loading on IVDs have been reported in recent years; however, the direction of loading by MAD in the optimal range is unclear. This study revealed that in slow running in the general population, the correlation between the change in T2 times and the direction of loading was different for each IVD level. In particular, in the anteroposterior direction, there were ROI with high adjusted R-square values at all IVDs levels, explaining approximately 50% of the changes in T2 times.

Previous studies have suggested that for IVDs, activities such as walking provide a “probable healthy” range of load based on the size, speed, direction, and type [15, 16]. Acute effects of running have been reported, with an average decrease of 1 mm in lumbar IVDs height after 1 h of running and a decrease in IVDs height and volume after moderate-intensity running [17, 18]. However, middle-aged men with >10 years of running experience showed differences in IVDs status, suggesting that long-term running may delay IDD [10]. These results suggest that the running load temporarily pushes water out of the IVDs; however, repeated loading in the optimal range may produce a favorable response in IVDs as a chronic change.

The IVDs in adults suggest that pumping by motion (fluid flow or convection) plays a role in transporting larger nutrients [19]. Thus, the transport of nutrients in IVDs requires motor loading, which may affect water movement. In this study, an optimal speed range of 8 km/h and running time were determined based on the reports of Chokan et al. and Bellaby et al. [8, 9]. Chokan et al. showed that in the nucleus pulposus, the T2 times significantly decreased after exercise and returned to the pre-exercise rest time. The results showed that there was a change in T2 times for trunk exercises (flexion-extension/rotation) and no change in T2 times for running, despite a similar exercise duration as in the present study, suggesting that T2 times are related to trunk range of motion and cause differences.

While T2 times may be affected by the range of motion, the direction of exercise loading at a running speed of 8 km/h suggests that loading in the anteroposterior direction of the Z-axis may affect T2 times by nearly 50%, based on the adjusted R-square value. Previous recommendations for beneficially modulating the IVDs include [15]: (a) use dynamic loading; (b) emphasize axial loading; and (c) choose exercises that load the IVDs in the range of 0.2–0.8 MPa, corresponding to an intradiscal pressure of 0.3–1.2 MPa [16]. In trunk motion, the greatest impact on intradiscal pressure is in flexion-extension [20], a movement that produces a range of 0.3-1.2 MPa, the optimal pressure for IVDs during the movement. Meadows et al. also observed that in their investigation of disc mechanics in disc strain, disc anterior-posterior shear, and mechanics due to flexion were related to the nucleus pulposus T2 values [21].

Therefore, IVDs may have been more susceptible to anteroposterior loading, and L3/L4 and L4/L5 [22], with greater flexion range of motion positively correlated with the application of optimal intradiscal pressure. In contrast, L5/S1 had greater posterior shear forces due to flexion than the other levels [21]. Shear forces may induce IDD [23, 24], and a negative correlation may be observed between different ranges of trunk motion and applied internal pressure at different IVDs levels.

Interestingly, L3/L4 was also affected by influences other than the Z-axis; L3/L4 had greater lateral flexion and rotation than the bottom two, and the X-axis was negatively correlated. This characteristic of cranial IVDs may be due to the high incidence of IDD in gymnasts with high rotational stress, which suggests that changes in T2 times are influenced by range of motion[25, 26].

Based on the present results, the hypothesis that slow running improves the T2 times of the nucleus pulposus was rejected because it did not significantly change the T2 times. However, the results partially support the hypothesis because it was improved by loading in the anteroposterior direction. This study had several limitations. In this study, a treadmill was used to maintain a constant running speed within the optimal loading range of the IVDs. However, the similarity in the triaxial acceleration data between treadmill and ground running suggests that the results for the loading direction were unaffected [27]. Although only male subjects were included in the study, sex was standardized because sex differences have been observed in pelvic movements during running [28], which may have affected the loading pattern. Finally, although there is a correlation, other factors may also have an effect. Further investigation is needed on factors other than speed that may contribute to IDD. The next investigation should be a slow-running motion analysis to determine the reasons for the high anteroposterior loading. Despite these limitations, the evidence suggests that IVDs are a beneficial response to certain types of loading and may have public health implications. IDD is an important factor in spinal pain. Understanding how IVDs respond to certain types of stress will lead to improved exercise guidelines for the prevention and management of spinal pain.

## Conclusion

The results of this study suggest that changes in the T2 times of the IVDs due to running at 8 km/h correlate with acceleration, suggesting that anteroposterior loading may contribute to the important response of the IVDs.

## Data Availability

All relevant data are within the manuscript and its Supporting Information files.

## Acknowledgements

The authors would like to thank all members of the research team and all participants of the study. This study was supported by the Japan Society for the Promotion of Science (, ) KAKENHI (30434153).

*   Received May 30, 2024.
*   Revision received May 30, 2024.
*   Accepted May 31, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Huang YC, Urban JP, Luk KD. Intervertebral disc regeneration: do nutrients lead the way? Nat Rev Rheumatol. 2014;10(9):561–6. Epub 20140610. doi: 10.1038/nrrheum.2014.91. PubMed PMID: 24914695.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrrheum.2014.91&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24914695&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

2.  2.Urban JP, Holm S, Maroudas A, Nachemson A. Nutrition of the intervertebral disk. An in vivo study of solute transport. Clin Orthop Relat Res. 1977;(129):101–14. PubMed PMID: 608268.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=608268&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

3.  3.Urban JP, Smith S, Fairbank JC. Nutrition of the intervertebral disc. Spine (Phila Pa 1976). 2004;29(23):2700-9. doi: 10.1097/01.brs.0000146499.97948.52. PubMed PMID: 15564919.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/01.brs.0000146499.97948.52&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15564919&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225481300013&link_type=ISI) 

4.  4.Grunhagen T, Wilde G, Soukane DM, Shirazi-Adl SA, Urban JP. Nutrient supply and intervertebral disc metabolism. J Bone Joint Surg Am. 2006;88 Suppl 2:30–5. doi: 10.2106/JBJS.E.01290. PubMed PMID: 16595440.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiamJqc2FtIjtzOjU6InJlc2lkIjtzOjEzOiI4OC9zdXBwbF8yLzMwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDUvMzEvMjAyNC4wNS4zMC4yNDMwODIwNC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

5.  5.Huang L, Liu Y, Ding Y, Wu X, Zhang N, Lai Q, et al. Quantitative evaluation of lumbar intervertebral disc degeneration by axial T2∗ mapping. Medicine. 2017;96(51).
    
    

6.  6.Marinelli NL, Haughton VM, Anderson PA. T2 relaxation times correlated with stage of lumbar intervertebral disk degeneration and patient age. AJNR Am J Neuroradiol. 2010;31(7):1278–82. Epub 20100401. doi: 10.3174/ajnr.A2080. PubMed PMID: 20360340; PubMed Central PMCID: PMCPMC7965459.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo5OiIzMS83LzEyNzgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wNS8zMS8yMDI0LjA1LjMwLjI0MzA4MjA0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

7.  7.Stelzeneder D, Welsch GH, Kovacs BK, Goed S, Paternostro-Sluga T, Vlychou M, et al. Quantitative T2 evaluation at 3.0T compared to morphological grading of the lumbar intervertebral disc: a standardized evaluation approach in patients with low back pain. Eur J Radiol. 2012;81(2):324–30. Epub 20110218. doi: 10.1016/j.ejrad.2010.12.093. PubMed PMID: 21315527.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejrad.2010.12.093&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21315527&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

8.  8.Chokan K, Murakami H, Endo H, Mimata Y, Yamabe D, Tsukimura I, et al. Evaluation of Water Retention in Lumbar Intervertebral Disks Before and After Exercise Stress With T2 Mapping. Spine (Phila Pa 1976). 2016;41(7):E430–6. doi: 10.1097/BRS.0000000000001283. PubMed PMID: 27018906.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/BRS.0000000000001283&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27018906&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

9.  9.Belavy DL, Quittner MJ, Ridgers N, Ling Y, Connell D, Rantalainen T. Running exercise strengthens the intervertebral disc. Sci Rep. 2017;7:45975. Epub 20170419. doi: 10.1038/srep45975. PubMed PMID: 28422125; PubMed Central PMCID: PMCPMC5396190.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/srep45975&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28422125&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

10. 10.Mitchell UH, Bowden JA, Larson RE, Belavy DL, Owen PJ. Long-term running in middle-aged men and intervertebral disc health, a cross-sectional pilot study. PLoS One. 2020;15(2):e0229457. Epub 20200221. doi: 10.1371/journal.pone.0229457. PubMed PMID: 32084224; PubMed Central PMCID: PMCPMC7034897.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0229457&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

11. 11.Takatalo J, Karppinen J, Nayha S, Taimela S, Niinimaki J, Blanco Sequeiros R, et al. Association between adolescent sport activities and lumbar disk degeneration among young adults. Scand J Med Sci Sports. 2017;27(12):1993–2001. Epub 20170319. doi: 10.1111/sms.12840. PubMed PMID: 28075521.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/sms.12840&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28075521&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

12. 12.Tyrrell A, Reilly T, Troup J. Circadian variation in stature and the effects of spinal loading. Spine. 1985;10(2):161–4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00007632-198503000-00011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=4002039&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1985AEV1600011&link_type=ISI) 

13. 13.Belavy DL, Brisby H, Douglas B, Hebelka H, Quittner MJ, Owen PJ, et al. Characterization of Intervertebral Disc Changes in Asymptomatic Individuals with Distinct Physical Activity Histories Using Three Different Quantitative MRI Techniques. J Clin Med. 2020;9(6). Epub 20200612. doi: 10.3390/jcm9061841. PubMed PMID: 32545639; PubMed Central PMCID: PMCPMC7357091.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/jcm9061841&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32545639&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

14. 14.Pfirrmann CW, Metzdorf A, Zanetti M, Hodler J, Boos N. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine (Phila Pa 1976). 2001;26(17):1873–8. doi: 10.1097/00007632-200109010-00011. PubMed PMID: 11568697.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00007632-200109010-00011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11568697&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170914000008&link_type=ISI) 

15. 15.Belavy DL, Albracht K, Bruggemann GP, Vergroesen PP, van Dieen JH. Can Exercise Positively Influence the Intervertebral Disc? Sports Med. 2016;46(4):473–85. doi: 10.1007/s40279-015-0444-2. PubMed PMID: 26666742.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s40279-015-0444-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26666742&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

16. 16.Chan SC, Ferguson SJ, Gantenbein-Ritter B. The effects of dynamic loading on the intervertebral disc. Eur Spine J. 2011;20(11):1796–812. Epub 20110504. doi: 10.1007/s00586-011-1827-1. PubMed PMID: 21541667; PubMed Central PMCID: PMCPMC3207351.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00586-011-1827-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21541667&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

17. 17.Dimitriadis AT, Papagelopoulos PJ, Smith FW, Mavrogenis AF, Pope MH, Karantanas AH, et al. Intervertebral disc changes after 1 h of running: a study on athletes. J Int Med Res. 2011;39(2):569–79. doi: 10.1177/147323001103900226. PubMed PMID: 21672362.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/147323001103900226&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21672362&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

18. 18.Kingsley MI, D’Silva LA, Jennings C, Humphries B, Dalbo VJ, Scanlan AT. Moderate-intensity running causes intervertebral disc compression in young adults. Med Sci Sports Exerc. 2012;44(11):2199–204. doi: 10.1249/MSS.0b013e318260dbc1. PubMed PMID: 22648342.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/MSS.0b013e318260dbc1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22648342&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

19. 19.Ferguson SJ, Ito K, Nolte LP. Fluid flow and convective transport of solutes within the intervertebral disc. J Biomech. 2004;37(2):213–21. doi: 10.1016/s0021-9290(03)00250-1. PubMed PMID: 14706324.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0021-9290(03)00250-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14706324&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000188498700007&link_type=ISI) 

20. 20.Wilke H, Neef P, Hinz B, Seidel H, Claes L. Intradiscal pressure together with anthropometric data--a data set for the validation of models. Clin Biomech (Bristol, Avon). 2001;16 Suppl 1:S111–26. doi: 10.1016/s0268-0033(00)00103-0. PubMed PMID: 11275349.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0268-0033(00)00103-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11275349&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

21. 21.Meadows KD, Peloquin JM, Newman HR, Cauchy PJK, Vresilovic EJ, Elliott DM. MRI-based measurement of in vivo disc mechanics in a young population due to flexion, extension, and diurnal loading. JOR Spine. 2023;6(1):e1243. Epub 20230109. doi: 10.1002/jsp2.1243. PubMed PMID: 36994458; PubMed Central PMCID: PMCPMC10041375.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jsp2.1243&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36994458&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

22. 22.Panjabi MM. The stabilizing system of the spine. Part II. Neutral zone and instability hypothesis. J Spinal Disord. 1992;5(4):390–6; discussion 7. doi: 10.1097/00002517-199212000-00002. PubMed PMID: 1490035.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00002517-199212000-00002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1490035&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992KC00800002&link_type=ISI) 

23. 23.Kim J, Yang SJ, Kim H, Kim Y, Park JB, Dubose C, et al. Effect of shear force on intervertebral disc (IVD) degeneration: an in vivo rat study. Ann Biomed Eng. 2012;40(9):1996–2004. Epub 20120420. doi: 10.1007/s10439-012-0570-z. PubMed PMID: 22527013.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10439-012-0570-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22527013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

24. 24.Xia DD, Lin SL, Wang XY, Wang YL, Xu HM, Zhou F, et al. Effects of shear force on intervertebral disc: an in vivo rabbit study. Eur Spine J. 2015;24(8):1711–9. Epub 20150318. doi: 10.1007/s00586-015-3816-2. PubMed PMID: 25784595.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00586-015-3816-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25784595&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

25. 25.Koyama K, Nakazato K, Min SK, Gushiken K, Hatakeda Y, Seo K, et al. Anterior Limbus Vertebra and Intervertebral Disk Degeneration in Japanese Collegiate Gymnasts. Orthop J Sports Med. 2013;1(3):2325967113500222. Epub 20130821. doi: 10.1177/2325967113500222. PubMed PMID: 26535240; PubMed Central PMCID: PMCPMC4555487.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/2325967113500222&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26535240&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

26. 26.Sward L, Hellstrom M, Jacobsson B, Nyman R, Peterson L. Disc degeneration and associated abnormalities of the spine in elite gymnasts. A magnetic resonance imaging study. Spine (Phila Pa 1976). 1991;16(4):437–43. doi: 10.1097/00007632-199104000-00009. PubMed PMID: 1828629.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00007632-199104000-00009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1828629&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991FG08900009&link_type=ISI) 

27. 27.Vanhelst J, Zunquin G, Theunynck D, Mikulovic J, Bui-Xuan G, Beghin L. Equivalence of accelerometer data for walking and running: treadmill versus on land. J Sports Sci. 2009;27(7):669–75. doi: 10.1080/02640410802680580. PubMed PMID: 19424900.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/02640410802680580&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19424900&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F05%2F31%2F2024.05.30.24308204.atom) 

28. 28.Perpiñá-Martínez S, Arguisuelas-Martínez MD, Pérez-Domínguez B, Nacher-Moltó I, Martínez-Gramage J. Differences between Sexes and Speed Levels in Pelvic 3D Kinematic Patterns during Running Using an Inertial Measurement Unit (IMU). International Journal of Environmental Research and Public Health. 2023;20(4):3631.

 [1]: /embed/inline-graphic-1.gif
 [2]: /embed/inline-graphic-2.gif