Skip to main content

Alterations in articular cartilage T2 star relaxation time following mechanical disorders: in vivo canine supraspinatus tendon resection models

Abstract

Background

The role of altered joint mechanics on cartilage degeneration in in vivo models has not been studied successfully due to a lack of pre-injury information. We aimed 1) to develop an accurate in vivo canine model to measure the changes in joint loading and T2 star (T2*) relaxation time before and after unilateral supraspinatus tendon resections, and 2) to find the relationship between regional variations in articular cartilage loading patterns and T2* relaxation time distributions.

Methods

Rigid markers were implanted in the scapula and humerus of tested dogs. The movement of the shoulder bones were measured by a motion tracking system during normal gaits. In vivo cartilage contact strain was measured by aligning 3D shoulder models with the motion tracking data. Articular cartilage T2* relaxation times were measured by quantitative MRI scans. Articular cartilage contact strain and T2* relaxation time were compared in the shoulders before and 3 months after the supraspinatus tendon resections.

Results

Excellent accuracy and reproducibility were found in our in vivo contact strain measurements with less than 1% errors. Changes in articular cartilage contact strain exhibited similar patterns with the changes in the T2* relaxation time after resection surgeries. Regional changes in the articular cartilage T2* relaxation time exhibited positive correlations with regional contact strain variations 3 months after the supraspinatus resection surgeries.

Conclusion

This is the first study to measure in vivo articular cartilage contact strains with high accuracy and reproducibility. Positive correlations between contact strain and T2* relaxation time suggest that the articular cartilage extracellular matrix may responds to mechanical changes in local areas.

Peer Review reports

Background

Articular cartilage is known to be modulated by mechanical loading. Changes in loading conditions have been reported to alter the structural and compositional characteristics of articular cartilage, and may initiate osteoarthritis (OA) development [1,2,3]. Previous studies have shown that the biochemical composition of articular cartilage varies across different regions and is related to regional loading history [4,5,6]. Weight-bearing regions of articular cartilage exhibit greater proteoglycan content compared to non-weight-bearing articular cartilage, whereas collagen content has been found to be greater in non-weight-bearing areas [5,6,7,8]. Because proteoglycan and collagen networks are crucial for maintaining the mechanical strength of articular cartilage, the conditions of proteoglycan content and collagen networks can provide valuable information for estimating the degree of articular cartilage degeneration [9,10,11].

Traditional measurements of articular cartilage compositional properties are based on destructive biochemical and histological procedures that are not applicable for clinical diagnostic purposes. Recently, T1 and T2 relaxation properties of quantitative magnetic resonance imaging (qMRI) have been proposed for estimating articular cartilage biochemical compositions. Because qMRI utilizes the magnetic relaxation characteristics of water components and macromolecules, magnetic relaxation parameters (T1, T1ρ, and T2) are believed to reflect the tissue hydration, collagen, and proteoglycan compositions of articular cartilage [12,13,14,15,16]. Among these relaxation values, T2 relaxation time has been widely used to examine the degeneration of articular cartilage tissues [12, 17,18,19]. Prolonged T2 relaxation time is believed to be a sign of degenerative articular cartilage, which indicates damage to collagen networks and increased water content in articular cartilage extracellular matrices [18,19,20].

Many studies have investigated the relationship between T2 relaxation time and loading patterns [21,22,23]. Significant unloading of knee joints resulted in increased T2 relaxation time, but T2 relaxation time returned to the normal level when the knee resumed natural weight-bearing activities [21]. Mild exercises with low-level shear and compressive loadings have been reported to effectively reduce the T2 relaxation time of articular cartilage in the early stages of OA [22, 23]. These studies indicate that mild mechanical loadings are important for maintaining healthy articular cartilage whereas a large magnitude of mechanical loading or unloading is detrimental to articular cartilage. Articular cartilage T2 relaxation time measurements have also been used to estimate articular cartilage conditions following joint injuries in various clinical studies. The cartilage T2 relaxation times in knees with anterior cruciate ligament (ACL) tears or reconstructed knees were found to be greater than the times in intact knees [24]. Because ACL injuries significantly alter joint biomechanics [25, 26], changes in T2 relaxation time may be the consequence of mechanically induced articular cartilage degradation. Recently, T2 star (T2*) relaxation time has been widely examined as an alternative method of T2 relaxation time measurements because T2* relaxation time have been reported to positively correlate with T2 relaxation time values on the estimation of articular cartilage degeneration [27, 28]. Additionally, the short acquisition time, high signal-to-noise ratio (SNR), and high out-of-plane resolutions of T2* imaging procedures should be significant advantages in clinical imaging situations [16, 27].

In previous clinical studies, mechanically induced articular cartilage degeneration and T2 relaxation time has been investigated by comparing stress/strain patterns and T2 relaxation time distributions between the joints with injuries and the contralateral intact joints, while the cartilage conditions prior to injuries were not available in most cases [29, 30]. However, it is important to study the changes in joint mechanics and T2 relaxation time before and after injuries in a joint to completely understand the role of mechanical alterations on articular cartilage degenerations because articular cartilage biochemical and biomechanical properties might be different between contralateral joints even in the same individual [31]. In this study, we aimed to develop an accurate in vivo canine model to measure the changes in joint loading and T2* relaxation time before and after supraspinatus tendon tears. We measured regional variations in articular cartilage loading patterns and T2* relaxation time distributions after supraspinatus resections. We hypothesized that supraspinatus tendon tears would alter the loading patterns and T2* relaxation time distributions in the shoulder articular cartilages. We also hypothesized that spatial variations in joint loading and T2* relaxation time would be correlated. 

