Long-term running in middle-aged men and intervertebral disc health, a cross-sectional pilot study

Abstract

Purpose

To measure intervertebral disc (IVD) health parameters in middle-aged long-term runners compared to matched non-physically active controls.

Methods

Seventeen males aged 44-62yr were included in the study: 9 runners with a running history of >10yr, averaging >50km/week, and eight matched non-physically active controls, the data from one participant had to be excluded. T2-relaxometry, diffusion weighted imaging, T1- and T2-weighted MR scanning, as well as T2 time mapping were performed. Morphological data relating to IVD were extrapolated.

Results

Compared to controls on average, runners had 20% greater IVD height (p = 0.002) and seven percentage points greater IVD-vertebral body height ratio (p = 0.001). No significant differences were observed between groups for mean(SD) IVD hydration status, as indicated by similar T2-times (runners: 94.4(11.1)ms, controls: 88.6(23.6)ms), or apparent diffusion coefficients (runners: 249.0(175.2)mm2/s, controls: 202.3(149.5)mm2/s). Average Pfirrmann score for the L5-S1 IVD was 2.2(0.7) for runners and 3.3(1.0) for controls (p = 0.026), average scores for all lumbar levels (L2-S1) were 1.9(0.2) and 2.5(0.7), respectively (p = 0.036). Anterior annulus T2-time and overall average lumbar level Pfirrmann grades were strongly correlated (r = 0.787, p = 0.021 and r = -0.704, p = 0.034, respectively) with greater distances run per week. Average lumbar level Pfirrmann grades were also strongly correlated (r = -0.823, p = 0.006) to total years of running.

Conclusion

Middle-aged long-term endurance runners exhibit less age-related decline in their lumbar IVDs. In addition, the measures of IVD morphology appeared to be better in those who had been running for a greater number of years, as well as in those who ran a greater distance per week.

About the Authors:

Ulrike H. Mitchell

Contributed equally to this work with: Ulrike H. Mitchell, Jennifer A. Bowden, Robert E. Larson, Daniel L. Belavy, Patrick J. Owen

Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Department of Exercise Sciences, Brigham Young University, Provo, Utah, United States of America

ORCID logo http://orcid.org/0000-0002-2298-8756

Jennifer A. Bowden

Contributed equally to this work with: Ulrike H. Mitchell, Jennifer A. Bowden, Robert E. Larson, Daniel L. Belavy, Patrick J. Owen

Roles Methodology, Writing – review & editing

Affiliation: Department of Exercise Sciences, Brigham Young University, Provo, Utah, United States of America

ORCID logo http://orcid.org/0000-0003-4959-675X

Robert E. Larson

Contributed equally to this work with: Ulrike H. Mitchell, Jennifer A. Bowden, Robert E. Larson, Daniel L. Belavy, Patrick J. Owen

Roles Data curation, Writing – review & editing

Affiliation: Department of Exercise Sciences, Brigham Young University, Provo, Utah, United States of America

Daniel L. Belavy

Contributed equally to this work with: Ulrike H. Mitchell, Jennifer A. Bowden, Robert E. Larson, Daniel L. Belavy, Patrick J. Owen

Roles Formal analysis, Software, Writing – review & editing

Affiliation: Deakin University, Institute for Physical Activity and Nutrition, School of Exercise and Nutrition Sciences, Geelong, Victoria, Australia

Patrick J. Owen

Contributed equally to this work with: Ulrike H. Mitchell, Jennifer A. Bowden, Robert E. Larson, Daniel L. Belavy, Patrick J. Owen

Roles Formal analysis, Software, Writing – review & editing

Affiliation: Deakin University, Institute for Physical Activity and Nutrition, School of Exercise and Nutrition Sciences, Geelong, Victoria, Australia

ORCID logo http://orcid.org/0000-0003-3924-9375

Introduction

Intervertebral disc (IVD) health is an important factor in overall spine health and the incidence of low back pain. It is strongly associated with nutrient flow within the disc [1]. In other tissues, nutrients are distributed to the cells through the cardiovascular system, however, IVD tissue is largely avascular in nature. While the outer annulus of the disc receives nutrients from the surrounding vasculature, the inner annulus and nucleus receive nutrients via bulk fluid flow within the disc tissue and diffusion through the vertebral endplates [2]. Multiple factors have been hypothesized to affect nutrient flow and thus disc health. For example, physical activity level and other lifestyle factors, such as high amounts of sedentary time, both in occupational and recreational time, influence disc nutrient flow [2–4]. Moreover, loading and resulting compression and deformation of the disc leads to expulsion of bulk fluid, while distraction encourages uptake of fluid and nutrition [5, 6]. As the IVD ages, its morphology changes and it loses its ability to attract water, rendering it more susceptible to injury [7]. After the age of 30 there is an increased rate of normal age-related IVD degeneration and the peak for these degenerative changes occurs around the age of 50 [8].

