Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Primary kidney disease modifies the effect of comorbidities on kidney replacement therapy patients’ survival

  • Jaakko Helve ,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Jaakko.helve@helsinki.fi

    Affiliations Finnish Registry for Kidney Diseases, Helsinki, Finland, Abdominal Center Nephrology, University of Helsinki and Helsinki University Hospital, Helsinki, Finland

  • Mikko Haapio,

    Roles Conceptualization, Validation, Writing – review & editing

    Affiliation Abdominal Center Nephrology, University of Helsinki and Helsinki University Hospital, Helsinki, Finland

  • Per-Henrik Groop,

    Roles Conceptualization, Writing – review & editing

    Affiliations Finnish Registry for Kidney Diseases, Helsinki, Finland, Abdominal Center Nephrology, University of Helsinki and Helsinki University Hospital, Helsinki, Finland, Folkhälsan Institute of Genetics, Folkhälsan Research Center Biomedicum Helsinki, Helsinki, Finland, Department of Diabetes, Central Clinical School, Monash University, Melbourne, Australia

  • Patrik Finne

    Roles Conceptualization, Data curation, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing

    Affiliations Finnish Registry for Kidney Diseases, Helsinki, Finland, Abdominal Center Nephrology, University of Helsinki and Helsinki University Hospital, Helsinki, Finland

Abstract

Background

Comorbidities are associated with increased mortality among patients receiving long-term kidney replacement therapy (KRT). However, it is not known whether primary kidney disease modifies the effect of comorbidities on KRT patients’ survival.

Methods

An incident cohort of all patients (n = 8696) entering chronic KRT in Finland in 2000–2017 was followed until death or end of 2017. All data were obtained from the Finnish Registry for Kidney Diseases. Information on comorbidities (coronary artery disease, peripheral vascular disease, left ventricular hypertrophy, heart failure, cerebrovascular disease, malignancy, obesity, underweight, and hypertension) was collected at the start of KRT. The main outcome measure was relative risk of death according to comorbidities analyzed in six groups of primary kidney disease: type 2 diabetes, type 1 diabetes, glomerulonephritis (GN), polycystic kidney disease (PKD), nephrosclerosis, and other or unknown diagnoses. Kaplan-Meier estimates and Cox regression were used for survival analyses.

Results

In the multivariable model, heart failure increased the risk of death threefold among PKD and GN patients, whereas in patients with other kidney diagnoses the increased risk was less than twofold. Obesity was associated with worse survival only among GN patients. Presence of three or more comorbidities increased the age- and sex-adjusted relative risk of death 4.5-fold in GN and PKD patients, but the increase was only 2.5-fold in patients in other diagnosis groups.

Conclusions

Primary kidney disease should be considered when assessing the effect of comorbidities on survival of KRT patients as it varies significantly according to type of primary kidney disease.

Introduction

Many of the patients entering long-term kidney replacement therapy (KRT) have other medical conditions as well [1,2]. In addition, cause of kidney failure is associated with the frequency of comorbidities; for instance, KRT patients with diabetic nephropathy have more comorbidities than patients with other causes of kidney disease [3,4].

Several studies have shown that the effect of both comorbidities and of primary kidney disease is significant on KRT patients’ survival [5,6]. However, Van Manen et al. [7] reported that adding comorbidities to a multivariable model on top of primary kidney disease and general patient characteristics had only a small impact on hazard ratios of death. Nevertheless, all published models predicting risk of death for patients starting dialysis assume that the effect of comorbidities on survival is constant despite primary kidney disease [8].

Knowledge on whether primary kidney disease modifies the effect of comorbidities on survival is scarce or nonexistent. This study shows how comorbidities are associated with mortality according to primary kidney disease.

Materials and methods

All data were obtained from the Finnish Registry for Kidney Diseases. This registry has 97–99% coverage of all Finnish patients on chronic KRT since 1965. The registry is maintained by the Finnish Liver and Kidney Association and financed by the Finnish government. All patients aged 18 years or older who entered chronic KRT in an outpatient setting from 1 January 2000 to 31 December 2017 were included in the study (n = 8696). Information on age at start of KRT, sex, cause of kidney disease coded using ICD-10 codes, type of KRT (hemodialysis, peritoneal dialysis, or kidney transplantation), and time and cause of death were collected. Data on comorbidities (coronary artery disease, peripheral vascular disease, left ventricular hypertrophy, heart failure, cerebrovascular disease, and malignancy) had been collected at start of KRT as well as data on blood pressure and body mass index (BMI). The cut-off value for obesity was set at BMI >30 kg/m2 and for underweight at BMI <20 kg/m2.

