Daria Pašalić
Department of Medical Chemistry, Biochemistry and Clinical Chemistry
Zagreb University School of Medicine
Šalata ul 2.
10 000 Zagreb, Croatia
Phone +385 (1) 4590 205; +385 (1) 4566 940
E-mail: dariapasalic [at] gmail [dot] com

Useful links

Original papers:


Mirjana Fuček*1, Živka Dika2, Sandra Karanović2, Ivana Vuković Brinar2, Vedran Premužić2, Jelena Kos2, Ante Cvitković3, Maja Mišić4, Josip Samardžić4, Dunja Rogić1, Bojan Jelaković2. Reliability of CKD-EPI predictive equation in estimating chronic kidney disease prevalence in the Croatian endemic nephropathy area.
Biochemia Medica 2018;28(1):010701. 
1Department of Laboratory Diagnostics, University Hospital Center Zagreb, School of Medicine University of Zagreb, Zagreb, Croatia
2Department of Nephrology, Hypertension and Dialysis, University Hospital Center Zagreb, School of Medicine University of Zagreb, Zagreb, Croatia
3Institute for Public Health, Brodsko Posavska County, Slavonski Brod, Croatia
4Department of PathologyGeneral Hospital "Dr.Josip Benčević ", Slavonski Brod, Croatia
*Corresponding author: mfucek [at] kbc-zagreb [dot] hr




Introduction: Chronic kidney disease (CKD) is a significant public health problem and it is not possible to precisely predict its progression to terminal renal failure. According to current guidelines, CKD stages are classified based on the estimated glomerular filtration rate (eGFR) and albuminuria. Aims of this study were to determine the reliability of predictive equation in estimation of CKD prevalence in Croatian areas with endemic nephropathy (EN), compare the results with non-endemic areas, and to determine if the prevalence of CKD stages 3-5 was increased in subjects with EN.

Materials and methods: A total of 1573 inhabitants of the Croatian Posavina rural area from 6 endemic and 3 non-endemic villages were enrolled. Participants were classified according to the modified criteria of the World Health Organization for EN. Estimated GFR was calculated using Chronic Kidney Disease Epidemiology Collaboration equation (CKD-EPI).

Results: The results showed a very high CKD prevalence in the Croatian rural area (19%). CKD prevalence was significantly higher in EN then in non EN villages with the lowest eGFR value in diseased subgroup.

Conclusions: eGFR correlated significantly with the diagnosis of EN. Kidney function assessment using CKD-EPI predictive equation proved to be a good marker in differentiating the study subgroups, remained as one of the diagnostic criteria for EN.

Key words: chronic kidney disease; estimated glomerular filtration rate; predictive equation; endemic nephropathy


Received: February 27, 2017                                   Accepted:October 08, 2017             Online publication date: November 24, 2017                




Over the past three decades the prevalence of chronic kidney disease (CKD) is on the rise worldwide (1). The Kidney Disease Outcomes Quality Initiative (KDOQI) consequently published guidelines for diagnosis of CKD and postulated the five stages for disease classification (CKD1-5), based on disease severity. Glomerular filtration rate (GFR) was set as the most important criteria for disease severity determination according to the 2002 guidelines (2).

In order to better assess disease progression, in 2012 the Kidney Disease: Improving Global Outcomes (KDIGO) group revised the CKD staging system and issued the new KDIGO CKD guidelines (3). This new disease staging system was based on five eGFR stages (the third stage (G3) was furthermore stratified into substages G3a and G3b) and three stages of albuminuria determined on the basis of the albumin-to-creatinine ratio (ACR). However the definition of CKD did not change, i.e. CKD was defined as an eGFR < 60 mL/min/1.73m2 during the period of 3 months (4). Based on these indicators and following the results of a meta-analysis that confirmed the independent value of estimated GFR (eGFR) and albuminuria for predicting mortality due to cardiovascular events, acute kidney injury, progressive CKD and end-stage renal disease (ESRD), a “heat map” was generated according to which patients could be classified into four risk categories (depending on renal and cardiovascular outcomes): low (stages G1-A1, G2-A1), moderate (stages G1-A2, G2-A2, G3a-A1), high (G1-A3, G2-A3, G3a-A2, G3b-A1) and very high risk (G3a-A3, G3b-A2-3, all G4 and G5) (5).

