NDT Advance Access published online on June 10, 2009
Nephrology Dialysis Transplantation, doi:10.1093/ndt/gfp273
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
A scoring system to predict renal outcome in IgA nephropathy: a nationwide 10-year prospective cohort study
1 Kyoto University Health Service, Kyoto 2 Department of Preventive Medicine/Biostatistics and Medical Decision Making, Nagoya University Graduate School of Medicine, Nagoya 3 Division of Nephrology and Metabolism, Department of Internal Medicine, Tokai University School of Medicine, Isehara 4 Division of Nephrology, Department of Internal Medicine, Juntendo University School of Medicine, Tokyo, Japan
Correspondence and offprint requests to: Masashi Goto; E-mail: goto{at}msa.biglobe.ne.jp
| Abstract |
|---|
|
|
|---|
Background. Immunoglobulin A nephropathy (IgAN) is the most common form of glomerulonephritis, and a substantial number of patients succumb to end-stage renal disease (ESRD). However, prediction of the renal outcome in individual patients remains difficult. We have already published a scoring system using the data in a prospective cohort of IgAN patients followed up from 1995 to 2002.
Methods. The cohort was further followed up until 2005 in 97 clinical units in Japan. The data from 2283 patients were analysed by Cox regression to determine the predictors of ESRD in IgAN, and their β-coefficients were converted into scores to estimate ESRD risk within 10 years.
Results. During the follow-up (median, 87 months), 252 patients developed ESRD. Male sex, age less than 30 years, family histories of chronic renal failure and chronic glomerulonephritis, hypertension, proteinuria, mild haematuria, hypoalbuminaemia, low glomerular filtration rate and a high histological grade at initial renal biopsy were associated with the risk of ESRD in the multivariable analysis. A scoring system was framed to estimate the 10-year ESRD risk using eight variables significant in both univariable and multivariable models. This prognostic score accurately classified patients by risk: patients with estimates of 0–4.9, 5.0–19.9, 20.0–49.9 and 50.0–100% had an observed incidence of 1.7, 8.3, 36.7 and 85.5%, respectively. The corresponding area under the receiver-operating characteristic curve was 0.942 (95% confidence interval, 0.925–0.958).
Conclusion. This validated scoring system to quantitatively estimate ESRD risk during the 10-year follow-up of IgAN patients will serve as a useful prognostic tool in clinical practice.
Keywords: cohort studies; IgA nephropathy; prognosis; renal dialysis; risk factors
| Introduction |
|---|
|
|
|---|
Immunoglobulin A nephropathy (IgAN) was described as a new clinical entity in 1968 by Berger and Hinglais [1], and is now the most common cause of idiopathic glomerulonephritis [2–4]. Long-term outcomes and prognostic factors of patients with IgAN have been evaluated in many studies. Although this disorder is thought to follow a benign course, many patients are at a risk of at least slow progression. Furthermore, end-stage renal disease (ESRD) developed in
15% of cases within 10 years [5]. Several tools have been available to estimate the renal prognosis using various clinical and pathological factors [6–9]. However, one was based on a relatively small size cohort [6], another used surrogate endpoints such as an increase in serum creatinine [7] and the others did not provide quantitative estimates of the ESRD risks [8,9].
We earlier proposed a valid scoring system to quantitatively predict renal outcomes based on the 7-year follow-up data involving more than 2000 patients with biopsy-proven IgAN [10]. However, renal function deteriorates so slowly in the early stage of IgAN that the longer the follow-up period, the more useful in the clinical settings is the constructed prediction model [7]. We, therefore, extended the follow-up of this cohort for 3 more years and refined the scoring system to predict renal outcomes at 10 years.
| Methods |
|---|
|
|
|---|
The details of the methods used here were described in our earlier 7-year follow-up study [10]. Briefly, 2450 patients with biopsy-proven IgAN from 97 clinical units were followed up from 1995 when a nationwide survey on IgAN was jointly conducted by the two research committees for the specified intractable diseases organized by the Japanese Government. Follow-up mail surveys to collect information on outcomes such as death, ESRD and serum creatinine were performed in 1997, 1999 and 2002 with response rates of 82.5, 95.7 and 93.3%, respectively. An additional survey was carried out in 2005 (response rate, 82.7%).
The baseline data of the patients were obtained by reviewing medical records in the nationwide survey in 1995. The data included sex, age, family history of chronic renal failure and chronic glomerulonephritis, initial clinical manifestations, year of diagnostic renal biopsy, systolic and diastolic blood pressure, urine protein and blood, serum total protein and albumin and serum creatinine. Proteinuria was semiquantified with a standard urine dipstick with (–), (±), (+), (++) and (+++) corresponding to <10, 10–29, 30–99, 100–299 and
300 mg/dL of urine albumin, respectively. For estimation of the glomerular filtration rate (GFR), the estimation equation for Japanese patients with chronic kidney disease was applied [11]. This equation calculates the GFR from serum creatinine, age and gender. Histological grade at initial renal biopsy was assessed using the criteria from the joint committee of one of the aforementioned governmental research committees and the Japanese Society of Nephrology (Table 1) [12]. Information on therapy was collected in 1997, 2 years after the beginning of the follow-up.
|
Statistical methods
The study endpoint was ESRD defined as the initiation of chronic haemodialysis. Since there was little influence of competing risks in the current study, the complement of a Kaplan–Meier survival estimate was referred to as the 10-year cumulative incidence (risk) of ESRD [13]. Hazard ratios of potential prognostic factors were estimated using the Cox proportional hazards regression model. The independent effect of each variable was assessed by multivariable analysis. Systolic or diastolic blood pressure and serum total protein or albumin were included in the multivariable model because of their close relationship with each other. Akaike's Information Criterion (AIC) and the Bayesian Information Criterion (BIC) were used to select the most appropriate covariance structure.
Development of the scoring system
A scoring system to predict ESRD in individual patients with IgAN was composed based on the proportional hazards model including the significant variables in the aforementioned multivariable analysis. To enhance parsimony, only variables that possessed significant univariate relationships with ESRD were included in the final scoring system [14]. To generate a simple integer-based point score for each predictor variable, scores were given by multiplying the β-coefficient by 10 and rounding up or down to the nearest integer. The overall risk score for each patient was calculated by summing the scores of all components. The baseline survivor function was estimated by the Kaplan–Meier product–limit estimate.
Evaluation of the scoring system
To examine the goodness of fit of the scoring system to the data, we divided the patients into four groups according to the predicted 10-year risk of ESRD, that is, low (0–4.9%), moderate (5.0–19.9%), high (20.0–49.9%) and very high (50.0–100%). The renal survival curve of ESRD was then drawn in each group using the Kaplan–Meier method. To further assess the utility of the score, we used the area under the receiver-operating characteristic (ROC) curve for the 10-year risk of ESRD. The area and its 95% confidence interval (CI) were estimated by the nonparametric method.
As an additional analysis, one-third of the participants were randomly allocated to a validation sample and the remainder to a derivation sample. The prognostic score was developed in the derivation sample, and the actual 10-year cumulative incidence of ESRD was computed by the predicted risk in the validation sample. The area under the ROC curve was also estimated in the validation group. Considering the sampling error, we repeated this procedure in 100 different validation sets.
Statistical analyses were performed using the Stata 10.1 software (Stata Corporation, College Station, TX, USA). All tests of significance were two-tailed, and P-values <0.05 were considered statistically significant. This investigation was approved by the Ethics Committee of Kyoto University Graduate School of Medicine and the Ethics Committee of the Juntendo University School of Medicine.
| Results |
|---|
|
|
|---|
Participants demographic and clinical characteristics
Out of the 2450 patients tracked from 1995, 165 patients with unknown outcomes and 2 patients with erroneous baseline serum creatinine levels were excluded, leaving 2283 patients in the current analysis. The median age of the included patients at baseline was 32.1 years [interquartile range (IQR), 20.7–46.9], and 51.3% were women. During the follow-up of 14 975 person-years [median follow-up period, 87 months (IQR, 42–122)], 252 patients (11.0%) developed ESRD. The renal survival rate at 10 years was 85.0% (95% CI, 83.1–86.7). Twenty-one deaths without ESRD were also reported: six from circulatory diseases, six from neoplasms, two from other causes and seven from unknown causes. By 1997, 34.5%, 10.6% and 28.2% of the patients had received corticosteroids, immunosuppressive agents and angiotensin-converting enzyme inhibitors (ACEi), respectively.
Table 2 summarizes the 10-year cumulative incidence of ESRD by demographic and clinical factors and their hazard ratios. Male sex, advanced age, earlier renal biopsy, higher systolic and diastolic blood pressure, more severe proteinuria, lower serum total protein and albumin, lower GFR and higher histological grade were significantly associated with the risk of ESRD. On the other hand, patients whose initial clinical manifestation was macrohaematuria had favourable renal outcomes. Mild haematuria, 1–29 red blood cells per high-power field (RBCs/HPF), showed a higher risk than severe haematuria,
30 RBCs/HPF (18.2% versus 12.3%). There were no remarkable differences in risk between histological grade I and II (4.7% versus 6.2%).
|
Multivariable analysis and development of the scoring system
Results of the multivariable analysis for the risk of ESRD by factors are summarized in Table 3. Family histories of chronic renal failure and chronic glomerulonephritis were significant factors for predicting renal outcome in the multivariable model. Instead, initial manifestation of macrohaematuria and earlier renal biopsy were no longer significant. Male sex, higher systolic blood pressure, more severe proteinuria, mild haematuria, hypoalbuminaemia, lower GFR and higher histological grade were associated with the development of ESRD. Patients aged less than 30 years were at higher risk in the multivariable model, whereas an upward trend in the hazard ratio with advancing age was found in the univariate analysis.
|
Based on these analyses, we composed a scoring system to estimate a 10-year risk of ESRD. Because the presence of family histories of chronic renal failure and chronic glomerulonephritis had no significant univariate relationships with ESRD and the area under the ROC curve of the model was rather smaller than that without these variables, they were removed from the final prediction scheme. Table 4(a) lists the scores of individual prognostic factors. The sum of individual scores can then be converted into the corresponding estimated risk using those shown in Table 4(b). The baseline survivor function was estimated as 0.9993504 at 10 years.
|
Evaluation of the scoring system
The renal survival curves according to the estimated 10-year risk are shown in Figure 1. The prognostic score successfully classified the patients by risk. Those with an estimated risk of 0–4.9% (score, 0–43), 5.0–19.9% (44–58), 20.0–49.9% (59–70) and 50.0–100% (71 or more) had an actual cumulative incidence of ESRD in 10 years of 1.7% (95% CI, 0.86–3.5), 8.3% (5.4–12.5), 36.7% (28.4–46.5) and 85.5% (78.0–91.5), respectively. This indicated a good agreement between the estimated and observed risks. Furthermore, the area under the ROC curve of this prediction model reached 0.942 (95% CI, 0.925–0.958) (Figure 2).
|
|
Even when the prognostic scores were developed using derivation samples randomly selected from all participants, the estimated 10-year cumulative incidence of ESRD well predicted the observed ones in the remaining validation sample. The median values of observed 10-year incidence were 2.2% (IQR, 1.5–3.0%), 9.2% (7.4–11.8), 34.3% (30.4–38.6) and 83.4% (76.7–87.1) in patients with an estimated risk of 0–4.9, 5.0–19.9, 20.0–49.9 and 50.0–100%, respectively. The median of the corresponding area under the ROC curve (0.935; IQR, 0.924–0.944) was comparable with the area in the full dataset (0.942).
| Discussion |
|---|
|
|
|---|
Based on a large-scale cohort study, we described the prognostic indicators for IgAN and developed a scoring system for estimating the ESRD risk within 10 years. Male sex, age less than 30 years, the presence of family histories of chronic renal failure and chronic glomerulonephritis, higher systolic blood pressure, more severe proteinuria, mild haematuria, hypoalbuminaemia, lower GFR and higher histological grade were related to the risk. The prognostic score comprising eight variables significant in both univariable and multivariable models successfully classified patients according to their ESRD risk, and the accuracy in predicting the ESRD was excellent.
We could establish a much simpler scoring system than the former one based on the 7-year follow-up data [10] by reducing the number of choices of each scoring item. Nevertheless, we did not compromise with the predictability; the areas under the ROC curve estimating ESRD risk rather increased by 0.003 from 0.939 (95% CI, 0.921–0.958) [10].
Male sex was a significant risk factor for ESRD in the current model, whereas it made a favourable contribution in the former scoring. It would be attributable to the overestimation of the women's kidney function based only on serum creatinine in our previous analysis. We could estimate GFR more elaborately based on serum creatinine, as well as age and gender in the current analysis, which made our understanding of the relationships between each predictor much simpler. One of the reasons for the increase in the relative weight of the age variable in the current model compared to that in the previous one is also this overestimation of the baseline renal function among the older patients. Contrary to the age variable, the prognostic value of the histological grade at initial renal biopsy declined by extending the follow-up period.
We selected serum albumin in the current multivariable model, instead of serum total protein in the prior one, because both AIC and BIC of the model were substantially decreased by replacing serum total protein with serum albumin. This replacement seems rational from the pathophysiological viewpoint that glomerular proteinuria is mainly composed of albumin.
In the multivariable analysis, we found that the presence of family histories of chronic renal failure and chronic glomerulonephritis were associated with the development of ESRD. Several genetic factors are considered to be associated with the prognosis of IgAN. For example, some genes related to the renin–angiotensin system are suggested to have prognostical importance [15–17]. A report from China suggested an association between the polymorphism of the megsin gene, a kidney-specific serine protease inhibitor, and the progression of IgAN [18]. Furthermore, familial IgAN was reported to have a worse prognosis in an observational epidemiological study [19]. Nevertheless, the influences of genetic inheritance are still in controversy and have not reached a consensus [16,20–22].
In virtue of the current scoring system, clinicians can quantitatively estimate ESRD risks during 10 years of follow-up of IgAN patients by summing up only eight integers yielded from the routine clinico-pathological data. The validated, multivariable regression-based model should aid medical decision making about patients with biopsy-proved IgAN. Patients judged to be at higher risk can receive more frequent monitoring and more intensive treatment, whereas patients estimated to be at lower risk can be reassured and given less intensive care.
Given the nature of the observational design and the insufficient information regarding treatment, a clear evaluation of the influence of therapy is beyond the scope of the current analysis. Even though the optimal approach to the treatment of IgAN is uncertain, some strategies have been proved effective in reducing proteinuria and the rate of disease progression through the accumulation of evidence during the last decade. According to both observational studies [23,24] and randomized trials [25–27], ACEi and angiotensin II receptor blockers have a definite role in treating IgAN, particularly of the hypertensive and proteinuric forms. A review supports the use of corticosteroids and other immunosuppressive agents in reducing proteinuria or preventing progression to ESRD [28]. Increasing attention has been paid to the role of tonsillectomy for the long-term prognosis of IgAN [29,30]. The estimated risk of ESRD by the current scoring system, therefore, might be attenuated for patients properly treated in the latter portion of the observation interval.
Some other potential limitations of the current analysis must be acknowledged. First, because the data on 24-h urinary protein excretion were not available for two-thirds of the patients and those on urinary creatinine were not collected, we had to assess proteinuria with simple dipstick urine test results. However, the semiquantified proteinuria was reproducible and considerably correlated with the 24-h urinary excretion of protein [10,31], and the predictability of the current scoring system was excellent. Feasibility in clinical settings would override theoretical propriety. From a pragmatic standpoint, clinicians should average or take the median of the results of several urinalyses when applying this scoring system. Second, the current scoring system was developed among Japanese patients, and the applicability of the results to other races was not verified. Another study adjusting for racial differences in serum creatinine, which was a main determinant of GFR in the estimation equation, would assure the generalizability of the current scoring system.
In summary, the present study revealed that male sex, age less than 30 years, family histories of chronic renal failure and chronic glomerulonephritis, hypertension, proteinuria, mild haematuria, hypoalbuminaemia, low GFR and advanced histological changes increased the risk of ESRD in IgAN patients. The ESRD prediction score based on a multivariable model was sufficiently valid and will serve as a useful tool for clinicians treating IgAN patients.
| Acknowledgments |
|---|
The authors express their sincere appreciation to the physicians who participated in this study. This study was supported in part by a Grant-in-Aid for the Research Group on Progressive Renal Diseases and the Research Committee on Epidemiology of Intractable Diseases from the Ministry of Health, Labor and Welfare, Japan.
Conflict of interest statement. Seven-year follow-up data from this cohort have already been published in Nephrology Dialysis Transplantation, October, 2006. The authors declare no conflict of interest.
| References |
|---|
|
|
|---|
- Berger J, Hinglais N. Intercapillary deposits of IgA-IgG. J Urol Nephrol (Paris) (1968) 74:694–695.[Medline]
- Li LS, Liu ZH. Epidemiologic data of renal diseases from a single unit in China: analysis based on 13 519 renal biopsies. Kidney Int (2004) 66:920–923.[CrossRef][Web of Science][Medline]
- Simon P, Ramee MP, Boulahrouz R, et al. Epidemiologic data of primary glomerular diseases in western France. Kidney Int (2004) 66:905–908.[CrossRef][Web of Science][Medline]
- Barratt J, Feehally J. IgA nephropathy. J Am Soc Nephrol (2005) 16:2088–2097.
[Free Full Text] - DAmico G. Natural history of idiopathic IgA nephropathy: role of clinical and histological prognostic factors. Am J Kidney Dis (2000) 36:227–237.[Web of Science][Medline]
- Beukhof JR, Kardaun O, Schaafsma W, et al. Toward individual prognosis of IgA nephropathy. Kidney Int (1986) 29:549–556.[Web of Science][Medline]
- Alamartine E, Sabatier JC, Guerin C, et al. Prognostic factors in mesangial IgA glomerulonephritis: an extensive study with univariate and multivariate analyses. Am J Kidney Dis (1991) 18:12–19.[Web of Science][Medline]
- DAmico G, Colasanti G, Barbiano di Belgioioso G, et al. Long-term follow-up of IgA mesangial nephropathy: clinico-histological study in 374 patients. Semin Nephrol (1987) 7:355–358.[Web of Science][Medline]
- Frimat L, Briançon S, Hestin D, et al. IgA nephropathy: prognostic classification of end-stage renal failure. Nephrol Dial Transplant (1997) 12:2569–2575.
[Abstract/Free Full Text] - Wakai K, Kawamura T, Endoh M, et al. A scoring system to predict renal outcome in IgA nephropathy: from a nationwide prospective study. Nephrol Dial Transplant (2006) 21:2800–2808.
[Abstract/Free Full Text] - Matsuo S, Imai E, Horio M, et al. Revised equations for estimated GFR from serum creatinine in Japan. Am J Kidney Dis (2009) 53:982–992.[CrossRef][Medline]
- Sakai H, Abe K, Kobayashi Y, et al. Clinical guidelines of IgA nephropathy. Jpn J Nephrol (1995) 37:417–421.
- Gooley TA, Leisenring W, Crowley J, et al. Statist Med (1999) 18:695–706.[CrossRef]
- Lemeshow S, Teres D, Klar J, et al. Mortality probability models (MPM II) based on an international cohort of intensive care unit patients. JAMA (1993) 270:2478–2486.
[Abstract/Free Full Text] - Harden PN, Geddes C, Rowe PA, et al. Polymorphisms in angiotensin-converting-enzyme gene and progression of IgA nephropathy. Lancet (1995) 345:1540–1542.[CrossRef][Web of Science][Medline]
- Hunley TE, Julian BA, Phillips JA III, et al. Angiotensin converting enzyme gene polymorphism: potential silencer motif and impact on progression in IgA nephropathy. Kidney Int (1996) 49:571–577.[Web of Science][Medline]
- Pei Y, Scholey J, Thai K, et al. Association of angiotensinogen gene T235 variant with progression of immunoglobin A nephropathy in Caucasian patients. J Clin Invest (1997) 100:814–820.[Web of Science][Medline]
- Xia YF, Huang S, Li X, et al. A family-based association study of megsin A23167G polymorphism with susceptibility and progression of IgA nephropathy in a Chinese population. Clin Nephrol (2006) 65:153–159.[Web of Science][Medline]
- Schena FP, Cerullo G, Rossini M, et al. Increased risk of end-stage renal disease in familial IgA nephropathy. J Am Soc Nephrol (2002) 13:453–460.
[Abstract/Free Full Text] - Frimat L, Philippe C, Maghakian MN, et al. Polymorphism of angiotensin converting enzyme, angiotensinogen, and angiotensin II type 1 receptor genes and end-stage renal failure in IgA nephropathy: IGARAS—a study of 274 men. J Am Soc Nephrol (2000) 11:2062–2067.
[Abstract/Free Full Text] - Schena FP, DAltri C, Cerullo G, et al. ACE gene polymorphism and IgA nephropathy: an ethnically homogeneous study and a meta-analysis. Kidney Int (2001) 60:732–740.[CrossRef][Web of Science][Medline]
- Izzi C, Ravani P, Torres D, et al. IgA nephropathy: the presence of familial disease does not confer an increased risk for progression. Am J Kidney Dis (2006) 47:761–769.[CrossRef][Web of Science][Medline]
- Cattran DC, Greenwood C, Ritchie S. Long-term benefits of angiotensin-converting enzyme inhibitor therapy in patients with severe immunoglobulin a nephropathy: a comparison to patients receiving treatment with other antihypertensive agents and to patients receiving no therapy. Am J Kidney Dis (1994) 23:247–254.[Web of Science][Medline]
- Kanno Y, Okada H, Yamaji Y, et al. Angiotensin-converting-enzyme inhibitors slow renal decline in IgA nephropathy, independent of tubulointerstitial fibrosis at presentation. Q J Med (2005) 98:199–203.[Web of Science]
- Maschio G, Cagnoli L, Claroni F, et al. ACE inhibition reduces proteinuria in normotensive patients with IgA nephropathy: a multicentre, randomized, placebo-controlled study. Nephrol Dial Transplant (1994) 9:265–269.
[Abstract/Free Full Text] - Praga M, Gutierrez E, Gonzalez E, et al. Treatment of IgA nephropathy with ACE inhibitors: a randomized and controlled trial. J Am Soc Nephrol (2003) 14:1578–1583.
[Abstract/Free Full Text] - Li PK, Leung CB, Chow KM, et al. Hong Kong study using valsartan in IgA nephropathy (HKVIN): a double-blind, randomized, placebo-controlled study. Am J Kidney Dis (2006) 47:751–760.[CrossRef][Web of Science][Medline]
- Samuels JA, Strippoli GF, Craig JC, et al. Immunosuppressive agents for treating IgA nephropathy. Cochrane Database Syst Rev (2003) 4:CD003965.[Medline]
- Hotta O, Miyazaki M, Furuta T, et al. Tonsillectomy and steroid pulse therapy significantly impact on clinical remission in patients with IgA nephropathy. Am J Kidney Dis (2001) 38:736–743.[Web of Science][Medline]
- Xie Y, Nishi S, Ueno M, et al. The efficacy of tonsillectomy on long-term renal survival in patients with IgA nephropathy. Kidney Int (2003) 63:1861–1867.[CrossRef][Web of Science][Medline]
- Kawamura T, Ohta T, Ohno Y, et al. Significance of urinalysis for subsequent kidney and urinary tract disorders in mass screening of adults. Intern Med (1995) 34:475–480.[Web of Science][Medline]
Accepted in revised form: 15. 5.09
![]()
CiteULike
Connotea
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

