Skip to content

Original Research

Frontiers in Nephrology ·

Volume 6 - 2026 | https://doi.org/10.3389/fneph.2026.1826486

Development and validation of a routine urinalysis-based risk stratification model for IgA nephropathy in a multicenter real-world hospital population

  • Chi Wang1,2,3
  • Wan He1,2,3
  • QingYuan Zheng1,2,3
  • Yuan Zhong1,2,3
  • Meng Tang1,2,3
  • ChunYing Zhang1,2,3*corresponding author
  • Yong He1,2,3*corresponding author
  1. 1Department of Laboratory Medicine, West China Hospital, Sichuan University, Chengdu, China
  2. 2Sichuan Clinical Research Center for Laboratory Medicine, Chengdu, China
  3. 3Clinical Laboratory Medicine Research Center of West China Hospital, Chengdu, China
Article metrics0 viewsView details

Abstract

IgA nephropathy (IgAN) is the most common primary glomerular disease worldwide. Routine urinalysis is universally accessible but underutilized for IgAN risk stratification. This study aimed to develop and validate a urinalysis-based risk stratification model for IgAN in a multicenter real-world hospital population. This retrospective multicenter study was based on routine urinalysis records from 8 hospitals. A total of 354 patients had biopsy-confirmed primary IgAN, and 18,002 controls had no documented IgAN diagnosis after electronic medical record (EMR) review. The dataset was randomly divided into a training cohort (70%) and a hold-out validation cohort (30%). LASSO penalized logistic regression was used for variable selection. Model performance was evaluated by discrimination, calibration, and decision curve analysis (DCA). Bootstrap validation with 1000 resamplings was performed. An 8-variable simplified model was additionally developed. LASSO regression selected 13 variables for the primary full model. The primary model achieved an AUC of 0.849 (95% CI: 0.824–0.873) in the training cohort and 0.861 (95% CI: 0.833–0.890) in the validation cohort. Calibration was good in the training cohort (Brier score = 0.017), while moderate calibration deviation was observed in the validation cohort (Brier score = 0.020). The model yielded a high negative predictive value (>0.99). The 8-variable simplified model retained favorable discrimination with a validation AUC of 0.840. This routine urinalysis-based model may serve as an accessible laboratory-based risk stratification tool for identifying patients with relatively low or elevated predicted probability of IgAN in real-world hospital settings. It should be interpreted as an adjunct to clinical evaluation rather than as a stand-alone diagnostic test or replacement for renal biopsy.

Keywords

1 Introduction

Immunoglobulin A nephropathy (IgAN), first described by Berger and Hinglais in 1968, is the most prevalent primary glomerular disease worldwide and a leading cause of end-stage renal disease (ESRD) among young and middle-aged adults (1, 2). Characterized by predominant IgA1 deposition in the glomerular mesangium with mesangial cell proliferation and matrix expansion, IgAN exhibits substantial clinical and pathological heterogeneity (3, 4). Approximately 20-30% of patients progress to ESRD within 10–20 years following diagnosis, highlighting the importance of early clinical recognition and risk assessment (5).

The gold standard for IgAN confirmation remains renal biopsy, a procedure that is invasive, technically demanding, and associated with potential complications (e.g., bleeding), particularly in patients with bleeding tendencies or skin infections (6, 7). Additionally, while galactose-deficient IgA1 (Gd-IgA1)-specific antiglycan antibodies show promise as diagnostic markers, their clinical application is not yet widespread (8). The Kidney Disease Improving Global Outcomes (KDIGO) guidelines acknowledge the lack of reliable serum or urinary screening tools for IgAN, highlighting an unmet clinical need for non-invasive, accessible diagnostic tools (9, 10).

Urinalysis represents one of the most commonly performed laboratory examinations in clinical medicine. Unlike other biological fluids regulated by homeostatic mechanisms, urine serves as a reservoir that reflects and accumulates alterations in internal homeostasis, making it particularly suitable for detecting early disease-related changes (11, 12). Urine contains abundant metabolic products that provide valuable information regarding physiological status and disease manifestation, facilitating screening of renal diseases, urinary tract infections, and metabolic disorders (12, 13). Despite its widespread use, the composite predictive value of routine urinalysis parameters for IgAN risk stratification remains underexplored. Previous studies have largely focused on isolated markers or limited clinical cohorts, lacking the scale and methodological rigor required for real-world implementation, highlighting the importance of early clinical recognition and risk assessment (14, 15).