Methods

Ethical statement and study design

This study was approved by the Animal Care and Use Committee (approval number: BA1507–180/047–01) and carried out in accordance with the Guide for the Care and Use of Laboratory Animals of the Seoul National University College of Medicine (Republic of Korea). Total six shoulders from three adolescent mongrel dogs (approximately one year old and weighed 20 kg) were examined. Dogs were obtained from Kukje Laboratory Animal Center (Republic of Korea). After the completion of the study, the animals were euthanized with intraventricular administration of potassium chloride (2 mmol/kg) under deep anesthesia by following AVMA (American Veterinary Medical Association) guidelines.

Experimental procedure and experimental animals

Gait analysis

In this study, we measured articular cartilage deformation patterns by using kinematics-based 3 dimensional (3D) shoulder models. We decided to use a marker-based motion tracking system (eight infrared cameras with a sampling frequency of 120 frames/sec, Eagle Digital real-time system, MotionAnalysis Co., USA). The large tracking areas of motion tracking systems are suitable for normal walking activities in in vivo canine models.

We decided to implant custom markers into shoulder bones rigidly to minimize skin artifacts in the motion tracking data. A custom marker was designed with three conventional reflective balls (diameter of 8 mm) (Fig. 1a). Two custom markers were installed in the humerus and scapula using orthopedic stainless surgical threaded pins. Scapular pins were inserted in two locations along the scapular spine and humeral pins were inserted into the distal and proximal ends of the humerus along the lateral side of the diaphysis (Fig. 1b). After four to 5 days of recovery following the marker insertion surgeries, we recorded the movements of the shoulder bones by tracking the ball markers as the dogs walked freely within our motion laboratory. Only one shoulder was tested in a single motion analysis examination to minimize the loosening of surgical pins. We found that the pins in shoulder bones were easily loosen when the tested dog laid the pin-inserted shoulder on the floor. Thus, the contralateral shoulder remained intact to allow the tested dog to use the intact shoulder in resting positions. After 2 weeks of recovery time, gait analysis of the contralateral shoulder was followed thorough the identical marker installation and motion tracking procedures.

Fig. 1
figure 1

A custom marker and its installation to shoulder bones. a A custom marker was designed with three conventional reflective balls. b Two custom markers were inserted along the scapular spine and other two markers were inserted into the distal and proximal ends of the humerus along the lateral side of the diaphysis

3D shoulder models

The 3D shoulder models were created from computed tomography (CT) (matrix: 512 × 512, field of view (FOV): 400 mm, slice thickness: 2.2 mm, slice spacing: 1 mm, Brilliance CT 64-channel, Philips, Netherlands) and magnetic resonance (MR) (T2 weighted fast field echo (FFE) sequence, repetition time (TR): 575 ms, echo time (TE): 11.51 ms, flip angle: 20°, matrix: 512 × 512, FOV: 130 mm, slice thickness: 2 mm, slice spacing: 2.2 mm, Achieva 3.0 T TX, Philips, Netherlands) images of the shoulder bones of tested dogs. We tried to maintain the location and orientation of shoulder joints in each imaging step to minimize potential artifacts in MR images which might be generated from the magic angle effect [32]. The shoulder bones and custom markers were automatically reconstructed from the CT images using the OSIRIX (Pixmeo, Switzerland) and GeoMagic (Research Triangle Park, NC, USA) software, while the articular cartilage layers were manually segmented from the MR images using a custom MATLAB code. The 3D bone and cartilage models were then combined by aligning the subchondral bone profiles of each model (best-fit alignment in GeoMagic software). The final shoulder models included scapular and humeral bones, articular cartilage, and custom markers (Fig. 2a).

Fig. 2
figure 2

Registration of 3D shoulder models to motion analysis data. a Image-based 3D shoulder models with scapular and humeral bones, articular cartilage, and custom markers were aligned with the marker positions from the motion tracking data. b The relative positions between the centers of the humeral head and glenoid cavity were calculated with six degrees of freedom motions (an example of the rotational motion between the Y axis of the humeral head and glenoid cavity)

Supraspinatus tendon resection model

After completing articular cartilage contact strain and T2* relaxation time maps for intact shoulders, we completely resected the supraspinatus tendon in one of the shoulders while the contralateral shoulder remained intact. The resection area was covered with penrose drain tubes to prevent a recovery of the resected muscle. The dogs were allowed to perform normal activities in our animal research laboratory for 3 months following the resection surgeries (cage type: SUS304 stainless steel frame cage with fiber reinforced plastic (FRP) bottom plate (3.4 m × 1.3 m × 2.4 m), one dog per each cage with automatic watering system, temperature: 20 ± 2 °C, humidity: 50 ± 10%, light cycle (12 h): 7 am to 7 pm). Measurements of the articular cartilage contact strain and T2* relaxation time patterns were then repeated for both shoulders in the supraspinatus-resected dogs. All surgical processes were conducted under general anesthesia with an intramuscular injection of 1) satropine sulfate (0.1 ml/kg of body weight, DAI HAN PHARM. Co., Ltd., Seoul, Korea) and 2) the mixture (0.2 ml/kg of body weight) of xylazine (Rompun, Bayer Korea, Seoul, Korea) and tiletamine-zolazepam (Zoletil 50, Virbac, Carros, France).

