A non-invasive diagnostic model of immunoglobulin A nephropathy and serological markers for evaluating disease severity.

BACKGROUND: Immunoglobulin A nephropathy (IgAN) is the most common pathological type of glomerular disease. Kidney biopsy, the gold standard for IgAN diagnosis, has not been routinely applied in hospitals worldwide due to its invasion nature. Thus, we aim to establish a non-invasive diagnostic model and determine markers to evaluate disease severity by analyzing the serological parameters and pathological stages of patients with IgAN.
METHODS: A total of 272 biopsy-diagnosed IgAN inpatients and 518 non-IgA nephropathy inpatients from the Department of Nephrology of Chinese People's Liberation Army General Hospital were recruited for this study. Routine blood examination, blood coagulation testing, immunoglobulin-complement testing, and clinical biochemistry testing were conducted and pathological stages were analyzed according to Lee grading system. The serological parameters and pathological stages were analyzed. The receiver operating characteristic (ROC) analysis was performed to estimate the diagnostic value of the clinical factors. Logistic regression was used to establish the diagnostic model.
RESULTS: There were 15 significantly different serological parameters between the IgAN and non-IgAN groups (all P < 0.05). The ROC analysis was performed to measure the diagnostic value for IgAN of these parameters and the results showed that the area under the ROC curve (AUC) of total protein (TP), total cholesterol (TC), fibrinogen (FIB), D-dimer (D2), immunoglobulin A (IgA), and immunoglobulin G (IgG) were more than 0.70. The AUC of the "TC + FIB + D2 + IgA + age" combination was 0.86, with a sensitivity of 85.98% and a specificity of 73.85%. Pathological grades of I, II, III, IV, and V accounted for 2.21%, 17.65%, 62.50%, 11.76%, and 5.88%, respectively, with grade III being the most prevalent. The levels of urea nitrogen (UN) (13.57 ± 5.95 vs. 6.06 ± 3.63, 5.92 ± 2.97, 5.41 ± 1.73, and 8.41 ± 3.72 mmol/L, respectively) and creatinine (Cr) (292.19 ± 162.21 vs. 80.42 ± 24.75, 103.79 ± 72.72, 96.41 ± 33.79, and 163.04 ± 47.51 μmol/L, respectively) were significantly higher in grade V than in the other grades, and the levels of TP (64.45 ± 7.56, 67.16 ± 6.94, 63.22 ± 8.56, and 61.41 ± 10.86 vs. 37.47 ± 5.6 mg/d, respectively), direct bilirubin (DB) (2.34 ± 1.23, 2.58 ± 1.40, 1.91 ± 0.97, and 1.81 ± 1.44 vs. 0.74 ± 0.57 μmol/L, respectively), and IgA (310.35 ± 103.78, 318.48 ± 107.54, 292.58 ± 81.85, and 323.29 ± 181.67 vs. 227.17 ± 68.12 g/L, respectively) were significantly increased in grades II-V compared with grade I (all P < 0.05).
CONCLUSIONS: The established diagnostic model that combined multiple factors (TC, FIB, D2, IgA, and age) might be used for IgAN non-invasive diagnosis. TP, DB, IgA, Cr, and UN have the potential to be used to evaluate IgAN disease severity.

Introduction

Immunoglobulin A nephropathy (IgAN), first discovered by Berger in 1968, is the most common pathological type of glomerular disease. This disease is featured by the deposition of immunoglobulin A (IgA) in the mesangium. The clinical spectrum of IgAN contains a wide range of features, ranging from asymptomatic microscopic haematuria with/without proteinuria to gross haematuria, and few patients present with clinical signs of nephrotic or nephritic syndrome. In recent years, the proportion of IgAN in kidney disease has exceeded 50% in the world.[,,] Approximately 20% to 40% patients with IgAN could develop end-stage renal disease (ESRD) within 20 years after first detection. Therefore, early identification of risk factors and predicting IgAN prognosis are of great value. Nowadays, the detection of IgAN mainly depends on the microscopic examination of renal biopsies, with immunocytochemical techniques providing further confirmation of the diagnosis results. Biopsy, which is essential for validating the diagnosis, may cause severe complications. Therefore, it is important to develop a non-invasive diagnostic model for IgAN and to find markers for evaluating the disease severity and progress of IgAN.