Therefore, this multicenter hospital laboratory-based study aimed to develop and validate a routine urinalysis-based risk stratification model for IgAN. The model was designed to support laboratory-based identification of patients with relatively elevated predicted probability of IgAN in real-world hospital urinalysis settings, rather than to provide a stand-alone diagnosis or replace renal biopsy.

2 Materials and methods

2.1 Study design and population

This was a multicenter cross-sectional study conducted among outpatients and inpatients from 8 hospitals (West China Hospital of Sichuan University, Shanghai Fengxian District Central Hospital, Shanghai General Hospital, The First Hospital of Jilin University, University-Town Affiliated Hospital of Chongqing Medical University, West China Second University Hospital, Sichuan University, The Second Affiliated Hospital of Kunming Medical University and The First Affiliated Hospital of Henan University of CM). The study period was from October 2023 to December 2025. During routine clinical practice, urine specimens from outpatients and inpatients were submitted to the laboratory for standard urinalysis, and the corresponding results were recorded in the laboratory information system.

Urinalysis records were linked with EMRs using unique patient identifiers. IgAN cases were defined as patients with biopsy-confirmed primary IgAN, as determined by renal pathology reports and EMR review. Controls were defined as patients who underwent routine urinalysis during the same period and had no documented diagnosis of IgAN after review of EMRs, diagnostic records, nephrology records, and available pathology reports. When multiple urinalysis records were available for the same patient, one eligible index record was retained to avoid repeated measurements from the same individual. When available, urinalysis records obtained before disease-specific treatment were preferentially selected for biopsy-confirmed IgAN cases. This study was approved by the institutional review board, and the requirement for informed consent was waived.

2.2 Data collection and laboratory measurements

Urine samples were collected from eligible participants during outpatient visits or hospitalization. All samples underwent standardized routine urinalysis using the CA-500 Automatic Urine Analyzer and MA-500 Automated Urine Sediment Analyzer (MEDICSIDE Corp, China), ensuring methodological consistency across all participating centers. For samples yielding inconclusive automated results, manual microscopic examination was performed by experienced laboratory technicians to confirm findings. Patients receiving renin-angiotensin-aldosterone system inhibitors (RAASi) or immunosuppressive agents prior to urinalysis were excluded to prevent confounding from treatment-modified urinalysis profiles. A total of 50 routine urine test variables were recorded, including:

Dipstick parameters: Urobilinogen (UBG), bilirubin (BIL), ketone bodies (KET), occult blood (BLD), protein (PRO), nitrite (NIT), leukocyte esterase (LEU), glucose (GLU), pH, vitamin C (VC), microalbumin (MALB), creatinine (Cr), urinary calcium (CA), albumin/creatinine ratio (A:C), specific gravity (SG), conductivity (Cond), osmolality (Osm).

Urinary sediment quantitative parameters: Red blood cells (RBC), normal RBC, microcytes, acanthocytes, ghost cells, other abnormal RBCs, leukocytes, leukocyte clusters, squamous epithelial cells, non-squamous epithelial cells, renal tubular epithelial cells, caudate epithelial cells, small round epithelial cells, hyaline casts, pathological casts, granular casts, waxy casts, broad casts, other casts, bacilli, total yeasts, yeasts, candida, crystals, calcium oxalate monohydrate crystals, calcium oxalate dihydrate crystals, uric acid crystals, magnesium ammonium phosphate crystals, amorphous salt crystals, other crystals, sperm, mucus threads, suspected cocci.

Laboratory measurements were performed according to CLSI guidelines and the operational manuals of the analyzers. Dipstick semi-quantitative results were converted to ordinal variables for statistical analysis (0 = negative, 1 = 1+, 2 = 2+, 3 = 3+, 4 = 4+, 5 = 5+), while urinary sediment parameters were used directly as continuous variables (reported as counts per high-power field [HPF] or low-power field [LPF]).