Epidemiological data on the prevalence of CKD in Croatia are still missing despite the increasing awareness of public health problems associated with CKD. There is only one preliminary report on prevalence of CKD in Croatian rural population which showed that CKD is more frequently present in this area comparing to other regions (6). One specific cause of CKD in Croatia is endemic nephropathy (EN). It was reported in Posavina and comprises 14 villages situated west of the city of Slavonski Brod (Figure 1) with a population of 12,686 inhabitants according to the census from 2011 (7).


Figure 1. Endemic focus in Croatia with respective epidemiological data. The average prevalence of endemic nephropathy in endemic villages from 1980 to 1991 is presented as percentage (20). EN-villages enrolled in this study are rounded with a full black line; other EN villages (not included in this study ) are rounded with dashed black line; non-EN villages are rounded with a dotted line.


EN is a chronic tubulointerstitial nephropathy with subtle onset and characteristically gradual progression to ESRD with no significant gender difference  - some reports in which EN is more frequent in women than in men are probably incorrect because before developing a kidney disease many men die of some cardiovascular disease (8).

The cause of EN is chronic aristolochic acid poisoning through contaminated food in genetically predisposed individuals (9). The adducts of aristolochic acid were isolated from the renal cortex of EN patients only, and not from patients with CKD of different etiology. Furthermore, a fingerprint mutation of the p53 gene (an AT:TA transversion) was confirmed in EN patients (9). Today, EN is considered as an environmental form of aristolochic acid nephropathy. Clinical course and disease progression of EN do not differ from other tubulointerstitial nephropathies (10). For now, there are no specific biomarkers for EN diagnosis and diagnostic criteria have not been precisely established. This means that various countries used different criteria. However, all these criteria have serious disadvantages (outdated, uneven and not in accordance with the new recommendations for the classification of renal diseases). Recently, a consensus document with new diagnostic criteria was accepted and it is hoped that researcher from all countries will start to use these compilation of criteria (11).

In EN, GFR is, just as in other CKD, considered a marker of overall renal function (12). Furthermore, individual patients’ risk estimation, prognosis of the disease clinical course and treatment decisions are based upon patients’ GFR. At the population level, it is important to accurately determine the prevalence of CKD to enable patient care planning and to estimate the need for renal replacement therapy as this is essential for planning public health budget, especially in areas with high incidence of CKD such as the Croatian areas with EN.

Because of higher exposure to aristolochic acid in the past, prevalence of CKD is higher in Croatian endemic area than in non-endemic villages. Our aims were to determine the prevalence of CKD using Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation and to analyse diagnostic value of eGFR in diagnosis and stratification of subjects from EN area. The results obtained contribute to the detection of differences in the prevalence of CKD stages and in their frequency in endemic and non-endemic villages.


Materials and methods



This cross-sectional observational study is based on data collected by field research between 2008 and 2010 as a part of the Croatian Ministry of Science research project entitled “Endemic Nephropathy in Croatia: Epidemiology, Diagnosis, Etiology” (108-0000000329). The survey was approved by the Ethical Committee of the Croatian National Institute of Public Health, the School of Medicine University of Zagreb and General Hospital ‘Dr. Josip Benčević’ Slavonski Brod.

A total of 1573 adult inhabitants from six endemic villages (Pričac, Slavonski Kobaš, Živike, Šumeće, Zbjeg and Dobočac) and three non-endemic villages (Klakar, Donja Bebrina, Rastušje) were invited to participate on a door-to-door basis and were enrolled after signing informed consent. The study included significant and comparable percentages of participants from each village over previous researches (Table 1) (13).