Statistics have been applied to determine significant predictors for the diagnosis or classification of various diseases. Subsequently, different statistical algorithms, biological datasets, and parameters may result in different outputs.[,,] More importantly, multicollinearity is always used in medical laboratory parameters, which may also cause variability and instability in a statistical model. Thus, it is of great importance to choose appropriate variable for multiparameter analyses. In recent years, modeling techniques have been widely applied in medicine to help monitor disease progression and predict treatment outcomes. However, no such approach exists for IgAN.[,] On the other hand, in the past decades, considerable efforts have been dedicated to detecting useful clinical markers for the non-invasive diagnosis of IgAN, and several markers have been proposed.[,] To date, several biomarkers have been suggested, but their validity has not been demonstrated in clinical practice.[,,] Some studies have found several risk factors consistently associated with the progression of IgAN, such as proteinuria >1 g/day, arterial hypertension, reduced renal function at diagnosis and glomerular sclerosis or tubulointerstitial scarring at renal biopsy.[,,] However, the influence of other factors such as older age, male sex, overweight, obesity, hypertriglyceridemia or hyperuricemia on IgAN development remains controversial.[,,]

Thus, our study aims to establish a non-invasive diagnostic model and determine markers for evaluating disease severity by analyzing serological parameters and pathological stages in patients with IgAN.

Methods

Ethics approval

The study was approved by the Ethics Committee of the Chinese People's Liberation Army General Hospital. All patients provided informed written consent for study sample collection as well as permission for their use in research.

Study design

A total of 790 inpatients from the Department of Nephrology of Chinese People's Liberation Army General Hospital were enrolled in our study at Chinese People's Liberation Army General Hospital from January 1, 2013, to March 15, 2016. Among them, a total of 272 patients were confirmed to have IgAN and 518 patients were confirmed to be non-IgAN (other types of nephropathy) by renal biopsy. These patients did not undergo any treatments before they were diagnosed and the patients’ complete history, clinical information, and pathological data were collected.

The inclusion criteria were as follows: (1) the patient accepted a renal biopsy during their hospitalization and (2) receiving no renal biopsy prior to pathological diagnosis at our hospital. The exclusion criteria used for the final selection of cases were as follows: (1) the renal biopsy was not conducted; (2) the patient received immunosuppression treatment or renal replacement therapy; (3) the pathological results indicated that the patient had secondary kidney disease, including diabetic nephropathy, lupus nephritis, and hepatitis-related nephropathy; (4) the clinical data were incomplete; and (5) patients with inflammatory and infectious diseases that may cause secondary IgA nephropathy. Based on the exclusion criteria, 790 cases were finally selected from the Chinese People's Liberation Army General Hospital.

Samples and biological parameters

Blood samples of 790 patients were collected for routine blood examination, blood coagulation testing, immunoglobulin-complement testing, and clinical biochemistry testing. Pathological stages were analyzed according to Lee's grading system. Serum samples were obtained by collecting venous blood into Vacutainer serum separator tubes (Greiner Bio-One, Frickenhausen, Germany), and the samples were then centrifuged at 1500 × g at 4°C for 10 min in a centrifuge (Sigma-Aldrich, St. Louis, MO, USA) within 20 min of collection. The blood biochemical and routine blood examinations were performed using enzymatic assays (Roche Products Ltd., Basel, Switzerland) on a fully automatic biochemical autoanalyzer (Cobas 8000; Roche Products Ltd., Basel, Switzerland).

Statistical analysis

The normally distributed data are expressed as the mean ± standard deviation (SD) and were compared using unpaired Student's t tests. The non-normally distributed data are expressed as medians with the corresponding 25th and 75th percentiles (interquartile range) and compared using Mann-Whitney U tests. The categorical variables were analyzed using the Chi-square tests. Logistic regression analysis was employed to establish the diagnostic model. Receiver operating characteristic (ROC) analysis was performed to measure the diagnostic value of the clinical factors and the area under the ROC curve (AUC) more than 0.70 were considered to have a good specificity. A value of P < 0.05 was considered to indicate a significant difference. All statistical analyses were performed using SPSS software version 19.0 (SPSS Inc., Chicago, IL, USA) and GraphPad Prism software vision 6 (GraphPad Software Inc., San Diego, CA, USA).