2.3 Statistical analysis

2.3.1 Data preprocessing and standardization

The dataset was randomly divided into a training cohort (70%) and a hold-out validation cohort (30%). Continuous variables were standardized using z-scores to improve coefficient interpretability and numerical stability. The mean and standard deviation used for standardization were estimated exclusively from the training cohort and then applied unchanged to the hold-out validation cohort to avoid information leakage.

2.3.2 Missing data handling

Missing values were observed only for “coccilike bacteria” (Supplementary Table 2), with 35 missing observations (0.19% of this variable, ~0.004% of all data cells). Median imputation was applied. This variable was not retained by LASSO and was not included in the final model.

2.3.3 LASSO variable selection

LASSO-penalized logistic regression with 10-fold cross-validation was employed (16). The optimal penalty parameter λ was selected using the one-standard-error (1-SE) rule (λ = 0.002675191). All variables with non-zero coefficients at the selected λ were retained for the final model.

2.3.4 Model construction

Variables selected by LASSO were entered into a multivariable logistic regression model. P values were reported for descriptive purposes only; no additional variable elimination was performed based on post-LASSO P-values. Variance inflation factor (VIF) was calculated (threshold <5). To improve clinical applicability, an 8-variable simplified model using only dry chemistry parameters was additionally constructed as an implementation-oriented secondary model.

2.3.5 Model validation and performance evaluation

Model performance was evaluated by discrimination (AUC), calibration (Hosmer-Lemeshow test, Brier score, calibration intercept and slope), and clinical utility (DCA). Diagnostic metrics including sensitivity, specificity, PPV, NPV, F1 score, and balanced accuracy were calculated. Bootstrap validation with 1000 resamplings was performed to estimate AUC stability. All analyses were performed using R software (version 4.3.1).

3 Results

3.1 Baseline characteristics

Study participants were enrolled and processed as outlined in Figure 1. A total of 18,356 hospital-based patients undergoing routine urinalysis were included, comprising 354 biopsy-confirmed IgAN cases and 18,002 controls without documented IgAN diagnosis after EMR review. The training cohort included 12,850 patients, of whom 230 had biopsy-confirmed IgAN, whereas the hold-out validation cohort included 5,506 patients, of whom 124 had biopsy-confirmed IgAN. Key urinalysis characteristics related to the final models are summarized in Table 1, and the full distribution of all 50 candidate urinalysis variables is provided in Supplementary Table 1.

Study flow diagram of the multicenter hospital laboratory-based routine urinalysis cohort. Routine urinalysis records were collected from outpatients and inpatients at eight participating hospitals. Urinalysis records were linked with electronic medical records and available renal pathology reports. IgAN cases were defined as patients with biopsy-confirmed primary IgAN, whereas controls were defined as patients without documented IgAN diagnosis after EMR review. The final cohort included 354 biopsy-confirmed IgAN cases and 18,002 controls without documented IgAN diagnosis. The cohort was randomly divided into a training cohort and a hold-out validation cohort. EMR, electronic medical record; IgAN, IgA nephropathy; LASSO, least absolute shrinkage and selection operator.
Figure 1. Study flow diagram of the multicenter hospital laboratory-based routine urinalysis cohort. Routine urinalysis records were collected from outpatients and inpatients at eight participating hospitals. Urinalysis records were linked with electronic medical records and available renal pathology reports. IgAN cases were defined as patients with biopsy-confirmed primary IgAN, whereas controls were defined as patients without documented IgAN diagnosis after EMR review. The final cohort included 354 biopsy-confirmed IgAN cases and 18,002 controls without documented IgAN diagnosis. The cohort was randomly divided into a training cohort and a hold-out validation cohort. EMR, electronic medical record; IgAN, IgA nephropathy; LASSO, least absolute shrinkage and selection operator.

Table 1. Key urinalysis characteristics related to the final models according to IgAN status.