Table 1. The percentage of the population enrolled in the study


A physician, nurse and technician went to the village where an extensive epidemiological questionnaire was completed and clinical examination (including measurement of blood pressure, height, weight) and sampling for laboratory analyses of blood and urine were performed. Venous blood (4 mL in K2EDTA blood collection tube for red blood cell count and 8.5 mL in tube with Silica Clot Activator-SST for biochemistry tests; all BD Diagnostic, Sparks, USA) and 50 mL of second morning urine sample (Sterile Specimen Collection Cup of 120 mL without preservative; BD Diagnostic, Sparks, USA) were collected from each subject. Blood samples were centrifuged 10 minutes at 3500 rpm within 2 hours of collection. All samples were transported the same day (in conditions that ensure the stability of the analyte) to the Department of Laboratory Diagnostics, University Hospital Center Zagreb, and testing was performed immediately upon acceptance.

According to the modified World Health Organization (WHO) criteria, which have been used in Croatia for the last 40 years, participants in endemic villages were classified into diseased, suspected of having EN, those at risk of developing EN and others, farmers from villages who were completely unrelated to EN (14). Subjects from endemic villages with missing key data for classification were not included in the study. Inclusion criteria were based on the following data: (a) positive family history of EN and/or living > 20 years in EN village; (b) low molecular weight proteinuria (alpha1-microglobulin > 10 mg/L or alpha1-microglobulin/creatinine > 14 mg/g); (c) serum creatinine > 132.6 µmol/L; (d) anaemia (haemoglobin (Hb) < 120 g/L if male, Hb < 113 g/L if female), and (e) exclusion of other renal diseases (including diabetes).

Participants were considered “diseased” if they were positive for criteria listed as (a), (b), (c), (d) and (e); or (b), (c), (d) and (e); and/or (a), (b), (d) and (e). Patients with positive criteria listed as (a) and (b) or (b) and (d) were “suspected of having EN” and “at risk” if they were from a family with EN. To exclude other renal diseases, laboratory data and imagining techniques were used when necessary. The control group was comprised of subjects who did not have EN according to the previously listed criteria (subjects from non-endemic villages).

Due to the fact that diabetes and hypertension are major causes of CKD, participants were additionally divided into subjects with hypertension according to the European Society of Hypertension classification (i.e. systolic blood pressure of ≥ 140 mmHg and diastolic blood pressure of ≥ 90 mmHg, measured with an automated device (Omron HEM-907, Matsusaka, Japan); and/or taking any antihypertensive agent) and subjects with diabetes (DM) according to the American Diabetes Association guidelines (fasting blood glucose ≥ 7.0 mmol/L measured by enzymatic UV method with hexokinase and/or taking drugs to lower blood glucose) (15,16). Furthermore, participants were divided into three age groups: group I (age 18 - 42); group II (age 43 - 66), and group III (age ≥ 67).



Blood and urine samples were analyzed using standard laboratory methods. Red blood cell count (RBC) was determined using laser light scattering technology on the Cell Dyn Sapphire (Abbott Diagnostic, Illinois, USA) and Sysmex XE 5000 (Sysmex Europe GmbH, Norderstedt, Germany) analyzers. Serum (SCr) and urine creatinine (UCr) were measured on the Olympus AU 2700 analyzer (Beckman Coulter, California, USA) using the Jaffé kinetic uncompensated method with continuous measurement. Calibration was performed using proprietary calibrators traceable to isotope dilution mass spectrometry (IDMS) method and Standard Reference Material (SRM) 909B of the National Institute of Standards and Technology (NIST, USA) for serum and the NIST reference material SRM 914a (substance creatinine purity of 99.7 ± 0.3%) for urine samples. Urine alpha1-microglobulin (alpha1-MG) was measured by immunonephelometric assay (on the BN II nephelometer, Siemens Healthcare Diagnostics Inc, Tarrytown, USA). Continuous internal quality control was performed throughout the study using quality control materials provided by the respective manufacturers. For CKD staging, urine albumin (U-Alb) was measured in the second morning sample by immunonephelometric assay (BN II nephelometer, Siemens Healthcare Diagnostics Inc, Tarrytown, USA) standardized using primary ERM-DA470 calibrators with a method sensitivity of 3.0 mg/L. The same tests were performed in the endemic and the control group.

