Ethical statement
This study was approved by the ethics committee of Keio University School of Medicine (approval number 20150050), and written informed consent was obtained from all participants. Experiments were conducted in accordance with the Ministry of Education, Culture, Sports, Science and Technology of Japan, and Ministry of Health, Labor and Welfare of Japan guidelines.
Subjects
The present study, which examines the 10-year longitudinal changes of the thoracic paravertebral muscles, was conducted as part of a 20-year longitudinal MRI study of the cervical spine.8–12 (Fig. 1). Whole-spine MRI images were obtained from a total of 223 asymptomatic subjects during the period from 2005 to 2007 in the initial study. For this study, we contacted the original subjects by mail and requested their participation in a 10-year follow-up study taking place between 2015 and 2017. Among the 223 original subjects, 103 agreed to participate in the present study (follow-up rate of 46.2%), and MRI as well as physical examinations were conducted at eight of the original 11 participating institutions. The omission of required imaging sequences or the lack of image quality that allowed for clear visualization of the borders of the posterior extensor muscle led to the exclusion of 16 cases, leaving 85 cases as the subjects of this study (final follow-up rate of 38.1%) (Fig. 1).
There were 50 male and 35 female subjects, with a mean age of 44.8 ± 11.5 years at the initial study conducted 10 years ago (distribution shown in Table 1) and a mean interval of 9.9 ± 0.8 years between the MRI studies. All subjects filled questionnaires related to clinical symptoms and lifestyle habits and underwent physical examinations by spine surgeons.
MRI Examination
In the initial study, MRI were acquired using a 1.5-Tesla (T) (SIGNA Excite HD 1.5 T, General Electronic, WI, USA) superconducting MRI scanner, using the imaging protocol detailed in our earlier report.9 In the present study, at one principal institution, where about half of the subjects in this study under- went MRI, the MR images were taken by a fast spin-echo technique using a 1.5 T superconducting scanner and phased array coils (Excelart Vantage, TOSHIBA) with the following sequence: a T2-weighted axial images (TR/ TE, 3600/120; ETL, 23; 4-mm slice thickness; FOV 16 cm; matrix size, 256 × 192; number of excitations NEX, 2). At the other institutions, basically the same protocol was used as at the principal institution, i.e., a fast spin-echo sequence was employed. T2-weighted axial images parallel to each vertebral discs were analyzed in this study.
Measurements of Paraspinal Muscle Cross-sectional Area
The CSAs of the transversospinalis muscles (rotatores, multifidus, semispinalis), erector spinae muscles (spinalis, longissimus, iliocostalis), and total extensor muscles at each thoracic intervertebral level were measured on T2-weighted axial images after scale calibration. Using Image J 1.52a, a Java-based version of the public domain NIH software, the fascial border of each muscle was manually traced (Fig. 2), yielding the area of the muscle cross section. Each measurement was taken twice, and the values were averaged.
The measurements for each thoracic intervertebral level from T1/2 to T11/12 were compared between the initial and 10-year follow-up MRI studies. The change in muscle CSAs was calculated as follows:
CSAs change (%) = (CSAs at follow-up − CSAs at initial study) / CSAs at initial study x 100.
Evaluation of Fat Infiltration of the Muscle
The degree of fat infiltration was evaluated by the luminance of the cross section of each extensor muscle group at each thoracic intervertebral level, using the luminance of pure muscle tissue as reference. Since the luminance of fat tissue is higher than muscle, increased fat infiltration of muscle will be reflected as increased luminance. Fat infiltration rate was calculated as follows:
Fat infiltration rate (%) = luminance in CSAs of whole muscle / luminance in pure muscle luminance x 100.15,16
The changes in the fat infiltration rate were compared between the initial and follow-up study.
Evaluation of Intervertebral Disc Degeneration
Degeneration of the intervertebral disc was evaluated based on the following five findings on MRI: (1) decrease in the signal intensity of the intervertebral discs, (2) anterior compression of the dura and the spinal cord, (3) posterior disc protrusion, (4) disc-space narrowing, and (5) foraminal stenosis. MR findings were scored with several modifications in the Matsumoto classification8 in Okada's study9. The scoring system for the different MR findings is shown in Supplement 1. An increase of at least 1 grade in any intervertebral level was considered to be the progress of degeneration.
Two experienced neuroradiologists blinded to the subjects independently read the MRI images and graded the intervertebral discs. The final results used the findings of one of the two neuroradiologists. The interobserver reliability was good to excellent, as reported previously.11
Statistical Analysis
Student's t-tests were used to compare the mean CSAs of the study muscles in the initial study with those in the 10-year follow-up study. A Poisson regression analysis was applied to evaluate the associations of age, sex, medical history, lifestyle, and disc degeneration with the changes in posterior extensor muscle CSAs and luminance.
A P-value of < 0.05 was considered to be indicative of statistical significance. All statistical analyses were performed using Dr. SPSS Ⅱ for Windows (SPSS Japan Inc., Tokyo, Japan) and Stata15 for Windows (Stata Corporation, College Station, TX, USA).
For the interobserver error assessment, 20 randomly chosen MR images were measured independently by the first investigator (H.U.) and second investigator (H.I.). For the intraobserver error assessment, 20 images were randomly selected, and the images were measured two times with an interval of 3 months. The values measured by the first investigator (H.U.) were used for this study.