CharacteristicTotal cohort (N = 18,356)Controls without documented IgAN (N = 18,002)Biopsy-confirmed IgAN (N = 354)P value
Outcome group
Participants, n (%)18,356 (100.00)18,002 (98.07)354 (1.93)—
Dipstick-based predictors
Occult blood, grade, median [IQR]0 [0; 1]0 [0; 1]1 [0; 3]<0.001
Protein, grade, median [IQR]0 [0; 2]0 [0; 2]3 [1; 3]<0.001
Leukocyte esterase, grade, median [IQR]0 [0; 1]0 [0; 1]0 [0; 0]<0.001
Glucose, grade, median [IQR]0 [0; 0]0 [0; 0]0 [0; 3]<0.001
Ketone bodies, grade, median [IQR]0 [0; 0]0 [0; 0]0 [0; 0]0.005
pH, median [IQR]5.5 [5.5; 6.0]5.5 [5.5; 6.5]5.5 [5.0; 6.0]<0.001
Vitamin C, grade, median [IQR]0 [0; 1]0 [0; 1]1 [0; 1]<0.001
Specific gravity, median [IQR]1.02 [1.01; 1.02]1.02 [1.01; 1.02]1.02 [1.01; 1.02]0.436
Additional predictors in the primary 13-variable model
Microalbumin, grade, median [IQR]1 [0; 3]1 [0; 3]3 [3; 3]<0.001
Creatinine, grade, median [IQR]2 [1; 3]2 [1; 3]3 [2; 3]<0.001
Squamous epithelial cells, median [IQR]1.00 [0.00; 2.00]1.00 [0.00; 2.00]1.00 [0.00; 3.00]0.437
Conductivity, mS/cm, median [IQR]16.60 [11.90; 21.70]16.60 [11.90; 21.80]13.90 [9.70; 18.85]<0.001
Acanthocytes, median [IQR]0.00 [0.00; 0.00]0.00 [0.00; 0.00]0.00 [0.00; 1.00]<0.001

Values are presented as median [interquartile range] unless otherwise indicated. Dipstick semi-quantitative parameters were coded as ordinal grades. P values compare biopsy-confirmed IgAN cases and controls without documented IgAN diagnosis. Mann-Whitney U test was used for continuous or ordinal variables, and chi-square test was used for categorical variables. IgAN, IgA nephropathy; IQR, interquartile range.

3.2 Predictor selection

Among the 50 candidate variables, LASSO regression retained 13 variables for the primary full urinalysis model (Figure 2; Table 2). The events per variable (EPV) was 17.7 for the 13-variable model and 28.8 for the 8-variable model.

Predictor selection using LASSO penalized logistic regression. **(A)** LASSO coefficient profiles of the 50 candidate routine urinalysis variables. **(B)** Ten-fold cross-validation curve for selection of the penalty parameter λ. Thirteen variables with non-zero coefficients at the selected λ (λ = 0.002675191) were retained for the primary full urinalysis model. The vertical dotted lines indicate λmin and λ1se. LASSO, least absolute shrinkage and selection operator.
Figure 2. Predictor selection using LASSO penalized logistic regression. (A) LASSO coefficient profiles of the 50 candidate routine urinalysis variables. (B) Ten-fold cross-validation curve for selection of the penalty parameter λ. Thirteen variables with non-zero coefficients at the selected λ (λ = 0.002675191) were retained for the primary full urinalysis model. The vertical dotted lines indicate λmin and λ1se. LASSO, least absolute shrinkage and selection operator.

Table 2. Post-LASSO logistic regression coefficients of the primary 13-variable full urinalysis model.

VariableβSEOR95% CIP valueVIF
Intercept-5.5130.760.0040.001–0.018<0.001–
Microalbumin0.6290.0961.8761.555–2.264<0.0011.676
Creatinine0.4260.0781.5311.313–1.784<0.0011.544
Occult blood0.3690.0461.4461.321–1.584<0.0011.23
Vitamin C0.3320.0771.3941.199–1.620<0.0011.077
Glucose0.2390.0431.271.168–1.382<0.0011.462
Squamous epithelial cells0.1290.041.1371.052–1.2290.0011.037
Acanthocytes0.0180.021.0180.978–1.0600.3871.045
Protein0.0270.0651.0270.904–1.1670.6831.813
Ketone bodies-0.8430.270.430.253–0.7310.0021.04
Leukocyte esterase-0.6240.0950.5360.445–0.645<0.0011.096
Specific gravity-0.4420.1030.6430.526–0.786<0.0011.803
pH-0.2530.1160.7760.619–0.9730.0281.076
Conductivity-0.1520.0850.8590.728–1.0140.0731.366