Glomerular filtration rate (eGFR in mL/min/1.73m2) was calculated by applying SCr to CKD-EPI equation (17). The equations are presented below (Eq. 1 for women and Eq. 2 for men).

(Eq. 1)

if SCr £ 62 μmol/L
2) = 144 x (SCr (μmol/L) / 62) - 0.329 x (0.993)age (years)

if SCr > 62 μmol/L
2) = 144 x (SCr (μmol/L) / 62) - 1.209 x (0.993)age (years)

(Eq. 2)

if SCr £ 80 μmol/L
GFR (mL/min/1.73m
2) = 141 x (SCr (μmol/L)/80) - 0.411 x (0.993)age (years)

if SCr > 80 μmol/L
GFR (mL/min/1.73m
2) = 141 x (SCr (μmol/L) / 80) – 1.209 x (0.993)age (years)

Because the CKD-EPI equation was evaluated according to the SCr values measured by compensated Jaffé method, SCr values of all subjects were corrected before inclusion in the equation using the linear regression equation: y = 0.95 x + 0.44 (y = compensated Jaffè; x = uncompensated Jaffè) (18).

The revised KDIGO classification was used to define CKD stages which includes five stages of eGFR: G1 (≥ 90 mL/min/1.73m2), G2 (60 - 89 mL/min/1.73m2), G3a (45 - 59 mL/min/1.73m2), G3b (30 - 44 mL/min/1.73m2), G4 (15 - 29 mL/min/1.73m2), G5 (< 15 mL/min/1.73m2) and three levels of albuminuria based on the albumin-to-creatinine ratio (ACR): A1 (< 30 mg/g), A2 (30 – 300 mg/g) and A3 (> 300 mg/g).

CKD prevalence was assessed using the eGFR calculated with CKD-EPI equation. We evaluated the overall prevalence of CKD and the prevalence of stages 3-5 (CKD3-5) in all subgroups depending on gender, age and place of recruitment. CKD was defined as eGFR < 60mL/min/1.73m2.


Statistical analysis

The SPSS software (SPSS Version 18, SPSS Inc., Chicago, USA)   was used for statistical analyses, and P-value below 0.05 was considered statistically significant. Normality of data distribution was tested with Shapiro-Wilks test. Student t-test, analysis of variance and χ2 post hoc test were performed to

test the differences between groups. In case of normal distribution Mann-Whitney test was used. The correlation between eGFR and age subgroups was

analyzed using Pearson correlation. Differences between three or more groups were analyzed in the case of normally distributed variables with ANOVA

analysis of variance and Tukey-Kramer post hoc test.




A total of 1573 participants were enrolled, 1229 from EN villages and 344 from non-EN villages. In the group of endemic villages, 33 farmers were classified as having EN (2.7%), 86 as suspect of having EN (6.9%), 268 as those at risk (22%) and 842 were completely unrelated to EN – so called “others” (68%). We failed to find difference in gender prevalence within EN- and control villages. The average age of all participants was 52 years (range 18 - 90), with statistically significant difference in age between males and females in EN villages. In terms of body weight and SCr, there were statistically significant differences between males and females regardless of their place of recruitment, but also within EN- and control villages. No statistically significant differences regarding diabetes and hypertension were found between subjects from EN and control villages (Table 2).


Table 2. Basic participant characteristics in EN and non-EN villages


Prevalence of CKD in the Croatian rural area using CKD-EPI equation was very high (19%). Also, CKD prevalence was higher in females than in males (21% vs. 18%, P < 0.001) and significantly higher in EN then in non-EN villages in both sexes (P < 0.001) (Table 3).