Experimental outcomes

Gait pattern and articular cartilage strain measurement

The completed 3D shoulder models were aligned with the marker positions from the motion tracking data to determine the locations of shoulder bones in each gait frame (Fig. 2a). The relative positions between the centers of the humeral head and glenoid cavity were calculated with six degrees of freedom (DOFs) (three translational and three rotational motions in a Cartesian coordinate system) (Fig. 2b). An example of the rotational motion between the Y axis of the humeral head and glenoid cavity is shown in Fig. 2b. We removed any abnormal gait cycles that were outside of the mean ± one standard deviation (red shaded area in Fig. 2b) from all six DOF components. The remaining gait cycles were then averaged to generate a representative gait pattern for each dog (solid red line in Fig. 2b).

Articular cartilage thickness was determined by calculating the perpendicular distance from the subchondral bone interface to the articular cartilage surface for both the scapular and humeral cartilage. Articular cartilage contact strain was defined as the ratio between the undeformed articular cartilage thickness and the thickness of the overlapping areas between the scapular and humeral articular cartilage in each gait fame (Fig. 3a) [33, 34]. The articular cartilage contact areas were then defined as the cartilage areas in which the overlapping articular cartilage thickness was greater than 0.25 mm (in-plane resolution of the MR images). We created a 3D articular cartilage contact strain map for each gait frame from an entire representative gait cycle and combined all of these strain maps to generate a cumulative contact strain distribution. In this cumulative strain map, a cumulative contact area was the combination of individual contact areas from each gait frame.

Fig. 3
figure 3

3D maps of articular cartilage contact strain and T2* relaxation time. a Articular cartilage contact strain was defined as the ratio between the undeformed articular cartilage thickness and the thickness of the overlapping areas between the scapular and humeral articular cartilage. b The T2* relaxation time of each pixel in the MR images was calculated and each articular cartilage surface was divided into multiple anterior-to-posterior 20° regions

Articular cartilage T2 star (T2*) relaxation time measurement

A qMRI scan of each tested dog was performed on the same day as the morphological CT and MR imaging. A T2 weighted multi-echo fast-field sequence (TR: 700 ms, TE: 3.83/9.37/14.91/20.46/26.01 ms, flip Angle: 25°, matrix: 768 × 768, FOV: 150 mm, slice thickness: 2 mm, slice spacing: 2.2 mm) was used in a 3 T clinical MRI scanner (Achieva 3.0 T TX, Philips, Netherlands). The T2* relaxation time of each pixel in the MR images was calculated using mono-exponential least squares fitting [35]. T2* relaxation time values over 50 ms were removed to minimize the partial effects of synovial fluid around articular cartilage boundaries because T2 relaxation time of synovial fluid is known to be approximately 100 ms [36], and the T2* relaxation time is known to be around 50% of the T2 relaxation time values [27]. We divided each articular cartilage surface into multiple 20° regions from the posterior to the anterior ends (Fig. 3b) because our contact strain maps indicated anterior-to-posterior directional variations after the resection surgeries, but the changes in the medial-to-lateral direction were minimal.

Accuracy and reproducibility of contact strain measurements

Intra- and inter-observer reliability of manually segmented articular cartilage models was tested. Segmentation of articular cartilage surface and subchondral bone interface in a shoulder joint was done by three independent researchers in three difference days. Multiple 3 dimensional articular cartilage models were generated from various segmentation results for a shoulder joint. Articular cartilage thickness values at 70 random locations were compared in different cartilage models. Intra-class correlation coefficients (ICC) with 95% confident interval (CI) and average root-mean-square (RMS) differences in the thickness measurements among different segmentation results were calculated. We also calculated variations in average T2* relaxation time when articular cartilage thickness was increased or decreased within the average RMS difference in the thickness measurements.

Because articular cartilage contact strains were directly determined based on the location of the scapular and humeral bones, we decided to measure the accuracy and reproducibility of our motion-tracking-based 3D shoulder models using a custom plastic phantom model. The phantom consisted of two rectangular blocks with a plastic ball attached at one end of each block. The phantom blocks were positioned in three different configurations (relative angles of 0°, 30°, 60°). Two custom markers were installed in each phantom block to emulate the in vivo experiments (Fig. 4). We moved the phantom model in the motion laboratory in random directions for 10 min while the motion tracking system continuously recorded the positions of the markers. We created a 3D model of the phantom blocks and measured 1) the center distance between the two phantom balls and 2) the angle between the two phantom blocks for 2000 randomly selected frames. The bias (average differences in the motion-tracking-based measurements from the physical measurements of the phantom model) and precision (variations in the motion-tracking-based measurements) of the center distance and angle measurements were then calculated to find systemic errors in motion tracking system. Finally, corresponding bias and precision of articular cartilage contact strain and contact area ratio (ratio between contact area and total cartilage surface area) in in vivo shoulder models were estimated by changing the position of scapular and humeral bones within the range of bias and precision values measured in the phantom experiments at each configuration (0°, 30°, 60° angles between scapula and humerus).