Continuous variables were standardized using z-scores before model fitting; odds ratios for standardized continuous variables correspond to a 1-standard-deviation increase. Semi-quantitative dipstick variables were entered as ordinal variables. P values are provided for descriptive purposes only and were not used for additional variable elimination after LASSO selection. VIF, variance inflation factor. The LASSO penalty parameter used was λ = 0.002675191.

3.3 Performance of the primary 13-variable model

3.3.1 Discrimination

The primary 13-variable model showed an AUC of 0.849 (95% CI: 0.824–0.873) in the training cohort and 0.861 (95% CI: 0.833–0.890) in the hold-out validation cohort (Figure 3). All VIF values were below 5.0.

Receiver operating characteristic curves of the primary 13-variable full urinalysis model. **(A)** Training cohort: AUC = 0.849, 95% CI: 0.824–0.873. **(B)** Hold-out validation cohort: AUC = 0.861, 95% CI: 0.833–0.890. AUC, area under the receiver operating characteristic curve; CI, confidence interval.
Figure 3. Receiver operating characteristic curves of the primary 13-variable full urinalysis model. (A) Training cohort: AUC = 0.849, 95% CI: 0.824–0.873. (B) Hold-out validation cohort: AUC = 0.861, 95% CI: 0.833–0.890. AUC, area under the receiver operating characteristic curve; CI, confidence interval.

3.3.2 Calibration

Calibration was good in the training cohort (Brier score = 0.017; Hosmer-Lemeshow P = 0.649; intercept = 0.000; slope = 1.000). In the validation cohort, the Brier score remained low (0.020), but the Hosmer-Lemeshow test reached statistical significance (P = 0.017), indicating that the model is more suitable for relative risk stratification than for precise estimation of absolute probability (Figure 4).

Calibration curves of the primary 13-variable full urinalysis model. **(A)** Training cohort. **(B)** Hold-out validation cohort. The diagonal line represents ideal calibration. The primary model showed excellent calibration in the training cohort, whereas some calibration deviation was observed in the validation cohort despite a low Brier score.
Figure 4. Calibration curves of the primary 13-variable full urinalysis model. (A) Training cohort. (B) Hold-out validation cohort. The diagonal line represents ideal calibration. The primary model showed excellent calibration in the training cohort, whereas some calibration deviation was observed in the validation cohort despite a low Brier score.

3.3.3 Clinical utility

Decision curve analysis (DCA) showed that the primary 13-variable model provided a positive net benefit compared with the treat-all and treat-none strategies mainly within the low-threshold probability range from 0 to 0.20 (Figure 5).

Decision curve analysis of the primary 13-variable full urinalysis model. **(A)** Training cohort. **(B)** Hold-out validation cohort. Decision curves were displayed across threshold probabilities from 0 to 0.20 because IgAN prevalence was low in this hospital-based urinalysis cohort and low threshold probabilities were considered more relevant for laboratory-based risk stratification. DCA, decision curve analysis; IgAN, IgA nephropathy.
Figure 5. Decision curve analysis of the primary 13-variable full urinalysis model. (A) Training cohort. (B) Hold-out validation cohort. Decision curves were displayed across threshold probabilities from 0 to 0.20 because IgAN prevalence was low in this hospital-based urinalysis cohort and low threshold probabilities were considered more relevant for laboratory-based risk stratification. DCA, decision curve analysis; IgAN, IgA nephropathy.

3.3.4 Diagnostic performance

The model achieved a high negative predictive value (NPV > 0.99) in both cohorts, while the positive predictive value (PPV) was low (~6-8%) due to the low population prevalence of IgAN (Table 3).

Table 3. Performance metrics of the primary 13-variable full model.

