Relationship Between Endplate Defects, Modic Change, Facet Joint Degeneration, and Disc Degeneration of Cervical Spine.

OBJECTIVE: The ''disc degeneration precedes facet joint osteoarthritis'' hypothesis and multidimensional analysis were actively discussed in lumbar spine. However, in cervical spine degeneration, the multifactorial analyzes of disc degeneration (DD), Modic changes (Mcs), facet degeneration, and endplate degeneration (ED) is still limited. In this cross-sectional study, we aimed to analyze the prevalence and interrelationship of cervical DD parameters.
METHODS: We retrospectively recruited 62 patients aged between 60 and 70 years. The disc height, segmental angle, ossified posterior longitudinal ligament (OPLL), ED, facet joint degeneration (FD), uncovertebral joint degeneration (UD), DD, spinal stenosis (SS), Mc, and cord signal change (CS) were evaluated using a previously well-known grading system.
RESULTS: The prevalence of cervical degenerative parameters were DD (grade 1, 1.2%; grade 2, 13.3%; grade 3, 54.8%; grade 4, 19.0%; grade 5, 11.7%), OPLL (26.2%), SS (grade 0, 7.7%; grade 1, 42.3%; grade 2, 26.2%; grade 3, 23.8%), UD (39.1%), ED (normal, 69.0%; focal defect, 9.7%; corner defect, 11.7%; erosion, 6.9%; sclerosis, 2.8%), and FD (normal, 48.8%; narrowing, 27.0%; hypertrophied, 24.2%). The interrelationship of degenerative parameters showed close relation between UD, SS, DD, OPLL, Mc. ED, and CS has partial relation with degenerative finding. FD only has relation with UD, and Mc.
CONCLUSION: Our results may indicate that FD is a degeneration that occurs independently, rather than as a result of other degenerative factors.

INTRODUCTION

Neck pain is a multidimensional phenomenon that affects as much as two-thirds or more of the general population [1]. Among the many causes of neck pain, degeneration of the cervical disc is a common condition [2]. However, the veritable cause of the cervical disc degeneration (DD) remains unclear. In analogy to this, the main opinion on the initial process of spinal degeneration is that degeneration begins in the intervertebral disc and facet joint degeneration (FD) occurs following DD [3-5]. However, other authors point out that DD does not necessarily antecedent of facet degeneration [6-8]. In addition, in the study of lumbar intervertebral DD, a hypothesis about endplate degeneration (ED) driven DD was also presented [9-11], Furthermore, the multifactorial analyzes of DD, Modic changes (Mcs), facet degeneration, ED, and pain are also recently evaluated [12,13]. However, this approach is still limited in cervical spine degeneration studies. In particular, the cervical spine differs from the lumbar spine, such as uncovertebral joint, ossified posterior longitudinal ligament (OPLL), and a higher range of motion (ROM). In this regard, multifactor analysis thereof is necessary for the cervical spine.

In this cross-sectional study, we aimed to analyze the prevalence and interrelationship of segmental ROM, facet degeneration, endplate defect, OPLL, uncovertebral joint degeneration (UD), Mc, and cervical DD.

MATERIALS AND METHODS

1. Patient Population

We retrospectively recruited patients aged 60 and 70 years, who underwent all C-spine X-ray series (upright lateral images on neutral, flexion, and extension position), conventional magnetic resonance imaging (MRI), and computed tomography (CT) at our hospital, between January 2018 to December 2018. All patients with a history of fracture, tumor, infectious spondylitis, deformity, or previous cervical operation were excluded. Finally, 62 patients were included. (average age, 63.60±2.94; 39 males [62.9%] and 23 females [37.1%]) A total 62 patients, 310 segments (from C3/4 to C7/T1) were evaluated, but C7/T1 has an invisibility endplate margin (16/62 segments, 25.8%), so, we only excluded this level, final enrolled segments were 244, from C3/4 to C6/7.

2. Radiologic Measurements (Table 1)

All from C3/4 to C7/T1 level were measured following degeneration-related parameters.

1) C-spine X-ray series