IVD health can be assessed using magnetic resonance imaging (MRI) [9]. T1 and T2 imaging techniques are useful for observing morphological changes to the IVD. In addition, the images can be used to grade and categorize the level of disc degeneration (Pfirrmann grading system [9]). Marinelli et al. [10] showed that T2 relaxation time, an intrinsic property of tissue that can be measured with MRI, significantly and strongly correlated with water content in the IVD. Belavý et al. [11] used T2 relaxation times to demonstrate that chronic runners aged 25-35yr exhibited better IVD tissue quality (i.e. hydration) compared to non-physically active controls. Diffusion-weighted imaging, a more advanced MRI application, has allowed clinicians and researchers to evaluate real time water movement within the IVD [12]. It gives a snapshot of the overall fluid movement within the disc and provides an earlier indication of changes in disc tissue than T2 relaxation time. Diffusion-weighted imaging can also be used to calculate the apparent diffusion coefficient (ADC), a value that reflects the average diffusion rate of water molecules [13].

Running is an activity that is of particular interest in examining disc health, as there are substantial compression and rotational forces placed on the lumbar disc with running, which likely impact disc health. Depending on the magnitude, frequency and duration of the forces, they can have an anabolic effect on the IVD [14]; however prolonged dynamic loading has been shown to bring about signs of degeneration [15]. Dimiatriadis et al. [16] showed that runners undergo an average disc height decrease of almost 1mm in each of the lumbar IVD after 1h of running. This could be interpreted as being detrimental to the structure of the disc, or as a benefit to the metabolic processes. Therefore, two competing theories concerning loading of the spine and its consequences to IVD health were considered in this study: 1) cyclic loading is beneficial for IVD tissue and leads to hypertrophic changes that make the disc stronger [11], and 2) mechanical overload produces localized trauma and tissue damage, which outpaces the ability of a disc to repair itself and leads to accelerated degradation [14]. Both theories are based on robust previous research, but the implications are conflicting.

The aim of this pilot project was to assess if IVD parameters obtained from MRI scans of middle-aged long-term runners and matched non-physically active controls are different and as such invite further studies. The particular IVD parameters were: T2-relaxation time, ADC, disc height (relative to vertebral body height, an internal control to body size [2]), and disc degeneration (Pfirrmann grade).

Methods

The study was approved by the institutional internal review board of Brigham Young University and participants provided their written informed consent prior to participation. A cross-sectional study was conducted from March to September 2018. We hypothesized that we will find different IVD parameters obtained from MRI scans for each of our two groups.

Description of participants

This study included men aged 44-62yr. By this age, IVD tissue alterations due to normal age-related IVD degeneration should be visible and discernible on MRI.

Inclusion criteria.

Male volunteers aged 44-62yr; the members in the running group had to have a history of running for >10yr, and had to have run on average >50km/wk. Members in the comparison group (control) had to have performed less than 150min/wk moderate physical activity, walk less than 15min to or from their place of work and not have performed regular sport or exercise training in the past 10yr.

Exclusion criteria.

Participants were excluded if they had current spinal pain greater than 4 on the digital pain scale, have a history of spinal surgery, history of traumatic injury to the spine, known scoliosis or kyphosis for which prior medical consultation was sought, and being a current or prior smoker. Women were excluded as they have a greater predisposition to back-related problems compared to men because of sex-specific events. These include accelerated post-menopausal disc degeneration [17] and increased propensity to vertebral fracture [18] in post-menopausal years. Participants were also excluded based on contraindications for MRI, such as known claustrophobia, metal object in the body, pacemaker or implantable cardioverter defibrillator.

Recruitment.

Participants were recruited by word of mouth. After running group participants were enrolled, we recruited an age-, sex- and height-matched control group. The controls were matched in order of priority to height (within 5cm), age (within 3yr), weight (within 3kg) and BMI (within 1kg/m2). Each subject had to match in at least two of the categories. All data were collected between noon and 3pm to account for diurnal changes in disc height and hydration status.