Fig. 4
figure 4

A phantom model for the accuracy and reproducibility evaluation. The phantom consisted of two rectangular blocks with a plastic ball attached at one end of each block (three different block angles of 0°, 30°, 60°)

Statistical methods

Intra- and inter-observer ICC with 95% confidence interval (CI) of segmentations were analyzed by using SPSS software (SPSS statistics 25, SPSS Inc., IL, USA). Pearson correlation coefficients between articular cartilage contact strain and T2* relaxation time were calculated by using MATLAB software (R2016a, MathWorks, MA, USA). A p-value less than 0.05 was considered to be a statistically significant difference.

Results

Accuracy and reproducibility of contact strain measurements

Excellent intra- and inter-observer reliability of cartilage thickness measurements in different segmentations were found with ICCs over 0.850. Average RMS differences in the thickness measurements between cartilage models from 3 different days and observers were 0.097 mm and 0.125 mm respectively (Table 1). Variations in the average T2* relaxation time were 0.496 ± 0.202 ms when the articular cartilage thickness varied within 0.125 mm which was the maximum value among intra- and inter-observer RMS differences in cartilage thickness measurements.

Table 1 Intra- and inter-observer reliability of cartilage thickness measurements

Intra-class correlation coefficients (ICC) with 95% confident interval (CI) and average root-mean-square (RMS) differences in the thickness measurements between different segmentation results were calculated.

The bias and precision of motion-tracking-guided measurements in the phantom were less than 0.050 mm for the center distance and 0.692 degree for the angle measurements (Table 2). Estimated errors (bias) in the articular cartilage contact strain and contact area ratio were found to be less than 0.349 and 0.570% respectively. Variations in the measurements (precision) were less than 0.696% for average contact strains and 1.020% for contact area ratios (Table 2).

Table 2 Bias and precision of motion-tracking-guided measurements

The bias and precision of the center distance and angle measurements in the phantom model are listed, and corresponding bias and precision of articular cartilage contact strain and contact area ratio in in vivo shoulder models are also estimated.

Variations in articular cartilage contact strain and T2* relaxation time after supraspinatus resection

In the 3D shoulder models of all tested dogs, average and maximum thickness values of articular cartilage were measured at 0.819 ± 0.086 mm and 1.227 ± 0.166 mm respectively. Unilateral supraspinatus resection surgeries were found to alter articular cartilage contact strain patterns in both shoulders 3 months after the supraspinatus resection surgeries. An example of the articular cartilage cumulative contact strain patterns for one of the tested shoulders is presented in Fig. 5. The average articular cartilage cumulative contact strain was reduced in the supraspinatus-resected shoulders by approximately 2.145%. However, the average cumulative contact strain in the contralateral shoulders increased by approximately 2.243% (Fig. 6).

Fig. 5
figure 5

Changes in cumulative contact strain and T2* relaxation time distributions after a supraspinatus resection surgery. An example of the articular cartilage cumulative contact strain and T2* relaxation time patterns before and after the supraspinatus resection surgery for one of the tested shoulders are presented

Fig. 6
figure 6

Changes in average cumulative contact strain and T2* relaxation time after supraspinatus resections. Changes in the average articular cartilage cumulative contact strain and in the average T2* relaxation time are shown after the supraspinatus resection surgeries

Average T2* relaxation time variations after the supraspinatus resection surgery exhibited similar patterns with the patterns in the average cumulative contact strain. The average articular cartilage T2* relaxation time in the supraspinatus-resected shoulders decreased, whereas the T2* relaxation time in the intact shoulders increased over 2.316 ms after the resection surgeries (Fig. 6).

Regional changes in contact strain and T2* relaxation time

Because cartilage contact strain and T2* relaxation time were found to concentrate in local contact areas (Fig. 5), we decided to calculate regional correlations between the changes in articular cartilage contact strain and T2* relaxation time following the supraspinatus resections in each of the anterior-to-posterior 20° regions for all tested shoulders. The regional variations in articular cartilage contact strain and T2* relaxation time exhibited a positive linear correlation with a Pearson correlation coefficient (r) of 0.726 (p < 0.001, Fig. 7).

Fig. 7
figure 7

Correlation between regional changes in contact strain and T2* relaxation time following the supraspinatus resections. The regional variations in articular cartilage contact strain and T2* relaxation time exhibited a linear correlation in both the intact and supraspinatus-resected shoulders

Discussion

Injury-related articular cartilage degeneration is believed to be the consequence of mechanical alterations. However, the role of altered joint mechanics in articular cartilage degeneration in in vivo models has not been studied successfully due to the lack of information regarding pre-injury articular cartilage conditions. This is the first in vivo study to compare articular cartilage contact strains and T2 star (T2*) relaxation times before and after a joint injury. We successfully developed a technique for measuring in vivo articular cartilage contact strains using a motion tracking system and image-based 3D shoulder models with high accuracy and reproducibility. Unilateral supraspinatus tendon resection altered articular cartilage contact patterns in both shoulders. Contact strains in supraspinatus-resected shoulders were found to decrease, whereas the contact strains in the contralateral intact shoulders increased. Interestingly, these changes in articular cartilage contact strain exhibited similar patterns with the changes in the articular cartilage T2* relaxation time following the supraspinatus resection surgeries throughout all tested shoulders (Fig. 6). It was further supported by positive correlations between regional variations in articular cartilage contact strain and T2* relaxation time (Fig. 7). These results are consistent with previous clinical studies which reported increased T2* relaxation time on the regions of increased articular cartilage contact loads [37,38,39]. Although those clinical measurements must be the result of long-term response to altered joint loading, our relatively short three-month study suggests that the changes in the articular cartilage extracellular matrix may be initiated at the early stage of load variations.

