Renal Tubular Damage Marker, Urinary N-acetyl-β-D-Glucosaminidase, as a Predictive Marker of Hepatic Fibrosis in Type 2 Diabetes Mellitus
Article information
Abstract
Background
Non-alcoholic steatohepatitis is closely associated with the progression of diabetic kidney disease (DKD) in type 2 diabetes mellitus (T2DM). We investigated whether urinary N-acetyl-β-D-glucosaminidase (u-NAG), an early renal tubular damage biomarker in DKD, could be related to the degree of hepatic fibrosis in patients with T2DM.
Methods
A total of 300 patients with T2DM were enrolled in this study. Hepatic steatosis and fibrosis were determined using transient elastography. The levels of urinary biomarkers, including u-NAG, albumin, protein, and creatinine, and glucometabolic parameters were measured.
Results
Based on the median value of the u-NAG to creatinine ratio (u-NCR), subjects were divided into low and high u-NCR groups. The high u-NCR group showed a significantly longer duration of diabetes, worsened hyperglycemia, and a more enhanced hepatic fibrosis index. A higher u-NCR was associated with a greater odds ratio for the risk of higher hepatic fibrosis stage (F2: odds ratio, 1.99; 95% confidence interval [CI], 1.04 to 3.82). Also, u-NCR was an independent predictive marker for more advanced hepatic fibrosis, even after adjusting for several confounding factors (β=1.58, P<0.01).
Conclusion
The elevation of u-NAG was independently associated with a higher degree of hepatic fibrosis in patients with T2DM. Considering the common metabolic milieu of renal and hepatic fibrosis in T2DM, the potential use of u-NAG as an effective urinary biomarker reflecting hepatic fibrosis in T2DM needs to be validated in the future.
INTRODUCTION
Type 2 diabetes mellitus (T2DM) is associated with an increased risk of nonalcoholic fatty liver disease (NAFLD) [1]. The prevalence of NAFLD is estimated to be 45% to 75% in patients with T2DM, more than two times higher than in the general population [2,3]. Also, the overall prevalence of biopsy-proven nonalcoholic steatohepatitis (NASH) and advanced fibrosis in patients with T2DM is approximately 37.33% and 17.02%, respectively [3]. The coexistence of T2DM and NAFLD significantly increases the risk for the development of NASH and cirrhosis compared to that observed in cases of NAFLD without T2DM [1]. Diabetic kidney disease (DKD) is prevalent in approximately 30% to 50% of the population with T2DM [4], and NAFLD is associated with increased incidence and prevalence of nephropathy in T2DM, a phenomenon that is independent of traditional risk factors [5]. This strong relationship between NAFLD and DKD in T2DM may be mediated by the common pathogenic metabolic mechanisms underlying NAFLD and DKD, but the systemic release of several pro-inflammatory mediators in NAFLD may also putatively contribute to their relationship [6].
Due to the high prevalence of NAFLD and NASH in T2DM patients, an effective and non-invasive biomarker to predict the degree of NAFLD and NASH in T2DM is urgently needed. However, the prevalence of histologically proven NASH in patients with T2DM who exhibit normal liver enzyme levels is 56% [7], suggesting that serum liver enzyme levels are not sufficient to predict the presence of steatohepatitis and fibrosis in patients with T2DM and NAFLD. Considering the close association between NAFLD and DKD, it would be reasonable to investigate whether a urinary biomarker that originally reflects the presence of DKD could simultaneously reflect the development of NASH.
Urinary N-acetyl-β-D-glucosaminidase (u-NAG), a well-established renal tubular damage marker, is a potential candidate. NAG is present in the lysosomes of proximal tubule epithelial cells, and its urinary excretion is increased in response to renal proximal tubular cell injury [8]. Since profibrotic changes in renal tubular cells are the cardinal features associated with diabetic nephropathy [8,9], u-NAG levels may reflect the degree of renal fibrosis. Considering that advanced liver fibrosis is independently associated with chronic kidney disease (CKD) among individuals with NAFLD [10], u-NAG, as a potential renal fibrosis marker, may be associated with the degree of NASH and vice versa. Furthermore, the enzymatic activities of serum NAG and u-NAG were elevated in cirrhotic and fibrotic liver disorders, indicating that NAG can be used as a potential systemic biomarker that can identify fibrosis in both the liver and kidney [11,12].
In this study, we investigated, using transient liver elastography in patients with T2DM, whether u-NAG, an early biomarker of renal tubulopathy, is associated with hepatic fibrosis.
METHODS
Study population
In a retrospective cross-sectional study, we enrolled patients aged ≥19 years with T2DM who had been simultaneously evaluated for transient elastography, FibroScan (Echosens, Paris, France), as well as blood and urinary markers, including u-NAG, albumin, protein, and creatinine at the Severance Hospital between April 2017 and July 2020. T2DM patients were defined as (1) patients taking insulin or an oral hypoglycemic or (2) patients with glycosylated hemoglobin (HbA1c) ≥6.5% (47.5 mmol/mol) according to the Korean Diabetes Association clinical practice guidelines [13]. Patients were excluded if they presented any of the following: age <19 years, type 1 diabetes mellitus, pregnancy, hepatic diseases including various forms of hepatitis, alcoholic liver disease and hepatocellular carcinoma, end-stage renal disease (defined as estimated glomerular filtration rate [eGFR] <15 mL/min/1.73 m2), and renal replacement therapy, including renal transplantation and dialysis.
Demographic data such as age, sex, weight, height, waist circumference, blood pressure, duration of diabetes, and information about current medications were obtained. Body mass index (BMI) was calculated as body weight divided by height squared (kg/m2). Hypertension was defined as a systolic blood pressure of ≥130 mm Hg or diastolic blood pressure of ≥85 mm Hg, or the current use of antihypertensive medications according to the guidelines proposed in 2005 by the American Heart Association and the National Heart, Lung, and Blood Institute for Asian populations [14]. The study protocol was reviewed and approved by the Institutional Review Board of Yonsei University College of Medicine, Seoul, Republic of Korea (approval No. 4-2020-0490); it waived the need for informed consent because of the retrospective nature of the study and the lack of personal information.
Measurement of blood and urinary parameters
Following an overnight fast (≥8 hours), morning spot urine samples were collected for measuring u-NAG, albumin, protein, and creatinine levels, and blood samples were collected to measure complete blood count, fasting and postprandial glucose, lipid profile, platelet count, and levels of insulin, C-peptide, HbA1c, glycated albumin, creatinine, albumin, and liver enzymes. Low-density lipoprotein cholesterol (LDL-C) levels were calculated using the Friedewald equation in the absence of actual LDL-C measurement [15]. The equation was not used if the patient’s triglyceride level exceeded 400 mg/dL. Pancreatic β-cell function and insulin sensitivity were assessed in non-insulin users using the following indices: homeostasis model assessment-β-cell (HOMA-β)=[20×fasting insulin (μIU/mL)]/[fasting plasma glucose (mmol/L)–3.5] and homeostasis model assessment-insulin resistance (HOMA-IR)= [fasting insulin (μIU/mL)×fasting plasma glucose (mmol)]/22.5 [16]. The eGFR was calculated using the Chronic Kidney Disease Epidemiology Collaboration equation [17]. U-NAG, albumin, and protein levels were adjusted as u-NAG to creatinine ratio (u-NCR), urinary albumin to creatinine ratio (u-ACR), and urinary protein to creatinine ratio (u-PCR), respectively, to minimize the influence of kidney function. Proteinuria was defined as 150 mg/g creatinine, based on the most conservative reported normal value for urinary protein excretion (<150 mg/day) [18]. Albuminuria was defined as 30 mg/g creatinine according to the Kidney Disease: Improving Global Outcomes recommendations [19]. CKD was defined as CKD stage 3 (eGFR <60 mL/min/1.73 m2). A u-NCR was categorized into two groups: at or below median u-NCR (≤6.77 U/g creatinine) and above the median u-NCR (>6.77 U/g creatinine) as low u-NCR and high u-NCR groups, respectively.
Measurement of hepatic steatosis and fibrosis
Transient elastography, FibroScan, was performed by a welltrained physician who was blinded to the participants’ clinical details; the final value was obtained using a standard procedure. Hepatic steatosis was classified into the following stages based on a controlled attenuation parameter (CAP) value (decibels per meter [dB/m]): steatosis stage 0 (S0) <238 dB/m, 238 dB/m≤ steatosis stage 1 (S1) <260 dB/m, 260 dB/m≤ steatosis stage 2 (S2) <293 dB/m, and steatosis stage 3 (S3) ≥293 dB/m [20]. The stages of hepatic fibrosis were further divided based on the liver stiffness measurement (LSM) value (kPa): no fibrosis and mild fibrosis (F0–1) <7 kPa, significant fibrosis (F2) ≥7 kPa, advanced fibrosis (F3) ≥8.7 kPa, and cirrhosis (F4) ≥11.5 kPa [21].
Using a previously validated fatty liver prediction model, NAFLD was defined as a liver fat score (NLFS) of >–0.64 and Framingham steatosis index (FSI). The fibrosis-4 calculator (FIB-4 index) and NALFD fibrosis score (NFS) were used to determine the extent of hepatic fibrosis, and significant fibrosis was defined as either FIB-4 >2.67, or NFS >0.676, as previously validated [2].
Statistical analysis
The data are presented as mean±standard deviation for normally distributed continuous variables and as medians (interquartile range) for non-normally distributed continuous variables. Categorical data are expressed as numbers and percentages. The baseline characteristics of the study population were analyzed based on the u-NCR level using the two-sample Student’s t-test for continuous variables and the chi-square test for categorical variables. A normality test was performed for all the continuous variables. To correct the skewed distributions, non-normally distributed continuous variables, such as u-NCR, u-ACR, and u-PCR were logarithmically transformed before statistical analysis. Outliers of u-NCR ≥23.275 U/g creatinine (n=25) were excluded from the data using Tukey’s fences method [22].
Correlations between u-NAG and other parameters were analyzed using Pearson’s correlation analysis. Multiple logistic regression analyses were used to determine the independent association between u-NAG and the hepatic fibrosis stages. Receiver operating characteristic (ROC) analysis with Youden’s index was conducted to calculate the cut-off value of u-NCR to predict moderate and advanced hepatic fibrosis (stage ≥F2 and F3) [23].
Stepwise multiple linear regression analyses were performed to evaluate the relationship between LSM and various parameters using u-NCR and u-ACR. Statistical analyses were performed using IBM SPSS statistical software for Windows version 25.0 (IBM, Armonk, NY, USA). Adjusted odds ratio (aOR) and 95% confidence interval (CI) were determined. Statistical significance was set at P<0.05.
RESULTS
Baseline characteristics of the study population according to urinary NAG
A total of 300 patients (204 men and 96 women) were enrolled in this study. The mean age and median duration of diabetes were 57.1 and 3 years, respectively. The median value of u-NCR was 6.77 U/g creatinine (U/g Cr). The patients were divided into low and high u-NCR groups based on the median u-NCR value (Table 1). The high u-NCR group was characterized by the following: a higher percentage of men, a significantly longer duration of diabetes, and a higher prevalence of hypertension and CKD. Blood pressure, hepatic enzymes, lipid profile except for total cholesterol, and the use of lipid medication did not differ between the two groups. The use of antidiabetic medication, including sulfonylurea and insulin, and antihypertensive medications was more common in the high NCR group than in the low NCR group (Supplementary Table 1). The high NCR group exhibited significantly higher levels of fasting and postprandial glucose, HbA1c, and glycated albumin (all P<0.001), and a significantly lower HOMA-β level (P=0.009). HOMA-IR scores did not differ between the two groups (P=0.097). In the high NCR group, levels of u-ACR, albuminuria, and u-PCR were higher, but eGFR was lower than that in the low NCR group. The higher NCR group had a significantly higher degree of fibrosis as determined by LSM (5.4% vs. 6.3%, P=0.001) and NFS (–1.11 vs. –0.63, P=0.032). Also, more advanced stages of hepatic fibrosis were observed in the high NCR group. However, no significant difference in hepatic steatosis determined by CAP, NLFS, and FSI was observed between the two groups.
Correlations between urinary NAG and other parameters
The u-NCR level was positively and significantly correlated with the duration of diabetes and the levels of fasting and postprandial glucose, HbA1c, and glycated albumin (Table 2). Also, the values of u-ACR and u-PCR were positively correlated with that of u-NCR. However, HOMA-β and eGFR were negatively correlated with u-NCR. With respect to hepatic fibrosis indices, the FIB-4 index (r=0.34, P<0.001) and NFS (r=0.35, P<0.001) were positively correlated with u-NCR. However, the hepatic steatosis indices, NLFS and FSI, were not significantly correlated with u-NCR. BMI, waist circumference, and lipid profile, except for total cholesterol, were not associated with u-NCR.
Risk of hepatic steatosis/fibrosis stages according to urinary NAG
The association between u-NAG and the stages of hepatic steatosis and fibrosis in relation to potential confounders was analyzed using logistic regression analysis. Insulin non-users (n=249) were analyzed to include HOMA indices because of their significance in previous correlation results and clinical settings (Supplementary Tables 2 and 3). Even after adjusting for multiple confounding parameters including demographic factors, the consumption of hypertension and lipid medication, metabolic parameters, and glycemic parameters in Model 4, a higher u-NCR was associated with a greater odds ratio for the risk of higher stage hepatic fibrosis, as the fibrosis stage increased (F2 [OR, 1.99; 95% CI, 1.04 to 3.82; P=0.039], F3 and F4 [OR, 2.40; 95% CI, 1.52 to 3.80; P<0.001]) (Fig. 1B, Supplementary Table 3). However, no significant association was observed between u-NCR and any stage of hepatic steatosis (Fig. 1A, Supplementary Table 2). When the median value of u-NCR in each hepatic fibrosis stage was compared (Fig. 2A), u-NCR in the F4 stage showed the highest value across the fibrosis stages compared with those observed in the F0 and F1 stages (9.4 [5.04 to 14.3] vs. 5.9 [3.8 to 9.3], P=0.001); u-ACR did not exhibit a significant trend across the stages of hepatic fibrosis (Fig. 2B). Similar results were obtained when the total population including insulin users (n=300) was analyzed (Supplementary Tables 4 and 5). In this analysis, HOMA indices were excluded owing to the inclusion of insulin users.
The cut-off values of u-NCR showed the best equilibrium between sensitivity and specificity at 6.31 and 7.84 U/g creatinine to predict the hepatic fibrosis stage above F2 and F3, respectively. The area under the ROC curve values using each cut-off value for F2 and F3 were 0.654 and 0.651, respectively.
Multiple linear regression analysis of factors affecting the liver stiffness measurement
To determine the predictive factors for hepatic fibrosis, multiple linear regression analysis was performed with non-insulin users (n=249) (Table 3). Age, sex, BMI, duration of diabetes, history of hypertension, use of lipid medication, HOMA-IR, HOMA-β, u-NCR, and levels of triglyceride, alanine aminotransferase (ALT), and HbA1c, were included as independent determinants. In Model 1, when age, sex, BMI, and duration of diabetes were adjusted, u-NCR was identified as a significant independent factor for LSM (standardized β [STD β]=1.58, P<0.001). Hypertension, lipid medication, and metabolic parameters such as triglyceride and ALT levels and HOMA-IR were included in Models 2 and 3, showing that u-NCR was independently associated with LSM. After Model 4 was further adjusted for glycemic factors, including HbA1c level and HOMA-β, u-NCR remained a significant coefficient value for a higher degree of hepatic fibrosis across the models (STD β=1.58, P<0.001). Similar results were obtained when the total population including insulin users (n=300) was analyzed (Supplementary Table 6). In this analysis, HOMA indices were excluded owing to the inclusion of insulin users. The efficacy of u-ACR as a predictive marker for hepatic fibrosis was also investigated and compared with that of u-NCR (Supplementary Tables 7 and 8). U-ACR was not associated with LSM in Model 4 in both non-insulin users and the total population, including insulin users (STD β=0.21, P=0.33 and STD β=0.17, P=0.43).
DISCUSSION
Accumulating evidence suggests that the development and progression of NAFLD in T2DM could be associated with DKD characterized by glomerular and tubular damage [5,6,24]. Among the various urinary proteins involved in DKD, urinary albumin primarily constitutes proteinuria, but various non-albumin proteins are also pathophysiologically filtered or secreted. Among the renal tubular secretory proteins, NAG is a well-established renal tubular damage marker [25]. We assumed that systemic inflammation and oxidative stress in the metabolic milieu of T2DM may simultaneously induce both hepatic parenchymal and renal tubular fibrosis and that u-NAG may also be proportionally elevated due to hepatic fibrosis in response to renal fibrosis. Based on this hypothesis, we investigated the relationship between renal tubular damage and hepatic fibrosis in patients with T2DM by measuring u-NAG and liver fibrosis using transient elastography and demonstrated that the increase in u-NAG was significantly associated with a higher degree of hepatic fibrosis in patients with T2DM (Fig. 3). This study has three main findings. First, an increased level of u-NAG is associated with nephropathic indices, poor glycemic control, and the use of antihypertensive medication. Second, even before the development of overt albuminuria, increased u-NAG level is associated with an increased risk for a higher degree of hepatic fibrosis in patients with T2DM. Third, the degree of hepatic steatosis is not associated with the levels of u-NAG and albuminuria in patients with T2DM.
Serum and urinary levels of NAG may be distinguished due to their different origins. NAG, a lysosomal enzyme, is located in the intracellular lesions of several cells, including renal tubular cells [26]. In addition to the urine, NAG is also present in serum, and serum NAG is mainly secreted from the liver [26]. Therefore, increased serum NAG levels may be associated with hepatic damage, including fibrosis. In the present study, increased urinary excretion of NAG may have directly resulted from damage to renal tubular cells rather than an increased glomerular filtration of serum NAG considering its relatively high molecular weight (>130 kDa) [27]. However, despite their different sources, the urinary concentration of NAG may be closely associated with the serum concentration potentially due to the common metabolic and systemic milieu [28,29]. For example, both serum and urinary concentrations of NAG increased during endotoxemia in vivo [28], possibly by the simultaneous damage to the liver and kidney. Also, serum and urinary concentrations of NAG were commonly elevated in patients with type 1 diabetes mellitus and microangiopathy [29]. In the present study, the high u-NCR group had a longer duration of diabetes and worsened hyperglycemia compared to the low u-NCR group. Thus, both kidney and liver of the subjects from the high u-NCR group may have been chronically exposed to systemic metabolic stress, and these common pathophysiological mechanisms may have resulted in a simultaneous elevation of NAG in both urine and serum. As the serum NAG concentration was not measured, we could not demonstrate a correlation between serum and u-NAG concentrations. Further studies with a complete measurement of serum and u-NAG should be performed to validate their significant associations.
An association between non-albumin proteinuria (NAP) and indices of hepatic fibrosis was previously tested and showed a significant correlation [30]. NAP is mainly related to renal tubular damage [31], and it encompasses a variety of proteins such as α-1 microglobulin, β-2 macroglobulin, cystatin C, kidney injury molecule-1, and NAG [31,32]. Therefore, a definite protein that contributed to the level of NAP in association with hepatic fibrosis could not be identified in the previous study [30]. We selected a specific urinary biomarker, NAG, among several components of NAP, to investigate its association with hepatic fibrosis. To our knowledge, this is the first study to analyze the association between this specific renal tubular marker, NAG, and the degree of hepatic fibrosis. Furthermore, the current data make a distinct point compared with the previous study [30], given the use of a different modality to determine NAFLD. Instead of using hepatic fibrosis indices [30], NAFLD was determined using FibroScan, one of the most widely used and well-validated non-invasive methods to assess hepatic fibrosis [33]. The diagnostic accuracy of FibroScan for hepatic fibrosis is superior to that of fibrosis indices [34]. In the present study, we supported our thesis based on compelling evidence generated by FibroScan.
In another previous report [12], the urinary activity of NAG was elevated in patients with liver cirrhosis in Child-Pugh class C compared with class A and B. However, our data differ from the previous data due to the different characteristics of the study subjects. Patients with diverse liver etiologies but with the fibrosis stage restricted to liver cirrhosis were included in the previous study [12]. In our study, patients with a varying spectrum of hepatic fibrosis attributed to NAFLD were enrolled. Owing to a unified liver etiology with NALFD and a diverse range of hepatic fibrosis, a quantitative association between urinary NAG and hepatic fibrosis in a comprehensive range of NAFLD could be investigated in this study. As a result, a higher u-NAG level was associated with an increased risk of hepatic fibrosis, particularly advanced fibrosis. In contrast, u-NAG was not associated with hepatic steatosis determined using indices or transient elastography. In short, u-NAG may be associated with the development and progression of hepatic fibrosis rather than hepatic steatosis. Our data suggests that a concrete urinary biomarker, NAG, could be used to predict and detect NAFLD in reversible stages. Although further validation studies should be conducted, the adoption of our findings in the clinical field would contribute to a decrease in NAFLD-related mortality and morbidity by early detection and monitoring of progress during NAFLD treatment.
U-NAG constitutes isolated NAP which is related to renal tubulointerstitial fibrosis [35]. In a previous study [36], u-NAG level was proportional to urinary NAP levels increased in T2DM patients with DKD even before onset of overt albuminuria. Therefore, u-NAG along with NAP could be used as early biomarkers of DKD in patients with T2DM in normoalbuminuric stage. Early DKD with renal tubulopathy is aggravated in systemic response to uncontrolled hyperglycemia, resulting in pro-inflammatory changes and oxidative stress in the whole body [25,30,36]. Considering the localized effect of hyperglycemia on the kidney, increased glucose reabsorption in the proximal tubule causes tubulointerstitial hypoxia, oxidative stress, and fibrosis as well [37,38]. In the current study, increased u-NAG levels were observed in patients with poorly controlled diabetes, β-cell dysfunction, and long-standing diabetes. Thus, chronically enhanced systemic inflammation due to hyperglycemia and insulin depletion in patients increased NAG levels may have resulted in increased pro-inflammatory cytokine production and oxidative stress that consequently induced hepatic fibrosis [30,39]. According to the two-hit theory of NASH, the progression to hepatic fibrosis involves triglyceride deposition in hepatocytes, followed by the formation of reactive oxygen species and free radicals [30,39]. Such development of liver fibrosis may have caused systemic release of proinflammatory mediators that augmented u-NAG elevation in return. Given the correlation of u-NAG with glucometabolic parameters, including HbA1c, HOMA-β, and duration of diabetes that reflect the deterioration of hyperglycemia-associated chronic systemic inflammation, we postulate that u-NAG might have the potential to be a novel biomarker effectively reflecting both renal and hepatic fibrosis. Further investigations with a prospective design for detailed mechanisms about the association between u-NAG levels and fibrotic burden in the kidney and liver are warranted in the future.
Concerning the association between urinary albumin and NAFLD, albuminuria is related to insulin resistance that underlies the mechanism of metabolic syndrome and plays a role in the development of NAFLD [30,40]. However, in this study, the level of albuminuria was not associated with the degree of hepatic steatosis or fibrosis. These results could be explained by the characteristics of study population. The study subjects had normal kidney function with a mean eGFR of 92.1 mL/min/1.73 m2. Also, the median u-ACR was 12.5 mg/g creatinine (5.7 to 45.7), and albuminuria was prevalent in 32.3% of the total subjects. Thus, the study population had a normal or relatively minimal albuminuria that might have mitigated the association between u-ACR and NAFLD. Meanwhile, u-NAG still significantly reflected the degree of hepatic fibrosis even when albuminuria had not been fully developed in the present study.
The current study has several limitations. First, a causal relationship between increased u-NAG and hepatic fibrosis could not be determined due to the cross-sectional study design. Long-term, large prospective studies are needed to validate the relationship between elevated u-NAG levels and hepatic fibrosis in the future. Second, liver biopsy data were absent in the current study, although it is the gold standard for diagnosing and predicting the progression of NAFLD by detecting inflammation and the ballooning degeneration of hepatocytes [41,42]. The lack of biopsy data in the current study poses a potential pitfall for underestimating the severity and progression of NAFLD [42]. Third, we could not investigate the self-reported history of alcohol consumption, and patients with a history of alcohol-related liver disease were identified using the International Classification of Diseases Clinical Modification code. Thus, confounding by alcohol consumption remains when interpreting the analyses for the association between u-NAG and hepatic fibrosis.
To our knowledge, this is the first study demonstrating that increased u-NAG levels are independently associated with an increased risk for a higher degree of hepatic fibrosis in patients with T2DM, even before the development of albuminuria. Enhanced systemic pro-inflammatory response and oxidative stress caused by uncontrolled hyperglycemia may have created a common milieu that explains the quantitative association between u-NAG levels and the degree of hepatic fibrosis. The potential use of u-NAG as an effective clinical biomarker of hepatic fibrosis in T2DM needs to be validated in the future.
Supplementary Materials
Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2020.0273.
Notes
CONFLICTS OF INTEREST
No potential conflict of interest relevant to this article was reported.
AUTHOR CONTRIBUTIONS
Conception or design: B.W.L
Acquisition, analysis, or interpretation of data: H.K.K., M.L., B.W.L.
Drafting the work or revising: H.K.K, M.L., Y.L., E.S.K., B.S.C., B.W.L.
Final approval of the manuscript: B.W.L.
FUNDING
This work was supported by the Technology Innovation Program funded by the Ministry of Trade, Industry & Energy (MOTIE, Korea) (20004103) and by Severance Hospital Research Fund for Clinical excellence (SHRC C 2019-0032).
Acknowledgements
This study was supported by the ‘SENTINEL (Severance ENdocrinology daTa scIeNcE pLatform)’ program funded by the 2020 Research fund of Department of Internal Medicine, Severance Hospital, Seoul, Korea and Sung-Kil Lim Research Award (4-2018-1215; DUCD000002) in statistical analyses.
The authors thank Medical Illustration & Design, part of the Medical Research Support Services of Yonsei University College of Medicine, for all artistic support related to this work.