Advertisement
Research Article| Volume 22, ISSUE 11, P1840-1850, November 2014

Download started.

Ok

Biochemical cartilage alteration and unexpected signal recovery in T2* mapping observed in ankle joints with mobile MRI during a transcontinental multistage footrace over 4486 km

Open ArchivePublished:August 12, 2014DOI:https://doi.org/10.1016/j.joca.2014.08.001

      Summary

      Objective

      The effect of ultra-long distance running on the ankle cartilage with regard to biochemical changes, thickness and lesions is examined in the progress of a transcontinental ultramarathon over 4486 km.

      Method

      In an observational field study, repeated follow-up scanning of 22 participants of the TransEurope FootRace (TEFR) with a 1.5 T MRI mounted on a mobile unit was performed. For quantitative biochemical and structural evaluation of cartilage a fast low angle shot (FLASH) T2* weighted gradient-echo (GRE)-, a turbo-inversion-recovery-magnitude (TIRM)- and a fat-saturated proton density (PD)-weighted sequence were utilized. Statistical analysis of cartilage T2* and thickness changes was obtained on the 13 finishers (12 male, mean age 45.4 years, BMI 23.5 kg/m²). None of the nine non-finisher (eight male, mean age 53.8 years, BMI 23.4 kg/m²) stopped the race due to ankle problems.

      Results

      From a mean of 17.0 ms for tibial plafond and 18.0 ms for talar dome articular cartilage at baseline, nearly all observed regions of interest (ROIs) of the ankle joint cartilage showed a significant T2*-signal increase (25.6% in mean), with standard error ranging from 19% to 33% within the first 2500 km of the ultra-marathon. This initial signal behavior was followed by a signal decrease. This signal recovery (30.6% of initial increase) showed a large effect size. No significant morphological or cartilage thickness changes (at baseline 2.9 mm) were observed.

      Conclusion

      After initial T2*-increase during the first 2000–2500 km, a subsequent T2*-decrease indicates the ability of the normal cartilage matrix to partially regenerate under ongoing multistage ultramarathon burden in the ankle joints.

      Keywords

      Introduction

      Societal life style changes in the last 20 years have created a distinct divide with one group reaching for higher and higher levels of fitness and the other resigning to a more sedentary daily life. There is abundant research on consequences of obesity and sedentary lifestyles for the human body and lower extremity joints in particular
      • Booth F.W.
      • Roberts C.K.
      • Laye M.J.
      Lack of exercise is a major cause of chronic diseases.
      . Little or contradictory information is available about the impact of high levels of endurance running exercise on the joints of the legs
      • Schueller-Weidekamm C.
      • Schueller G.
      • Uffmann M.
      • Bader T.
      Incidence of chronic knee lesions in long-distance runners based on training level: findings at MRI.
      • Krampla W.W.
      • Newrkla S.P.
      • Kroener A.H.
      • Hruby W.F.
      Changes on magnetic resonance tomography in the knee joints of marathon runners: a 10-year longitudinal study.
      . The 4486 km long ultramarathon TransEurope FootRace (TEFR)
      • Schütz U.H.
      • Schmidt-Trucksäss A.
      • Knechtle B.
      • Machann J.
      • Wiedelbach H.
      • Ehrhardt M.
      • et al.
      The TransEurope Footrace Project: longitudinal data acquisition in a cluster randomized mobile MRI observational cohort study on 44 endurance runners at a 64-stage 4,486 km transcontinental ultramarathon.
      where the participants were continuously monitored using a mobile magnetic resonance imaging (MRI) traveling along on a semitrailer truck offered a unique, once in a life time opportunity for investigating the physiological responses to extreme contiguous 64 day exercise without any day of rest. How much exercise is beneficial for the human body and where is the limit for physiological adaptation of the ankle joint, if there is one, when human beings are running more than 100 marathons, one after another? Only very few data are published concerning the abnormalities and functional adaptation of the ankle joint (tibiotalar joint; TTJ) cartilage to mechanical loadings
      • El-Khoury G.Y.
      • Alliman K.J.
      • Lundberg H.J.
      • Rudert M.J.
      • Brown T.D.
      • Saltzman C.L.
      Cartilage thickness in cadaveric ankles: measurement with double-contrast multi-detector row CT arthrography versus MR imaging.
      • Trattnig S.
      • Breitenseher M.J.
      • Huber M.
      • Zettl R.
      • Rottmann B.
      • Haller J.
      • et al.
      Determination of cartilage thickness in the ankle joint. An MRT (1.5)-anatomical comparative study.
      , nothing is known about the consequences of ultra-long distance running for the cartilage of the lower extremity joints. In this setting, a strategic decision had to be made for just one of the “quantitative cartilage mapping techniques” to fit the severe time constraints, as whole body, brain and cardiovascular studies also were performed. The ankle joint is uniquely designed for repetitive and high biomechanical loading on a much smaller surface area than the knee with a “mortise” type fit. Furthermore, the articular cartilage is much thinner than in the knee
      • Koepp H.
      • Eger W.
      • Muehleman C.
      • Valdellon A.
      • Buckwalter J.A.
      • Kuettner K.E.
      • et al.
      Prevalence of articular cartilage degeneration in the ankle and knee joints of human organ donors.
      . Therefore, to obtain reasonable signal to noise and the necessary resolution to image a thin, curved small structure in a short acquisition, T2* relaxation time measurements were chosen
      • Hesper T.
      • Hosalkar H.S.
      • Bittersohl D.
      • Welsch G.H.
      • Krauspe R.
      • Zilkens C.
      • et al.
      T2* mapping for articular cartilage assessment: principles, current applications, and future prospects.
      . This measurement offers inherently higher signal-to-noise ratio and robustness when compared to other options such as T2, T1rho, chemical exchange saturation transfer (CEST), delayed gadolinium-enhanced MRI of cartilage (dGEMRIC), Na measurements etc
      • Koepp H.
      • Eger W.
      • Muehleman C.
      • Valdellon A.
      • Buckwalter J.A.
      • Kuettner K.E.
      • et al.
      Prevalence of articular cartilage degeneration in the ankle and knee joints of human organ donors.
      . As such, it was considered to be the only viable option. We hypothesize that the ankle joint will reveal overall superior resilience to degeneration
      • Cole A.A.
      • Margulis A.
      • Kuettner K.E.
      Distinguishing ankle and knee articular cartilage.
      , unless a pre-existing injury is present. Our goal in this study was to explore the physiological adaptive capacity of the TTJ cartilage under conditions of extreme wear and address the question how much is too much?

      Method

      Of the 67 runners in the 4486 km ultramarathon, 44 volunteers who were officially accepted as TEFR participants by the organizers
      • Schütz U.H.
      • Schmidt-Trucksäss A.
      • Knechtle B.
      • Machann J.
      • Wiedelbach H.
      • Ehrhardt M.
      • et al.
      The TransEurope Footrace Project: longitudinal data acquisition in a cluster randomized mobile MRI observational cohort study on 44 endurance runners at a 64-stage 4,486 km transcontinental ultramarathon.
      , did not have MRI contraindications and gave written informed consent
      • Shellock F.G.
      • Spinazzi A.
      MRI safety update 2008: part 2, screening patients for MRI.
      were included in the scientific TEFR project. The local ethics committee approved the study in accordance to the Declaration of Helsinki. The project included multiple MRI studies, such as brain, cardiac, vascular and musculoskeletal imaging at baseline and four consecutive time points approximately every 900 km in the race. Musculoskeletal imaging was randomized to 22 runners each (20 male, two female) in a knee MRI (reference) and an ankle/foot MRI group, results reported here, respectively
      • Schütz U.H.
      • Schmidt-Trucksäss A.
      • Knechtle B.
      • Machann J.
      • Wiedelbach H.
      • Ehrhardt M.
      • et al.
      The TransEurope Footrace Project: longitudinal data acquisition in a cluster randomized mobile MRI observational cohort study on 44 endurance runners at a 64-stage 4,486 km transcontinental ultramarathon.
      . All 22 randomly selected volunteers in the ankle cohort met the inclusion criteria which were official acceptance as a participant at the TEFR by the organizers
      • Schütz U.H.
      • Schmidt-Trucksäss A.
      • Knechtle B.
      • Machann J.
      • Wiedelbach H.
      • Ehrhardt M.
      • et al.
      The TransEurope Footrace Project: longitudinal data acquisition in a cluster randomized mobile MRI observational cohort study on 44 endurance runners at a 64-stage 4,486 km transcontinental ultramarathon.
      and absence of common contraindications against native MRI
      • Shellock F.G.
      • Spinazzi A.
      MRI safety update 2008: part 2, screening patients for MRI.
      scanning. Subsequently, 13 (59.1%) of the 22 participants in the ankle/foot MRI study group finished the race with a baseline and all four consecutive MRI data time points necessary for the data analysis. Main reason for prematurely voluntary withdraw from race were stress reaction to the lower extremity soft tissues including the muscles, tendons, bones and subcutaneous tissues unrelated to the ankle or foot joints (6). The data of the study subjects (13 finisher/9 non-finisher) at baseline were: mean age 45.4/53.8 years (SD 10.7/11.3, range 27–62/36–68), male 12/8 (92.7/88.9), body mass (BM) 73.0/67.9 kg (SD 11.3/11.3, range 51.9–94.2/49.2–81.8), body mass index (BMI) 23.4/23.5 kg/m² (SD 2.5/3.0, range 20.5–29.1/19.2–28.3). Age, gender, BMI and clinical alignment of the lower extremity joints showed no relevant differences between study subjects and non-participants
      • Schütz U.H.
      • Schmidt-Trucksäss A.
      • Knechtle B.
      • Machann J.
      • Wiedelbach H.
      • Ehrhardt M.
      • et al.
      The TransEurope Footrace Project: longitudinal data acquisition in a cluster randomized mobile MRI observational cohort study on 44 endurance runners at a 64-stage 4,486 km transcontinental ultramarathon.
      .
      MRI data were acquired with a mobile 1.5 T MR scanner (Avanto™, Siemens Ltd., Erlangen), which was mounted on a MRI-semitrailer truck traveling with the runners throughout whole TEFR. Their ankles were scanned consecutively using a table fixed, 8-channel ankle and foot array coil with a boot-like design ensuring a standardized foot position. MRI scanning was planned at baseline (t0) within the last 4 days before start and roughly every 900 km (±211.5 km) measurement interval (MI: t1–t3) during the race and t4 representing the endpoint.
      For T2*-mapping a sagittal fast low angle shot (FLASH) T2* gradient-echo (GRE) sequence was used: flip angle (FA) 60°, echo times (TE) 4.5/12.2/19.9/27.7/35.4 ms, repetition time (TR) 1010 ms, slice thickness (ST) 2.5 mm, field of view (FOV) 182.25 cm², pixel size (PS) 0.178 mm² (in plane resolution). For quantitative biochemical cartilage analysis T2*-relaxation times were obtained from online reconstructed T2*-maps by using a pixel wise, monoexponential nonnegative least squares fit analysis (syngo™ MapIt; Siemens Ltd.)
      • Apprich S.
      • Welsch G.H.
      • Mamisch T.C.
      • Szomolanyi P.
      • Mayerhoefer M.
      • Pinker K.
      • et al.
      Detection of degenerative cartilage disease: comparison of high-resolution morphological MR and quantitative T2 mapping at 3.0 Tesla.
      . For detection of soft tissue pathology a turbo inversion recovery (TIRM) sequence was utilized in the sagittal plane: FA 140°, TE 60 ms, inversion time 120 ms, ST 2 mm, FOV 900 cm², PS 0.343 mm². Additional proton density weighted fat saturated (PDfs) sequence was obtained if a suspected subchondral signal abnormality was detected: FA 150°, TE 32 ms, TR 5830 ms, ST 3 mm, FOV 256 cm², PS 0.172 mm².
      The articular cartilage of the ankle joint was defined as normal if the cartilage thickness was preserved, no cartilage signal alterations and no superficial and deep cartilage defects or fissures were present on the respective morphological imaging sequences (TIRM, PDfs, T2*-GRE).
      For image post processing three sagittal slices through TTJ centered between the medial and lateral margins of the talar dome and two additional slices 8.2 mm medially and laterally, respectively. Six regions of interest (ROIs) for analysis of mean T2*-values were manually drawn on the three slices to cover the entire tibial plafond and talar dome cartilage: tibial-medial, tibial-central, tibial-lateral, talar-medial, talar-central, talar-lateral [Fig. 1(B)]. Cartilage thickness was measured in the center of the anterior, mid and posterior third of each of the three sagittal slices [Fig. 1(C)]. All measurements were done by three experienced musculoskeletal radiologists (UHS, DS, BC).
      Figure thumbnail gr1
      Fig. 1MR-image post-processing for quantification of thickness and T2* relaxation time of TTJ cartilage: A: sagittal FLASH T2*w GRE: (1) tibia, (2) talus, (3) calcaneus, (4) navicular, (5) cuboid, (6) ankle joint, (7) subtalar joint, (8) talonavicular joint, (9) calcaneocuboidal joint. B: fused colored T2* GRE map (syngo™ MapIt fusion technique): colored visualization of T2* in tibial plafond and talar dome cartilage of the ankle joint. C: cartilage thickness [mm] measurement of ankle joint (anterior, central and posterior in each slice).
      All osteochondral lesions were independently classified according to the modified Outerbridge MRI grading system
      • Mosher T.J.
      MRI of osteochondral injuries of the knee and ankle in the athlete.
      . In lesions greater Outerbridge grade 1 the ROI for quantitative T2*-analysis was drawn in the subjacent cartilage slice.
      For data documentation and statistical analysis Microsoft™ Office Excel™ (Microsoft Inc.) and SPSS™ (IBM™ Statistics, SPSS Inc.) were utilized, respectively.
      The absolute thickness values are presented (Table I) as means and standard deviation (SD) with min and max. The T2*-values are presented as absolute values and relative differences to baseline for the finisher cohort.
      Table IThickness measurements (at baseline) and throughout TEFR (t0–t4): mean (SD), range (n = 13)
      MILateral sliceCentral sliceMedial slice
      AnteriorCentralPosteriorAnteriorCentralPosteriorAnteriorCentralPosterior
      t02.38 (0.39) 1.7–3.12.87 (0.40) 2.1–3.83.44 (0.56) 2.3–5.02.62 (0.48) 1.8–3.52.87 (0.44) 1.7–3.52.98 (0.42) 2.2–3.93.32 (0.48) 2.4–4.53.12 (0.40) 2.2–3.72.94 (0.48) 2.0–4.2
      t12.5 (0.34) 2–3.12.88 (0.31) 2.3–3.63.5 (0.45) 2.9–4.82.69 (0.4) 2–3.42.91 (0.4) 2.1–3.63.01 (0.46) 2.3–4.13.4 (0.47) 2.4–4.63.2 (0.39) 2.3–3.72.92 (0.39) 2.4–3.5
      t22.42 (0.37) 1.8–3.12.79 (0.33) 2.2–3.53.4 (0.46) 2.5–4.72.7 (0.42) 2–3.42.85 (0.35) 2.2–3.52.96 (0.45) 2.3–3.93.4 (0.56) 2.3–4.63.13 (0.39) 2.1–3.72.96 (0.46) 2.1–4.1
      t32.53 (0.37) 1.8–3.22.81 (0.33) 2.3–3.63.38 (0.43) 2.4–4.42.74 (0.38) 2.1–3.32.9 (0.36) 2.1–3.62.92 (0.41) 2.1–3.73.42 (0.46) 2.6–4.53.17 (0.33) 2.2–3.82.97 (0.44) 2.3–4.1
      t42.48 (0.35) 1.8–3.12.76 (0.36) 2.3–3.63.33 (0.45) 2.4–4.52.68 (0.35) 2.1–3.32.83 (0.39) 2.1–3.62.85 (0.40) 2.1–3.83.4 (0.45) 2.4–4.43.14 (0.28) 2.4–3.72.98 (0.43) 2.4–4
      For determination of intraobserver precision of ROIs, mean T2*- and thickness-values, all measurements of the right side were performed at baseline (t0) and at a mean interval of 3 weeks by one investigator (UHS). For interrater reliability baseline measurements of two radiologists (UHS and DS) were compared. For precision calculation the 95% limits of agreement (LOA: mean difference ± 1.96 SD)
      • Bland J.M.
      • Altman D.G.
      Statistical methods for assessing agreement between two methods of clinical measurement.
      , for calculation of correlation coefficients on rater reliabilities lambda as proposed by Jepsen et al.
      • Jepsen J.R.
      • Laursen L.H.
      • Hagert C.G.
      • Kreiner S.
      • Larsen A.I.
      Diagnostic accuracy of the neurological upper limb examination I: inter-rater reproducibility of selected findings and patterns.
      was utilized.
      A tests P-value of 0.05 indicated significance. For testing on differences between T2*-relaxation times of tibial plafond and talar dome ROIs within the same slice a t test for independent variables was used. Differences between T2*-relaxation time of ROIs within the same cartilage layer at the same time point were calculated by a 1-way analysis of variance (ANOVA) without repeated measurements. To analyze significant changes of T2*- and thickness-values between t0 and MI (t1–t4) during TEFR a 1-way ANOVA for repeated measurements utilized. Given the longitudinal nature of the test data, a general linear model for repeated measurements was applied. For correction of accumulation of the alpha level due to multiple testing a Bonferroni-procedure was applied. The precondition sphericity was proven by the Mauchly-Test. To determine a trend of T2*-value curves in progress of TEFR significance of innersubject effects of the ANOVA were calculated. To determine significant value differences at the end of TEFR (t4) to any maximal change (peak) during exercise (t1–t3) compared to baseline, a paired (two-tailed) samples t test with calculation of the effect size according to Cohen
      • Cohen J.
      Statistical Power Analysis for the Behavioral Sciences.
      was done.
      If relevant changes in the course of TEFR were found, relationship to influencing cofactors of running burden (Table I) and BM was proven with a specific regression analysis using a linear mixed model for fixed effects.

      Results

      Combined tibial plafond and talar dome mean cartilage thickness ranged from 2.4 to 3.5 mm and there was no significant change in articular cartilage thickness during the race (Table I).
      Drawn ROI sizes ranged from 32.3 mm² (lateral tibial) to 38.0 mm² (medial tibial). Intra- and interrater analysis showed high agreement for drawn ROI sizes with 7.5% (2.7 mm²) and 8.3% (3.0 mm²), for cartilage thickness with 5.5% (0.16 mm) and 6.3% (0.19 mm), and for T2* values with 3.1% (0.5 ms) and 3.3% (0.6 ms) respectively. The intra-and interrater reliabilities (λ) were high for all ROI sizes ranging from 0.962 to 0.991 and 0.965 to 0.992, for measured thickness ranging from 0.975 to 0.991 and 0.969 to 0.990 and for T2* values ranging from 0.995 to 0.998 and 0.994 to 0.998 at baseline respectively.
      In the finisher cohort the mean baseline T2*-values for the right and left ankle were 16.1 ms (SD 2.4) and 17.8 ms (SD 2.7) for tibial plafond and 17.5 ms (SD 3.1) and 18.4 ms (SD 3.3) for talar dome cartilage, respectively. Figure 2 demonstrates the mean, confidence interval (CI) and percentiles for each ROI, the talar dome revealed a higher T2*-value (1.1% in mean) compared to the tibial cartilage at baseline, but this was only significant in the central ROIs at t1 and t2, in the medial ROIs of the left side throughout total TEFR (t1–t4) and on the right side only at t2 (Table II). Tibial and talar cartilage T2*-values combined were significantly different on the left side at t2, t3, and t4 and on the right only at t2.
      Figure thumbnail gr2
      Fig. 2T2*-mapping of TTJ (relaxation time): Absolute values of measures of variation (nF = 13).
      Table IIP-values of mean T2* values [ms] of articular ROIs over time
      SideRegionP-value
      t0t1t2t3t4
      Differences over time within all tibial and all talar ROIs:
      (One-way) univariate ANOVA.
      RightTibial plafond0.1680.5640.0260.2300.325
      Talar dome0.2450.5590.846
      Variance homogenity not given (Levene test < 0.1).
      0.9220.611
      LeftTibial plafond0.1210.0130.0130.0070.295
      Talar dome0.3020.8550.4670.9900.757
      Differences over time of combined tibial and talar ROIs within a slice:
      Independent t test.
      RightLateral0.1180.1530.1770.5270.278
      Central0.3930.013<0.001
      Variance homogenity not given (Levene test < 0.1).
      0.1090.147
      Medial0.5680.2560.0080.0540.013
      All ROIs0.2120.0560.0020.0840.058
      LeftLateral0.2630.9060.9700.4660.590
      Central0.3180.005
      Variance homogenity not given (Levene test < 0.1).
      0.001
      Variance homogenity not given (Levene test < 0.1).
      0.0020.052
      Medial0.0350.0110.0010.0010.015a
      All ROIs0.5960.0210.0030.0020.070
      ROIs: regions of interest (tibial-medial, tibial-central, tibial-lateral, talar-medial, talar-central, talar-lateral).
      Bold fonts show significant differences (P-values).
      (One-way) univariate ANOVA.
      Independent t test.
      Variance homogenity not given (Levene test < 0.1).
      The relative changes of T2*-values compared to baseline during TEFR revealed a T2*-increase between t0 and t2 in all regions. This increase amounted 27.7% (tibial 26.1%, talar 29.5%) one the right and 23.5% (tibial 18.4%, talar 28.8%) on the left in mean. Total standard error of the signal increase ranged from 22.3 to 33.0% on the right and 18.9 to 28.0% on the left (Fig. 3) and was significant for all ROIs throughout TEFR, with exception of the left central tibial ROI (Table III). Subsequently, T2*-signal decreases in all ROIs: mean 10.1% (10.8% tibial, 9.4% talar) on the right and 6.2% (6.0% tibial, 6.4% talar) on the left (Fig. 3). With the exception of the right medial talar ROI–(linear) decrease, all other ROI's showed quadratic decrease. All tibial and talar ROIs, as well as combined ROI's show at least a medium to high effects size of T2*-value decrease (Table III).
      Figure thumbnail gr3
      Fig. 3T2*-mapping of TTJ: Relative changes of T2* relaxation times of single and aggregated segments compared to start (finisher group; nF = 13).
      Table IIIChanges of intrachondral T2* values [ms] in the course of TEFR
      SideROIs and total jointnMauchly-testTest on T2*-signal changes ANOVA
      “Greenhouse-Geisser” correction procedure was used.
      Paired t test on secondary T2*-decreaseTest on quadratic innersubject effects
      P-valueP-valueTest powerP-valueEffect sizeP-valueTest power
      RightTibial-lateral120.157<0.0010.9990.0130.95
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      0.0030.929
      Significant quadratic trend of T2*-signal curve shows high test power.
      Tibial-central0.390<0.0011.0000.0020.85
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      0.0020.949
      Significant quadratic trend of T2*-signal curve shows high test power.
      Tibial-medial0.4540.0080.871<0.0011.41
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      0.0060.864
      Significant quadratic trend of T2*-signal curve shows high test power.
      Talar-lateral0.563<0.0010.9870.0010.93
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      0.0030.923
      Significant quadratic trend of T2*-signal curve shows high test power.
      Talar-central0.860<0.0011.0000.0011.04
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      <0.0011.000
      Significant quadratic trend of T2*-signal curve shows high test power.
      Talar-medial0.650<0.0010.9900.0020.79
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      0.0280.632
      Tibial all0.508<0.0011.000<0.0011.08
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      <0.0010.995
      Significant quadratic trend of T2*-signal curve shows high test power.
      Talar all0.668<0.0011.0000.0030.95
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      <0.0010.993
      Significant quadratic trend of T2*-signal curve shows high test power.
      Total joint0.659<0.0011.0000.0011.02
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      <0.0010.995
      Significant quadratic trend of T2*-signal curve shows high test power.
      LeftTibial-lateral110.1440.0020.9370.0030.86
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      0.0050.895
      Significant quadratic trend of T2*-signal curve shows high test power.
      Tibial-central0.8760.2120.4300.0150.580.2820.178
      Tibial-medial0.528<0.0011.000<0.0010.85
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      <0.0010.996
      Significant quadratic trend of T2*-signal curve shows high test power.
      Talar-lateral0.068<0.0011.0000.0010.85
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      0.0050.897
      Significant quadratic trend of T2*-signal curve shows high test power.
      Talar-central0.0080.004
      “Greenhouse-Geisser” correction procedure was used.
      0.9000.0010.81
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      0.0030.929
      Significant quadratic trend of T2*-signal curve shows high test power.
      Talar-medial0.0020.001
      “Greenhouse-Geisser” correction procedure was used.
      0.9780.0020.60<0.0010.999
      Significant quadratic trend of T2*-signal curve shows high test power.
      Tibial all0.188<0.0010.9930.0040.750.0020.960
      Significant quadratic trend of T2*-signal curve shows high test power.
      Talar all0.017<0.001
      “Greenhouse-Geisser” correction procedure was used.
      0.9960.0030.65<0.0010.997
      Significant quadratic trend of T2*-signal curve shows high test power.
      Total joint0.040<0.001
      “Greenhouse-Geisser” correction procedure was used.
      0.9940.0140.59<0.0010.996
      Significant quadratic trend of T2*-signal curve shows high test power.
      ROIs: regions of interest (tibial-medial, tibial-central, tibial-lateral, talar-medial, talar-central, talar-lateral).
      Bold fonts show significance (P-values).
      “Greenhouse-Geisser” correction procedure was used.
      Significant T2*-relaxation time decrease after initial increase with large effect size (Cohen's d > 0.8) at the end of TEFR related to max. (peak) in the first part of TEFR.
      Significant quadratic trend of T2*-signal curve shows high test power.
      On the right a significant dependency between total distance run, total run time and number of finished stages till time of MRI scanning could be shown for all combined tibial and talar ROI's and all talar ROIS's (Table IV). On the left such dependencies were only detectable for the medial tibial ROI, whereas the all combined tibial and talar ROI's were dependent on number of finished stages till time of MRI also (Table IV). The distance run of the stage at the day of mobile MRI had relevant influence on T2*-value measurements of the right combined tibial and talar ROI's, all talar ROI's and the talar-central ROI (Table IV). On the left side this was only detectable for the tibial-lateral ROI (Table IV). BM and outcome parameters cartilage thickness and T2* showed no significant relationship at any measurement.
      Table IVAnalysis on significant relationships between T2* values [ms] and time of mobile MRI measurements (test of fixed effects – linear mixed models)
      SideROIs and total jointStage no. at time of MRI [n]Total distance run [km]Total run time [hrs]Distance run since last MRI [km]Run time since last MRI [hrs]Stage distance at day of MRI [km]Stage run time at day of MRI [hrs]Time gap between stage finish and MRI [hrs]
      RightTibial-lateral0.3360.3720.3640.2660.9050.2070.4160.681
      Tibial-central0.2860.2700.2070.2730.6690.2850.5710.483
      Tibial-medial0.0290.0280.0390.1420.2950.1660.0940.016
      Talar-lateral0.0440.0420.0600.0830.7520.0720.1920.143
      Talar-central0.0180.0070.0160.0360.2460.0030.0260.061
      Talar-medial0.5500.5520.4500.8900.1630.1530.1670.268
      Tibial all0.0720.0690.0680.1490.8840.1170.1860.118
      Talar all0.0490.0480.0490.1170.2110.0040.0430.060
      Total joint0.0490.0490.0490.1230.4300.0280.0660.071
      LeftTibial-lateral0.0790.0780.1050.0110.1520.0290.0870.490
      Tibial-central0.9570.9700.9740.6140.5590.6600.2970.528
      Tibial-medial0.0120.0120.0090.0170.1500.1500.2390.914
      Talar-lateral0.3140.2970.3470.1900.6700.7070.8840.716
      Talar-central0.1280.1270.1580.1130.1450.1410.0940.051
      Talar-medial0.2060.2170.1770.1930.0670.1710.2820.370
      Tibial all0.0660.0470.0800.0360.1060.1450.3840.268
      Talar all0.0860.0930.0890.1070.1140.1910.2810.162
      Total joint0.0340.1060.0740.0490.0520.1300.2740.115
      ROIs: regions of interest (tibial-medial, tibial-central, tibial-lateral, talar-medial, talar-central, talar-lateral).
      Bold fonts show significance (P-values).
      Subchondral bone marrow abnormalities were identified in two finishers (49 and 60 years). An ill defined area of subchondral posterior tibial bone marrow edema as evidenced by T2-elevation with an Outerbridge grade 1 lesion of the overlying articular cartilage. The osteochondral lesion did not change in size, however, there was a steady increase of the T2* value over time observed. The bone marrow edema showed no relevant changes throughout TEFR [Fig. 4(A) + (B)]. In three subjects (13.6%) a small unilateral subchondral tibial area of bone marrow edema was visible at baseline with no associated cartilage alterations (grade 0) [Fig. 4(C)]. A prominent posterior process of the talus in two subjects remained unchanged during whole TEFR [Fig. 4(D)].
      Figure thumbnail gr4
      Fig. 4Osteochondral lesions (yellow arrow) in the ankle joint. Sagittal TIRM slices through TTJ: A: of a 59-years-old male finisher (1: baseline, 2: stage 32/2176 km run, 3: stage 53/3669 km run, 4: 8 months after TEFR) B: of a 49-years-old male non-finisher (1: baseline, 2: stage 19/1260 km run, 3:stage 31/2131 km run, 4: stage 53/3669 km run) C: of a 46-years-old male finisher (1: baseline, 2: stage 52/3609 km run) and a 68-years-old female non-finisher (3: baseline, 4: stage 15/1003 km run) D: of a 34-years-old male finisher (1: baseline, 2: stage 18/1192 km run, 3:stage 29/1985 km run, 4: stage 54/3763 km run).

      Discussion

      This study provides a unique insight into the response of the TTJ to extreme physical strains of a contiguous 64-day 4486 km ultramarathon. Regarding our first hypothesis, that the normal ankle joint can withstand extreme physical loading, our study revealed three significant observations.
      First, upon morphologic evaluation of the MR images of the ankle no significant pathology was detected in the runners, which had no pre-existing abnormalities. Specifically, the cartilage of the talar dome and the tibial plateau remained intact over the entire time course of the race; there was no evidence of subchondral bone marrow edema or significant joint effusion. None of the runners were eliminated from the study because of subjective or radiographic evidence of ankle pathology. This observation is in keeping with statistics on prevalence of ankle osteoarthritis (OA), which is about 8–10 times less frequent than knee OA
      • Hesper T.
      • Hosalkar H.S.
      • Bittersohl D.
      • Welsch G.H.
      • Krauspe R.
      • Zilkens C.
      • et al.
      T2* mapping for articular cartilage assessment: principles, current applications, and future prospects.
      . While all joints are prone to develop OA over time, differences in joint kinetics, cartilage thickness and mechanical properties may result in different degrees of resilience to degeneration. The relative resilience of the native ankle mortise is attributed to its high degree of congruency and associated stability, the thin articular cartilage of the TTJ was found to be withstanding much higher forces of loading than in the knee
      • Koepp H.
      • Eger W.
      • Muehleman C.
      • Valdellon A.
      • Buckwalter J.A.
      • Kuettner K.E.
      • et al.
      Prevalence of articular cartilage degeneration in the ankle and knee joints of human organ donors.
      • Hesper T.
      • Hosalkar H.S.
      • Bittersohl D.
      • Welsch G.H.
      • Krauspe R.
      • Zilkens C.
      • et al.
      T2* mapping for articular cartilage assessment: principles, current applications, and future prospects.
      .
      Second, the morphometric measurements of mean cartilage thickness revealed no significant change over the course of the race. Measurements of cartilage thickness with MRI are precise and reproducible
      • Trattnig S.
      • Breitenseher M.J.
      • Huber M.
      • Zettl R.
      • Rottmann B.
      • Haller J.
      • et al.
      Determination of cartilage thickness in the ankle joint. An MRT (1.5)-anatomical comparative study.
      • Eckstein F.
      • Charles H.C.
      • Buck R.J.
      • Kraus V.B.
      • Remmers A.E.
      • Hudelmaier M.
      • et al.
      Accuracy and precision of quantitative assessment of cartilage morphology by magnetic resonance imaging at 3.0 T.
      . Al-Ali et al.
      • Al-Ali D.
      • Graichen H.
      • Faber S.
      • Englmeier K.H.
      • Reiser M.
      • Eckstein F.
      Quantitative cartilage imaging of the human hind foot: precision and inter-subject variability.
      found a similar reproducibility and interrater precision of cartilage thickness measurements. Our measurements agree with findings in the literature
      • Trattnig S.
      • Breitenseher M.J.
      • Huber M.
      • Zettl R.
      • Rottmann B.
      • Haller J.
      • et al.
      Determination of cartilage thickness in the ankle joint. An MRT (1.5)-anatomical comparative study.
      and confirm, that the articular cartilage of the ankle is thin, when compared to the knee
      • Hesper T.
      • Hosalkar H.S.
      • Bittersohl D.
      • Welsch G.H.
      • Krauspe R.
      • Zilkens C.
      • et al.
      T2* mapping for articular cartilage assessment: principles, current applications, and future prospects.
      . These findings provide further evidence of specific anatomical and biomechanical properties in the ankle as a primary rolling joint with congruent surfaces making the TTJ less susceptible to osteoarthritic changes and more resilient than other joints
      • Cole A.A.
      • Margulis A.
      • Kuettner K.E.
      Distinguishing ankle and knee articular cartilage.
      • Aurich M.
      • Squires G.R.
      • Reiner A.
      • Mollenhauer J.A.
      • Kuettner K.E.
      • Poole A.R.
      • et al.
      Differential matrix degradation and turnover in early cartilage lesions of human knee and ankle joints.
      . Ankle cartilage is stiffer than cartilage in the knee
      • Treppo S.
      • Koepp H.
      • Quan E.C.
      • Cole A.A.
      • Kuettner K.E.
      • Grodzinsky A.J.
      Comparison of biomechanical and biochemical properties of cartilage from human knee and ankle pairs.
      or hip
      • Kempson G.E.
      Age-related changes in the tensile properties of human articular cartilage: a comparative study between the femoral head of the hip joint and the talus of the ankle joint.
      , which is inversely correlated to cartilage thickness
      • Buettner O.
      • Leumann A.
      • Lehner R.
      • Dell-Kuster S.
      • Rosenthal R.
      • Mueller-Gerbl M.
      • et al.
      Histomorphometric, CT arthrographic, and biomechanical mapping of the human ankle.
      . This is attributed to its special molecular composition with higher glycosaminoglycan (GAG) and less water content
      • Treppo S.
      • Koepp H.
      • Quan E.C.
      • Cole A.A.
      • Kuettner K.E.
      • Grodzinsky A.J.
      Comparison of biomechanical and biochemical properties of cartilage from human knee and ankle pairs.
      when compared to the knee joint.
      Third, while conventional MR based morphological assessment and measurements revealed stability over the entire race, quantitative T2* measurements allowed insight into the adaptive capacities of the most important structure for load bearing and load dissipation, the articular cartilage. It is well established, that mechanical loading of the articular cartilage regulates homeostasis and normal tissues remodeling. The two components that enable the tissue to withstand compressive stress are a liquid and a multicomponent solid collagen and hydrophilic proteoglycan (PG) phase. Water, the most abundant component in conjunction with free mobile cations greatly influence the mechanical behavior of articular cartilage
      • Urban J.P.
      • Maroudas A.
      • Bayliss M.T.
      • Dillon J.
      Swelling pressures of proteoglycans at the concentrations found in cartilaginous tissues.
      . MRI is uniquely suited to probe the interaction between water, collagen and PGs
      • Burstein D.
      • Gray M.
      • Mosher T.
      • Dardzinski B.
      Measures of molecular composition and structure in osteoarthritis.
      , which have an important function in regulating the structural organization of the extracellular matrix and its swelling properties
      • Bader D.L.
      • Salter D.M.
      • Chowdhury T.T.
      Biomechanical influence of cartilage homeostasis in health and disease.
      via the Donnan osmotic equilibrium. T2* measurement is most sensitive to changes in free water content and collagen orientation with an additional component introduced by susceptibility changes otherwise refocused in T2-mapping
      • Koepp H.
      • Eger W.
      • Muehleman C.
      • Valdellon A.
      • Buckwalter J.A.
      • Kuettner K.E.
      • et al.
      Prevalence of articular cartilage degeneration in the ankle and knee joints of human organ donors.
      . Most of the water occupies the interfibrillar space of the extracellular matrix and it is believed that approximately 70% of the water is free to move when loaded by compressive forces
      • Mow V.C.
      • Guo X.E.
      Mechano-electrochemical properties of articular cartilage: their inhomogeneities and anisotropies.
      . This interstitial fluid movement is important in controlling cartilage mechanical behavior
      • Buettner O.
      • Leumann A.
      • Lehner R.
      • Dell-Kuster S.
      • Rosenthal R.
      • Mueller-Gerbl M.
      • et al.
      Histomorphometric, CT arthrographic, and biomechanical mapping of the human ankle.
      . The swelling pressures of articular cartilage
      • Urban J.P.
      • Maroudas A.
      • Bayliss M.T.
      • Dillon J.
      Swelling pressures of proteoglycans at the concentrations found in cartilaginous tissues.
      is, in turn resisted and balanced by tension developed in the collagen network according to Starling's law. Permeability decreased exponentially as function of both increasing compressive strain as well as increasing fluid pressure
      • Mansour J.M.
      • Mow V.C.
      The permeability of articular cartilage under compressive strain and at high pressures.
      ; an immediate adaptive mechanism with loading and unloading, by which decreased permeability allows an overall increase of free water within the articular cartilage as long as the collagen network and PGs can structurally withstand the increasing fluid pressure.

      Immediate T2/T2*-decrease and recovery

      With cyclical loading of cartilage, as it occurs in running, tissue deformation likely produces an increase in anisotropy of superficial collagen fibers and a concomitant decrease in free cartilage water which leads to lower T2-values on the post-exercise studies. The results of Mosher et al.
      • Mosher T.J.
      • Liu Y.
      • Torok C.M.
      Functional cartilage MRI T2 mapping: evaluating the effect of age and training on knee cartilage response to running.
      and other investigators indicated that this immediate T2-response to running or other physical load is not influenced by age or level of physical activity. During unloading of chondrocytes complete recovery of all structural deformation was observed after 30 min
      • Kääb M.J.
      • Richards R.G.
      • Ito K.
      • ap Gwynn I.
      • Nötzli H.P.
      Deformation of chondrocytes in articular cartilage under compressive load: a morphological study.
      . Volumetric MRI-analyses of the ankle cartilage showed significant initial talar cartilage volume reduction, which was restored within 30 min
      • Van Ginckel A.
      • Roosen P.
      • Almqvist K.F.
      • Verstraete K.
      • Witvrouw E.
      Effects of in vivo exercise on ankle cartilage deformation and recovery in healthy volunteers: an experimental study.
      . In our investigation such immediate T2-decreases could not be detected since T2* data were acquired at least 45 min after unloading and first data points were obtained after a mean distance run of 1000 km when the effects of long time running masked the initial T2-response.

      Long distance T2/T2*-increase

      After 1000 km of running our data revealed a significant mean T2*-elevation up to 23.5–27.7% in all areas of the ankle cartilage (Fig. 3) with the exception of the left tibial central ROI (Table III), where the baseline T2*-value was already high. Forty-eight hours after a marathon Luke et al.
      • Luke A.C.
      • Stehling C.
      • Stahl R.
      • Li X.
      • Kay T.
      • Takamoto S.
      • et al.
      High-field magnetic resonance imaging assessment of articular cartilage before and after marathon running: does long-distance running lead to cartilage damage?.
      found not only significantly higher T2-values but also increased T1ρ values in all articular cartilage of the knee (P < 0.01) except the lateral compartment, postulating a relative decrease of GAGs leading to more free water in the matrix
      • Regatte R.R.
      • Akella S.V.
      • Wheaton A.J.
      • Lech G.
      • Borthakur A.
      • Kneeland J.B.
      • et al.
      3D-T1rho-relaxation mapping of articular cartilage: in vivo assessment of early degenerative changes in symptomatic osteoarthritic subjects.
      . Increased T2 can be seen after marathon running as well as in early OA
      • Golditz T.
      • Steib S.
      • Pfeifer K.
      • Uder M.
      • Gelse K.
      • Janka R.
      • et al.
      Functional ankle instability as a risk factor for osteoarthritis: using T2-mapping to analyze early cartilage degeneration in the ankle joint of young athletes.
      , for different reasons. If the overall environment of the chondrocytes is impaired, the excessive fluid pressure will lead to structural changes of cartilage fibrillation, PG loss, and loss of collagen network organization which has been shown in a long-term (1-year) program of running exercise (up to 40 km/day) in the knee and humeral head cartilage of young dogs
      • Arokoski J.
      • Kiviranta I.
      • Jurvelin J.
      • Tammi M.
      • Helminen H.J.
      Long-distance running causes site-dependent decrease of cartilage glycosaminoglycan content in the knee joints of beagle dogs.
      • Kiviranta I.
      • Tammi M.
      • Jurvelin J.
      • Säämänen A.M.
      • Helminen H.J.
      Moderate running exercise augments glycosaminoglycans and thickness of articular cartilage in the knee joint of young beagle dogs.
      . Since no structural and morphometric changes were observed, the slow, but steady increase of T2 in our data over the course of the first 2000–2500 km, could be explained by changes in permeability and respective ion flow reaching a new equilibrium after every run. As increasing hydrostatic pressure up-regulates PG and type II collagen mRNA expression
      • Toyoda T.
      • Seedhom B.B.
      • Kirkham J.
      • Bonass W.A.
      Upregulation of aggrecan and type II collagen mRNA expression in bovine chondrocytes by the application of hydrostatic pressure.
      de-novo syntheses of PGs will be initiated.

      Ultra-long distance T2/T2*-decrease (“recovery”)

      MRI has not yet been utilized before to probe human cartilage properties beyond a single marathon. A significant T2-increase within the knee cartilage was observed 48 h after a single marathon, which recovered to baseline levels after 3 months
      • Luke A.C.
      • Stehling C.
      • Stahl R.
      • Li X.
      • Kay T.
      • Takamoto S.
      • et al.
      High-field magnetic resonance imaging assessment of articular cartilage before and after marathon running: does long-distance running lead to cartilage damage?.
      . Surprisingly, when ultra-marathon runners in our study continued beyond 2000–2500 km a significant T2*-decrease of about 6–10% occurred, which approximates one third of the primary signal increase. Animal experiments have shown that load bearing exercises minimize the development of OA, increases PG content and cartilage thickness in rodent models
      • Galois L.
      • Etienne S.
      • Grossin L.
      • Cournil C.
      • Pinzano A.
      • Netter P.
      • et al.
      Moderate-impact exercise is associated with decreased severity of experimental osteoarthritis in rats.
      and dogs
      • Kiviranta I.
      • Tammi M.
      • Jurvelin J.
      • Säämänen A.M.
      • Helminen H.J.
      Moderate running exercise augments glycosaminoglycans and thickness of articular cartilage in the knee joint of young beagle dogs.
      . The positive effect of joint loading on chondrocyte function with increased PG and collagen synthesis was found in multiple studies
      • Grodzinsky A.J.
      • Levenston M.E.
      • Jin M.
      • Frank E.H.
      Cartilage tissue remodeling in response to mechanical forces.
      . Hydrostatic pressure up-regulates PG and type II collagen mRNA expression in bovine chondrocytes
      • Toyoda T.
      • Seedhom B.B.
      • Kirkham J.
      • Bonass W.A.
      Upregulation of aggrecan and type II collagen mRNA expression in bovine chondrocytes by the application of hydrostatic pressure.
      .
      After moderate running (1 hr/4 km per day more than 15 weeks) an increase of compression stiffness and GAG-content was observed in young canine patellae and lateral femoral condyles
      • Säämämen A.M.
      • Kiviranta I.
      • Jurvelin J.
      • Helminen H.J.
      • Tammi M.
      Proteoglycan and collagen alterations in canine knee articular cartilage following 20 km daily running exercise for 15 weeks.
      . In several other animal studies it has been shown that moderate exercise may protect against cartilage degeneration
      • Galois L.
      • Etienne S.
      • Grossin L.
      • Cournil C.
      • Pinzano A.
      • Netter P.
      • et al.
      Moderate-impact exercise is associated with decreased severity of experimental osteoarthritis in rats.
      • Otterness I.G.
      • Eskra J.D.
      • Bliven M.L.
      • Shay A.K.
      • Pelletier J.P.
      • Milici A.J.
      Exercise protects against articular cartilage degeneration in the hamster.
      . In a comparison study between lifelong trained canines (4 km per day for more than 10 years) with normally active canines Newton et al.
      • Newton P.M.
      • Mow V.C.
      • Gardner T.R.
      • Buckwalter J.A.
      • Albright J.P.
      Winner of the 1996 Cabaud Award. The effect of lifelong exercise on canine articular cartilage.
      observed no differences of the biochemical and mechanical qualities of the joint cartilage between both groups, the experiments of Lapvetelainen et al.
      • Lapvetelainen T.
      • Hyttinen M.
      • Lindblom J.
      • Langsio T.K.
      • Sironen R.
      • Li S.W.
      • et al.
      More knee joint osteoarthritis (OA) in mice after inactivation of one allele of type II procollagen gene but less OA after lifelong voluntary wheel running exercise.
      with long distance running found less knee OA in running knockout mice than sedentary mice. They demonstrated that ultra-long physical activity does not predispose normal mice to OA
      • Lapvetelainen T.
      • Hyttinen M.
      • Lindblom J.
      • Langsio T.K.
      • Sironen R.
      • Li S.W.
      • et al.
      More knee joint osteoarthritis (OA) in mice after inactivation of one allele of type II procollagen gene but less OA after lifelong voluntary wheel running exercise.
      .
      In the human knee cartilage Roos and Dahlberg
      • Roos E.M.
      • Dahlberg L.
      Positive effects of moderate exercise on glycosaminoglycan content in knee cartilage: a four-month, randomized, controlled trial in patients at risk of osteoarthritis.
      measured an increase of GAG in the weight-bearing posterior medial femoral condyle following moderate exercise (1 hr exercise, 3 times weekly for 4 months). Tiderius et al.
      • Tiderius C.J.
      • Svensson J.
      • Leander P.
      • Ola T.
      • Dahlberg L.
      dGEMRIC (delayed gadolinium-enhanced MRI of cartilage) indicates adaptive capacity of human knee cartilage.
      found in a cross-sectional study with elite runners and untrained volunteers, that human knee cartilage adapts to exercise by increasing the GAG-content using dGEMRIC. The increased concentration of PGs has been shown to impede hydraulic fluid flow
      • Mow V.C.
      • Holmes M.H.
      • Lai W.M.
      Fluid transport and mechanical properties of articular cartilage: a review.
      . This mechanism shields the collagen-PG matrix from high stresses.

      Ultra-long distance T2/T2*-increase without recovery

      With respect to the second question, how much exercise is too much, our observations are significant, when considered in the context of the existing data. Development of ankle OA is intimately related to instability and incongruity caused by trauma
      • Fitzpatrick D.C.
      • Otto J.K.
      • McKinley T.O.
      • Marsh J.L.
      • Brown T.D.
      Kinematic and contact stress analysis of posterior malleolus fractures of the ankle.
      . The impact of traumatic injuries to the tibial plafond, fibula fractures and mainly ankle sprains involving the anterior talofibular ligament (ATFL) or high ankle ligaments leading to OA in up to 70–78% of the patient's
      • Saltzman C.L.
      • Salamon M.L.
      • Blanchard G.M.
      • Huff T.
      • Hayes A.
      • Buckwalter J.A.
      • et al.
      Epidemiology of ankle arthritis: report of a consecutive series of 639 patients from a tertiary orthopaedic center.
      • van Dijk C.N.
      • Reilingh M.L.
      • Zengerink M.
      • van Bergen C.J.
      The natural history of osteochondral lesions in the ankle.
      is established. T2-increase could be detected after ATFL injury in trochlear cartilage
      • Lee S.
      • Yoon Y.C.
      • Kim J.H.
      T2 mapping of the articular cartilage in the ankle: correlation to the status of anterior talofibular ligament.
      . However, overall ankle injuries are rare in ultramarathon runners
      • Bishop G.W.
      • Fallon K.E.
      Musculoskeletal injuries in a six-day track race: ultramarathoner's ankle.
      . Within the cohort of runners in this study, there were only two preexisting low grade osteochondral lesions, which on conventional MR imaging did not change appearance, however revealed a steady increase in T2* values, likely evidence of increased swelling, lacking the recovery phase observed in the ankles without pre-existing injuries.

      Limitations

      There are major limitations in our study, the small number of runners, the lack of complementary quantitative cartilage mapping sequences, which could not be applied due to severe time and study constrains. Specifically, the selective quantitative information on GAG content is missing for more definitive interpretation. However, negative charge density measurements using dGEMRIC, Na or GAG-CEST have their known significant limitations. T1rho (adiabatic or spin-lock) allow detection of PG contributions, but are not selective
      • Regatte R.R.
      • Akella S.V.
      • Wheaton A.J.
      • Lech G.
      • Borthakur A.
      • Kneeland J.B.
      • et al.
      3D-T1rho-relaxation mapping of articular cartilage: in vivo assessment of early degenerative changes in symptomatic osteoarthritic subjects.
      . Depth dependence of T2* could not be analyzed in this study due to the thin ankle cartilage. The stage distance run (in km) at the day of MRI measurements influenced T2*-values in some ROIs. It must be assumed that the immediate T2*-decrease based on an increase in collagen fiber anisotropy and the concomitant free cartilage water
      • Mosher T.J.
      • Liu Y.
      • Torok C.M.
      Functional cartilage MRI T2 mapping: evaluating the effect of age and training on knee cartilage response to running.
      • Souza R.B.
      • Stehling C.
      • Wyman B.T.
      • Hellio Le Graverand M.P.
      • Li X.
      • Link T.M.
      • et al.
      The effects of acute loading on T1rho and T2 relaxation times of tibiofemoral articular cartilage.
      exchange were present in parallel to the mechanism of long distance T2-increase. This “masking” effect on the slower T2*-increase is difficult to estimate quantitatively. The time elapsed between stage finish and MRI had no correlation with the data. Other confounders, which were related to performance, had only limited and inconsistent influence on the results (Table IV), they were not adjusted to confounders. T2* data were adjusted to age, regional variation, injuries as confounding factors. This study is subject to a possible selection bias regarding experience in ultra-running before TEFR, however, the influence on results is difficult to assess. Another bias regarding generalizability of our results may be the fact, that non-finisher were 8 years older than finisher.
      In conclusion, there are at least three distinct time domains, in the complex adaptation to loading of the cartilage; immediately during and after mechanical loading free water will be trapped or released from the articular cartilage via fast permeability changes. As no structural and morphometric changes were observed, the slow, but steady increase of T2 in our data over the course of the first 2500 km, could be explained by changes in permeability and respective ion flow reaching a new equilibrium after every run. As increasing hydrostatic pressure up-regulates PG and type II collagen mRNA expression
      • Toyoda T.
      • Seedhom B.B.
      • Kirkham J.
      • Bonass W.A.
      Upregulation of aggrecan and type II collagen mRNA expression in bovine chondrocytes by the application of hydrostatic pressure.
      subsequent higher PG content shifts the free fluid back into the PG and collagen bound compartment, the ultra-slow adaptive regime, where T2* will slowly decrease and recover. In cartilage regions with underlying pre-existing osteochondral abnormalities the recovery phase could not be observed in our study, to the contrary, T2* continued to increase over the entire race. As running was an important means for survival from an evolutionary point of view, our findings indicate that the human ankle indeed appears to be resilient to ultra-long distance running, unless a pre-existing injury exists.

      Contributions

      All authors of this manuscript had substantial contribution to conception and design or acquisition, analysis and interpretation of data; all revised it critically for important intellectual content and did final approval of the version to be published:
      Schütz UHW conceived the study (conception and design), implemented the project (including administrative, technical and logistical support), obtained the funding, participated in the data collection (MRI measurements), data assembly and data evaluation with statistical analysis and drafted the article. Ellermann J gave critical revision of the article for important intellectual content and did final approval of the article. Schoss D participated in the data evaluation. Wiedelbach H participated mainly in the data collection (MRI measurements) and in the implementation of the project. Beer M gave critical revision of the article for important intellectual content and did final approval of the article. All authors read and approved the final manuscript. Billich C participated in the implementation of the project (including administrative, technical and logistical support), in the data collection (MRI measurements) and in the implementation of the project (including administrative, technical and logistical support).
      Schütz UHW ([email protected]) takes responsibility for the integrity of the work as a whole, from inception to finished article.

      Role of funding source

      This work is supported in part by the German Research Association (DFG: “Deutsche Forschungsgemeinschaft”), under Grants SCHU 2514/1-1 and SCHU 2514/1-2. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. No additional external funding was received for this study.

      Competing interests

      The authors declare that they have no financial or non-financial competing interests. There are no financial or non-financial competing interests of other people or organizations influencing our interpretation of data or presentation of information.

      Acknowledgments

      We would like to thank all the athletes of TEFR who took part at this project. Considering their immense physical and mental stresses they showed an extraordinary compliance on every day of the race.

      References

        • Booth F.W.
        • Roberts C.K.
        • Laye M.J.
        Lack of exercise is a major cause of chronic diseases.
        Compr Physiol. 2012; 2: 1143-1211
        • Schueller-Weidekamm C.
        • Schueller G.
        • Uffmann M.
        • Bader T.
        Incidence of chronic knee lesions in long-distance runners based on training level: findings at MRI.
        Eur J Radiol. 2006; 58: 286-293
        • Krampla W.W.
        • Newrkla S.P.
        • Kroener A.H.
        • Hruby W.F.
        Changes on magnetic resonance tomography in the knee joints of marathon runners: a 10-year longitudinal study.
        Skelet Radiol. 2008; 37: 619-626
        • Schütz U.H.
        • Schmidt-Trucksäss A.
        • Knechtle B.
        • Machann J.
        • Wiedelbach H.
        • Ehrhardt M.
        • et al.
        The TransEurope Footrace Project: longitudinal data acquisition in a cluster randomized mobile MRI observational cohort study on 44 endurance runners at a 64-stage 4,486 km transcontinental ultramarathon.
        BMC Med. 2012; 10: 78
        • El-Khoury G.Y.
        • Alliman K.J.
        • Lundberg H.J.
        • Rudert M.J.
        • Brown T.D.
        • Saltzman C.L.
        Cartilage thickness in cadaveric ankles: measurement with double-contrast multi-detector row CT arthrography versus MR imaging.
        Radiology. 2004; 233: 768-773
        • Trattnig S.
        • Breitenseher M.J.
        • Huber M.
        • Zettl R.
        • Rottmann B.
        • Haller J.
        • et al.
        Determination of cartilage thickness in the ankle joint. An MRT (1.5)-anatomical comparative study.
        RöFo: Fortschritte auf dem Gebiete der Röntgenstahlen und Nuklearmedizin. 1997; 166: 303-306
        • Koepp H.
        • Eger W.
        • Muehleman C.
        • Valdellon A.
        • Buckwalter J.A.
        • Kuettner K.E.
        • et al.
        Prevalence of articular cartilage degeneration in the ankle and knee joints of human organ donors.
        J Orthop Sci. 1999; 4: 407-412
        • Hesper T.
        • Hosalkar H.S.
        • Bittersohl D.
        • Welsch G.H.
        • Krauspe R.
        • Zilkens C.
        • et al.
        T2* mapping for articular cartilage assessment: principles, current applications, and future prospects.
        Skelet Radiol. 2014; 43: 1429-1445
        • Cole A.A.
        • Margulis A.
        • Kuettner K.E.
        Distinguishing ankle and knee articular cartilage.
        Foot Ankle Clin. 2003; 8: 305-316
        • Shellock F.G.
        • Spinazzi A.
        MRI safety update 2008: part 2, screening patients for MRI.
        AJR Am J Roentgenol. 2008; 191: 1140-1149
        • Apprich S.
        • Welsch G.H.
        • Mamisch T.C.
        • Szomolanyi P.
        • Mayerhoefer M.
        • Pinker K.
        • et al.
        Detection of degenerative cartilage disease: comparison of high-resolution morphological MR and quantitative T2 mapping at 3.0 Tesla.
        Osteoarthritis Cartilage. 2010; 18: 1211-1217
        • Mosher T.J.
        MRI of osteochondral injuries of the knee and ankle in the athlete.
        Clin Sports Med. 2006; 25: 843-866
        • Bland J.M.
        • Altman D.G.
        Statistical methods for assessing agreement between two methods of clinical measurement.
        Lancet. 1986; 1: 307-310
        • Jepsen J.R.
        • Laursen L.H.
        • Hagert C.G.
        • Kreiner S.
        • Larsen A.I.
        Diagnostic accuracy of the neurological upper limb examination I: inter-rater reproducibility of selected findings and patterns.
        BMC Neurol. 2006; 6: 8
        • Cohen J.
        Statistical Power Analysis for the Behavioral Sciences.
        2nd edn. Lawrence Erlbaum Associates, , Hillsdale1988
        • Eckstein F.
        • Charles H.C.
        • Buck R.J.
        • Kraus V.B.
        • Remmers A.E.
        • Hudelmaier M.
        • et al.
        Accuracy and precision of quantitative assessment of cartilage morphology by magnetic resonance imaging at 3.0 T.
        Arthritis Rheum. 2005; 52: 3132-3136
        • Al-Ali D.
        • Graichen H.
        • Faber S.
        • Englmeier K.H.
        • Reiser M.
        • Eckstein F.
        Quantitative cartilage imaging of the human hind foot: precision and inter-subject variability.
        J Orthop Res. 2002; 20: 249-256
        • Aurich M.
        • Squires G.R.
        • Reiner A.
        • Mollenhauer J.A.
        • Kuettner K.E.
        • Poole A.R.
        • et al.
        Differential matrix degradation and turnover in early cartilage lesions of human knee and ankle joints.
        Arthritis Rheum. 2005; 52: 112-119
        • Treppo S.
        • Koepp H.
        • Quan E.C.
        • Cole A.A.
        • Kuettner K.E.
        • Grodzinsky A.J.
        Comparison of biomechanical and biochemical properties of cartilage from human knee and ankle pairs.
        J Orthop Res. 2000; 18: 739-748
        • Kempson G.E.
        Age-related changes in the tensile properties of human articular cartilage: a comparative study between the femoral head of the hip joint and the talus of the ankle joint.
        Biochim Biophys Acta. 1991; 1075: 223-230
        • Buettner O.
        • Leumann A.
        • Lehner R.
        • Dell-Kuster S.
        • Rosenthal R.
        • Mueller-Gerbl M.
        • et al.
        Histomorphometric, CT arthrographic, and biomechanical mapping of the human ankle.
        Foot Ankle Int. 2013; 34: 1025-1034
        • Urban J.P.
        • Maroudas A.
        • Bayliss M.T.
        • Dillon J.
        Swelling pressures of proteoglycans at the concentrations found in cartilaginous tissues.
        Biorheology. 1979; 16: 447-464
        • Burstein D.
        • Gray M.
        • Mosher T.
        • Dardzinski B.
        Measures of molecular composition and structure in osteoarthritis.
        Radiol Clin North Am. 2009; 47: 675-686
        • Bader D.L.
        • Salter D.M.
        • Chowdhury T.T.
        Biomechanical influence of cartilage homeostasis in health and disease.
        Arthritis. 2011; 2011: 979032
        • Mow V.C.
        • Guo X.E.
        Mechano-electrochemical properties of articular cartilage: their inhomogeneities and anisotropies.
        Annu Rev Biomed Eng. 2002; 4: 175-209
        • Mansour J.M.
        • Mow V.C.
        The permeability of articular cartilage under compressive strain and at high pressures.
        J Bone Joint Surg Am. 1976; 58: 509-516
        • Mosher T.J.
        • Liu Y.
        • Torok C.M.
        Functional cartilage MRI T2 mapping: evaluating the effect of age and training on knee cartilage response to running.
        Osteoarthritis Cartilage. 2010; 18: 358-364
        • Kääb M.J.
        • Richards R.G.
        • Ito K.
        • ap Gwynn I.
        • Nötzli H.P.
        Deformation of chondrocytes in articular cartilage under compressive load: a morphological study.
        Cells Tissues Organs. 2003; 175: 133-139
        • Van Ginckel A.
        • Roosen P.
        • Almqvist K.F.
        • Verstraete K.
        • Witvrouw E.
        Effects of in vivo exercise on ankle cartilage deformation and recovery in healthy volunteers: an experimental study.
        Osteoarthritis Cartilage. 2011; 19: 1123-1131
        • Luke A.C.
        • Stehling C.
        • Stahl R.
        • Li X.
        • Kay T.
        • Takamoto S.
        • et al.
        High-field magnetic resonance imaging assessment of articular cartilage before and after marathon running: does long-distance running lead to cartilage damage?.
        Am J Sports Med. 2010; 38: 2273-2280
        • Regatte R.R.
        • Akella S.V.
        • Wheaton A.J.
        • Lech G.
        • Borthakur A.
        • Kneeland J.B.
        • et al.
        3D-T1rho-relaxation mapping of articular cartilage: in vivo assessment of early degenerative changes in symptomatic osteoarthritic subjects.
        Acad Radiol. 2004; 11: 741-749
        • Golditz T.
        • Steib S.
        • Pfeifer K.
        • Uder M.
        • Gelse K.
        • Janka R.
        • et al.
        Functional ankle instability as a risk factor for osteoarthritis: using T2-mapping to analyze early cartilage degeneration in the ankle joint of young athletes.
        Osteoarthritis Cartilage. May 2014; ([Epub ahead of print])
        • Arokoski J.
        • Kiviranta I.
        • Jurvelin J.
        • Tammi M.
        • Helminen H.J.
        Long-distance running causes site-dependent decrease of cartilage glycosaminoglycan content in the knee joints of beagle dogs.
        Arthritis Rheum. 1993; 36: 1451-1459
        • Kiviranta I.
        • Tammi M.
        • Jurvelin J.
        • Säämänen A.M.
        • Helminen H.J.
        Moderate running exercise augments glycosaminoglycans and thickness of articular cartilage in the knee joint of young beagle dogs.
        J Orthop Res. 1988; 6: 188-195
        • Toyoda T.
        • Seedhom B.B.
        • Kirkham J.
        • Bonass W.A.
        Upregulation of aggrecan and type II collagen mRNA expression in bovine chondrocytes by the application of hydrostatic pressure.
        Biorheology. 2003; 40: 79-85
        • Galois L.
        • Etienne S.
        • Grossin L.
        • Cournil C.
        • Pinzano A.
        • Netter P.
        • et al.
        Moderate-impact exercise is associated with decreased severity of experimental osteoarthritis in rats.
        Rheumatol (Oxford). 2003; 42: 692-693
        • Grodzinsky A.J.
        • Levenston M.E.
        • Jin M.
        • Frank E.H.
        Cartilage tissue remodeling in response to mechanical forces.
        Annu Rev Biomed Eng. 2000; 2: 691-713
        • Säämämen A.M.
        • Kiviranta I.
        • Jurvelin J.
        • Helminen H.J.
        • Tammi M.
        Proteoglycan and collagen alterations in canine knee articular cartilage following 20 km daily running exercise for 15 weeks.
        Connect Tissue Res. 1994; 30: 191-201
        • Otterness I.G.
        • Eskra J.D.
        • Bliven M.L.
        • Shay A.K.
        • Pelletier J.P.
        • Milici A.J.
        Exercise protects against articular cartilage degeneration in the hamster.
        Arthritis Rheum. 1998; 41: 2068-2076
        • Newton P.M.
        • Mow V.C.
        • Gardner T.R.
        • Buckwalter J.A.
        • Albright J.P.
        Winner of the 1996 Cabaud Award. The effect of lifelong exercise on canine articular cartilage.
        Am J Sports Med. 1997; 25: 282-287
        • Lapvetelainen T.
        • Hyttinen M.
        • Lindblom J.
        • Langsio T.K.
        • Sironen R.
        • Li S.W.
        • et al.
        More knee joint osteoarthritis (OA) in mice after inactivation of one allele of type II procollagen gene but less OA after lifelong voluntary wheel running exercise.
        Osteoarthritis Cartilage. 2001; 9: 152-160
        • Roos E.M.
        • Dahlberg L.
        Positive effects of moderate exercise on glycosaminoglycan content in knee cartilage: a four-month, randomized, controlled trial in patients at risk of osteoarthritis.
        Arthritis Rheum. 2005; 52: 3507-3514
        • Tiderius C.J.
        • Svensson J.
        • Leander P.
        • Ola T.
        • Dahlberg L.
        dGEMRIC (delayed gadolinium-enhanced MRI of cartilage) indicates adaptive capacity of human knee cartilage.
        Magn Reson Med. 2004; 51: 286-290
        • Mow V.C.
        • Holmes M.H.
        • Lai W.M.
        Fluid transport and mechanical properties of articular cartilage: a review.
        J Biomech. 1984; 17: 377-394
        • Fitzpatrick D.C.
        • Otto J.K.
        • McKinley T.O.
        • Marsh J.L.
        • Brown T.D.
        Kinematic and contact stress analysis of posterior malleolus fractures of the ankle.
        J Orthop Trauma. 2004; 18: 271-278
        • Saltzman C.L.
        • Salamon M.L.
        • Blanchard G.M.
        • Huff T.
        • Hayes A.
        • Buckwalter J.A.
        • et al.
        Epidemiology of ankle arthritis: report of a consecutive series of 639 patients from a tertiary orthopaedic center.
        Iowa Orthop J. 2005; 25: 44-46
        • van Dijk C.N.
        • Reilingh M.L.
        • Zengerink M.
        • van Bergen C.J.
        The natural history of osteochondral lesions in the ankle.
        Instr Course Lect. 2010; 59: 375-386
        • Lee S.
        • Yoon Y.C.
        • Kim J.H.
        T2 mapping of the articular cartilage in the ankle: correlation to the status of anterior talofibular ligament.
        Clin Radiol. 2013; 68: e355-e361
        • Bishop G.W.
        • Fallon K.E.
        Musculoskeletal injuries in a six-day track race: ultramarathoner's ankle.
        Clin J Sport Med. 1999; 9: 216-220
        • Souza R.B.
        • Stehling C.
        • Wyman B.T.
        • Hellio Le Graverand M.P.
        • Li X.
        • Link T.M.
        • et al.
        The effects of acute loading on T1rho and T2 relaxation times of tibiofemoral articular cartilage.
        Osteoarthritis Cartilage. 2010; 18: 1557-1563