Because articular cartilage contact strains were calculated directly from position measurements of the shoulder bones, we attempted to estimate the accuracy and reproducibility of the position measurements of the proposed motion-tracking-derived 3D models using a plastic phantom. The errors of the position measurements were less than 50 μm in terms of translation and less than 0.7° in terms of angular displacement for all three configurations. The errors in the phantom experiments were expected to generate approximately 0.4% errors in the average contact strain and 0.6% errors in the contact area ratio measurements in in vivo shoulder models. These errors were a combination of systemic errors in the motion tracking device, reconstruction errors in the image-based 3D models, and alignment errors between the motion tracking information and 3D models. The results demonstrate that our technique is as accurate as biplane X-ray radiostereometric analysis, which is typically used for minimized and limited joint motions (accuracy of approximately 10–250 μm in terms of translation and 0.03–0.6° in terms of angular displacement) [40, 41]. The accuracy of our method was found to be better than the accuracy of the popular 2D/3D registration method of biplane fluoroscopy for real-time joint motion analysis (accuracy of approximately 200–700 μm in terms of translation and 0.1–0.9° in terms of angular displacement) [42,43,44]. Furthermore, our motion-tracking-based method has a wide measurement area, whereas X-ray-based methods have limited measurement space between X-ray sources and detectors. In our in vivo animal study, a large measurement space was crucial to allow natural joint motions without any pre-training.

Unilateral supraspinatus tendon resections altered the articular cartilage contact strain patterns in both shoulders. Decreases in the mean contact strain were found in the supraspinatus-resected shoulders. However, the magnitude of the mean contact strain in the contralateral intact joints markedly increased. This may be the result of compensatory weight shifting behaviors by the dogs to avoid discomfort without losing joint stability following the supraspinatus resection surgeries [45]. T2* relaxation time in the joints with increased contact strain was also increased by approximately 2.5 ms (Fig. 6). Previously, T2* relaxation time in osteoarthritic articular cartilage were reported to be approximately 2 to 5 ms greater than the values in healthy articular cartilage [28, 35, 37]. Although our 2.5 ms increases in T2* relaxation time were smaller than the increases in osteoarthritic cartilages, the increased T2* relaxation time in the region of increased contact strain may indicate the early onset of articular cartilage degenerative processes, which have been identified in abnormally overloaded cartilage areas in various studies [46, 47]. Because articular cartilage T2* relaxation time is related to collagen network integrity and water content in cartilage matrices [12, 16, 48], prominent positive correlations between regional variations in articular cartilage contact strain and T2* relaxation time (Fig. 7) suggests that the articular cartilage extracellular matrices responded to the altered loading patterns, resulting in structural and compositional changes within a relatively short three-month testing period.

This study has a few limitations that must be discussed to contextualize the presented results. First, we only completed six in vivo canine shoulder models from three dogs, despite our original intent of using five dogs for our experiments. One dog died from complications during the anesthesia process for CT scanning. We could not collect motion tracking data for a second dog because the dog did not want to walk in the motion laboratory following the supraspinatus resection surgery. This limited sample size did not allow us to establish various statistical examinations between average contact strain and T2* relaxation time before and after the supraspinatus resections. However, the patterns of increased T2* relaxation time in the shoulders with increased contact strain and decreased T2* relaxation time in the shoulders with decreased strain were consistent for all tested shoulders. The linear correlations between regional contact strain and T2* relaxation time changes in both the intact and supraspinatus-resected shoulders were statistically significant, which supports our hypothesis that alterations in joint mechanics produce important structural and compositional changes in articular cartilage extracellular matrices. Second, the gait patterns recorded by our motion tracking system were not consistent for each gait cycle because joint angles and walking speed varied during normal free walking episodes. Previous studies have shown that a treadmill can produce consistent gait patterns [49]. However, we chose to record the gait cycles of free-walking subjects using a motion tracking system because 1) treadmill walking requires many training sessions for test subjects to become familiar with the system and 2) we wished to avoid any trained gaits using the treadmill, which could have created biased joint mechanics. In this study, we only selected gait cycles from dogs walking in straight lines. We then excluded gait cycles outside of the mean ± one standard deviation range from the selected gaits. The remaining gait cycles were averaged to generate a consistent representative gait pattern. We feel that this final gait pattern successfully described a representative gait motion for each tested dog by including more than 70% of the recorded gait cycles.

Conclusion