The patients were followed from the start of KRT until death (55.1%), loss to follow-up (0.1%), moving abroad (0.2%), recovery of kidney function (1.9%), or end of follow-up on 31 December 2017 (42.7%). We divided the patients into six groups according to their primary kidney disease: type 2 diabetes, type 1 diabetes, glomerulonephritis (GN), polycystic kidney disease (PKD), nephrosclerosis, and other or unknown diagnoses. Information on comorbidities was available for 89% (left ventricular hypertrophy) to 97% (coronary artery disease) of patients. Information on all comorbidities was available for 85% of patients. Multiple imputation was used for sensitivity analyses.

The Finnish Registry for Kidney Diseases collects written informed consent from all patients for use of data for research purposes. The nephrological centers in Finland provide data on all patients at start of chronic KRT and once a year thereafter to the Finnish Registry for Kidney Diseases and these patient records are stored in the registry database. This study was approved by the Board of the Finnish Registry for Kidney Diseases and conducted in accordance with the Data Protection Law. The patient records used here were pseudonymized, the data were encrypted, and only specified persons of the research group had access to the database. Thus, individuals cannot be identified from the results, and the study was based entirely on previously collected registry data. The dataset for this study was retrieved from the database of the Finnish Registry for Kidney Diseases in November 2018.

Statistical methods

Between-group comparisons were performed using χ2 test for categorical variables and the Kruskal-Wallis test for continuous variables. Survival probabilities were calculated using Kaplan-Meier curves, and differences between groups were assessed using the log-rank test. Relative risks of death as a function of comorbidity were estimated using Cox proportional hazards regression with and without adjustment for potential confounding factors. Death was the event and patients were censored at the last day of follow-up or at the latest on 31 December 2017. Patients were not censored at kidney transplantation in the main analyses. Interactions between diagnosis group and comorbidities were assessed by including interaction terms in Cox regression with and without adjustment to determine whether the effect of comorbidity on survival differs according to diagnosis group. Two-sided p values lower than 0.05 were considered significant.

Statistical analyses were performed using SPSS Statistics, version 25.0.

Results

Baseline characteristics

Of the 8696 patients (41 793 patient-years), 4793 died during follow-up. Table 1 presents baseline characteristics, follow-up, and survival by diagnosis group. Median age at start of KRT was the lowest, 47.3 years, among patients with type 1 diabetes and the highest, 72.7 years, among patients with nephrosclerosis. Median survival time of all patients was 5.4 years, ranging from 3.3 years in patients with type 2 diabetes to 15.4 years in patients with PKD. The most common cause of death was cardiovascular disease, accounting for 42% of the deceased patients. Hemodialysis was the initial treatment modality in 75% and peritoneal dialysis in 25% of patients. Of all patients, 2% started in home hemodialysis and 33% received a kidney transplant during follow-up.

thumbnail
Table 1. Baseline characteristics, follow-up, and survival according to primary kidney disease.

https://doi.org/10.1371/journal.pone.0256522.t001

More than half of the patients had at least one comorbidity (coronary artery disease, peripheral vascular disease, left ventricular hypertrophy, heart failure, cerebrovascular disease, or malignancy). Left ventricular hypertrophy (in 33%) and coronary artery disease (in 25%) were the most common comorbidities. Patients with type 2 diabetes or nephrosclerosis as the cause of kidney failure had significantly more comorbidities than other patients, and the number of comorbidities was the lowest in patients with PKD. Obesity was more common in those with type 2 diabetes than in other groups. Systolic blood pressure was elevated in two-thirds and diastolic blood pressure in one-quarter of patients despite the high prevalence of medication for hypertension (Table 2).

thumbnail
Table 2. Proportion of patients with comorbidities and medication at start of KRT according to primary kidney disease.