MetricTraining cohortValidation cohort
Discrimination
AUC0.8490.861
95% CI (DeLong)0.824–0.8730.833–0.890
Bootstrap AUC (1000 resamples)0.8410.857
Calibration
Brier score0.0170.020
Calibration intercept0.0000.465
Calibration slope1.0001.060
Hosmer-Lemeshow P value0.6490.017
Diagnostic performance
Sensitivity0.7650.790
Specificity0.7870.794
Positive predictive value0.0610.081
Negative predictive value0.9950.994

3.4 Simplified 8-variable model

To further address clinical applicability, we developed a simplified dipstick-based model using eight routinely available parameters: occult blood, protein, leukocyte esterase, glucose, ketone bodies, pH, vitamin C, and specific gravity. The simplified model showed favorable discrimination, with an AUC of 0.827 in the training cohort and 0.840 in the validation cohort (Supplementary Figure 2; Supplementary Table 4). Calibration was good in the training cohort (Brier = 0.017, Hosmer-Lemeshow P = 0.3135), whereas less consistent calibration was observed in the validation cohort despite a low Brier score (Brier = 0.021, Hosmer-Lemeshow P = 0.0004). Therefore, this model was considered an implementation-oriented secondary model rather than a replacement for the primary full model (Table 4; Supplementary Figure 3).

Table 4. Performance metrics of the simplified 8-variable model.

MetricTraining cohortValidation cohort
Discrimination
AUC0.8270.840
95% CI (DeLong)0.802–0.8520.811–0.869
Bootstrap AUC (1000 resamples)0.8210.840
Calibration
Brier score0.0170.021
Hosmer-Lemeshow P value0.31350.0004
Diagnostic performance
Sensitivity0.7700.750
Specificity0.7670.775
Positive predictive value0.0570.071
Negative predictive value0.9950.993
F1 score0.1060.130
Balanced accuracy0.7680.762

3.5 Bootstrap internal validation

Bootstrap validation with 1000 resamplings yielded comparable AUC estimates (0.841 in training, 0.857 in validation), suggesting no obvious evidence of severe overfitting (Supplementary Figure 1; Supplementary Table 2).

4 Discussion

4.1 Principal findings

This multicenter laboratory-based urinalysis study developed a 13-variable full model and an 8-variable simplified implementation model for IgAN risk stratification in a real-world hospital setting.

4.2 Clinical positioning and intended use

It is imperative to emphasize that this model should be used for laboratory-based risk stratification rather than stand-alone diagnosis or rule-out decision-making. Due to the low prevalence of IgAN, the model yields a low PPV but a high NPV (>0.99). The high NPV should be interpreted in the context of low disease prevalence. It may help identify patients with relatively low predicted probability, but it should not be used alone to exclude IgAN when clinical suspicion remains high.

4.3 Interpretation of predictors

Classic IgAN-related markers (occult blood, microalbumin) reflect glomerular injury (17, 18). Non-classical variables (vitamin C, ketone bodies, specific gravity) may reflect urine concentration, metabolic status, tubular handling, or comorbid conditions, and their inclusion should be interpreted primarily as incremental predictive information rather than direct mechanistic evidence (6, 19).

4.4 Clinical applicability of the simplified model

The simplified model retained favorable discrimination and may support implementation in settings where only commonly available dipstick urinalysis parameters are available. However, validation calibration was less consistent, so absolute probabilities should be interpreted cautiously.

4.5 Limitations

First, controls were defined by the absence of documented IgAN rather than universal biopsy verification, which inherently introduces verification bias and may include occult IgAN cases. Second, the low prevalence of IgAN resulted in a low PPV. Third, the model cannot be used to independently rule out IgAN. Fourth, statistically significant calibration deviation was observed in the validation cohort for both models. Fifth, clinical variables (age, sex, BP, eGFR) were not incorporated. Sixth, exclusion of patients receiving RAAS inhibitors or immunosuppressive therapy may reduce confounding from treatment-modified urinalysis profiles, but it may also limit generalizability to treated populations. Although a hold-out validation cohort and bootstrap validation were used, both were derived from the same multicenter retrospective dataset; therefore, the present study does not constitute prospective external validation. Prospective external validation with standardized diagnostic follow-up is warranted.