Table 3. CKD prevalence and distribution of CKD stages in EN and non-EN villages


No differences were observed in the prevalence of CKD in hypertension and diabetes subgroups of participants between EN and non-EN villages (P = 0.704 and P = 0.945, respectively). The prevalence of CKD rose with patient’s age regardless of used equation. In EN diseased, CKD prevalence was the highest (88%).

There was no statistically significant difference of eGFR between age groups in EN and non-EN villages. Compared according to gender, significantly higher eGFR were observed in males than females. This difference was more pronounced in subjects from EN villages (Table 4).


Table 4. eGFR according to CKD-EPI equation in EN and non-EN villages


The lowest eGFR value was expected in diseased (eGFR in all subgroups are presented in Table 5).


Table 5. eGFR in EN and control subgroups according to the CKD-EPI equation


There was no statistically significant difference between subjects with or without diabetes in eGFR. Between subjects with or without hypertension, no statistically significant difference was found in eGFR, with lower values in hypertensive subjects.




Diagnosis and staging of CKD depends on the correct calculation of GFR. We tested differences in eGFR among subjects from EN and non-EN areas by applying serum creatinine values to recommended CKD-EPI predictive equation.

CKD-EPI equation differentiated healthy from diseased subjects well. Participants from EN villages had higher prevalence of CKD than participants from non-EN villages. There are no new cases of EN today, and participants who were classified into diseased subgroup thirty years ago, today have more CKD3-5 stages than others in EN villages. No statistically significant difference was observed in the eGFR between patients with or without diabetes, although CKD-EPI equation underestimated their GFR (rather consistently with literature data) (19). In the hypertension group, significantly more G3a, G3b, G4 and G5 CKD stages were observed (20% vs. 7.7% in the non hypertension group). High prevalence of CKD3-5 stages in that subgroup could be a result of subjects lifestyle – high prevalence of obesity and hypertension, poor control of salt intake and untreated hypertension.

We are aware that this study has several limitations: recommended methods for serum creatinine measurements were not used in the study (enzymatic method has been applied in the Department of Laboratory Diagnostics, University Hospital Center Zagreb since 2014 and study is based on data collected between 2008 and 2010) and creatinine was determined only once so we do not know if there was a clinically significant difference in re-determining (based on the analytical and biological variation of serum creatinine the clinically significant difference for eGFR is about 14%). Also, eGFR in our study group, were not compared with the GFR measured using the gold standard  with an ideal filtration marker.

However, the value of our study is that we found differences in the CKD prevalence and the frequency of certain stages of CKD in EN and non-EN villages. The prevalence of EN diseased (2.7%) and suspected (6.9%) patients is in accordance with the previously described incidence decline trend. This percentage of diseased and suspected patients explain the difference in CKD prevalence and it is an important indicator that in the coming years, regardless of if there will be new cases of EN, the number of new patients to start with renal replacement therapy will be higher in EN villages than in control. These results confirm that the assessment of kidney function with predictive equations proved to be a good marker in differentiating the study subgroups in the EN areas and it should therefore remain the diagnostic criteria for EN, and that the optimal equation for GFR estimation is CKD-EPI equation.




The study was part of the Croatian Ministry of Science research project entitled “Endemic Nephropathy in Croatia: Epidemiology, Diagnosis, Etiology” (108-0000000329).

Potential conflict of interest

None declared.