The segmental angle (SA) of each level from C3–4 to C7–T1 at neutral, flexion, and extension position, and disc height of the mid vertebral body at the neutral position were measured. Segmental ROM (SA ROM) was the difference of SA between flexion and extension (Fig. 1)

Fig. 1.

Measurement methods for X-ray parameters: (A) segmental angle at flexion position, (B) disc height and segmental angle at neutral position, and (C) segmental angle at extension position. SA, segmental angle.

2) CT of C-spine

Using sagittal images, we measured OPLL (present or not) and OPLL type (segmental, continuous, circumscribed, and mixed) [12,14-17], and ED (normal, focal defect, corner defect, erosion, and sclerosis) [18]. Using coronal images, UD was evaluated (normal, hypertrophied). On axial images, FD (normal, narrowing, and hypertrophied) [19,20], and facet tropism (bilateral orientation difference > 10°) [21] were measured.

3) MRI of C-spine

Using sagittal T2-weight sequences, we measured DD grade (using Pfirrmann grade—grade 1: disc is homogeneous with bright hyper-intense white signal intensity [SI] and normal disc height, nucleus and annulus are clearly differentiated, and disc height is normal; grade 2: disc is inhomogeneous, but keeping the hyper-intense white signal, and a grey horizontal band could be present; grade 3: disc is inhomogeneous with an intermittent grey SI, distinction between nucleus and annulus is unclear, disc height is normal or slightly decreased; grade 4: disc is inhomogeneous with a hypo-intense dark grey SI, there is no more distinction between the nucleus and annulus, disc height is slightly or moderately decreased; grade 5: disc is inhomogeneous with a hypo-intense black SI, and the disc space is collapsed) [22], spinal stenosis grade (SS; using Muhle grade—grade 0: normal, grade 1: partial obliteration of the anterior or posterior subarachnoid space or of both, grade 2: complete obliteration, grade 3: complete obliteration with cervical cord compression of displacement) [23], Mc (only T2 sequence was measured), and cervical spinal cord signal change (CS: only T2 sequence was measured).

3. Statistical Analysis

All statistical analyses were evaluated using IBM SPSS Statistics ver. 21.0 (IBM Co., Armonk, NY, USA). Some ordinary parameters were reclassified into binary categories according to degeneration or not; SS: grade 0+1 (control) vs. grade 2+3 (degeneration), ED: normal vs focal defect, corner defect, erosion, sclerosis (degeneration), FD: grade 1 (control) vs. grade 2+3 (degeneration), DD: grade 1+2 (relative lower degeneration group) vs. grade 4+5 (relativ e high degeneration group). In Table 2, generalized estimating equation was performed to reduce within-person bias that measured the binary degenerative variable incidence (degeneration or not) for multiple levels (from C3/4 to C6/7). In Tables 3 to 5, it was assumed that the degeneration degree of each level was patient-independent and level-independent, and chi-square or Fisher exact test was used. For continuous variables, presented as mean±standard deviation, and the independent t-test or 1-way analysis of variance (ANOVA) was adopted to assess the difference of each degenerative parameters, or levels. For categorical or ordinal variables, the Spearman correlation analysis was performed to detect the interrelationship between degenerative parameters.

Table 2.

Degenerative change according to level

