Study participants
Patients were prospectively enrolled from the general and orthopedic outpatient clinics of the 960th Hospital of the Joint Logistics Support Force of the PLA between January 2023 and July 2023.
Approval for the study was granted by the Institutional Research Board of the 960th Hospital of the Joint Logistics Support Force of the PLA under Grant No. (2023) Research Ethics Review No. (020). Written informed consent was obtained from all study participants, and the experiments adhered to relevant guidelines and regulations.
The inclusion criteria comprised patients with discordant knee pain lasting ≥ 3 months in the past 12 months or pain for most days in the past month, along with definite radiographic evidence of osteoarthritis in at least one compartment of at least one knee, and an age exceeding 40 years. Exclusion criteria included individuals with a prior history of other forms of joint disease, such as bone implants, osteonecrosis, or CT scans revealing bone cysts. Patients taking medications affecting bone metabolism (e.g., bisphosphonates) were also excluded.
Participant assessment
Body mass index (BMI) was calculated by dividing weight in kilograms by height in meters squared. To assess pain intensity at the affected knee joint, the pain section of the Western Ontario McMasters Osteoarthritis Index (WOMAC) was utilized [15]. This involved measuring pain on a 10-point scale (ranging from 0 for no pain to 9 for severe pain) across five categories: walking on a flat surface, climbing or descending stairs, pain at night, sitting or lying, and standing upright (representing the most severe pain). The individual element pain scores were then summed to derive a WOMAC pain score within the range of 0 to 45. During knee X-ray scans, the knee joint was observed in a non-flexed position for standing images, while a flexed position of 135 degrees was applied for lateral images. The severity of knee osteoarthritis (KOA) was determined using Kellgren-Lawrence (KL) grading [16], with grading ranging from 0 to 4. Knee alignment was categorized as varus, valgus, or neutral based on the angle between the femoral and tibial axes. Varus alignment was defined as an angle less than 176 degrees, valgus alignment as an angle greater than 180 degrees, and neutral alignment as 176–180 degrees[7]. The diagnosis of cysts relies on CT image reconstructions using bone and soft tissue windows, where cysts are identified as well-defined, round or oval, low-density lesions with a diameter exceeding two voxels. KL grading and cyst diagnosis were conducted using a double-blind method by two diagnostic physicians, each possessing over ten years of work experience.
Participant classification
A total of 53 patients diagnosed with knee osteoarthritis (KOA) participated in the study, with 11 patients excluded due to the presence of subchondral cysts revealed in CT scans. The final cohort comprised 42 patients, consisting of 14 men and 28 women, with a mean age of 55.7 ± 9.0 years and a BMI of 25.9 ± 3.0. This group contributed data for a comprehensive analysis of 84 knees. Relevant participant characteristics, including Kellgren-Lawrence (KL) grading and knee alignment in weight-bearing situations, are presented in Table 1.
Patients received separate WOMAC scores for their right and left knee joints based on their perception of knee pain. Subsequently, patients' right and left knee joints were divided into two groups based on their respective WOMAC scores. The "moderate-severe pain" group was characterized by a relatively high WOMAC score, while the "mild pain" group had a relatively low WOMAC score. A control group was established, comprising participants with no pain according to WOMAC and a Kellgren-Lawrence (KL) grade of 0.
CT scan parameters
CT scanning parameters included a tube voltage of 120 kVp, automatic tube current modulation, axial scanning plane, a slice thickness of 0.625 mm, and an in-plane pixel size of 0.625 mm × 0.625 mm. The effective radiation dose per scan was estimated to be less than 0.16 mSv using software (CT-DOSE, Neusoft Medical System, Shenyang, China), corresponding to a natural background radiation exposure of fewer than 16 days in a comparable time frame [17].
Quantitative CT (QCT) acquisition
A 256-slice spiral CT instrument (NeuViz Glory, Neusoft Medical, Shenyang, China) was employed for bone imaging. Each participant underwent CT scans with a solid QCT reference phantom (FT/HK-2000, Huake Testing New Technology Development Institute, Chengdu, China) beneath them. Phantoms facilitated the conversion of grayscale CT Hounsfield Units (HU) to apparent BMD (mg/cm³ K₂HPO₄). Previous human studies validated that QCT density measurements accurately reflect true BMD. Simultaneous CT scans of both legs were performed with the participant in a supine position, focusing on the distal femur, patella, and proximal tibia. However, only the proximal tibia was analyzed in this study.
Equivalent volume BMD values are determined by converting grayscale Hounsfield units (HU) to BMD, a linearly proportional measure to the CT value of bone tissue and its BMD on the same CT image frame. This relationship is expressed as:
$$\:\text{B}\text{M}\text{D}=\frac{{\text{H}\bullet\:({\mu\:}_{x}/{\mu\:}_{w})\bullet\:(\text{k}\text{V}\text{p}/\text{V})}^{\text{n}}}{\text{C}\text{T}\:\text{v}\text{a}\text{l}\text{u}\text{e}\:\text{i}\text{n}\text{d}\text{e}\text{x}\text{i}\text{n}\text{g}\:\text{f}\text{a}\text{c}\text{t}\text{o}\text{r}\bullet\:{\rho\:}_{w}\bullet\:{\mu\:}_{x}\bullet\:{\rho\:}_{x}}$$
Here, H is the CT value in Hu, \(\:{\mu\:}_{x}\) and \(\:{\mu\:}_{w}\) are the linear attenuation coefficients of the measured tissue and water, respectively in \(\:{cm}^{-1}\); \(\:\text{k}\text{V}\text{p}/\text{V}\) is the ratio of the mass attenuation coefficient associated with the CT tube's voltage value (a constant for the same device), n is the specified CT value indexing factor (usually 1000) in HU. \(\:{\rho\:}_{w}\) is the density of water and \(\:{\rho\:}_{x}\) is the mass attenuation coefficient and mass density of the measured tissue, in \(\:{cm}^{2}/g\) and \(\:{g/cm}^{3}\).
To mitigate the impact of uncontrollable fluctuations in CT equipment's kV values on CT values, the body model is synchronized with bone tissue during scanning. The 3D imaging and precise indexing of CT values enable the Quantitative CT (QCT) technique to detect changes as subtle as several milligrams of BMD. This technique demonstrates a repeatability of 0.8% over 10 scans at any interface and maintains an accuracy error within ± 1.5%.
Image normalization
Image normalization procedures were conducted for each case before measurement. Utilizing AVW 2.0 software (Neusoft Medical, Shenyang, China) for postprocessing, a bone window was selected for Multi-Planar Reconstruction (MPR). CT image postprocessing at the workstation involved 3D BMD measurement. In the sagittal position, the localization line's longitudinal axis aligned parallel to the knee joint's long axis. For coronal positioning, the line connecting the lowest points of the inner and outer condyles of the tibial plateau was identified, and the line was placed in proximity to the upper edge of the tibial cortex. In the transverse position, this layer was defined as the subchondral level.
Regional analysis
Delineation of the medial and lateral tibial compartments: To delineate the medial and lateral tibial compartments, a continuous line is drawn at the widest point of the tibial plateau. The length of this line is measured, and it is then divided into three segments in proportions of 40%, 20%, and 40%. Subsequently, a plumb line is drawn at the two nodes representing 40%, creating a clear demarcation. The areas on either side of the plumb line, extending to the edge of the tibia, are defined as the medial and lateral tibial compartments.
Delineation of the anterior and posterior compartments: The anterior and posterior compartments of the tibia are delineated by dividing the anterior and posterior dimensions of each plateau into two subregions with equal spacing. These subregions are labeled as L1, L2, L3 and L4, as illustrated in Fig. 1.
Region of Interest (ROI) Definition: The ROI is manually drawn by modifying the ROI boundary of the medial and lateral platforms while ensuring the exclusion of cortical bone. This meticulous delineation process, conducted by operators who have received proper training to identify and distinguish different types of bone tissue, allows for accurate assessment and analysis within the designated compartments. Following the definition of the ROI, a second check by an experienced expert is performed to ensure that no cortical bone has been mistakenly included.
All assessments were performed at four subchondral surface depths, each corresponding to specific anatomical structures: (1) 0-2.5 mm: This depth aligns closely with the density of the subchondral endplate and cortical bone [18]; (2) 2. 5–5.0 mm: Proximity to the density of the subchondral trabeculae characterizes this depth [10]; (3) 5.0-7.5 mm: This range corresponds to the density of the area between the subchondral trabecular and the proximal epiphyseal line trabeculae.. (4) 7.5–10.0 mm: This depth mirrors the density of the proximal epiphyseal plate trabeculae.
Statistical analysis
BMD was determined by calculating the mean from three CT scans. Quantitative variables were summarized using descriptive statistics (mean ± SD). To affirm the correlation between body mass index (BMI) and Western Ontario McMasters Osteoarthritis Index (WOMAC), the Spearman correlation coefficient was utilized. Additionally, a chi-square test was conducted to verify the correlation between Kellgren-Lawrence (KL) grading and WOMAC. Exploring the variability in BMD between the "moderate-severe pain" group and the "mild pain" group involved the application of a paired t-test. Statistical analysis was carried out using SPSS 25 software (SPSS Inc.), and statistical significance was defined at p < 0.05.