- Open Access
Heterogeneity of clinical indices among the older dialysis population—a study on Japanese dialysis population
Renal Replacement Therapyvolume 3, Article number: 1 (2017)
The older dialysis population is increasing, and several studies have investigated the effect of actual age on the relationship between clinical parameters or clinical practice and prognosis. There is concern that biological age does not necessarily reflect the actual effects of aging, and it is important to take heterogeneity into account when considering the older population. In this study, we investigated whether such heterogeneities exist in a large number of clinical parameters and how any disparities might be affected by the domains of clinical indices.
We investigated the published, summarized results of the Japanese Society for Dialysis Therapy Renal Data Registry, a nationwide survey of the Japanese dialysis population, for the year 2013.
The potential population comprised 306,925 dialysis patients. We investigated coefficient of variations (CVs) of 25 clinical indices across age groups and compared them with the CVs of the 45–59 years age group as a reference group, using tests for homogeneity of variances. Almost all variables showed some heterogeneity by age. The domain for muscle mass or visceral protein exhibited greater heterogeneity in older patients, but surprisingly, the domain for body mass or physique exhibited less heterogeneity. The mean values of most parameters declined in the older population.
This study demonstrated that variations exist in clinical indices in the older dialysis population. Apart from biological age, we should take these heterogeneities into account when interpreting the findings of clinical studies involving older dialysis patients and in their individualized management.
The dialysis population has become much older in recent years, and the proportion of elderly patients is increasing in many registries [1, 2]. There are several clinical problems to overcome in the older hemodialysis population. One of the most important involves malnutrition-related syndromes, such as sarcopenia , protein energy wasting , and frailty . These syndromes are in fact associated with a worse prognosis in the dialysis population as a whole [6, 7]; these phenotypes are considered to be most important in the management of the older population.
Strategies against complications for older dialysis patients sometimes differ from those for their younger counterparts . However, the definition of “older” patients varies in different reports [9–11]. The reason for this may be due to the heterogeneity among the elderly population in terms of phenotype. Some patients look younger than their actual age, while some look older. Moss et al. reported that the “surprise” question, “Would I be surprised if this patient died in the next year?” correlated well with 1 year survival , indicating that the phenotype, presumably wasting or malnutrition, seems important for survival.
We speculated that such phenotype of the “older” population does matter and that it is important to take the heterogeneities of such patients into account, which might affect the older phenotype. However, there have been no reports on the heterogeneity of the elderly population. Therefore, here, we focused on the heterogeneity of the older dialysis population and investigated which clinical indices exhibited large heterogeneity and were potentially related to a poor prognosis, using national registry data compiled by the Japanese Society for Dialysis Therapy (JSDT). It is anticipated that our results will help the clinician to undertake individualized management of older dialysis patients.
Data from the annual Japanese Society for Dialysis Therapy Renal Data Registry (JRDR) survey run by JSDT were used for this analysis. The details of this registry have been described previously . We analyzed the data collected in 2013. The summarized tables are recorded in CD-ROM or are reported on the internet. In this analysis, we used data that appears on JSDT’s website .
The potential population for this study was the entire dialysis population as of December 31, 2013, namely 306,925 patients. The population was divided by 15-year age groups. We excluded patients with missing data on age (n = 12) and patients who were younger than 15 years (n = 116).
The following clinical indices were investigated: (1) mineral and bone disorder markers (corrected calcium and phosphate); (2) malnutrition-related factors (total cholesterol, high-density lipoprotein (HDL) cholesterol, non-HDL cholesterol, percent creatinine generation rate (%CGR) , normalized protein catabolic rate (nPCR) , and albumin); (3) body mass and volume status (body height, postdialytic body weight (BW), body mass index (BMI), BW change, and fractional BW change); (4) dialysis prescription and its consequences (dialysis duration, spKt/V , eqKt/V , blood flow rate, predialytic blood urea, and creatinine); (5) circulatory factors (systolic blood pressure, diastolic blood pressure, and pulse rate); and (6) other factors (hemoglobin (Hb) level, glycoalbumin, and HbA1c). Postdialytic BW, dialysis duration, and blood flow rate were investigated among the patients who were treated by hemodialysis, hemodiafiltration, or hemoperfusion. The values of%CGR, nPCR, BW change, fractional BW change, spKt/V, eqKt/V, and creatinine were investigated only in patients who were treated by hemodialysis, hemodiafiltration, or hemoperfusion three times weekly with a dialysis history of ≥2 years. Glycoalbumin and HbA1c were investigated among those who had history of diabetes or were on glucose-lowering therapies. Parameters other than those described above were studied in the entire dialysis population.
This analysis consists of two parts: one is to determine the heterogeneities of age populations, and the other is to compare them with the standard of 45–59 year-old group.
The heterogeneity was measured by coefficient of variations (CV) in this analysis, because standard deviation (SD) does not necessarily express the degree of heterogeneity when the averages vary. The tables of the JRDR summarized data give the number of patients and mean values and SD for each parameter in all age groups. CVs were calculated by dividing the SD by the mean for each age category. The larger CV is more heterogenic or heteroscedastic the population can be considered.
Relative CVs were calculated by dividing the CVs of each age group by that of the 45–59 years age group. If the value is larger than one, the index is more heterogenic in older population than the younger counterpart. Relative means were also calculated as the ratio of the mean of each age group to that of the 45–59 years age group, similarly to relative CVs.
Finally, we statistically compared the heterogeneities across age groups by investigating the homogeneity of variances. In details, we calculated variances by squaring standard deviations. The tests for homogeneity of variances were performed on the squared relative CVs using Bartlett’s test among overall groups and F tests for each pair of age groups. The former investigated the homogeneity of the variances by age groups as a whole. When the test indicated statistical significance, it shows that the clinical index is heterogenic across age groups. On the other hand, the latter tested statistical significance by comparing the CV of each age group with that of the reference group of ages 45–59. When a test was significant, it indicates that the clinical index is more heterogenic in that age group compared to that in the 45–59-year-old age group. Post hoc adjustment for multiple comparisons was done by Bonferroni correction. The p values less than 0.05 were considered significant. For multiple comparisons, p values less than 0.01 were considered significant because five pairs of comparisons were performed with the 45–59 years age group as the reference group. All statistical calculations and analyses were performed using Microsoft Excel 2013 (Microsoft Corporation, Redmond, WA).
This study is the secondary analyses of the summarized and published data and does not deal with any personal information. Therefore, we consider that this study is deemed exempt from the Institutional Review Board or Ethics Committee approval. This study adhered to the Declaration of Helsinki.
The demographics and characteristics of the total population who underwent dialysis therapy in 2013 are shown in Table 1 . Mean age was 67.21 ± 12.51 years (mean ± SD). The age groups 15–29, 30–44, 45–59, 60–74, 75–89, and ≥90 years accounted for 0.37, 4.91, 18.56, 45.87, 28.72, and 1.54% of the total population, respectively.
Table 2 shows the number of patients and their proportion of sex, vintage, and primary diagnoses for each age category. In the older age groups, there was a higher proportion of female patients and a higher proportion of patients with the primary diagnosis of nephrosclerosis; less predominant were the proportion of patients with a long history of dialysis and the proportion with diabetes as the primary diagnosis.
Table 3 shows the SDs and means of the variables investigated as well as the number of the patients. Table 4 shows the CVs relative to those of the 45–59 years age group, which is used as the reference group. The overall differences of variances, including chi-square values and p values calculated with the Bartlett test, are shown, as well as p values for the differences of variances of the age groups, taking the variance of the 45–59 years age group as the references. For almost all variables, the variances differed by age group. However, the trends themselves differed by index. We can categorize these indices according to the relative CV for the 75–89 years age group into four categories: (1) relative CVs of more than 1.1 (indicating the most heteroscedastic indices), namely for %CGR, albumin, creatinine, fractional BW reduction, blood urea nitrogen, and diastolic blood pressure; (2) relative CVs of 1.0–1.1 (mildly heterogeneous), namely for phosphate, total cholesterol, nPCR, body height, BW reduction, blood flow rate, systolic blood pressure, pulse rate, and Hb; (3) relative CVs of 0.9–1.0 (less heterogeneous), namely for corrected calcium, non-HDL, HDL, session duration, spKt/V, and eqKt/V; and (4) relative CVs of less than 0.9 (least heterogeneous), namely for BW, BMI, glycoalbumin, and HbA1c.
Figure 1 indicates that the disparities in representative indices between several domains differ significantly across the age groups. The relative CVs of indicators relating to the muscle mass or visceral proteins increased steadily with age, becoming more heteroscedastic. The parameters related to dietary intake showed a modest relationship across the age groups. On the other hand, indices of body mass demonstrated that the dispersion progressively decreased with the age of the population.
Figure 1 was designated to demonstrate the CVs by themselves so that the heteroscedasticity can be recognized easily. On the other hand, the CVs were calculated from the SDs and means. The means are also subjected to change by age groups, and the relative means indicate the trend of the distributions of the absolute values as a whole by age groups. Therefore, in Fig. 2, we added the information about the relative means to the results shown in Fig. 1 so that information about whether the absolute values tended to become smaller among the older population as well as the information about the heteroscedasticity. This figure indicates that the mean values of all indices included in Fig. 1 were lower in the older age groups, irrespective of relative CVs. This fact might indicate that the normal aging process accompanies universal declines in clinical indices.
Heterogeneity as determined by CV increased with age in some domains, especially that of muscle mass or visceral proteins. However, in some other domains such as body mass, the disparities seem to be diminished in the older population. This is the first study investigating the heteroscedasticity of clinical indices, using nationwide registry data.
There were differences among the domains of indices; factors related to sarcopenia or wasting were most prominently dispersed. Moreover, the mean values of these parameters declined with age. Therefore, muscle mass and albumin level generally decreases with age, but their speed of decline varies among individuals. The decrease in these parameters, as discussed above, reportedly relates to poor survival relating to wasting [6, 7]. Therefore, the heterogeneity of these clinical indices might reflect the diversity of the older population in terms of “old” phenotypes and relate to a worse prognosis. Moreover, the mean values of CGR and albumin decreased less in older population, compared to other indices as shown in Fig. 2. This fact might reflect that CGR and albumin are less affected by aging and they could be good markers of aging phenotype. In this analysis, we investigated the disparity in muscle mass in terms of CGR. However, muscle strength is reported to relate to survival to a greater degree . The concept of dynapenia has been proposed [18, 19]. Future studies investigating the heterogeneity of muscular strength and its relation to survival are required.
Albumin has been used as a marker of malnutrition, but many other factors also affect albumin levels . Dietary restriction only does not decrease the albumin level and it is now regarded as a factor of wasting . The 2013 annual JRDR survey, the source of the data analyzed in this study, also investigated C-reactive protein (CRP). However, the distribution of CRP does not follow a normal distribution, and its heterogeneity cannot be assessed by SD or CV as we performed on other variables because the variables that do not follow normal distribution cannot be summarized by SDs or means. Therefore, we could not investigate the heterogeneity of CPR, and we cannot draw the conclusion that the disparity of albumin relates to inflammation status.
Dietary intake is reduced in the elderly population because of reduced appetite . Moreover, the degree of appetite itself relates to prognosis . In this study, we used nPCR as a surrogate of dietary intake. On the other hand, phosphate is contained in foods that contain much protein . Phosphate binders and dialysis dosage, especially session length, also affect the level of phosphate . It is interesting to note then that nPCR and phosphate exhibited similar trends. Intradialytic BW reduction can be related to sodium intake during interdialytic periods , which also followed a similar trend to that of nPCR and phosphate levels. Therefore, these clinical indices can be considered as surrogate markers of dietary intake in patients.
Surprisingly, the distribution of body mass tended to be smaller in the older population even though the disparities of muscle mass become wide among the elderly population. An inverse relationship can be observed in the dialysis population ; patients with larger BMI experience better survival . However, in the Asian population, such relationship might be weaker than in the white or black populations . Female patients tend to have better survival among the Japanese population, and the average age of female patients is higher than their male counterparts . Actually, Table 2 shows that the proportion of female patients was larger among the older age groups. Obviously, female patients have a smaller physique than male patients. Therefore, the present results might be attributed to the predominance of female patients in the older age groups. On the other hand, it is possible that the findings of the smaller and relatively homogeneous physique in older patients merely reflect the trend in the general population because this study did not examine variances in the general population.
Moreover, there was a discrepancy in terms of heteroscedasticity between the indices regarding body mass such as BMI or BW and the indices regarding muscle mass such as %CGR or creatinine. The former showed to be less heterogenic, while the latter showed to be more heterogenic in older population. BMI is one of the simplest indices to describe the physiques. However, it cannot distinguish fat mass from muscle mass. The dialysis patients are shown to reduce muscle mass, while they often lose less amount of fat or even gain it during the course of dialysis therapy . Therefore, BMI can underestimate the changes in muscle mass . This point might cause the discrepancy of heteroscedasticity between body mass and muscle mass. On the other hand, body height can contribute to BMI. However, as Table 4 shows, the heterogeneity of height remained small even in the older population. Moreover, the mean values of body height were less in older population. Several recent studies demonstrated that taller dialysis patients might experience worse survival [28, 29]. Thus, the smaller and less heterogenic older population in terms of their body height might be due to the survival effects partly.
Nonetheless, we should pay attention to the clinical indices which exhibited large heterogeneity, such as creatinine generation rate, creatinine, and albumin, when we perform clinical practice on the older population. The deterioration of these indices might reflect the degree of the morbid conditions among older population. Future study will be required to investigate the relationship to the clinical outcomes, such as mortality, quality of life (QOL), and activity of daily living (ADL). After that, the interventions to the indices are warranted to improve the outcomes of the patients.
This study has several limitations. The first is its cross-sectional and observational nature. Therefore, we cannot draw any conclusions about the time course or cause-effect relationship. The second is that this study was based on results that appeared in summarized tables, and therefore, detailed investigations, such as further stratification by sex or primary diagnoses, could not be performed. The third is that the data that we used did not contain the data for clinical outcomes. Therefore, we could not investigate the effect of heterogeneity observed on mortality, QOL, or ADL. The fourth is that the JRDR data were almost exclusively derived from Japanese people. Therefore, the applicability of our findings to other ethnicities is unclear. Lastly, data was missing on some parameters in some patients. Although the number of the patients was sufficiently large, which is the most important advantage of this study, it is possible that errors exist.
Although there are several limitations in this study, the strengths of this study are the large number of patients and parameters that were examined. These results may contribute to the individualized management of the older dialysis population. Further studies on the cause of disparities and the relationship to the prognosis of the patients are required. We anticipate that the findings of such studies will provide a better understanding of phenotypes in the older dialysis population.
This is the first study to investigate heterogeneity in the older dialysis population. There were some disparities in heterogeneity across the parameters. Clinical parameters related to muscle mass or wasting exhibited larger heterogeneity, while those related to body mass or BW exhibited less heterogeneity. Individualized management for older dialysis patients with these heterogeneities might contribute to better outcomes in this population.
Percent creatinine generation rate
Activity of daily living
Body mass index
Compact disk read-only memory
Coefficient of variation
The Japanese Society for Dialysis Therapy Renal Data Registry
The Japanese Society for Dialysis Therapy
Normalized protein catabolic rate
Quality of life
Single pool Kt/V
Masakane I, Nakai S, Ogata S, Kimata N, Hanafusa N, Hamano T, Wakai K, Wada A, Nitta K. An overview of regular dialysis treatment in Japan (as of 31 December 2013). Ther Apher Dial. 2015;19:540–74.
Saran R, Li Y, Robinson B, Ayanian J, Balkrishnan R, Bragg-Gresham J, Chen JT, Cope E, Gipson D, He K, et al. US renal data system 2014 annual data report: epidemiology of kidney disease in the United States. Am J Kidney Dis. 2015;66:Svii. S1-305.
Fielding RA, Vellas B, Evans WJ, Bhasin S, Morley JE, Newman AB, Abellan van Kan G, Andrieu S, Bauer J, Breuille D, et al. Sarcopenia: an undiagnosed condition in older adults. Current consensus definition: prevalence, etiology, and consequences. International working group on sarcopenia. J Am Med Dir Assoc. 2011;12:249–56.
Fouque D, Kalantar-Zadeh K, Kopple J, Cano N, Chauveau P, Cuppari L, Franch H, Guarnieri G, Ikizler TA, Kaysen G, et al. A proposed nomenclature and diagnostic criteria for protein-energy wasting in acute and chronic kidney disease. Kidney Int. 2008;73:391–8.
Fried LP, Tangen CM, Walston J, Newman AB, Hirsch C, Gottdiener J, Seeman T, Tracy R, Kop WJ, Burke G, et al. Frailty in older adults: evidence for a phenotype. J Gerontol A Biol Sci Med Sci. 2001;56:M146–56.
Isoyama N, Qureshi AR, Avesani CM, Lindholm B, Barany P, Heimburger O, Cederholm T, Stenvinkel P, Carrero JJ. Comparative associations of muscle mass and muscle strength with mortality in dialysis patients. Clin J Am Soc Nephrol. 2014;9:1720–8.
Johansen KL, Chertow GM, Jin C, Kutner NG. Significance of frailty among dialysis patients. J Am Soc Nephrol. 2007;18:2960–7.
Hanafusa N, Nomura T, Hasegawa T, Nangaku M. Age and anemia management: relationship of hemoglobin levels with mortality might differ between elderly and nonelderly hemodialysis patients. Nephrol Dial Transplant. 2014;29:2316–26.
Locatelli F, Pisoni RL, Combe C, Bommer J, Andreucci VE, Piera L, Greenwood R, Feldman HI, Port FK, Held PJ. Anaemia in haemodialysis patients of five European countries: association with morbidity and mortality in the Dialysis Outcomes and Practice Patterns Study (DOPPS). Nephrol Dial Transplant. 2004;19:121–32.
Canaud B, Tong L, Tentori F, Akiba T, Karaboyas A, Gillespie B, Akizawa T, Pisoni RL, Bommer J, Port FK. Clinical practices and outcomes in elderly hemodialysis patients: results from the Dialysis Outcomes and Practice Patterns Study (DOPPS). Clin J Am Soc Nephrol. 2011;6:1651–62.
Foote C, Ninomiya T, Gallagher M, Perkovic V, Cass A, McDonald SP, Jardine M. Survival of elderly dialysis patients is predicted by both patient and practice characteristics. Nephrol Dial Transplant. 2012;27:3581–7.
Moss AH, Ganjoo J, Sharma S, Gansor J, Senft S, Weaner B, Dalton C, MacKay K, Pellegrino B, Anantharaman P, Schmidt R. Utility of the “surprise” question to identify dialysis patients with high mortality. Clin J Am Soc Nephrol. 2008;3:1379–84.
Hanafusa N, Nakai S, Iseki K, Tsubakihara Y. Japanese society for dialysis therapy renal data registry—a window through which we can view the details of Japanese dialysis population. Kidney Int Suppl. 2015;5:15–22.
An Overview of Regular Dialysis Treatment in Japan [http://member.jsdt.or.jp/member/contents/cdrom/2013/main.html]. Accessed 17 Apr 2016.
Shinzato T, Nakai S, Miwa M, Iwayama N, Takai I, Matsumoto Y, Morita H, Maeda K. New method to calculate creatinine generation rate using pre- and postdialysis creatinine concentrations. Artif Organs. 1997;21:864–72.
Shinzato T, Nakai S, Fujita Y, Takai I, Morita H, Nakane K, Maeda K. Determination of Kt/V and protein catabolic rate using pre- and postdialysis blood urea nitrogen concentrations. Nephron. 1994;67:280–90.
Daugirdas JT. Simplified equations for monitoring Kt/V, PCRn, eKt/V, and ePCRn. Adv Ren Replace Ther. 1995;2:295–304.
Clark BC, Manini TM. Sarcopenia =/= dynapenia. J Gerontol A Biol Sci Med Sci. 2008;63:829–34.
Kim JC, Kalantar-Zadeh K, Kopple JD. Frailty and protein-energy wasting in elderly patients with end stage kidney disease. J Am Soc Nephrol. 2013;24:337–51.
Friedman AN, Fadem SZ. Reassessment of albumin as a nutritional marker in kidney disease. J Am Soc Nephrol. 2010;21:223–30.
Lopes AA, Elder SJ, Ginsberg N, Andreucci VE, Cruz JM, Fukuhara S, Mapes DL, Saito A, Pisoni RL, Saran R, Port FK. Lack of appetite in haemodialysis patients—associations with patient characteristics, indicators of nutritional status and outcomes in the international DOPPS. Nephrol Dial Transplant. 2007;22:3538–46.
Kalantar-Zadeh K, Gutekunst L, Mehrotra R, Kovesdy CP, Bross R, Shinaberger CS, Noori N, Hirschberg R, Benner D, Nissenson AR, Kopple JD. Understanding sources of dietary phosphorus in the treatment of patients with chronic kidney disease. Clin J Am Soc Nephrol. 2010;5:519–30.
Levin NW, Kotanko P, Eckardt KU, Kasiske BL, Chazot C, Cheung AK, Redon J, Wheeler DC, Zoccali C, London GM. Blood pressure in chronic kidney disease stage 5D-report from a kidney disease: improving global outcomes controversies conference. Kidney Int. 2010;77:273–84.
Kalantar-Zadeh K, Block G, Humphreys MH, Kopple JD. Reverse epidemiology of cardiovascular risk factors in maintenance dialysis patients. Kidney Int. 2003;63:793–808.
Kalantar-Zadeh K, Kopple JD, Kilpatrick RD, McAllister CJ, Shinaberger CS, Gjertson DW, Greenland S. Association of morbid obesity and weight change over time with cardiovascular survival in hemodialysis population. Am J Kidney Dis. 2005;46:489–500.
Marcelli D, Brand K, Ponce P, Milkowski A, Marelli C, Ok E, Merello Godino JI, Gurevich K, Jirka T, Rosenberger J, et al. Longitudinal changes in body composition in patients after initiation of hemodialysis therapy: results from an international cohort. J Ren Nutr. 2016;26:72–80.
Abramowitz MK, Sharma D, Folkert VW. Hidden obesity in dialysis patients: clinical implications. Semin Dial. 2016;29:391–5.
Shapiro BB, Streja E, Ravel VA, Kalantar-Zadeh K, Kopple JD. Association of height with mortality in patients undergoing maintenance hemodialysis. Clin J Am Soc Nephrol. 2015;10:965–74.
Elsayed ME, Ferguson JP, Stack AG. Association of height with elevated mortality risk in ESRD: variation by race and gender. J Am Soc Nephrol. 2016;27:580–93.
The data reported here were provided by the JSDT. The interpretation and reporting of these data are the responsibility of the authors and do not reflect the official interpretation or views of the JSDT.
The authors did not have specific funding source to be disclosed.
Availability of data and materials
The data we used for the analysis have been already publicized by the JSDT on the web pages.
NH conceived of the study design, analyzed the data, and wrote the manuscript. SS and MN made the critical revisions to the manuscript and helped to finalize the article. All authors read and approved the final manuscript.
NH has received research consulting fees from Kyowa Hakko Kirin Co., Ltd.; lecture fees from Bayer Yakuhin, Ltd., Chugai Pharmaceutical Co., Ltd., Kyowa Hakko Kirin Co., Ltd., and Nikkiso Co., Ltd.; and grant support from Chugai Pharmaceutical Co., Ltd. and Kyowa Hakko Kirin Co., Ltd. NH was also a member of a division that is funded by Terumo Corporation. The other authors declare that they have no competing interests.
Consent for publication
This study does not include any individual person’s data. Therefore, this point is not applicable.
Ethics approval and consent to participate
Although this study employs the data originally derived from the patients, the data is widely available from the JSDT web pages and the data have already been anonymized. We the authors recognized that this study was not within the scope of the Ethical Guidelines for Medical and Health Research Involving Human Subject issued by the Ministry of Health, Labour and Welfare of Japan in December 2014, pursuant to the chapter 1 of part 3 of the guideline. Therefore, the ethics approval by the institutional review boards is exempted.