- 1Department of Orthopedic Surgery, Soonchunhyang University Gumi Hospital, Gumi, South Korea
- 2Division of Nephrology, Department of Internal Medicine, CHA Gumi Medical Center, CHA University, Gumi, South Korea
- 3Division of Nephrology, Department of Internal Medicine, College of Medicine, Yeungnam University, Daegu, South Korea
Introduction: This study aimed to evaluate the association between the psoas muscle (PM) index with or without fatty infiltration and various indices associated with muscle mass in hemodialysis (HD) patients.
Methods: We included stable HD patients (n = 83). The collected data included subjective global assessment (SGA) score, ASM/Ht2 (appendicular skeletal muscle mass divided by the squared height in meters), gait speed (GS; m/s), and handgrip strength (HGS; kg). The abdominal computed tomography (CT) image was obtained using a CT scanner. The PM and normal-density PM (NPM) indices (mm2/m2) were calculated using the whole PM area or the area with an attenuation range of 35–100 HU divided by the squared height in meters.
Results: Correlation coefficients for the SGA score, ASM/Ht2, HGS, and GS were greater for the NPM index than for the PM index. The linear regression analysis showed that, on multivariate analysis, the NPM index was significantly associated with the SGA score, ASM/Ht2, and GS. However, the PM index was significantly associated with the SGA score and ASM/Ht2 but not with HGS or GS. For calculating the low GS, the area under the receiver operating characteristic curve area was significantly greater for the NPM index than for the PM index (P = 0.012).
Conclusion: The present study suggested that the NPM index excluding fatty infiltration may be an early and useful indicator for detecting muscle strength and physical performance among HD patients.
Introduction
Hemodialysis (HD) is one of most commonly used dialysis modalities in end-stage renal disease patients requiring renal replacement therapy. With advances in HD technologies and cares, the survival of HD patients has improved. However, the prevalence of long-term complications increased over time with treatment and as the dialysis population ages (1). Decreased muscle mass is one complication that can develop in long-term HD patients (2). This is associated with decreased quality of life and increased morbidity or mortality in HD patients (3). Therefore, accurate muscle mass measurement is important for predicting prognosis or physical performance in HD patients.
Computed tomography (CT) and dual X-ray absorptiometry (DXA) are commonly used for measuring muscle mass in HD patients (4, 5). Previous studies have shown that, among various muscle mass indices, the psoas muscle (PM) index measured on CT is associated with clinical outcomes in HD patients (6, 7). However, HD patients are prone to insulin resistance caused by various conditions such as diabetes mellitus (DM), chronic inflammation, or uremia (8). Insulin resistance can be associated with fatty changes in the muscle mass in HD patients and lead to overestimated muscle mass measurements. Both PM including or excluding fatty infiltration can predict the total muscle mass index and nutritional index. The two indices are positively correlated with total muscle mass. Both muscle mass and fatty infiltration can be increased by improving nutrition in dialysis patients. However, the predictability of muscle strength or physical performance differs between the two indices. Although the total amount of muscle mass was to be the same, muscle mass with fatty infiltration is associated with a greater decrease in muscle strength or physical performance than that without fatty infiltration. These reveal that muscle mass without fatty changes may be more closely associated with muscle strength or physical performance than total muscle mass per se. However, there are few studies of the association between muscle without fatty changes and indicators associated with muscle function in HD patients. This study aimed to evaluate the association between the PM index with or without fatty infiltration and various indices associated with muscle mass in HD patients.
Materials and Methods
Study Design and Participants
Our study was a cross-sectional study and it was performed between September 2012 and March 2015 in a medical center. We included HD patients aged ≥20 years and with duration of dialysis ≥6 months, ability to ambulate without an assistive device, ability to communicate with the interviewer, and no hospitalization within the 3 months before enrollment. This study was approved by the institutional review board of a local medical center (approval no. 12-07). Informed consent was obtained from all participants before enrollment: a total of 83 patients were enrolled. None of the participants took medications associated with physical activity such as opioids, antihistamines, or antidepressants.
Baseline Variables
The collected baseline data included sex, age, presence of DM, dialysis vintage, hemoglobin (g/dL), C-reactive protein (CRP; mg/dL), blood urea nitrogen (mg/dL), creatinine (mg/dL), calcium (mg/dL), phosphorus (mg/dL), sodium (mEq/L), potassium (mEq/L), chloride (mEq/L), intact parathyroid hormone (i-PTH; pg/mL), total cholesterol (mg/dL), albumin (g/dL), and single pool measurement of solute removal during HD that focuses on urea (spKt/Vurea). All laboratory tests were performed prior to the HD sessions and repeated thrice in the following 3 weeks. The mean of three values was recorded for each parameter. DM was defined as a patient-reported history and medical record of a DM diagnosis or use of antiglycemic medication. The spKt/Vurea values were calculated using the described as described previously (9).
Assessment of the Nutritional Status, Muscle Mass, Strength, and Gait Speed
The subjective global assessment (SGA) score was calculated based on scores from seven items (weight loss, dietary intake, gastrointestinal symptoms, functional capacity, comorbidity, decreased fat, and decreased muscle mass) as described previously (10). Muscle mass was evaluated using whole-body DXA (GE Medical Systems Lunar, Madison, WI, USA). Appendicular skeletal muscle mass (ASM; kg) was calculated using the sum of lean mass of both extremities. Finally, ASM/Ht2 was calculated as ASM divided by the squared height in meters. Hand grip strength (HGS) was measured and each patient performed three trials with the dominant hand using a manual hydraulic dynamometer (Jamar® Hydraulic hand dynamometer; Sammons Preston, Chicago, IL, USA). Maximum values among the three trials were selected. Gait speed (GS; m/s) was evaluated using the timed 4-m walking test (11). A low GS group was defined based on their values <0.8 m/s, according to a previous guideline (11). We also defined the low SGA (<6 of total SGA score), low ASM/Ht2 (<7.0 kg/m2 for males and <5.4 kg/m2 for females), and low HGS (<26 kg for males and <18 kg for females) groups based on the same guidelines (11).
PM Evaluation
The abdominal CT image was obtained using a 320-slice CT scanner (Aquilion ONE; Toshiba Medical Systems Corp., Tokyo, Japan). A single axial image was obtained at the lower border of the L3 vertebral level. The images were analyzed using an image analysis software (ImageJ 1.45S; National Institutes of Health, Bethesda, MD, USA). In our study, Hounsfield units (HU) of muscle were defined using a modified version described previously (12). Muscle attenuations at very low, low, or normal density were −29 to −1, 0–34, and 35–100 HU, respectively. A quantitative assessment of the whole PM was obtained and PM area (mm2) was defined as the area with an attenuation range of −29 to 100 HU among the sum of the right and left whole PM areas, which includes very low-, low-, and normal-density muscle. The PM index (mm2/m2) was calculated using the PM area divided by the squared height in meters. Normal-density PM (NPM) area was defined as the area with an attenuation range of 35–100 HU among the sum of the right and left whole PM areas including normal-density muscle alone. The NPM index (mm2/m2) was calculated using the NPM area divided by the squared height in meters. Three measurements were taken of both PMs and the average was used for the analysis.
Statistical Analysis
The data were analyzed using IBM SPSS Statistics version 25 (SPSS Inc., Chicago, IL, USA). Categorical variables are expressed as counts (percentages). The distribution of the continuous variables was evaluated using the Kolmogorov-Smirnov test. Continuous variables with normal distribution are expressed as means ± standard deviations. Those with non-normal distributions are expressed as medians (interquartile ranges). Comparisons between two groups were performed using t-tests. The correlations between two continuous variables were assessed using Pearson's correlation analysis. For linear regression analyses, the dependent variables were the SGA, ASM/Ht2, HGS, and GS, and the independent variables were the PM and NPM indices. The multivariate analysis was adjusted for age, sex, and presence of DM and performed using the backward conditional method. Discrimination, which is the ability of the model to differentiate between participants who have a low GS and those who do not, was examined using the area under the receiver operating characteristic curve (AUROC). The statistical significance between the AUROCs was calculated using the DeLong method. SPSS version 25 did not compare the AUROC curves using the DeLong method. Therefore, we also used other software to compare and evaluate the statistical significance between the ROC curves. The AUROC was calculated using MedCalc version 11.6.1.0 (MedCalc, Mariakerke, Belgium). The level of statistical significance was set at P < 0.05.
Results
The clinical characteristics of the 83 patients are shown in Table 1. The mean age was 56.7 ± 11.9 years; 44 (53%) of the participants were male. The dialysis duration was 3.0 (5.1) years. The mean SGA score, HGS, GS, and ASM/Ht2 values were 6 (2), 26.0 (9) kg, 0.92 ± 0.20 m/s, and 6.58 ± 0.98 kg/m2, respectively. The proportions of very low-, low-, or normal-density muscle among the entire PM area were 5.8 ± 3.1%, 17.2 ± 7.2%, and 77.0 ± 9.6%, respectively. The mean PM and NPM index values were 448.7 ± 138.7 mm2/m2 and 348.1 ± 126.3 mm2/m2, respectively.
Correlation coefficients for the PM index were 0.309 for the SGA score (P = 0.005), 0.580 for ASM/Ht2 (P < 0.001), 0.474 for HGS (P < 0.001), and 0.295 for GS (P = 0.007). Those for the NPM index were 0.338 for the SGA score (P = 0.002), 0.604 for ASM/Ht2 (P < 0.001), 0.537 for HGS (P < 0.001), and 0.384 for GS (P < 0.001). Because age is an important confounding factor, we also performed subgroup analyses according to 30–49, 50–59, or 60–79 years of age (Table 2). For all the groups and variables except ASM/Ht2 in the 30–49 age group, the correlation coefficients were greater in the NPM index than the PM index.
Table 2. Correlation coefficients between the PM indices and various variables according to age groups.
Comparisons between the low and normal groups are presented in Table 3. The PM and NPM indices in the low GS group were lower than those in the normal GS group. There were no significant differences in the PM and NPM indices between the low and normal groups for HGS, SGA, and ASM/Ht2, but the trends of the mean differences between the two groups were consistently greater in the NPM index than the PM index. In addition, the difference between the low and normal groups was relatively greater in HGS or GS than in ASM/Ht2.
The linear regression analysis showed that, on multivariate analysis, the NPM index was significantly associated with the SGA score, ASM/Ht2, and GS (Table 4). However, the PM index was significantly associated with the SGA score and ASM/Ht2 but not with HGS or GS. Multivariate linear regression analysis showed that statistical significance was not obtained the association between the two indices and HGS (P = 0.251 for PM index and P = 0.053 for NPM index). Non-significance between the NPM index and HGS may be associated with small sample size. However, considering lower P-values for NPM index compared to that for PM index, analyses with a larger sample size would lead to statistical significance between the NPM and HGS. Figure 1 shows the AUROCs for predicting the low GS of the two indices. For calculating the low GS, the AUROC values of the PM index were 0.675 [95% confidence interval (CI), 0.562–0.774; P = 0.007]. The AUROC values of the NPM index were 0.725 (95% CI, 0.616–0.818; P < 0.001). The AUROC area was significantly greater for the NPM index than for the PM index (P = 0.012).
Figure 1. The area under the receiver operating characteristic curve for the PM or NPM indices for predicting low gait speed. NPM, normal density-psoas muscle; PM, psoas muscle.
Discussion
Previous studies have shown the importance of muscle quantity and quality in populations prone to the risk of sarcopenia and/or insulin resistance (13, 14). Muscle quality is generally associated with fat infiltration within or between the muscles. Fatty changes in the muscle can result from various conditions, such as age, obesity, DM, or surgery, that decrease muscle attenuation on CT images (15). A range of −190 HU and −30 HU is commonly used to define the fat mass on CT images. A previous study defined HU ranges of very low-, low-, or normal-density muscle of −29 to −1, 0–34, and 35–100, respectively (12). However, there were inconsistencies regarding HU attenuation cut-off values or muscle mass on CT. Low cut-off values for muscle mass ranged from −29 to 35 HU, while high cut-off values ranged from 100 to 150 HU. The wide range of HU values of muscle mass can be useful for predicting whole muscle mass area regardless of the fatty infiltration status, but a narrow range of HU values of muscle mass would be more useful for measuring the mass of muscles without fatty changes. Therefore, in our study, we used normal-density muscle mass alone to represent muscle mass without fatty infiltration.
Previous studies evaluated the association between mean PM density and prognosis in patients with various malignancies, such as melanoma, adrenocortical carcinoma, prostate cancer, and biliary or gastrointestinal tract cancer (16–20). Trotter et al. evaluated patients who underwent emergency laparotomy and showed the positive association between sarcopenia, assessed using mean PM density, and prognosis (21). Lindström et al. showed the association between PM quality and prognosis in patients treated for abdominal aortic aneurysms (22). Three studies showed a significant association between PM quality and prognosis in patients who underwent cardiac surgery (23–25). HD patients are at risk of insulin resistance. A previous study calculated the fatty changes of 49 HD patients using mean whole PM values and showed that intramuscular fat infiltration was associated with GS and 6-min walk test performance (26).
Physical performance such as GS is associated with muscle strength and/or mass (27, 28). Physical performance is also associated with various factors beyond muscle mass or strength, such as motor coordination, excitation-contraction coupling, or skeletal integrity. Therefore, decreases in physical performance can be addressed sooner than decreased muscle mass or strength. In addition, decreases in muscle strength would occur sooner than decreases in muscle mass. Although the associations among muscle mass, strength, and physical performance are complex, muscle mass with fatty changes would lead to lower strength and/or physical performance than that without fatty changes despite the same total muscle mass.
In our study, the SGA score as a nutritional index and ASM/Ht2 as a total skeletal muscle mass index were associated with the PM index and NPM index. ASM/Ht2 determined using DXA simply evaluates the total lean mass regardless of fatty changes within the muscle mass. The SGA score is a useful indicator for predicting nutritional status, but it does not discriminate whether the dominantly increased mass is muscle or fat. GS was associated with the NPM index alone. Statistical significance was not obtained, but trend showed a positive association between NPM index and HGS. Strength or physical performance can be influenced by muscle quantity and quality. Our results revealed that the NPM index as an indicator of muscle quality is more closely associated with strength or physical performance. In addition, the correlation coefficients for the SGA score, ASM/Ht2, HGS, and GS were greater for the NPM index than for the PM index.
Our study has inherent limitations, including a cross-sectional design with a single-center setting and a small number of patients. Our study was performed without sample size calculation according to the hypothesis. In addition, our study had a relatively too small a sample size to draw definitive conclusions from. A cross-sectional study with a relatively small sample size, such as ours, would be useful in suggesting issues regarding the associations between the variables rather than drawing definitive conclusions. Second, we did not perform multivariate analyses with adjustment for all confounding factors, such as the presence of exercise, calcium, phosphorus, i-PTH, CRP, or 25-(OH) vitamin D levels. Multivariate analyses using all confounding factors may be optimal for evaluating the independent effects of a variable. However, our study had a relatively small sample size, and we did not perform multivariate analyses adjusting for all variables. Therefore, we selected three covariates (age, sex, and the presence of DM) that can be significantly influenced by outcome measurements and performed multivariate analyses including three covariates and one variable (PM or NPM index). A prospective longitudinal study including a large number of patients and additional indicators is warranted to overcome these limitations.
In conclusion, the present study suggested that the NPM index excluding fatty infiltration may be applicable as an early and useful indicator for detecting muscle strength and physical performance among HD patients. However, limitations regarding the small sample size and the cross-sectional design without sample size calculation according to the hypothesis should be noted. Prospective studies with a larger sample size are needed to identify the definitive associations between the NPM index and muscle strength or physical performance in HD patients.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by IRB of CHA Gumi medical center. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
SK conceptualized and designed the study, performed the data analysis and interpretation, and wrote the manuscript. BK and JK generated and collected the data. JD and JK drafted and revised the manuscript. All authors approved the final version of the manuscript.
Funding
This work was supported by the Medical Research Center Program (2015R1A5A2009124) through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning. The funder had no role in the study design, data collection, analysis, and interpretation, manuscript writing, or decision to submit the article for publication.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. US Renal Data System USRDS 2020 Annual Data Report. Bethesda, MD: Atlas of Chronic Kidney Disease in the United States, National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases (2020).
2. Sabatino A, Cuppari L, Stenvinkel P, Lindholm B, Avesani CM. Sarcopenia in chronic kidney disease: what have we learned so far? J Nephrol. (2020). doi: 10.1007/s40620-020-00840-y. [Epub ahead of print].
3. Giglio J, Kamimura MA, Lamarca F, Rodrigues J, Santin F, Avesani CM. Association of sarcopenia with nutritional parameters, quality of life, hospitalization, and mortality rates of elderly patients on hemodialysis. J Ren Nutr. (2018) 28:197–207. doi: 10.1053/j.jrn.2017.12.003
4. Ikizler TA, Burrowes JD, Byham-Gray LD, Campbell KL, Carrero JJ, Chan W, et al. KDOQI clinical practice guideline for nutrition in CKD: 2020 update. Am J Kidney Dis. (2020) 76(3 Suppl. 1):S1–107. doi: 10.1053/j.ajkd.2020.05.006
5. Kang SH, Lee HS, Lee S, Cho JH, Kim JC. Comparison of muscle mass indices using computed tomography or dual X-ray absorptiometry for predicting physical performance in hemodialysis patients. Kidney Blood Press Res. (2017) 42:1119–27. doi: 10.1159/000485779
6. Kurumisawa S, Kawahito K. The psoas muscle index as a predictor of long-term survival after cardiac surgery for hemodialysis-dependent patients. J Artif Organs. (2019) 22:214–21. doi: 10.1007/s10047-019-01108-4
7. Takata T, Motoe A, Tanida K, Taniguchi S, Ida A, Yamada K, et al. Feasibility of computed tomography-based assessment of skeletal muscle mass in hemodialysis patients. J Nephrol. (2021) 34:465–71. doi: 10.1007/s40620-020-00871-5
8. da Costa JA, Ikizler TA. Inflammation and insulin resistance as novel mechanisms of wasting in chronic dialysis patients. Semin Dial. (2009) 22:652–7. doi: 10.1111/j.1525-139X.2009.00664.x
9. Daugirdas JT. Second generation logarithmic estimates of single-pool variable volume Kt/V: an analysis of error. J Am Soc Nephrol. (1993) 4:1205–13. doi: 10.1681/ASN.V451205
10. Steiber AL, Kalantar-Zadeh K, Secker D, McCarthy M, Sehgal A, McCann L. Subjective Global Assessment in chronic kidney disease: a review. J Ren Nutr. (2004) 14:191–200. doi: 10.1053/j.jrn.2004.08.004
11. Chen LK, Liu LK, Woo J, Assantachai P, Auyeung TW, Bahyah KS, et al. Sarcopenia in Asia: consensus report of the Asian Working Group for Sarcopenia. J Am Med Dir Assoc. (2014) 15:95–101. doi: 10.1016/j.jamda.2013.11.025
12. Long DE, Villasante Tezanos AG, Wise JN, Kern PA, Bamman MM, Peterson CA, et al. A guide for using NIH Image J for single slice cross-sectional area and composition analysis of the thigh from computed tomography. PLoS ONE. (2019) 14:e0211629. doi: 10.1371/journal.pone.0211629
13. Goodpaster BH, Krishnaswami S, Harris TB, Katsiaras A, Kritchevsky SB, Simonsick EM, et al. Obesity, regional body fat distribution, and the metabolic syndrome in older men and women. Arch Intern Med. (2005) 165:777–83. doi: 10.1001/archinte.165.7.777
14. Albu JB, Kovera AJ, Allen L, Wainwright M, Berk E, Raja-Khan N, et al. Independent association of insulin resistance with larger amounts of intermuscular adipose tissue and a greater acute insulin response to glucose in African American than in white nondiabetic women. Am J Clin Nutr. (2005) 82:1210–7. doi: 10.1093/ajcn/82.6.1210
15. Goodpaster BH, Kelley DE, Thaete FL, He J, Ross R. Skeletal muscle attenuation determined by computed tomography is associated with skeletal muscle lipid content. J Appl Physiol. (2000). 89:104–10. doi: 10.1152/jappl.2000.89.1.104
16. Sabel MS, Lee J, Cai S, Englesbe MJ, Holcombe S, Wang S. Sarcopenia as a prognostic factor among patients with stage III melanoma. Ann Surg Oncol. (2011) 18:3579–85. doi: 10.1245/s10434-011-1976-9
17. Miller BS, Ignatoski KM, Daignault S, Lindland C, Doherty M, Gauger PG, et al. Worsening central sarcopenia and increasing intra-abdominal fat correlate with decreased survival in patients with adrenocortical carcinoma. World J Surg. (2012) 36:1509–16. doi: 10.1007/s00268-012-1581-5
18. McDonald AM, Swain TA, Mayhew DL, Cardan RA, Baker CB, Harris DM, et al. CT measures of bone mineral density and muscle mass can be used to predict noncancer death in men with prostate cancer. Radiology. (2017) 282:475–483. doi: 10.1148/radiol.2016160626
19. Chakedis J, Spolverato G, Beal EW, Woelfel I, Bagante F, Merath K, et al. Pre-operative sarcopenia identifies patients at risk for poor survival after resection of biliary tract cancers. J Gastrointest Surg. (2018) 22:1697–708. doi: 10.1007/s11605-018-3802-1
20. Dohzono S, Sasaoka R, Takamatsu K, Hoshino M, Nakamura H. Prognostic value of paravertebral muscle density in patients with spinal metastases from gastrointestinal cancer. Support Care Cancer. (2019) 27:1207–13. doi: 10.1007/s00520-018-4465-x
21. Trotter J, Johnston J, Ng A, Gatt M, MacFie J, McNaught C. Is sarcopenia a useful predictor of outcome in patients after emergency laparotomy? A study using the NELA database. Ann R Coll Surg Engl. (2018) 100:377–81. doi: 10.1308/rcsann.2017.0230
22. Lindström I, Khan N, Vänttinen T, Peltokangas M, Sillanpää N, Oksala N. Psoas muscle area and quality are independent predictors of survival in patients treated for abdominal aortic aneurysms. Ann Vasc Surg. (2019) 56:183–93. doi: 10.1016/j.avsg.2018.08.096
23. Tzeng YH, Wei J, Tsao TP, Lee YT, Lee KC, Liou HR, et al. Computed tomography-determined muscle quality rather than muscle quantity is a better determinant of prolonged hospital length of stay in patients undergoing transcatheter aortic valve implantation. Acad Radiol. (2020) 27:381–8. doi: 10.1016/j.acra.2019.05.007
24. Yamashita M, Kamiya K, Matsunaga A, Kitamura T, Hamazaki N, Matsuzawa R, et al. Prognostic value of psoas muscle area and density in patients who undergo cardiovascular surgery. Can J Cardiol. (2017) 33:1652–59. doi: 10.1016/j.cjca.2017.10.009
25. Foldyna B, Troschel FM, Addison D, Fintelmann FJ, Elmariah S, Furman D, et al. Computed tomography-based fat and muscle characteristics are associated with mortality after transcatheter aortic valve replacement. J Cardiovasc Comput Tomogr. (2018) 12:223–8. doi: 10.1016/j.jcct.2018.03.007
26. Cheema B, Abas H, Smith B, O'Sullivan AJ, Chan M, Patwardhan A, et al. Investigation of skeletal muscle quantity and quality in end-stage renal disease. Nephrology. (2010) 15:454–63. doi: 10.1111/j.1440-1797.2009.01261.x
27. Tietz NW, Shuey DF, Wekstein DR. Laboratory values in fit aging individuals–sexagenarians through centenarians. Clin Chem. (1992) 38:1167–85. doi: 10.1093/clinchem/38.6.1167
Keywords: psoas muscle, hemodialysis, muscle strength, gait speed, fatty change
Citation: Kwack BH, Kim JC, Do JY and Kang SH (2021) Association Between the Normal-Density Psoas Muscle Index and Handgrip Strength or Gait Speed in Maintenance Hemodialysis Patients. Front. Med. 8:712497. doi: 10.3389/fmed.2021.712497
Received: 20 May 2021; Accepted: 05 July 2021;
Published: 28 July 2021.
Edited by:
Ming Yang, Sichuan University, ChinaReviewed by:
Irma Ruslina Defi, Dr. Hasan Sadikin General Hospital, IndonesiaTuyen Van Duong, Taipei Medical University, Taiwan
Copyright © 2021 Kwack, Kim, Do and Kang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Seok Hui Kang, kangkang@ynu.ac.kr
†These authors have contributed equally to this work and share first authorship