5 Conclusion

In this multicenter hospital laboratory-based routine urinalysis cohort, we developed and validated a 13-variable LASSO-selected full urinalysis model for IgAN risk stratification. The model showed favorable discrimination and low overall prediction error, although calibration deviation was observed in the validation cohort. An 8-variable simplified dipstick-based model retained acceptable discrimination and may support implementation in settings where comprehensive urinalysis parameters are unavailable. These models should be interpreted as adjunctive risk stratification tools and should not replace clinical evaluation or renal biopsy.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneph.2026.1826486/full#supplementary-material

Supplementary material Supplementary material at the source

Supplementary material Supplementary material at the source

Supplementary material Supplementary material at the source

Supplementary material Supplementary material at the source

Supplementary material Supplementary material at the source

Supplementary material Supplementary material at the source

Supplementary material Supplementary material at the source

Footnotes

Edited by: Vladimir Tesar, Charles University, Czechia

Reviewed by: Irene De Lourdes Noronha, University of São Paulo, Brazil

Soumita Bagchi, All India Institute of Medical Sciences, India

Statements

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the clinical research project commissioned by Zhejiang Medicside Medical Technology Co., Ltd., entitled 'Development and Multicenter Validation of Rules for AI-Assisted Automated Review of Urinalysis and Clinical Disease Alerts' (Project No. HX-H2309304), and by the Science and Technology Agency of Sichuan Province (Grant No. 2024YFFK0144).

Data availability

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics approval

The study was approved by the ethics committee of West China Hospital of Sichuan University (Approval number: NO. 20231272). The studies were conducted in accordance with the local legislation and institutional requirements. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

CW: Validation, Writing – review & editing, Methodology, Investigation, Writing – original draft. WH: Writing – review & editing, Data curation, Investigation, Visualization. QZ: Software, Writing – review & editing, Resources. YZ: Data curation, Project administration, Methodology, Writing – review & editing. MT: Writing – review & editing, Supervision, Resources. CZ: Investigation, Supervision, Writing – review & editing, Methodology, Data curation. YH: Writing – review & editing, Conceptualization, Resources, Funding acquisition, Visualization.

Acknowledgements

We thank all participating hospitals and laboratory staff for their support in data collection and quality control.

