Quantitative dynamic contrast-enhanced magnetic resonance imaging for renal perfusion measurement in autosomal dominant polycystic kidney disease
Introduction
Autosomal dominant polycystic kidney disease (ADPKD) is a progressive, inherited disorder and a leading cause of end-stage renal disease. It is characterized by the gradual accumulation of renal cysts, leading to increased total kidney volume (TKV), loss of renal function, and, ultimately, kidney failure (1). Height-adjusted TKV (htTKV), an Food and Drug Administration (FDA)-approved imaging biomarker, has been widely used to predict renal outcomes in ADPKD (2,3). However, despite its clinical utility, htTKV exhibits only a moderate correlation with long-term decline in glomerular filtration rate (GFR). Besides, htTKV accounted for only 35% and 22% of individual GFR decline in the CRISP and HALT-PKD cohorts, respectively (4,5), which suggests that htTKV captures only a subset of the complex pathological mechanisms driving kidney function loss in ADPKD. Furthermore, the htTKV trajectory often fluctuates and may require several years of follow-up to detect meaningful change, limiting its responsiveness as a short-term surrogate endpoint. Therefore, reliable, quantitative biomarkers are urgently needed to monitor early disease progression, guide individualized therapy, and optimize patient selection for clinical trials.
Alternative imaging biomarkers have been proposed, such as renal texture analysis, but these are highly sensitive to scanner-specific parameters, image quality, and post-processing techniques (4). Texture metrics are also not directly linked to underlying pathophysiological processes, raising concerns about their biological relevance. Emerging magnetic resonance imaging (MRI) techniques, such as blood oxygen level-dependent (BOLD) imaging, magnetization transfer imaging (MTI), intravoxel incoherent motion (IVIM), and chemical exchange saturation transfer (CEST), offer insights into renal physiology (6-9). Nonetheless, each method has its own limitations related to sensitivity, reproducibility, or physiological specificity, and none currently provides a fully validated quantitative biomarker for monitoring disease progression in ADPKD.
Given these limitations, intrarenal perfusion has emerged as a promising, biologically meaningful biomarker of ADPKD progression. Patients with ADPKD exhibit both endothelial and vascular smooth-muscle dysfunction, which contributes to impaired renal blood flow (RBF), a known independent predictor of renal outcomes in ADPKD (10,11). Notably, RBF improves existing multivariable models for GFR slope prediction, even when adjusted for htTKV and PKD genotype. Contrast-enhanced ultrasound (CEUS) has been proposed for measuring RBF; however, its relatively narrow field of view compared with CT or MRI and its dependence on operator expertise present limitations (12,13). Moreover, obtaining absolute RBF values is challenging, and CEUS is typically limited to semi-quantitative assessment (14). Dynamic CT can measure RBF (15), but concerns regarding radiation exposure and the use of nephrotoxic iodinated contrast agents remain. Arterial spin labelling (ASL) MRI enables non-contrast, quantitative assessment of kidney perfusion but is limited by low signal-to-noise ratio, motion sensitivity, and variable reproducibility (16).
Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) provides a noninvasive means of quantifying renal perfusion by tracking the kinetics of contrast agent uptake (17). This technique provides quantitative pharmacokinetic (PK) parameters such as the volume transfer constant (Ktrans), which reflects capillary-level tissue perfusion (18-21). When small (<1 kDa), extracellular gadolinium-based contrast agents are used, Ktrans provides a reliable estimation of perfusion. Moreover, macrocyclic gadolinium agents are considered safe even in patients with advanced chronic kidney disease, including those with stage 4 or 5 disease (22). In a study of over 6,000 patients receiving up to triple doses of gadoteridol, no serious adverse events were observed (23).
However, a major challenge in the broader application of quantitative DCE-MRI is the variability in measurements at different times due to hardware instability (24-26) and across scanners due to differences in vendor-specific hardware, pulse sequences, and image reconstruction pipelines (27). External calibration phantoms containing known concentrations of contrast agents have been proposed to detect and correct scanner-specific errors. Yet, static phantoms do not accurately model the dynamic signal changes produced by contrast agent movement through living tissue (28,29).
We developed a point-of-care portable perfusion phantom (P4) to detect and correct scanner-specific bias in quantitative DCE-MRI measurements (30,31). The P4 phantom can be scanned alongside the patient without modification to standard imaging protocols, enabling real-time quality assurance. In a recent multicenter study, P4-based correction improved data agreement [i.e., intraclass correlation coefficient (ICC)] across three MRI scanners for Ktrans measurements, increasing it from 0.38 to 0.99 (31). The P4 device is low-cost, user-friendly, and compatible with any commercial MRI system, making it well-suited for multi-institutional studies.
This study aimed to verify that intrarenal Ktrans measurements obtained using quantitative DCE-MRI are highly reproducible when P4-based error correction is applied, and that the corrected Ktrans serves as a reliable biomarker for assessing disease severity and progression in patients with ADPKD. Additionally, other PK parameters were evaluated and compared to Ktrans in terms of their diagnostic accuracy for ADPKD prognosis. We present this article in accordance with the STROBE reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-1764/rc).
Methods
Human subjects
The protocols comply with current regulatory requirements for studies involving human subjects and are approved by the IRB of the University of Alabama at Birmingham (IRB-300011141). This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. All participants were fully informed about the experimental nature of the study and provided written informed consent before the scans. Five healthy subjects were recruited from October 2020 to November 2020, and twenty ADPKD patients were recruited from March 2024 to March 2025 at the University of Alabama at Birmingham (UAB). Table 1 summarizes the clinical profiles of all participants. The healthy subject cohort consisted of five individuals [all females; 3 White, 2 Black; aged 23–41 years; estimated glomerular filtration rate (eGFR) >90 mL/min/1.73 m2]. The patient cohort consisted of twenty individuals (4 males, 16 females; 15 White, 4 Black, 1 Asian; aged 20–46 years), 10 with mild ADPKD (eGFR ≥60 mL/min/1.73 m2 and htTKV ≤750 mL/m) and 10 with severe ADPKD (eGFR <60 mL/min/1.73 m2 or htTKV >750 mL/m). Mayo Imaging Classification (MIC), a method for assessing ADPKD severity, incorporates age as a factor (1). As a result, one patient in the mild ADPKD group was classified as 1D. In healthy subjects, three imaging sessions were conducted on three different 3T MRI scanners (GE Signa, Chicago, USA; Siemens Prisma, Erlangen, Germany; and Philips Ingenia, Best, The Netherlands) within a week to assess the reproducibility of PK parameter measurements across scanners. For patients, two imaging sessions were conducted on a 3T GE Signa MRI scanner, spaced 4.2±2.5 days apart (range, 1–9 days). However, the second DCE-MRI scan for a patient with mild ADPKD (subject 11 in Table 1) was unsuccessful due to improper catheter placement. In addition, the first DCE-MRI scan for a patient with severe ADPKD (subject 16 in Table 1) was terminated at the subject’s request because of the cold sensation from the contrast injection. Participants fasted for at least 12 hours before imaging to minimize the effects of salt intake on RBF (32) and to reduce abdominal peristaltic motion artifacts in MRI. Also, the participants were refrained from consuming caffeine and alcohol for 24 hours before imaging, as these substances can affect kidney perfusion.
Table 1
| Group | Subject ID | Age (years) | Sex | Race | BW (kg) | Height (m) | eGFR (mL/min/1.73 m2) | htTKV (mL/m) |
TCV (mL) |
MIC |
|---|---|---|---|---|---|---|---|---|---|---|
| Healthy subjects | 1 | 34 | Female | White | 76 | 1.68 | 96 | 191 | NA | NA |
| 2 | 29 | Female | White | 81 | 1.65 | 125 | 214 | NA | NA | |
| 3 | 23 | Female | Black | 68 | 1.47 | >90 | 184 | NA | NA | |
| 4 | 26 | Female | White | 59 | 1.6 | >90 | 144 | NA | NA | |
| 5 | 41 | Female | Black | 66 | 1.63 | 90 | 192 | NA | NA | |
| Mild ADPKD | 6 | 46 | Male | White | 92 | 1.85 | 76 | 598 | 617 | 1C |
| 7 | 46 | Female | White | 108 | 1.65 | 108 | 545 | 207 | 1B | |
| 8 | 39 | Female | White | 48 | 1.57 | 96 | 713 | 548 | 1C | |
| 9 | 43 | Female | White | 103 | 1.68 | 94 | 325 | 151 | 1B | |
| 10 | 23 | Male | Black | 74 | 1.8 | 108 | 371 | 144 | 1C | |
| 11 | 22 | Female | White | 62 | 1.63 | 93 | 228 | 14 | 1B | |
| 12 | 41 | Female | White | 112 | 1.88 | 82 | 526 | 378 | 1C | |
| 13 | 30 | Female | White | 75 | 1.66 | 102 | 252 | 43 | 1B | |
| 14 | 20 | Female | White | 85 | 1.7 | 94 | 370 | 140 | 1D | |
| 15 | 26 | Female | White | 104 | 1.65 | 127 | 272 | 26 | 1B | |
| Severe ADPKD | 16 | 44 | Male | Asian | 129 | 1.8 | 41 | 1,644 | 2,294 | 1D |
| 17 | 40 | Female | White | 111 | 1.63 | 59 | 2,038 | 2,910 | 1E | |
| 18 | 40 | Female | White | 80 | 1.7 | 36 | 737 | 677 | 1C | |
| 19 | 41 | Female | Black | 85 | 1.75 | 48 | 1,711 | 2,374 | 1E | |
| 20 | 36 | Male | White | 104 | 1.75 | 53 | 2,252 | 2,841 | 1E | |
| 21 | 40 | Female | Black | 84 | 1.65 | 66 | 1,517 | 1,704 | 1D | |
| 22 | 33 | Female | White | 98 | 1.68 | 117 | 999 | 827 | 1D | |
| 23 | 38 | Female | White | 109 | 1.7 | 59 | 658 | 668 | 1C | |
| 24 | 24 | Female | Black | 53 | 1.6 | 92 | 1,207 | 943 | 1E | |
| 25 | 45 | Female | White | 102 | 1.7 | 52 | 1,222 | 1,704 | 1D |
ADPKD, autosomal-dominant polycystic kidney disease; BW, body weight; eGFR, estimated glomerular filtration rate; htTKV, height-adjusted total kidney volume; MIC, Mayo imaging classification; NA, not applicable; TCV, total cyst volume.
MRI parameters
Participants were positioned on the phantom package containing three P4s with a torso-phased array coil placed around the abdomen as demonstrated in a previous study (31). For healthy subjects, the MRI parameters and protocols are the same as described in a previous study (31). For kidney and cyst volume measurement of patient cohort, T2-weighted (T2W) images were acquired with single shot fast spin echo sequence (SSFSE) in the coronal plane with the following parameters: field of view =44–48 cm, matrix size =512×512, slice thickness =4–5 mm, number of slices =40–57, repetition time (TR) =605–944 ms, echo time (TE) =70–90 ms, and number of averages =0.6–0.7. Additionally, AIR Recon DL (GE Healthcare, Waukesha, WI, USA) was employed for denoising of T2W images (33). Prior to DCE-MRI, T1-weighted (T1W) images were obtained using a fast spoiled gradient echo sequence (FSPGR) with four flip angles (FAs) (2°, 5°, and 10°) for T1 mapping (34). The imaging protocol included a field of view of 40×40 cm, with phase and frequency encoding set to 160 and 192, respectively. A total of 12 slices were acquired using a matrix of 256×256, with a slice thickness of 5 mm and no interslice gap. The acquisition was performed with a TR of 3.9 ms and a TE of 2.1 ms. The scan utilized a single average and a SENSE acceleration factor of 2, achieving a temporal resolution of 2.9 seconds. T1-weighted images were acquired over a 30-second free-breathing period, during which images captured in the expiratory phase were automatically identified and averaged for each FA. DCE-MRI was also performed in free-breathing mode using the same imaging sequence and parameters as for T1W imaging, with a FA of 15°. DCE-MRI was performed over approximately 9 minutes. Intravenous injection of gadoteridol (0.1 mmol/kg), a macrocyclic contrast agent approved by the FDA, was initiated 30 seconds after the start of imaging, followed by a 20-mL saline flush delivered at a steady rate of 2 mL/s. For phantom imaging, each of the three P4 phantoms received a 4-mL injection of 100 mM gadoteridol at a rate of 0.24 mL/s, beginning 15 seconds after image acquisition commenced, using a syringe pump (NE-1600).
Image processing
The automatic segmentation of kidneys and cysts was performed using a lab-trained nnUNet model (35), trained on 604 3D volumes from the CRISP dataset, tested on our study subjects, and compared against semi-automatically generated ground truth. Training incorporated extensive data augmentation, including random rotations, scaling, intensity shifts, and elastic deformations, to improve robustness to anatomical variability and imaging artifacts. The model was optimized using a combination of weighted Dice and Jaccard losses over 1,000 epochs with a learning rate schedule. The kidney and cyst regions in ten patients with mild (n=7) and severe (n=3) ADPKD were segmented using a semi-automatic method described in a previous study (35), which served as the ground truth for evaluating the segmentation accuracy of the nnUNet model. Model performance was assessed using the Dice similarity coefficient (DSC). The processing of DCE-MRI data to generate PK maps involved seven sequential steps. First, images were spatially aligned using a B-spline registration method based on the expiratory phase (36). Second, local FA variations were quantified via B1 mapping. Third, T1 relaxation maps were generated using the variable flip angle (VFA) method (34). Fourth, a look-up table was constructed to relate the reference contrast enhancement coefficients (CECs) of P4 phantoms to those measured by MRI (30). Fifth, maps of contrast agent concentration (CAC) in tissues were computed (37), with FA inhomogeneity corrected using the B1 maps from step two. Sixth, CAC maps were adjusted using the established look-up table from step four. Lastly, intrarenal PK parameter maps were derived using a population-based arterial input function (38) to standardize outputs and reduce variability. All image processing steps were performed using a lab-made software package developed in LabVIEW v17.0, with the SSM-based PK mapping module implemented in MATLAB v2020a.
PK parameter mapping
Eleven PK maps were generated using three modeling approaches: the extended Tofts model (ETM), which provided Ktrans, kep, ve, and vp (39); the Tofts model (TM), which yielded Ktrans, kep, and ve (40); and the Shutter Speed Model (SSM), which produced Ktrans, kep, ve, and τi (41). The detailed equations for each model are provided in a previous publication (27). In these models, Ktrans represents the volume transfer constant, kep is the flux rate constant, ve denotes the extravascular extracellular volume fraction, vp is the fractional plasma volume, and τi is the mean intracellular water lifetime.
Statistical analysis
The reproducibility (or repeatability) of PK parameter measurement was assessed using the within-subject coefficient of variation (wCV) or ICC. Comparisons of intrarenal PK parameters between mild and severe ADPKD groups were performed using analysis of variance (ANOVA). Correlations between clinical indicators (e.g., htTKV, total cyst volume (TCV), and eGFR) and intrarenal PK metrics were evaluated using Pearson correlation coefficient. The sensitivity, specificity, and overall accuracy of each PK parameter in differentiating mild ADPKD from severe ADPKD were assessed using receiver operating characteristic (ROC) curve analysis (42). Optimal cutoff values were determined using the Youden index (43). All data in this manuscript are presented as mean ± standard deviation (SD). Statistical significance was defined as a P value of less than 0.05. Statistical analyses were performed using SAS software (version 9.4; SAS Institute Inc., Cary, NC, USA).
Results
AI-generated segmentations exhibited excellent agreement with the semi-automatically derived ground truth, yielding DSCs of 91% for kidney regions and 98% for cyst regions (Figure 1). The reproducibility (wCV) of renal Ktrans measurement across three scanners was 50–80%, while the repeatability measured on a single scanner was 20–25% before P4-based error correction (Table 2). However, after applying P4-based error correction, reproducibility and repeatability improved to comparable levels. The P4-based error correction enhanced the repeatability of ETM-derived Ktrans measurements by approximately threefold, while the repeatability of ve and vp remained largely unchanged across all PK models. Figure 2A shows representative ETM-based Ktrans maps of a patient with ADPKD, acquired on days 0 and 2, illustrating the visual impact of P4-based correction. Figure 2B shows the scatter plot of ETM-based Ktrans values from two repeated measurements (Scan 1 and Scan 2) in ADPKD patients, demonstrating improved agreement after P4-based correction. TM-derived Ktrans repeatability improved twofold following correction, with minimal changes observed for ve and kep. In contrast, SSM-derived repeatability metrics for Ktrans and kep showed no meaningful improvement, and ve and τi remained stable pre- and post-correction.
Table 2
| wCV | ETM | TM | SSM | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Ktrans | ve | kep | vp | Ktrans | ve | kep | Ktrans | ve | kep | τi | |||
| Reproducibility (healthy subjects), % | |||||||||||||
| Before correction | 52 | 38 | 36 | 70 | 58 | 36 | 30 | 78 | 64 | 42 | 74 | ||
| After correction | 11 | 14 | 16 | 39 | 16 | 14 | 16 | 52 | 21 | 42 | 27 | ||
| Repeatability (ADPKD patients), % | |||||||||||||
| Before correction | 21 | 11 | 12 | 48 | 22 | 11 | 14 | 25 | 9 | 24 | 23 | ||
| After correction | 8 | 9 | 9 | 49 | 10 | 8 | 11 | 36 | 8 | 35 | 17 | ||
kep, blood influx rate; Ktrans, volume transfer constant; ve, extravascular extracellular space; vp, fractional plasma volume; τi, intracellular lifetime of water. ADPKD, autosomal-dominant polycystic kidney disease; DCE-MRI, dynamic contrast-enhanced magnetic resonance imaging; ETM, extended Tofts model; PK, pharmacokinetic; SSM, Shutter Speed Model; TM, Tofts model; wCV, within-subject coefficient of variation.
Following P4-based error correction, Ktrans and ve values were significantly higher in the mild ADPKD group compared with the severe group, regardless of the PK model used (Table 3). The vp and τi values also remained significantly higher in the mild group. However, kep derived from the TM and ETM models showed no significant difference between groups after correction. The renal Ktrans values of healthy subjects were approximately 25–60% higher than those of patients with mild ADPKD, although the difference did not reach statistical significance (P>0.05). In Figure 3, the representative ETM-based PK maps (Ktrans, ve, and vp) from a healthy subject and two patients with mild and severe ADPKD illustrate these differences.
Table 3
| PK parameters | ETM | TM | SSM | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Ktrans (min−1) | ve | kep (min−1) | vp | Ktrans (min−1) | ve | kep (min−1) | Ktrans (min−1) | ve | kep (min−1) | τi (s) | |||
| Healthy subjects | |||||||||||||
| Before correction | 0.33±0.15 | 0.73±0.23 | 0.49±0.27 | 0.07±0.05 | 0.41±0.21 | 0.74±0.22 | 0.53±0.14 | 0.48±0.37 | 0.40±0.22 | 1.06±0.56 | 0.19±0.13 | ||
| After correction | 0.21±0.05 | 0.57±0.14 | 0.38±0.07 | 0.05±0.02 | 0.27±0.07 | 0.60±0.15 | 0.46±0.08 | 0.68±0.38 | 0.51±0.14 | 1.28±0.61 | 0.31±0.13 | ||
| ADPKD patients | |||||||||||||
| Before correction | |||||||||||||
| Mild | 0.41±0.14 | 0.79±0.10 | 0.51±0.13 | 0.06±0.04 | 0.49±0.18 | 0.80±0.11 | 0.60±0.17 | 0.55±0.17 | 0.79±0.09 | 0.69±0.23 | 0.16±0.06 | ||
| Severe | 0.29±0.07 | 0.63±0.09 | 0.46±0.06 | 0.03±0.01 | 0.33±0.08 | 0.63±0.09 | 0.51±0.07 | 0.39±0.11 | 0.65±0.07 | 0.59±0.15 | 0.26±0.09 | ||
| P value | 0.003 | <0.001 | 0.175 | 0.001 | 0.001 | <0.001 | 0.033 | 0.002 | <0.001 | 0.107 | 0.001 | ||
| After correction | |||||||||||||
| Mild | 0.17±0.04 | 0.46±0.11 | 0.38±0.09 | 0.03±0.01 | 0.20±0.05 | 0.47±0.11 | 0.44±0.11 | 0.43±0.19 | 0.54±0.10 | 0.82±0.40 | 0.38±0.08 | ||
| Severe | 0.09±0.02 | 0.25±0.05 | 0.39±0.05 | 0.01±0.01 | 0.10±0.02 | 0.25±0.05 | 0.42±0.06 | 0.19±0.08 | 0.32±0.07 | 0.60±0.23 | 0.54±0.10 | ||
| P value | <0.001 | <0.001 | 0.733 | <0.001 | <0.001 | <0.001 | 0.66 | <0.001 | <0.001 | 0.049 | <0.001 | ||
Data are presented as mean ± standard deviation. kep, blood influx rate; Ktrans, volume transfer constant; ve, extravascular extracellular space; vp, fractional plasma volume; τi, intracellular lifetime of water. ADPKD, autosomal-dominant polycystic kidney disease; DCE-MRI, dynamic contrast-enhanced magnetic resonance imaging; ETM, extended Tofts model; PK, pharmacokinetic; SSM, Shutter Speed Model; TM, Tofts model.
ETM-based Ktrans achieved the highest accuracy (95%) in distinguishing between mild and severe ADPKD following error correction, compared with only 50% before correction (Table 4). Accuracy for ve improved from 70–80% to 85–90% across all PK models after correction. TM-based and SSM-based Ktrans accuracy increased from 60% and 55% to 90% and 70%, respectively, whereas kep accuracy remained unchanged. Accuracy for vp rose from 50% to 70%, and for τi from 60% to 75% after correction.
Table 4
| Performance | ETM | TM | SSM | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Ktran | ve | kep | vp | Ktrans | ve | kep | Ktrans | ve | kep | τi | |||
| Before correction | |||||||||||||
| Threshold | 0.35 | 0.73 | 0.51 | 0.05 | 0.42 | 0.73 | 0.6 | 0.47 | 0.74 | 0.75 | 0.19 | ||
| Sensitivity, % | 70 | 90 | 70 | 90 | 80 | 80 | 80 | 80 | 90 | 80 | 70 | ||
| Specificity, % | 60 | 60 | 30 | 10 | 40 | 60 | 30 | 30 | 70 | 20 | 50 | ||
| Accuracy, % | 50 | 75 | 50 | 50 | 60 | 70 | 55 | 55 | 80 | 50 | 60 | ||
| After correction | |||||||||||||
| Threshold | 0.13 | 0.31 | 0.39 | 0.02 | 0.15 | 0.35 | 0.45 | 0.34 | 0.45 | 0.73 | 0.37 | ||
| Sensitivity, % | 100 | 80 | 40 | 100 | 100 | 100 | 80 | 100 | 100 | 80 | 100 | ||
| Specificity, % | 90 | 90 | 50 | 40 | 80 | 80 | 30 | 40 | 90 | 20 | 50 | ||
| Accuracy, % | 95 | 85 | 45 | 70 | 90 | 90 | 55 | 70 | 95 | 50 | 75 | ||
kep, blood influx rate; Ktrans, volume transfer constant; ve, extravascular extracellular space; vp, fractional plasma volume; τi, intracellular lifetime of water. ADPKD, autosomal-dominant polycystic kidney disease; DCE-MRI, dynamic contrast-enhanced magnetic resonance imaging; ETM, extended Tofts model; SSM, Shutter Speed Model; TM, Tofts model.
Significant correlations between PK parameters (Ktrans, ve, vp, and τi) and htTKV, TCV, and eGFR were observed across all PK models, even before correction (Table 5). These correlations were further strengthened for Ktrans, vp, and τi after correction. The most notable improvements were observed for ETM-based Ktrans correlations with htTKV (r=−0.68, P=0.001 before vs. r=−0.79, P<0.001 after), TCV (r=−0.69, P=0.001 before vs. r=−0.77, P<0.001 after), and eGFR (r=0.64, P=0.002 before vs. r=0.69, P<0.001 after). Of note, the ETM-based Ktrans after P4-based error correction was significantly correlated with htTKV in healthy subjects and mild ADPKD patients (n=15; r=−0.56, P=0.029) and severe ADPKD patients (n=10; r=−0.64, P=0.048).
Table 5
| Correlation | ETM | TM | SSM | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Ktrans | ve | kep | vp | Ktrans | ve | kep | Ktrans | ve | kep | τi | |||
| Before correction | |||||||||||||
| htTKV | |||||||||||||
| r | −0.68 | −0.85 | −0.42 | −0.72 | −0.70 | −0.85 | −0.53 | −0.69 | −0.83 | −0.40 | 0.77 | ||
| P | 0.001 | <0.001 | 0.062 | <0.001 | 0.001 | <0.001 | 0.016 | 0.001 | <0.001 | 0.082 | <0.001 | ||
| TCV | |||||||||||||
| r | −0.69 | −0.85 | −0.45 | −0.69 | −0.70 | −0.85 | −0.54 | −0.69 | −0.82 | −0.42 | 0.79 | ||
| P | 0.001 | <0.001 | 0.046 | 0.001 | 0.001 | <0.001 | 0.014 | 0.001 | <0.001 | 0.069 | <0.001 | ||
| eGFR | |||||||||||||
| r | 0.64 | 0.70 | 0.46 | 0.69 | 0.67 | 0.71 | 0.55 | 0.65 | 0.70 | 0.42 | −0.66 | ||
| P | 0.002 | <0.001 | 0.040 | <0.001 | 0.001 | <0.001 | 0.012 | 0.002 | <0.001 | 0.064 | 0.002 | ||
| After correction | |||||||||||||
| htTKV | |||||||||||||
| r | −0.79 | −0.74 | −0.08 | −0.73 | −0.78 | −0.75 | −0.21 | −0.72 | −0.77 | −0.46 | 0.75 | ||
| P | <0.001 | <0.001 | 0.733 | <0.001 | <0.001 | <0.001 | 0.381 | <0.001 | <0.001 | 0.041 | <0.001 | ||
| TCV | |||||||||||||
| r | −0.77 | −0.70 | −0.14 | −0.68 | −0.76 | −0.72 | −0.25 | −0.70 | −0.73 | −0.48 | 0.74 | ||
| P | <0.001 | <0.001 | 0.549 | 0.001 | <0.001 | <0.001 | 0.279 | 0.001 | <0.001 | 0.032 | <0.001 | ||
| eGFR | |||||||||||||
| r | 0.69 | 0.62 | 0.17 | 0.65 | 0.70 | 0.63 | 0.29 | 0.61 | 0.62 | 0.45 | −0.67 | ||
| P | <0.001 | 0.004 | 0.463 | 0.002 | <0.001 | 0.003 | 0.208 | 0.005 | 0.003 | 0.049 | 0.001 | ||
kep, blood influx rate; Ktrans, volume transfer constant; ve, extravascular extracellular space; vp, fractional plasma volume; τi, intracellular lifetime of water. DCE-MRI, dynamic contrast-enhanced magnetic resonance imaging; eGFR, estimated glomerular filtration rate; ETM, extended Tofts model; htTKV, height-adjusted total kidney volume; SSM, Shutter Speed Model; TCV, total cyst volume; TM, Tofts model.
Discussion
Based on the PK models, Ktrans measurements from quantitative DCE-MRI scans distinguished between mild and severe ADPKD and were strongly correlated with htTKV, TCV, and eGFR in both groups. As the disease progresses, microcysts form within the renal parenchyma, which are often undetectable by conventional MRI due to its limited spatial resolution. These microcysts reduce ve and concurrently diminish the vascular fraction (vp), which delivers oxygen and nutrients. This reduction in vascular support results in decreased tissue perfusion, as reflected by lower Ktrans values. Notably, because both ve and Ktrans decrease in tandem, the flux rate constant, kep=Ktrans/ve, remains relatively stable throughout disease progression. A similar trend was observed between healthy subjects and mild ADPKD patients, although the difference did not reach statistical significance with the current small sample size. These findings suggest that the proposed method may be useful for monitoring early disease progression at the individual level.
Mild ADPKD patients (eGFR ≥60 mL/min/1.73 m2 and htTKV ≤750 mL/m) exhibited a shorter intracellular water lifetime (τi) compared to severe ADPKD patients (eGFR <60 mL/min/1.73 m2 or htTKV >750 mL/m). Recent investigations have identified τᵢ as a surrogate marker of cellular metabolic activity, showing an inverse association with Na+-K+-ATPase (NKA) pump activity. NKA, driven by phosphate derived from ATP hydrolysis, is essential for maintaining both intracellular and extracellular ion gradients, particularly those of potassium and sodium (44-46). Our results support the notion that as ADPKD severity and microcyst burden increase, the amount of functional renal parenchyma decreases, thereby lowering ATP consumption. This reduction likely diminishes NKA activity, thereby prolonging τᵢ.
Our findings also indicate that implementing the P4 method significantly enhanced the repeatability of Ktrans measurements by up to 3-fold, improving differentiation between mild and severe ADPKD. In contrast, the repeatability of kep and ve remained lower than that of Ktrans, even prior to applying the P4-based correction, which aligns with previous observations (47,48). This discrepancy is attributable to the physiological nature of these parameters: kep, representing the rate constant for the return of contrast agent from the extracellular space to plasma, is largely governed by the washout phase, during which contrast concentration changes are relatively limited. In comparison, Ktrans reflects the extravasation rate during the uptake phase, where broader fluctuations in contrast concentration occur, particularly in high-field MRI systems, leading to increased susceptibility to quantification errors. By correcting these fluctuations, the P4 approach effectively reduced variability and improved the reliability of Ktrans measurements. The limited repeatability of vp was mainly attributed to its inherently low signal-to-noise ratio.
Notably, the Ktrans values derived from the SSM were approximately twice as high as those calculated using the TM or ETM after applying P4-based error correction. This discrepancy arises because TM and ETM assume instantaneous water exchange between the intracellular and extracellular compartments, whereas the SSM explicitly accounts for finite transmembrane water-exchange kinetics (49). As a result, TM and ETM tend to underestimate perfusion-related parameters, particularly in renal parenchyma with microcysts. Therefore, SSM-derived Ktrans may more accurately reflect renal perfusion in this patient population.
In this study, although the reproducibility of renal PK parameter measurements in healthy subjects was evaluated across three MRI scanners, the data from ADPKD patients were acquired on a single scanner. Consequently, while the results are promising, they should be interpreted with caution and validated in larger, multicenter studies. Expanding the investigation to include diverse clinical settings would enhance the generalizability of the findings and help identify potential challenges in implementing the P4-phantom protocol across various MRI platforms and institutional workflows.
Conclusions
In conclusion, this study highlights the potential of quantitative DCE-MRI, enhanced by P4-based error correction, to accurately characterize the severity of ADPKD by detecting the changes in renal perfusion caused by microstructural and functional alterations in the renal parenchyma as cysts develop and grow in patients with ADPKD. The P4 method substantially improved the repeatability of PK parameter measurements in the TM and ETM, although its impact was less pronounced for parameters derived from the SSM. Future studies over a longer timeframe, including larger patient cohorts across multiple institutions and different scanners, may further support the use of Ktrans in tandem with htTKV for ADPKD risk stratification.
Acknowledgments
The authors thank Ms. Quenteeria Mooney and Ms. Deja Cunningham for patient recruitment, and Ms. Nicole Haynes and Mr. Nicholas Hatfield for image acquisition.
Footnote
Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-1764/rc
Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-1764/dss
Funding: This study 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-1764/coif). M.M. reports research support from Otsuka, Sanofi, Vertex, AbbVie, and has been a member of advisory panels/boards for Sanofi, Santa Barbara Nutrients, and PKD Foundation, and has provided consultancy for Otsuka, Sanofi, Vertex, AbbVie, and Regulus. H.K. reports funding from the PKD Foundation and patents on the P4-based error correction strategy (US 10,578,702 B2 &US Provisional Application No. 63/394836). The other authors have no 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. All experiments comply with current regulatory requirements, and are approved by the IRB of the University of Alabama at Birmingham (IRB-300011141), and abide by the laws of the United States of America. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. All participants were fully informed about the experimental nature of the study and provided written informed consent before the scans.
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
- Irazabal MV, Rangel LJ, Bergstralh EJ, Osborn SL, Harmon AJ, Sundsbak JL, Bae KT, Chapman AB, Grantham JJ, Mrug M, Hogan MC, El-Zoghby ZM, Harris PC, Erickson BJ, King BF, Torres VE. Imaging classification of autosomal dominant polycystic kidney disease: a simple model for selecting patients for clinical trials. J Am Soc Nephrol 2015;26:160-72. [Crossref] [PubMed]
- Grantham JJ, Torres VE, Chapman AB, Guay-Woodford LM, Bae KT, King BF Jr, Wetzel LH, Baumgarten DA, Kenney PJ, Harris PC, Klahr S, Bennett WM, Hirschman GN, Meyers CM, Zhang X, Zhu F, Miller JP. Volume progression in polycystic kidney disease. N Engl J Med 2006;354:2122-30. [Crossref] [PubMed]
- Chapman AB, Bost JE, Torres VE, Guay-Woodford L, Bae KT, Landsittel D, Li J, King BF, Martin D, Wetzel LH, Lockhart ME, Harris PC, Moxey-Mims M, Flessner M, Bennett WM, Grantham JJ. Kidney volume and functional outcomes in autosomal dominant polycystic kidney disease. Clin J Am Soc Nephrol 2012;7:479-86. [Crossref] [PubMed]
- Kline TL, Korfiatis P, Edwards ME, Bae KT, Yu A, Chapman AB, Mrug M, Grantham JJ, Landsittel D, Bennett WM, King BF, Harris PC, Torres VE, Erickson BJ. Image texture features predict renal function decline in patients with autosomal dominant polycystic kidney disease. Kidney Int 2017;92:1206-16. [Crossref] [PubMed]
- Chapman AB, Torres VE, Perrone RD, Steinman TI, Bae KT, Miller JP, Miskulin DC, Rahbari Oskoui F, Masoumi A, Hogan MC, Winklhofer FT, Braun W, Thompson PA, Meyers CM, Kelleher C, Schrier RW. The HALT polycystic kidney disease trials: design and implementation. Clin J Am Soc Nephrol 2010;5:102-9. [Crossref] [PubMed]
- Kline TL, Edwards ME, Garg I, Irazabal MV, Korfiatis P, Harris PC, King BF, Torres VE, Venkatesh SK, Erickson BJ. Quantitative MRI of kidneys in renal disease. Abdom Radiol (NY) 2018;43:629-38. [Crossref] [PubMed]
- Kline TL, Irazabal MV, Ebrahimi B, Hopp K, Udoji KN, Warner JD, Korfiatis P, Mishra PK, Macura SI, Venkatesh SK, Lerman LO, Harris PC, Torres VE, King BF, Erickson BJ. Utilizing magnetization transfer imaging to investigate tissue remodeling in a murine model of autosomal dominant polycystic kidney disease. Magn Reson Med 2016;75:1466-73. [Crossref] [PubMed]
- Caroli A, Schneider M, Friedli I, Ljimani A, De Seigneux S, Boor P, Gullapudi L, Kazmi I, Mendichovszky IA, Notohamiprodjo M, Selby NM, Thoeny HC, Grenier N, Vallée JP. Diffusion-weighted magnetic resonance imaging to assess diffuse renal pathology: a systematic review and statement paper. Nephrol Dial Transplant 2018;33:ii29-40. [Crossref] [PubMed]
- Wang F, Kopylov D, Zu Z, Takahashi K, Wang S, Quarles CC, Gore JC, Harris RC, Takahashi T. Mapping murine diabetic kidney disease using chemical exchange saturation transfer MRI. Magn Reson Med 2016;76:1531-41. [Crossref] [PubMed]
- Rahbari-Oskoui F, Williams O, Chapman A. Mechanisms and management of hypertension in autosomal dominant polycystic kidney disease. Nephrol Dial Transplant 2014;29:2194-201. [Crossref] [PubMed]
- Dambreville S, Chapman AB, Torres VE, King BF, Wallin AK, Frakes DH, Yoganathan AP, Wijayawardana SR, Easley K, Bae KT, Brummer ME. Renal arterial blood flow measurement by breath-held MRI: Accuracy in phantom scans and reproducibility in healthy subjects. Magn Reson Med 2010;63:940-50. [Crossref] [PubMed]
- David E, Del Gaudio G, Drudi FM, Dolcetti V, Pacini P, Granata A, Pretagostini R, Garofalo M, Basile A, Bellini MI, D'Andrea V, Scaglione M, Barr R, Cantisani V. Contrast Enhanced Ultrasound Compared with MRI and CT in the Evaluation of Post-Renal Transplant Complications. Tomography 2022;8:1704-15. [Crossref] [PubMed]
- Yaming G, Zuoguan C, Yongjun L, Wayne W Z. Renal cortex blood perfusion contrast-enhanced ultrasound: Assisting diagnosis and treatment of renal artery disease. Aging Med (Milton) 2024;7:438-42. [Crossref] [PubMed]
- Kogan P, Johnson KA, Feingold S, Garrett N, Guracar I, Arendshorst WJ, Dayton PA. Validation of dynamic contrast-enhanced ultrasound in rodent kidneys as an absolute quantitative method for measuring blood perfusion. Ultrasound Med Biol 2011;37:900-8. [Crossref] [PubMed]
- Sah A, Gupta A, Garg S, Yadav N, Khan MA, Das CJ. Can quantitative perfusion CT-based biomarkers predict renal cell carcinoma subtypes? Abdom Radiol (NY) 2025;50:2586-94. [Crossref] [PubMed]
- Odudu A, Nery F, Harteveld AA, Evans RG, Pendse D, Buchanan CE, Francis ST, Fernández-Seara MA. Arterial spin labelling MRI to measure renal perfusion: a systematic review and statement paper. Nephrol Dial Transplant 2018;33:ii15-21. [Crossref] [PubMed]
- Shukla-Dave A, Obuchowski NA, Chenevert TL, Jambawalikar S, Schwartz LH, Malyarenko D, Huang W, Noworolski SM, Young RJ, Shiroishi MS, Kim H, Coolens C, Laue H, Chung C, Rosen M, Boss M, Jackson EF. Quantitative imaging biomarkers alliance (QIBA) recommendations for improved precision of DWI and DCE-MRI derived biomarkers in multicenter oncology trials. J Magn Reson Imaging 2019;49:e101-21. [Crossref] [PubMed]
- Barnes SL, Whisenant JG, Loveless ME, Yankeelov TE. Practical dynamic contrast enhanced MRI in small animal models of cancer: data acquisition, data analysis, and interpretation. Pharmaceutics 2012;4:442-78. [Crossref] [PubMed]
- Zhang CC, Yan Z, Giddabasappa A, Lappin PB, Painter CL, Zhang Q, Li G, Goodman J, Simmons B, Pascual B, Lee J, Levkoff T, Nichols T, Xie Z. Comparison of dynamic contrast-enhanced MR, ultrasound and optical imaging modalities to evaluate the antiangiogenic effect of PF-03084014 and sunitinib. Cancer Med 2014;3:462-71. [Crossref] [PubMed]
- Craciunescu OI, Blackwell KL, Jones EL, Macfall JR, Yu D, Vujaskovic Z, Wong TZ, Liotcheva V, Rosen EL, Prosnitz LR, Samulski TV, Dewhirst MW. DCE-MRI parameters have potential to predict response of locally advanced breast cancer patients to neoadjuvant chemotherapy and hyperthermia: a pilot study. Int J Hyperthermia 2009;25:405-15. [Crossref] [PubMed]
- Kim H, Folks KD, Guo L, Sellers JC, Fineberg NS, Stockard CR, Grizzle WE, Buchsbaum DJ, Morgan DE, George JF, Zinn KR. Early therapy evaluation of combined cetuximab and irinotecan in orthotopic pancreatic tumor xenografts by dynamic contrast-enhanced magnetic resonance imaging. Mol Imaging 2011;10:153-67.
- Woolen SA, Shankar PR, Gagnier JJ, MacEachern MP, Singer L, Davenport MS. Risk of Nephrogenic Systemic Fibrosis in Patients With Stage 4 or 5 Chronic Kidney Disease Receiving a Group II Gadolinium-Based Contrast Agent: A Systematic Review and Meta-analysis. JAMA Intern Med 2020;180:223-30. [Crossref] [PubMed]
- Cho SB, Lee AL, Chang HW, Kim KA, Yoo WJ, Yeom JA, Rho MH, Kim SJ, Lim YJ, Han M. Prospective Multicenter Study of the Safety of Gadoteridol in 6163 Patients. J Magn Reson Imaging 2020;51:861-8. [Crossref] [PubMed]
- Eikefjord E, Andersen E, Hodneland E, Hanson EA, Sourbron S, Svarstad E, Lundervold A, Rørvik JT. Dynamic contrast-enhanced MRI measurement of renal function in healthy participants. Acta Radiol 2017;58:748-57. [Crossref] [PubMed]
- Rata M, Collins DJ, Darcy J, Messiou C, Tunariu N, Desouza N, Young H, Leach MO, Orton MR. Assessment of repeatability and treatment response in early phase clinical trials using DCE-MRI: comparison of parametric analysis using MR- and CT-derived arterial input functions. Eur Radiol 2016;26:1991-8. [Crossref] [PubMed]
- Waterton JC, Ho M, Nordenmark LH, Jenkins M, DiCarlo J, Guillard G, Roberts C, Buonaccorsi G, Parker GJM, Bowes MA, Peterfy C, Kellner H, Taylor PC. Repeatability and response to therapy of dynamic contrast-enhanced magnetic resonance imaging biomarkers in rheumatoid arthritis in a large multicentre trial setting. Eur Radiol 2017;27:3662-8. [Crossref] [PubMed]
- Kim H. Variability in Quantitative DCE-MRI: Sources and Solutions. J Nat Sci 2018;4:e484.
- O'Callaghan J, Wells J, Richardson S, Holmes H, Yu Y, Walker-Samuel S, Siow B, Lythgoe MF. Is your system calibrated? MRI gradient system calibration for pre-clinical, high-resolution imaging. PLoS One 2014;9:e96568. [Crossref] [PubMed]
- jackson ef, gupta sn, a. Rm, et al. Qiba dce-mri technical committee update: phantom studies and first dcemri profile. Proceedings of the 96th scientific assembly and annual meeting of the radiological society of north america, chicago, ill, usa, december 2010. 2010.
- Kim H, Mousa M, Schexnailder P, Hergenrother R, Bolding M, Ntsikoussalabongui B, Thomas V, Morgan DE. Portable perfusion phantom for quantitative DCE-MRI of the abdomen. Med Phys 2017;44:5198-209. [Crossref] [PubMed]
- Holland MD, Morales A, Simmons S, Smith B, Misko SR, Jiang X, Hormuth DA, Christenson C, Koomullil RP, Morgan DE, Li Y, Xu J, Yankeelov TE, Kim H. Disposable point-of-care portable perfusion phantom for quantitative DCE-MRI. Med Phys 2022;49:271-81. [Crossref] [PubMed]
- Kurtz TW, DiCarlo SE, Pravenec M, Schmidlin O, Tanaka M, Morris RC Jr. An alternative hypothesis to the widely held view that renal excretion of sodium accounts for resistance to salt-induced hypertension. Kidney Int 2016;90:965-73. [Crossref] [PubMed]
- Brendel JM, Jacoby J, Dehdab R, Herrmann J, Ursprung S, Werner S, Gassenmaier S, Nickel D, Nikolaou K, Afat S, Almansour H. Deep learning reconstruction for accelerated high-resolution upper abdominal MRI improves lesion detection without time penalty. Diagn Interv Imaging 2025;106:85-92. [Crossref] [PubMed]
- Liberman G, Louzoun Y, Ben Bashat D. T1 mapping using variable flip angle SPGR data with flip angle correction. J Magn Reson Imaging 2014;40:171-80. [Crossref] [PubMed]
- Krishnan C, Schmidt E, Onuoha E, Mrug M, Cardenas CE, Kim H. nnUNet for Automatic Kidney and Cyst Segmentation in Autosomal Dominant Polycystic Kidney Disease. Curr Med Imaging 2024;20:e15734056272767. [Crossref] [PubMed]
- Li Z, Tielbeek JAW, Caan MWA, Puylaert CAJ, Ziech MLW, Nio CY, Stoker J, van Vliet LJ, Vos FM. Expiration-phase template-based motion correction of free-breathing abdominal dynamic contrast enhanced MRI. IEEE Trans Biomed Eng 2015;62:1215-25. [Crossref] [PubMed]
- Bokacheva L, Rusinek H, Chen Q, Oesingmann N, Prince C, Kaur M, Kramer E, Lee VS. Quantitative determination of Gd-DTPA concentration in T1-weighted MR renography studies. Magn Reson Med 2007;57:1012-8. [Crossref] [PubMed]
- Parker GJ, Roberts C, Macdonald A, Buonaccorsi GA, Cheung S, Buckley DL, Jackson A, Watson Y, Davies K, Jayson GC. Experimentally-derived functional form for a population-averaged high-temporal-resolution arterial input function for dynamic contrast-enhanced MRI. Magn Reson Med 2006;56:993-1000. [Crossref] [PubMed]
- Tofts PS, Brix G, Buckley DL, Evelhoch JL, Henderson E, Knopp MV, Larsson HB, Lee TY, Mayr NA, Parker GJ, Port RE, Taylor J, Weisskoff RM. Estimating kinetic parameters from dynamic contrast-enhanced T(1)-weighted MRI of a diffusable tracer: standardized quantities and symbols. J Magn Reson Imaging 1999;10:223-32. [Crossref] [PubMed]
- Tofts PS. Modeling tracer kinetics in dynamic Gd-DTPA MR imaging. J Magn Reson Imaging 1997;7:91-101. [Crossref] [PubMed]
- Li X, Huang W, Yankeelov TE, Tudorica A, Rooney WD, Springer CS Jr. Shutter-speed analysis of contrast reagent bolus-tracking data: Preliminary observations in benign and malignant breast disease. Magn Reson Med 2005;53:724-9. [Crossref] [PubMed]
- Fawcett T. An introduction to roc analysis. Pattern Recognition Letters 2006;27:861-74.
- Liu A, Schisterman EF, Zhu Y. On linear combinations of biomarkers to improve diagnostic accuracy. Stat Med 2005;24:37-47. [Crossref] [PubMed]
- Ruggiero MR, Baroni S, Pezzana S, Ferrante G, Geninatti Crich S, Aime S. Evidence for the Role of Intracellular Water Lifetime as a Tumour Biomarker Obtained by In Vivo Field-Cycling Relaxometry. Angew Chem Int Ed Engl 2018;57:7468-72. [Crossref] [PubMed]
- Springer CS Jr, Li X, Tudorica LA, Oh KY, Roy N, Chui SY, Naik AM, Holtorf ML, Afzal A, Rooney WD, Huang W. Intratumor mapping of intracellular water lifetime: metabolic images of breast cancer? NMR Biomed 2014;27:760-73. [Crossref] [PubMed]
- Zhang Y, Poirier-Quinot M, Springer CS Jr, Balschi JA. Active trans-plasma membrane water cycling in yeast is revealed by NMR. Biophys J 2011;101:2833-42. [Crossref] [PubMed]
- Huang W, Chen Y, Fedorov A, Li X, Jajamovich GH, Malyarenko DI, et al. The Impact of Arterial Input Function Determination Variations on Prostate Dynamic Contrast-Enhanced Magnetic Resonance Imaging Pharmacokinetic Modeling: A Multicenter Data Analysis Challenge. Tomography 2016;2:56-66. [Crossref] [PubMed]
- Huang W, Chen Y, Fedorov A, Li X, Jajamovich GH, Malyarenko DI, et al. The Impact of Arterial Input Function Determination Variations on Prostate Dynamic Contrast-Enhanced Magnetic Resonance Imaging Pharmacokinetic Modeling: A Multicenter Data Analysis Challenge, Part II. Tomography 2019;5:99-109. [Crossref] [PubMed]
- Huang W, Li X, Chen Y, Li X, Chang MC, Oborski MJ, et al. Variations of dynamic contrast-enhanced magnetic resonance imaging in evaluation of breast cancer therapy response: a multicenter data analysis challenge. Transl Oncol 2014;7:153-66. [Crossref] [PubMed]

