T1ρ magnetic resonance imaging quantification of early articular cartilage degeneration in a rabbit model

Background Osteoarthritis (OA) is a serious problem in the recent aging society, and early diagnosis and intervention of articular cartilage degeneration are very important for the onset of OA. Therefore, development of newer MRI techniques is necessary and expected for detection of early articular cartilage degeneration. Methods 24 rabbits were randomly divided into four equal experimental groups (Group A, B, C, D) to establish articular cartilage models in different grades of early degeneration by injecting papain into the left knee joint cavity. Another 8 rabbits were considered as blank control (Group E), and then randomized into four subgroups (EA, EB, EC, ED). T1ρ and T2-weighted images of the bilateral knee joints were obtained for rabbits by using 3.0 T MRI. Group A, B, C, and D were imaged respectively at 1, 2, 3, and 4 weeks post-operation, and EA, EB, EC, ED underwent the same period imaging. Rabbits were sacrificed after scanning and the femoral condyle cartilage (FCC) was histological examined. T1ρ values of the femoral condyle cartilage were measured and statistically analyzed, and contrasted with the histologic results. Results T1ρ values of the left side in experimental groups were significantly higher than the right side (P < 0.05), and which increased gradually with the passage of post-operation time (P < 0.05). Histological examination demonstrated the proteoglycan content of the left side decreased, and indicated the occurrence of early degeneration. Conclusions T1ρ MRI can sensitively and quantitatively reflect the change in proteoglycans prior to the morphologic alterations of articular cartilage, and T1ρ value is gradually increased with a decrease in proteoglycan content, therefore that T1ρ could potentially act as a reliable tool to identify early cartilage degeneration.


