Non-invasive risk stratification of immunoglobulin A nephropathy and crescent formation using contrast-enhanced ultrasound combined with serological markers
Introduction
Immunoglobulin A nephropathy (IgAN) is the most prevalent primary glomerulonephritis globally and a leading cause of end-stage renal disease (1,2). Its clinical manifestations are highly heterogeneous, ranging from isolated hematuria with preserved renal function to progressive proteinuria, renal insufficiency, and even end-stage renal disease. This heterogeneity highlights the importance of early identification of pathological activity and disease severity (3).
Among the pathological features of IgAN, crescent formation has special clinical significance. Crescents reflect active glomerular injury and are usually associated with a more aggressive disease course, higher inflammatory burden, and poorer renal prognosis. Therefore, the ability to identify patients at increased risk of crescent formation is crucial for risk stratification, treatment selection, and follow-up intensity.
Although renal biopsy remains the diagnostic gold standard, its invasiveness and sampling variability restrict repeated use in routine practice. Non-invasive tools for identifying patients with a high probability of IgAN and for stratifying crescent risk may therefore support earlier risk assessment and clinical decision-making.
Contrast-enhanced ultrasound (CEUS) provides real-time quantitative assessment of renal microcirculation without ionizing radiation or nephrotoxic contrast agents (4,5). Compared with traditional gray-scale ultrasound, CEUS can provide dynamic functional information about renal perfusion. Parameters derived from time-intensity curves (TICs), such as peak enhancement (PE), wash-in rate (WiR), and wash-in perfusion index (WiP), have shown correlations with histopathological severity (6-10). Integrating these with serological markers like estimated glomerular filtration rate (eGFR) and urine albumin-to-creatinine ratio (UACR) could enhance diagnostic and prognostic accuracy.
In this study, we aimed to (I) examine the associations between CEUS parameters, serological biomarkers, and Oxford MEST-C (mesangial hypercellularity, endocapillary hypercellularity, segmental sclerosis, tubular atrophy/interstitial fibrosis, crescent formation) scores, with a focus on crescents; (II) develop predictive models for IgAN diagnosis and crescent formation; and (III) determine whether combining CEUS with clinical biomarkers improves diagnostic performance compared with clinical-only or CEUS-only models. We present this article in accordance with the TRIPOD reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0025/rc).
Methods
Patient selection and ethics approval
This retrospective study included patients who visited the Department of Nephrology, Huashan Hospital, Fudan University, for persistent proteinuria, hematuria, or unexplained renal impairment between February 28, 2023, and May 9, 2024, and who were considered by clinicians to have clear indications for renal biopsy.
Inclusion criteria were as follows: (I) age ≥18 years; (II) chronic kidney disease (CKD) diagnosed according to Kidney Disease: Improving Global Outcomes (KDIGO) guidelines, defined as evidence of kidney injury (including proteinuria, hematuria, or imaging abnormalities) and/or eGFR <60 mL/min/1.73 m2 lasting ≥3 months; (III) CKD stage 1–5 with a biopsy-confirmed pathological diagnosis; and (IV) CEUS examination completed within 7 days before renal biopsy, with stable clinical status and no acute change in renal function during the examination period.
Exclusion criteria were as follows: (I) acute kidney injury at enrollment; (II) marked fluctuation in renal function within 1 month before examination (eGFR change >30%); (III) an interval of more than 7 days between CEUS examination and renal biopsy; and (IV) poor CEUS image quality that precluded quantitative analysis.
A total of 184 eligible patients were included. The sample size was based on the number of eligible patients available during the study period. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. This retrospective study was approved by the Ethics Committee of Huashan Hospital, Fudan University (Approval No. KY2024-730), and all patients provided informed consent. Among the 184 patients who completed renal biopsy, 94 were diagnosed with IgAN and 90 were diagnosed with non-IgAN kidney diseases.
Among the 94 patients with IgAN, crescent status was defined according to the Oxford classification C score. Patients with C1 or C2 lesions were assigned to the crescent subgroup, whereas patients with C0 lesions were assigned to the non-crescent subgroup. Each subgroup included 47 patients.
Clinical and laboratory data collection
Demographic information (age and gender) and laboratory values, including eGFR calculated using the CKD-EPI equation, serum creatinine, urea, albumin, and UACR, were recorded within 7 days before biopsy.
Contrast-enhanced ultrasound examination
Examinations were performed using a Siemens ACUSON Sequoia (Siemens Healthineers) color Doppler ultrasound system with a 1–6 MHz convex array probe for CEUS at low mechanical index. After conventional B-mode imaging, 1.0 mL of SonoVue (Bracco, Milan, Italy) was injected intravenously via the antecubital vein, followed by a 5 mL normal saline flush. Valid 90-second CEUS videos were recorded and saved in Digital Imaging and Communications in Medicine (DICOM) format.
CEUS videos in DICOM format were analyzed using VueBox 7.0 software with cortex-to-medulla normalization. CEUS analysis was performed independently without knowledge of pathological results. The regions of interest (ROI) consisted of one elliptical area (approximately 4 mm × 3 mm) in the renal cortex and one in the renal medulla, placed in the mid renal cortex (red) and mid renal medulla (blue) closest to the probe, respectively (Figure 1A). The software generated TICs (Figure 1B), which reflected dynamic changes in renal cortical and medullary microcirculatory perfusion. ROI placement was performed by an experienced radiologist and reviewed by a second investigator.
CEUS TIC analysis generated a series of renal perfusion parameters, including rise time (RT), time to peak (TTP), fall time (FT), mean transit time (mTT), PE, wash-in area under the curve (WiAUC), wash-out AUC (WoAUC), WiR, WiP, wash-out rate (WoR), and wash-in and wash-out AUC (WiWoAUC). The definitions of the key parameters are shown in Figure 1B: PE was defined as the maximum peak intensity above baseline (dB), RT as the time required for the signal to increase from 10% to 100% of PE(s), FT as the time required for the signal to decrease from PE to baseline(s), TTP as the time from contrast agent arrival to PE(s), WiR as the slope of the ascending segment (dB/s), WiP as PE/RT (dB/s), and WoR as the slope of the descending segment (dB/s). mTT, WiAUC, WoAUC, and WiWoAUC were also recorded. Only curves with a goodness of fit ≥85% were included.
The cortex-to-medulla ratio was used as a normalization index to quantify renal microcirculation. This approach was adopted for two reasons. First, CEUS parameters are affected by technical factors, including contrast-agent dose, injection speed, machine gain, imaging depth, and probe angle, as well as individual factors such as blood pressure, cardiac function, cardiac output, and respiratory motion. Second, renal blood flow is distributed unevenly, with the cortex receiving approximately 90% of renal blood flow and the medulla approximately 10%. Decreased cortical microcirculation is an early and important feature of CKD, whereas medullary perfusion impairment tends to occur later and may be less pronounced. Cortex-to-medulla normalization may therefore reduce technical and physiological variability and improve data stability and comparability. Each subject was measured twice using the same method, and the average value was used for statistical analysis.
Renal biopsy and histopathological evaluation
All patients underwent ultrasound-guided percutaneous renal biopsy under local anesthesia, and renal tissue samples were obtained using a 16 G biopsy needle. The samples were fixed in 10% neutral buffered formalin, embedded in paraffin, and examined by light microscopy and immunofluorescence. Pathological assessment was performed by experienced pathologists who were blinded to CEUS findings. All patients with IgAN were newly diagnosed, untreated cases without recurrence.
The diagnosis of IgAN was based on renal biopsy findings, characterized by predominant IgA deposition in the glomerular mesangium on immunofluorescence together with compatible light microscopic findings. According to renal biopsy pathology, CKD patients were divided into the IgAN group and the non-IgAN group. Pathological lesions in IgAN were evaluated using the Oxford MEST-C scoring system, including mesangial hypercellularity (M), endocapillary hypercellularity (E), segmental sclerosis (S), tubular atrophy/interstitial fibrosis (T), and crescent formation (C).
According to the MEST-C scoring system in the Oxford classification, crescents in IgAN patients were graded: C0: no crescents; C1: crescents involving <25% of glomeruli; C2: crescents involving ≥25% of glomeruli. In this study, C1 + C2 grades were defined as the crescent formation group (1,2).
Statistical analysis
SPSS 25.0 (IBM, USA) and R 4.1.0 (R Foundation, Vienna) were used for statistical analysis. All statistical tests were two-tailed, and P<0.05 was considered statistically significant. Normality was assessed using the Shapiro-Wilk test. Normally distributed continuous variables were expressed as mean ± standard deviation and compared using independent-samples t-tests. Non-normally distributed variables were expressed as median (interquartile range) and compared using the Mann-Whitney U test. Categorical variables were expressed as n (%) and compared using the Chi-squared test or Fisher’s exact test, as appropriate. Spearman correlation analysis was used to evaluate associations with MEST-C scores.
Univariate logistic regression was used to screen potential predictors. Variables with P<0.05 in univariate analysis were entered into multivariable logistic regression. Before multivariable analysis, multicollinearity was assessed using the variance inflation factor (VIF), with VIF <5 indicating no significant multicollinearity. Final predictors were retained according to their multivariable associations with each endpoint, which explains why different CEUS parameters entered different models. Independent predictors of IgAN diagnosis and crescent formation were identified, and odds ratios (ORs) with 95% confidence intervals (CIs) were calculated. Diagnostic models were constructed based on multivariable logistic regression results and visualized as nomograms.
Model performance was evaluated using receiver operating characteristic (ROC) curves, and AUCs, 95% CIs, sensitivity, and specificity were calculated. Positive predictive value (PPV) and negative predictive value (NPV) were derived from sensitivity, specificity, and the observed prevalence of each endpoint in the study cohort. Calibration curves were generated to evaluate agreement between predicted and observed probabilities. Decision curve analysis (DCA) was used to assess clinical net benefit. Optimal cut-off values were determined using the Youden index. DeLong tests were used for pairwise comparisons between ROC curves. Bootstrap resampling was used for internal validation.
Results
Baseline characteristics and CEUS metrics
Among the 184 patients, 94 were diagnosed with IgAN. Compared with non-IgAN patients, IgAN patients were significantly younger (40.3±12.8 vs. 49.4±15.1 years; P<0.001) and had lower eGFR (78.0±33.4 vs. 90.3±47.5 mL/min/1.73 m2; P=0.045), lower UACR, and higher albumin levels (P=0.021; Table 1). Among non-IgAN cases, membranous nephropathy (15.2%), lupus nephritis (7.1%), focal segmental glomerulosclerosis (6.0%), Henoch-Schonlein purpura nephritis (3.8%), diabetic nephropathy (3.3%), and other pathological types (13.6%) were confirmed. Crescent formation varied across diagnoses, with the highest prevalence in Henoch-Schonlein purpura nephritis (71.4%) and lupus nephritis (23.1%). No crescents were observed in membranous nephropathy, focal segmental glomerulosclerosis, or diabetic nephropathy (Table 2). CEUS showed that PE, WiR, and WiP were significantly lower in the IgAN group (all P<0.05; Table 3). The inclusion of several non-IgAN pathological entities allowed further evaluation of the crescent prediction model in the entire cohort.
Table 1
| Variables | Non-IgAN (n=90) | IgAN (n=94) | P |
|---|---|---|---|
| Gender | 0.369 | ||
| Female | 40 (44.44) | 48 (51.06) | |
| Male | 50 (55.56) | 46 (48.94) | |
| Age (years) | 49.42±15.05 | 40.29±12.75 | <0.001 |
| eGFR (mL/min/1.73 m2) | 90.31±47.51 | 78.04±33.44 | 0.045 |
| Urea (mmol/L) | 8.46±5.31 | 7.51±3.30 | 0.149 |
| SCr (μmol/L) | 86.5 (68.0, 115.0) | 92.0 (74.0, 128.0) | 0.857 |
| ALB (g/L) | 35.26±7.94 | 37.64±5.73 | 0.021 |
| UACR (mg/g) | 1,850 (890, 4200) | 1,120 (640, 2350) | 0.005 |
Data are expressed as mean ± standard deviation, median (interquartile range), or n (%) unless otherwise stated. ALB, albumin; eGFR, estimated glomerular filtration rate; IgAN, immunoglobulin A nephropathy; SCr, serum creatinine; UACR, urine albumin-to-creatinine ratio.
Table 2
| Pathological diagnosis | Total | Crescent present | Crescent absent |
|---|---|---|---|
| IgAN | 94 (51.1) | 47 (50.0) | 47 (50.0) |
| Membranous nephropathy | 28 (15.2) | 0 (0.0) | 28 (100.0) |
| Lupus nephritis | 13 (7.1) | 3 (23.1) | 10 (76.9) |
| Focal segmental glomerulosclerosis | 11 (6.0) | 0 (0.0) | 11 (100.0) |
| Henoch-Schönlein purpura nephritis | 7 (3.8) | 5 (71.4) | 2 (28.6) |
| Diabetic nephropathy | 6 (3.3) | 0 (0.0) | 6 (100.0) |
| Other pathological types | 25 (13.6) | 1 (4.0) | 24 (96.0) |
Data are presented as n (%). IgAN, immunoglobulin A nephropathy.
Table 3
| Variables | Non-IgAN (n=90) | IgAN (n=94) | P |
|---|---|---|---|
| PE | 2.09±0.75 | 1.81±0.62 | 0.007 |
| RT | 0.38 (0.21, 0.65) | 0.45 (0.28, 0.72) | 0.080 |
| mTT | 0.51 (0.35, 0.78) | 0.54 (0.38, 0.75) | 0.599 |
| TTP | 0.68 (0.45, 0.95) | 0.72 (0.48, 0.98) | 0.352 |
| FT | 0.62 (0.48, 0.85) | 0.56 (0.41, 0.78) | 0.039 |
| WiR | 2.45 (1.85, 3.20) | 2.05 (1.45, 2.75) | 0.002 |
| WiP | 2.32±0.83 | 2.07±0.82 | 0.046 |
| WoR | 2.81±1.33 | 2.79±1.49 | 0.920 |
| WiAUC | 1.42 (1.15, 1.75) | 1.40 (1.10, 1.78) | 0.856 |
| WoAUC | 0.52 (0.35, 0.78) | 0.54 (0.38, 0.75) | 0.956 |
| WiWoAUC | 0.72 (0.55, 0.98) | 0.75 (0.58, 1.02) | 0.525 |
Data are expressed as mean ± standard deviation or median (interquartile range). All CEUS parameters are presented as dimensionless cortex-to-medulla ratios. CEUS, contrast-enhanced ultrasound; FT, fall time; IgAN, immunoglobulin A nephropathy; mTT, mean transit time; PE, peak enhancement; RT, rise time; TTP, time to peak; WiAUC, wash-in area under the curve; WiP, wash-in perfusion index; WiR, wash-in rate; WiWoAUC, wash-in and wash-out area under the curve; WoAUC, wash-out area under the curve; WoR, wash-out rate.
Correlation analysis
Spearman correlation analysis demonstrated partial significant associations between Oxford classification scores [mesangial hypercellularity (M), endocapillary hypercellularity (E), segmental sclerosis (S), tubular atrophy/interstitial fibrosis (T), and crescent formation (C)] and both serological markers (eGFR and serum creatinine) and CEUS parameters (PE, WiR, and WiP) (P<0.05), indicating that CEUS metrics reflect histopathological severity (Table 4).
Table 4
| Variables | Age | eGFR | Urea | SCr | ALB | UACR | PE | RT | mTT | TTP | FT | WiR | WiP | WoR | WiAUC | WoAUC | WiWoAUC |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| M (r) | –0.166 | –0.112 | 0.008 | 0.065 | –0.163 | 0.247 | –0.058 | 0.054 | 0.139 | 0.069 | 0.081 | –0.071 | –0.127 | –0.063 | –0.105 | 0.063 | 0.034 |
| M (P) | 0.111 | 0.281 | 0.939 | 0.532 | 0.117 | 0.017 | 0.576 | 0.603 | 0.183 | 0.507 | 0.436 | 0.494 | 0.224 | 0.546 | 0.312 | 0.549 | 0.748 |
| E (r) | –0.095 | 0.171 | –0.052 | –0.134 | –0.178 | 0.336 | –0.002 | –0.035 | –0.078 | –0.106 | –0.122 | 0.091 | –0.007 | 0.046 | –0.026 | 0.004 | 0.077 |
| E (P) | 0.364 | 0.1 | 0.616 | 0.198 | 0.086 | 0.001 | 0.988 | 0.736 | 0.458 | 0.311 | 0.241 | 0.383 | 0.943 | 0.662 | 0.807 | 0.972 | 0.463 |
| S (r) | –0.251 | 0.011 | –0.004 | –0.068 | –0.019 | 0.006 | –0.004 | 0.111 | –0.008 | 0.055 | 0.021 | –0.015 | –0.085 | –0.026 | –0.051 | –0.148 | –0.013 |
| S (P) | 0.015 | 0.918 | 0.968 | 0.515 | 0.853 | 0.954 | 0.967 | 0.288 | 0.941 | 0.601 | 0.839 | 0.883 | 0.414 | 0.802 | 0.624 | 0.153 | 0.902 |
| T (r) | –0.043 | –0.325 | 0.164 | 0.283 | –0.121 | 0.227 | –0.267 | 0.107 | 0.133 | 0.182 | –0.004 | –0.285 | –0.314 | –0.195 | –0.387 | –0.215 | –0.128 |
| T (P) | 0.68 | 0.001 | 0.115 | 0.006 | 0.244 | 0.028 | 0.009 | 0.305 | 0.2 | 0.08 | 0.972 | 0.005 | 0.002 | 0.06 | <0.001 | 0.038 | 0.22 |
| C (r) | –0.319 | 0.035 | –0.049 | –0.025 | –0.033 | 0.132 | –0.074 | 0.109 | 0.09 | 0.178 | –0.011 | 0.029 | –0.212 | –0.192 | 0.038 | –0.042 | –0.074 |
| C (P) | 0.002 | 0.735 | 0.638 | 0.809 | 0.753 | 0.204 | 0.479 | 0.296 | 0.388 | 0.086 | 0.916 | 0.781 | 0.041 | 0.063 | 0.715 | 0.689 | 0.48 |
MEST-C: Oxford classification [mesangial hypercellularity (M), endocapillary hypercellularity (E), segmental sclerosis (S), tubular atrophy/interstitial fibrosis (T), and crescent formation (C)]. ALB, albumin; CEUS, contrast-enhanced ultrasound; eGFR, estimated glomerular filtration rate; FT, fall time; mTT, mean transit time; PE, peak enhancement; RT, rise time; SCr, serum creatinine; TTP, time to peak; UACR, urine albumin-to-creatinine ratio; WiAUC, wash-in area under the curve; WiP, wash-in perfusion index; WiR, wash-in rate; WiWoAUC, wash-in and wash-out area under the curve; WoAUC, wash-out area under the curve; WoR, wash-out rate.
Construction of IgAN diagnosis model and nomogram
The diagnostic model for IgAN identified age, eGFR, UACR, and WiR as independent predictors (P<0.05; Table 5). A corresponding nomogram was developed based on these variables to estimate the probability of IgAN (Figure 2A). The model showed moderate discrimination, with an AUC of 0.78 (95% CI: 0.72–0.85; Figure 2B). Calibration and decision curve analyses showed acceptable agreement between predicted and observed probabilities and suggested potential clinical utility (Figure 2C,2D).
Table 5
| Variables | Univariate | Multivariate | |||||
|---|---|---|---|---|---|---|---|
| β | P | OR (95% CI) | β | P | OR (95% CI) | ||
| Age | −0.05 | <0.001 | 0.95 (0.93–0.98) | −0.06 | <0.001 | 0.94 (0.91–0.96) | |
| eGFR | −0.01 | 0.046 | 0.99 (0.99–0.99) | −0.01 | 0.005 | 0.99 (0.98–0.99) | |
| Urea | −0.05 | 0.150 | 0.95 (0.89–1.02) | ||||
| SCr | −0.00 | 0.856 | 1.00 (1.00–1.00) | ||||
| ALB | 0.05 | 0.023 | 1.05 (1.01–1.10) | ||||
| UACR | −0.01 | 0.007 | 0.99 (0.99–0.99) | −0.01 | 0.008 | 0.99 (0.99–0.99) | |
| PE | −0.59 | 0.008 | 0.56 (0.36–0.86) | ||||
| RT | 0.83 | 0.081 | 2.30 (0.90–5.87) | ||||
| mTT | 0.29 | 0.595 | 1.34 (0.46–3.94) | ||||
| TTP | 0.45 | 0.350 | 1.58 (0.61–4.08) | ||||
| FT | −1.34 | 0.041 | 0.26 (0.07–0.94) | ||||
| WiR | −0.48 | 0.002 | 0.62 (0.46–0.84) | −0.49 | 0.004 | 0.61 (0.44–0.86) | |
| WiP | −0.36 | 0.048 | 0.70 (0.49–0.99) | ||||
| WoR | −0.01 | 0.919 | 0.99 (0.81–1.22) | ||||
| WiAUC | −0.07 | 0.855 | 0.94 (0.47–1.88) | ||||
| WoAUC | −0.03 | 0.955 | 0.97 (0.36–2.64) | ||||
| WiWoAUC | 0.32 | 0.523 | 1.37 (0.52–3.63) | ||||
All CEUS parameters are presented as dimensionless cortex-to-medulla ratios. ALB, albumin; CEUS, contrast-enhanced ultrasound; CI, confidence interval; eGFR, estimated glomerular filtration rate; FT, fall time; IgAN, immunoglobulin A nephropathy; mTT, mean transit time; OR, odds ratio; PE, peak enhancement; RT, rise time; SCr, serum creatinine; TTP, time to peak; UACR, urine albumin-to-creatinine ratio; WiAUC, wash-in area under the curve; WiP, wash-in perfusion index; WiR, wash-in rate; WiWoAUC, wash-in and wash-out area under the curve; WoAUC, wash-out area under the curve; WoR, wash-out rate.
The combined model achieved an AUC of 0.78 (95% CI: 0.72–0.85), which was higher than that of the clinical model (AUC 0.64, 95% CI: 0.55–0.72) and the CEUS-only model (AUC 0.63, 95% CI: 0.55–0.71).
Pairwise comparisons using DeLong tests are shown in Table S1. The combined model was significantly superior to both the clinical model and the CEUS-only model, whereas no significant difference was observed between the clinical model and the CEUS-only model. These findings suggest that CEUS alone should not replace conventional clinical variables, but that it may provide incremental information when combined with clinical variables.
Construction of crescent formation diagnosis model and nomogram in IgAN, and subgroup analysis by renal function
For crescent formation prediction in IgAN, age and WiP were identified as independent predictors (P<0.05; Figure 3A). A nomogram was constructed to assess the risk of crescent formation in patients with IgAN. ROC analysis yielded an AUC of 0.74 (95% CI: 0.64–0.85; Figure 3B; Table 6), indicating moderate discrimination. Calibration and decision curve analyses showed acceptable agreement between predicted and observed probabilities and suggested potential clinical utility (Figure 3C,3D).
Table 6
| Variables | Univariate | Multivariate | |||||
|---|---|---|---|---|---|---|---|
| β | P | OR (95% CI) | β | P | OR (95% CI) | ||
| Age | −0.06 | 0.003 | 0.94 (0.91–0.98) | −0.07 | 0.001 | 0.93 (0.89–0.97) | |
| eGFR | 0.00 | 0.732 | 1.00 (0.99–1.01) | ||||
| Urea | −0.03 | 0.635 | 0.97 (0.85–1.10) | ||||
| SCr | −0.00 | 0.807 | 1.00 (0.99–1.01) | ||||
| ALB | −0.01 | 0.750 | 0.99 (0.92–1.06) | ||||
| UACR | 0.00 | 0.221 | 1.00 (1.00–1.00) | ||||
| PE | −0.25 | 0.475 | 0.78 (0.39–1.54) | ||||
| RT | 0.71 | 0.293 | 2.03 (0.54–7.64) | ||||
| mTT | 0.75 | 0.384 | 2.11 (0.39–11.35) | ||||
| TTP | 1.27 | 0.090 | 3.58 (0.82–15.59) | ||||
| FT | −0.09 | 0.915 | 0.91 (0.17–4.90) | ||||
| WiR | 0.06 | 0.779 | 1.07 (0.68–1.68) | ||||
| WiP | −0.56 | 0.043 | 0.57 (0.33–0.98) | −0.74 | 0.014 | 0.48 (0.26–0.86) | |
| WoR | −0.30 | 0.069 | 0.74 (0.54–1.02) | ||||
| WiAUC | 0.17 | 0.712 | 1.18 (0.48–2.91) | ||||
| WoAUC | −0.33 | 0.686 | 0.72 (0.14–3.59) | ||||
| WiWoAUC | −0.54 | 0.476 | 0.58 (0.13–2.56) | ||||
All CEUS parameters are presented as dimensionless cortex-to-medulla ratios. ALB, albumin; CEUS, contrast-enhanced ultrasound; CI, confidence interval; eGFR, estimated glomerular filtration rate; FT, fall time; IgAN, immunoglobulin A nephropathy; mTT, mean transit time; OR, odds ratio; PE, peak enhancement; RT, rise time; SCr, serum creatinine; TTP, time to peak; UACR, urine albumin-to-creatinine ratio; WiAUC, wash-in area under the curve; WiP, wash-in perfusion index; WiR, wash-in rate; WiWoAUC, wash-in and wash-out area under the curve; WoAUC, wash-out area under the curve; WoR, wash-out rate.
The AUC of the combined model was 0.74 (95% CI: 0.64–0.85), while the AUC of the clinical model was 0.64 (95% CI: 0.52–0.75) and that of the CEUS-only model was 0.62 (95% CI: 0.51–0.74). The CEUS-only model showed relatively higher sensitivity, while the clinical model had higher specificity. The combined model provided a more balanced performance profile, indicating that these two types of variables provide complementary information.
According to Table S1, the combined model was significantly superior to both the clinical model and the CEUS-only model. No significant difference was observed between the clinical model and the CEUS-only model, suggesting that the added value of CEUS lies mainly in its combination with clinical markers rather than in its independent use.
Subgroup analysis showed that nomogram performance was higher in patients with preserved renal function (eGFR ≥60 mL/min/1.73 m2) and lower in patients with advanced CKD. This pattern may be explained by diffuse and homogeneous microvascular impairment in advanced CKD, which reduces perfusion heterogeneity detectable by CEUS.
Construction of crescent formation diagnosis model and nomogram in the entire cohort
For crescent formation prediction in the entire cohort, age and RT were identified as independent predictors (P<0.05; Table 7; Figure 4A). A nomogram was constructed to assess the risk of crescent formation across the entire cohort. ROC analysis yielded an AUC of 0.72 (95% CI: 0.64–0.81; Figure 4B). Calibration and decision curve analyses further confirmed good agreement between predicted and observed probabilities and supported the clinical utility of the model (Figure 4C,4D).
Table 7
| Variables | Univariate | Multivariate | |||||
|---|---|---|---|---|---|---|---|
| β | P | OR (95% CI) | β | P | OR (95% CI) | ||
| Age | −0.05 | <0.001 | 0.95 (0.92–0.97) | −0.05 | <0.001 | 0.95 (0.92–0.97) | |
| eGFR | −0.00 | 0.725 | 1.00 (0.99–1.01) | ||||
| Urea | −0.08 | 0.100 | 0.92 (0.83–1.02) | ||||
| SCr | −0.00 | 0.748 | 1.00 (1.00–1.00) | ||||
| ALB | 0.02 | 0.348 | 1.02 (0.97–1.08) | ||||
| UACR | −0.00 | 0.246 | 1.00 (1.00–1.00) | ||||
| PE | −0.37 | 0.150 | 0.69 (0.41–1.14) | ||||
| RT | 1.09 | 0.041 | 2.97 (1.04–8.43) | 1.09 | 0.045 | 2.98 (1.03–8.45) | |
| mTT | −0.21 | 0.736 | 0.81 (0.23–2.79) | ||||
| TTP | 0.91 | 0.112 | 2.48 (0.81–7.58) | ||||
| FT | −0.96 | 0.198 | 0.38 (0.09–1.66) | ||||
| WiR | −0.14 | 0.405 | 0.87 (0.62–1.21) | ||||
| WiP | −0.38 | 0.075 | 0.68 (0.45–1.04) | ||||
| WoR | −0.21 | 0.129 | 0.81 (0.62–1.06) | ||||
| WiAUC | 0.25 | 0.549 | 1.28 (0.57–2.86) | ||||
| WoAUC | −0.40 | 0.507 | 0.67 (0.21–2.18) | ||||
| WiWoAUC | −0.18 | 0.754 | 0.84 (0.28–2.54) | ||||
All CEUS parameters are presented as dimensionless cortex-to-medulla ratios. ALB, albumin; CEUS, contrast-enhanced ultrasound; CI, confidence interval; eGFR, estimated glomerular filtration rate; FT, fall time; mTT, mean transit time; OR, odds ratio; PE, peak enhancement; RT, rise time; SCr, serum creatinine; TTP, time to peak; UACR, urine albumin-to-creatinine ratio; WiAUC, wash-in area under the curve; WiP, wash-in perfusion index; WiR, wash-in rate; WiWoAUC, wash-in and wash-out area under the curve; WoAUC, wash-out area under the curve; WoR, wash-out rate.
The AUC of the combined model was 0.72 (95% CI: 0.64–0.81), which was higher than that of the clinical model (0.64, 95% CI: 0.55–0.73) and the CEUS-only model (0.60, 95% CI: 0.51–0.69). Compared with the IgAN subgroup, the overall discrimination of the entire cohort was slightly lower. This reduction in performance is likely attributable to the heterogeneous pathophysiological mechanisms underlying crescent formation across different renal diseases. Nevertheless, the combined model remained the best-performing method.
DeLong tests showed that the combined model was significantly superior to the CEUS-only model and also outperformed the clinical model, with statistical significance or marginal significance depending on the comparison (Table S1). These results suggest that CEUS may add information when crescent formation is evaluated across a broader disease spectrum, although the magnitude of improvement remained moderate.
Comparison of model performance at each endpoint
Table 8 compares the diagnostic performance of the clinical, CEUS-only, and combined models across all three endpoints, including AUC, sensitivity, specificity, PPV, and NPV. Overall, the combined model consistently achieved the highest AUC than the clinical-only and CEUS-only models, supporting the incremental value of the integrated modeling strategy. ROC curves for the three models at each endpoint are displayed in Figure 5A-5C, with curves labeled as Combined, Clinical, and CEUS models. The renal-function subgroup results are summarized in Table 9, and the corresponding ROC curves are shown in Figure 6.
Table 8
| Endpoint | Model | AUC (95% CI) | Sensitivity | Specificity | PPV | NPV |
|---|---|---|---|---|---|---|
| IgAN diagnosis | Combined model (age + eGFR + UACR + WiR) | 0.78 (0.72–0.85) | 0.59 | 0.86 | 0.81 | 0.69 |
| IgAN diagnosis | Clinical model (age + eGFR + UACR) | 0.64 (0.55–0.72) | 0.77 | 0.54 | 0.64 | 0.69 |
| IgAN diagnosis | CEUS model (WiR) | 0.63 (0.55–0.71) | 0.55 | 0.72 | 0.67 | 0.61 |
| Crescent in IgAN | Combined model (age + WiP) | 0.74 (0.64–0.85) | 0.77 | 0.62 | 0.69 | 0.74 |
| Crescent in IgAN | Clinical model (age) | 0.64 (0.52–0.75) | 0.49 | 0.75 | 0.66 | 0.60 |
| Crescent in IgAN | CEUS model (WiP) | 0.62 (0.51–0.74) | 0.78 | 0.49 | 0.61 | 0.69 |
| Crescent in entire cohort | Combined model | 0.72 (0.64–0.81) | 0.74 | 0.66 | 0.49 | 0.85 |
| Crescent in entire cohort | Clinical model | 0.64 (0.55–0.73) | 0.47 | 0.77 | 0.47 | 0.77 |
| Crescent in entire cohort | CEUS model | 0.60 (0.51–0.69) | 0.64 | 0.56 | 0.39 | 0.78 |
AUC, area under the curve; CEUS, contrast-enhanced ultrasound; CI, confidence interval; eGFR, estimated glomerular filtration rate; IgAN, immunoglobulin A nephropathy; NPV, negative predictive value; PPV, positive predictive value; UACR, urine albumin-to-creatinine ratio; WiP, wash-in perfusion index; WiR, wash-in rate.
Table 9
| Subgroup | AUC | 95% CI | P value |
|---|---|---|---|
| eGFR ≥60 (mL/min/1.73 m2) | 0.788 | 0.664–0.913 | <0.001 |
| eGFR <60 (mL/min/1.73 m2) | 0.694 | 0.510–0.879 | 0.049 |
AUC, area under the curve; CI, confidence interval; eGFR, estimated glomerular filtration rate.
The corresponding pairwise statistical comparisons are summarized in Table S1. In summary, these analyses support the conclusion that CEUS provides additional information beyond conventional clinical variables at multiple clinically relevant endpoints, but the moderate AUC values indicate that the models should be considered supportive rather than definitive diagnostic tools.
Discussion
The main finding of this study is that CEUS may add diagnostic information when integrated with conventional clinical variables. Quantitative CEUS parameters (PE, WiR, and WiP) were associated with Oxford MEST-C scores, suggesting that CEUS can capture microvascular perfusion changes related to IgAN, particularly in fibrotic (T) and crescentic (C) lesions (6,7,11). By integrating these parameters with serological markers, we developed nomograms for IgAN diagnosis (AUC 0.78) and crescent prediction (AUC 0.74). These AUC values indicate moderate discrimination rather than high diagnostic accuracy, and the models should therefore be interpreted as exploratory supportive tools. Differences in the selected predictors across models may reflect endpoint-specific pathophysiological mechanisms and were determined by multivariable regression. The pre-test probability of IgAN in this cohort was relatively high (51.1%), exceeding that reported in many Western populations. This difference may reflect regional epidemiology and biopsy indications, and it may limit the generalizability of the proposed nomogram in populations with a lower prevalence of IgAN. Future multicenter studies should include more diverse populations for external validation. Reduced WiP in crescentic cases may reflect capillary disruption and consequent ischemia, which can contribute to progressive fibrosis. The normalized CEUS metrics used in this study may reduce inter-individual variability and strengthen links with histopathology compared with non-normalized measurements (4,5,11). Previous studies have also developed radiomic or clinical nomograms for predicting crescent formation in IgAN (12,13). Similarly, a CEUS-based nomogram has been reported for evaluating glomerulosclerosis in transplanted kidneys, supporting the feasibility of integrating quantitative CEUS metrics for non-invasive renal pathological assessment (14).
Another important consideration is the heterogeneity of the non-IgAN control group. Crescent formation across different renal diseases may arise through different pathological mechanisms, and renal perfusion patterns may also vary by disease type. This heterogeneity makes biological interpretation more complex and may partly explain why the crescent model performed slightly less well in the entire cohort than in the IgAN subgroup. Nevertheless, the combined model remained the best-performing strategy, suggesting that CEUS may have broader exploratory utility beyond a single pathological subtype, although disease-specific models may provide higher discrimination.
From a clinical perspective, the proposed nomograms may be useful in selected settings, but several practical constraints should be considered. The models are not intended to replace renal biopsy, which remains the diagnostic gold standard. Instead, they may assist pre-biopsy risk stratification and patient selection. Compared with radiomics, magnetic resonance imaging (MRI)-based approaches, and established clinical prediction models, CEUS has the advantages of real-time microvascular assessment and avoidance of nephrotoxic contrast agents. However, it also requires specific equipment, trained operators, standardized acquisition, and post-processing expertise. Cost, availability, and operator dependence may limit routine implementation, particularly outside centers with CEUS experience. For IgAN diagnosis, the model may help identify CKD patients with suspected glomerular disease who have a higher probability of biopsy-confirmed IgAN. For crescent prediction, the models may help identify patients at increased risk of active pathological injury. CEUS may also have potential value for longitudinal assessment or relapse monitoring in selected patients, although biopsy remains indispensable when diagnostic confirmation is required.
The subgroup analysis by renal function provided clinically relevant information. Better performance in patients with preserved renal function suggests that CEUS may be particularly useful in earlier disease stages, when microvascular changes remain heterogeneous enough to distinguish pathological states. In advanced CKD, diffuse and overlapping perfusion impairment may reduce the specificity of CEUS-derived parameters. This pattern supports the potential role of CEUS as an auxiliary tool in earlier rather than later CKD stages. However, blood pressure, disease duration, and other factors affecting renal microcirculation may still have influenced CEUS parameters and should be considered in future prospective studies.
This study has several strengths, including a diverse CKD cohort, standardized CEUS quantification, and internal validation. Several limitations should also be acknowledged. First, this was a retrospective single-center study, which may introduce selection bias and limit generalizability because patient selection and biopsy indications may differ between hospitals. Second, the sample size was relatively small for prediction-model development, especially for crescent and renal-function subgroup analyses. Third, no external validation was performed, and the nomograms require independent confirmation in larger multicenter cohorts. Fourth, formal interobserver or intraobserver reproducibility analysis for ROI placement and CEUS parameter extraction was not performed, although ROI placement was reviewed by a second investigator. Fifth, missing data were handled using complete-case analysis, which may introduce bias if missingness was not random. Sixth, because of the cross-sectional design, the prognostic value of CEUS for long-term renal outcomes could not be evaluated. Future prospective studies should include external validation, reproducibility assessment, detailed calibration statistics, optimism-corrected performance estimates, and prespecified clinical threshold analyses.
Conclusions
CEUS provides microvascular perfusion information that may complement conventional clinical biomarkers. In this exploratory single-center study, models combining CEUS and clinical variables showed better performance than clinical-only and CEUS-only models for diagnosing IgAN and predicting crescent formation, although discrimination was moderate. These findings suggest that CEUS may serve as a promising supportive tool for pre-biopsy risk stratification and crescent formation prediction in CKD, but external validation is required before routine clinical use.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the TRIPOD reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0025/rc
Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0025/dss
Funding: None.
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-0025/coif). The 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. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. This study was approved by the Ethics Committee of Huashan Hospital, Fudan University (Approval No. KY2024-730), and all patients provided informed consent.
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
- Roberts IS, Cook HT, Troyanov S, Alpers CE, Amore A, et al. The Oxford classification of IgA nephropathy: pathology definitions, correlations, and reproducibility. Kidney Int 2009;76:546-56.
- Cattran DC, Coppo R, Cook HT, Feehally J, Roberts IS, et al. The Oxford classification of IgA nephropathy: rationale, clinicopathological correlations, and classification. Kidney Int 2009;76:534-45.
- Trimarchi H, Barratt J, Cattran DC, Cook HT, Coppo R, Haas M, Liu ZH, Roberts IS, Yuzawa Y, Zhang H, Feehally J. IgAN Classification Working Group of the International IgA Nephropathy Network and the Renal Pathology Society; Conference Participants. Oxford Classification of IgA nephropathy 2016: an update from the IgA Nephropathy Classification Working Group. Kidney Int 2017;91:1014-21. [Crossref] [PubMed]
- Piscaglia F, Nolsøe C, Dietrich CF, Cosgrove DO, Gilja OH, Bachmann Nielsen M, et al. The EFSUMB Guidelines and Recommendations on the Clinical Practice of Contrast Enhanced Ultrasound (CEUS): update 2011 on non-hepatic applications. Ultraschall Med 2012;33:33-59. [Crossref] [PubMed]
- Sidhu PS, Cantisani V, Dietrich CF, Gilja OH, Saftoiu A, Bartels E, et al. The EFSUMB Guidelines and Recommendations for the Clinical Practice of Contrast-Enhanced Ultrasound (CEUS) in Non-Hepatic Applications: Update 2017 (Long Version). Ultraschall Med 2018;39:e2-e44. [Crossref] [PubMed]
- Wang L, Wu J, Cheng JF, Liu XY, Ma F, Guo LH, Xu JM, Wu T, Mohan C, Peng A, Xu HX, Song YX. Diagnostic value of quantitative contrast-enhanced ultrasound (CEUS) for early detection of renal hyperperfusion in diabetic kidney disease. J Nephrol 2015;28:669-78. [Crossref] [PubMed]
- Cao W, Cui S, Yang L, Wu C, Liu J, Yang F, Liu Y, Bin J, Hou FF. Contrast-Enhanced Ultrasound for Assessing Renal Perfusion Impairment and Predicting Acute Kidney Injury to Chronic Kidney Disease Progression. Antioxid Redox Signal 2017;27:1397-411. [Crossref] [PubMed]
- Greis C. Quantitative evaluation of microvascular blood flow by contrast-enhanced ultrasound (CEUS). Clin Hemorheol Microcirc 2011;49:137-49. [Crossref] [PubMed]
- Tranquart F, Mercier L, Frinking P, Gaud E, Arditi M. Perfusion quantification in contrast-enhanced ultrasound (CEUS)--ready for research projects and routine clinical use. Ultraschall Med 2012;33:S31-8. [Crossref] [PubMed]
- Platz Batista da Silva N, Jung EM, Jung F, Schlitt HJ, Hornung M. VueBox® perfusion analysis of contrast-enhanced ultrasound (CEUS) examinations in patients with primary hyperparathyroidism for preoperative detection of parathyroid gland adenoma. Clin Hemorheol Microcirc 2018;70:423-31. [Crossref] [PubMed]
- Han BH, Park SB. Usefulness of Contrast-enhanced Ultrasound in the Evaluation of Chronic Kidney Disease. Curr Med Imaging 2021;17:1003-9. [Crossref] [PubMed]
- Qin X, Xia L, Hu X, Xiao W, Huaming X, Xisheng X, Zhang C. A novel clinical-radiomic nomogram for the crescent status in IgA nephropathy. Front Endocrinol (Lausanne) 2023;14:1093452. [Crossref] [PubMed]
- Lin Z, Feng L, Zeng H, Lin X, Lin Q, Lu F, Wang L, Mai J, Fang P, Liu X, Tan Q, Zou C. Nomogram for the prediction of crescent formation in IgA nephropathy patients: a retrospective study. BMC Nephrol 2023;24:262. [Crossref] [PubMed]
- Xu N, Wang D, Hong Y, Huang P, Xu Q, Sun H, Cai L, Yin J, Zhang L, Yang B. A nomogram based on contrast-enhanced ultrasound for evaluating the glomerulosclerosis rate in transplanted kidneys. Quant Imaging Med Surg 2024;14:3060-74. [Crossref] [PubMed]

