$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
This study was conducted in accordance with the principles of the Declaration of Helsinki. Ethical approval was obtained from the Institutional Ethics Committee of Panzhihua Hospital of Integrated Traditional Chinese and Western Medicine (The Affiliated Hospital of Panzhihua University), Panzhihua, China (Approval No. 2024ZD-S-4). As this was a retrospective study based on anonymized clinical and imaging records, the requirement for written informed consent was waived by the ethics committee. All patient data were handled confidentially and anonymized prior to analysis.
Study design and setting
This retrospective cross-sectional diagnostic accuracy study was conducted at a tertiary care center with a dedicated hemophilia treatment unit. The study followed the standards for reporting diagnostic accuracy studies (STARD) guidelines. Data were obtained from hospital records and the radiology database after both the index test, MSK-US, and the reference standard (magnetic resonance imaging [MRI]) had been performed as part of routine clinical care between January 2020 and December 2025.
The primary diagnostic endpoint was the accuracy of MSK-US for detecting hemophilic arthropathy, using MRI as the reference standard. MRI examinations were requested by treating hematologists or orthopedic specialists when clinically indicated, including persistent joint pain, recurrent hemarthrosis, suspected progression of target-joint disease, discordance between clinical findings and previous imaging results, or the need for detailed structural assessment before treatment modification. MSK-US examinations were performed as part of routine joint evaluation and monitoring within the hemophilia care program.
Given the retrospective design, there is a potential risk of selection bias; this was minimized by including a consecutive series of eligible patients identified through predefined inclusion and exclusion criteria. However, because MRI was performed according to clinical indications rather than as a systematic screening procedure, the study population may have had a higher pre-test probability of joint pathology than the general hemophilia population. This potential referral bias was considered during the interpretation of the diagnostic accuracy findings and is addressed in the limitations section.
Participants
Eligibility criteria and selection
Patients with a confirmed diagnosis of hemophilia A or B who underwent both MSK-US and MRI of the same joint within a maximum interval of three months were eligible for inclusion. Potentially eligible participants were identified through electronic medical records and the radiology database at the study center. Both pediatric and adult patients were eligible for inclusion.
Patients were excluded if they had a history of prior joint surgery, coexisting inflammatory or degenerative joint diseases (such as rheumatoid arthritis or advanced osteoarthritis), or incomplete or poor-quality imaging studies. Patients who underwent major clinical interventions between MSK-US and MRI examinations, including joint surgery or initiation of new disease-modifying treatment, were also excluded to minimize potential bias arising from interval changes in joint status.
A consecutive series of eligible patients meeting these criteria during the study period was included to minimize selection bias and improve representativeness. Baseline demographic and clinical characteristics, including age, hemophilia subtype (A or B), disease severity, prophylaxis status, inhibitor status, and joint type examined, were recorded to characterize the study population and assess its representativeness.
Unit of analysis and joint selection
To avoid statistical dependence arising from multiple joints within the same individual, a single joint per patient was included in the primary analysis. The representative joint was defined as the joint that underwent both MSK-US and MRI within the specified interval and was clinically indicated for imaging (e.g., symptomatic or designated target joint as documented in clinical records).
In cases where more than one joint met these criteria, the joint with the most complete imaging dataset and highest clinical relevance was selected according to predefined rules. Although this approach minimizes clustering bias and prevents overestimation of diagnostic performance, it may underestimate the burden of multi-joint involvement in hemophilia. To address this limitation, joints from different anatomical sites were included, and subgroup analyses were performed according to joint type (knee, ankle, and elbow).
Index test (musculoskeletal ultrasound)
MSK-US examinations were performed using high-frequency linear transducers (7–15 MHz) following a standardized scanning protocol for hemophilic joints, including the knee, ankle, and elbow14. These examinations were performed as formal diagnostic MSK-US studies by experienced musculoskeletal radiologists rather than point-of-care ultrasound (POCUS) assessments. Each joint was systematically evaluated in longitudinal and transverse planes, covering anterior, medial, lateral, and posterior recesses where applicable.
The following ultrasound features were assessed:
Synovial hypertrophy: non-compressible hypoechoic intra-articular tissue, with or without Doppler signal.
Joint effusion: compressible anechoic or hypoechoic intra-articular fluid without a Doppler signal. Joint effusion was recorded as a binary variable (present/absent). Any intra-articular fluid collection exceeding the physiological amount expected for the examined joint was considered positive. Because ultrasound cannot reliably distinguish simple effusion from hemarthrosis, fluid collections were interpreted in conjunction with clinical and MRI findings.
Cartilage damage: thinning, irregularity, or loss of the anechoic cartilage layer.
Bone erosion: cortical discontinuity visible in two perpendicular planes.
Ultrasound findings were recorded as both binary variables (present/absent) and ordinal variables using the HEAD-US scoring system11. According to the internationally accepted HEAD-US methodology, synovial hypertrophy was graded from 0 to 2, cartilage damage from 0 to 4, and bone changes from 0 to 2, yielding a composite HEAD-US score ranging from 0 to 8. Higher scores reflected greater severity of hemophilic arthropathy. Joint effusion was recorded separately and was not included in the composite HEAD-US score. ROC analyses and threshold-based assessments were performed using the composite HEAD-US score.
Ultrasound examinations were performed by experienced radiologists with at least 5 years of expertise in musculoskeletal imaging. Image interpretation for diagnostic accuracy analysis was conducted by a primary reader who was blinded to MRI findings. To assess reproducibility, a random subset of 40 patients (26.7%) was independently re-evaluated by a second radiologist who was blinded to both the initial ultrasound interpretation and MRI results. Interobserver agreement analyses were performed using the independent pre-consensus interpretations of both readers. Discrepancies were subsequently resolved by consensus discussion to establish an adjudicated imaging classification for descriptive purposes. For the remaining cases, the primary reader’s interpretation was used in the final diagnostic accuracy analysis.
Detailed MSK-US acquisition protocol
All examinations were performed using high-frequency linear-array transducers (7–15 MHz) with patients positioned comfortably to permit full access to the target joint. Grayscale imaging was performed before Doppler assessment. Images were acquired in longitudinal and transverse planes, and abnormalities were confirmed in at least two orthogonal planes before being recorded.
Knee joint:
The patient was positioned supine with the knee maintained in slight flexion (20–30°) using a cushion for support. The suprapatellar recess was examined in both longitudinal and transverse planes to assess for synovial hypertrophy and joint effusion. The medial and lateral parapatellar recesses were then scanned in longitudinal and transverse orientations. The femoral trochlear cartilage was evaluated with the knee in maximal flexion, and the cortical bone surfaces were assessed for erosive changes. Representative grayscale images of the suprapatellar recess, cartilage surface, and any abnormal findings were obtained and saved for documentation.
Ankle joint:
The patient was positioned supine with the ankle maintained in a neutral position. The anterior tibiotalar recess was examined in both longitudinal and transverse planes. The medial and lateral joint recesses were assessed for the presence of synovial hypertrophy and joint effusion. Accessible articular cartilage and cortical bone surfaces were evaluated for structural abnormalities. Representative images of the anterior recess and any abnormal findings were obtained and saved for documentation.
Elbow joint:
The patient was positioned either seated or supine with the elbow flexed to approximately 90°. The anterior humeroradial and humeroulnar recesses were assessed in both longitudinal and transverse planes. The posterior olecranon recess was then examined for evidence of synovial hypertrophy and joint effusion. Accessible cartilage surfaces and cortical bone margins were evaluated for structural abnormalities. Representative images of the anterior and posterior recesses, as well as any abnormal findings, were obtained and saved for documentation.
Doppler acquisition protocol:
Power Doppler imaging was performed when synovial hypertrophy was identified on grayscale imaging. Doppler settings were standardized across examinations, including low wall-filter settings, pulse repetition frequency of approximately 500–750 Hz, and gain adjusted to the highest level below background noise. Doppler signals were considered positive only when reproducible intra-synovial flow was identified on consecutive image acquisitions.
Image quality requirements and scoring workflow
For each joint, image acquisition included both longitudinal and transverse views of all accessible recesses. Bone erosions were required to be visible in two perpendicular planes. Images affected by substantial motion, poor acoustic contact, or inadequate visualization of target structures were excluded.
The scoring workflow consisted of: (1) assessment of synovial hypertrophy; (2) assessment of joint effusion; (3) evaluation of cartilage integrity; (4) evaluation of cortical bone surfaces; and (5) assignment of the composite HEAD-US score according to established criteria.
Reader training and standardization
Prior to image evaluation, participating radiologists underwent structured standardization sessions that included review of representative cases, consensus discussions on imaging definitions, and calibration using reference images based on the internationally accepted HEAD-US scoring system. These sessions were conducted to ensure uniform application of imaging criteria and to reduce interobserver variability before formal image interpretation.
Doppler assessment standardization
When Doppler imaging was used to assess synovial vascularity, machine settings (including gain, pulse repetition frequency, and wall filter) were standardized across examinations according to institutional protocol. Doppler findings were interpreted in conjunction with grayscale imaging, and only consistent intra-synovial signals were considered positive. The Doppler signal was not used as a standalone criterion for defining disease.
Interobserver agreement (ultrasound)
Interobserver agreement was assessed in a randomly selected subset of 40 patients (26.7%), which was independently evaluated by a second radiologist blinded to both the primary ultrasound interpretation and MRI findings. Agreement was quantified using Cohen’s kappa (κ) for binary variables and weighted kappa for ordinal HEAD-US scores. Kappa values were interpreted as follows: <0.20 poor agreement, 0.21–0.40 fair agreement, 0.41–0.60 moderate agreement, 0.61–0.80 good agreement, and >0.80 excellent agreement according to the criteria proposed by Landis and Koch15.
Agreement statistics were calculated using the independent initial interpretations of the two readers before consensus review. Discrepant cases were subsequently reviewed jointly and resolved by consensus for adjudication purposes, but were not used in the calculation of kappa statistics.
Reference standard (MRI)
MRI examinations were performed using standard joint imaging protocols and served as the reference standard because of their superior ability to detect both early synovial changes and deep structural abnormalities. MRI examinations were performed using 1.5-T or 3.0-T scanners equipped with dedicated joint coils appropriate for the anatomical site examined. Standardized imaging protocols included T1-weighted, T2-weighted, proton-density (PD) fat-suppressed, and gradient-echo sequences for the detection of hemosiderin deposition. Images were acquired in sagittal, coronal, and axial planes with slice thicknesses of approximately 3–4 mm.
The following MRI findings were recorded:
Synovial proliferation: thickened synovial tissue demonstrating characteristic signal intensity and expansion of the synovial recesses.
Hemarthrosis: intra-articular fluid demonstrating signal characteristics consistent with blood products and/or hemosiderin deposition on gradient-echo sequences.
Osteochondral damage: cartilage loss, surface irregularity, subchondral cyst formation, bone erosion, or other osteochondral abnormalities consistent with hemophilic arthropathy.
All MRI examinations underwent quality review before interpretation to ensure adequate visualization of synovial, cartilaginous, and osseous structures. Synovial proliferation was assessed based on synovial thickening and expansion of joint recesses. Hemarthrosis was identified using characteristic signal intensities of blood products and hemosiderin deposition, particularly on gradient-echo sequences. Osteochondral damage was evaluated by assessing cartilage thinning or loss, surface irregularity, subchondral cyst formation, bone erosions, and other structural abnormalities consistent with hemophilic arthropathy.
Where applicable, findings were assessed using established scoring systems such as the International Prophylaxis Study Group (IPSG) MRI scale16. When available, lesion severity was additionally characterized using the IPSG scoring system. MRI images were interpreted by a primary experienced radiologist who was blinded to MSK-US findings.
Interobserver agreement (MRI)
To assess reproducibility, a random subset of 40 patients (26.7%) was independently evaluated by a second radiologist who was blinded to both the primary MRI interpretation and MSK-US findings. Interobserver agreement was quantified using Cohen’s kappa (κ). Kappa values were interpreted according to the criteria proposed by Landis and Koch (<0.20 poor, 0.21–0.40 fair, 0.41–0.60 moderate, 0.61–0.80 good, and >0.80 excellent agreement).
Agreement statistics were calculated using the independent initial interpretations of both readers before consensus review. Cases with discrepant interpretations were subsequently adjudicated by consensus; however, consensus readings were not used for the calculation of interobserver agreement statistics.
Definition of test positivity
For diagnostic accuracy analyses, three complementary ultrasound-based definitions were evaluated.
Primary binary MSK-US definition:
Ultrasound results were dichotomized as positive or negative. The primary definition of a positive ultrasound result was the presence of at least one characteristic feature of hemophilic arthropathy, including synovial hypertrophy, joint effusion, cartilage damage, or bone erosion. This definition was used for the primary diagnostic endpoint and for the calculation of sensitivity, specificity, positive predictive value, negative predictive value, likelihood ratios, and overall diagnostic accuracy.
HEAD-US–based definition
The composite HEAD-US score was analyzed separately for ROC and threshold analyses. Consistent with the standard HEAD-US methodology, the score incorporated synovial hypertrophy, cartilage damage, and bone changes but excluded joint effusion. HEAD-US scores ranged from 0 to 8, with higher scores indicating greater arthropathy severity. To evaluate the robustness of the findings and the trade-off between sensitivity and specificity, pre-specified secondary analyses were performed using more stringent criteria, including:
Presence of ≥2 abnormal ultrasound features
Higher HEAD-US score thresholds (e.g., ≥2 and ≥3)
Feature-level definitions
Individual ultrasound features were also analyzed separately to assess lesion-specific diagnostic performance relative to MRI. For lesion-level analyses, synovial hypertrophy detected by MSK-US was compared with synovial proliferation on MRI; joint effusion on MSK-US was compared with MRI-detected hemarthrosis or intra-articular fluid abnormalities; cartilage damage on MSK-US was compared with MRI-detected osteochondral cartilage abnormalities; and bone erosion on MSK-US was compared with MRI-detected osseous abnormalities. These comparisons were based on the closest corresponding pathological features identifiable by both imaging modalities. MRI was considered positive if one or more characteristic abnormalities (synovial proliferation, hemarthrosis, or osteochondral damage) were present and served as the reference standard for all diagnostic accuracy analyses.
Blinding
Radiologists interpreting MSK-US images were blinded to MRI findings and detailed clinical information, except for the known diagnosis of hemophilia. MRI readers were similarly blinded to MSK-US results. All image interpretations were performed independently to minimize observer bias.
Time interval between tests
The interval between MSK-US and MRI examinations was restricted to a maximum of 3months to reflect real-world clinical practice. To reduce the potential for bias due to disease progression during this interval: (i) The time interval between examinations was recorded and analyzed. (ii) Patients undergoing major clinical interventions (e.g., surgery or initiation of new disease-modifying therapy) during this period were excluded. (iii) A pre-specified sensitivity analysis was performed for patients with an interval of ≤30 days.
The mean interval between MSK-US and MRI examinations was subsequently reported in the Results section to allow assessment of the potential influence of temporal changes in joint status on diagnostic performance. The distribution of imaging intervals was evaluated descriptively to assess the potential influence of disease progression between examinations. In addition to the mean interval, the median interval and range were recorded. Clinical records were reviewed for documented bleeding episodes, initiation or modification of prophylactic regimens, and other major interventions occurring between MSK-US and MRI examinations. Patients undergoing major clinical interventions during this period were excluded from the study. Patients with documented interval events that could substantially alter joint status were reviewed individually, and those meeting exclusion criteria were not included in the final analysis.
Statistical analysis
Statistical analyses were performed using SPSS Statistics for Windows (version 26.0) and R software (version 4.3.1). The primary analysis evaluated the diagnostic accuracy of MSK-US for detecting hemophilic arthropathy, using MRI as the reference standard. Sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and positive and negative likelihood ratios (LR+ and LR–) were calculated from 2 × 2 contingency tables, with corresponding 95% confidence intervals (CI). In addition to overall diagnostic accuracy, feature-level sensitivity, specificity, and agreement analyses were performed for synovial hypertrophy, joint effusion, cartilage damage, and bone erosion.
Receiver operating characteristic (ROC) curves were generated using ordinal HEAD-US scores, and the area under the curve (AUC) was calculated to assess overall diagnostic performance. Optimal cut-off values were determined using the Youden index.
Agreement between MSK-US and MRI findings, as well as interobserver agreement, was assessed using Cohen’s kappa (κ) for binary variables and weighted kappa for ordinal variables. Kappa values were interpreted according to the criteria proposed by Landis and Koch: <0.20 poor agreement, 0.21–0.40 fair agreement, 0.41–0.60 moderate agreement, 0.61–0.80 good agreement, and >0.80 excellent agreement.
Handling of indeterminate and missing data
All included MSK-US and MRI examinations were of sufficient diagnostic quality for interpretation. No indeterminate imaging findings remained after initial review, and no cases required exclusion because of unresolved diagnostic uncertainty following adjudication. Indeterminate findings identified during image review were independently evaluated by a second radiologist. Any discrepancies were resolved by consensus discussion. No cases remained unresolved after adjudication. Patients with missing or incomplete data for either MSK-US or MRI were excluded before analysis, and no imputation of missing data was performed.
Subgroup and exploratory analyses
Pre-specified subgroup analyses were conducted according to joint type (knee, ankle, and elbow) and type of pathology (synovial abnormalities versus structural abnormalities). For joint-specific analyses, sensitivity, specificity, and overall diagnostic performance were calculated separately for knee, ankle, and elbow joints.
Exploratory analyses evaluated the relationship between disease severity, as reflected by HEAD-US scores, and diagnostic performance. Threshold analyses were performed using increasing HEAD-US cut-off values to assess the trade-off between sensitivity and specificity across different levels of disease severity.
Sample size
No formal prospective sample-size calculation was performed because of the retrospective study design. However, sample precision was evaluated post hoc. Among the 150 included patients, 112 were MRI-positive, and 38 were MRI-negative. Assuming an expected sensitivity of approximately 90% and specificity of approximately 80%, this sample size provided 95% confidence intervals with widths of approximately ±5%–8% for sensitivity estimates and ±10%–15% for specificity estimates. The wider confidence intervals for specificity reflect the smaller number of MRI-negative cases available for analysis. Consequently, sensitivity estimates are considered more precise than specificity estimates and should be interpreted accordingly.