Results

Patient characteristics

Clinical characteristics of the participants in this study are shown in Table 1. The differences in the serological parameters between the two groups are shown in Table 2. There were 15 significantly different serological parameters between the IgAN and non-IgAN groups, including total protein (TP), total bilirubin (TB), direct bilirubin (DB), creatinine (Cr), uric acid (Ua), total cholesterol (TC), triglyceride (TG), lactate dehydrogenase (LDH), high density lipoprotein cholesterol (HDL), low density lipoprotein cholesterol (LDL), fibrinogen (FIB), D-dimer (D2), immunoglobulin A (IgA), immunoglobulin G (IgG), and immunoglobulin E (IgE) (all P < 0.05). However, urea nitrogen (UN), creatine kinase (CK), immunoglobulin M (IgM), complement 3 (C3), and complement 4 (C4) were not statistically different between the 2 groups.

Table 1

Clinical characteristics of the participants in this study, n = 790.

Table 2

Differences in the serological parameters between the IgAN and non-IgAN groups.

ROC analysis of related characteristics

ROC analysis was performed to measure the diagnostic value for IgAN of the 15 different serological parameters between the IgAN and non-IgAN groups. The area under the curves (AUCs) of TP (AUC = 0.78), TC (AUC = 0.73), FIB (AUC = 0.74), D2 (AUC = 0.71), IgA (AUC = 0.74), and IgG (AUC = 0.71) were more than 0.70 [Table 3].

Table 3

Receiver operating characteristic analysis for the diagnostic value of the 15 clinical parameters for immunoglobulin A nephropathy.

Diagnostic model based on logistic regression analysis

Logistic regression analysis was used to establish the diagnostic model for IgAN patients. The results showed that the high level of IgA, low level of TC, FIB, D2, and young age were risk indicators for IgAN [Table 4]. The predicted probabilities (PRE) were calculated based on logistic regression analysis. The classification equation for predicting IgAN was as follows: PRE = 1/(1+e–[2.389–0.237×TC–0.457×FIB–0.397×D2+0.008×IgA–0.043×age]).

Table 4

Multivariate logistic regression analysis for immunoglobulin A nephropathy.

Based on the PRE of the logistic regression analysis, the AUC of the “TC + FIB + D2 + IgA + age” combination was 0.86 (P < 0.001, 95% CI: 0.83–0.89), with a sensitivity and specificity of 85.98% and 73.85%, respectively [Figure 1].

Figure 1

The diagnostic value of the “TC + FIB + D2 + IgA + age” combination for immunoglobulin A nephropathy as detected using receiver operating characteristic analysis. The AUC was 0.86, with a sensitivity of 85.98% and a specificity of 73.85%. AUC: area under the receiver operating characteristic curve; D2: D-dimer; FIB: Fibrinogen; IgA: Immunoglobulin A; ROC: Receiver operating characteristic; TC: Total cholesterol.

Relationship between the pathological grades and the laboratory test results

We collected several representative clinical parameters to analyze the relationship between the pathological grades and laboratory test results. As shown in Table 5, the levels of UN and Cr were significantly higher in grade V than in grades II–V, respectively (all P < 0.05). And the levels of TP, DB, and IgA were significantly increased in grades II–V compared with grade I (all P < 0.05).

Table 5

Pathological grades and laboratory test results for the immunoglobulin A nephropathy patients.

Discussion

IgAN has a high incidence worldwide and a large proportion of IgAN patients will progress to ESRD. Therefore, it is of great value to evaluate the disease severity and development of IgAN. Currently, although several IgAN biomarkers have been extensively researched, none have been applied for its screening in clinical practice. Moreover, a number of computational studies have been performed on many types of kidney disease, but none have focused on noninvasively diagnosing IgAN. In the present study, we used retrospective data to analyze the serological parameters and pathological stages of patients with IgAN and to establish a noninvasive diagnostic model for IgAN.

