Research | Open | Published:
Coronary artery disease screening and prognosis in incident dialysis patients
Renal Replacement Therapyvolume 2, Article number: 22 (2016)
Guidelines have recommended that physicians screen for coronary artery disease (CAD) at the initiation of dialysis. The purpose of this study was to examine the effects of CAD screening at the initiation of dialysis on prognosis after starting dialysis
This retrospective cohort study involved 224 consecutive incident dialysis patients without advanced cardiac disease. Pharmacological stress myocardial perfusion imaging (MPI) was performed to routinely screen incident dialysis patients within 3 months of starting dialysis therapy. The influence of MPI on prognosis after starting dialysis was explored by the propensity score adjustment method. The study outcome was defined as all-cause mortality and cardiac death.
Of the 224 patients, 164 (73.2 %) underwent MPI screening at initiation of dialysis. During the median follow-up period of 5.5 years, 77 patients died, among whom 17 (22.1 %) cardiac deaths were observed. The patients who received MPI had a significantly higher survival rate for all-cause mortality and cardiac death compared with the patients without MPI (log-rank test, p < 0.001, p = 0.021, respectively). In multiple Cox regression analysis, MPI screening continued to have a statistically significant negative association with all-cause mortality after adjusting for confounding variables (hazard ratio 0.581, p = 0.034) and propensity score (hazard ratio 0.568, p = 0.025).
MPI CAD screening at dialysis initiation could improve prognosis after starting dialysis.
Although the underlying mechanism of the cardio-renal axis in patients with chronic kidney disease (CKD) has been determined in considerable detail over the past decade , cardiovascular disease is still a leading cause of death in maintenance dialysis patients worldwide [2, 3] and also in Japan . The accelerated progression of atherosclerosis during the conservative phase of CKD, before starting dialysis therapy , is thought to be one of the causes of the high incidence of cardiovascular events in the dialysis phase [6, 7]. Therefore, it is well known that a high prevalence of atherosclerotic cardiovascular diseases is observed at the initiation of dialysis. For example, approximately 50 % of incident asymptomatic dialysis patients have been observed to suffer from advanced coronary artery disease (CAD) [8, 9] and show poor prognosis after starting dialysis [6, 10]. Based on this background, guidelines created by several organizations [11–13] have recommended that physicians screen for CAD at the initiation of dialysis in order to find it earlier and improve prognosis. However, it has not yet been established whether screening for CAD plays a crucial role in improving the prognosis for dialysis patients.
Pharmacological stress myocardial perfusion imaging (MPI), even in end-stage kidney disease (ESKD) patients, is one of the most frequently performed procedures for identifying patients with myocardial ischemia and CAD, and for stratifying the patients at high risk for poor prognosis . Recently, the severity of abnormal MPI results has begun to have a great impact on the prediction of cardiac events in the conservative phase of advanced CKD patients without cardiovascular disease . Therefore, MPI is becoming a useful screening procedure for patients at any stage of kidney disease as a non-invasive procedure that is not harmful to kidney function.
Our hypothesis for the present study was that screening for CAD at the initiation of dialysis by using MPI may lead to improved prognosis after starting dialysis. Therefore, this study examined the effect of MPI screening for CAD on the prognosis of incident dialysis patients.
Study design and patients
This was a retrospective cohort study that involved 314 consecutive patients who started maintenance dialysis for ESKD at the Toho University Ohashi Medical Center between January 2001 and December 2010. In order to precisely evaluate the effect of MPI screening on prognosis, ESKD patients who already had cardiac disease before starting dialysis were excluded. As shown in Fig. 1, 74 patients met the exclusion criteria as follows: (1) history of New York Heart Association (NYHA) level III or IV heart failure, (2) coronary heart disease such as myocardial infarction and/or coronary revascularization, or (3) left ventricular dysfunction defined as ejection fraction (EF) <50 % evaluated by echocardiography. Additionally, 16 patients who expired within 3 months of starting dialysis were excluded. The remaining 224 patients were enrolled into this study, among whom 164 underwent MPI as screening for CAD. The influence of MPI on prognosis after starting dialysis was explored by crude cohort analysis and by the propensity score adjustment method.
All patients provided informed consent to participate in the study, and the Ethics Committee for Clinical Research at Toho University Ohashi Medical Center approved the study protocol [Permission no. 13-52, 13-61]. This study adhered to the Declaration of Helsinki.
Clinical information from all patients was recorded before the first hemodialysis (HD) session. All patients were initially interviewed to determine age, sex, type of primary renal disease, previous hospitalization, and history of hypertension and cardiac disease. We also collected information about the medications the patients were receiving after starting dialysis. Blood pressure was recorded from the supine position, and blood samples were withdrawn immediately before the first HD session. Serum calcium (Ca) concentrations were corrected by albumin levels as follows: corrected Ca (cCa) = measured Ca (mg/dL) + (4—measured albumin level (g/dL)). The estimated glomerular filtration rate (eGFR) was calculated using the Japanese equation as follows : eGFR (mL/min/1.73 m2) = 194 × Cr − 1.094 × Age − 0.287 (for women × 0.739). Body mass index (BMI) at optimal weight was calculated as weight (kg) divided by height (m) squared. All participants were assessed by echocardiography at optimum weight after a dialysis session or on an intermittent day of hemodialysis before discharge. Ejection fraction was measured as a marker of left ventricular (LV) dysfunction.
MPI as a screening method for coronary artery disease
From January 2001, ESKD patients at our hospital have been routinely screened for CAD within 3 months of starting HD therapy, regardless of any suspicion for ischemic heart disease (IHD). Pharmacological stress MPI by using thallium-201 was performed for screening in 164 of the 224 ESKD study patients. Clinically significant CAD was defined when myocardial perfusion defects were reversible in patients examined by pharmacological stress thallium-201 MPI.
The outcome of this study was defined as all-cause mortality and cardiac death. Cardiac death was defined as death from heart failure, myocardial infarction, or sudden death. Information about death was obtained from a questionnaire survey sent to the maintenance hemodialysis facilities or the hospital medical records. The cohort was followed until the end of 2012.
Data are summarized as number, prevalence, arithmetic mean ± standard deviation, or median (interquartile range), as appropriate. Comparisons of prevalence and values between groups were performed by Student’s t test and the χ 2 test. In the prospective analysis, the Kaplan-Meier method was used with the log-rank test. Prognostic variables were examined by Cox proportional hazard models, and hazard ratios and 95 % confidence intervals are reported. p values <0.05 were considered significant. To adjust for differences in baseline characteristics between the two patient groups (with and without screening MPI), a propensity score analysis was performed using multivariate logistic regression modeling, including all baseline covariates such as age, sex, diabetes, BMI, systolic blood pressure (SBP), LV EF, albumin, hemoglobin (Hb), total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-C), eGFR, cCa, phosphate (P), intact-parathyroid hormone (iPTH), c-reactive protein (CRP), and usage of three medications: renin angiotensin aldosterone system inhibitor (RAS-I), beta blocker (BB), and statin. The score was then incorporated into a Cox proportional hazards model as a covariate. All statistical analyses were performed using SPSS for Windows version 20 (IBM, Tokyo, Japan).
Of the 224 study patients, 164 (73.2 %) patients underwent MPI screening at the initiation of dialysis (Table 1). Significantly higher incidence of diabetes, higher diastolic pressure, and higher creatinine concentration (lower eGFR) were observed in the MPI screening group compared with the non-MPI group. Although the MPI group had a younger age, more males, higher SBP, lower Hb, higher HDL-C, lower cCa, and higher P than the non-MPI group, no statistically significant differences were found. Among the 164 patients who underwent the screening, 27 patients (16.5 %) showed abnormal MPI. The usage of three medications (statin, RAS-I, and BB) did not differ significantly between the two groups.
The median follow-up period was 5.5 years. During the follow-up period, 77 patients died, including 17 cardiac deaths (22.1 %), as shown in Table 2. A total of 31 patients were lost to follow-up, and two patients received renal transplantation and were censored at the time of the last observation, which was the operation day for renal transplantation (Fig. 1). At the end of the follow-up period on December 31, 2012, a total of 114 patients were alive on maintenance hemodialysis and had completed follow-up.
As shown in Fig. 2a, the Kaplan-Meier curves indicated that the patients with MPI screening had a significantly higher survival rate than the patients who did not undergo MPI screening (log-rank test, p < 0.001). The cumulative survival rates at 5 years after starting dialysis with and without MPI were 74.9 and 49.1 %, respectively. Regarding cardiac death in Fig. 2b, the cardiac death event-free 5-year survival rate was 94.8 % in patients with the screening and 86.6 % in patients without the screening (log-rank test, p = 0.021). After adjusting by propensity score, the Kaplan-Meier curves still showed higher survival in patients with MPI than in those without MPI (dotted line in Fig. 2a, b).
Figure 3 depicts the survival curves in three patient groups: 27 patients with CAD, 137 patients without CAD, and 60 patients without screening. Interestingly, the 5-year survival rate was the same for patients with and without CAD (74.9 vs. 74.9 %). After 5 years, the survival curve of patients with CAD steeply decreased and became close to the same level as the patients without the screening at 10 years (26.7 vs. 27.4 %) (Fig. 3a). Similar findings were confirmed in the survival curves for cardiac death (Fig. 3b).
Coronary angiography and revascularization
Coronary angiography was performed within 180 days in 9 (33.3 %) of 27 with abnormal MPI findings. Seven of nine with coronary angiography had significant CAD that was treated with coronary revascularization therapy (CRT) (one coronary artery bypass surgery and six percutaneous transluminal coronary angioplasty).
Cox proportional hazards model for the factors associated with death
In the univariate Cox proportional hazards model, age, SBP, serum albumin, Hb, creatinine, eGFR, cCa, CRP, usage of RAS-I, and MPI screening showed significant positive or negative associations with death (Table 3). CRT did not show any significant association with death (Table 3). In the multiple Cox regression analysis, MPI screening was significantly associated with all-cause mortality after adjusting for age, SBP, albumin, hemoglobin, eGFR, calcium, CRP, and usage of RAS-I (Table 4). Furthermore, a significant association with all-cause mortality was observed for MPI screening after adjusting the propensity score (Table 4).
Age, triglyceride, creatinine, eGFR, corrected Ca, CRP, and MPI screening showed significant positive or negative associations with cardiac death in the univariate Cox proportional hazards model. CRT did not show any significant association with cardiac death. In multiple Cox regression analysis, MPI screening was negatively associated with cardiac death after adjusting for age, eGFR and, CRP (Table 4). Furthermore, a negative association with cardiac death was observed for MPI screening after adjusting the propensity score (Table 4). However, neither of these negative associations was statistically significant.
Main findings of this study
Several studies have observed a high prevalence of CAD at the initiation of dialysis [8, 9], with severe CAD having a large influence on poor prognosis after starting dialysis . Therefore, international [11, 12] and Japanese  guidelines have recommended that physicians screen patients for coronary heart disease at the initiation of dialysis in order to improve their prognosis. However, to date, no evidence has been obtained to prove the effect of screening on prognosis in dialysis patients. The present study demonstrated that MPI screening for CAD at the initiation of dialysis may lead to improved prognosis in incident dialysis patients without cardiac disease. An approximately 50 % risk reduction for death produced by the screening was confirmed by using crude cohort analysis, and also by using propensity score adjustment analysis. Our findings may serve as evidence to support the international guidelines [11, 12].
Effect of CAD screening on mortality
Several previous studies have shown the effect of screening on prognosis in the field of cardiovascular disease. Thompson et al. reported the elegant results of their randomized controlled study, in which ultrasound screening for abdominal aortic aneurysm (AAA) resulted in a reduction in all-cause mortality and was beneficial for AAA-related mortality in a population-based sample of men aged 65–74 years . A relative 42 % risk reduction for AAA-related death and 3 % risk reduction in all-cause mortality were confirmed in the screening group compared with the non-screening group. Similarly, Lindholt et al. conducted a randomized single-center mass screening trial for AAA . Screening reduced the number of emergency operations by 68 % and lowered the mortality of AAA by 73 % over a 10-year period in the screening group compared with the control group. Although these reports demonstrate the benefit of screening for AAA in the asymptomatic general population, which is different from the aim of the present study, it is not surprising that the same beneficial findings are confirmed by screening for CAD even in ESKD patients without cardiac disease.
Conversely, several reports have demonstrated that screening for disease had no beneficial effects on mortality. In a randomized controlled trial , ultrasound screening for AAA was found to be not effective in men aged 65–83 years and did not reduce overall death rates. Another study  showed that echocardiographic screening for structural and valvular heart disease in the general population provided no benefit for mortality or for the risk of myocardial infarction or stroke. The reasons for the contrary outcomes between those studies remain unclear; however, some potential underlying reasons might exist. The prevalence of an implicit target disease for screening in a study cohort would be an important factor affecting the outcome for screening success. If the prevalence of a target cardiovascular disease is very low in a study cohort, there is a high possibility that the screening will not show any beneficial effect if the study cohort is not enough large. Indeed the clinical guideline in American college of physicians has not recommended screening for cardiac disease in asymptomatic, low-risk adults with resting or stress electrocardiography, stress echocardiography, or stress myocardial perfusion imaging . The accuracy of a screening test would also play an important role in the study outcome. A lack of sufficient sensitivity and specificity to capture a target disease would crucially influence the results. It is well known that a high prevalence (40–50 %) of underlying CAD is observed in asymptomatic ESKD patients [8, 9, 22]. Even in ESKD patients, the accuracy of the pharmacological stress MPI for detecting CAD is reported to be sufficiently high, with a sensitivity of 80 % and specificity of 73 % observed in 42 chronic dialysis patients  and a sensitivity of 86 % and specificity of 79 % in 80 type I diabetes ESKD patients . Also, the event rate after screening influences the study outcome. Young et al. assessed whether routine screening for CAD by MPI affects cardiac outcomes in asymptomatic type 2 diabetic patients . These investigators concluded that cardiac event rates were not significantly reduced by MPI screening over a follow-up period of 4.8 years due at least in part to a low cardiac event rate of 2.9 %. In the present study, during the median follow-up period of 5.5 years, the cardiac death rate was 22.1 %, which is seven times higher than that in the Young study. This is also one of the reasons for the beneficial effect of screening. Therefore, although we used a small sample of ESKD patients in a single-center study, we believe that the findings of the present study provide valuable data.
Greater beneficial effect for all-cause mortality than for cardiac death
In general, a greater beneficial effect is expected for cardiac death than for all-cause mortality by MPI screening. In the present study, MPI screening had a negative impact on cardiac death; however, this did not reach a statistically significant level in multivariate Cox proportional hazards analysis. Few cardiac death events during the follow-up period may lead to insufficient statistical power for analysis. However, a possible explanation can be raised for this issue. The first possible explanation is that the patients with MPI screening may be a special sub-cohort who pays more attention to their own health, compared to the patients without the screening. Although it is a matter of speculation, screening for CAD at the initiation of dialysis could be a trigger to alert patients to be vigilant regarding their overall health. This would be a reason why the screening might exert a more beneficial effect on improving prognosis for all-cause mortality than for cardiac death.
The second possible explanation may relate to the results of our previous study . In that study, we analyzed the survival rates for major adverse cardiac events in incident dialysis patients without CAD and found that even in patients without significant CAD, cardiac events occurred after 1 year of starting dialysis. We have proposed to administer a second screening for CAD at least 2 years after starting dialysis, even if no significant CAD was confirmed at the initiation of dialysis. It is conceivable that a higher number of physician visits may be made by patients who underwent a first screening compared with the patients without screening; consequently, the mortality rate may decrease in the patients who undergo screening for CAD.
The third possible explanation is attempts to lower the rate of cardiac death. Twenty-seven patients among the 164 (16.4 %) who received screening suffered from CAD in the present study. Physicians may pay more attention and add intensive treatment to those patients to improve prognosis. In fact, it was found that the 5-year survival rate of the 27 ESKD patients with CAD was 74.9 %, which is not inferior to the 5-year survival rate of 74.9 % in the 137 ESKD patients without CAD. The 1-year survival rate showed a large difference (26.7 vs. 64.5 %) between patients with and without CAD. On the other hand, in the 60 ESKD patients without the screening, approximately 10 patients (about 16.5 %) were calculated to have suffered from CAD if the same prevalence in the 164 ESKD patients with the screening had occurred. Failure to pay attention and insufficient management in those 10 patients may have played a role, at least in part, in poor prognosis for all-cause mortality. However, the possible advantages of screening on prognosis could not be fully explained based only on the difference in management of CAD.
The fourth possible explanation of greater benefit for all-cause mortality than for cardiac death is the potential for selection bias. The design of our study was an observational cohort study, not a randomized prospective study. It cannot be denied that relatively low-risk patients were selected and underwent the screening by chance. In order to minimize the selection bias, according to the analytical technique used in the present study, we calculated the propensity score by using all background factors to adjust for differences in baseline characteristics. As shown in Table 4, an approximately 50 % risk reduction for death engendered by the screening appeared even after adjusting the propensity score. Therefore, it is more plausible that screening at the initiation of dialysis could lead to a better prognosis for incident dialysis than ever before.
The timing of screening in CKD patients
To improve the prognosis of hemodialysis patients, the best goal we should reach is to start dialysis without any cardiovascular complication. Therefore, we think that the screening should be performed in earlier predialytic phase of CKD such as CKD 3 or 4. The initiation of HD as timing for screening may be too late. This would be a reason why the beneficial effect of MPI screening at the initiation of dialysis on cardiac mortality in this study showed marginal.
The present study had some limitations. In order to accurately evaluate the effects of screening on prognosis, a prospective randomized controlled study is the gold standard. However, based on the recommendations of international guidelines [11, 12], it is impossible to plan a study in which the patients are randomly divided according to whether screening is performing. We felt constrained to use the propensity score adjustment method to prove our hypothesis. Another limitation was that this was a single-center observational study with a relatively small number of patients. A larger-scale study should be conducted in the future.
Ronco C, McCullough P, Anker SD, Anand I, Aspromonte N, Bagshaw SM, et al. Cardio-renal syndromes: report from the consensus conference of the acute dialysis quality initiative. Eur Heart J. 2010;31:703–11.
Investigators ET, Chertow GM, Block GA, Correa-Rotter R, Drueke TB, Floege J, et al. Effect of cinacalcet on cardiovascular disease in patients undergoing dialysis. N Engl J Med. 2012;367:2482–94.
Wheeler DC, London GM, Parfrey PS, Block GA, Correa-Rotter R, Dehmel B, et al. Effects of cinacalcet on atherosclerotic and nonatherosclerotic cardiovascular events in patients receiving hemodialysis: the EValuation Of Cinacalcet HCl Therapy to Lower CardioVascular Events (EVOLVE) trial. J Am Heart Assoc. 2014;3:e001363.
Moroi M, Tamaki N, Nishimura M, Haze K, Nishimura T, Kusano E, et al. Association between abnormal myocardial fatty acid metabolism and cardiac-derived death among patients undergoing hemodialysis: results from a cohort study in Japan. Am J Kidney Dis. 2013;61:466–75.
Shoji T, Emoto M, Tabata T, Kimoto E, Shinohara K, Maekawa K, et al. Advanced atherosclerosis in predialysis patients with chronic renal failure. Kidney Int. 2002;61:2187–92.
Hase H, Tsunoda T, Tanaka Y, Takahashi Y, Imamura Y, Ishikawa H, et al. Risk factors for de novo acute cardiac events in patients initiating hemodialysis with no previous cardiac symptom. Kidney Int. 2006;70:1142–8.
Tanaka Y, Joki N, Hase H. History of acute coronary events during the predialysis phase of chronic kidney disease is a strong risk factor for major adverse cardiac events in patients initiating haemodialysis. Nephrol Dial Transplant. 2007;22:2917–23.
Joki N, Hase H, Nakamura R, Yamaguchi T. Onset of coronary artery disease prior to initiation of haemodialysis in patients with end-stage renal disease. Nephrol Dial Transplant. 1997;12:718–23.
Ohtake T, Kobayashi S, Moriya H, Negishi K, Okamoto K, Maesato K, et al. High prevalence of occult coronary artery stenosis in patients with chronic kidney disease at the initiation of renal replacement therapy: an angiographic examination. J Am Soc Nephrol. 2005;16:1141–8.
Joki N, Hase H, Takahashi Y, Ishikawa H, Nakamura R, Imamura Y, et al. Angiographical severity of coronary atherosclerosis predicts death in the first year of hemodialysis. Int Urol Nephrol. 2003;35:289–97.
K/DOQI. Clinical practice guidelines for cardiovascular disease in dialysis patients. Am J Kidney Dis. 2005;45:S1–S153.
Herzog CA, Asinger RW, Berger AK, Charytan DM, Diez J, Hart RG, et al. Cardiovascular disease in chronic kidney disease. A clinical update from Kidney Disease: Improving Global Outcomes (KDIGO). Kidney Int. 2011;80:572–86.
Hirakata H, Nitta K, Inaba M, Shoji T, Fujii H, Kobayashi S, et al. Japanese Society for Dialysis Therapy guidelines for management of cardiovascular diseases in patients on chronic hemodialysis. Ther Apher Dial. 2012;16:387–435.
Hase H, Joki N, Ishikawa H, Fukuda H, Imamura Y, Saijyo T, et al. Prognostic value of stress myocardial perfusion imaging using adenosine triphosphate at the beginning of haemodialysis treatment in patients with end-stage renal disease. Nephrol Dial Transplant. 2004;19:1161–7.
Joki N, Hase H, Kawano Y, Nakamura S, Nakajima K, Hatta T, et al. Myocardial perfusion imaging for predicting cardiac events in Japanese patients with advanced chronic kidney disease: 1-year interim report of the J-ACCESS 3 investigation. Eur J Nucl Med Mol Imaging. 2014;41:1701–9.
Matsuo S, Imai E, Horio M, Yasuda Y, Tomita K, Nitta K, et al. Revised equations for estimated GFR from serum creatinine in Japan. Am J Kidney Dis. 2009;53:982–92.
Thompson SG, Ashton HA, Gao L, Buxton MJ, Scott RA, Multicentre Aneurysm Screening Study G. Final follow-up of the Multicentre Aneurysm Screening Study (MASS) randomized trial of abdominal aortic aneurysm screening. Br J Surg. 2012;99:1649–56.
Lindholt JS, Juul S, Fasting H, Henneberg EW. Preliminary ten year results from a randomised single centre mass screening trial for abdominal aortic aneurysm. Eur J Vasc Endovasc Surg. 2006;32:608–14.
Norman PE, Jamrozik K, Lawrence-Brown MM, Le MT, Spencer CA, Tuohy RJ, et al. Population based randomised controlled trial on impact of screening on mortality from abdominal aortic aneurysm. BMJ. 2004;329:1259.
Lindekleiv H, Lochen ML, Mathiesen EB, Njolstad I, Wilsgaard T, Schirmer H. Echocardiographic screening of the general population and long-term survival: a randomized clinical study. JAMA Intern Med. 2013;173:1592–8.
Chou R, High Value Care Task Force of the American College of P. Cardiac screening with electrocardiography, stress echocardiography, or myocardial perfusion imaging: advice for high-value care from the american college of physicians. Ann Intern Med. 2015;162:438–47.
Charytan D, Kuntz RE, Mauri L, DeFilippi C. Distribution of coronary artery disease and relation to mortality in asymptomatic hemodialysis patients. Am J Kidney Dis. 2007;49:409–16.
Dahan M, Legallicier B, Himbert D, Faraggi M, Aubry N, Siohan P, et al. Diagnostic value of myocardial thallium stress scintigraphy in the detection of coronary artery disease in patients undergoing chronic hemodialysis. Arch Mal Coeur Vaiss. 1995;88:1121–3.
Boudreau RJ, Strony JT, duCret RP, Kuni CC, Wang Y, Wilson RF, et al. Perfusion thallium imaging of type I diabetes patients with end stage renal disease: comparison of oral and intravenous dipyridamole administration. Radiology. 1990;175:103–5.
Young LH, Wackers FJ, Chyun DA, Davey JA, Barrett EJ, Taillefer R, et al. Cardiac outcomes after screening for asymptomatic coronary artery disease in patients with type 2 diabetes: the DIAD study: a randomized controlled trial. JAMA. 2009;301:1547–55.
Joki N, Tanaka Y, Ishikawa H, Takahashi Y, Iwakura Y, Masuda H, et al. Optimum second screening point for detection of coronary artery disease in hemodialysis patients without advanced coronary artery disease. Am J Nephrol. 2009;29:420–5.
The authors greatly appreciate the staff of the Cardiac Nuclear Medicine Department at Toho University Ohashi Medical Center for their evaluations of the myocardial perfusion imaging studies.
NJ received lecture fee from CHUGAI Pharmaceutical and KYOWA KIRIN and received research funding from Roche Diagnostics. The other authors declare that they have no competing interests.
HH and NJ had full access to all the data in the study and take responsibility for the integrity of the data and analytical accuracy. HH, NJ, and YT contributed to the study concept and design. MI, SK, TA, YT, YI, AM, MH, HN, and KH participated in the data acquisition. YT, TH, HN, and NJ performed the data analysis and interpretation. YT and NJ drafted the manuscript. YT, TH, and NJ were involved in the statistical analysis. KH gave administrative, technical, and material support. All authors read and approved the final manuscript.