- Open Access
Prediction of Methionine and Homocysteine levels in Zucker diabetic fatty (ZDF) rats as a T2DM animal model after consumption of a Methionine-rich diet
Nutrition & Metabolism volume 15, Article number: 14 (2018)
Although alterations in the methionine metabolism cycle (MMC) have been associated with vascular complications of diabetes, there have not been consistent results about the levels of methionine and homocysteine in type 2 diabetes mellitus (T2DM). The aim of the current study was to predict changes in plasma methionine and homocysteine concentrations after simulated consumption of methionine-rich foods, following the development of a mathematical model for MMC in Zucker Diabetic Fatty (ZDF) rats, as a representative T2DM animal model.
The model building and simulation were performed using NONMEM® (ver. 7.3.0) assisted by Perl-Speaks-NONMEM (PsN, ver. 4.3.0). Model parameters were derived using first-order conditional estimation method with interactions permitted among the parameters (FOCE-INTER). NCA was conducted using Phoenix (ver. 6.4.0). For all tests, we considered a P-value < 0.05 to reflect statistical significance.
Our model featured seven compartments that considered all parts of the cycle by applying non-linear mixed effects model. Conversion of S-adenosyl-L-homocysteine (SAH) to homocysteine increased and the metabolism of homocysteine was reduced under diabetic conditions, and consequently homocysteine accumulated in the elimination phase.
Using our model, we performed simulations to compare the changes in plasma methionine and homocysteine concentrations between ZDF and normal rats, by multiple administrations of the methionine-rich diet of 1 mmol/kg, daily for 60 days. The levels of methionine and homocysteine were elevated approximately two- and three-fold, respectively, in ZDF rats, while there were no changes observed in the normal control rats.
These results can be interpreted to mean that both methionine and homocysteine will accumulate in patients with T2DM, who regularly consume high-methionine foods.
The methionine metabolism cycle (MMC) plays important roles in regulating both methionine and cysteine homeostasis , and is related with functions in terms of cell growth and development under normal conditions [2, 3]. As a universal key intermediate in the MMC, homocysteine is not obtained from the diet, but is remethylated to methionine, or converted to cysteine by the trans-sulfuration pathway. It is well-known that elevated homocysteine levels are associated with the early development of cardiovascular diseases , especially increasing the risk in patients with diabetes mellitus (DM) [5, 6]. A study performed in type 1 DM (T1DM) patients revealed an association between homocysteine levels and diabetic complications , and recent systematic review with a meta-analysis studies concluded that high levels of homocysteine may be associated with type 2 DM (T2DM) progression  and the development of vascular complications of diabetes . This is because high blood glucose status, caused by insulin-resistance conditions like T2DM, induces decrements of methionine transmethylation, homocysteine transsulfuration, and clearance in type 2 diabetic subjects , leading to increased concentrations of plasma homocysteine [5, 8]. On the other hand, some studies reported that the levels of homocysteine were either reduced or normal under diabetic conditions [11, 12]. These inconsistent results seem to be attributable to varying extents of hepatic and renal dysfunction [11, 13], the methionine content of the diet [14, 15], and the clinical status of diabetic disease [12, 16, 17]. Here, we focused on the fact that a high-fat diet affects the metabolism of both methionine and homocysteine in a diabetic rat model . Thus, 0we hypothesized that long-term administration of a high-methionine diet to T2DM rats with normal renal function may create metabolic changes, culminating in an elevated level of circulating homocysteine. To verify this hypothesis, mathematical modeling of methionine metabolism was required to predict the levels of homocysteine derived from given amounts of methionine.
Although some earlier studies had developed several mechanistic mathematical modeling approaches in an effort to understand the MMC, these models specially focused on and explained the systematic changes in cellular metabolism of methionine and homocysteine [18,19,20], making it difficult to apply the theoretical predictions of such models to real in vivo situations. Therefore, the objective of the current study was to develop a mathematical model of the MMC using in vivo data from T2DM rats. And we aimed to predict plasma methionine and homocysteine concentrations using the model, following the ingestion of a methionine-rich diet.
A single dose of 0.8 mmol/kg methionine (a mixture of L-methionine 0.6 mmol/kg and L-methionine-d4 0.2 mmol/kg) was given intravenously to both the Zucker diabetic fatty (ZDF) rats (ZDF/Gmi fa/fa)  and the controls (ZDF/Gmi fa/?). Blood samples were collected at 0, 10, 30, 60, 120, 210, 300, and 420 min after administration. Methionine and homocysteine concentrations were determined using LC-ESI/MS/MS; the details of the study design were described earlier .
Mechanistic modeling of the methionine metabolism cycle (MMC)
All parameters of the MMC model were transformed to in vivo scales, with consideration of physiological aspects, including liver volume and systemic circulation of methionine. The rate constants of distribution between compartments were simplified, by reference to prior data, suggesting that biological processes followed generally first-order kinetics. The MMC model featured a total of seven compartments associated with the five compartments of methionine, homocysteine, S-adenosylmethionine (SAM), S-adenosyl- L-homocysteine (SAH), cystathione, and the other two compartments addressing systemic distribution of methionine and homocysteine (Fig. 1). The movements of each component among compartments were described by the following differential equations (Eq):
K CM was the rate constant for the movement of methionine from compartment 1 to 2, K MC was for movement from compartment 2 to 1, and K MS from compartment 2 to 3 (correspondence of methionine adenosyl transferase [MAT]). K EL represented the rate of methionine elimination, K SH was the transformation rate constant from SAH to homocysteine (correspondence of SAH hydrolase [SAHH]), K HM was the transformation rate constant from homocysteine to methionine (corresponding to the activity of betaine-homocysteine methyltransferase [BHMT]), and K HC was the transformation rate constant from homocysteine to cysteine (corresponding to the activity of cystathionine beta-synthase [CβS]).
Inter-individual variability (activity differences among subjects) was explained by the exponential relationships between the transformation rate constant from homocysteine to methionine (K HM , Eq. 8) and the distribution rate constant of homocysteine from the central to the peripheral compartment (K HP , Eq. 9). In addition, ZDF status was added to the MMC model as dichotomous covariate (0 or 1), to quantify diabetic effects on rate constants in the ZDF rats. Equations corresponding to the above explanations were as follows:
Residual variability was explained by a combined error model using both additive and proportional error (equation was not shown).
The performance of the MMC model was assessed using graphical and numerical diagnostic tools. The numerical analysis showed that the reduction in the objective function value (OFV) was at least 10.83 (P value < 0.001 for one degree of freedom). Diagnostic plots, including standard goodness-of-fit plots and graphical assessment of the distribution of conditional weighted residuals (CWRES), were also employed to establish whether the MMC model was adequate for describing the current dataset. In addition, the final model was evaluated via visual predictive check (VPC), where we graphically compared the original datasets and the simulated results, with 95% confidence intervals (CIs), at the 5th, 50th, and 95th percentiles.
Use of the MMC model to simulate a methionine-rich diet
We assumed that methionine and homocysteine metabolisms were altered when T2DM patients ingested a methionine-rich diet for 60 days. To estimate the changes in two component levels after administration of a methionine-rich diet, 1000 virtual datasets were generated with the aid of Monte Carlo simulations. All simulated subjects were given 1 mmol/kg methionine divided into three times daily. For methionine and homocysteine concentrations, the simulated predicted accumulation ratio, total area under the curves (AUC) to total duration (AUCtotal), the AUC to 8 h (AUC0-8h), and the AUC to the last time point (AUClast) values, were compared between ZDF and control rats. The predicted accumulation ratio was the ratio of the AUC8-16h to the AUC0-8h (on day 1). The maximum concentration (C max ), time to maximum concentration (T max ), and AUC values were calculated via non-compartmental analysis (NCA).
The MMC model building and simulation were performed using NONMEM® (ver. 7.3.0) assisted by Perl-Speaks-NONMEM (PsN, ver. 4.3.0). Model parameters were derived using first-order conditional estimation method with interactions permitted among the parameters (FOCE-INTER). NCA was conducted using Phoenix (ver. 6.4.0). For all tests, we considered a P-value < 0.05 to reflect statistical significance.
A mathematical MMC model for ZDF rats
A new mathematical MMC model was successfully developed using plasma measurements of both methionine and homocysteine, and was characterized using seven compartments. From the prior study , we found that the plasma concentration kinetics of methionine were similar in the ZDF and control groups, while homocysteine exhibited a different kinetic pattern. Furthermore, since the hepatic concentrations of SAM and SAH did not differ between control and ZDF rats, the K CM , K MC , K MS , K SS , K EL , and V c parameters were predicted without consideration of disease status. Changes in homocysteine levels in the ZDF rat, the left shift of T max , and the increment in C max , are all explained by an increment in K SH and decrements in both K HM and K HC . The K SH value of ZDF rats was almost 16% greater than that of controls (ZDF versus control, 13.5 versus 11.6 h− 1). Although the homocysteine elimination slope was greater in ZDF rats in the early phase of elimination, such elimination was delayed after 210 h compared with that of controls, explained by decreases in both K HM and K HC at this time. The final model and kinetic parameters are shown in Table 1.
The goodness-of-fit plot for the final model showed that the population predicted and individual predicted concentrations were in good agreement with the observed value. The CWRES versus time plot was unbiased (data not shown). The VPC plot successfully explained the plasma concentrations of methionine and homocysteine; the 95% CIs of the 5th, 50th, and 95th percentiles of the observational data were well-covered in the MMC model (Fig. 2).
Simulation of a methionine-rich diet
In the virtual scenario, both the methionine and homocysteine levels steadily increased in ZDF rats after high-dose methionine administration (Fig. 3). The concentrations of both components steadily increased over this time period and, after 60 days of administration, the level of methionine in the ZDF group was twice that of controls, triggering significant changes in the homocysteine level. Although the accumulation ratios of methionine and homocysteine at a dose of 0.8 mmol/kg methionine did not differ significantly between ZDF rats and controls (methionine in ZDF rats and controls, 1.22 and 1.10, P = 0.059; homocysteine in ZDF rats and controls, 1.23 and 2.78, P = 0.162), the ratio of homocysteine at a dose of 1 mmol/kg in ZDF rats was significantly greater than that in controls (ZDF rats and controls, 1.21 and 1.39, P = 0.003). At 1350 h, the AUC8h of homocysteine was 307.84% greater in ZDF rats than controls (ZDF vs. control, 47.31 ± 6.06 vs. 145.64 ± 13.67; P < 0.001).
In the current study, we aimed to develop an MMC model using in vivo data under diabetic conditions, and to predict the changes in methionine and homocysteine levels after high-dose methionine administration. The final model successfully reproduced the kinetic profiles of methionine and homocysteine mathematically, and the elevated homocysteine levels were predicted for repeated high methionine ingestion. Thus, we confirmed the clinical importance of diet regulation in diabetic patients, for maintaining homocysteine levels within the normal range.
A prior study showed that T max of homocysteine decreased and C max increased in ZDF rats relative to normal rats, which may have resulted in homocysteine levels in the terminal phase being lower in ZDF rats than in the controls . It seemed that the metabolism of homocysteine was upregulated under diabetic conditions . More rapid transformation of methionine to homocysteine was related to rapid systemic distribution of homocysteine in patients of high-glucose status. The increased distribution of homocysteine during the initial phase was characterized by increased SAAH enzyme activity , as well as by higher protein levels in ZDF rats than in the controls . The K SH value of ZDF rats was higher than that of controls, which meant that upregulation of the insulin-signaling pathway induced the increased SAHH activity . In addition, our model showed that elevated levels of enzymes catalyzing homocysteine production and recycling to methionine, corresponding with previous studies conducted using diabetic in vivo rat models [13, 24]. Among all catalyzing enzymes, BHMT played a significant role in catalyzing remethylation of homocysteine to methionine. In the early phase, BHMT activity was elevated under diabetic conditions, as a result of increased levels of mRNA encoding the enzyme. Thus, homocysteine was rapidly reduced, accompanied by remarkable reductions in the hepatic concentrations of betaine . However, the catalyzing rate of homocysteine gradually decreased as betaine was depleted by the activation of BHMT in terminal phase . Consequently, an accumulation of homocysteine occurred. Previously, we described how K HM and K HC value were lower in ZDF rats than in controls. On the basis of the calculated accumulation factor, we predicted that blood homocysteine levels would rise about three-fold after administration of methionine on three-times per day over 60 days (a total of 180 occasions). Furthermore, the decline in K HC was associated with the decrease in the level of CβS, which is predominantly associated with the trans-sulfuration pathway . A previous study revealed that hydrogen sulfide (H2S) signaling was impaired under diabetic conditions, which can contribute to mediating a reduction in CβS activity . From these MMC modeling results, we confirmed that the metabolisms of homocysteine under diabetic conditions were different than those under normal conditions.
Subsequently, we explored the influence of a high-methionine diet on homocysteine levels in T2DM rats, on the basis of our model. The methionine dose was determined based on previous reports that have evaluated a high-dose methionine diet in rats. Yagisawa M et al. injected with methionine at a dose of 0.67 mmol/kg as a high-dose methionine diet in rats for estimating the influence of betaine on cardiovascular diseases. In the study of Cao Y et al., rats were treated with saline and intravenous high-dose methionine (0.8 mmol/kg) to generate hyperhomocysteinemia model. The methionine dose was equivalent to the human dose of 0.12 g/kg by scaling the body size, which was higher dose used in Haulrik N et al. study . The data obtained after repeated methionine administrations over 2 months showed that homocysteine concentrations were higher in those with diabetes, regardless of individual differences in the enzymes of the MMC. The greater accumulation ratio of homocysteine with high-dose methionine led to a greater increase in the AUC of homocysteine, in the ZDF rats compared to the controls. This may imply an increased risk of hyperhomocysteinemia in diabetic patients. These results are consistent with those of previous studies that used protein-rich or high-fat diets. High amount of methionine intake induced more efficient homocysteine metabolism, explained by the adaptive activation of enzymes involved in the methionine cycle [27, 28]. Another study found that protein-rich breakfasts and dinners triggered acute diphasic changes in plasma methionine and homocysteine concentrations . Thus, in diabetic patients, a high-methionine diet can induce significant variation in methionine metabolism, which is associated with an increased risk of vascular complications.
This mathematical approach is both useful and necessary when attempting to identify modes of regulation of complex metabolic systems and to predict changes in methionine and homocysteine profiles affected by various factors that cannot be investigated in the lab. Nevertheless, our study had certain limitations. First, the model did not evaluate clinical factors other than disease status. As the data were obtained in a previous well-controlled experiment using a rat model, clinical diversity was low and inter- and intra-subject variabilities could not be considered. In addition, the simulation results were analyzed by direct extrapolation from in vivo rat data to humans. However, it has been demonstrated that the several enzymes influencing plasma homocysteine profiles are similar in the rat and human .
In conclusion, we found that methionine-rich foods elevated both methionine and homocysteine levels in ZDF model, regardless of individual differences in enzymes of the MMC cycle. It is well known that high homocysteine levels are significantly associated with both atherosclerosis and vascular disease, by inhibiting vascular endothelial cell growth , and hyperhomocysteinemia is considered to be a biomarker of atherosclerosis in mammals [31, 32], especially in patients with T2DM than those without . Thus, it is important that T2DM patients avoid high-methionine diets to regulate their homocysteine levels within the normal range, and to prevent cardiovascular disease and mortality.
Area under the concentration curve
- Cmax :
Conditional weighted residuals
First-order conditional estimation method with interactions
Liquid chromatography-electrospray ionization/mass/mass
Methionine metabolism cycle
Non-linear mixed effect modeling
Objective function value
Type 2 diabetic mellitus
- Tmax :
Time to maximum concentration
Visual predictive check
Zucker Diabetic Fatty
Finkelstein JD. Methionine metabolism in mammals. J Nutr Biochem. 1990;1(5):228–37.
Refsum H, Ueland PM, Nygard O, Vollset SE. Homocysteine and cardiovascular disease. Annu Rev Med. 1998;49:31–62.
Kang SS, Wong PW, Malinow MR. Hyperhomocyst(e)inemia as a risk factor for occlusive vascular disease. Annu Rev Nutr. 1992;12:279–98.
Kernan WN, Ovbiagele B, Black HR, Bravata DM, Chimowitz MI, Ezekowitz MD, Fang MC, Fisher M, Furie KL, Heck DV, Johnston SC, Kasner SE, Kittner SJ, Mitchell PH, Rich MW, Richardson D, Schwamm LH, Wilson JA, C.O.C. American Heart Association stroke council, C.O.C.C. Stroke nursing, D. Council on peripheral vascular, guidelines for the prevention of stroke in patients with stroke and transient ischemic attack: a guideline for healthcare professionals from the American Heart Association/American Stroke Association, stroke. a journal of cerebral circulation. 2014;45(7):2160–236.
Hoogeveen EK, Kostense PJ, Jakobs C, Dekker JM, Nijpels G, Heine RJ, Bouter LM, Stehouwer CD. Hyperhomocysteinemia increases risk of death, especially in type 2 diabetes : 5-year follow-up of the Hoorn study. Circulation. 2000;101(13):1506–11.
Okada E, Oida K, Tada H, Asazuma K, Eguchi K, Tohda G, Kosaka S, Takahashi S, Miyamori I. Hyperhomocysteinemia is a risk factor for coronary arteriosclerosis in Japanese patients with type 2 diabetes. Diabetes Care. 1999;22(3):484–90.
Agullo-Ortuno MT, Albaladejo MD, Parra S, Rodriguez-Manotas M, Fenollar M, Ruiz-Espejo F, Tebar J, Martinez P. Plasmatic homocysteine concentration and its relationship with complications associated to diabetes mellitus, Clinica chimica acta. int j clin chem. 2002;326(1–2):105–12.
Huang T, Ren J, Huang J, Li D. Association of homocysteine with type 2 diabetes: a meta-analysis implementing Mendelian randomization approach. BMC Genomics. 2013;14:867.
Hoogeveen EK, Kostense PJ, Beks PJ, Mackaay AJ, Jakobs C, Bouter LM, Heine RJ, Stehouwer CD. Hyperhomocysteinemia is associated with an increased risk of cardiovascular disease, especially in non-insulin-dependent diabetes mellitus: a population-based study. Arterioscler Thromb Vasc Biol. 1998;18(1):133–8.
Tessari P, Coracina A, Kiwanuka E, Vedovato M, Vettore M, Valerio A, Zaramella M, Garibotto G. Effects of insulin on methionine and homocysteine kinetics in type 2 diabetes with nephropathy. Diabetes. 2005;54(10):2968–76.
Davies L, Wilmshurst EG, McElduff A, Gunton J, Clifton-Bligh P, Fulcher GR. The relationship among homocysteine, creatinine clearance, and albuminuria in patients with type 2 diabetes. Diabetes Care. 2001;24(10):1805–9.
Yun KU, Ryu CS, Lee JY, Noh JR, Lee CH, Lee HS, Kang JS, Park SK, Kim BH, Kim SK. Hepatic metabolism of sulfur amino acids in db/db mice, Food and chemical toxicology. int j pub British Industrial Bio R Assoc. 2013;53:180–6.
Stead LM, Brosnan ME, Brosnan JT. Characterization of homocysteine metabolism in the rat liver. Biochem j. 2000;350(Pt 3):685–92.
Kwak HC, Kim YM, Oh SJ, Kim SK. Sulfur amino acid metabolism in Zucker diabetic fatty rats. Biochem Pharmacol. 2015;96(3):256–66.
Yun KU, Ryu CS, Oh JM, Kim CH, Lee KS, Lee CH, Lee HS, Kim BH, Kim SK. Plasma homocysteine level and hepatic sulfur amino acid metabolism in mice fed a high-fat diet. Eur J Nutr. 2013;52(1):127–34.
Abbasi F, Facchini F, Humphreys MH, Reaven GM. Plasma homocysteine concentrations in healthy volunteers are not related to differences in insulin-mediated glucose disposal. Atherosclerosis. 1999;146(1):175–8.
Meigs JB, Jacques PF, Selhub J, Singer DE, Nathan DM, Rifai N, D'Agostino RB Sr, Wilson PW, Framingham S. Offspring, fasting plasma homocysteine levels in the insulin resistance syndrome: the Framingham offspring study. Diabetes Care. 2001;24(8):1403–10.
Martinov MV, Vitvitsky VM, Mosharov EV, Banerjee R, Ataullakhanov FI. A substrate switch: a new mode of regulation in the methionine metabolic pathway. J Theor Biol. 2000;204(4):521–32.
Lee EH, Lee K, Kwak GH, Park YS, Lee KJ, Hwang KY, Kim HY. Evidence for the dimerization-mediated catalysis of methionine sulfoxide reductase a from clostridium oremlandii. PLoS One. 2015;10(6):e0131523.
Reed MC, Nijhout HF, Sparks R, Ulrich CM. A mathematical model of the methionine cycle. J Theor Biol. 2004;226(1):33–43.
Wang H, Yoshizumi M, Lai K, Tsai JC, Perrella MA, Haber E, Lee ME. Inhibition of growth and p21ras methylation in vascular endothelial cells by homocysteine but not cysteine. J Biol Chem. 1997;272(40):25380–5.
Chiang EP, Wang YC, Chen WW, Tang FY. Effects of insulin and glucose on cellular metabolic fluxes in homocysteine transsulfuration, remethylation, S-adenosylmethionine synthesis, and global deoxyribonucleic acid methylation. J Clin Endocrinol Metab. 2009;94(3):1017–25.
Depuydt G, Xie F, Petyuk VA, Smolders A, Brewer HM, Camp DG 2nd, Smith RD, Braeckman BP. LC-MS proteomics analysis of the insulin/IGF-1-deficient Caenorhabditis Elegans daf-2(e1370) mutant reveals extensive restructuring of intermediary metabolism. J Proteome Res. 2014;13(4):1938–56.
Wijekoon EP, Hall B, Ratnam S, Brosnan ME, Zeisel SH, Brosnan JT. Homocysteine metabolism in ZDF (type 2) diabetic rats. Diabetes. 2005;54(11):3245–51.
Imbard A, Benoist JF, Esse R, Gupta S, Lebon S, De Vriese AS, Ogier De Baulny H, Kruger W, Schiff M, Blom HJ. High homocysteine induces betaine depletion. Biosci Rep. 2015;35(4):e00222.
Manna P, Gungor N, McVie R, Jain SK. Decreased cystathionine-gamma-lyase (CSE) activity in livers of type 1 diabetic rats and peripheral blood mononuclear cells (PBMC) of type 1 diabetic patients. J Biol Chem. 2014;289(17):11767–78.
Haulrik N, Toubro S, Dyerberg J, Stender S, Skov AR, Astrup A. Effect of protein and methionine intakes on plasma homocysteine concentrations: a 6-mo randomized controlled trial in overweight subjects. Am J Clin Nutr. 2002;76(6):1202–6.
Loehrer FM, Schwab R, Angst CP, Haefeli WE, Fowler B. Influence of oral S-adenosylmethionine on plasma 5-methyltetrahydrofolate, S-adenosylhomocysteine, homocysteine and methionine in healthy humans. J Pharmacol Exp Ther. 1997;282(2):845–50.
Guttormsen AB, Schneede J, Fiskerstrand T, Ueland PM, Refsum HM. Plasma concentrations of homocysteine and other aminothiol compounds are related to food intake in healthy human subjects. J Nutr. 1994;124(10):1934–41.
Garrido F, Gasset M, Sanz-Aparicio J, Alfonso C, Pajares MA. Rat liver betaine-homocysteine S-methyltransferase equilibrium unfolding: insights into intermediate structure through tryptophan substitutions. Biochem j. 2005;391(Pt 3):589–99.
Castro R, Rivera I, Struys EA, Jansen EE, Ravasco P, Camilo ME, Blom HJ, Jakobs C, Tavares de Almeida I. Increased homocysteine and S-adenosylhomocysteine concentrations and DNA hypomethylation in vascular disease. Clin Chem. 2003;49(8):1292–6.
Liu C, Wang Q, Guo H, Xia M, Yuan Q. Y. Hu, H. Zhu, M. Hou, J. Ma, Z. Tang, W. Ling, plasma S-adenosylhomocysteine is a better biomarker of atherosclerosis than homocysteine in apolipoprotein E-deficient mice fed high dietary methionine. J Nutr. 2008;138(2):311–5.
This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Science, ICT & Future Planning (Grant number 2009–0093815 and 2014R1A1A1006006) and supported by a grand of the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (Grant number: HI17C0927). This work was also supported by research funds from Chungnam National University.
Availability of data and materials
Please, contact corresponding author for data requests.
Ethics approval and consent to participate
All animal experiments were approved by the Institutional Animal Care and Use Committee in Chungnam National University (CNU), and were performed under institutional guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Nayoung Han and Jung-woo Chae contributed equally to this work as the first author.
Sang Kyum Kim and Hwi-yeol Yun contributed equally to this work as the corresponding author.