Variable	C3–4	C4–5	C5–6	C6–7	p-value[†]
OPLL	12 (19.4)	23 (37.1)	21 (33.9)	9 (14.5)	< 0.001
No	50/62 (80.6)	39/62 (62.9)	41/62 (66.1)	53/62 (85.5)
Segmental	11/62 (17.7)	17/62 (27.4)	17/62 (27.4)	7/62 (11.3)
Mixed	1/62 (1.6)	3/62 (4.8)	3/62 (4.8)	2/62 (3.2)
Continuous	0 (0)	0 (0)	0 (0)	0 (0)
Circumscribed	0 (0)	3/62 (4.8)	1/62 (1.6)	0 (0)
ED					0.010[*]
Normal	47 (75.8)	42 (67.7)	34 (54.8)	48 (77.4)
Focal defect	3 (4.8)	4 (6.5)	10 (16.1)	7 (11.3)
Corner defect	6 (9.7)	11 (17.7)	9 (14.5)	3 (4.8)
Erosion	2 (3.2)	4 (6.5)	8 (12.9)	3 (4.8)
Sclerosis	4 (6.5)	1 (1.6)	1 (1.6)	1 (1.6)
UD	13 (21.0)	18 (29.0)	36 (58.1)	30 (48.4)	< 0.001[*]
FD					0.591
Normal	32 (51.6)	32 (51.6)	27 (43.5)	30 (48.4)
Narrowed	15 (24.2)	15 (19.4)	20 (32.3)	20 (32.3)
Hypertrophied	15 (24.2)	18 (29.0)	15 (24.2)	12 (19.4)
DD					< 0.001[*]
Grade 1	1 (1.6)	1 (1.6)	1 (1.6)	0 (0)
Grade 2	7 (11.3)	3 (4.8)	6 (9.7)	17 (27.4)
Grade 3	37 (59.7)	38 (61.3)	33 (53.2)	28 (45.2)
Grade 4	13 (21.0)	15 (24.2)	6 (9.7)	13 (21.0)
Grade 5	4 (6.5)	5 (8.1)	16 (25.8)	4 (6.5)
Mc	15 (24.2)	15 (24.2)	17 (27.4)	12 (19.4)	0.674
CS	8 (12.9)	13 (21.0)	12 (19.4)	4 (6.5)	0.058
SS					< 0.001[*]
Grade 0	11 (17.7)	4 (6.5)	1 (1.6)	3 (4.8)
Grade 1	28 (45.2)	27 (43.5)	20 (32.3)	30 (48.4)
Grade 2	14 (22.6)	11 (17.7)	19 (30.6)	21 (33.9)
Grade 3	9 (14.5)	20 (32.3)	22 (35.5)	8 (12.9)

Values are presented as number (%).

OPLL, ossified posterior longitudinal ligament; ED, endplate degeneration; UD, uncovertebral joint degeneration; FD, facet joint degeneration; DD, disc degeneration; Mc, Modic change; CS, cord signal change; SS, spinal stenosis.

p < 0.05, statistical significance.

p-value was the results of generalized estimating equation for binary degenerative variable (degeneration or not).

Table 3.

Interrelationship between polynomial degenerative parameters

Variable	UD	DD	OPLL	Mc	SS	ED	CS	FD
Spearman correlation test
DD (1–5)	< 0.001[*]		0.003[*]	0.030[*]	< 0.001[*]	0.009[*]	< 0.001[*]	0.803
SS (0–3)	< 0.001[*]	< 0.001[*]	< 0.001[*]	0.006[*]		0.126	< 0.001[*]	0.624
ED (0–4)	0.002[*]	0.009[*]	0.053	< 0.001[*]	0.126		0.796	0.736
FD (1-3)	0.108	0.803	0.505	0.004[*]	0.624	0.918	0.736

UD, uncovertebral joint degeneration; DD, disc degeneration; OPLL, ossified posterior longitudinal ligament; Mc, Modic change; SS, spinal stenosis; ED, endplate degeneration; CS, cord signal change; FD, facet joint degeneration.

p < 0.05, statistical significance.

Table 5.

Disc height, segmental angle (SA), segmental range of motion (ROM), and degenerative parameters

Parameter	Category	Disc height	p-value	SA	p-value	SA ROM	p-value
Sex	Male	4.88 ± 1.18	0.001[*]	3.70 ± 3.14	0.048[*]	2.23 ± 4.67	0.402
	Female	4.40 ± 1.03		2.94 ± 2.36		1.69 ± 5.19
OPLL	No	4.75 ± 1.10	0.301	3.70 ± 2.96	0.009[*]	2.08 ± 4.77	0.785
	Yes	4.58 ± 1.28		2.60 ± 2.54		1.89 ± 5.17
UD	No	4.87 ± 1.13	0.001[*]	3.53 ± 3.01	0.360	2.17 ± 4.96	0.518
	Yes	4.33 ± 1.12		2.16 ± 2.63		1.73 ± 4.68