https://doi.org/10.1371/journal.pone.0256522.t002

Comorbidities and relative risk of death

Without adjustment, all comorbidities, obesity, and underweight were associated with increased risk of death, and only elevated blood pressure was associated with decreased risk of death regardless of whether patients were on medication for high blood pressure. The relative risk of death that was associated with each comorbidity varied significantly according to primary kidney disease. Relative risk of death associated with presence of heart failure varied the most, from 1.91 (95% CI 1.66–2.19) in patients with type 2 diabetes to 6.61 (95% CI 3.57–12.22) in patients with PKD, when the risk of death was compared with patients who did not have heart failure within the same diagnosis group (S1 Table).

When adjusted for age and sex, primary kidney disease significantly modified the effect of coronary artery disease, heart failure, and malignancy on death, as all of these increased risk of death more in PKD and GN patients than in other patients (S2 Table). Obesity was associated with increased risk of death only in patients with GN. Underweight was also associated with worse prognosis, especially among patients with type 1 diabetes or nephrosclerosis.

In the multivariable model (Table 3), coronary artery disease, heart failure, and higher age had the strongest effect on worse prognosis in patients with GN or PKD. Malignancy had the most marked worsening effect on prognosis in patients with PKD. Similar to the less adjusted model, in our multivariable model the risk of death that was associated with underweight was the highest in patients with type 1 diabetes, and obesity was associated with worse prognosis only in GN patients.

thumbnail
Table 3. Multivariable model of comorbidities’ effect on relative risk of death according to primary kidney disease.

https://doi.org/10.1371/journal.pone.0256522.t003

Patients for whom we had all comorbidity data had better survival than patients with non-complete comorbidity data, and median survival times were 5.57 (95% CI 5.32–5.82) years and 4.63 (95% CI 4.25–5.01) years, respectively. Multiple imputation of missing comorbidity data did not change our results. When patients were censored at time of kidney transplantation, the differences in comorbidities’ association with risk of death according to primary kidney disease remained similar to results without censoring at kidney transplantation.

Number of comorbidities and risk of death

Six comorbidities (coronary artery disease, peripheral vascular disease, left ventricular hypertrophy, heart failure, cerebrovascular disease, and malignancy) were included in the survival analyses on number of comorbidities. Increasing number of comorbidities was associated with higher risk of death in all primary kidney disease groups, but after adjustment for age and sex, the increase was more significant among patients with GN or PKD. Three or more comorbidities, compared with no comorbidity, increased risk of death 4.62-fold in patients with GN and 4.36-fold in patients with PKD, whereas relative risk varied from 2.25 to 2.63 in other primary kidney disease groups (S3 Table). Furthermore, survival time varied markedly in primary kidney disease groups according to number of comorbidities. In patients without comorbidities, median survival ranged from 4.9 years (patients with type 2 diabetes) to more than 18 years (PKD patients). However, if patients had three or more comorbidities the variation in median survival time was much smaller: from 1.5 years in GN patients to 2.9 years in PKD patients (Fig 1). This also indicates that the number of comorbidities is more relevant to the survival prognosis in patients with PKD and GN than in other diagnostic groups.

thumbnail
Fig 1. Effect of the number of comorbidities on unadjusted survival probability according to primary kidney disease.

https://doi.org/10.1371/journal.pone.0256522.g001

Discussion

Our study shows that comorbidities affect KRT patients’ risk of death differently depending on the type of primary kidney disease. This study included all adult patients who entered long-term KRT during 2000–2017 in Finland. A significantly greater increase in risk of death was seen in PKD and GN patients with heart failure, coronary artery disease, or malignancy than with other diagnoses, when comparing patients in the same diagnosis group with or without these comorbidities. PKD or GN patients diagnosed with heart failure before start of KRT showed a fivefold age- and sex-adjusted risk of death relative to those without this comorbidity, whereas among patients with other primary kidney diseases the increase was only twofold. PKD and GN patients had less comorbidities and better survival than other KRT patients. However, if these patients had multiple comorbidities their survival was similar to other patients with multiple comorbidities due to the greater effect of comorbidities in the survival prognosis in patients with PKD or GN than in patients with other diagnoses. To our knowledge, no previously published study has assessed this primary kidney disease-related variation in comorbidities’ effect on mortality.

