Higher Muscle Mass Protects Women with Gestational Diabetes Mellitus from Progression to Type 2 Diabetes Mellitus
Article information
Abstract
Background
We evaluated whether postpartum muscle mass affects the risk of type 2 diabetes mellitus (T2DM) in Korean women with gestational diabetes mellitus (GDM).
Methods
A total of 305 women with GDM (mean age, 34.9 years) was prospectively evaluated for incident prediabetes and T2DM from 2 months after delivery and annually thereafter. Appendicular skeletal muscle mass (ASM) was assessed with bioelectrical impedance analysis at the initial postpartum visit, and ASM, either divided by body mass index (BMI) or squared height, and the absolute ASM were used as muscle mass indices. The risk of incident prediabetes and T2DM was assessed according to tertiles of these indices using a logistic regression model.
Results
After a mean follow-up duration of 3.3 years, the highest ASM/BMI tertile group had a 61% lower risk of incident prediabetes and T2DM compared to the lowest tertile group, and this remained significant after we adjusted for covariates (adjusted odds ratio, 0.37; 95% confidence interval [CI], 0.15 to 0.92; P=0.032). Equivalent findings were observed in normal weight women (BMI <23 kg/m2), but this association was not significant for overweight women (BMI ≥23 kg/m2). Absolute ASM or ASM/height2 was not associated with the risk of postpartum T2DM.
Conclusion
A higher muscle mass, as defined by the ASM/BMI index, was associated with a lower risk of postpartum prediabetes and T2DM in Korean women with GDM.
INTRODUCTION
Women with gestational diabetes mellitus (GDM) have an increased risk of developing type 2 diabetes mellitus (T2DM) compared to women with normoglycemic pregnancy [1,2]. Several risk factors for postpartum T2DM have been identified which include maternal age, body mass index (BMI), family history of diabetes, and the degree of hyperglycemia during pregnancy [1,2]. These risk factors also contribute to postpartum diabetes in the Korean population [3-5], but their attributable risk seems to vary according to ethnic groups [6-8].
Asians develop diabetes mellitus at a lower BMI than Westerners [9,10]. A previous multiethnic study identified that Asian women had a more than three times higher prevalence of GDM in the same BMI range compared to non-Hispanic white or African American women [11]. Furthermore, following GDM, Asian women had a higher prevalence of postpartum T2DM and metabolic syndrome compared to Europeans [7]. For a given BMI, there are differences in the body composition of fat and muscle mass across the ethnic groups: Asians have more body fat at a lower BMI compared to Caucasians [12], and adiposity is well known to elevate the risk of T2DM [13]. Recent reports have suggested that low muscle mass can increase insulin resistance and diabetes risk [14,15]. Since skeletal muscle is a metabolically active organ that accounts for approximately 30% of the resting metabolic rate, reduced skeletal muscle mass can impair glucose utilization and decrease reservoir capacity [16,17]. A 9-year prospective cohort study in Korea demonstrated that middle-aged and older adults with the lowest muscle mass tertile had a more than 2-fold increased risk of T2DM after adjusting for covariates, including obesity parameters [14]. Therefore, we hypothesized that assessment of muscle mass could provide important information regarding diabetic risk especially in Asian women.
Here, we set out to investigate if low muscle mass contributes to an increased incidence of prediabetes and T2DM in women with GDM, with the goal of offering additional information beyond BMI.
METHODS
Study design
This prospective observational study of women with GDM was conducted in Seoul National University Bundang Hospital. All subjects participated voluntarily, and informed consent was obtained from each subject. The Institutional Review Board of Seoul National University Bundang Hospital approved the study protocol (IRB Number: B-0503/018-001) and retrospective data analysis for this study (IRB Number: B-1909/565-106). This study was conducted according to the Declaration of Helsinki.
Subjects were recruited from 2009 to 2018 and followed until the end of 2019. The initial postpartum follow-up visit was performed at 6 to 8 weeks postpartum, and annual follow-up visits were made thereafter to evaluate the postpartum maternal glycemic status. Subjects were censored at the onset of diabetes or the last visit of examination (or when lost to follow-up). The mean follow-up duration was 3.3 years (Supplementary Fig. 1). A total of 378 women with GDM were screened for eligibility; 305 women were eligible to participate according to the inclusion and exclusion criteria listed below. GDM was diagnosed with a 100-g oral glucose tolerance test (OGTT) at 24 to 28 weeks of gestation according to the criteria recommended by the American Diabetes Association [18]. A diagnosis of GDM was made if two or more of the following were present: fasting plasma glucose ≥95 mg/dL, postprandial (PP) 1 hour ≥180 mg/dL, PP 2 hours ≥155 mg/dL, and/or PP 3 hours ≥140 mg/dL. Subjects with a history of type 1 diabetes mellitus or T2DM, recurrent GDM, or T2DM at postpartum 2 months were excluded, and only newly diagnosed GDM subjects were included. We excluded those with an underlying medical condition, including cardiovascular disease, thyroid disease, or any condition that required concurrent medication (e.g., steroids, lipid-lowering, or anti-hypertensive agents).
Body composition analysis
Body composition was assessed with a direct segmental eight-point multifrequency bioelectrical impedance analysis (BIA, InBody770, Biospace Inc., Seoul, Korea) at the initial postpartum evaluation. This analysis was previously shown to have acceptable agreement with dual-energy X-ray absorptiometry (DXA) and magnetic resonance imaging [19,20]. Total fat mass and visceral fat area (VFA) were measured by BIA. VFA measured by BIA was previously reported to correlate well with that measured by computed tomography (CT), especially in Korean women with age <50 years and BMI <25 kg/m2 [21]. We employed a formula developed to approximate the VFA measured by BIA (VFA-BIA) to that measured by CT (VFACT) [21]. Estimated VFA, which was used for further analysis, is calculated as: (–142.77)+(1.40×VFA-BIA)+(–1.29×age)+(–0.98×BMI)+(2.14×waist circumference)+(–0.03×VFABIA×BMI)+(0.07×age×BMI).
The primary muscle quantity was calculated as the appendicular skeletal muscle mass (ASM) divided by BMI, and participants were assorted into tertiles as follows: lowest, 0.51–0.66; middle, 0.67–0.73; and highest, 0.74–0.94. As secondary muscle mass measures, we examined the absolute ASM and the ASM normalized with the squared height (ASM/Ht2). These were also divided into tertiles, as follows: for ASM, lowest 11.3–15.1, middle 15.2–16.7, and highest 16.8–23.5 (kg); and for ASM/Ht2, lowest 4.60–5.93, middle 5.94–6.40, and highest 6.41-8.33 (kg/m2).
Anthropometric and clinical parameters
Antepartum parameters were collected at 24 to 28 weeks of gestation; these included glucose levels at OGTT and the levels of glycosylated hemoglobin (HbA1c), insulin, and C-peptide. The parameters evaluated at initial postpartum follow-up, which occurred at 6 to 8 weeks after delivery, are presented in Table 1; they included baseline anthropometric, clinical, and biochemical measurements such as maternal age, height, body weight, waist circumference, ASM, fat mass, blood pressure, and lipid profiles. Systolic and diastolic blood pressures were measured with subjects in a seated position using an electronic sphygmomanometer (UA-1020 device, A&D Co., Tokyo, Japan). We assessed physical activity after delivery. The physical activity was classified into mild, moderate, and vigorous-intensity, and we defined exercise as >150 minutes per week of moderate-intensity activity or >75 minutes per week of vigorous-intensity activity following the World Health Organization (WHO) recommendations [22]. Postpartum BMI change was defined as the difference between the BMI at the initial postpartum follow-up and that at the last visit of examination.
Metabolic assessment
Postpartum maternal glycemic status was evaluated with a 75-g OGTT given at 6 to 8 weeks after delivery and annually thereafter. Participants were classified into the following groups: normoglycemia; prediabetes, including impaired glucose or impaired glucose tolerance; and T2DM based on American Diabetes Association criteria [18]. HbA1c was measured using a Bio-Rad Variant II Turbo HPLC analyzer (BioRad, Hercules, CA, USA). Glucose concentrations were analyzed using the hexokinase method. Plasma concentrations of insulin and C-peptide were measured by radioimmunoassay (IZOTOP, HoilBioMed Co., Seoul, Korea). The insulinogenic index was calculated from OGTT parameters as follows: insulin (30 min)–insulin (0 min)/glucose (30 min)–glucose (0 min). Homeostasis model assessment of insulin resistance (HOMA-IR) and β-cell function (HOMA-β) were calculated [23]. Triglyceride concentrations were measured by the glycerol-3-phosphate oxidase peroxide method, and total cholesterol, high-density lipoprotein, and low-density lipoprotein were measured with a Hitachi 747 chemistry analyzer (Hitachi, Tokyo, Japan).
Statistical analysis
All data are expressed as the mean±standard deviation or as numbers and percentages. Analysis of variance (ANOVA) was used to compare baseline characteristics according to ASM/BMI tertiles and post hoc analysis was performed with the Bonferroni method among different groups. The incidence of prediabetes and T2DM was compared using the chi-square test and statistical significance was assessed using linear by linear association.
The risk of incident ‘prediabetes and T2DM’ was assessed using a binary logistic regression model, and ‘prediabetes’ and ‘T2DM’ were assessed using a multinomial logistic regression model, according to the tertiles of ASM/BMI. Three models were fitted; they sequentially included the covariates from the initial postpartum follow-up, as follows. Model 1 was unadjusted; Model 2 was adjusted for maternal age, family history of diabetes, systolic blood pressure, total fat mass, and waist circumference; and Model 3 was adjusted for the variables listed in Model 2 plus breast feeding, exercise after delivery, postpartum BMI change, and the antepartum fasting plasma glucose concentration. The lowest tertile was set as the reference category. Unadjusted and adjusted odds ratios (ORs) for incident prediabetes and T2DM were calculated. The participants were stratified by their initial postpartum BMI (normal weight group, BMI <23 kg/m2; overweight group, BMI ≥23 kg/m2) for sensitivity analysis. The secondary muscle mass measures (absolute ASM and ASM/Ht2) were examined in separate models to verify the extents of association with incident prediabetes and diabetes and assess the level of consistency across the different muscle indices. All statistical analyses were performed using SPSS version 20.0 (IBM Corp., Armonk, NY, USA).
RESULTS
Clinical characteristics of the participants at the antepartum and postpartum periods
A total of 305 women (mean age, 34.9±3.7 years) was prospectively evaluated for incident prediabetes and T2DM, with the mean follow-up duration of 3.3 years. At the 6 to 8 weeks postpartum visit, the mean BMI was 23.1±3.5 kg/m2. The ASM/BMI tertiles comprised 103, 111, and 91 women with GDM from the lowest to the highest tertiles, respectively. As shown in Table 1, the highest tertile group had lower antepartum and postpartum HbA1c as well as lower baseline postpartum BMI, waist circumference, body fat mass, blood pressures, and levels of fasting plasma glucose, insulin, C-peptide, and triglycerides, compared to the lowest tertile group.
Incidence of prediabetes and T2DM according to ASM/BMI
The composite incidence of prediabetes and T2DM was significantly lower in the highest ASM/BMI group, being found at 57.3%, 55.9%, and 36.3% in the lowest, middle, and highest groups, respectively (P for trend=0.004) (Table 2, Supplementary Table 1). After we adjusted for covariates such as age, family history of diabetes, blood pressure, total fat mass, waist circumference, antepartum fasting glucose concentration, breast feeding, exercise after delivery, and postpartum BMI change, the highest tertile had a 63% lower risk of developing prediabetes and T2DM (OR, 0.37; 95% confidence interval [CI], 0.15 to 0.92; P=0.032). The separate incidence of prediabetes also showed a significant association, with the risk being significantly lower in women of the highest tertile (OR, 0.36; 95% CI, 0.14 to 0.94; P=0.036). The incidence of T2DM was 72% lower in the highest group in univariate analysis (OR, 0.28; 95% CI, 0.10 to 0.82; P=0.020), but this association was attenuated after we adjusted for confounders (OR, 0.44; 95% CI, 0.06 to 30.7; P=0.437).
We assessed the contribution of visceral fat, which is a known risk factor for T2DM [24]. Subjects with more visceral fat had a higher prevalence of postpartum prediabetes and T2DM (55.3% in the upper half vs. 45.8% in the lower half). Next, we performed multivariate logistic analysis including estimated VFA as a covariate. After we adjusted for VFA and other covariates, subjects in the highest tertile of ASM/BMI still had lower risks of developing postpartum prediabetes and T2DM (OR, 0.32; 95% CI, 0.14 to 0.72; P=0.007) (Supplementary Table 2).
The risk of incident prediabetes and T2DM according to ASM/BMI after stratification by BMI
For further analysis, we stratified the participants by BMI (normal weight, BMI <23 kg/m2 and overweight BMI ≥23 kg/m2) (Table 3). In the normal weight group, prediabetes and T2DM occurred at rates of 51.6%, 48.4%, and 32.9% in the lowest, middle, and highest tertiles, respectively. In a multivariate adjusted model, the highest muscle group had a 79% lower risk for developing prediabetes and T2DM (OR, 0.21; 95% CI, 0.06 to 0.77; P=0.019). However, this association was not significant for the overweight group (OR, 0.56; 95% CI, 0.10 to 3.26; P=0.522).
Postpartum changes of HOMA-IR and HOMA-β according to ASM/BMI
There was no significant difference in HOMA-IR across the tertiles during the antepartum period (Table 4). However, postpartum insulin resistance was lowest in the highest ASM/BMI group. HOMA-IR significantly decreased in the middle and highest tertiles from the antepartum to postpartum periods, but no statistically significant decline was seen for HOMA-IR in the lowest tertile.
Within each tertile, HOMA-β was decreased significantly from the antepartum to postpartum periods. However, there was no statistically significant difference in the antepartum to postpartum change in HOMA-β or postpartum insulinogenic index among the different tertiles.
Incident prediabetes and T2DM according to the absolute ASM and ASM/Ht2
Next, we evaluated the association of absolute ASM and ASM normalized by height (ASM/Ht2) with incident prediabetes and T2DM according to the absolute ASM tertiles (Supplementary Tables 3 and 4). Of note, the highest ASM tertile group was older and had higher BMI, waist circumference, body fat mass, blood pressures, and levels of fasting plasma glucose, insulin, C-peptide, and triglycerides at baseline (Supplementary Table 3). Antepartum fasting plasma glucose and HbA1c were higher for the highest tertile compared to the lowest tertile. These metabolic phenotypes therefore showed trends opposite to those observed for ASM/BMI.
In a crude analysis, the absolute ASM was not associated with incident prediabetes and T2DM (Supplementary Table 3). No significant associations between the absolute ASM with incident prediabetes and T2DM were found after further adjusting for covariates. In the normal weight subgroup (BMI <23 kg/m2), intriguingly, the highest ASM group had a lower risk of incident prediabetes and T2DM in a univariate analysis (OR, 0.20; 95% CI, 0.06 to 0.74; P=0.015), and this association showed marginal significance after we adjusted for covariates (OR, 0.21; 95% CI, 0.04 to 1.16; P=0.072). No such trend was observed in women with overweight or obese (BMI >23 kg/m2).
The trends in baseline characteristics by ASM/Ht2 were similar to those seen for the absolute ASM groups: Women in the highest ASM/Ht2 tertile were older and had higher BMI, waist circumference, ASM, body fat mass, blood pressures, and levels of fasting plasma glucose, HbA1c, insulin, C-peptide, and triglycerides at baseline (Supplementary Table 4). Antepartum fasting plasma glucose and HbA1c were also increased for the highest group compared to the lowest group. Intriguingly, the incidence of prediabetes and T2DM was 46.5%, 43.0%, and 62.2% in the lowest, middle, and highest groups, respectively, which was opposite the relationship observed for ASM/BMI.
In univariate analysis, the highest ASM/Ht2 group had an approximately 2-fold increased risk for incident prediabetes and T2DM (OR, 1.95; 95% CI, 1.10 to 3.44; P=0.022), but the association was lost after we adjusted for covariates (OR, 2.00; 95% CI, 0.73 to 5.48; P=0.176). In addition, we did not observe any significant association between ASM/Ht2 and incident prediabetes and T2DM after we adjusted for covariates following BMI stratification.
DISCUSSION
This study investigated the association between postpartum muscle mass and the risk of prediabetes and T2DM in women with a previous history of GDM. Higher muscle mass defined by the ASM/BMI index was associated with lower insulin resistance after delivery and a reduced risk of prediabetes and T2DM. In a sensitivity analysis, this association was found specifically in the normal weight group (BMI <23 kg/m2). We compared different muscle indices, including ASM/BMI, absolute ASM, and ASM/Ht2, and found that maternal body parameters (e.g., body weight) highly affected the risk of postpartum prediabetes and T2DM. Our results emphasize the need to use an appropriate muscle mass index to predict metabolic outcomes in postpartum women. From our findings, we conclude that evaluation of muscle mass and adopting ASM/BMI in the postpartum period is useful for predicting future T2DM in women who had GDM, especially in a non-obese population.
GDM is characterized by underlying maternal defects in the β-cell response to insulin resistance during pregnancy [25]. We previously reported that several genetic variants related to GDM, including CDKN2A/2B, HHEX, and CDKAL1, are associated with β-cell function and an increased risk for postpartum T2DM [26,27]. Diabetes in Asians is characterized by an underlying functional β-cell deficit that occurs years before the development of T2DM, with such individuals failing to show a proper compensatory increase in insulin secretion upon insulin resistance [28]. This partly explains the higher incidence of GDM and postpartum T2DM in Asians compared to other ethnicities at a given BMI [11,29].
In addition to these non-modifiable factors, there are modifiable risk factors that contribute to postpartum T2DM; they include prepregnancy BMI and postpartum weight change, waist circumference, exercise, and lactation [5,27,30-32], all of which are associated with insulin resistance. Our previous study showed that women who lost weight had improvements in insulin sensitivity and the lowest incidence of T2DM, while women who gained weight showed the opposite associations during 4 years of follow-up [4]. Given that skeletal muscle is the metabolic organ that utilizes most of the glucose in the body, it is important to evaluate the association between muscle mass and progression to T2DM in GDM women.
Muscle mass has been widely studied in association with insulin resistance and the incidence of T2DM in middle-aged and elderly populations. Low muscle mass is associated with an increased risk of T2DM and higher insulin resistance in both Koreans and Caucasians [14,15]. However, the effect of muscle mass on T2DM has been less well studied in young populations. The present study included young women at their peak of muscle mass, which begins to decrease after age 40 [33], and who were free from other comorbidities. The significant improvement in insulin sensitivity after delivery in the highest muscle mass group is consistent with our hypothesis that more muscle mass may confer a protective effect after GDM. We demonstrate that, even for a young population, relatively low muscle mass manifests as a risk for T2DM.
Higher muscle mass is expected to improve insulin sensitivity, but women with polycystic ovary syndrome (PCOS) typically exhibit increased insulin resistance while having more muscle mass compared to women without PCOS [34,35]. To better understand the relationship between muscle mass and the risk of T2DM, we need to clarify the factors that contribute to increased muscle mass in women with PCOS (e.g., androgens, higher body weight, and/or fat mass). A recent meta-analysis reported that higher total and trunk lean tissue mass were observed in a subgroup of women with PCOS who are BMI ≥25 kg/m2, which were attributed to overweight and obesity, but not androgen or insulin resistance [36]. On the other hand, increases in total or trunk lean mass were not evident in women with BMI <25 kg/m2. This suggests that increased body weight, rather than PCOS itself, may contribute to the higher muscle mass observed in women with PCOS. In the present study, which did not exclude women with PCOS, high ASM/BMI was shown to be protective for postpartum prediabetes and T2DM after we adjusted for total fat or body weight. Nevertheless, the possible contribution of PCOS to muscle mass or insulin resistance, as well as its linking mechanisms, should be studied further.
Asian populations develop T2DM at lower BMI than Western populations [9,10]. A previous multiethnic study including GDM women of European, Asian, and African ethnicities found an increased prevalence of impaired glucose tolerance and metabolic syndrome in Asians compared to Europeans at postpartum 20 months [7]. In addition, the relative contributions of conventional risk factors for GDM differ by ethnicity as advanced maternal age or overweight/obesity affects Chinese to a lesser extent than other ethnic groups [6]. Therefore, it is important to distinguish clinical risk factors for postpartum diabetes in Asian GDM women with relatively lower BMI, as these factors may differ from those in other ethnicities.
Here, we stratified subjects by BMI 23 kg/m2, and found that the normal weight women (BMI <23 kg/m2) in the highest ASM/BMI group developed less prediabetes and T2DM, whereas no association or risk reduction was found in the overweight group (BMI ≥23 kg/m2). As in women with PCOS, overweight women have more muscle mass, but this excess muscle mass is considered metabolically inactive because these women are insulin resistant. Therefore, the muscle mass in overweight women was not associated with the development of postpartum glucose intolerance. The results of a previous Korean nationwide study support our finding that subjects with higher muscle mass showed improved insulin sensitivity only among lean (low fat) subjects, but this protective association of muscle mass and insulin sensitivity was attenuated in subjects with higher proportions of fat [37]. This suggests that the contribution of muscle mass to glucose metabolism may be more important in women with relatively lower BMI, but further validation studies including women with other ethnicities are required.
Since Baumgartner first proposed the height adjusted skeletal muscle mass as a measure for representing muscle quantity [38], a number of working groups have adopted the ASM/Ht2 as a standard muscle mass index [39,40]. Meanwhile, the Foundation for the National Institutes of Health proposed ASM adjusted with BMI as having the best clinical implications for immobility and mortality [41,42]. ASM/BMI is better correlated with the age pattern in Korean women than other indices [33], but debate remains regarding the best measure for muscle mass with respect to clinical implications. Here, we compared different muscle indices to better understand the contribution of muscle mass to the development of postpartum diabetes in association with maternal body weight or BMI.
In the present study, higher ASM/BMI was found to be a protective predictor for the composite postpartum prediabetes and T2DM. On the other hand, no association was found between the absolute ASM or ASM/Ht2 and T2DM, suggesting that there is a wide divergence between the different muscle indices. Similarly discordant findings were reported in a recent cohort study, in which absolute ASM was positively associated with incident T2DM [43]. Thus, future studies including other ethnicities are required to investigate the disparities observed in how different muscle indices are associated with metabolic outcomes, including T2DM.
This study has several limitations. First, we used the BIA method in assessing the body composition. This may not fully reflect the actual skeletal muscle mass. However, it is a non-invasive and valid technique that provides accurate estimates of skeletal muscle mass that are closely correlated with measurements obtained using DXA and magnetic resonance imaging across ranges of age, volume status, and BMI [19,20]. Assessment of muscle mass as well as VFA using magnetic resonance imaging or CT would provide additional insights into predicting metabolic outcomes according to body composition. Second, we were unable to evaluate the contribution of muscle strength or quality to the development of postpartum glucose intolerance. The clinical implications of muscle strength have been thoroughly assessed in older adults, but future studies are needed to clarify the link between muscle strength and metabolic risk in younger populations. Third, the dietary pattern, which has been associated with the risk of postpartum diabetes, was not assessed in this study. As an alternative, we included postpartum BMI change as a covariate to reflect nutritional status, since maternal BMI has been reported to attenuate the beneficial effect of a healthy dietary pattern postpartum [44]. Fourth, due to the observational nature of this study, the follow-up duration was not uniform among the study subjects and we cannot confirm causal relationships or exclude confounding effects. The sample size was small, especially for the sensitivity analyses, which limited further assessment.
In conclusion, the muscle mass defined by the ASM/BMI index was associated with incident prediabetes and T2DM in Korean women with GDM, independent of previously known risk factors including maternal age, glycemia during pregnancy, body fat mass, postpartum BMI change, breastfeeding, and exercise. When seeking to quantify the overall diabetes risk in women after GDM pregnancy, clinicians may find that assessing muscle quantity will help identify women at high risk for developing T2DM. Pre- and post-pregnancy aerobic and resistance exercise may help prevent postpartum diabetes by increasing muscle mass and reducing body fat, which should be further investigated through intervention studies.
SUPPLEMENTARY MATERIALS
Supplementary materials related to this article can be found online at https://doi.org/10.4093/dmj.2021.0334.
Notes
CONFLICTS OF INTEREST
Tae Jung Oh was editorial board member of the Diabetes & Metabolism Journal from 2020 to 2021. Sung Hee Choi has been editorial board member of the Diabetes & Metabolism Journal since 2021. They were not involved in the review process of this article. Otherwise, there was no conflict of interest.
AUTHOR CONTRIBUTIONS
Conception or design: Y.S., H.C.J
Acquisition, analysis, or interpretation of data: Y.S., J.H.M., T.J.O., H.C.J.
Drafting the work or revising: Y.S., J.H.M., H.C.J.
Final approval of the manuscript: Y.S., J.H.M., T.J.O., C.H.A., J.H.M., S.H.C., H.C.J
FUNDING
This research was supported by the grant from Seoul National University Bundang Hospital to Hak Chul Jang (grant No. 02-2019-015).
Acknowledgements
The authors thank Division of Statistics in Medical Research Collaborating Center at Seoul National University Bundang Hospital for statistical supports.