1.James MT, Hemmelgarn BR, Tonelli M. Early recognition and prevention of chronic kidney disease. Lancet 2010;375:1296–309.

 2. National Kidney Foundation. K/DOQI Clinical Practice Guidelines for Chronic Kidney Disease: Evaluation, Classification and Stratification. Am J Kidney Dis 2002;39:S1-S266.

 3. Kidney Disease: Improving Global Outcomes (KDIGO) CKD Work Group. KDIGO 2012 Clinical Practice Guideline for the Evaluation and Management of Chronic Kidney Disease. Kidney Inter 2013;Suppl.3:1–150.

 4. Levey AS, de Jong PE, Coresh J, El Nahas M, Astor BC, Matsushita K, et al. The definition, classification and prognosis of chronic kidney disease: a KDIGO controversies conference report. Kidney Int 2011;80:17-28.

 5. Matsushita K, van der Velde M, Astor BC, Woodward M, Levey AS, de Jong PE, et al. Association of estimated glomerular filtration rate and albuminuria with all-cause and cardiovascular mortality in general population cohorts: a collaborative meta-analysis. Lancet 2010;375:2073-81.

 6. Premužić V, Vrdoljak A, Fuček M, Karanović S, Vuković-Lela I, Kos J, et al. Surprisingly high prevalence of cronic kidney disease not only in Croatian endemic nephropathy area but also in non endemic villages. 50th ERA-EDTA CONGRES Nephrol Dial Transplant 2013;28:i140-i54.

 7. Croatian bureau of statistics. Census of Population, Households and Dwellings 2011. Available at: . Accessed May 22th 2017. (in Croatian)

 8. Jelaković B, Dika Ž, Karanović S, Vuković LI. Endemska nefropatija u Hrvatskoj. Liječ Vjesn 2015;137:100–8. (in Croatian)

 9. Jadot I, Declèves AE, Nortier J, Caron N. An Integrated view of aristolochic acid nephropathy: Update of the Literature. Int J Mol Sci 2017;18:297.

10. Bamias G, Boletis J. Balkan nephropathy: evolution of our knowledge. Am J Kidney Dis 2008;52:606-16.

11. Jelaković B, Nikolić J, Radovanović Z, Nortier J, Cosyns JP, Grollman PA, et al. Consensus statement on screening, diagnosis, classification and treatment of endemic (Balkan) nephropathy. Nephrol Dial Transplant 2014;29:2020–7.

12. Gifford FJ, Gifford RM, Eddleston M, Dhaun N. Endemic nephropathy around the world. Kidney Int Rep 2017;2:282-92.

13. Čeović S, Miletić-Medved M. Epidemiological features of endemic nephropathy in focal area of Brodska Posavina, Croatia. In: Čvorišćec S, Čeović S, Stavljenić-Rukavina A, eds. Endemic nephropathy in Croatia. Zagreb: Academia Croatica scientiarum medicarum;1996. p. 7-21.

14. World Health Organization Memorandum. The endemic nephropathy of South-Eastern Europe. Bull World Health Organ 1965;32:441–8.

15. Mancia G, De Backer G, Dominiczak A, Cifkova R, Fagard R, Germano G,et al. 2007 Guidelines for the management of arterial hypertension: The Task Force for the Management of Arterial Hypertension of the European Society of Hypertension (ESH) and of the European Society of Cardiology (ESC). Eur Heart J 2007;28:1462-536.

16. American Diabetes Association (ADA). Diagnosis and classification of diabetes mellitus. Diabetes Care 2008;31:S55-60.

17. Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF 3rd, Feldman HI, et al. A new equation to estimate glomerular filration rate. Ann Intern Med 2009;150:604-12.

18. Wuyts B, Bernard D, Van Den Noortgate N, Van De Walle J, Van Vlem B, De Smet R, et al. Reevaluation of formulas for predicting creatinine clearance in adults and children using compensated creatinine methods. Clin Chem 2003;49:1011-4.

19. Rigalleau V, Lasseur C, Raffaitin C, Perlemoine C, Barthe N, Chauveau P, et al. The Mayo Clinic quadratic equation improves the prediction of glomerular filtration rate in diabetic subjects. Nephrol Dial Transplant 2007;22:813-8.

20. Miletić-Medved M, Domijan AM, Peraica M. Recent data on endemic nephropathy and related urothelial tumors in Croatia. Wien Klin Wochenschr 2005;117:604-9.