This nationwide study included all patients starting KRT in Finland during the 18-years study period. Data on mortality were complete, the number of patients lost to follow-up was very low, data were collected in a similar manner throughout the study period, and data on the most important comorbidities were available for a large proportion of patients. Therefore, the possibility of selection or information bias was minimal. However, a few limitations should be noted. The population in Finland is genetically quite homogeneous and almost entirely white, which may reduce the generalizability of the findings. Another potential limitation is that comorbidity data were collected mainly from pre-existing information, that is, patients were not examined separately for registry purposes. This could cause reporting bias since more severe cases were more likely reported. We also had no information on severity of comorbidities. Results of a study from the Netherlands, however, showed that adding severity grading of several comorbid conditions did not lead to improved prognostic power [9]. Although we did not have complete information on comorbidities, the data we used were significantly more comprehensive in this respect than the data used in some previous registry studies [5,10,11].

Generally, when comparing the results of our study with results of earlier ones from other countries, it is crucial that diagnoses of primary kidney disease groups and classification of comorbid conditions are comparable. Diabetes is the leading cause of end-stage kidney disease (ESKD) in most of the countries, but incidence of ESKD caused by type 1 diabetes is higher in Finland than in other countries. The ERA-EDTA Registry, the United States Renal Data System, and previous studies have reported higher proportions of incident KRT patients with hypertension or renal vascular disease as the cause of kidney disease than we report here [5,7,12,13]. Some of the differences may emerge from differences between populations, but also classification of primary kidney disease may differ. Definitions of comorbidities and methods of data collection are likely to have an impact on reported numbers of comorbidities. Previous studies have described a higher incidence of peripheral vascular disease and heart failure at the start of KRT, whereas incidence of coronary artery disease, cerebrovascular disease, and malignancies has been similar to our study [2,4,5,7,14].

All six major comorbidities reported here have also earlier been shown to be associated with increased mortality [2,6,1416]. Of note, van Manen et al. [7] showed that the influence of comorbidities may be less important than expected, explaining only 1.9% of the variance in mortality. Underweight has previously been established to be associated with an increased risk of death, but the association between obesity and mortality among GN patients shown here, has not been published earlier [1720]. Although treatment of hypertension reduces mortality among patients on hemodialysis [21], hypertension has been demonstrated to have an inverse association with mortality in epidemiological studies [6,22], which is in line with our results. This may be due to an increased risk of death associated with intradialytic hypotension [23]. Taken together, our results on the entire KRT population resemble those of previous studies, and this supports the generalizability of our findings on how primary kidney disease modifies the effect of comorbidities on survival.

How could the results of our study be interpreted? Lower age in GN and PKD patients could be one explanation. If malignancy is diagnosed at a younger age, it is more likely to be more aggressive than among older patients, thus worsening the prognosis more in a young population. In addition, conservative treatment for ESKD may have been chosen more often for elderly patients who have malignancy with a poor prognosis. However, although patients with type 1 diabetes in our study were younger than GN and PKD patients, their risk of death associated with malignancy was lower. Lower risk of malignancy-related mortality in patients with diabetes than in patients with GN has been reported previously and may be caused by a higher risk of cardiovascular death in patients with diabetes as a competing factor [24]. Cardiovascular morbidity is high in patients with diabetes and in elderly patients and some of these patients may have undiagnosed coronary artery disease and heart failure. Consequently, had some of these patients been misclassified as not having these comorbidities, this would have diluted the association with comorbidity and mortality and could explain the stronger association with mortality in GN and PKD patients. In addition, higher overall risk of death in patients with diabetes and in elderly patients may reduce the additional effect of comorbidities on mortality. Furthermore, severity and type of comorbid disease at time of diagnosis and the different treatments for these comorbidities may have influenced our results.

In conclusion, comorbidities affect KRT patients’ survival differently depending on the type of primary kidney disease. When assessing prognosis of patients with GN or PKD, especially cardiological problems and malignancies should be taken into account. Our findings should be considered when building prognostic models that include comorbidities, because otherwise these models are not accurate, and the estimates calculated with these models may lead to wrong conclusions for some patients. Future studies could shed more light on potential differences in severity or nature of comorbidities according to primary kidney disease. To avoid known shortcomings of an observational study like ours, a randomized study could reveal whether systematic screening and treatment of comorbidities in patients starting KRT would improve their survival.

