Capability of arterial spin labeling combined with morphological characteristics to distinguish early- and mid-stage chronic kidney disease: a prospective study on 5 T MRI
Introduction
Chronic kidney disease (CKD), as a major global public health issue, requires optimized early diagnosis and management strategies to critically improve patient outcomes (1-3). Early interventions can significantly delay progression to end-stage renal disease and reduce complications such as cardiovascular events (4,5). The Kidney Disease: Improving Global Outcomes (KDIGO) 2024 Clinical Practice Guideline emphasizes that early-stage CKD (G1–G2 stages) should focus on identifying primary etiologies and preserving renal function, whereas mid-stage CKD (G3A–3B stages) prioritizes slowing the decline rate of estimated glomerular filtration rate (eGFR) and managing complications (6). Previous research demonstrated that the standardized application of guideline-recommended therapies (e.g., renin–angiotensin system inhibitors) in mid-stage CKD patients significantly improves annual eGFR decline rates during a two-year follow-up period (7). This further underscores the clinical value of early screening and timely intervention.
Current CKD diagnosis primarily relies on serum creatinine (Scr), urinary albumin-to-creatinine ratio (UACR), and eGFR, supplemented by renal biopsy histopathological assessment. However, Scr levels are significantly influenced by age, sex, and muscle mass, limiting its sensitivity in detecting early renal impairment (6). Although renal biopsy remains the gold standard, its invasive nature introduces risks such as bleeding and infection, and it is unsuitable for dynamic monitoring (8). Consequently, exploring non-invasive renal function assessment methods with high sensitivity represents a pressing clinical need.
Magnetic resonance imaging (MRI) demonstrates unique advantages in assessing CKD morphology and function through its multi-parametric, multi-sequence imaging capabilities (9-11). Currently, the clinical assessment of renal morphology primarily relies on ultrasound, which is susceptible to interference from intestinal gas and surrounding tissues and has strong operator dependence, thus compromising the reproducibility of morphological evaluations (12). In contrast, MRI enables three-dimensional (3D) visualization of renal anatomy with its multi-parametric and multi-planar imaging capabilities. Particularly at 5 T MRI, higher spatial resolution allows for precise identification of subtle structural alterations in the renal parenchyma and heterogeneity in blood flow perfusion (13,14). Pathological studies have confirmed that CKD progression involves glomerulosclerosis, interstitial fibrosis, and microvascular rarefaction, radiologically manifested as renal volume reduction, cortical thinning, and blurring of corticomedullary differentiation (CMD). Francis et al. (9) and Noda et al. (15) have established that renal volume and cortical thickness serve as independent imaging predictors of disease progression. Arterial spin labeling (ASL) technology non-invasively evaluates renal perfusion by quantifying renal blood flow (RBF) (16,17). Multiple studies have confirmed that RBF is significantly lower in CKD patients than in healthy individuals and is significantly positively correlated with eGFR (18,19). Using a 1.5 T scanner, Rossi et al. demonstrated that RBF was significantly higher in healthy volunteers compared to CKD patients (20). Mao et al., utilizing a 3 T system, further demonstrated a significant positive correlation between RBF and eGFR in CKD patients (21).
Notably, while existing studies predominantly utilize 3 T or 1.5 T MRI, the clinical application of 5 T MRI has enhanced the image quality in renal imaging, demonstrating anatomical and functional abdominal image quality comparable to or better than 3 T (13,14). Its superior spatial resolution enables detailed visualization of subtle renal tissue structures, while also significantly improving the signal-to-noise ratio (SNR) and quantitative reproducibility of ASL imaging. However, to our knowledge, no previous study has specifically focused on the application of 5 T MRI for CKD staging. This study innovatively integrates 5 T ASL perfusion imaging with morphological characteristics. It validates the clinical feasibility of ASL technology at 5 T and then applies it to systematically evaluate the ability of a morpho-functional multiparametric model to distinguish early-stage (G1–G2) CKD from mid-stage (G3A–G3B) CKD. We present this article in accordance with the STARD reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0296/rc).
Methods
Participants
The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. This study was approved by the Ethics Committee of Chongqing University Three Gorges Hospital [No. 2024-KS-(149)]. Written informed consent was provided by all participants prior to enrollment. CKD was defined as the presence of kidney damage (blood or urine composition abnormalities, kidney biopsy findings, radiographic abnormalities) or an eGFR lower than 60 mL/min/1.73 m2, persisting for 3 or more months. From March to December 2024, 103 patients meeting this definition (69 males and 34 females; mean age, 47.29±15.53 years) were prospectively enrolled at Chongqing University Three Gorges Hospital. The study flowchart is shown in Figure 1. The inclusion criteria were as follows: (I) confirmed CKD diagnosis per KDIGO 2024 guidelines (6) with eGFR ≥30 mL/min/1.73 m2; (II) valid SCr and UACR measurements within 7 days before enrollment; (III) no contraindications for MRI examination; and (IV) good compliance with the examination process. The exclusion criteria were as follows: (I) significant kidney asymmetry (>10% inter-kidney parameter discrepancy); (II) renal neoplasms or cysts (>10 mm diameter); (III) history of renal replacement therapy (dialysis/transplantation); (IV) congenital renal anomalies; (V) decompensated heart failure; and (VI) suboptimal imaging quality for quantitative analysis.
The Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation was used to calculate eGFR. The participants were categorized into two groups according to KDIGO 2024 staging criteria (6): the early-stage group (G1–G2) with an eGFR of 60–120 mL/min/1.73 m2, and the mid-stage group (G3A–G3B) with an eGFR of 30–59 mL/min/1.73 m2. Baseline demographic and laboratory parameters are summarized in Table 1.
Table 1
| Characteristics | Early-stage group | Mid-stage group | t/Z/χ2 | P value |
|---|---|---|---|---|
| Number of patients | CKD G1: 36, CKD G2: 30 | CKD G3A–G3B: 37 | – | – |
| Age (years) | 46.00 (32.50, 52.50) | 55.00 (49.00, 67.75) | −4.52 | <0.001 |
| Sex (male/female) | 45/21 | 24/13 | 0.12 | 0.73 |
| Systolic pressure (mmHg) | 125.00 (116.75, 132.00) | 129.00 (117.00, 139.50) | −0.95 | 0.34 |
| Diastolic pressure (mmHg) | 76.50 (70.00, 82.25) | 72.00 (67.50, 82.50) | −0.71 | 0.48 |
| Cause of CKD | 8.45 | 0.07 | ||
| Glomerulonephritis | 32 (48.48) | 13 (35.14) | ||
| Diabetic kidney disease | 4 (6.06) | 3 (8.11) | ||
| Hypertensive nephropathy | 4 (6.06) | 9 (24.32) | ||
| Other | 9 (13.64) | 6 (16.22) | ||
| Unknown etiology | 17 (25.76) | 6 (16.22) | ||
| Scr (μmol/L) | 82.86±18.59 | 131.16±19.98 | −12.32 | <0.001 |
| eGFR (mL/min/1.73 m2) | 93.00 (80.15, 102.60) | 49.65 (42.68, 53.60) | −8.21 | <0.001 |
| UACR (mg/mmol) | 110.22±26.45 | 102.49±27.02 | 0.28 | 0.60 |
Data are presented as n, n (%), median (Q1, Q3) or mean ± SD. CKD, chronic kidney disease; eGFR, estimated glomerular filtration rate; Q1, first quartile; Q3, third quartile; Scr, serum creatinine; SD, standard deviation; UACR, urinary albumin-to-creatinine ratio.
MRI protocols
MRI examinations were performed using a 5 Tesla (5 T) MRI system (uMR Jupiter, United Imaging Healthcare, Shanghai, China). An 8-channel parallel volumetric transmit coil was used to mitigate dielectric artifacts through independent channel control. A custom-built 24-channel body coil was used for all studies at 5 T using local B1 ± shimming for B1 ± optimization. All patients were required to discontinue medications that may affect kidney hemodynamics, such as angiotensin-converting enzyme inhibitors (ACEIs) and angiotensin receptor blockers (ARBs), for at least 12 hours prior to the examination, and to abstain from food and water for at least 6 hours. Coronal breath-hold T1-weighted imaging (T1WI) inversion recovery-fast spoiled gradient echo (IR-GRE-FSP) sequence, sagittal breath-hold T2-weighted imaging (T2WI) single-shot fast spin echo (SSFSE) sequence, and transversal free-breathing T2WI fast spin echo (FSE) sequence with fat suppression were employed to exclude organic renal lesions and to evaluate morphological characteristics. Kidney perfusion was measured using a free-breathing, two-dimensional (2D) echo-planar sequence with pulsed ASL and flow-sensitive alternating inversion recovery (PASL_FAIR) scheme. Retrospective registration of the image volumes was performed prior to averaging to reduce the influence of breathing motion on the perfusion imaging. All the relevant scan parameters are summarized in Table 2.
Table 2
| Parameters | T1WI | T2WI | T2WI | ASL |
|---|---|---|---|---|
| Sequence type | IR-GRE-FSP | SSFSE | FSE | EPI PASL_FAIR |
| Orientation | Coronal | Sagittal | Transversal | Coronal |
| Respiratory control | Breath-holding | Breath-holding | Free-breathing | Free-breathing |
| TR/TE/TI/PLD (ms) | 5.4/2.2/1,100/– | 191/42/–/– | 3,686/100.52/–/– | 4,000/19.3/–/1,800 |
| Voxel size (mm3) | 1.25×1.00×5 | 1.93×1.35×5 | 1.25×1.25×5 | 2.08×2.08×8 |
| FOV (mm2) | 320×320 | 260×200 | 380×380 | 200×300 |
| Matrix | 256×320 | 135×148 | 304×304 | 96×144 |
| Phase oversampling | 0% | 30% | 5% | 10% |
| Slices | 15 | 35 | 28 | 1 |
| FA | 9° | 90° | 90° | 90° |
| Fat suppression | – | – | FS | FS |
| Bandwidth (HZ/pixel) | 300 | 500 | 400 | 2000 |
| Acceleration factor | – | 3 (GRAPPA) | 4 (GRAPPA) | 3 (GRAPPA) |
| Acquisition time | 54.3 s (18.1 s ×3) | 15.8 s | 3 min | 3 min |
ASL, arterial spin labeling; EPI, echo planar imaging; FA, flip angle; FOV, field of view; FS, fat saturation; FSE, fast spin echo; GRAPPA, generalized autocalibration partially parallel acquisition; IR-GRE-FSP, inversion recovery-fast spoiled gradient echo imaging; MRI, magnetic resonance imaging; PASL_FAIR, pulsed arterial spin labeling with flow-sensitive alternating inversion recovery; PLD, post labeling delay; SSFSE, single shot fast spin echo; TE, echo time; TI, inversion time; TR, repetition time; T1WI, T1-weighted imaging; T2WI, T2-weighted imaging.
Image analysis
All image analyses were conducted by two radiologists (J.K., with over 20-years’ experience of abdominal MRI; and W.Z., with more than 30-years’ experience of abdominal MRI) blindly and independently.
Renal morphological assessment
Quantitative image analysis
(I) 3D renal dimensions: according to ultrasonography standards (22,23), renal length was measured at the hilum level on T2WI sagittal images (Figure 2A). Thickness and width were assessed at the maximum transverse section and hilum level on T2WI axial images (Figure 2B,2C). All other morphological evaluations were performed on T1WI coronal images, except for the 3D renal dimensions. (II) Cortical thickness: three-point measurements were taken at the upper, middle, and lower poles of the kidney at the hilum level, and the average cortical thickness was calculated (Figure 2D). (III) Cortical signal heterogeneity: large contiguous regions of interest (ROIs) (110–250 pixels) were manually delineated along the cortical margin at the hilum level, avoiding renal sinus and vascular structures (Figure 2E). Heterogeneity index was calculated as: .
Qualitative image analysis
Given the preserved visibility of CMD in early-to-mid stage CKD patients on T1WI, CMD contrast was categorized into two types: good (CMD is maintained and clearly visible in all segments of the kidney); moderate (CMD is present but not distinct in some parts of the kidney) (24) (Figure 2F,2G).
ASL processing and assessment
Raw Digital Imaging and Communications in Medicine (DICOM) datasets were processed using Matlab R2018b (MathWorks, Natick, MA, USA). The post-processing pipeline consisted of the following steps: (I) initial contour delineation: the user manually outlined the left and right kidney boundaries on an averaged display of all 20 control and label images to define a baseline anatomical reference. (II) Rigid motion correction: All individual control and label images were co-registered to this reference mask using an automated rigid-body registration algorithm to correct for respiratory motion. (III) Iterative mask validation: a quality control loop allowed the user to inspect the kidney contour overlaid on each registered image; if respiratory motion caused the kidney to fall outside the mask, the user corrected the mask for the affected image indices, iterating until all 20 masks were confirmed as anatomically accurate. (IV) Signal averaging: the mean signal intensity within each kidney ROI was calculated for the control, label, and proton density (PD)-weighted images (M0) by averaging across the 20 repetitions to enhance the SNR. (V) RBF calculation and export: The mean label, control, and M0 images obtained from the scan were used to calculate RBF voxel-wise using the following formula, yielding a color-coded RBF parametric map (25):
Based on previous studies and ASL guidelines, the tissue–blood partition coefficient (λ) was set at 0.9 mL/100 g, the labelling efficiency (α) was 0.95, and SI with the respective subscripts represents the signal intensities of the control, label, and PD images. The magnetic resonance (MR) longitudinal relaxation time of arterial blood (T1,blood) was set at 1.65 s. Inflowing suppression of arterial blood and background suppression were both performed. The post-labeling delay (TI) and the bolus duration (TI_1) were both set at 1,800 ms.
Two ROIs (each 110–250 pixels) covering the entire renal cortex for each kidney were manually drawn on the RBF maps, while avoiding areas with cysts, artifacts, the collecting system, or vascular structures (Figure 2H). The MRI parameter values from both kidneys were averaged for subsequent statistical analysis, as CKD typically affects both kidneys without bias and eGFR reflects the overall function of both kidneys.
Statistical analysis
Statistical analysis was performed using the software SPSS 26.0 (IBM Corp., Armonk, NY, USA) and MedCalc version 23.1 (MedCalc, Ostend, Belgium). Continuous variables were first assessed for normality using the Shapiro-Wilk test. Normally distributed data were presented as mean ± standard deviation (SD) and compared between the two groups using independent samples t-test. Non-normally distributed data were presented as median with interquartile range (IQR) [median (Q1, Q3)] and analyzed with Mann-Whitney U test. Categorical variables were presented as frequencies (%) and compared via Pearson’s χ2 test. Positive imaging indicators demonstrating significant intergroup differences (P<0.05) were included in a multivariate logistic regression analysis to identify independent diagnostic markers for CKD. Receiver operating characteristic (ROC) curve analysis was performed to evaluate diagnostic performance of individual characteristics and combined RBF-morphology model, with the area under the curve (AUC) calculated for each. The DeLong test was applied to compare AUC differences among individual characteristics and combined models. Spearman correlation analysis was used to evaluate the correlations between eGFR and imaging characteristics. Intraclass correlation coefficient (ICC) and Cohen’s Kappa were used to assess the interobserver agreement between the two radiologists (a coefficient <0.60 indicated poor agreement, 0.61–0.80 indicated moderate agreement, and >0.80 indicated excellent agreement). A two-tailed P value <0.05 was considered statistically significant.
Results
Patient characteristics
Table 1 summarizes the clinical characteristics of the enrolled CKD patients. A total of 103 patients were included, with 66 (64.08%) classified as early-stage CKD and 37 (35.92%) as mid-stage CKD. The etiologies were as follows: glomerulonephritis in 45 patients, diabetic kidney disease in 7, hypertensive nephropathy in 13, other causes in 15, and unknown etiology in 23. Significant differences were observed between the two groups in SCr and eGFR (both P<0.001), whereas UACR did not differ significantly (P=0.60).
Interobserver agreement of imaging characteristics
Quantitative measurements demonstrated excellent interobserver agreement (Table 3): renal length [ICC =0.938, 95% confidence interval (CI): 0.874–0.970], width (ICC =0.915, 95% CI: 0.829–0.959), thickness (ICC =0.906, 95% CI: 0.813–0.954), cortical thickness (ICC =0.931, 95% CI: 0.862–0.967), cortical signal heterogeneity (ICC =0.893, 95% CI: 0.788–0.948), and RBF values (ICC =0.905, 95% CI: 0.812–0.954). Qualitative assessment of CMD contrast also showed excellent agreement (Kappa =0.870).
Table 3
| Characteristics | Observer_1 | Observer_2 | ICC (95% CI)/Kappa value |
|---|---|---|---|
| Renal length (mm) | 97.53±9.75 | 97.71±10.03 | 0.938 (0.874–0.970) |
| Renal width (mm) | 59.50 (54.75, 63.25) | 59.33 (55.00, 63.00) | 0.915 (0.829–0.959) |
| Renal thickness (mm) | 46.50 (44.75, 49.00) | 46.50 (44.23, 49.63) | 0.906 (0.813–0.954) |
| Cortical thickness (mm) | 5.45 (4.85, 5.90) | 5.50 (4.80, 6.03) | 0.931 (0.862–0.967) |
| Cortical signal heterogeneity (%) | 46.25 (39.19, 55.06) | 46.11 (40.74, 53.88) | 0.893 (0.788–0.948) |
| RBF values (mL/100 g/min) | 206.00 (183.00, 228.60) | 206.65 (190.85, 231.19) | 0.905 (0.812–0.954) |
| CMD contrast | 0.870 | ||
| Good | 26 (86.67) | 25 (83.33) | |
| Moderate | 4 (13.33) | 5 (16.67) |
Data are presented as n (%), median (Q1, Q3) or mean ± SD. CI, confidence interval; CMD, corticomedullary differentiation; ICC, intraclass correlation coefficient; Q1, first quartile; Q3, third quartile; RBF, renal blood flow; SD, standard deviation.
Comparison of imaging characteristics between early- and mid-stage CKD groups
The renal length, width, cortical thickness, cortical signal heterogeneity, CMD contrast, and RBF values showed statistically significant differences between groups. Table 4 shows that renal length decreased significantly in mid-stage CKD (t=2.62, P<0.001), with concomitant reductions in renal width (Z=−2.61, P<0.001) and cortical thickness (Z=−3.69, P<0.001). Cortical signal heterogeneity increased significantly (Z=−4.55, P<0.001), with significant degradation of CMD contrast (χ2=25.97, P<0.001). RBF values were markedly lower in mid-stage patients (Z=−5.07, P<0.001).
Table 4
| Characteristics | Early-stage group | Mid-stage group | t/Z/χ2 | P value |
|---|---|---|---|---|
| Renal length (mm) | 98.64±9.50 | 93.76±8.23 | 2.62 | <0.001 |
| Renal width (mm) | 60.50 (55.50, 65.00) | 58.00 (54.50, 60.88) | −2.61 | <0.001 |
| Renal thickness (mm) | 47.00 (44.32, 50.43) | 46.08 (42.00, 53.00) | −0.65 | 0.52 |
| Cortical thickness (mm) | 5.80 (5.40, 6.10) | 5.20 (4.71, 5.64) | −3.69 | <0.001 |
| CMD contrast | 25.97 | <0.001 | ||
| Good | 63 (95.45) | 20 (54.05) | ||
| Moderate | 3 (4.55) | 17 (45.95) | ||
| Cortical signal heterogeneity (%) | 42.93 (37.27, 47.68) | 51.18 (44.69, 58.58) | −4.55 | <0.001 |
| RBF values (mL/100 g/min) | 217.70 (203.78, 258.53) | 187.53 (170.71, 205.69) | −5.07 | <0.001 |
Data are presented as n (%), median (Q1, Q3) or mean ± SD. CKD, chronic kidney disease; CMD, corticomedullary differentiation; Q1, first quartile; Q3, third quartile; RBF, renal blood flow; SD, standard deviation.
Multivariate logistic regression analysis identified the independent predictors for the diagnosis of early-stage CKD
The results of the multivariate logistic regression analysis are shown in Table 5. Renal length (B=−0.133, P=0.004), CMD contrast (B=2.351, P=0.010), cortical signal heterogeneity (B=0.093, P=0.022), and RBF values (B=−0.038, P=0.002) were identified as independent predictors for distinguishing between early and mid-stage CKD.
Table 5
| Characteristics | B | SE | Walds χ2 | OR (95% CI) | P value |
|---|---|---|---|---|---|
| Renal length | −0.133 | 0.046 | 8.247 | 0.875 (0.799–0.959) | 0.004 |
| CMD contrast | 2.351 | 0.912 | 6.652 | 10.497 (1.758–62.656) | 0.010 |
| Cortical signal heterogeneity | 0.093 | 0.041 | 5.255 | 1.098 (1.014–1.189) | 0.022 |
| RBF values | −0.038 | 0.032 | 9.605 | 0.963 (0.940–0.986) | 0.002 |
CI, confidence interval; CKD, chronic kidney disease; CMD, corticomedullary differentiation; OR, odds ratio; RBF, renal blood flow; SE, standard error.
Correlations between eGFR and the independent predictors
A positive correlation was observed between renal length (r=0.191, P=0.054), CMD contrast (r=0.484, P<0.001), and RBF values (r=0.519, P<0.001), whereas eGFR was significantly negatively correlated with cortical signal heterogeneity (r=−0.339, P<0.001).
Diagnostic efficacy for early-stage CKD
The results of the ROC curve analysis are presented in Table 6. The AUC values for renal length, CMD contrast, cortical signal heterogeneity, RBF values, and the combined model were 0.629 (95% CI: 0.528–0.722), 0.707 (95% CI: 0.609–0.793), 0.771 (95% CI: 0.678–0.848), 0.802 (95% CI: 0.712–0.874), and 0.916 (95% CI: 0.844–0.961), respectively. Compared to renal length (Z=4.828, P<0.001), cortical signal heterogeneity (Z=3.270, P=0.001), CMD contrast (Z=5.202, P<0.001), and RBF values alone (Z=2.778, P=0.005), the combined RBF-morphology model demonstrated the highest diagnostic efficacy, with an AUC of 0.916 (95% CI: 0.844–0.961), a specificity of 87.88%, and a sensitivity of 81.08% at the optimal cut-off value of 44.0% (Figure 3).
Table 6
| Characteristics | AUC (95% CI) | P value | Sensitivity (%) | Specificity (%) | Cut-off value | Youden’s index |
|---|---|---|---|---|---|---|
| Renal length | 0.629 (0.528–0.722) | 0.024 | 32.43 | 89.39 | 88.15 | 0.218 |
| CMD contrast | 0.707 (0.609–0.793) | <0.001 | 45.95 | 95.45 | 2.500 | 0.414 |
| Cortical signal heterogeneity | 0.771 (0.678–0.848) | <0.001 | 62.16 | 83.33 | 49.972 | 0.455 |
| RBF values | 0.802 (0.712–0.874) | <0.001 | 75.68 | 78.79 | 201.75 | 0.545 |
| RBF-morphology model | 0.916 (0.844–0.961) | <0.001 | 81.08 | 87.88 | 0.440 | 0.690 |
AUC, area under the curve; CI, confidence interval; CKD, chronic kidney disease; CMD, corticomedullary differentiation; RBF, renal blood flow.
Discussion
This study utilized 5 T MRI to achieve imaging of subtle renal tissue structures and combined it with ASL perfusion imaging to differentiate early- from mid-stage CKD from both morphological and functional perspectives. The results demonstrate that renal length, CMD contrast, cortical signal heterogeneity, and RBF values are independent factors for differentiating CKD stages. With declining eGFR, these parameters progressively decrease. The combined RBF-morphology model achieved an AUC of 0.916 for differentiating early- from mid-stage CKD, significantly outperforming single-characteristic assessments (P<0.05), demonstrating superior diagnostic value. These findings suggest that perfusion and morphological biomarkers derived from 5 T MRI may serve as complementary imaging indicators for assessing renal structural and functional alterations in CKD, and could help to propel the clinical translation of 5 T MRI.
Morphological assessment
In this study, renal length decreased significantly in mid-stage CKD, with concomitant reductions in renal width and cortical thickness. During CKD progression, pathological changes including glomerulosclerosis, tubular atrophy, and interstitial fibrosis lead to progressive renal parenchymal loss, manifesting as cortical thinning and medullary fibrosis, ultimately causing a 3D reduction (9,10). Renal length, as a key morphological indicator, is commonly used to characterize kidney size (26). Beland et al. demonstrated a significant correlation between renal length and eGFR (P=0.003), suggesting its potential as a renal functional biomarker (27). At the early-stage, despite the initial structural damage, compensatory mechanisms such as glomerular hyperfiltration and tubular solute load regulation in the residual nephrons help to maintain relatively stable renal function. This compensatory effect causes morphological changes (renal length reduction) to precede functional decline (eGFR reduction).
Our study revealed that cortical signal heterogeneity increased significantly, with significant degradation of CMD contrast. Mid-stage CKD patients frequently exhibit characteristic “mosaic-like” signal heterogeneity on T1WI, attributable to spatial heterogeneity of renal injury. In early-stage diseases, fibrotic foci demonstrate prolonged T1 relaxation times (hypointensity); meanwhile, adjacent edematous areas also show prolonged T1 relaxation times (hypointensity) (19). However, normal renal cortex exhibits relatively shorter T1 values, appearing hyperintense compared to fibrotic or edematous tissue. This spatial intermingling between pathology and normal tissue results in T1WI signal heterogeneity. Inflammation-driven edema in renal impairment increases cortical free water and lengthens T1, attenuating cortical signal and CMD (15,19,28), while fibrosis further prolongs T1 and deepens CMD loss by restricting diffusion (19,29-31). Marotti et al. (29) and Kettritz et al. (12) emphasized that a reduction or complete loss of CMD contrast is a sensitive but non-specific indicator. In contrast to these findings, our study demonstrated a relatively high specificity (95.45%) for CMD contrast. This discrepancy may be attributed to the subjective assessment of CMD, introducing potential observer variability, and selection bias within the study population.
Functional assessment
ASL is a non-invasive imaging technique that provides a reliable method for quantifying RBF (32-34). Under normal physiological conditions, RBF accounts for 20–25% of cardiac output and exhibits regional distribution differences. Over 90% of blood flow is distributed to the renal cortex with its high metabolic demand, whereas only a small portion reaches the medulla (33,35,36). This anatomical distribution results in the cortex exhibiting significantly high-perfusion characteristics. During the progression of CKD, renal hemodynamic disturbances lead to progressive perfusion deficits in the peritubular capillary network and microvascular loss, ultimately resulting in a significant reduction in cortical blood flow (37,38). This study found that RBF is an independent predictor for diagnosing early- and mid-stage CKD [odds ratio (OR) =0.963, 95% CI: 0.940–0.986] and shows a significantly positive correlation with eGFR (r=0.519, P<0.001). These findings are consistent with those reported by Mao et al. (21), which were confirmed through ROC curve analysis that RBF values exhibited good discriminative ability for early renal injury (AUC =0.827, 95% CI: 0.687–0.922) and showed a moderate correlation with eGFR (r=0.586, P<0.05), suggesting a close association between the degree of renal ischemia and renal functional decompensation. Compared to traditional eGFR assessments based on Scr, ASL directly quantifies microcirculatory perfusion in the renal cortex. This approach not only overcomes the limitations of creatinine metabolism, which is influenced by muscle mass and diet, but also provides a direct visualization of the dynamic progression from compensatory hyperperfusion to decompensated hypoperfusion in nephrons (25,39).
The combined RBF-morphology model assessment
Morphological changes and RBF values in CKD can reflect the pathological progression of the disease and are closely related to eGFR. Currently, the clinical assessment of renal morphology primarily relies on ultrasound, which is susceptible to interference from intestinal gas and surrounding tissues and is strongly operator-dependent, thus compromising the reproducibility of morphological evaluations (23,26). In contrast, MRI with its multi-parametric and multi-planar imaging capabilities, enables 3D visualization of renal anatomy. Notably, 5 T MRI provides significantly superior SNR and contrast-to-noise ratio (CNR) compared to 3 T, allowing clearer delineation of fine anatomical details and more sensitive detection of subtle lesions (13,40). Consistent with this, Zheng et al. (13) demonstrated that 5 T MRI yields morphological and functional renal image quality comparable or superior to that of 3 T, with the enhanced SNR and CNR improving intra-lesional structural visualization and pathological diagnostic accuracy. Furthermore, the increase in SNR and the longer T1 relaxation time of blood at higher field strengths may improve ASL imaging, potentially allowing for a shorter scan time, higher spatial resolution, and increased sensitivity to low perfusion levels (41-43). Supporting this, Tong et al. (44) showed that, compared to 3 T, 5 T yielded a 20–50% increase in SNR and a 60–70% increase in CNR, with excellent repeatability and reproducibility of RBF values (ICC =0.855–0.973). Capitalizing on these advantages, our study employed 5 T MRI for renal morphological and ASL perfusion imaging. ROC curve analysis confirmed that the combined RBF-morphology model achieved an AUC of 0.916 (95% CI: 0.844–0.961) for distinguishing between early- and mid-stage CKD, significantly outperforming any single-characteristic assessment (P<0.05). These findings indicate that the RBF-morphology combination has good diagnostic value for differentiating early-stage from mid-stage CKD, and may serve as a complementary imaging indicator for assessing renal structural and functional alterations in CKD.
Despite providing valuable insights, this study has several limitations that warrant consideration. First, ROIs were manually delineated, which may have introduced subjective bias and affected the accuracy of the results. Future studies could employ deep learning-based automatic ROI segmentation algorithms to improve measurement reproducibility. Second, the lack of patient oxygenation level assessment limits deeper pathophysiological analysis of the hemodynamic disturbances. Third, the single-center design and small sample size may restrict the clinical applicability of the conclusions. Future research should expand the sample size and conduct multi-center studies to enhance the reliability of the findings. Fourth, CKD staging in this study was based on eGFR rather than directly measured glomerular filtration rate (GFR), which may introduce classification bias due to inherent inaccuracies of estimation equations. Future studies incorporating measured GFR could help validate the imaging biomarkers against a more definitive reference standard.
Conclusions
The combination of ASL and morphological characteristics achieves significantly better diagnostic performance between early- and mid-stage CKD than any single imaging characteristic, suggesting that 5 T MRI derived perfusion and morphological biomarkers may serve as complementary imaging indicators for assessing renal structural and functional alterations in CKD.
Acknowledgments
We would especially like to thank Dr. Zhichao Feng for his valuable suggestions on the revision of this manuscript.
Footnote
Reporting Checklist: The authors have completed the STARD reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0296/rc
Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0296/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-2026-1-0296/coif). S.X. and L.Z. are currently employed by United Imaging Healthcare. 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. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. This study was approved by the Ethics Committee of Chongqing University Three Gorges Hospital [No. 2024-KS-(149)]. Written informed consent was provided by all participants prior to enrollment.
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
- Global, regional, and national burden of chronic kidney disease, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet 2020;395:709-33.
- Francis A, Harhay MN, Ong ACM, Tummalapalli SL, Ortiz A, Fogo AB, et al. Chronic kidney disease and the global public health agenda: an international consensus. Nat Rev Nephrol 2024;20:473-85. [Crossref] [PubMed]
- Song J, Pan T, Xu Z, Yi W, Pan R, Cheng J, Hu W, Su H. A systematic analysis of chronic kidney disease burden attributable to lead exposure based on the global burden of disease study 2019. Sci Total Environ 2024;908:168189. [Crossref] [PubMed]
- Levin A, Okpechi IG, Caskey FJ, Yang CW, Tonelli M, Jha V. Perspectives on early detection of chronic kidney disease: the facts, the questions, and a proposed framework for 2023 and beyond. Kidney Int 2023;103:1004-8. [Crossref] [PubMed]
- Shlipak MG, Tummalapalli SL, Boulware LE, Grams ME, Ix JH, Jha V, Kengne AP, Madero M, Mihaylova B, Tangri N, Cheung M, Jadoul M, Winkelmayer WC, Zoungas S. The case for early identification and intervention of chronic kidney disease: conclusions from a Kidney Disease: Improving Global Outcomes (KDIGO) Controversies Conference. Kidney Int 2021;99:34-47. [Crossref] [PubMed]
- KDIGO 2024 Clinical Practice Guideline for the Evaluation and Management of Chronic Kidney Disease. Kidney Int 2024;105:S117-314.
- Nie S, Zhou S, Cabrera C, Chen S, Jia M, Zhang S, Su L, Gao Q, Tangri N, Hou FF. Characteristics of patients with undiagnosed stage 3 chronic kidney disease: results from an observational study (REVEAL-CKD) in China. Lancet Reg Health West Pac 2025;54:101275. [Crossref] [PubMed]
- Luciano RL, Moeckel GW. Update on the Native Kidney Biopsy: Core Curriculum 2019. Am J Kidney Dis 2019;73:404-15. [Crossref] [PubMed]
- Francis ST, Selby NM, Taal MW. Magnetic Resonance Imaging to Evaluate Kidney Structure, Function, and Pathology: Moving Toward Clinical Application. Am J Kidney Dis 2023;82:491-504. [Crossref] [PubMed]
- Mao W, Ding X, Ding Y, Cao B, Fu C, Kuehn B, Benkert T, Grimm R, Nickel D, Zhou J, Zeng M. Evaluation of interstitial fibrosis in chronic kidney disease by multiparametric functional MRI and histopathologic analysis. Eur Radiol 2023;33:4138-47. [Crossref] [PubMed]
- Wei C, Liu R, Li X, Zhou C, Ma Q, Xu Y, Zhu Y, Shen J, Zeng Y, Song K, Jiang Z. Evaluation of multiparametric MRI and clinical indicators for renal fibrosis in chronic kidney disease. Quant Imaging Med Surg 2025;15:6200-16. [Crossref] [PubMed]
- Kettritz U, Semelka RC, Brown ED, Sharp TJ, Lawing WL, Colindres RE. MR findings in diffuse renal parenchymal disease. J Magn Reson Imaging 1996;6:136-44. [Crossref] [PubMed]
- Zheng L, Yang C, Sheng R, Dai Y, Zeng M. Renal imaging at 5 T versus 3 T: a comparison study. Insights Imaging 2022;13:155. [Crossref] [PubMed]
- Zhang Y, Yang C, Liang L, Shi Z, Zhu S, Chen C, Dai Y, Zeng M. Preliminary Experience of 5.0 T Higher Field Abdominal Diffusion-Weighted MRI: Agreement of Apparent Diffusion Coefficient With 3.0 T Imaging. J Magn Reson Imaging 2022;56:1009-17. [Crossref] [PubMed]
- Noda Y, Ito K, Kanki A, Tamada T, Yamamoto A, Kazuya Y, Higaki A. Measurement of renal cortical thickness using noncontrast-enhanced steady-state free precession MRI with spatially selective inversion recovery pulse: Association with renal function. J Magn Reson Imaging 2015;41:1615-21. [Crossref] [PubMed]
- Hueper K, Gutberlet M, Rong S, Hartung D, Mengel M, Lu X, Haller H, Wacker F, Meier M, Gueler F. Acute kidney injury: arterial spin labeling to monitor renal perfusion impairment in mice-comparison with histopathologic results and renal function. Radiology 2014;270:117-24. [Crossref] [PubMed]
- Wu WC, Su MY, Chang CC, Tseng WY, Liu KL. Renal perfusion 3-T MR imaging: a comparative study of arterial spin labeling and dynamic contrast-enhanced techniques. Radiology 2011;261:845-53. [Crossref] [PubMed]
- Khatir DS, Pedersen M, Jespersen B, Buus NH. Reproducibility of MRI renal artery blood flow and BOLD measurements in patients with chronic kidney disease and healthy controls. J Magn Reson Imaging 2014;40:1091-8. [Crossref] [PubMed]
- Buchanan CE, Mahmoud H, Cox EF, McCulloch T, Prestwich BL, Taal MW, Selby NM, Francis ST. Quantitative assessment of renal structural and functional changes in chronic kidney disease using multi-parametric magnetic resonance imaging. Nephrol Dial Transplant 2020;35:955-64. [Crossref] [PubMed]
- Rossi C, Artunc F, Martirosian P, Schlemmer HP, Schick F, Boss A. Histogram analysis of renal arterial spin labeling perfusion data reveals differences between volunteers and patients with mild chronic kidney disease. Invest Radiol 2012;47:490-6. [Crossref] [PubMed]
- Mao W, Ding Y, Ding X, Fu C, Cao B, Kuehn B, Benkert T, Grimm R, Zhou J, Zeng M. Capability of arterial spin labeling and intravoxel incoherent motion diffusion-weighted imaging to detect early kidney injury in chronic kidney disease. Eur Radiol 2023;33:3286-94. [Crossref] [PubMed]
- Thakur V, Watkins T, McCarthy K, Beidl T, Underwood N, Barnes K, Cook ME. Is kidney length a good predictor of kidney volume? Am J Med Sci 1997;313:85-9. [Crossref] [PubMed]
- Bakker J, Olree M, Kaatee R, de Lange EE, Beek FJ. In vitro measurement of kidney size: comparison of ultrasonography and MRI. Ultrasound Med Biol 1998;24:683-8. [Crossref] [PubMed]
- Lee VS, Kaur M, Bokacheva L, Chen Q, Rusinek H, Thakur R, Moses D, Nazzaro C, Kramer EL. What causes diminished corticomedullary differentiation in renal insufficiency? J Magn Reson Imaging 2007;25:790-5. [Crossref] [PubMed]
- Nery F, Buchanan CE, Harteveld AA, Odudu A, Bane O, Cox EF, et al. Consensus-based technical recommendations for clinical translation of renal ASL MRI. MAGMA 2020;33:141-61. [Crossref] [PubMed]
- Cheong B, Muthupillai R, Rubin MF, Flamm SD. Normal values for renal length and volume as measured by magnetic resonance imaging. Clin J Am Soc Nephrol 2007;2:38-45. [Crossref] [PubMed]
- Beland MD, Walle NL, Machan JT, Cronan JJ. Renal cortical thickness measured at ultrasound: is it better than renal length as an indicator of renal function in chronic kidney disease? AJR Am J Roentgenol 2010;195:W146-9. [Crossref] [PubMed]
- Suo S, Lin N, Wang H, Zhang L, Wang R, Zhang S, Hua J, Xu J. Intravoxel incoherent motion diffusion-weighted MR imaging of breast cancer at 3.0 tesla: Comparison of different curve-fitting methods. J Magn Reson Imaging 2015;42:362-70. [Crossref] [PubMed]
- Marotti M, Hricak H, Terrier F, McAninch JW, Thuroff JW. MR in renal disease: importance of cortical-medullary distinction. Magn Reson Med 1987;5:160-72. [Crossref] [PubMed]
- Friedli I, Crowe LA, Berchtold L, Moll S, Hadaya K, de Perrot T, Vesin C, Martin PY, de Seigneux S, Vallée JP. New Magnetic Resonance Imaging Index for Renal Fibrosis Assessment: A Comparison between Diffusion-Weighted Imaging and T1 Mapping with Histological Validation. Sci Rep 2016;6:30088. [Crossref] [PubMed]
- Friedli I, Crowe LA, de Perrot T, Berchtold L, Martin PY, de Seigneux S, Vallée JP. Comparison of readout-segmented and conventional single-shot for echo-planar diffusion-weighted imaging in the assessment of kidney interstitial fibrosis. J Magn Reson Imaging 2017;46:1631-40. [Crossref] [PubMed]
- Martirosian P, Boss A, Schraml C, Schwenzer NF, Graf H, Claussen CD, Schick F. Magnetic resonance perfusion imaging without contrast media. Eur J Nucl Med Mol Imaging 2010;37:S52-64. [Crossref] [PubMed]
- Chowdhury AH, Cox EF, Francis ST, Lobo DN. A randomized, controlled, double-blind crossover study on the effects of 1-L infusions of 6% hydroxyethyl starch suspended in 0.9% saline (voluven) and a balanced solution (Plasma Volume Redibag) on blood volume, renal blood flow velocity, and renal cortical tissue perfusion in healthy volunteers. Ann Surg 2014;259:881-7. [Crossref] [PubMed]
- Wan J, Jin M, Li J, Ma J, Que C, Jiang B, Tian Y, Hu L, Yu Y, Hu C, Wang J, Zhu M. Magnetic resonance diffusion tensor imaging is superior to arterial spin labeling in detecting renal allograft fibrosis. Quant Imaging Med Surg 2025;15:3211-21. [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]
- Pallone TL, Edwards A, Mattson DL. Renal medullary circulation. Compr Physiol 2012;2:97-140. [Crossref] [PubMed]
- Li LP, Tan H, Thacker JM, Li W, Zhou Y, Kohn O, Sprague SM, Prasad PV. Evaluation of Renal Blood Flow in Chronic Kidney Disease Using Arterial Spin Labeling Perfusion Magnetic Resonance Imaging. Kidney Int Rep 2017;2:36-43. [Crossref] [PubMed]
- Singh P, Ricksten SE, Bragadottir G, Redfors B, Nordquist L. Renal oxygenation and haemodynamics in acute kidney injury and chronic kidney disease. Clin Exp Pharmacol Physiol 2013;40:138-47. [Crossref] [PubMed]
- Zhang X, Ye C, Lu F, Yang J, Xu Y, Wang C. Evaluation of renal oxygenation and perfusion in patients with chronic kidney disease: a preliminary prospective study based on functional magnetic resonance. Ren Fail 2024;46:2428337. [Crossref] [PubMed]
- Lu H, Miao X, Wang D, Zheng X, Zhang S, Wang R, Li G, Qian X, Chen Y, Hu C, Jin H, Zeng M. Feasibility and Clinical Application of 5-T Noncontrast Dixon Whole-Heart Coronary MR Angiography: A Prospective Study. Radiology 2024;313:e240389. [Crossref] [PubMed]
- Wang K, Ma SJ, Shao X, Zhao C, Shou Q, Yan L, Wang DJJ. Optimization of pseudo-continuous arterial spin labeling at 7T with parallel transmission B1 shimming. Magn Reson Med 2022;87:249-62. [Crossref] [PubMed]
- Zhao C, Shao X, Shou Q, Ma SJ, Gokyar S, Graf C, Stollberger R, Wang DJ. Whole-Cerebrum distortion-free three-dimensional pseudo-continuous arterial spin labeling at 7T. Neuroimage 2023;277:120251. [Crossref] [PubMed]
- Hernandez-Garcia L, Aramendía-Vidaurreta V, Bolar DS, Dai W, Fernández-Seara MA, Guo J, et al. Recent Technical Developments in ASL: A Review of the State of the Art. Magn Reson Med 2022;88:2021-42. [Crossref] [PubMed]
- Tong X, Chen C, Ning Z, Shen R, Yang N, Zhang S, Chen S, He ZX, Zhao X. Saturated Multi-Delay Renal Arterial Spin Labeling Technique at 5 T MR: A Comparative Study With 3 T. J Magn Reson Imaging 2026;63:130-9. [Crossref] [PubMed]
(English Language Editor: J. Jones)