Background
Osteoarthritis (OA) is the most common joint disorder worldwide, which is the major cause of mobility impaired and disabled in humans [1,2]. It affects approximately 75 % of the population over 70 years of age, while a quarter of people aged over 55 have an episode of persistent knee pain, and there is a rapid increase in the number of the cases due to the aging of the population and the obesity epidemic [2][3][4][5]. This suggests that OA will continue to be a large and growing public health problem in the future [5,6].
The most recent accomplishments have indicated that articular cartilage degeneration is highly likely to play a key role in the pathogenesis of OA, and early diagnosis and intervention of cartilage degeneration before the onset of irreversible changes, are very important for the prevention and treatment of OA [2,7,8]. For this purpose, considerable efforts have recently been in search of non-invasive techniques to sensitively and accurately assess the prophase of cartilage degeneration [9]. Articular cartilage starts early to degenerate from extracellular matrix (ECM) metabolism imbalance, which mainly involves a gradual loss of proteoglycan (PG), collagen damage and increased water content, and these macromolecular (biochemical) changes are not accompanied by significant structural disturbances in the tissue [10][11][12][13]. Among them, decreased PG content is an initiating factor and an important marker of early cartilage degeneration. Presently, it is believed that cartilage destruction results from large PG aggrecan degradation, which precedes the damage to collagen fibrillar network and is responsible for the whole changes [12][13][14][15]. In short, early cartilage changes are characterized by PGs depletion, and detection of PG content can identify early cartilage degeneration. Although conventional magnetic resonance imaging (MRI) could provide excellent detection of morphological changes in the thickness and structure of articular cartilage, it has some difficulties in discriminating minimal changes associated with the early stage of degeneration before morphological or clinical alterations [2,12,13]. Given as such, newer MRI techniques have emerged for investigation of the physiology and biochemistry of cartilage.
T 1ρ MRI is one of the most promising techniques focused on PG quantification [2,12,16]. The longitudinal relaxation time that occurs during the application of T 1ρ imaging is referred to as spin-lattice relaxation in the rotating frame, or T 1ρ relaxation time. T 1ρ relaxation time (T 1ρ value) can be measured to describe the slow-motion interactions between the macromolecule protons and bulk water protons at high static fields in vivo, and therefore indirectly reflect PG content [17][18][19][20]. As a consequence, T 1ρ imaging is mainly used for the evaluation of tissues composed partly of slow-frequency macromolecules, such as brain, articular cartilage, and intervertebral disc [21][22][23]. In articular cartilage, several studies have reported that a decrease in PG content is the major factor in change of T 1ρ value, and which has been shown to detect cartilage changes [13,23,24]. However, T 1ρ is still available on a research basic, and further studies are required to evaluate its value in early degeneration, and optimize the scanning sequence to create a reliable quantitative scale [2,[25][26][27]. Therefore, we designed the present study to confirm the sensibility and reliablity of T 1ρ in quantification of early cartilage degeneration with histology.

Ethics statement
The experimental protocol was approved by the Animal Care and Use Committee of Jinan University, Guangzhou, China. All procedures performed in this study involving animals were in accordance with the revised Animals (Scientific Procedures) Act 1986 in the UK and Directive 2010/63/EU in Europe.

Preliminary experiment
In order to investigate the feasibility of this study, including training operations and sample size estimation, two healthy New Zealand white rabbits were performed the preliminary experiment first. Sample size was determined by the calculation of the pre-test results and the relevant literatures [12,28].

Animal model
Thirty-two healthy New Zealand white rabbits (16 male and 16 female, with a mean age of 8 months and weighting from 2.0 to 3.0 kg) were used for the formal experiment. Rabbits were randomly numbered 1-32, and divided into 5 groups. Numbers 1 to 6 were included in Group A, 7 to 12 in Group B, 13 to 18 in Group C, 19 to 24 in Group D, and numbers 25 to 32 were included in Group E (among the 8 rabbits in Group E, 2 rabbits were randomly selected for each subgroup E A , E B , E C , and E D , which were considered as the blank controls of Group A, B, C, and D, respectively). The left knees of experimental groups (Group A, B, C, D) were served as the experimental side to establish articular cartilage models at different stages of early degeneration with papain, which is a known and validated mediator of inducing PG catabolism [28][29][30][31], and the right knees acted as the control side. On the 1st, 4th, and 7th day of the   , the twenty-four rabbits of experimental groups were administered an injection of 0.5 ml of 1.6 % papain solution into the joint cavity of the left knee by using a 16-gauge hypodermic needle, and a control injection of the same amount of solution without papain in the right knee. The rabbits were sedated via ear vein injection of Ketamine (10 mg/kg) and Midazolam (1 mg/kg) while being monitored by an anesthesiologist for the period of the surgical procedure. Postoperatively, the rabbits were housed in individual cages without restriction of joint movement.

MRI Protocol
T 1ρ and T 2 -weighted images of the bilateral knee joints were obtained for 32 rabbits by using a clinical 3.0 T MRI scanner (GE discovery MR750; the First Affiliated Hospital of Jinan University, Guangzhou, China) and human 8ch knee array coil. Group A, B, C, and D were scanned in sequence at 1, 2, 3 and 4 weeks post-operation (from the last injection of papain) respectively, and rabbits in the subgroups (E A , E B , E C , and E D ) underwent MRI in the same period for comparison. For the duration of the imaging session (approximately 15 mins), each rabbit was under general anesthesia with intravenous Ketamine (8 mg/kg) and Midazolam (0.8 mg/kg), and the protocol included 2 sequences: the first, conventional T 2 -weighted imaging (T 2 WI) for morphological evaluation were acquired using a fast spin-echo (FSE) imaging on sagittal view, with the following parameters: repetition time/echo time (TR/TE) = 2000 ms/85 ms, slice thickness/spacing (thn/spa) = 2 mm/0 mm, field of view (FOV) = 12 × 12 cm, acquisition matrix = 352 × 224, NEX = 2; the second MRI pulse sequence used for a series of T 1ρ -weighted sagittal images was 3D Magnetization-Prepared Angle-Modulated Partitioned k-Space Spoiled Gradient Echo Snapshots (3D MAPSS) with a spin-lock pulse amplitude of 500 Hz, under the following parameters: time of spin-lock pulse (TSL) = 0, 10, 40, or 80 ms, TR/TE = 7.8 ms/3.7 ms, thn/ spa = 12 mm/0 mm, FOV = 12 × 12 cm, matrix = 256 × 128, NEX = 1, scan time = 3 min 33 s, and the B1 of spinlock radio frequency (RF) was 0.145 Gauss.

Histology
Immediately after imaging, rabbits in each group were euthanized by intravenous pentobarbital (100 mg/kg), and both knees were harvested, formalin-fixed, and decalcified. Tissue blocks of 0.5 cm in thickness were harvested from the medial and lateral femoral condyle cartilage (MFCC and LFCC) on the sagittal plane, and made into paraffin-embedded tissue sections. Morphology of the chondrocytes was observed after hematoxylin and eosin (H&E) staining, and the PG content was determined by Safranin-O and Fast Green staining [12,32].

T 1ρ quantification
Sun Advantage Workstation 4.5 (GE) was used for the post-processing of T 1ρ images. With an exponential decay, T 1ρ maps were reconstructed on a pixel-by-pixel basis by fitting a linear regression of the image intensity data to the following equation: where S is the signal intensity on T 1ρ images with the given TSL.
On the intermediate layer of the T 1ρ sagittal image, three regions of interest (ROIs) were placed from anterior to posterior on the images of MFCC and LFCC for   measurement of the T 1ρ values. Each ROI was drawn using a circular tool by the radiologist with more than 10 years of the professional experience of MRI, and in order to reduce the influence of partial volume effect on results, the ROIs were not allowed to be set on the adjacent tissue. Average values of the three ROIs were calculated and recorded as the T 1ρ values of MFCC and LFCC, which were statistically analyzed.

Image interpretation
The interpretations of images acquired were determined by two radiologists with more than 10 years of experience, which blinded to all experimental information and the final histological diagnosis, and independently readed the T1ρ original and T2-weighted images of the 32 rabbits. Each observer performed the same interpretation twice, with a delay of two weeks, to assess intraobserver agreement. The image interpretation results were determined according to the agreement and disagreement between the two observers as follows: the same results and the final results after discussion.

Statistical analysis
Interobserver and intraobserver reliability of image interpretation was evaluated with the interclass and intraclass correlation coefficient (ICC): poor (< 0.40), fair (0.40 -0.59), good (0.60 -0.74), and excellent (0.75 -1.00). The mean ± standard deviation of measurements were calculated for the T 1ρ values of MFCC and LFCC. A normal distribution of our data was shown by using the normality test, and then the differences in the T 1 ρ values of LFCC and MFCC within experimental groups were compared using one-way analysis of variance (ANOVA). Comparison of the T 1 ρ values in the control side between Group E and experimental groups was conducted with independent sample T-test. A P value of < 0.05 was considered statistically significant. Statistical analyses were performed with SPSS statistics software package, version 10.0 (SPSS Inc., Chicago, Illinois, USA).

Analysis of images
Two observers agreed perfectly that there was no imaging change of the femoral condyle cartilage (FCC) in the acquired T 1 ρ original and T 2 -weighted images of the 32 rabbits (Figures 1 and 2), in addition, there was no difference in images analysis between the first and the second interpretation by the same observer. Interclass and intraclass correlation coefficient were excellent, with k values of 1.00. T 1ρ pseudo-color images showed relatively significant chromatic aberration in the experimental sides compared with the control sides and Group E, and the color of the experimental sides were darker. Progressive changes from pale green to dark green were observed in the experimental sides of Group A-D (Fig. 3), which indicated the T 1ρ value increased gradually, and the process was consistent with the histological finding (Figs. 4 and 5).

T 1ρ value differences of intra-and inter-group
T 1ρ values of LFCC and MFCC in the bilateral knees of rabbits are shown in  Fig. 6). Fig. 7 showed that T 1ρ values of the experimental side increased gradually with the passage of post-operation time (P < 0.05), and there was no significant difference in T 1ρ values of the control side among Group A, B, and C (P > 0.05), while the T 1ρ values in Group D were obviously greater than in Group A, B, and C (P < 0.05). There was no difference in T 1ρ values on control side between Group A, B, C and the subgroup E A , E B , and E C (P > 0.05, Table 1), while the T 1ρ values of LFCC and MFCC on the control side in Group D were obviously greater than in subgroup E D (P < 0.05, Table 1) . There was no significant difference in T 1ρ values between the LFCC and MFCC, within the bilateral knees of Group A-D (P > 0.05, Table 1). Finally, No difference of T 1ρ values existed in Group E (P > 0.05).