Supporting information

S1 Table. Unadjusted effect of comorbidities on relative risk of death according to primary kidney disease.

https://doi.org/10.1371/journal.pone.0256522.s001

(DOCX)

S2 Table. Age- and sex-adjusted effect of comorbidities on relative risk of death according to primary kidney disease.

https://doi.org/10.1371/journal.pone.0256522.s002

(DOCX)

S3 Table. Age- and sex-adjusted relative risk of death according to the number of comorbidities in various groups of primary kidney disease.

https://doi.org/10.1371/journal.pone.0256522.s003

(DOCX)

Acknowledgments

The authors acknowledge support from the Board of the Finnish Registry for Kidney Diseases and all the nephrologists and staff in all Finnish central hospitals that have reported to the Finnish Registry for Kidney Diseases: Helsinki University Central Hospital, Turku University Central Hospital, Satakunta Central Hospital, Kanta-Häme Central Hospital, Tampere University Central Hospital, Päijät-Häme Central Hospital, Kymenlaakso Central Hospital, Etelä-Karjala Central Hospital, Mikkeli Central Hospital, Itä-Savo Central Hospital, Pohjois-Karjala Central Hospital, Kuopio University Central Hospital, Keski-Suomi Central Hospital, Etelä-Pohjanmaa Central Hospital, Vaasa Central Hospital, Keski-Pohjanmaa Central Hospital, Oulu University Central Hospital, Kainuu Central Hospital, Länsi-Pohja Central Hospital, Lappi Central Hospital, and Åland Central Hospital.