Measures

Magnetic resonance imaging (MRI), image processing and analysis.

We used a 3-Tesla Siemens MAGNETOM Tim Trio (Siemens Healthineers, Erlangen Germany) scanner with a 4-channel flexible coil. To quantify bone and IVD morphology, a sagittal T1-weighted scan was taken. For determination of IVD degeneration level (Pfirrmann grade) we took a sagittal T2-weighted turbo-spin echo sequence. A T2 mapping multi-echo spin echo (MESE) sequence was used to quantify T2 time, and a transverse two-dimensional diffusion weighted image was used to evaluate ADC for the L5-S1 IVD. Specific parameters for all MRI imaging are included in Table 1.

[Figure omitted. See PDF.]

Table 1. MRI parameters.

https://doi.org/10.1371/journal.pone.0229457.t001

After obtaining images in DICOM format, they were loaded and viewed on OsiriX (Pixmeo, Geneva, Switzerland). A mid-sagittal plane view was used and the IVDs of interest were manually segmented from the 2D T1 structural scan. T2 relaxation times and ADC values were determined for each IVD of interest.

Two investigators (UHM and REL) determined Pfirrmann grades of each lumbar IVD on sagittal T2-weighted images independently and in a blinded fashion. When the two scores did not coincide, the assessors discussed the cases and came to an agreement.