References

  1. 1 Knoppova B, Reily C, King RG, Julian BA, Novak J, Green TJ. Pathogenesis of IgA nephropathy: Current understanding and implications for development of disease-specific treatment. J Clin Med. (2021) 10: 4501. doi: doi:10.3390/jcm10194501
  2. 2 Jiawen H, Tangli X, Jun Z, Yan L, Liangyu Y, Chun Z. Ambrisentan as adjunctive therapy for IgA nephropathy: a retrospective single-center analysis. J Nephrol. (2026) 39 (2): 335–43. doi: doi:10.1093/joneph/aajaf052
  3. 3 Pattrapornpisut P, Avila-Casado C, Reich HN. IgA nephropathy: Core curriculum 2021. Am J Kidney Dis. (2021) 78: 429–41. doi: doi:10.1053/j.ajkd.2021.01.024
  4. 4 Fang Z, Yang Y, Fei T, Yebei L, Daijin R, Chengmin C. Efficacy and safety of telitacicept combined with immunosuppressive therapy for IgA nephropathy: a retrospective multicenter cohort study. Front Immunol. (2026) 17. doi: doi:10.3389/fimmu.2026.1740891
  5. 5 Ito T, Akamatsu K. Echocardiographic manifestations in end-stage renal disease. Heart Fail Rev. (2023) 29: 465–78. doi: doi:10.1007/s10741-023-10376-5
  6. 6 Xu Z, Zhan H, Zhang J, Li Z, Cheng L, Chen Q. New biomarkers in IgA nephropathy. Clin Immunol. (2025) 274: 110468. doi: doi:10.1016/j.clim.2025.110468
  7. 7 Dan X, Minjie Z, Weiwei L, Lijiang F, Feifei G. Risk factors of disease progression in IgA nephropathy: a systematic review and meta-analysis. Immun Inflammation Dis. (2026) 14: e70393. doi: doi:10.2215/cjn.01430307
  8. 8 Gleeson PJ, O'Shaughnessy MM, Barratt J. IgA nephropathy in adults-treatment standard. Nephrol Dial Transplant. (2023) 38: 2464–73. doi: doi:10.1093/ndt/gfad146
  9. 9 Floege J, Barratt J, Cook HT, Noronha IL, Reich HN, Suzuki Y. Executive summary of the KDIGO 2025 clinical practice guideline for the management of immunoglobulin A nephropathy (IgAN) and immunoglobulin A vasculitis (IgAV). Kidney Int. (2025) 108: 548–54. doi: doi:10.1016/j.kint.2025.04.003
  10. 10 Balázs S, Tibor V, Sadra S, Tibor József K. Prognostic significance of lactate dehydrogenase-to-albumin ratio and neutrophil percentage-to-albumin ratio in IgA nephropathy. Biomedicines. (2026) 14: 318. doi: doi:10.1093/ndt/11.5.906b
  11. 11 Haq K, Patel DM. Urinalysis: Interpretation and clinical correlations. Med Clin North Am. (2023) 107: 659–79. doi: doi:10.1016/j.mcna.2023.03.002
  12. 12 Pelin A, Ahmet A. Urinary system involvement in children with cystic fibrosis: a single-center retrospective cohort. Pediatr Pulmonol. (2026) 61: e71551. doi: doi:10.1002/ppul.71551
  13. 13 Vorperian SK, DeFelice BC, Buonomo JA, Chinchinian HJ, Gray IJ, Yan J. Deconvolution of human urine across the transcriptome and metabolome. Clin Chem. (2024) 70: 1344–54. doi: doi:10.1093/clinchem/hvae137
  14. 14 Xu L-L, Zhang D, Weng H-Y, Wang L-Z, Chen R-Y, Chen G. Machine learning in predicting T-score in the Oxford classification system of IgA nephropathy. Front Immunol. (2023) 14: 1224631. doi: doi:10.3389/fimmu.2023.1224631
  15. 15 Yang S, Lin Y, Chen Z. Predictive model integrating imaging and inflammatory biomarkers for preoperative urinary tract infection in patients with ureteral calculi and metabolic syndrome. Arch Esp Urol. (2026) 79: 105–13. doi: doi:10.56434/j.arch.esp.urol.20267901.13
  16. 16 Najafabadi AH, Ramspek CL, Dekker FW, Heus P, Hooft L, Moons KGM. TRIPOD statement: a preliminary pre-post analysis of reporting and methods of prediction models. BMJ Open. (2020) 10: e041537. doi: doi:10.1136/bmjopen-2020-041537
  17. 17 Kamal F, Kim J, Lafayette R. Current biomarkers of IgA nephropathy. Semin Nephrol. (2024) 44: 151572. doi: doi:10.1016/j.semnephrol.2025.151572
  18. 18 Ramsawak S, Cohen S, Linares A, Cavanaugh C. IgA nephropathy: Update on pathogenesis and treatment. Cleve Clin J Med. (2025) 92: 373–83. doi: doi:10.3949/ccjm.92a.24105
  19. 19 Novak J, Reily C, Steers NJ, Schumann T, Rizk DV, Julian BA. Emerging biochemical and immunologic mechanisms in the pathogenesis of IgA nephropathy. Semin Nephrol. (2024) 44: 151565. doi: doi:10.1016/j.semnephrol.2025.151565

Citation

Wang, C., et al. (2026). Development and validation of a routine urinalysis-based risk stratification model for IgA nephropathy in a multicenter real-world hospital population. Frontiers in Nephrology, Volume 6. https://doi.org/
Received
Accepted
Published
Volume
6 - 2026

Copyright

Copyright © 2026 Wang, He, Zheng, Zhong, Tang, Zhang and He.

License: https://creativecommons.org/licenses/by/4.0/

*Correspondence: ChunYing Zhang; Yong He

Article metrics

Views
0
Downloads
0
Citations
0
Shares
0