This is the first study to investigate the relationship between articular cartilage contact strain and T2* relaxation time before and after a joint injury in an in vivo animal model. We intentionally created unilateral supraspinatus tendon resections in in vivo canine models to alter joint mechanics, and measured articular cartilage contact strain and T2* relaxation time before and 3 months after the supraspinatus resections. Interestingly, patterns in the articular cartilage contact strain changes were similar to the patterns in the articular cartilage T2* relaxation time variations following supraspinatus resection surgeries. Regional comparisons between articular cartilage contact strain and T2* relaxation time revealed positive correlations in all tested shoulders. Because increased T2* relaxation time is believed to be early signs of articular cartilage degeneration, the degenerative process of articular cartilage may be initiated in areas with increased cartilage strain, potentially leading to the development of osteoarthritis [28, 35].

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ACL:

Anterior cruciate ligament

CI:

Confidential interval;

CT:

Computed tomography

DOF:

Degrees of freedom

FFE:

Fast-field-echo

FOV:

Field-of-view

FRP:

Fiber reinforced plastic

ICC:

Intraclass correlation coefficients

MR:

Magnetic resonance

OA:

Osteoarthritis

qMRI:

Quantitative magnetic resonance imaging

RMS:

Root-mean-square

SNR:

Signal to noise ratio

TE:

Time to echo

3D:

Three dimension

TR:

Time to repeat

T2*:

T2 star