Mc	No	4.85 ± 1.07	< 0.001[*]	3.59 ± 2.93	0.107	2.17 ± 5.00	0.432
	Yes	4.24 ± 1.29		2.89 ± 2.73		1.59 ± 4.45
CS	No	4.75 ± 1.15	0.132	3.46 ± 2.83	0.586	2.05 ± 4.89	0.926
	Yes	4.44 ± 1.16		3.18 ± 3.27		1.96 ± 4.84
ED	Gr 0	4.84 ± 1.12	0.010[†]	3.84 ± 3.17	0.014	2.30 ± 5.17	0.189
	Gr 1	4.31 ± 0.94		2.70 ± 1.68		1.19 ± 2.94
	Gr 2	4.79 ± 1.22		2.16 ± 1.73		0.38 ± 5.21
	Gr 3	4.00 ± 1.24		2.93 ± 2.53		3.49 ± 3.21
	Gr 4	4.17 ± 1.34		2.13 ± 1.31		1.74 ± 2.98
FD	Gr 0	4.82 ± 1.14	0.297	3.17 ± 2.93	0.240	1.74 ± 5.14	0.530
	Gr 1	4.59 ± 1.14		3.40 ± 2.68		2.04 ± 3.85
	Gr 2	4.60 ± 1.18		3.94 ± 3.03		2.62 ± 5.33
DD	Gr 1	4.99 ± 1.74	< 0.001[†]	3.73 ± 2.01	0.031	3.23 ± 5.08	0.146
	Gr 2	4.86 ± 1.01		4.42 ± 3.59		2.64 ± 5.22
	Gr 3	4.87 ± 1.15		3.32 ± 2.86		2.04 ± 4.92
	Gr 4	4.78 ± 1.07		3.80 ± 2.92		2.79 ± 4.80
	Gr 5	3.59 ± 0.75		2.14 ± 1.64		0.14 ± 3.99
SS	Gr 0	5.07 ± 0.87	0.035	4.19 ± 2.63	0.470	4.92 ± 5.31	0.005[†]
	Gr 1	4.84 ± 1.02		3.57 ± 3.10		2.56 ± 5.01
	Gr 2	4.68 ± 1.25		3.19 ± 2.61		1.50 ± 4.37
	Gr 3	4.37 ± 1.28		3.15 ± 2.90		0.75 ± 4.57

Values are presented as mean±standard deviation.

OPLL, ossified posterior longitudinal ligament; UD, uncovertebral joint degeneration; Mc, Modic change; CS, cord signal change; ED, endplate degeneration; Gr, grade; FD, facet joint degeneration; DD, disc degeneration; SS, spinal stenosis.

p < 0.05, statistical significance for independent t-test.

p < 0.05/number of categorical variables for analysis of variance test.

4. Ethics Statement

The study protocol was approved by the Institution Review Board (IRB) of the affiliated hospital, which waived the requirement for informed consent due to the retrospective nature of this study (IRB No. I-2020-249).

RESULTS

1. Prevalence of Degenerative Change

1) CT findings

OPLL (26.2%), OPLL type (normal, 74.6%; segmental, 20.2%; continuous, 0%; circumscribed, 1.6%; mixed, 3.6%), ED (normal, 69.0%; focal defect, 9.7%; corner defect, 11.7%; erosion; 6.9% and sclerosis 2.8%), UD (39.1%), FD (normal, 48.8%; narrowing 27.0%; hypertrophied, 24.2%). Only 8 segments (3.2%) had a facet tropism.

2) MR findings

DD (grade 1, 1.2%; grade 2, 13.3%; grade 3, 54.8%; grade 4, 19.0%; grade 5, 11.7%), SS (grade 0, 7.7%; grade 1, 42.3%; grade 2, 26.2%; grade 3, 23.8%), Mc (23.8%), and CS (14.9%).

2. Degenerative Change According to Level

The prevalence of cervical degenerative parameter was summarized in Table 2. There was a difference in the prevalence of degenerative parameter according to level, except FD, Mc, and CS.