Based on the statistical analyses and clinical experience, 15 out of 20 routine and useful parameters were selected as predictors of IgAN: TP, TB, DB, Cr, Ua, TC, TG, LDH, HDL, LDL, FIB, D2, IgA, IgG, and IgE [Table 2]. Compared with previous studies, this study included more characteristics, including fibrinogen, D-dimer, serum IgA, and C3, all of which are known biomarkers of kidney disease.[,] In addition, ROC analysis was performed to further assess the diagnostic value of the 15 parameters and the results showed that the AUCs of TP, TC, FIB, D2, IgA, and IgG were all more than 0.7 [Table 3]. Although serum IgA appears not to be a specific biomarker of IgAN, previous studies have reported that IgA levels are still statistically different and have differentially diagnostic value, especially when combined with other clinical parameters. Berthoux et al reported that IgG was a biomarker for the prediction of clinicopathologic recurrence events in IgAN. Our non-IgAN group consists mainly of membranous nephropathy, minimally pathological nephropathy, mesangial proliferative glomerulonephritis and other diseases characterized by nephrotic syndrome. Patients with nephrotic syndrome are often in a state of hypercoagulability, hyperfibrinolysis, hyperlipemia, and hypoproteinemia. This may explain why the relevant index levels for blood clotting and blood lipids were higher in the non-IgAN group than in the IgAN group, such as FIB, D2, and TC, while the TP levels were lower in the non-IgAN group than in the IgAN group.

The multivariate logistic regression analysis requires that each explanatory variable is independent. Based on our clinical experience, the level of TP is not independent of IgA and IgG levels. Therefore, we removed TP and selected the other 5 parameters for further analysis. Given that the incidence of IgAN varies from age to age, we added the age variable to the logistic regression model. The predicted probabilities were calculated based on logistic regression analysis [Table 4]. The AUC of the “TC + FIB + D2 + IgA + age” combination was 0.86, with a sensitivity of 85.98% and a specificity of 73.85% [Figure 1]. The established diagnostic model that combined multiple factors (TC, FIB, D2, IgA, and age) might be used for IgAN noninvasive diagnosis.

We found that the levels of UN and Cr were significantly higher in grade V patients than those in other grades. Furthermore, apparent increases in TP, DB, and IgA were observed in grades II–V compared with grade I [Table 5]. Many studies have shown that elevated Ua, serum Cr and other indicators are associated with an increased risk of IgAN.[,] However, there are limited reports on the relationship between these blood indicators and pathological grades or other indicators reflecting the severity of the disease.

Several strengths of our study should be stated. First, readily available clinical parameters such as patient demographics were applied. Second, all clinical characteristics were derived from biopsy-proven patients with IgAN. These patients were probably representative of patients with increased diagnostic uncertainty, which is the most challenging patients encountered in clinical practice. Lastly, our models were internally validated. However, this study has a few limitations that must be considered. First, this study was not a longitudinal investigation but rather a cross-sectional study. We are unable to determine the impact of these altered indicators on the pathogenesis of IgAN and whether this elevation is progressive or reversible. Second, because the study individuals were all of Chinese descent, we cannot ensure that our research results are applicable to individuals of other ethnic backgrounds. Third, this was a single-center study; further multicenter studies and large cohort studies should be conducted for validation.

In conclusion, the established diagnostic model that combined multiple factors (TC, FIB, D2, IgA, and age) could effectively distinguish IgAN patients from non-IgAN patients, with high sensitivity and specificity. TP, DB, IgA, Cr, and UN could be used to evaluate IgAN disease severity.

Funding

This work was supported by grants from the National Key R&D Program of China (No. 2016YFC1305500), the National Natural Science Foundation of China (Nos. 61471399, 61671479, and 81670663), the National Key Research and Development Program (No. 2016YFC1305404), and the Joint Funds of National Natural Science Foundation of China and Henan Province (No. U1604284).

Conflicts of interest

None.