ImageJ 1.38x (http://rsb.info.nih.gov/ij/) was used to perform T2 time analysis. All available IVD from L2-S1 were measured. After segmenting the IVD, a custom written ImageJ plugin (“ROI Analyzer”; https://github.com/tjrantal/RoiAnalyzer) was used to measure area, height, width and signal intensity of the IVD in its entirety as well as in five subregions from the anterior to posterior aspect of the disc (Fig 1). T2-time was calculated via linear fit to the natural logarithm of the image intensity in each of the available MR echoes.

[Figure omitted. See PDF.]

Fig 1. Example diagram showing IVD ROIs on one MR image.

https://doi.org/10.1371/journal.pone.0229457.g001

Statistical analyses

All analyses were conducted using Stata statistical software version 15 (College Station TX, USA). Measures were compared between groups by one-way analysis of variance. In an exploratory analysis, the strength and direction of associations between running activity (runners only) and outcomes of interest were assessed with Pearson correlation coefficient. Analyses considered outcomes averaged across all lumbar vertebral levels (L2-S1) or each individual level, except ADC which used L5-S1 only [19]. An alpha-level of 0.05 was adopted for all statistical tests.

Results

Seventeen participants were included in analyses (runners: n = 9; control: n = 8). One participant in the control group complained of claustrophobia and we ended the scan pre-maturely. We had to consequently exclude his data. Participant demographics, including level of disability and health status, as measured by the Oswestry Disability questionnaire [20], health-related quality of life domains (36-Item Short Form Survey Australian version; SF36V2 [21]), and running activity are shown in Table 2. On average, controls had 9.7% greater body mass index (BMI) than runners, despite no significant difference in height or weight. The range of running activity parameters were: 10-39yr, 56-129km/wk, 4-6d/wk and 3.9–4.8min/km. No differences were observed in subjective measures of physical function or health-related quality of life domains.

[Figure omitted. See PDF.]

Table 2. Participants demographics and running activity in runners (n = 9) and controls (n = 8).

https://doi.org/10.1371/journal.pone.0229457.t002

Average lumbar level morphology of the vertebral body and IVD and outcomes of IVD health are displayed in Table 3. Compared to controls runners had significantly lower (4.4%) vertebral body height, greater (20%) IVD height and greater (seven percentage points) IVD-vertebral body height ratio. No differences were observed between groups for T2-time or ADC. Average Pfirrmann score was approximately half a grade significantly lower in runners compared to controls (Fig 2). Notably, additional analyses using linear mixed models with random effects for within participant variance across all lumbar levels did not alter these results (Table 3).

[Figure omitted. See PDF.]

Fig 2. Average and individual lumbar IVD Pfirrmann grades.

*P<0.05 when compared to control Average Pfirrmann grades for the whole lumbar spine and the L5/S1 segment were approximately half a grade significantly lower in runners compared to controls.

https://doi.org/10.1371/journal.pone.0229457.g002

[Figure omitted. See PDF.]

Table 3. Vertebral body and intervertebral disc (IVD) morphology and IVD health outcomes in runners (n = 9) and controls (n = 8).

https://doi.org/10.1371/journal.pone.0229457.t003

Lumbar level-specific vertebral body and IVD morphology are shown in Table 4. Runners had 20% and 30% significantly greater IVD height, for L2-L3 and L3-L4, respectively, compared to controls. Moreover, IVD-vertebral body height ratio was significantly greater in runners than controls for L2-L3 (36%), L3-L4 (25%) and L5-S1 (24%). Pfirrmann score was also one grade significantly lower for level L5-S1 when comparing runners to controls.

[Figure omitted. See PDF.]

Table 4. Vertebral body and intervertebral disc (IVD) morphology and IVD health outcomes in runners (n = 9) and controls (n = 8) by lumbar level.

https://doi.org/10.1371/journal.pone.0229457.t004

Associations between running activity and outcomes across all lumbar levels are displayed in Table 5. Greater total years of running (P = 0.006), kilometres ran per week (P = 0.034) and minutes per kilometre (P = 0.044) were associated with lower Pfirrmann grades. Greater kilometres ran per week was also associated with longer IVD anterior annulus T2-times (P = 0.021).

[Figure omitted. See PDF.]

Table 5. Correlations between variables examined in runners (n = 9).

https://doi.org/10.1371/journal.pone.0229457.t005

Discussion

The results of this study demonstrate that MRI scans of middle-aged long term endurance runners are different to scans from matched non-runners. The main finding was that long-term running has a positive effect on IVD health. Specifically, the current study showed that runners had lower average lumbar spine Pfirrmann grade and greater IVD-to-vertebral body height ratio. Moreover, these measures of IVD health appeared to be better in those who had been running for a greater number of years, as well as those who ran a greater distance per week.

The IVD of runners demonstrated less age-related degeneration, as measured by Pfirrmann grades [9]. The 5-point grading system is based on T2-weighted MR images and the structure of the IVD is taken into consideration. An IVD with a grade 1 represents a healthy, non-degenerated disc, while a grade 5 represents a severely degenerated IVD. In our study, runners exhibited healthier IVD, as demonstrated by lower Pfirrmann grades, compared to the non-physically active control participants across all lumbar levels, as well as specifically at the L5-S1 IVD. This finding is of particular interest, because the lower lumbar levels (i.e. L4-L5 and L5-S1) have been found to be the most common sites for disc degeneration[22] and disc herniation[23]. In contrast to our findings, Belavý et al. [11] found that there was no difference in Pfirrmann grades between participants who were not regularly involved in sport compared to runners. Their population was an average of 20yr younger, and morphological changes due to age and activity may have not yet manifested themselves (at least not on MR images). Our findings suggest that long-term running may delay IVD degeneration or possibly result in beneficial modulation of this tissue.

Runners in the current study had a greater IVD-to-vertebral body height ratio compared to their non-physically active counterparts. The difference in ratios was a function of both a smaller vertebral body height and a larger IVD height in the runners. Belavý et al. [11] also found that long distance runners, but not the joggers, nor the participants in the non-sporting group, exhibited a greater IVD-to-vertebral body height ratio and referred to this as ‘IVD hypertrophy’. This finding is interesting because it is commonly accepted that lumbar IVDs decrease in height with advancing age [24] as a result of morphology changes and the reduced ability to bind water [7]. On the other hand, endurance running may just be delaying degeneration and not actually causing disc growth. Somewhat related to this are findings by Videman et al.[25] who reported that the heavier twin in weight discordant monozygotic twins exhibited slightly larger IVDs compared to their lighter twin. This was explained by the greater habitual loading the heavier twin experienced, a concept similarly proposed to explain why individuals with greater fat mass often have greater absolute muscle mass. In our study the vertebral body height was smaller in runners compared to their non-sporting counterparts. Therefore, long-term running (i.e. habitual dynamic axial loading of the spine) may have elicited a bone response with possible physiologic changes in bone architecture. We do not have data on bone architecture, nor density, of any of the vertebrae, thus our suggestions that the decrease in vertebral height was coupled with increased bone density remains speculative. Nonetheless, these morphological observations in the current study support that running is associated with IVD hypertrophy and thus better IVD health.

Both greater total years of running and total distance ran per week were associated with lower Pfirrmann grades, which indicates that greater amounts of running may delay IVD degeneration. Similarly, there was a positive association between running distance and T2-times, indicating greater IVD hydration, albeit this was only observed in the annulus. The annulus contains substantially fewer water binding proteoglycans compared to the nucleus, which makes this finding remarkable.

Surprisingly, parameters in the current study indicating hydration of the disc (i.e. T2-time and ADC) did not differ between groups, despite differences in Pfirrmann grade. Previously, Belavý et al. [11] found significantly longer T2-times in runners compared to non-athletic individuals. The discrepancy between our findings may be explained by the age difference of our participants. Our study used participants with an average age of 50yr (20yr older than comparator), when disc hydration had been found to be greatly reduced [26]. We intentionally recruited middle-aged participants, to assess if habitual cyclic loading of the disc was able to delay the desiccation process associated with ageing [7]. However, the data from our current study do not support the notion that running delays this process. Nevertheless, our data indicate longer average T2 times in all areas of the IVD and at all levels in the runners compared to the controls, except for the anterior annulus. This difference is not statistically significant, most likely because the data are likely underpowered and have high variability.

The results of this cross-sectional study confirmed previous findings [11] that there is an association between long-term long-distance running and better IVD health. They also support earlier speculations that the IVD is able to undergo anabolic adaptations within a human life span [11]. The question of an optimal loading pattern, however, remains elusive. The following information should be used to elucidate this matter in further studies: The relationship between habitual physical activity and risk of IVD degeneration has been described as U-shaped [27], indicating that either too much or insufficient loading of the spine is detrimental to IVD health. In addition, it has been suggested that there might be a ‘physiological range’, or ‘zone’ of dynamic axial compressive loads that lead to an anabolic response within the IVD [27]. This also implies that there might exist a ‘therapeutic window’ of optimal benefits one can receive from running. Establishing these parameters has potential public health implications as IVD degeneration is a leading cause of mechanical low back pain. Therefore, interventional exercise training studies examining the efficacy of running for IVD health are warranted.

Our study was strengthened by the blinded nature of MRI analyses, as well as the matched-nature of controls. Cross-sectional studies, however, have well known limitations. Despite observations that running was associated with better IVD health, we cannot draw a cause and effect conclusion (e.g. running is the cause for healthy IVDs). Similarly, we cannot dismiss the possibility that there is reverse causality (e.g. healthier discs allow people to continue running, while degenerated discs force people to stop). Another limitation was the variance in running activity (e.g. range of years: 10-39yr). Lastly, while we explained why only men were included in this study, these findings may not apply to women.

Conclusion

Middle-aged long-term endurance runners exhibit less age-related decline in their lumbar IVDs compared to matched non-runners. In addition, the measures of IVD morphology appeared to be better in those who had been running for a greater number of years, as well as in those who ran a greater distance per week. Our findings support those of the only other study that investigated the effect of long term running on the IVD [11] and broaden the implications to an older and more active population.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Supporting data spreadsheet.

https://doi.org/10.1371/journal.pone.0229457.s001

(XLSX)

Citation: Mitchell UH, Bowden JA, Larson RE, Belavy DL, Owen PJ (2020) Long-term running in middle-aged men and intervertebral disc health, a cross-sectional pilot study. PLoS ONE 15(2): e0229457. https://doi.org/10.1371/journal.pone.0229457

References

1. Guiot BH, Fessler RG. Molecular biology of degenerative disc disease. Neurosurgery. 2000;47(5):1034–40. pmid:11063096

2. 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.

3. Bowden JA, Bowden AE, Wang H, Hager RL, LeCheminant JD, Mitchell UH. In Vivo Correlates Between Daily Physical Activity and Intervertebral Disc Health. J Orthop Res. 2017.

4. Adams MA, Hutton WC. The effect of posture on the fluid content of lumbar intervertebral discs. Spine (Phila Pa 1976). 1983;8(6):665–71.

5. Gullbrand SE, Peterson J, Ahlborn J, Mastropolo R, Fricker A, Roberts TT, et al. ISSLS Prize Winner: Dynamic Loading-Induced Convective Transport Enhances Intervertebral Disc Nutrition. Spine (Phila Pa 1976). 2015;40(15):1158–64.

6. Mitchell UH, Beattie PF, Bowden J, Larson R, Wang H. Age-related differences in the response of the L5-S1 intervertebral disc to spinal traction. Musculoskelet Sci Pract. 2017;31:1–8. pmid:28624722

7. Urban JP, Roberts S. Degeneration of the intervertebral disc. Arthritis Res Ther. 2003;5(3):120–30. pmid:12723977

8. Boos N, Weissbach S, Rohrbach H, Weiler C, Spratt KF, Nerlich AG. Classification of age-related changes in lumbar intervertebral discs: 2002 Volvo Award in basic science. Spine (Phila Pa 1976). 2002;27(23):2631–44.

9. 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.

10. Marinelli NL, Haughton VM, Munoz A, Anderson PA. T2 relaxation times of intervertebral disc tissue correlated with water content and proteoglycan content. Spine (Phila Pa 1976). 2009;34(5):520–4.

11. Belavy DL, Quittner MJ, Ridgers N, Ling Y, Connell D, Rantalainen T. Running exercise strengthens the intervertebral disc. Sci Rep. 2017;7:45975. pmid:28422125

12. Zhang Z, Chan Q, Anthony MP, Samartzis D, Cheung KM, Khong PL, et al. Age-related diffusion patterns in human lumbar intervertebral discs: a pilot study in asymptomatic subjects. Magn Reson Imaging. 2012;30(2):181–8. pmid:22055854

13. Beattie PF, Morgan PS, Peters D. Diffusion-weighted magnetic resonance imaging of normal and degenerative lumbar intervertebral discs: a new method to potentially quantify the physiologic effect of physical therapy intervention. J Orthop Sports Phys Ther. 2008;38(2):42–9. pmid:18560192

14. Iatridis JC, MacLean JJ, Roughley PJ, Alini M. Effects of mechanical loading on intervertebral disc metabolism in vivo. J Bone Joint Surg Am. 2006;88 Suppl 2:41–6.

15. Fearing BV, Hernandez PA, Setton LA, Chahine NO. Mechanotransduction and cell biomechanics of the intervertebral disc. JOR Spine. 2018;1(3).

16. 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. pmid:21672362

17. Lou C, Chen H, Mei L, Yu W, Zhu K, Liu F, et al. Association between menopause and lumbar disc degeneration: an MRI study of 1,566 women and 1,382 men. Menopause. 2017.

18. Puche RC, Morosano M, Masoni A, Perez Jimeno N, Bertoluzzo SM, Podadera JC, et al. The natural history of kyphosis in postmenopausal women. Bone. 1995;17(3):239–46. pmid:8541136

19. Beattie PF, Donley JW, Arnot CF, Miller R. The Change in the Diffusion of Water in Normal and Degenerative Lumbar Intervertebral Discs Following Joint Mobilization Compared to Prone Lying. Journal of Orthopaedic & Sports Physical Therapy. 2009;39(1):4–11.

20. Fairbank JC, Couper J, Davies JB, O'Brien JP. The Oswestry low back pain disability questionnaire. Physiotherapy. 1980;66(8):271–3. pmid:6450426

21. Hawthorne G, Osborne RH, Taylor A, Sansoni J. The SF36 Version 2: critical analyses of population weights, scoring algorithms and population norms. Qual Life Res. 2007;16(4):661–73. pmid:17268926

22. Teraguchi M, Yoshimura N, Hashizume H, Muraki S, Yamada H, Minamide A, et al. Prevalence and distribution of intervertebral disc degeneration over the entire spine in a population-based cohort: the Wakayama Spine Study. Osteoarthritis and cartilage / OARS, Osteoarthritis Research Society. 2014;22(1):104–10.

23. Ma D, Liang Y, Wang D, Liu Z, Zhang W, Ma T, et al. Trend of the incidence of lumbar disc herniation: decreasing with aging in the elderly. Clin Interv Aging. 2013;8:1047–50. pmid:23966775

24. Ferguson SJ, Steffen T. Biomechanics of the aging spine. Eur Spine J. 2003;12 Suppl 2:S97–S103.

25. Videman T, Gibbons LE, Kaprio J, Battie MC. Challenging the cumulative injury model: positive effects of greater body mass on disc degeneration. Spine J. 2010;10(1):26–31. pmid:19926343

26. Oh CH, Yoon SH. Whole spine disc degeneration survey according to the ages and sex using Pfirrmann disc degeneration grades. Korean J Spine. 2017;14(4):148–54. pmid:29301175

27. 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. pmid:26666742