3. Interrelationship Between Degenerative Parameters

1) The Spearman correlation test

For the multiple or ordinary parameters (ED, DD, SS, and FD), the interrelationship between degenerative change was summarized in Table 3. DD had interrelationship with all degenerative parameter but FD. SS had interrelationship with all degenerative parameter but ED. ED related with UD, DD, OPLL, but had no relation with SS, CS, and FD. FD only had a relation with UD and Mc.

2) Interrelationship between degenerative parameters (binary analysis)

The interrelationship between binary degenerative parameter was summarized in Table 4. UD related to all degenerative parameters. Other parameters showed partial interrelationship with each other. However, FD had only 2 related factors: Mc and UD. Compared with the result of K independent samples test, DD had relatively low interrelationship compared with other parameters because of lower enrolled segments (grade 3 degeneration: 136/248 segments, 54.8% was excluded in binary analysis).

Table 4.

Interrelationship between binary degenerative parameters

Variable	UD	DD	OPLL	Mc	SS	ED	CS	FD
Binary categorical analysis
UD		OR = 3.046	1.912	2.720	2.367	2.339	3.042	1.767
		p = 0.014[*]	0.027[*]	0.001[*]	0.002[*]	0.003[*]	0.003[*]	0.042[*]
DD			6.787	1.716	2.404	2.411	5.667	1.171
			0.001[*]	0.276	0.040[*]	0.060	0.018[*]	0.840
OPLL				2.475	3.975	2.070	3.786	0.898
				0.006	< 0.001[*]	0.019[*]	< 0.001[*]	0.773
Mc					2.167	4.618	1.225	2.468
					0.017[*]	< 0.001[*]	0.676	0.004[*]
SS						1.302	15.237	1.524
						0.410	< 0.001[*]	0.127
ED							0.930	0.898
							> 0.999	0.784
CS								1.144
								0.725
FD

OPLL, ossified posterior longitudinal ligament; ED, endplate degeneration; UD, uncovertebral joint degeneration; FD, facet joint degeneration; DD, disc degeneration; Mc, Modic change; CS, cord signal change; SS, spinal stenosis.

Odd ratio (OR) and p-value were calculated by pairwise 2×2 contingency table with converting some polynomial factors into binomial factors.

Binary variable conversion: DD (no: grade 1+2 vs. yes: grade 4+5), SS (no: grade 0+1 vs. yes: grade 2+3), ED (no: normal vs. yes: defect, erosion, and sclerosis), FD (no: grade 0 vs. yes: grade 1+2).

p < 0.05, statistical significance.

4. Disc Height, SA, Segmental ROM, and Degenerative Parameters

Male patients had a higher disc height, and SA than female (p = 0.001 and p = 0.048, respectively). There was a difference in disc height, SA, or SA ROM according to degenerative parameters except for CS and FD (Fig. 2, Table 5).

Fig. 2.

Analysis of variance study for disc height, segmental angle (SA), SA range of motion (SA ROM), and spinal stenosis, disc degeneration, facet joint degeneration, and endplate degeneration grade.

DISCUSSION

Degenerative changes at the lower cervical spine can affect different anatomical structures: the intervertebral disc, uncovertebral joints, and facet joints [24]. Cervical radiculopathy is frequently secondary to osteoarthritic changes of the uncovertebral joints anteriorly or facet joints posteriorly rather than DD [25]. In this analogy to this, to understand the interrelationship between cervical degenerative parameters is important. There was some study about the relationship between age [26], FD [20], facet orientation/tropism [21], Mc [12,27,28], and DD, but there was no study including ED, OPLL, UD, and SS.

The cervical FD was classified by several authors. Rydman et al. [16] study found that the facet narrowing has relatively high reliability than osteophyte and articular surface irregularities (kappa value, 0.569, 0.313, 0.365, respectively). In that sense, we adopted the method used by Park et al. for the FD grade [20]. Grade 4, a bony fusion of the facet joint, was not identified in our study. Kettler et al. [26] classified the degree of FD through cadaveric studies. They also report that FD increases with age and evenly distributed at all levels. Kim et al. [29] evaluated the prevalence of asymptomatic cervical FD using CT imaging. They reported 33% of prevalence, and C6/7 were the most affected level. They also agreed on the age-related prevalence of cervical FD. As in the previous 2 studies, age has high relevance in FD.