Histodiagnosis
This study performed the injections of papain on rabbits to produce a decline in PG content for setting up the animal models of early articular cartilage degeneration, and the histodiagnosis substantiated the occurrence of biochemical changes in the FCC, which with little change in morphology. Figure 8 showed that the LFCC and MFCC of the specimens were smooth and devoid of morphologic alteration upon histological inspection. Figure 4 revealed early degenerative changes in FCC on the experimental side because of the gradually disappearing of normal chondrocytes and the fibrosis of ECM. And there was a visible reduction in histological staining of PG by Safranin-O for papain-treated FCC as compared with control FCC, and the loss of PG was proportional to the degree of cartilage degeneration, meanwhile the range of Fast Green staining of the collagen fibers gradually increased (Figs. 4 and 5). However, there were no significant change in the FCC on control side (except in Group D) and Group E; there was mild degeneration in the FCC on control side of Group D. Histological verification was determined by the observer with the professional experience of pathology.

Discussion
In the present study, the primary purpose was to investigate whether T 1ρ MRI could be uesd as a means to detect biochemical changes prior to the morphologic alterations in articular cartilage with histology. Out of the conventional MRI methods, T 2 WI is commonly used to depict the morphological changes in cartilage, such as fibrillation, fissure, and partial-or full-thickness defect, which are preceded by biochemical changes of ECM [12,33,34]. Even though it can clearly show morphological abnormalities of cartilage, and has relatively The FCC of the bilateral knees was observed to be smooth and with little change in morphology high detection rates for the late cartilage degeneration, which mostly indicate the progressed stage of OA with irreversible variation of the cartilage, it was limited in detecting ultrastructural changes associated with early cartilage degeneration [2,35,36]. In this study, T 2weighted images were used to assess whether the morphologic alterations had occurred in the FCC of rabbit models or not, and were taken as a representative of conventional MRI to contrast with T1ρ images on the capability for characterizing the micro-changes in the prophase of cartilage degeneration. The histological examination demonstrated that early cartilage degeneration had already existed in the experimental sides of Group A-D, however, T 2 WI showed there was no obvious imaging change of the cartilage, and T 1ρ was highly sensitive to reflecting the biochemical micro-changes due to the early stages of FCC degeneration. It had suggested that T 1ρ MRI was a physiological imaging technique for detecting the biochemical state and micro information of tissue at the molecular level, and can identify early lesions which cannot be identified by routine MRI sequences. By mean of an animal model, this study aimed to measure the T 1ρ value of FCC for quantitative reflection of PG depletion, which was the earliest finding of early cartilage degeneration. Toward this end, we successfully established a rabbit model of PG degradation inducing FCC degeneration. Histodiagnosis verified that the model highly replicated the pathological process of early articular cartilage degeneration, and Group A, B, C, and D were at different stages of early degeneration, which ensured the reliability of the results. In our experiment, to select the rabbit as a model was for the reason that there was a high degree of similarity between the biochemical properties of articular cartilage in a rabbit and a human, such as the main composition of ECM, the morphology and alignment of chondrocytes, and the pathological mechanism of cartilage degeneration in rabbits also resembled human being [37,38]. We chose to focus on the FCC mostly because that was readily identified by using surface coil for generating images, and easy to dissect from the knee joint thereby histological examination.
Histological sections showed that there was a visible reduction in staining intensity of PG in papain-treated cartilage as compared with control cartilage, and the red of staining in Group A-D gradually faded, while the normal morphology of chondrocytes gradually disappeared with the emergence of fibrosis in ECM over time. All of that indicated the degree of early cartilage degeneration progressively increased from Group A to Group D, and PG depletion was the characteristics of early cartilage degeneration, its loss was proportional to the degree of cartilage degeneration. With comprehensive analysis of T 1ρ imaging, a strong correlation was revealed between changes in PG content and T 1ρ value. The result of this study showed a obviously difference in T 1ρ values of FCC between experimental side and control side, the T 1ρ values on the experimental side were significantly higher than the control side. Moreover, there was a trend of gradual increase in T 1ρ values with the passage of postoperation time (after the last papain injection). This result was consistent with the loss of PG determined by the histodiagnosis, which suggested that a decreased PG content induced the increase of T 1ρ value, and the increased T 1ρ value reflected the decreased PG content. However, there was no significant difference in T 1ρ values between the control sides of Group A-C and the blank group (Group E), and histological examination showed that there was no cartilage degeneration in the control sides (except Group D) and Group E, which indicated the scientific validity of the case-control design and the objective presence of T 1ρ relaxation time extension in the experimental-side cartilage. Therefore, it can be clearly demonstrated based on the above analyses that the results of the present study were supported by the histological evidence, the increased T 1ρ value suggests early cartilage degeneration, and was positively correlated with the degree of cartilage degeneration. In addition, T 1ρ values increased in the control sides of Group D, and the phenomenon was consistent with the histological result. We speculated that this result was caused by the compensatory overload in the control-side cartilage due to cartilage damage of the experimental side, and such speculations surely warrants further investigations. Furthermore, we also found that 1 week after the last injection of papain, there was a significant difference in T 1ρ values of Group A between the experimental side and control side prior to cartilage degeneration shown by histological examination. This showed that the sensitivity of T 1ρ imaging may exceed that of histological examination for ultra-early low-level cartilage lesions, however, which should be verified by further research.

Conclusions
In conclusion, the main contribution of this study lay in the investigation of T 1ρ imaging to quantitatively reflect biochemical changes in the initial phase of cartilage degeneration. The results suggested that T 1ρ value was a sensitive indicator for accurately detecting early articular cartilage degeneration, its increasing was closely related to the degenerative change. Therefore, T 1ρ MRI could be potentially used for the prospective and quantitative evaluation of early cartilage degeneration, and might allow intervention of cartilage degeneration at the earliest stages. However, T 1ρ technique was still available with only a few types of MR scanners (minimum usable field strength of 3.0 T) with specialized modified pulse sequences on a research basic, and it had not been widely applied [2,25,26]. The present study is also a preliminary exploration of the T 1ρ imaging in the evaluation of articular cartilage, and there were some difficulties with the T 1ρ technique, such as the repeatability of T 1ρ imaging, the setting of scanning parameters, the design of the T 1ρ image reconstruction program and the quality of the imaging. But despite all this, with the improvement of T 1ρ technique and optimization of the scanning sequence, T 1ρ MRI may be a useful tool for detecting early cartilage degeneration in clinical practice.