Long-term association of pericardial adipose tissue with incident diabetes and prediabetes: the Coronary Artery Risk Development in Young Adults Study
Article information
Abstract
OBJECTIVES
We examined whether pericardial adipose tissue (PAT) is predictive of prediabetes and type 2 diabetes over time.
METHODS
In total, 2,570 adults without prediabetes/diabetes from the Coronary Artery Risk Development in Young Adults Study were followed up over 15 years. PAT volume was measured by computed tomography scans, and the new onset of prediabetes/diabetes was examined 5 years, 10 years, and 15 years after the PAT measurements. Multivariable Cox regression models were used to examine the association between the tertile of PAT and incident prediabetes/diabetes up to 15 years later. The predictive ability of PAT (vs. waist circumference [WC], body mass index [BMI], waist-to-height ratio [WHtR]) for prediabetes/diabetes was examined by comparing the area under the receiver operating characteristic curve (AUC).
RESULTS
The highest tertile of PAT was associated with a 1.56 times (95% confidence interval [CI], 1.03 to 2.34) higher rate of diabetes than the lowest tertile; however, no association was found between the highest tertile of PAT and prediabetes in the fully adjusted models, including additional adjustment for BMI or WC. In the fully adjusted models, the AUCs of WC, BMI, WHtR, and PAT for predicting diabetes were not significantly different, whereas the AUC of WC for predicting prediabetes was higher than that of PAT.
CONCLUSIONS
PAT may be a significant predictor of hyperglycemia, but this association might depend on the effect of BMI or WC. Additional work is warranted to examine whether novel adiposity indicators can suggest advanced and optimal information to supplement the established diagnosis for prediabetes/diabetes.
INTRODUCTION
Obesity is a primary etiology of insulin resistance and diabetes. Specifically, the well-established biological mechanism is that visceral, unlike subcutaneous, adipose tissue secretes free-fatty acids and has higher lipolytic activity, which results in insulin resistance, hyperinsulinemia, and/or diabetes [1-4]. Emerging evidence has also shown that adipose tissue surrounding the heart, which is metabolically akin to and highly correlated with visceral adipose tissue (VAT) [5], may be associated with a higher risk of type 2 diabetes (hereafter “diabetes”) [4,6]. Pericardial adipose tissue (PAT) is an ectopic depot surrounding the coronary arteries, the heart, and the aortic root. PAT secretes inflammatory cytokines into the coronary circulation and more free fatty acids than other adipose tissue depots [5]. As such, PAT is a known risk factor for hyperglycemia, which could potentially have deleterious cardiovascular effects, irrespective of obesity status [5,7].
The association between PAT and diabetes remains inconclusive. A few studies found that the volume of either epicardial adipose tissue (e.g., intra-pericardial adipose tissue) [8] or PAT [4,6] was significantly higher in individuals with prediabetes and diabetes than in normoglycemic individuals, independent of body mass index (BMI), waist circumference (WC), and/or visceral adipose tissue. Conversely, Liu et al. [9] reported that adjustment for VAT attenuated the cross-sectional association between PAT and diabetes. The findings of previous studies are inconsistent in terms of adjustments and limited by cross-sectional study designs. Research with a large sample using a longitudinal study design could clarify the temporality and independence of the associations.
Clinical measures of adiposity, such as WC, BMI, the waist-to-hip ratio, or the waist-to-height ratio (WHtR), have been thoroughly investigated as predictors of diabetes [10-13]. A systematic review and meta-analysis by Kodama et al. [13] concluded that anthropometric indicators of obesity/adiposity (e.g., BMI, WC, the waist-to-hip ratio, WHtR) were statistically significant in predicting diabetes, but WC and WHtR had greater predictive abilities than the other indicators. However, whether PAT predicts long-term incident prediabetes/diabetes and is comparable to the classic anthropometric measures of obesity remains unclear.
The purpose of this study was to examine whether PAT predicts future incident prediabetes/diabetes in individuals without prevalent prediabetes/diabetes. We hypothesized that higher PAT volume is associated with a higher incidence of prediabetes/diabetes over time. Further, we examined variation in the associations by sex, age, race, and/or BMI. Moreover, we aimed to examine the predictive ability of PAT for incident prediabetes/diabetes over 15 years of follow-up in comparison to the predictive abilities of WC, BMI, and WHtR. We hypothesized that the predictive ability of PAT in predicting incident prediabetes/diabetes is comparable to that of WC, BMI, and/or WHtR. Lastly, as an exploratory aim, we examined an optimal cut-off value of PAT for predicting incident prediabetes/diabetes.
MATERIALS AND METHODS
Study sample
The Coronary Artery Risk Development in Young Adults (CARDIA) study is a prospective cohort study investigating the trends and correlates of coronary artery disease risk factors and their clinical sequelae. The CARDIA study was initiated in 1985-1986 (year 0), with a total of 5,115 healthy Black and White male and female, aged 18 years to 30 years, who completed an in-person clinical examination and were enrolled from 4 United States metropolitan communities: Birmingham, AL; Minneapolis, MN; Chicago, IL; and Oakland, CA. In the present study, we set CARDIA year 15 (2000-2001) as the baseline because computed tomography (CT) scans were first performed in year 15; therefore, PAT data have been available since 2000-2001. From here on, to avoid any confusion on timeline between ours and other CARDIA papers, we will refer to our baseline as “year 15.” From the 3,671 participants (72% of year 0) who participated in the CARDIA year 15 exam (2000-2001), we additionally excluded participants who were missing data on PAT (n=778), had diabetes or prediabetes (n=125), or were missing data on covariates (n=198) at exam year 15. Thus, a total of 2,570 participants were included in the analysis (Supplementary Material 1).
Assessment of pericardial adipose tissue
The volume (cm3) of PAT at year 15 was quantified on 2.5-mm or 3.0-mm-thick scans of the chest obtained using 64-channel, ECG-gated, multidetector CT scanners (GE Healthcare, Milwaukee, WI at the Birmingham and Oakland centers; Siemens, Erlangen, Germany at the Chicago and Minneapolis centers) [6]. To quantify PAT volume, CT slices 15 mm above and 30 mm below the superior extent of the left main coronary artery were selected using an image processing workstation (OsiriX, Pixmeo, Geneva, Switzerland). The analysts manually segmented the images in this 45-mm block of images and then applied a threshold of -190 to -30 Hounsfield units to isolate adipose tissue. Reader reproducibility was calculated in 156 randomly selected re-reads (intrareader: 49 pairs; inter-reader: 107 pairs). The intra-reader variability was 2.0% and the inter-reader variability was 4.2% [6].
Assessment of prediabetes/diabetes
Diabetes status was defined based on a combination of (1) self-report of oral hypoglycemic medication or insulin use (at years 0, 7, 10, 15, 20, 25, and 30); (2) measured fasting glucose ≥ 126 mg/dL (at years 0, 7, 10, 15, 20, 25, and 30); (3) a glycated hemoglobin (HbA1c) ≥ 6.5% (at years 20 and 25); or (4) a 2-hour postload glucose ≥ 200 mg/dL during a 75-g oral glucose tolerance test (at years 10, 20, 25, and 30), and prediabetes was defined based on a combination of (1) measured fasting glucose 100-125 mg/dL; (2) a 2-hour postload glucose 140-199 mg/dL; or (3) HbA1c 5.7-6.4% [14,15]. Fasting blood samples drawn by venipuncture were processed at a central laboratory using a standard protocol (Molecular Epidemiology and Biomarker Research Laboratory, University of Minnesota, Minneapolis, MN, USA and Northwest Lipid Research Center, Seattle, WA, USA).
Other assessments
Sex, age (years), study center, race (Black, White), education (years), occupation status (full time; yes/no), alcohol consumption (mL/day), family history of diabetes (either mother or father; years 0, 5, 10, and 25 only), antihypertensive and blood lipid-lowering medication use, and menopausal status and hormone use (among female only) were assessed using a standard questionnaire. Smoking status (never, former, current) was assessed using a tobacco use questionnaire [16]. BMI was calculated as the measured weight in kilograms divided by height in meters squared. WC was measured to the nearest 0.5 cm in duplicate and averaged using anthropometric tape midway between the iliac crest and the lowest lateral portion of the rib cage. WHtR was calculated as WC divided by height in meters. Systolic and diastolic blood pressure were measured 3 times after 5 minutes of rest (1-minute rest between the last 2 readings) using a Hawksley random-zero sphygmomanometer (Hawksley, Lancing, UK) [17], with the first and fifth phase Korotkoff sounds recorded and the second and third measurements averaged. Diet quality was assessed using the A Priori Diet Quality Score at years 0, 7, and/or 20 only [18]. Diet quality scores were summed across the 10 components for a maximum score of 100, with low values reflecting a poor diet and high values reflecting a healthy diet. If a diet quality score was available from 1 year (i.e., year 0, 7, or 20), the score from the available year was used; otherwise, the average scores from 2 exam or 3 exam years were calculated to account for variation in diet across time and used for analysis [19]. Moderate-to-vigorous intensity physical activity (exercise units) was assessed using the interviewer-administered CARDIA Physical Activity History Questionnaire [20]. Plasma concentrations of total cholesterol and high-density lipoprotein-cholesterol (HDL-C) were determined enzymatically using a Hitachi 917 analyzer (Hitachi, Tokyo, Japan).
Statistical analysis
Descriptive statistics, including frequency distributions, as well as measures of central tendency and variability, were used to calculate and present the characteristics of participants from years 15 to 30 stratified by tertile of PAT at year 15. The differences between tertiles of PAT were examined using analysis of variance, the Kruskal-Wallis test, or the chi-square test as appropriate.
Person-time (years) was calculated from the number of participants multiplied by 5 years, 10 years, or 15 years because the observation cycle of the outcome variables was 5-year intervals (beginning at year 15) for all participants. The number of participants with prediabetes/diabetes per 1,000 person-years from year 15 by PAT tertile was calculated, and then the incidence rate was calculated as the total number of new cases of prediabetes/diabetes divided by the sum of the person-years. Multivariable-adjusted Cox proportional hazard regression models (PROC PHREG in SAS version 9.4) were used to estimate the hazard ratios (HRs) and 95% confidence intervals (CIs) for incident prediabetes/diabetes 5 years to 15 years later across PAT tertiles and a continuous PAT variable (per 10 cm3 increment) at year 15. The covariates were as follows. Model 1 adjusted for sex, race, center, age at year 15, education (years), and occupation status at year 30; and model 2 adjusted for the variables in model 1, plus smoking status, averages (exam years 15, 20, 25, and 30) of physical activity, alcohol consumption (binary; 0 vs. ≥ 0 mL/day), systolic/diastolic blood pressure, total cholesterol, and HDL-C, and diet quality score derived from years 0, 7, or 20, family history of diabetes at year 25, and antihypertensive and lipid-lowering medication use at year 15. Model 3 adjusted for the variables in model 2, plus average values (exam years 15, 20, 25, and 30) of BMI. We adjusted for clinically relevant covariates a priori based on knowledge from previous studies that they may confound the associations between adipose tissue and diabetes. Potential effect modification was examined using multiplicative interactions of sex, age (median split: 32-40, 41-51 years at year 15), race (Black, White), BMI (< 25.0, 25.0-29.9, ≥30.0 kg/m2), and WC (<102 vs. 102 cm ≤for male; <88 vs. 88 cm ≤ for female) at year 15 with PAT.
Unadjusted and multivariable-adjusted (model 2 of the Cox regression models) receiver operating characteristic (ROC) curve analysis was used to examine the diagnostic ability of PAT at year 15 for predicting prediabetes/diabetes 5 years to 15 years later by estimating sensitivity, specificity, and area under the ROC curve (AUC). To examine statistical differences in AUC for PAT versus WC, BMI, and WHtR for predicting prediabetes/diabetes, the algorithm developed by DeLong et al. [21] was used. In addition, separate fully-adjusted ROC curve models were analyzed by WC, BMI, and WHtR to examine whether the addition of PAT changed the AUC of each model. For reference, Spearman correlation coefficients between PAT, WC, BMI, and WHtR from year 15 were examined (Supplementary Material 2). To estimate an optimal cut-off value of PAT in predicting prediabetes/diabetes, we examined the balance of sensitivity and specificity with the highest Youden index [22].
In sensitivity analyses, first, we adjusted for WC instead of BMI in the fully adjusted Cox regression models to examine whether PAT predicts incident prediabetes/diabetes independent of a proxy measure of VAT. Second, we examined the association of PAT with hyperglycemia (prediabetes and diabetes). Third, we examined the association of PAT with prediabetes/diabetes defined by fasting glucose only across exam years. Fourth, we examined the association of PAT (continuous variable) with fasting glucose 5 years, 10 years, and 15 years later, separately, using multivariable linear regression. Lastly, we examined the association of PAT using quartile and quintile splits with incident prediabetes/diabetes. All statistical analyses were conducted using SAS version 9.4 (SAS Institute Inc., Cary, NC, USA). A 2-sided p-value < 0.05 was considered statistically significant.
Ethics statement
All participants provided written informed consent, and the CARDIA study protocols were annually approved by the institutional review board at each CARDIA center (Birmingham, AL; Minneapolis, MN; Chicago, IL; and Oakland, CA; USA).
RESULTS
The characteristics of the study cohort at year 15 across tertiles of PAT are presented in Table 1. The mean age of the 2,570 eligible participants was 40.3 years, 55.8% were female, and 44.2% were Black. Except for the rate of alcohol consumption and family history of diabetes, all variables were significantly different across the tertiles of PAT (all p<0.05). The correlations of PAT with BMI, WC, and WHtR were 0.59, 0.74, and 0.69, respectively (all p<0.001; Supplementary Material 2). Compared to the participants in the analytic sample, excluded participants were more likely to be Black, younger, and current smokers, reported lower education levels, had higher WC, BMI, WHtR, systolic/diastolic blood pressure, and PAT, and had lower HDL-C levels, diet quality scores, and physical activity scores (all p<0.05; Supplementary Material 3).
The adjusted HRs and 95% CIs of incident prediabetes/diabetes in 5 years to 15 years across tertiles of PAT at year 15 are presented in Table 2. The incidence rates of diabetes and prediabetes were higher in higher tertiles of PAT. In the fully adjusted model, the highest tertile of PAT was 1.56 times (95% CI, 1.03 to 2.34) more likely to have incident diabetes than the lowest tertile of PAT. Each 10-cm3 higher PAT volume was associated with a 10.14 times (95% CI, 10.08 to 10.20) higher incidence of diabetes in the fully adjusted model. There were no significant interactions of sex, age, and BMI with the tertile of PAT in the association with incident prediabetes/diabetes. A race interaction in the association with incident prediabetes was found; however, none of the associations were significant in the stratified analyses.
The predictive abilities (AUCs) of PAT, WC, BMI, and WHtR in predicting incident prediabetes/diabetes 5 years to 15 years later are presented in Figure 1. In the unadjusted model predicting diabetes, the AUCs of WC, BMI, and WHtR were significantly higher than that of PAT (p<0.05); however, none were significantly different in the fully adjusted model. The AUCs of WC in predicting incident prediabetes were significantly higher than that of PAT (p<0.05) in both unadjusted and fully adjusted models (the AUC of WHtR was significantly higher than that of PAT).
The AUCs of WC, BMI, and WHtR, after adjustment for covariates and with and without adjustment for PAT, in predicting incident prediabetes/diabetes at 5 years and 15 years are presented in Figure 2. The AUCs of WC, BMI, and WHtR in predicting incident diabetes improved after adjusting for PAT; however, none were statistically significant. Furthermore, none of the AUCs for predicting incident prediabetes significantly improved after adjusting for PAT.
The highest sensitivity, specificity, and Youden index to estimate the best cut-off values of PAT for predicting prediabetes/diabetes were 34.5 cm3/42.3 cm3 for female and 51.2 cm3/63.1 cm3 for male, respectively (Supplementary Material 4). In the fully adjusted model, the association between the highest tertile of PAT with incident diabetes was not significant after adjusting for WC instead of BMI (data not shown). The associations of PAT tertile with incident hyperglycemia (Supplementary Material 5) and fasting glucose-defined prediabetes/diabetes (Supplementary Material 6) were not significant in the fully adjusted models. PAT was significantly associated with high fasting glucose 10 years and 15 years later in the fully adjusted models (Supplementary Material 7). Similar to our findings using tertile splits for classifying PAT, the highest quartile (Supplementary Material 8) and the highest quintile (Supplementary Material 9) of PAT showed significant positive associations with incident diabetes in the fully adjusted models. Lastly, none of the associations of PAT tertiles with incident prediabetes/diabetes were significant for BMI (Supplementary Material 10; except for tertile 2 and incident prediabetes in the BMI category < 25 kg/m2) or WC (Supplementary Material 11) stratifications.
DISCUSSION
In this longitudinal cohort study, high PAT was associated with a higher incidence of diabetes and prediabetes 5 years to 15 years later, independent of the tested covariates; however, these associations were attenuated after additionally adjusting for BMI (for prediabetes) or WC (for diabetes). In addition, the predictive ability of PAT for incident diabetes was comparable to that of WC, BMI, and WHtR when accounting for all the tested covariates. However, the predictive ability of PAT for incident prediabetes was lower than that of WC.
Our findings on the prospective associations of PAT with incident prediabetes/diabetes fill a gap in the current literature that has been limited by cross-sectional observations on prevalent diabetes or higher fasting glucose [6,8,9]. Similar to previous cross-sectional findings, we found that high PAT was associated with a higher incidence of prediabetes/diabetes 5 years to 15 years later. However, after accounting for different adiposity variables in the association of PAT with diabetes, results were mixed across previous studies. For example, Alman et al. [6] found that the highest quartile of PAT was associated with a higher prevalence of diabetes after adjusting for all covariates, including BMI (odds ratio [OR], 2.57; 95% CI, 1.66 to 3.98) or VAT (OR, 2.08; 95% CI, 1.32 to 3.29). Furthermore, Wang et al. [8] found that the highest tertile of epicardial adipose tissue was associated with a higher prevalence of diabetes (OR, 4.82; 95% CI, 1.55 to 16.58) after accounting for all covariates, including WC. On the contrary, Liu et al. [9] found that PAT as a continuous variable was associated with a higher prevalence of diabetes (OR, 1.40; 95% CI, 1.01 to 1.94); however, this association was attenuated after adjusting for BMI or VAT. We found that both the highest tertile of PAT and PAT as a continuous variable were associated with a higher incidence of diabetes after adjusting for BMI (despite a null association when adjusting for WC), and adjusting for BMI attenuated the association between the highest tertile of PAT and prediabetes. Alman et al. [6] suggested that there might be a threshold effect of the highest quartile of PAT in the association with a higher prevalence of diabetes (null findings in quartiles 2 and 3 of PAT). Similarly, our results also appeared to show a threshold effect in the highest tertile of PAT for predicting diabetes (null finding in tertile 2, quartiles 1-3, and quintiles 1-4 of PAT; as shown in Table 3, Supplementary Materials 8 and 9). However, continuous PAT was also associated with both a high risk of incident diabetes and elevated fasting glucose. PAT may be associated with incident diabetes independent of confounders, including BMI; however, WC, as a proxy of visceral adiposity, may confound the association.
Additional longitudinal research is warranted because it is scientifically premature to directly compare our longitudinal results to the previous cross-sectional findings [6,8,9], especially with heterogeneity in measurement matrices (e.g., thickness) and with different approaches to PAT categorization. Bosy-Westphal et al. [23] examined both cross-sectional and longitudinal associations of PAT with insulin resistance and sensitivity in overweight female during approximately 14 weeks of a calorie-restricted diet intervention. PAT was correlated with insulin resistance at baseline (r= 0.46); however, after the 14-week short-term intervention, the longitudinal association between PAT and insulin resistance was not significant.
Although there is no clear consensus as to which adiposity indicators best predict diabetes, WC and WHtR have been recognized as stronger indicators than BMI [10,24,25]. We observed that the AUCs of WC, BMI, and WHtR seemed to be superior to that of PAT for predicting diabetes in crude models. However, interestingly, the AUC of PAT increased and became comparable to those of BMI, WC, and WHtR after adjusting for potential confounders. Furthermore, adding PAT into the BMI and fasting glucose diagnostic models for predicting incident diabetes improved the AUCs (although not to a statistically significant extent). PAT increases progressively with age [26] and is dependent on major health conditions and behaviors [5,27,28]. This suggests that PAT may also be a good adiposity measure, when accounting for socio-demographics and key health behaviors and characteristics, to understand the future risk of diabetes, compared to the wellknown indicators of diabetes [10,24,25]. However, we found that WC was a better indicator than PAT for incident prediabetes, even after adjusting for covariates.
The optimal cut-off values of PAT volume at exam year 15 for predicting incident prediabetes/diabetes 5 years to 15 years later were 34.5 cm3/42.3 cm3 for female and 51.2 cm3/63.1 cm3 for male, respectively. A small body of evidence has examined the optimal cut-off values of epicardial adipose tissue thickness in predicting cardiovascular events [29,30] or metabolic syndrome [31,32]. However, to our knowledge, there is no evidence suggesting optimal cut-off values of PAT in predicting prediabetes/diabetes. Additional research with more diverse study populations is needed to confirm our preliminary findings.
The health implications of prediabetes have been underestimated; however, prediabetes is a critical indicator of diabetes risk, and approximately 38% of United States adults ≥ 18 years have prediabetes [33]. As such, early and precise diagnostic actions for identifying prediabetes are important to prevent long-term complications [34]. Overall, we found that the magnitude of associations between PAT and hyperglycemia, and the predictive ability of PAT for hyperglycemia varied by prediabetes/diabetes status. This is possibly because the incidence rate of prediabetes is much higher than that of diabetes, and there may be pathophysiological differences between these conditions [34]. However, there is a lack of evidence about the differences in the predictive ability and the patterns of the association between diagnostic indicators (including PAT) and prediabetes, compared to the associations of those indicators with diabetes.
Strengths of this study include the large sample size with a high participant retention rate, long-term observations of study variables, thorough examinations of potential confounding and effect modifications, a sex/race-balanced cohort, standardized CT measurement of PAT, and repeated screening of glycemic status to define prediabetes/diabetes. However, our study has several limitations. First, the findings cannot be generalized to wider populations given that the study participants were midlife adults residing in the United States. Second, glycemic status was measured every 5 years; therefore, it was not possible to assess the exact timing of the new cases of prediabetes/diabetes. Third, the criteria to define prediabetes/diabetes varied across exam years and 2-hour glucose and HbA1c were not measured at year 15; therefore, it is likely that prediabetes/diabetes might have been underestimated. However, the findings were similar when prediabetes/diabetes was defined with the use of fasting glucose only. Fourth, not only WHtR but also the waist-to-hip ratio may be a good indicator [24,25]. A hip measure was not conducted at years 15 and 20 in CARDIA; therefore, we could not compare the predictive ability of the waist-to-hip ratio. However, WHtR may be superior to the waist-to-hip ratio in predicting prediabetes/diabetes [35]. Lastly, there may have been unmeasured residual confounders that could have interacted with the associations.
In conclusion, a high level of PAT was associated with a higher risk of incident prediabetes/diabetes 5 years to 15 years later, independent of major health characteristics, but this association was attenuated after accounting for BMI or WC. In addition, PAT was comparable to other known anthropometric measures, such as BMI, WC, or WHtR, for predicting incident diabetes. The significance of our findings suggests that PAT can be used as a diagnostic indicator to identify individuals who have prediabetes/diabetes or may be at risk for future diabetes. However, we also suggest that the association of PAT with prediabetes/diabetes may be dependent on BMI or WC. Furthermore, there are health implications of examining the association between PAT and prediabetes/diabetes in that individuals with diabetes are at high risk of cardiovascular events [36] and atherosclerosis [37], resulting from the adverse health consequences of excessive PAT deposition [38].
SUPPLEMENTARY MATERIALS
Supplementary materials are available at https://www.e-epih.org/.
Notes
CONFLICT OF INTEREST
The authors have no conflicts of interest to declare for this study.
FUNDING
The Coronary Artery Risk Development in Young Adults Study (CARDIA) is conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with the University of Alabama at Birmingham (HHSN268201800005I & HHSN268201800007I), Northwestern University (HHSN268201800003I), University of Minnesota (HHSN268201800006I), and Kaiser Foundation Research Institute (HHSN268201800004I).
AUTHOR CONTRIBUTIONS
Conceptualization: Oh M, Whitaker KM, Kim JY. Data curation: Oh M, Schreiner PJ. Formal analysis: Oh M, Cho W, Lee DH, Kim JY. Funding acquisition: Schreiner PJ. Methodology: Oh M, Cho W, Lee DH, Whitaker KM, Schreiner PJ, Kim JY. Project administration: Oh M, Kim JY. Visualization: Oh M, Cho W, Lee DH, Whitaker KM, Schreiner PJ, Terry JG, Kim JY. Writing – original draft: Oh M. Writing – review & editing: Oh M, Cho W, Lee DH, Whitaker KM, Schreiner PJ, Terry JG, Kim JY.
Acknowledgements
This manuscript has been reviewed by CARDIA for scientific content.