To exclude the effects of age, and to focus on interpreting the interrelationship between cervical degenerative parameters, we only included patients with age from 60 to 70. Walraevens also pointed out the lack of evidence about the relationship between the DD and facet degeneration in the cervical spine, which has been thoroughly investigated in lumbar spine degeneration [30-33]. They reported a spatial weak correlation between FD and DD [34]. In the current study, FD had a relationship with UD and Mc. But FD was not related with other factors. In that sense, FD occurs relatively independently of other cervical degeneration parameters. The prevalence of degenerative factors according to level is also in line with these findings. The prevalence of cervical spine degeneration was different according to the cervical level. In particular, the prevalence of OPLL, ED, UD, grade 4 DD, and grade 2–3 SS was higher in C5–6. However, FD, Mc, and CS did not show a statistically significant difference by level. These results also could be evidence that these exceptions factors have a different mechanism from others. We can also found that these results are on the opposite side of the “disc degeneration precedes facet joint osteoarthritis” hypothesis [30,35] in the cervical spine.

Recently, there were 2 studies of multifactorial analysis for lumbar spine degeneration including ED. Lv et al. [12] demonstrated that ED would accelerate the change of facet joint orientation and facet tropism in segments with MCs. Zehra et al. [13] found the prevalence of ED (67.5% of the subjects and in 13.5% of the endplates). ED has a significant relationship with disc degenerative scores, Mcs, and posterior disc displacement (p < 0.05). The cumulative endplate defect scores, which they devised, increased with narrow facet joints (p = 0.004). They also showed clinical correlation according to the size of endplate defect. Sahoo et al. [18] also found the association between ED and symptomatic lumbar disc herniation, and the presence of ED can increase the neurological deficit, and reduce the chance of recovery with conservative management. These studies emphasized the importance of ED and endplate morphology [10]. The classification of ED of the lumbar spine was defined several authors [12,15,18], but there was no method to measure ED of the cervical spine. We revised the ED classification for the lumbar spine to the cervical spine. However, the small size of cervical vertebrae, we cannot obtain a clear view of the endplate in magnetic resonance (MR) images. We found that Sahoo et al. [18] showed that the CT image could demonstrate endplate defects similar to MR. Because we could obtain more clear view for ED in CT images, we adopted ED classification using CT images.

We tried to suggest whether our analysis is meaningful not only for categorized variables (UD, FD, etc.) but also for radiological continuity variables. In that sense, the radiological factors related to segment, disc height, SA, and segment motion were measured. In particular, disc height reduction, segmental kyphosis, and loss of segmental motion are related with degeneration. OPLL, UD, Mc, and Cs are binary variables and can present all the information by Table 5, but ordinary variables such as ED, FD, DD, and SS cannot present all the information. So, we visualized it as Fig. 2. It demonstrates the change of disc height, SA, and segmental motion according to the degeneration grade. As the degeneration grade of SS, DD, ED progresses, the disc height, SA, and SA ROM decrease, but the FD shows a different pattern. These results also support that FD has a different mechanism from other degenerative factors such as SS, DD, and ED. In the case of ED, some statistical significance was obtained in ANOVA, but there was a low agreement in linear to linear combination than other factors (Fig. 2). In the grade 3, endplate erosion, SA and SA ROM tended to be relatively high than other degenerative grades, and in the grade 2, corner defects, disc height was relatively high than other degenerative grades. We suggested that these results are due to the accuracy limitations of ED classification. Previous studies have reported poor agreement of endplate sclerosis of the cervical spine using CT images (kappa value= 0.241, p = 0.092) [16]. In our study, endplate sclerosis was founded in 1 of 3 disc grade 1 segments. This shows that there is a limitation in measuring the endplate erosion using CT images. In a future study, it is necessary to re-evaluate grade 3 through a reliability test.