References

  1. 1. Khan IH, Catto GR, Edward N, Fleming LW, Henderson IS, MacLeod AM. Influence of coexisting disease on survival on renal-replacement therapy. Lancet. 1993; 341: 415–418. pmid:8094182
  2. 2. Goodkin DA, Bragg-Gresham JL, Koenig KG, Wolfe RA, Akiba T, Andreucci VE, et al. Association of Comorbid Conditions and Mortality in Hemodialysis Patients in Europe, Japan, and the United States: The Dialysis Outcomes and Practice Patterns Study (DOPPS). J Am Soc Nephrol. 2003; 14: 3270–3277. pmid:14638926
  3. 3. Helve J, Haapio M, Groop PH, Grönhagen-Riska C, Finne P. Comorbidities and survival of patients with type 1 diabetes on renal replacement therapy. Diabetologia. 2011; 54: 1663–1669. pmid:21465326
  4. 4. Ceretta ML, Noordzij M, Luxardo R, De Meester J, Abad Diez JM, Finne P, et al. Changes in co-morbidity patterns in patients starting renal replacement therapy in Europe–data from the ERA-EDTA Registry. Nephrol Dial Transplant. 2018; 33:1794–1804. pmid:29361126
  5. 5. Wagner M, Ansell D, Kent DM, Griffith JL, Naimark D, Wanner C, et al. Predicting mortality in incident dialysis patients: an analysis of the United Kingdom Renal Registry. Am J Kidney Dis. 2011; 57: 894–902. pmid:21489668
  6. 6. van Walraven C, Austin PC, Knoll G. Predicting potential survival benefit of renal transplantation in patients with chronic kidney disease. CMAJ. 2010; 182: 666–672. pmid:20351122
  7. 7. van Manen JG, van Dijk PC, Stel VS, Dekker FW, Cleries M, Conte F, et al. Confounding effect of comorbidity in survival studies in patients on renal replacement therapy. Nephrol Dial Transplant. 2007; 22: 187–195. pmid:16998216
  8. 8. Anderson RT, Cleek H, Pajouhi AS, Bellolio MF, Mayukha A, Hart A. Prediction of Risk of Death for Patients starting Dialysis: A Sistematic Review and Meta-Analysis. Clin J Am Soc Nephrol. 2019; 14: 1213–1227. pmid:31362990
  9. 9. van Manen JG, Korevaar JC, Dekker FW, Boeschoten EW, Bossuyt PM, Krediet RT; NECOSAD Study Group. Netherlands Co-operative Study on the Adequacy of Dialysis-2. How to adjust for comorbidity in survival studies in ESRD patients: a comparison of different indices. Am J Kidney Dis. 2002; 40: 82–89. pmid:12087565
  10. 10. Collier T, Steenkamp R, Tomson C, Caskey F, Ansell D, Roderick P, et al. Patterns and effects of missing comorbidity data for patients starting renal replacement therapy in England, Wales and Northern Ireland. Nephrol Dial Transplant. 2011; 26: 3651–3658. pmid:21436380
  11. 11. Rao A, Steenkamp R, Caskey F. UK Renal Registry 16th annual report: chapter 5 comorbidities and current smoking status amongst patients starting renal replacement therapy in England, Wales and Northern Ireland from 2011 to 2012. Nephron Clin Pract. 2013; 125: 99–110. pmid:24662169
  12. 12. ERA-EDTA Registry: ERA–EDTA Registry Annual Report 2016. Academic Medical Center, Department of Medical Informatics, Amsterdam, The Netherlands: 2018.
  13. 13. U.S. Renal Data System. 2018 USRDS annual data report: Epidemiology of kidney disease in the United States. National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases, Bethesda, MD, 2018.
  14. 14. Kotwal S, Webster AC, Cass A, Gallagher M. Comorbidity recording and predictive power of comorbidities in the Australia and New Zealand dialysis and transplant registry compared with administrative data: 2000–2010. Nephrology (Carlton). 2016; 21: 930–937. pmid:26636746
  15. 15. Fotheringham J, Jacques RM, Fogarty D, Tomson CR, El Nahas M, Campbell MJ. Variation in centre-specific survival in patients starting renal replacement therapy in England is explained by enhanced comorbidity information from hospitalization data. Nephrol Dial Transplant. 2014; 29: 422–430. pmid:24052459
  16. 16. Stack AG, Saran R. Clinical correlates and mortality impact of left ventricular hypertrophy among new ESRD patients in the United States. Am J Kidney Dis. 2002; 40: 1202–1210. pmid:12460039
  17. 17. Leavey SF, McCullough K, Hecking E, Goodkin D, Port FK, Young EW. Body mass index and mortality in ‘healthier’ as compared with ‘sicker’ haemodialysis patients: results from the Dialysis Outcomes and Practice Patterns Study (DOPPS). Nephrol Dial Transplant. 2001; 16: 2386–2394. pmid:11733631
  18. 18. Chazot C, Gassia JP, Di Benedetto A, Cesare S, Ponce P, Marcelli D. Is there any survival advantage of obesity in Southern European haemodialysis patients? Nephrol Dial Transplant. 2009; 24: 2871–2876. pmid:19369686
  19. 19. Miskulin DC, Meyer KB, Martin AA, Fink NE, Coresh J, Powe NR, et al. Comorbidity and its change predict survival in incident dialysis patients. Am J Kidney Dis. 2003; 41: 149–161. pmid:12500232
  20. 20. Ladhani M, Craig JC, Irving M, Clayton PA, Wong G. Obesity and the risk of cardiovascular and all–cause mortality in chronic kidney diasese: a systematic review and meta–analysis. Nephrol Dial Transplant. 2017; 32: 439–449. pmid:27190330
  21. 21. Aftab R, Sellappans R, Ming C, Shaik I. Taking a step further in identifying ideal blood pressure range among hemodialysis patients: A systematic review and meta-analysis. Front Pharmacol. 200; 11: 729. pmid:32528285
  22. 22. 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. pmid:12631061
  23. 23. Sars B, van der Sande FM, Kooman JP. Intradialytic hypotension: Mechanisms and outcome. Blood Purif. 2020; 49: 158–167. pmid:31851975
  24. 24. Vogelzang JL, van Stralen KJ, Noordzij M, Diez JA, Carrero JJ, Couchoud C, et al. Mortality from infections and malignancies in patients treated with renal replacement therapy: data from the ERA-EDTA Registry. Nephrol Dial Transplant. 2015; 30: 1028–1037. pmid:25637641