Feasibility and reliability of high-resolution phase-contrast quantitative-flow magnetic resonance imaging for assessing femoral head blood flow: a prospective study
Introduction
Ischemic necrosis of the femoral head (i.e., a loss of blood flow to the bone tissue) is caused by a myriad of factors, including corticosteroid use, excessive alcohol intake, and femoral neck fractures (1,2). These factors can precipitate circulatory disruptions, notably to the arterial blood supply and venous congestion (1). The early identification of aberrant blood supply patterns in the femoral head is crucial for initiating timely interventions to mitigate and prevent the progression of ischemic necrosis (3).
Various diagnostic techniques can be used to evaluate blood supply to the femoral head. For example, angiographic methods, such as digital subtraction angiography, high-resolution computed tomography angiography, and contrast-enhanced magnetic resonance (MR) angiography, can be used to visualize the blood vessels supplying the femoral head to a certain degree (4-8). However, these approaches predominantly delineate the morphologic abnormalities of the nutrient arteries. Difficulties may arise in analyzing early hemodynamic changes caused by corticosteroid use and excessive alcohol intake because these are not embolic infarctions (9). Bone scans and dynamic contrast-enhanced magnetic resonance imaging (MRI) can quantitatively detect perfusion anomalies but require the administration of radiotracers and gadolinium-based contrast agents (10). Ultrasonography may reveal deviations in blood circulation in the early phases of femoral head osteonecrosis, preceding the manifestations observed in plain radiography (11-13). Nevertheless, the accuracy and reproducibility of Doppler-based flow assessments are contingent on the operator, and the volumetric flow measurements are influenced by various factors, such as the angle dependence, variability in the estimation of the vessel area, and a lack of consideration of the spatial velocity profile (14-17).
The MR quantitative-flow (Q-flow) sequence, which uses the phase-contrast technique, has the capability to assess a range of hemodynamic parameters. This technique has been widely used in cardiovascular and cerebrovascular assessments (18). Q-flow allows for flow quantification without the need for contrast agent injection, offering a more precise estimation of flow volume compared to ultrasound by delivering a more accurate measurement of the vessel lumen and considering the spatial velocity across the actual vessel lumen (19). Thus, Q-flow is a potentially superior technique for precisely determining blood flow to the femoral head. However, the principal vessels supplying the femoral head [i.e., the medial femoral circumflex artery (MFCA) and the lateral femoral circumflex artery (LFCA), along with their branches] are characterized by small diameters, and the Q-flow faces challenges in accurately measuring small vessels due to errors induced by the partial volume effect (20,21). Consequently, Q-flow faces significant challenges in accurately quantifying the primary blood vessels feeding the femoral head, particularly the deep MFCA, which has a luminal diameter of approximately 1.6–1.8 mm (5).
Currently, research on the use of Q-flow to quantify blood supply to the femoral head is limited. Previous studies have shown the feasibility and accuracy of high-resolution Q-flow for quantifying blood flow in vessels with luminal diameters smaller than 2 mm, with reported errors of less than 6.2% in the phantom experiment (22,23). Therefore, we hypothesized that high-resolution phase-contrast Q-flow could be used to quantify blood flow to the femoral head. The main aims of this study were to use high-resolution phase-contrast Q-flow to quantify the feeding arteries and draining veins of the femoral head, and to ascertain its applicability and effectiveness in clinical settings. We present this article in accordance with the STARD reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-411/rc).
Methods
Participants
Between October 2023 and August 2024, this prospective study recruited a convenience sample from Deyang People’s Hospital, Deyang, comprising 10 healthy volunteers (4 men and 6 women), 10 women with systemic lupus erythematosus (SLE) on corticosteroid therapy, and 10 age-matched healthy women for the SLE patients, who served as controls. The characteristics of the study population are presented in Table 1.
Table 1
| Characteristic | Healthy volunteers (N=10) | SLE patients (N=10) | Healthy females (control group) (N=10) | P† |
|---|---|---|---|---|
| Female | 6 | 10 | 10 | – |
| Age (years) | 35±8.38 | 36.1±5.47 | 36.2±4.83 | 0.96 |
| Heart rate (beat/min) | 74±7.16 | 76.3±8.27 | 72.2±8.24 | 0.24 |
| Weight (kg) | 64 (50, 67.5) | 52.10±6.89 | 53.53±5.94 | 0.54 |
| Height (m) | 1.58 (1.5, 1.7) | 1.55±0.04 | 1.58±0.06 | 0.57 |
| BMI (kg/m2) | 22.22 (20.9, 25.09) | 21.54±2.40 | 21.52±1.64 | 0.49 |
Data are presented as mean ± standard deviation for normally distributed variables, and median (25th percentile, 75th percentile) for non-normally distributed variables. †, comparison between SLE patients and healthy females (control group). BMI, body mass index; SLE, systemic lupus erythematosus.
The healthy volunteers and controls were selected based on specific criteria that excluded those with a history of conditions affecting femoral head vascularization, such as hip trauma, diabetes mellitus, alcoholism, and steroid hormone use (7). Conventional MRI was performed on the SLE patients to rule out hips with imaging-based ischemic necrosis of the femoral head. Participants with contraindications to MRI, such as individuals with a pacemaker or claustrophobia, were excluded from the study.
The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Deyang People’s Hospital Ethics Committee (approval No. 2020-04-118-K01), and informed consent was obtained from all the participants before they underwent the Q-flow examination. This study has been registered in the China Medical Research Registration and Filing System (https://www.medicalresearch.org.cn/, No. MR-51-25-001006).
MRI protocol
All the participants were scanned using a 3T scanner (Ingenia; Philips Healthcare) equipped with a 32-channel torso coil covering both hip joints. The relaxation-enhanced angiography without contrast and triggering (REACT) technique was used to visualize the critical femoral vessels for the Q-flow sequence localization. The isotropic REACT technique employed in this study was a flow-independent, Dixon-based fat suppression method that effectively mitigates artifacts arising from both flow variability and inhomogeneous magnetic fields during fat suppression (24). This magnetization-prepared non-balanced dual-echo approach enables the simultaneous imaging of arteries and veins (24). It also enables the visualization of vital femoral nutrient arteries, including the MFCA and LFCA, and their deep branches, as well as corresponding draining veins such as the medial femoral circumflex vein (MFCV) and lateral femoral circumflex vein (LFCV) (25,26).
For the imaging procedure, after a half-hour of rest, the patients were positioned supine with both lower limbs abducted to approximately 30°, as abduction may affect femoral blood flow (12). After the REACT acquisition, the two-dimensional Q-flow slice with an in-plane resolution of 0.4×0.4 (interpolated to 0.2×0.2) mm2, oriented perpendicular to the target vessel, was determined based on the maximum intensity projection (MIP) and multi-planar reconstruction (MPR) of the REACT image for flow parameter measurement as depicted in Figure 1. The velocity encoding was set to 50 cm/s for the MFCA and 40 cm/s for the deep MFCA and ascending LFCA, guided by flow velocity measurements obtained by ultrasound (11). The imaging parameters are detailed in Table 2.
Table 2
| Parameters | REACT | Q-flow |
|---|---|---|
| Scan mode | 3D | 2D |
| Technique | FFE | FFE |
| Slice orientation | Coronal | Oblique |
| Scan duration (s) | 03:25 | 04:32 |
| TR/TE1/TE2 or TR/TE (ms) | 4.0/1.31/2.3 | 9.6/4.8 |
| TI delay (ms) | 55.5 | – |
| Phase-encoding direction | RL | RL/AP/RL |
| FOV (mm3) (read × phase × slice) | 350×448×135 | 200×139×5 |
| Acquisition voxel size (mm3) (read × phase × slice) | 1.5×1.5×1.5 | 0.4×0.4×5 |
| Reconstruction voxel size (mm3) | 0.75×0.75×0.75 | 0.2×0.2×5 |
| CS-SENSE | 4 | No |
| Echoes | 2 | 1 |
| Partial echo | No | Yes |
| Slice oversampling | 1 | – |
| Water-fat shift | 0.313 | 2 |
| Flip angle | 10° | 18° |
| No. of signals acquired | 1 | 1 |
| T2prep (ms) | 50 | No |
| Cardiac gating (synchronization/device/heart phases) | – | Retrospective/PPU/18 |
2D, two-dimensional; 3D, three-dimensional; AP, anterior-posterior; CS, compressed sensing technique; FFE, fast field echo; FOV, field of view; PPU, peripheral pulse unit; Q-flow, quantitative flow; REACT, relaxation-enhanced angiography without contrast and triggering; RL, right-left; SENSE, sensitivity encoding; TE, echo time; TI, inversion time; TR, repetition time.
The Q-flow measurement positions for the target vessels were meticulously determined by an imaging technologist with an associate title under the guidance of a deputy chief physician. In measuring the MFCA, the origin of the MFCA was identified anterior to the commencement of the inferior retinacular artery, which signified the blood supply of the femoral head through the superior retinacular artery and inferior retinacular artery (27,28). The position of the deep MFCA was pinpointed following the branching out of the inferior retinacular artery and the descending branch of the MFCA, which indicated the femoral blood supply through the superior retinacular artery (28). In addition, the location of the ascending LFCA was determined at the origin of the ascending LFCA, corresponding to the femoral blood supply through the anterior retinacular artery originating from the ascending LFCA (25,29). The associated veins (the one closest to the measured artery), specifically the lateral and medial circumflex femoral veins, which typically run parallel to their corresponding arteries, were simultaneously measured (26,30). The Q-flow measurement positions are shown in Figure 1.
For each participant, three Q-flow slices were acquired per hip to measure the MFCA, deep MFCA, and LFCA, totaling six Q-flow slices across both hips, with a total Q-flow scan time of approximately 30 minutes per participant. To assess the retest reliability, each healthy volunteer underwent two scans in a single session, separated by a 30-minute interval during which they got up from the MRI table and left the examination room. This allowed for independent positioning and the setup for the second scan. Further, a Q-flow scan was performed on the SLE patients and the controls to assess the efficacy of high-resolution phase-contrast Q-flow in quantifying femoral flow variations.
Quantitative measurements
Two readers, reader 1 with 11 years of experience and reader 2 with 6 years of experience in interpreting osteoarticular images, were blinded to the patient and volunteer specifics. The readers delineated the contours of the target arteries and the accompanying veins in the Q-flow images using the Q-flow analysis software package (Philips Healthcare) to quantify the blood flow parameters. The accuracy of vessel delineation was meticulously examined to ensure precise alignment with vessel boundaries on the magnitude, phase, and velocity images. When necessary, proper alignment with the velocity images was achieved through adjustments. Subsequently, the flow parameters, including stroke (the absolute difference between forward and backward flows), regurgitant fraction (RF; representing the fraction of backward to forward flow), mean flux (MF), stroke distance (SD; representing the net distance blood travels in the vessel per beat), mean velocity (MV), and peak systolic velocity (PSV), were automatically computed. The reliability of Q-flow in quantifying blood flow in the femoral head was evaluated based on established methods for reliability assessment as described previously (31,32), including retest, intra-, and inter-rater reliability as presented in Figure 2. Reader 1 measured the scan data of the SLE and the control groups to investigate early blood flow changes.
Statistical analysis
The retest reliability, and inter- and intra-rater reliability were assessed using the paired-samples t-test (for normal distributions) or the Wilcoxon rank-sum test (for non-normal distributions), intra-class correlation coefficients (ICCs), and a Bland-Altman analysis. An ICC value between 0.61 and 0.80 indicated good reliability, while a value exceeding 0.80 indicated excellent reliability; good or excellent reliability (ICC >0.6) was considered acceptable (33). The quantitative results between the SLE and control groups were analyzed using the Wilcoxon rank-sum test (for non-normal distributions) or the t-test (for normal distributions), with the Shapiro-Wilk test used to ascertain distribution normality. All the statistical analyses were conducted using Rstudio (version 1.4.1717; https://posit.co/).
Results
Study cohort
A total of 20 hips from 10 healthy volunteers (4 men and 6 women) were included in the study to assess the reliability. None of the SLE patients exhibited any ischemic findings on conventional MRI. Thus, 20 hips from 10 SLE patients (mean SLE duration: 0.95±0.44 years, range: 6 months to 1.5 years; mean corticosteroid use: 0.89±0.47 years, range: 3 months to 1.5 years), and 20 hips from 10 aged-matched controls were analyzed to evaluate the efficacy of assessing early blood flow changes. In the hip evaluation, the measurement of the deep MFCVs proved challenging, and had a success rate of only 66.7%. Conversely, nearly all the other vessels were successfully measured (see Table 3).
Table 3
| Vessel | Healthy volunteers | SLE patients | Healthy females (control group) |
|---|---|---|---|
| MFCA | 19 [1] | 20 | 19 [1] |
| LFCA | 20 | 20 | 20 |
| Deep MFCA | 18 [2] | 18 [2] | 18 [2] |
| MFCV | 19 [1] | 20 | 18 [2] |
| LFCV | 20 | 20 | 20 |
| Deep MFCV | 12 [8] | 14 [6] | 14 [6] |
Data are presented as number of successful measurements [unsuccessful measurements]. LFCA, lateral femoral circumflex artery; LFCV, lateral femoral circumflex vein; MFCA, medial femoral circumflex artery; MFCV, medial femoral circumflex vein; SLE, systemic lupus erythematosus.
Retest, intra-rater, and inter-rater reliability
No significant differences in the flow parameters were observed between the two scans in relation to retest reliability, nor between the measurements for intra-rater and inter-rater reliability [P=0.08–1/not applicable (NA)], except for the RF (P=0.01, P<0.01 and P=0.04 for the LFCA, deep MFCA, and LFCV for the retest reliability). The detailed comparative results are presented in Table 4.
Table 4
| Parameters | Data 1 | Data 2 | Measurement 1 | Measurement 2 | P† | P‡ | P§ |
|---|---|---|---|---|---|---|---|
| MFCA | |||||||
| Stroke (mL) | 0.18±0.11 | 0.17±0.1 | 0.18±0.11 | 0.18±0.11 | 0.09 | 1 | 0.38 |
| RF | 0.23% (0, 10%) | 0 (0, 3%) | 0.7% (0, 2%) | 0 (0, 1%) | 0.61 | 0.27 | 0.18 |
| MF (mL/s) | 0.22±0.12 | 0.2±0.1 | 0.22±0.11 | 0.22±0.12 | 0.28 | 0.26 | 0.37 |
| SD (cm) | 2.4±1.03 | 2.28±1.03 | 2.35±1.01 | 2.31±0.93 | 0.38 | 0.31 | 0.24 |
| MV (cm/s) | 2.85±1.05 | 2.81±1.45 | 2.78±1 | 2.75±0.94 | 0.87 | 0.33 | 0.08 |
| PSV (cm/s) | 26.87±7.11 | 25.57±5.62 | 26.87±7.11 | 26.87±7.11 | 0.16 | NA | 0.90 |
| LFCA | |||||||
| Stroke (mL) | 0.21±0.12 | 0.22±0.13 | 0.21±0.13 | 0.21±0.12 | 0.44 | 0.94 | 0.38 |
| RF | 0 (0, 0) | 0.06% (0, 1%) | 0 (0, 0) | 0 (0, 0) | 0.01 | 1 | 0.18 |
| MF (mL/s) | 0.27±0.15 | 0.24±0.14 | 0.27±0.16 | 0.27±0.15 | 0.31 | 0.27 | 0.37 |
| SD (cm) | 2.72±0.97 | 2.1 (1.53, 3.49) | 2.31 (1.97, 3.53) | 2.66±1 | 0.64 | 0.61 | 0.24 |
| MV (cm/s) | 3.32±1.04 | 3.08±1.23 | 3.27±1.06 | 3.23±1.09 | 0.29 | 0.3 | 0.08 |
| PSV (cm/s) | 23.08±7.64 | 23.24±7.53 | 22.88±7.66 | 23.08±7.63 | 0.86 | 0.33 | 0.90 |
| Deep MFCA | |||||||
| Stroke (mL) | 0.07±0.05 | 0.07 (0.04, 0.08) | 0.08±0.05 | 0.07±0.05 | 0.09 | 0.17 | 1 |
| RF | 0 (0, 0) | 15.7% (6%, 39%) | 0 (0, 1%) | 0 (0, 1%) | <0.01 | 0.53 | 0.59 |
| MF (mL/s) | 0.09±0.06 | 0.08 (0.06, 0.09) | 0.09±0.05 | 0.1 (0.05, 0.12) | 0.29 | 0.24 | 0.71 |
| SD (cm) | 1.95±0.96 | 1.76±0.76 | 1.99±1.13 | 1.93±1.13 | 0.12 | 0.28 | 0.89 |
| MV (cm/s) | 2.33±1.02 | 2.15±0.83 | 2.33±1.14 | 2.29±1.18 | 0.12 | 0.59 | 0.71 |
| PSV (cm/s) | 20.82±5.68 | 19.67±5.16 | 20.79±5.64 | 20.82±5.42 | 0.26 | 0.83 | 1 |
| MFCV | |||||||
| Stroke (mL) | 0.05 (0.02, 0.11) | 0.05 (0.02, 0.11) | 0.05 (0.02, 0.11) | 0.04 (0.02, 0.11) | 0.89 | 1 | 0.34 |
| RF | 0 (0, 5%) | 0 (0, 0) | 0 (0, 5%) | 0.03% (0, 5%) | 0.83 | 0.36 | 0.8 |
| MF (mL/s) | 0.06 (0.04, 0.11) | 0.09±0.06 | 0.06 (0.04, 0.13) | 0.06 (0.04, 0.13) | 0.75 | 1 | 0.36 |
| SD (cm) | 0.9 (0.7, 1.55) | 1.28 (1.08, 1.6) | 0.94 (0.62, 1.63) | 0.85 (0.59, 1.68) | 0.19 | 0.28 | 0.2 |
| MV (cm/s) | 1.29 (0.94, 2.01) | 1.4 (1.23, 1.87) | 1.17 (0.83, 1.89) | 1.41±0.9 | 0.34 | 0.3 | 0.2 |
| PSV (cm/s) | 14.84±3.14 | 15.58±4.22 | 15.48±3.27 | 15.62±3.55 | 0.37 | 0.63 | 0.09 |
| LFCV | |||||||
| Stroke (mL) | 0.07 (0.03, 0.2) | 0.09 (0.04, 0.19) | 0.08 (0.03, 0.18) | 0.08 (0.03, 0.2) | 0.95 | 1 | 0.28 |
| RF | 0 (0, 0) | 0 (0, 0) | 0 (0, 0) | 0 (0, 0) | 0.04 | 0.42 | 0.59 |
| MF (mL/s) | 0.09 (0.03, 0.26) | 0.125 (0.09, 0.25) | 0.11 (0.05, 0.26) | 0.11 (0.04, 0.26) | 0.11 | 0.36 | 0.5 |
| SD (cm) | 1.41 (0.58, 2.7) | 1.83±1.14 | 1.14 (0.58, 2.23) | 1.01 (0.59, 2.23) | 0.40 | 0.58 | 0.24 |
| MV (cm/s) | 1.76 (0.75, 3.41) | 2.52±1.64 | 1.42 (0.75, 3.01) | 1.23 (0.73, 3.01) | 0.12 | 0.1 | 0.31 |
| PSV (cm/s) | 9.98 (6.91, 14.46) | 11.3±3.78 | 11.4±4.16 | 11.64±4.03 | 0.99 | 0.33 | 0.21 |
| Deep MFCV | |||||||
| Stroke (mL) | 0.06±0.05 | 0.04 (0.01, 0.11) | 0.07±0.05 | 0.06±0.05 | 0.93 | 0.39 | 0.28 |
| RF | 0 (0, 29%) | 2.65% (2%, 3%) | 0 (0, 29%) | 0 (0, 15%) | 0.69 | 0.79 | 1 |
| MF (mL/s) | 0.06±0.06 | 0.08±0.07 | 0.07±0.05 | 0.07±0.06 | 0.23 | 0.36 | 0.10 |
| SD (cm) | 1.03±0.82 | 1.3±0.89 | 1.08±0.96 | 1.09±0.95 | 0.13 | 0.38 | 0.38 |
| MV (cm/s) | 1.18±0.93 | 1.63±1.15 | 1.22±1.02 | 1.26±1.1 | 0.11 | 0.25 | 0.34 |
| PSV (cm/s) | 15.03±3.04 | 14.485 (12.18, 19.65) | 14.76±3.23 | 15.11±3 | 0.34 | 0.07 | 0.21 |
Data are presented as mean ± standard deviation for normally distributed variables, and median (25th percentile, 75th percentile) for non-normally distributed variables. †, comparison between two scan data (data 1 and data 2) for evaluating retest reliability; ‡, comparison between two measurements (measurement 1 and measurement 2) by the same radiologist for evaluating intra-rater reliability; §, comparison between two radiologists (data 1 and measurement 1) for evaluating inter-rater reliability. LFCA, lateral femoral circumflex artery; LFCV, lateral femoral circumflex vein; MF, mean flux; MFCA, medial femoral circumflex artery; MFCV, medial femoral circumflex vein; MV, mean velocity; PSV, peak systolic velocities; RF, regurgitant fraction; SD, stroke distance.
The agreement for feeding arteries was not always acceptable for the RF (ICC: 0.01–0.45 and 0.53–0.97 for the retest and intra-rater reliability), but it was acceptable for all the other flow parameters in terms of the retest reliability (ICC: 0.92–0.99, 0.61–0.9, 0.67–0.83, 0.63–0.86, and 0.7–0.86 for the stroke, MF, SD, MV, and PSV, respectively), intra-rater reliability (ICC: 0.97–1, 0.98–0.99, 0.95–0.98, 0.96–0.98, and 0.99–1 for the stroke, MF, SD, MV, and PSV, respectively), and inter-rater reliability (ICC: 1, 0.99–1, 0.93–0.98, 0.92–0.97, and 0.99–1 for the stroke, MF, SD, MV, and PSV, respectively). For the draining veins, the RF, MV, and PSV did not consistently achieve acceptable agreement in terms of retest reliability (ICC: 0.04–0.9, 0.46–0.84, and 0.5–0.73, respectively), while the other parameters consistently exhibited acceptable agreement in terms of the retest reliability (ICC: 0.9–0.97, 0.63–0.75, and 0.75–0.92 for stroke, MF, and SD, respectively), intra-rater reliability (ICC: 0.84–1, 0.87–1 and 0.93–1 for stroke, MF and SD respectively), and inter-rater reliability (ICC: 0.9–0.99, 0.87–0.99 and 0.9–0.99 for the stroke, MF and SD, respectively). The ICC values are illustrated in Figure 3.
In the Bland-Altman analysis of retest reliability, the percentages of points outside the 95% agreement threshold were 0–10.53% for arteries and 0–16% for veins, with high RF values at 8.33–16.67%. In relation to inter-rater reliability, the values were 5–11.11% for arteries and 5–16.67% for veins, with a high of 16.67% for the deep MFCV in relation to RF. In relation to intra-rater reliability, the artery values were 0–11.11%, and the vein values were 5–15%, with a high of 15% for the LFCV in terms of the MV. The Bland-Altman analysis results are presented in Figures S1-S3.
Efficacy of assessing early blood flow changes
In relation to the feeding arteries, the SLE group had lower blood flow parameters than the control group. Specifically, significant decreases in the MFCA (in the stroke and SD, P=0.01 and P<0.01), LFCA (in the stroke, MF, SD, MV, and PSV, all P<0.01), and deep MFCA (in the stroke, MF, SD, MV, and PSV, all P<0.01) were observed. In relation to the draining veins, the SLE group had lower flow parameters than the control group for the LFCV (in the SD and MV, P<0.01 and P=0.02) and deep MFCV (in the stroke, MF, SD, and PSV, P<0.01, P=0.02, P=0.01 and P=0.03). These outcomes are detailed in Table 5 and Figure 4. Sample Q-flow images of healthy controls and SLE patients are presented in Figure 5.
Table 5
| Parameters | Controls | SLE patients | P value |
|---|---|---|---|
| MFCA | |||
| Stroke (mL) | 0.21±0.11 | 0.14±0.05 | 0.01 |
| RF | 0 (0, 1%) | 0 (0, 2%) | 0.48 |
| MF (mL/s) | 0.24±0.11 | 0.18±0.07 | 0.06 |
| SD (cm) | 2.73±0.84 | 2±0.63 | <0.01 |
| MV (cm/s) | 3.13±0.71 | 2.73±1 | 0.16 |
| PSV (cm/s) | 28.33±6.42 | 24.83±5.23 | 0.07 |
| LFCA | |||
| Stroke (mL) | 0.2±0.09 | 0.1±0.04 | <0.01 |
| RF | 0 (0, 0) | 0 (0, 0) | 0.33 |
| MF (mL/s) | 0.27 (0.18, 0.3) | 0.12 (0.08, 0.16) | <0.01 |
| SD (cm) | 2.62±0.76 | 1.64±0.67 | <0.01 |
| MV (cm/s) | 3.08±0.93 | 2.21±0.97 | <0.01 |
| PSV (cm/s) | 20.87 (19.96, 23.81) | 18.11 (14.98, 19.37) | <0.01 |
| Deep MFCA | |||
| Stroke (mL) | 0.08±0.04 | 0.04±0.02 | <0.01 |
| RF | 0 (0, 0) | 0 (0, 5%) | 0.14 |
| MF (mL/s) | 0.09±0.04 | 0.05±0.03 | <0.01 |
| SD (cm) | 1.94 (1.64, 2.64) | 1.15 (0.86, 1.6) | <0.01 |
| MV (cm/s) | 2.56±0.9 | 1.68±0.83 | <0.01 |
| PSV (cm/s) | 19.85 (18.8, 22.81) | 16.23 (13.35, 18.94) | <0.01 |
| MFCV | |||
| Stroke (mL) | 0.07 (0.02, 0.15) | 0.07 (0.05, 0.09) | 0.97 |
| RF | 0 (0, 7%) | 0.5% (0, 6%) | 0.70 |
| MF (mL/s) | 0.08 (0.03, 0.16) | 0.09 (0.06, 0.12) | 0.85 |
| SD (cm) | 0.91 (0.58, 1.55) | 0.75 (0.49, 1.03) | 0.16 |
| MV (cm/s) | 1.13 (0.76, 1.55) | 1.02 (0.59, 1.39) | 0.46 |
| PSV (cm/s) | 15.93±3.06 | 17.69±3.01 | 0.08 |
| LFCV | |||
| Stroke (mL) | 0.17±0.13 | 0.11±0.08 | 0.06 |
| RF | 0 (0, 0) | 0 (0, 0) | 0.20 |
| MF (mL/s) | 0.21±0.15 | 0.14±0.09 | 0.10 |
| SD (cm) | 1.7 (1.39, 2.7) | 1.08 (0.7, 1.54) | <0.01 |
| MV (cm/s) | 1.99 (1.6, 3.41) | 1.34 (0.96, 2.04) | 0.02 |
| PSV (cm/s) | 12.27 (9.26, 14.08) | 10.57 (9.2, 13.09) | 0.53 |
| Deep MFCV | |||
| Stroke (mL) | 0.05 (0.02, 0.09) | 0.02 (0, 0.03) | <0.01 |
| RF | 0 (0, 6%) | 8.4% (1%, 30%) | 0.08 |
| MF (mL/s) | 0.06±0.04 | 0.03±0.02 | 0.02 |
| SD (cm) | 1.32±0.78 | 0.63±0.48 | 0.01 |
| MV (cm/s) | 1.47±0.91 | 0.9±0.66 | 0.07 |
| PSV (cm/s) | 15.59±1.59 | 13.65±2.68 | 0.03 |
Data are presented as mean ± standard deviation for normally distributed variables, and median (25th percentile, 75th percentile) for non-normally distributed variables. LFCA, lateral femoral circumflex artery; LFCV, lateral femoral circumflex vein; MF, mean flux; MFCA, medial femoral circumflex artery; MFCV, medial femoral circumflex vein; MV, mean velocity; PSV, peak systolic velocities; RF, regurgitant fraction; SD, stroke distance; SLE, systemic lupus erythematosus.
Discussion
The evaluation of femoral head blood supply abnormalities is crucial for clinical diagnosis and treatment. This study successfully used the MR phase-contrast Q-flow technique to reliably quantify blood flow to the femoral head. Further, it validated the capacity of Q-flow to assess blood flow changes.
Although bone scans and dynamic contrast-enhanced MRI can monitor early blood supply abnormalities, they require the administration of radiotracers and gadolinium contrast agents (10). Color Doppler can detect changes in blood velocity in early osteonecrosis, but it faces challenges related to accuracy and repeatability, such as insonation angle problems, and an inability to account for spatial velocity profiles (14-17,19). The MR phase-contrast Q-flow technique offers an alternative approach for blood supply quantification, providing the benefits of contrast-free and accurate measurements. To successfully use Q-flow to reliably quantify blood flow to the femoral head, the MRI protocol was optimized. First, limitations related to the partial volume effect were addressed by enhancing the resolution. Previous studies have reported that a resolution of approximately one-quarter of the vessel diameter (pixel/vessel: 1.14/5) allows for the accurate measurement of blood flow parameters (34), and 16 voxels per vessel cross-section can achieve measurement accuracy within a 10% error margin (20,21). Therefore, this study adopted a similar imaging resolution strategy (pixel/vessel: 0.4/1.6, approximately 64 interpolated voxels per cross-section) to ensure the accuracy of the measurement results.
The imaging resolution of this study (0.4×0.4 mm2) is comparable to the resolutions of previous studies (0.35–0.5 mm) that successfully measured blood flow in vessels smaller than 2 mm in luminal diameter (22,23). In addition, unlike previous studies that focused only on the MFCA and LFCA (11,13), this study extended the measurements to include the deep MFCA (which supplies the femoral head directly) and ascending LFCA (which supplies the femoral head through the anterior retinacular artery). This refined protocol improved the accuracy of femoral blood supply assessment by omitting branches not serving the femoral head, including the descending branches of the MFCA, and the descending and transverse branches of the LFCA (25,27). The measurement positions, located away from the femoral head, may facilitate analysis in cases of femoral neck fracture and core decompression.
High-resolution Q-flow demonstrated acceptable reliability in measuring the femoral head vessels, yielding acceptable retest agreement in all the arterial flow parameters except the RF, along with strong inter- and intra-agreement in this study. However, a few MFCAs and their corresponding deep branches could not be detected by REACT because of the flow variability to the femoral head. For example, the absence of the MFCA may occur if the piriformis branch predominantly supplies blood to the femoral head (5,35). Despite these challenges, the MPR and MIP images of the REACT facilitated the precise localization of the Q-flow scanning plane, ensuring reliable repeatability for the feeding arteries (only the RF did not consistently achieve acceptable agreement). Although these vessels are diminutive in size, high-resolution Q-flow effectively illustrated the lumen, enabling the precise delineation of the regions of interest (ROIs) along the lumen. Consequently, this led to strong inter- and intra-agreement in measuring the target vessels, except for the RF. When arterial and venous vessels are closely situated, delineating ROIs may inadvertently incorporate peripheral pixels from accompanying reversed-flow vessels, which may lead to considerable variability in RF between measurements. Nonetheless, its effect on other parameters is limited because peripheral pixels typically exhibit lower flow rates and contribute less than the central pixel (36). Additionally, successfully measuring the deep MFCVs using Q-flow is a significant challenge due to their extremely sluggish blood flow.
This study further validated the utility of Q-flow in detecting flow reduction in the femoral head of SLE patients on corticosteroids. Corticosteroids induce a femoral intra-osseous compartment syndrome characterized by the hyperplasia of marrow fat cells, intra-osseous hypertension, vascular compression, and intravascular coagulation (9), ultimately resulting in hemodynamic changes in both blood-supplying arteries and draining veins. In addition, corticosteroid-induced osteoporosis may further contribute to reduced perfusion in the femoral head (37,38).
Previous studies have shown that the main blood supply from the superior retinacular arteries is significantly impaired in the early stage of corticosteroid use (39,40). These studies indicate that the vascular injury process begins with impaired blood flow in the true superior retinacular arteries. If perfusion is not restored, the femoral head enters a silent necrotic phase, during which radiographs may still appear normal. In some cases, revascularization may occur through small peripheral neovascularization extending into the ischemic region; however, this process is often disrupted by structural damage, ultimately leading to irreversible osteonecrosis. This sequence typically unfolds over approximately 13.2 months.
Using dynamic MRI, Nakamura et al. (41) showed that within four months of corticosteroid therapy, peak enhancement values in the femoral heads of adult patients recovered from 10.7% to 12.7%, but remained lower than the 20% observed in healthy adults. In addition, an early decreased blood supply was also found in steroid-induced osteonecrosis of the rabbit femoral head (42). Our study found comparable blood flow alterations in the SLE patients who had received hormone therapy for an average of less than 9 months, evidenced by reduced blood flow parameters in the LFCAs, MFCAs, and their deep branches, highlighting the ability of high-resolution Q-flow imaging to detect early changes in arterial blood flow to the femoral head. The early administration of hormones can significantly compromise the blood supply to the femoral head; however, this does not inevitably lead to femoral head necrosis, as the initial blood supply impairment is often followed by a subsequent process of revascularization (39-41). Further, this process may lead to a blood supply that exceeds the physiological blood supply levels of the femoral head (13,43). Consequently, Q-flow is a promising tool for dynamically assessing the recovery of femoral head blood flow.
This study had two limitations. First, it had a small sample size. Second, although we aimed to align the measurement plane perpendicular to the artery, it was not necessarily perpendicular to the accompanying vein, introducing potential inaccuracies in venous measurements. This challenge was at times exacerbated by the presence of multiple accompanying veins, complicating precise venous drainage quantification. Additionally, blood flow is influenced by factors such as age and cardiac function (40). Thus, relative blood flow parameters normalized by these factors may be more practical (e.g., target vessel flow relative to femoral artery flow or cardiac output). However, for longitudinal assessments of specific patients (e.g., for evaluating the effect of core decompression) or when using the healthy side as a reference (e.g., for evaluating unilateral femoral neck fractures), the use of measured raw parameters is also feasible (44).
Conclusions
The high-resolution MR phase-contrast Q-flow technique, which had acceptable reliability, was able to detect early changes in arterial blood flow to the femoral head, highlighting its potential clinical applicability.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the STARD reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-411/rc
Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-411/dss
Funding: This work was supported by
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-411/coif). X.Z. is an employee of Philips Healthcare. All authors report funding from the Sichuan Medical Association (Sichuan Province Medical Research Youth Innovation Project) (No. Q20020). The authors have no other conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The study was approved by the Deyang People’s Hospital Ethics Committee (approval No. 2020-04-118-K01), and informed consent was obtained from all the participants.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Mont MA, Salem HS, Piuzzi NS, Goodman SB, Jones LC. Nontraumatic Osteonecrosis of the Femoral Head: Where Do We Stand Today?: A 5-Year Update. J Bone Joint Surg Am 2020;102:1084-99. [Crossref] [PubMed]
- Zhao D, Zhang F, Wang B, Liu B, Li L, Kim SY, Goodman SB, Hernigou P, Cui Q, Lineaweaver WC, Xu J, Drescher WR, Qin L. Guidelines for clinical diagnosis and treatment of osteonecrosis of the femoral head in adults (2019 version). J Orthop Translat 2020;21:100-10. [Crossref] [PubMed]
- Liu LH, Zhang QY, Sun W, Li ZR, Gao FQ. Corticosteroid-induced Osteonecrosis of the Femoral Head: Detection, Diagnosis, and Treatment in Earlier Stages. Chin Med J (Engl) 2017;130:2601-7. [Crossref] [PubMed]
- Wang B, Li L, Wang Y, Wang Z, Li C, Fu W, Qiu X, Zhao D. Digital Subtraction Angiography and Magnetic Resonance Imaging-Based Staging of Circulatory Obstruction in the Femoral Head During Osteonecrosis of the Femoral Head Development. Ann Plast Surg 2020;85:677-84. [Crossref] [PubMed]
- Zlotorowicz M, Czubak J, Kozinski P, Boguslawska-Walecka R. Imaging the vascularisation of the femoral head by CT angiography. J Bone Joint Surg Br 2012;94:1176-9. [Crossref] [PubMed]
- Zlotorowicz M, Czubak J, Caban A, Kozinski P, Boguslawska-Walecka R. The blood supply to the femoral head after posterior fracture/dislocation of the hip, assessed by CT angiography. Bone Joint J 2013;95-B:1453-7. [Crossref] [PubMed]
- Liao Z, Bai Q, Ming B, Ma C, Wang Z, Gong T. Detection of vascularity of femoral head using sub-millimeter resolution steady-state magnetic resonance angiography-initial experience. Int Orthop 2020;44:1115-21. [Crossref] [PubMed]
- Liao Z, Liu C, Wu B, Ma C, Ming B, Zhou Q, Zhang X, Zhou S, Chen Y. Evaluating blood supply changes in the osteonecrosis of the femoral head using gadobutrol-based steady-state MR angiography. Eur Radiol 2023;33:8597-604. [Crossref] [PubMed]
- Hines JT, Jo WL, Cui Q, Mont MA, Koo KH, Cheng EY, et al. Osteonecrosis of the Femoral Head: an Updated Review of ARCO on Pathogenesis, Staging and Treatment. J Korean Med Sci 2021;36:e177. [Crossref] [PubMed]
- Dyke JP, Aaron RK. Noninvasive methods of measuring bone blood perfusion. Ann N Y Acad Sci 2010;1192:95-102. [Crossref] [PubMed]
- Graif M, Schweitzer ME, Nazarian L, Matteucci T, Goldberg BB. Color Doppler hemodynamic evaluation of flow to normal hip. J Ultrasound Med 1998;17:275-80; quiz 281-2.
- Yousefzadeh DK, Jaramillo D, Johnson N, Doerger K, Sullivan C. Biphasic threat to femoral head perfusion in abduction: arterial hypoperfusion and venous congestion. Pediatr Radiol 2010;40:1517-25. [Crossref] [PubMed]
- Hassan SZ, Abdalla MA, El-Azizi HM, Mohamed NN. Ultrasonographic evaluation of haemodynamic flow changes to the hip in systemic lupus erythematosus patients: Correlation with risk factors for osteonecrosis. The Egyptian Rheumatologist 2018;40:161-6.
- Krishnamurthy U, Yadav BK, Jella PK, Haacke EM, Hernandez-Andrade E, Mody S, Yeo L, Hassan SS, Romero R, Neelavalli J. Quantitative Flow Imaging in Human Umbilical Vessels In Utero Using Nongated 2D Phase Contrast MRI. J Magn Reson Imaging 2018;48:283-9. [Crossref] [PubMed]
- Du Y, Ding H, He L, Yiu BYS, Deng L, Yu ACH, Zhu L. Quantitative Blood Flow Measurements in the Common Carotid Artery: A Comparative Study of High-Frame-Rate Ultrasound Vector Flow Imaging, Pulsed Wave Doppler, and Phase Contrast Magnetic Resonance Imaging. Diagnostics (Basel) 2022;12:690. [Crossref] [PubMed]
- Applegate GR, Thaete FL, Meyers SP, Davis PL, Talagala SL, Recht M, Wozney P, Kanal E. Blood flow in the portal vein: velocity quantitation with phase-contrast MR angiography. Radiology 1993;187:253-6. [Crossref] [PubMed]
- Blanco P. Volumetric blood flow measurement using Doppler ultrasound: concerns about the technique. J Ultrasound 2015;18:201-4. [Crossref] [PubMed]
- Markl M, Schnell S, Wu C, Bollache E, Jarvis K, Barker AJ, Robinson JD, Rigsby CK. Advanced flow MRI: emerging techniques and applications. Clin Radiol 2016;71:779-95. [Crossref] [PubMed]
- Maier IL, Hofer S, Joseph AA, Merboldt KD, Tan Z, Schregel K, Knauth M, Bähr M, Psychogios MN, Liman J, Frahm J. Carotid artery flow as determined by real-time phase-contrast flow MRI and neurovascular ultrasound: A comparative study of healthy subjects. Eur J Radiol 2018;106:38-45. [Crossref] [PubMed]
- Tang C, Blatter DD, Parker DL. Accuracy of phase-contrast flow measurements in the presence of partial-volume effects. J Magn Reson Imaging 1993;3:377-85. [Crossref] [PubMed]
- Greil G, Geva T, Maier SE, Powell AJ. Effect of acquisition parameters on the accuracy of velocity encoded cine magnetic resonance imaging blood flow measurements. J Magn Reson Imaging 2002;15:47-54. [Crossref] [PubMed]
- Promelle V, Bouzerar R, Milazzo S, Balédent O. Quantification of blood flow in the superior ophthalmic vein using phase contrast magnetic resonance imaging. Exp Eye Res 2018;176:40-5. [Crossref] [PubMed]
- Ambarki K, Hallberg P, Jóhannesson G, Lindén C, Zarrinkoob L, Wåhlin A, Birgander R, Malm J, Eklund A. Blood flow of ophthalmic artery in healthy individuals determined by phase-contrast magnetic resonance imaging. Invest Ophthalmol Vis Sci 2013;54:2738-45. [Crossref] [PubMed]
- Yoneyama M, Zhang S, Hu HH, Chong LR, Bardo D, Miller JH, Toyonari N, Katahira K, Katsumata Y, Pokorney A, Ng CK, Kouwenhoven M, Van Cauteren M. Free-breathing non-contrast-enhanced flow-independent MR angiography using magnetization-prepared 3D non-balanced dual-echo Dixon method: A feasibility study at 3 Tesla. Magn Reson Imaging 2019;63:137-46. [Crossref] [PubMed]
- Seeley MA, Georgiadis AG, Sankar WN. Hip Vascularity: A Review of the Anatomy and Clinical Implications. J Am Acad Orthop Surg 2016;24:515-26. [Crossref] [PubMed]
- Xiao J, Yang XJ, Xiao XS. DSA observation of hemodynamic response of femoral head with femoral neck fracture during traction: a pilot study. J Orthop Trauma 2012;26:407-13. [Crossref] [PubMed]
- Zlotorowicz M, Szczodry M, Czubak J, Ciszek B. Anatomy of the medial femoral circumflex artery with respect to the vascularity of the femoral head. J Bone Joint Surg Br 2011;93:1471-4. [Crossref] [PubMed]
- Lazaro LE, Klinger CE, Sculco PK, Helfet DL, Lorich DG. The terminal branches of the medial femoral circumflex artery: the arterial supply of the femoral head. Bone Joint J 2015;97-B:1204-13. [Crossref] [PubMed]
- Netter FH. Atlas of human anatomy, 7th edn. Philadelphia, PA: Elsevier; 2019.
- Caggiati A, Bergan JJ, Gloviczki P, Eklof B, Allegra C, Partsch H. Nomenclature of the veins of the lower limb: extensions, refinements, and clinical application. J Vasc Surg 2005;41:719-24. [Crossref] [PubMed]
- Morgan AG, Thrippleton MJ, Stringer M, Jin N, Wardlaw JM, Marshall I. Repeatability and comparison of 2D and 4D flow MRI measurement of intracranial blood flow and pulsatility in healthy individuals and patients with cerebral small vessel disease. Front Psychol 2023;14:1125038. [Crossref] [PubMed]
- Chen X, Qin L, Pan D, Huang Y, Yan L, Wang G, Liu Y, Liang C, Liu Z. Liver diffusion-weighted MR imaging: reproducibility comparison of ADC measurements obtained with multiple breath-hold, free-breathing, respiratory-triggered, and navigator-triggered techniques. Radiology 2014;271:113-25. [Crossref] [PubMed]
- Chang KH, Lee YH, Chen CY, Lin MF, Lin YC, Chen JH, Chan WP. Inter- and Intra-Rater Reliability of Individual Cerebral Blood Flow Measured by Quantitative Vessel-Flow Phase-Contrast MRI. J Clin Med 2020;9:3099. [Crossref] [PubMed]
- Dyvorne HA, Knight-Greenfield A, Besa C, Cooper N, Garcia-Flores J, Schiano TD, Markl M, Taouli B. Quantification of hepatic blood flow using a high-resolution phase-contrast MRI sequence with compressed sensing acceleration. AJR Am J Roentgenol 2015;204:510-8. [Crossref] [PubMed]
- Zlotorowicz M, Czubak-Wrzosek M, Wrzosek P, Czubak J. The origin of the medial femoral circumflex artery, lateral femoral circumflex artery and obturator artery. Surg Radiol Anat 2018;40:515-20. [Crossref] [PubMed]
- Koerte I, Haberl C, Schmidt M, Pomschar A, Lee S, Rapp P, Steffinger D, Tain RW, Alperin N, Ertl-Wagner B. Inter- and intra-rater reliability of blood and cerebrospinal fluid flow quantification by phase-contrast MRI. J Magn Reson Imaging 2013;38:655-62. [Crossref] [PubMed]
- Weinstein RS. Glucocorticoid-induced osteoporosis and osteonecrosis. Endocrinol Metab Clin North Am 2012;41:595-611. [Crossref] [PubMed]
- Griffith JF, Yeung DK, Tsang PH, Choi KC, Kwok TC, Ahuja AT, Leung KS, Leung PC. Compromised bone marrow perfusion in osteoporosis. J Bone Miner Res 2008;23:1068-75. [Crossref] [PubMed]
- Atsumi T, Kuroki Y. Role of impairment of blood supply of the femoral head in the pathogenesis of idiopathic osteonecrosis. Clin Orthop Relat Res 1992;22-30.
- Atsumi T, Kuroki Y, Yamano K, Muraki M. Revascularization in nontraumatic osteonecrosis of the femoral head. Clin Orthop Relat Res 1996;168-73. [Crossref] [PubMed]
- Nakamura J, Ohtori S, Watanabe A, Nakagawa K, Inoue G, Kishida S, Harada Y, Suzuki M, Takahashi K. Recovery of the blood flow around the femoral head during early corticosteroid therapy: dynamic magnetic resonance imaging in systemic lupus erythematosus patients. Lupus 2012;21:264-70. [Crossref] [PubMed]
- Wang G, Zhang CQ, Sun Y, Feng Y, Chen SB, Cheng XG, Zeng BF. Changes in femoral head blood supply and vascular endothelial growth factor in rabbits with steroid-induced osteonecrosis. J Int Med Res 2010;38:1060-9. [Crossref] [PubMed]
- Bluemke DA, Petri M, Zerhouni EA. Femoral head perfusion and composition: MR imaging and spectroscopic evaluation of patients with systemic lupus erythematosus and at risk for avascular necrosis. Radiology 1995;197:433-8. [Crossref] [PubMed]
- Hirata T, Konishiike T, Kawai A, Sato T, Inoue H. Dynamic magnetic resonance imaging of femoral head perfusion in femoral neck fracture. Clin Orthop Relat Res 2001;294-301. [Crossref] [PubMed]