Some limitation of the present study should be acknowledged. First, our enrolled patients had a radiological selective bias. Only a few patients take both CT and MRI, and they are also more likely to have ventral lesions that require surgical treatment. In that sense, the statistical interpretations were radiologically dominated with; grade 3 DD with a higher prevalence of CS up to 14.9%, as well as a predilection to contain OPLL up to 26.2%.in current study group. The patients enrolled in our retrospective study cannot reflect the interrelationship between degenerative parameters of healthy volunteer or patients with mild degeneration. Second, the binary analysis was performed except grade 3 DD. We want to present the strength of the relationship (odds ratio) between degenerative parameters using binary analysis (degeneration or not). OPLL, UD, Mc, CS are originally binary variables, and FD and ED can be classified as normal or degeneration. The classification for SS also has significant differences between grades 0, 1 and grades 2, 3, so there are no problems with conversion to binary variables. However, our enrolled patients have the most prevalence of grade 3 DD (not too advanced, 54.8%). There were limitations in the binary analysis including grade 3 DD to either group (grade 1+2+3 vs. 4+5 or grade 1+2 vs. 3+4+5). We used nondegenerated group (grade 1+2) and degeneration group (grade 4+5) except grade 3 DD for more clarity binary analysis. Third, we did not validate the ED classification using CT images. In addition, our study data reported by one reader and does not include intraclass correlation (ICC) results. In this study, our goal is not to verify the accuracy or validity of the measurement method. In the future study, the ICC test and reliability test will be required for more accurate measurement. Fourth, 2 cross-sectional studies negated the association between neck pain and radiological studies [36,37]. In this study, we did not evaluate the correlation with clinical information. Since we had to present too much factors information, there was a limit to analyzing the clinical information together. Future research will further analyze the relevance to clinical information on the basis of the result of the current study.

CONCLUSION

We confirmed a close correlation between DD, ED, SS, Mc, OPLL, and UD. However, FD was only related to Mc and UD, not to other factors. These results indicate that FD is a degeneration that occurs independently, rather than as a result of other degenerative factors, and negates the “disc degeneration precedes facet joint osteoarthritis” hypothesis in the cervical spine.

Table 1.

Measuring method for radiological parameters

Modality	Parameter	Criteria	Definition
CT	(Sagittal) OPLL	No/yes
	(Sagittal) OPLL type	Segmental
		Continuous
		Circumscribed
		Mixed
	(Sagittal) Endplate degeneration	Normal
		Focal defect	Local hollow or discontinuity on EP
		Corner defect	Anterior or posterior corner of VB
		Erosion	Extensive alteration of EP
		Sclerosis	Subchondral sclerosis
	(Coronal) Uncovertebral joint degeneration	Normal
	(Axial) Facet joint degeneration	Hypertrophied	Blunting of the apex, narrowing, and sclerosing
		Normal	Facet joint space ≥ 2 mm
		Narrow	Facet joint space < 2 mm, cyst formation
		Hypertrophied	≥ 3mm lateral margin than normal level
		Bony fusion
	(Axial) Facet joint tropism	No/Yes	Bilateral facet orientation ≥ 10°
MR	(T2 sagittal) Disc degeneration	Grade 1	(NP) homogenous
		Grade 2	(NP) Inhomogeneous +white band
		Grade 3	(NP/AF demarcation) unclear
		Grade 4	(NP T2 signal) hypo-intense
		Grade 5	(disc height) collapsed
	(T2 sagittal) Spinal stenosis	Grade 0
		Grade 1	Partial obliteration
		Grade 2	Complete obliteration
		Grade 3	Cord compression or displacement
	(T2 sagittal) Modic change	No/Yes	Only T2 image was evaluated
	(T2 sagittal) Cord signal change	No/Yes	Only T2 image was evaluated

CT, computed tomography; MR, magnetic resonance; OPLL, ossified posterior longitudinal ligament; EP, endplate; VB, vertebral body; NP, nucleus pulposus; AF, annulus fibrosus.