References

  1. Maldonado M, Nam J. The role of changes in extracellular matrix of cartilage in the presence of inflammation on the pathology of osteoarthritis. Biomed Res Int. 2013;2013:284873.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Jorgensen AEM, Kjaer M, Heinemeier KM. The effect of aging and mechanical loading on the metabolism of articular cartilage. J Rheumatol. 2017;44(4):410–7.

    Article  PubMed  CAS  Google Scholar 

  3. Bader DL, Salter DM, Chowdhury TT. Biomechanical influence of cartilage homeostasis in health and disease. Arthritis. 2011;2011:979032.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Brama PA, Tekoppele JM, Bank RA, Karssenberg D, Barneveld A, van Weeren PR. Topographical mapping of biochemical properties of articular cartilage in the equine fetlock joint. Equine Vet J. 2000;32(1):19–26.

    Article  CAS  PubMed  Google Scholar 

  5. Rogers BA, Murphy CL, Cannon SR, Briggs TW. Topographical variation in glycosaminoglycan content in human articular cartilage. J Bone Joint Surg Br. 2006;88(12):1670–4.

    Article  CAS  PubMed  Google Scholar 

  6. Armstrong SJ, Read RA, Price R. Topographical variation within the articular cartilage and subchondral bone of the normal ovine knee joint: a histological approach. Osteoarthr Cartil. 1995;3(1):25–33.

    Article  CAS  PubMed  Google Scholar 

  7. Stenhamre H, Slynarski K, Petren C, Tallheden T, Lindahl A. Topographic variation in redifferentiation capacity of chondrocytes in the adult human knee joint. Osteoarthr Cartil. 2008;16(11):1356–62.

    Article  CAS  PubMed  Google Scholar 

  8. Brama PA, Tekoppele JM, Bank RA, Barneveld A, van Weeren PR. Functional adaptation of equine articular cartilage: the formation of regional biochemical characteristics up to age one year. Equine Vet J. 2000;32(3):217–21.

    Article  CAS  PubMed  Google Scholar 

  9. Han EH, Chen SS, Klisch SM, Sah RL. Contribution of proteoglycan osmotic swelling pressure to the compressive properties of articular cartilage. Biophys J. 2011;101(4):916–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Canal Guterl C, Hung CT, Ateshian GA. Electrostatic and non-electrostatic contributions of proteoglycans to the compressive equilibrium modulus of bovine articular cartilage. J Biomech. 2010;43(7):1343–50.

    Article  PubMed  Google Scholar 

  11. Stolz M, Gottardi R, Raiteri R, Miot S, Martin I, Imer R, et al. Early detection of aging cartilage and osteoarthritis in mice and patient samples using atomic force microscopy. Nat Nanotechnol. 2009;4(3):186–92.

    Article  CAS  PubMed  Google Scholar 

  12. Choi JA, Gold GE. MR imaging of articular cartilage physiology. Magn Reson Imaging Clin N Am. 2011;19(2):249–82.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Keenan KE, Besier TF, Pauly JM, Han E, Rosenberg J, Smith RL, et al. Prediction of glycosaminoglycan content in human cartilage by age, T1rho and T2 MRI. Osteoarthr Cartil. 2011;19(2):171–9.

    Article  CAS  PubMed  Google Scholar 

  14. Li X, Cheng J, Lin K, Saadat E, Bolbos RI, Jobke B, et al. Quantitative MRI using T1rho and T2 in human osteoarthritic cartilage specimens: correlation with biochemical measurements and histology. Magn Reson Imaging. 2011;29(3):324–34.

    Article  CAS  PubMed  Google Scholar 

  15. Blumenkrantz G, Majumdar S. Quantitative magnetic resonance imaging of articular cartilage in osteoarthritis. Eur Cell Mater. 2007;13:76–86.

    Article  CAS  PubMed  Google Scholar 

  16. Guermazi A, Alizai H, Crema MD, Trattnig S, Regatte RR, Roemer FW. Compositional MRI techniques for evaluation of cartilage degeneration in osteoarthritis. Osteoarthr Cartil. 2015;23(10):1639–53.

    Article  CAS  PubMed  Google Scholar 

  17. Baum T, Joseph GB, Karampinos DC, Jungmann PM, Link TM, Bauer JS. Cartilage and meniscal T2 relaxation time as non-invasive biomarker for knee osteoarthritis and cartilage repair procedures. Osteoarthr Cartil. 2013;21(10):1474–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Prasad AP, Nardo L, Schooler J, Joseph GB, Link TM. T (1) rho and T (2) relaxation times predict progression of knee osteoarthritis. Osteoarthr Cartil. 2013;21(1):69–76.

    Article  CAS  PubMed  Google Scholar 

  19. Apprich S, Mamisch TC, Welsch GH, Stelzeneder D, Albers C, Totzke U, et al. Quantitative T2 mapping of the patella at 3.0T is sensitive to early cartilage degeneration, but also to loading of the knee. Eur J Radiol. 2012;81(4):e438–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Lusse S, Claassen H, Gehrke T, Hassenpflug J, Schunke M, Heller M, et al. Evaluation of water content by spatially resolved transverse relaxation times of human articular cartilage. Magn Reson Imaging. 2000;18(4):423–30.

    Article  CAS  PubMed  Google Scholar 

  21. Souza RB, Baum T, Wu S, Feeley BT, Kadel N, Li X, et al. Effects of unloading on knee articular cartilage T1rho and T2 magnetic resonance imaging relaxation times: a case series. J Orthop Sports Phys Ther. 2012;42(6):511–20.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Hovis KK, Stehling C, Souza RB, Haughom BD, Baum T, Nevitt M, et al. Physical activity is associated with magnetic resonance imaging-based knee cartilage T2 measurements in asymptomatic subjects with and those without osteoarthritis risk factors. Arthritis Rheum. 2011;63(8):2248–56.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Munukka M, Waller B, Rantalainen T, Hakkinen A, Nieminen MT, Lammentausta E, et al. Efficacy of progressive aquatic resistance training for tibiofemoral cartilage in postmenopausal women with mild knee osteoarthritis: a randomised controlled trial. Osteoarthr Cartil. 2016;24(10):1708–17.

    Article  CAS  PubMed  Google Scholar 

  24. Kim CW, Hosseini A, Lin L, Wang Y, Torriani M, Gill T, et al. Quantitative analysis of T2 relaxation times of the patellofemoral joint cartilage 3 years after anterior cruciate ligament reconstruction. J Orthop Translat. 2018;12:85–92.

    Article  PubMed  Google Scholar 

  25. Gao B, Zheng NN. Alterations in three-dimensional joint kinematics of anterior cruciate ligament-deficient and -reconstructed knees during walking. Clin Biomech (Bristol, Avon). 2010;25(3):222–9.

    Article  Google Scholar 

  26. Hosseini A, Van de Velde S, Gill TJ, Li G. Tibiofemoral cartilage contact biomechanics in patients after reconstruction of a ruptured anterior cruciate ligament. J Orthop Res. 2012;30(11):1781–8.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Mamisch TC, Hughes T, Mosher TJ, Mueller C, Trattnig S, Boesch C, et al. T2 star relaxation times for assessment of articular cartilage at 3 T: a feasibility study. Skelet Radiol. 2012;41(3):287–92.

    Article  Google Scholar 

  28. Marik W, Apprich S, Welsch GH, Mamisch TC, Trattnig S. Biochemical evaluation of articular cartilage in patients with osteochondrosis dissecans by means of quantitative T2- and T2-mapping at 3T MRI: a feasibility study. Eur J Radiol. 2012;81(5):923–7.

    Article  CAS  PubMed  Google Scholar 

  29. Teng HL, Pedoia V, Link TM, Majumdar S, Souza RB. Local associations between knee cartilage T1rho and T2 relaxation times and patellofemoral joint stress during walking: a voxel-based relaxometry analysis. Knee. 2018;25(3):406–16.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Amano K, Pedoia V, Su F, Souza RB, Li X, Ma CB. Persistent biomechanical alterations after ACL reconstruction are associated with early cartilage matrix changes detected by quantitative MR. Orthop J Sports Med. 2016;4(4):2325967116644421.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Roberts S, Weightman B, Urban J, Chappell D. Mechanical and biochemical properties of human articular cartilage in osteoarthritic femoral heads and in autopsy specimens. J Bone Joint Surg Br. 1986;68(2):278–88.

    Article  CAS  PubMed  Google Scholar 

  32. Mosher TJ, Smith H, Dardzinski BJ, Schmithorst VJ, Smith MB. MR imaging and T2 mapping of femoral cartilage: in vivo determination of the magic angle effect. AJR Am J Roentgenol. 2001;177(3):665–9.

    Article  CAS  PubMed  Google Scholar 

  33. Liu F, Kozanek M, Hosseini A, Van de Velde SK, Gill TJ, Rubash HE, et al. In vivo tibiofemoral cartilage deformation during the stance phase of gait. J Biomech. 2010;43(4):658–65.

    Article  PubMed  Google Scholar 

  34. Carter TE, Taylor KA, Spritzer CE, Utturkar GM, Taylor DC, Moorman CT. 3rd, et al. In vivo cartilage strain increases following medial meniscal tear and correlates with synovial fluid matrix metalloproteinase activity. J Biomech. 2015;48(8):1461–8.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Tsai PH, Lee HS, Siow TY, Chang YC, Chou MC, Lin MH, et al. Sequential change in T2* values of cartilage, meniscus, and subchondral bone marrow in a rat model of knee osteoarthritis. PLoS One. 2013;8(10):e76658.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Souza RB, Kumar D, Calixto N, Singh J, Schooler J, Subburaj K, et al. Response of knee cartilage T1rho and T2 relaxation times to in vivo mechanical loading in individuals with and without knee osteoarthritis. Osteoarthr Cartil. 2014;22(10):1367–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Hu Y, Tao H, Qiao Y, Ma K, Hua Y, Yan X, et al. Evaluation of the Talar cartilage in chronic lateral ankle instability with lateral ligament injury using biochemical T2* mapping: correlation with clinical symptoms. Acad Radiol. 2018;25(11):1415–21.

    Article  PubMed  Google Scholar 

  38. Behzadi C, Welsch GH, Laqmani A, Henes FO, Kaul MG, Schoen G, et al. Comparison of T2* relaxation times of articular cartilage of the knee in elite professional football players and age-and BMI-matched amateur athletes. Eur J Radiol. 2017;86:105–11.

    Article  CAS  PubMed  Google Scholar 

  39. Behzadi C, Maas KJ, Welsch G, Kaul M, Schoen G, Laqmani A, et al. Quantitative T2 * relaxation time analysis of articular cartilage of the tibiotalar joint in professional football players and healthy volunteers at 3T MRI. J Magn Reson Imaging. 2018;47(2):372–9.

    Article  PubMed  Google Scholar 

  40. Karrholm J. Roentgen stereophotogrammetry. Review of orthopedic applications. Acta Orthop Scand. 1989;60(4):491–503.

    Article  CAS  PubMed  Google Scholar 

  41. Tashman S, Anderst W. In-vivo measurement of dynamic joint motion using high speed biplane radiography and CT: application to canine ACL deficiency. J Biomech Eng. 2003;125(2):238–45.

    Article  PubMed  Google Scholar 

  42. Bey MJ, Zauel R, Brock SK, Tashman S. Validation of a new model-based tracking technique for measuring three-dimensional, in vivo glenohumeral joint kinematics. J Biomech Eng. 2006;128(4):604–9.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Bey MJ, Kline SK, Tashman S, Zauel R. Accuracy of biplane x-ray imaging combined with model-based tracking for measuring in-vivo patellofemoral joint motion. J Orthop Surg Res. 2008;3:38.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Anderst W, Zauel R, Bishop J, Demps E, Tashman S. Validation of three-dimensional model-based tibio-femoral tracking during running. Med Eng Phys. 2009;31(1):10–6.

    Article  PubMed  Google Scholar 

  45. Mikola K, Piras A, Hakala L. Isolated avulsion of the tendon of insertion of the Infraspinatus and supraspinatus muscles in five juvenile Labrador retrievers. Vet Comp Orthop Traumatol. 2018;31(4):285–90.

    Article  PubMed  Google Scholar 

  46. Song Y, Carter DR, Giori NJ. Cartilage nominal strain correlates with shear modulus and glycosaminoglycans content in meniscectomized joints. J Biomech Eng. 2014;136(6):064503.

    Article  PubMed  Google Scholar 

  47. Arunakul M, Tochigi Y, Goetz JE, Diestelmeier BW, Heiner AD, Rudert J, et al. Replication of chronic abnormal cartilage loading by medial meniscus destabilization for modeling osteoarthritis in the rabbit knee in vivo. J Orthop Res. 2013;31(10):1555–60.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Andreisek G, Weiger M. T2* mapping of articular cartilage: current status of research and first clinical applications. Investig Radiol. 2014;49(1):57–62.

    Article  Google Scholar 

  49. Owen MR, Richards J, Clements DN, Drew ST, Bennett D, Carmichael S. Kinematics of the elbow and stifle joints in greyhounds during treadmill trotting - an investigation of familiarisation. Vet Comp Orthopaed. 2004;17(3):141–5.

    Article  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This work was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (NRF-2012R1A1A3010896 and NRF-2015R1C1A1A0203708).

Author information

Authors and Affiliations

Authors

Contributions

DL contributed to the acquisition, analysis and interpretation of data, and drafting the manuscript. KH participated in the acquisition and analysis of data. YHL, JSP and WK contributed to the surgery and the treatment of tested animals. TSL and EL contributed to the acquisition and interpretation of medical images. JHO contributed to the design of the surgery and the treatment of tested animals, drafting and revising the manuscript, and approval of the final submitted article. JAC contributed to the design of medical imaging procedures, drafting and revising the manuscript, and approval of the final submitted article. YS contributed to research design, interpretation of data, drafting and revising the manuscript, and approval of the final submitted article. All authors have read and approved the final submitted manuscript.

Corresponding author

Correspondence to Yongnam Song.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Animal Care and Use Committee (approval number: BA1507–180/047–01) and carried out in accordance with the Guide for the Care and Use of Laboratory Animals of the Seoul National University College of Medicine (Republic of Korea).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lee, D., Hong, KT., Lim, T.S. et al. Alterations in articular cartilage T2 star relaxation time following mechanical disorders: in vivo canine supraspinatus tendon resection models. BMC Musculoskelet Disord 21, 424 (2020). https://doi.org/10.1186/s12891-020-03447-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-020-03447-3

Keywords