We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
Association between the DNA methylations of POMC, MC4R, and HNF4A and metabolic profiles in the blood of children aged 7–9 years
BMC Pediatricsvolume 18, Article number: 121 (2018)
Proopiomelanocortin (POMC), melanocortin 4 receptor (MC4R), and hepatocyte nuclear factor 4 alpha (HNF4A) are closely associated with weight gain and metabolic traits. In a previous study, we demonstrated associations between the methylations of POMC, MC4R, and HNF4A and metabolic profiles at birth. However, little is known about these associations in obese children. To evaluate the clinical utility of epigenetic biomarkers, we investigated to determine whether an association exists between the methylations of POMC, MC4R, and HNF4A and metabolic profiles in blood of normal weight and overweight and obese children.
We selected 79 normal weight children and 41 overweight and obese children aged 7–9 years in the Ewha Birth and Growth Cohort study. POMC methylation levels at exon 3, and MC4R and HNF4A methylation levels in promoter regions were measured by pyrosequencing. Serum glucose, total cholesterol (TC), triglyceride, high-density lipoprotein cholesterol (HDL–c), and insulin levels were analyzed using a biochemical analyzer and an immunoradiometric assay. Partial correlation and multiple regression analysis were used to assess relationships between POMC, MC4R, and HNF4A methylation levels and metabolic profiles.
Significant correlations were found between POMC methylation and HDL–c levels, and between HNF4A methylation and both TC and HDL–c levels. Interestingly, associations were found between POMC methylation status and HDL–c levels, and between HNF4A methylation status and TC levels independent of body mass index.
These findings show that POMC, MC4R, and HNF4A methylation status in the blood of children are associated with metabolic profiles. Therefore, we suggest that the DNA methylation status might serve as a potential epigenetic biomarkers of metabolic syndrome.
The prevalence of obesity and its related disorders have been increasing worldwide over the last few decades . In particular, the rates of childhood obesity has increased dramatically in many countries, including South Korea [2, 3]. It has been noted that around 40% of overweight children continue to be overweight during adolescence and that up to 80% of obese adolescents become obese adults . Obesity in childhood can therefore severely impact future health outcomes .
Obesity is also an important risk factor for the developments of type 2 diabetes, metabolic syndrome (MetS), and cardiovascular disease . Although, the molecular mechanisms underlying the pathogenesis of obesity remain largely unknown, one study concluded that heritability could account for 85% of obesity , whereas in another it was estimated heredity accounted for < 2% . The diversity of these findings suggests that obesity is the result of complex interactions between genetic and environment factors.
Epigenetic modifications, especially DNA methylation, are associated with changes of gene expression through binding of transcription factor . The establishment and maintenance of methylation are known to be sensitive to various environmental factors and to be heritable, although the mechanism underlying this inheritance in mammals remains unclear . In rats, maternal diet during pregnancy has been shown to alter DNA methylation patterns and to affect phenotypic consequences . Similarly, human studies have shown that prenatal exposure to famine led to changes in the methylation patterns of genes involved in lipid metabolism .
We previously reported that the methylation levels of proopiomelanocortin (POMC), melanocortin 4 receptor (MC4R), and hepatocyte nuclear factor 4 alpha (HNF4A) in cord blood were associated with metabolic profiles [13, 14]. POMC and MC4R genes play critical roles in the regulation of appetite, body weight, and energy homeostasis . The disruptions of POMC or MC4R in humans have been linked to early onset obesity [16, 17]. POMC is cleaved into α-melanocortin stimulating hormones that acts on MC4R in the hypothalamus . MC4R deficiency has been shown to be associated with elevated insulin levels and accelerated linear growth in childhood . HNF4A encodes a transcription factor relevant to glucose and lipid metabolism and has two distinct promoters, P1 and P2, which are separated by more than 45.6 kb [18, 19]. Furthermore, HNF4A mutation has been reported to be associated with impaired insulin secretion, resulting in maturity-onset diabetes of the young , and the expression of genes involved in liver and beta-cell functions, such as, glucose transport and glycolysis . To date, most studies have focused on the associations between genetic variants of POMC, MC4R, and HNF4A and metabolic phenotypes [20,21,22].
In the present study, we analyzed methylation levels of POMC, MC4R, and HNF4A in the blood of normal weight and overweight and obese children aged 7–9 years. We also evaluated associations between these methylation status and metabolic profiles.
This study was performed on 79 normal weight and 41 overweight and obese children aged 7–9 years that were selected in the Ewha Birth and Growth Cohort study. Methodological details of the cohort study have been reported elsewhere [13, 23]. We choose those with an age- and gender-matched body mass index (BMI) ≥ 85th percentile as overweight and obese children, and those with BMI < 85th percentile as normal weight children based on the 2007 Korean Children and Adolescent Growth Standards . Written informed consent was obtained from parents or guardians with respect to participation in the study. The study was approved by the Institutional Review Board of the Ewha Womans University Hospital (ECT 13-01A-13).
Anthropometric measurements and blood sampling were carried out according to the standard procedures by trained researchers. Height and weight were measured to the nearest 0.1 cm and 0.1 kg using an automatic electronic scale (Dong Sahn Jenix Co. Ltd., Seoul, Korea). BMI was calculated by dividing weight by height squared (kg/m2).
Fasting blood samples were collected in 10 ml Vacutainer tubes containing EDTA or serum tubes (BD Biosciences, San Jose, CA). Blood samples were centrifuged at 3000 rpm for 10 min, and the serum so obtained was stored at − 80 °C until required for chemical analysis. Serum glucose, triglyceride (TG), total cholesterol (TC), and high-density lipoprotein cholesterol (HDL-c) were analyzed by an automatic analyzer (model 7180; Hitachi, Tokyo, Japan). Insulin levels were measured using an immunoradiometric assay kit (MyBiosource, San Diego, CA). Insulin resistance was determined by the commonly used homeostasis model assessment of insulin resistance (HOMA-IR) method, which was calculated as (plasma glucose [mmol/L] × insulin [μIU/mL]) /22.5.
DNA methylation analysis by pyrosequencing
Genomic DNA samples were extracted from blood samples of children using the DNeasy Blood and Tissue kit (Qiagen, Hilden, Germany), according to manufacturer’s protocols. Concentrations and purities of DNA were measured using a Nanodrop spectrophotometer (Nanodrop Technologies, Wilmington, DE). Methylation target regions of the POMC, MC4R, and HNF4A were amplified by primer sets designed using the PSQ Assay Design software (Biotage AB, Uppsala, Sweden) [13, 14]. In POMC, the 4 CpG sites are located in the exon 3 region up to 7485 to 7520 bp downstream of the transcription start site (TSS) (Fig. 1a, 13). In MC4R, the 3 CpG sites are located in the promoter region up to − 800 to − 783 bp upstream of the TSS (Fig. 1b, 14). The HNF4A gene has 2 promoters (P1 and P2). The 4 CpG sites in P1 promoter region are located up to − 119 to − 92 bp and 4 CpG sites in P2 promoter region are located up to − 44,924 to − 44,907 bp upstream of the TSS (Fig. 1c, 14). Of these, the CpG2 site in P2 promoter of HNF4A was almost fully methylated (> 99%), and thus, was excluded from this study. Details of pyrosequencing analysis have been described in our previous study [13, 14]. The PCR conditions were as follows; 94 °C for 10 min, followed by 45 amplification cycles (denaturation 94 °C for 30 s, annealing for 30 s, and extension 72 °C for 30 s). For the annealing step, the temperature were set at 52, 53, 54, or 56 °C for POMC, MC4R, HNF4A-P1, and HNF4A-P2, respectively. After amplification, all reaction were incubated at 72 °C for 10 min and then cooled to 4 °C.
Quantitative data were expressed as means ± standard deviations and were analyzed using the Student t-test or analysis of covariance, adjusted for age. TG and insulin levels were analyzed as log-transformed values to satisfy normality. Partial correlation was used to identify the relationship between POMC, MC4R, and HNF4A methylation levels and metabolic profiles after adjusting for age in all subjects. To control for the effects of BMI as a major risk factor of MetS, multiple linear regression analysis was used to assess the associations between DNA methylations of the three genes and metabolic profiles in all subjects. In addition, age and gender were adjusted in the analysis. The statistical analysis was conducted using SAS software (ver. 9.4; SAS Institute Inc., Cary, NC). All analyses were two-tailed, and p values of < 0.05 were considered as significant.
General characteristics of the study subjects
Seventy-nine normal weight and 41overweight and obese children were included in this study (Table 1). There were significant differences between the two groups in age, height, weight, BMI, and waist circumference (p < 0.001). After adjusting for age, average TG, TC, insulin, glucose levels and HOMA index were significantly higher in overweight and obese children (p < 0.05), whereas mean HDL-c level was significantly lower (p < 0.001).
Average DNA methylation levels
The following methylation levels were significantly less methylated in overweight and obese children (p < 0.05): POMC methylation at 4 CpG sites, MC4R methylation at CpG3, and HNF4A methylation levels at CpG2, 3, and 4 in P1. On the contrary, the following methylation levels were significantly more methylated in overweight and obese children (p < 0.05): MC4R methylation at CpG1 and 2, and HNF4A methylation at CpG1 in P1 and at CpG1 and 3 in P2, after adjusting for age (Table 2).
The relationship between DNA methylation and blood metabolic profiles
Table 3 shows significant correlations between hypermethylation of POMC at 4 CpG sites, and between hypermethylation of HNF4A at CpG4 in P1 and lower HDL-c levels (r = − 0.19 ~ 0.23, p < 0.05). In addition, hypomethylation of HNF4A at CpG4 in P1, and hypermethylation of HNF4A at CpG3 in P2 were significantly correlated with higher TC levels (r = − 0.32, p = 0.001, r = 0.30, p = 0.001, respectively). Hence, we focused on these 6 CpG sites and 2 metabolic profiles showing significant relation (p < 0.05).
The association between DNA methylation and blood metabolic profiles in all subjects
To control for the effect of BMI on metabolic profiles, we analyzed associations between the methylation status of POMC, MC4R, and HNF4A and metabolic profiles after adjusting for age, gender, and BMI in all subjects (Table 4). We observed a hypermethylation of POMC-CpG2 was significantly related to a lower HDL-c level (β = − 0.23, p = 0.048). Hypomethylation of HNF4A-CpG4 in P1 and hypermethylation of HNF4A-CpG3 in P2 were significantly associated with higher TC levels (β = − 2.79, p = 0.001; β = 5.70, p = 0.002).
Epigenetic variations in metabolism–related genes partly contributed to the development of obesity and MetS. Most studies have reported the associations between POMC, MC4R, and HNF4A variants and obesity and MetS-related phenotypes [16,17,18,19,20]. In this study, we identified significant correlations between the methylation status of POMC and HNF4A and metabolic profiles in the blood of normal weight and overweight and obese children. Moreover, BMI-independent associations were found between the methylation status of POMC and HNF4A, and TC and HDL–c levels in Korean children. These findings suggest that epigenetic changes in metabolism-related genes may influence metabolic profiles and could serve as potential epigenetic biomarkers for MetS.
In the present study, mean TG, TC, insulin, and glucose levels and HOMA index were significantly higher in overweight and obese children, whereas mean HDL–c levels had significantly lower than in normal weight children. Ba et al. also found that mean TG, TC, and insulin levels and HOMA index were significantly higher in obese children than in normal children . Similarly, Azab et al. reported that obese children had significantly higher TG and TC levels, and lower HDL-c levels than non-obese children . However, Woo et al. reported that mean TC levels were not different between obese and non-obese children, although mean TG and insulin levels were significantly different .
Furthermore, our results indicate average POMC, MC4R, and HNF4A methylation levels at CpG sites were significantly different in overweight and obese and normal weight children. Interestingly, our earlier study found that POMC, MC4R, and HNF4A methylation in specific CpG sites were significantly different in term and preterm infants .
In this study, the methylation levels at some CpG sites differed by < 1% in the two groups. Although these differences were quite marginal, they were found to be statistically significant. These results may be due to higher rate of overweight and obese children (34.2%) in this study rather than in the general population. A lower rate of overweight and obese children may have reduced the bias and possibly increased the difference in methylation levels between the two groups. After adjusting for age, we found a negative correlation between POMC methylation levels and HDL–c levels in children. In a previous study, we identified that high POMC methylation levels were marginally related to lower HDL–c levels . Although there is no significant association between MC4R methylation and insulin levels in this study, larger sample size and further study are needed to clarify this association. Martinelli et al. reported that MC4R–deficient children had considerably higher insulin level . Furthermore, a study on MC4R–null mice reported elevated insulin levels . In this study, negative correlations were found between HNF4A methylation levels in P1 and TC and HDL–c levels, whereas HNF4A methylation levels in P2 were positively correlated with TC levels. Hayhurst et al. found that HNF4A–null mice showed markedly lower TC and TG levels . In addition, we previously reported significant associations between the hypomethylation status of MC4R and HNF4A and higher TG levels in term infants and preterm infants . These findings suggest that the methylation status of POMC, MC4R, and HNF4A may contribute to the metabolic phenotypes of children.
Interestingly, we observed a significant association between hypermethylation of POMC–CpG2 in exon 3 and lower HDL–c levels after adjusting for age, gender, and BMI. Crujeiras et al. showed hypermethylation of the promoter region of POMC was associated with weight gain after an energy-restricted diet . In addition, hypermethylation of exon 3 of POMC has been shown to be highly related to childhood BMI, and to be a relatively stable region in different cell types . Biological actions of POMC–derived peptides are mediated by melanocortin receptors . Perez-Tilve et al. reported that inhibition of melanocortin signaling led to higher HDL-c levels independent of body weight . In the current study, hypermethylation of the promoter of MC4R-CpG1 tended to be related to higher insulin levels in children, although not statistically significant. In contrast, Widiker et al. found the hypomethylation status of MC4R was related with obesity in mice fed a high–fat diet . Likewise, the hypomethylation of MC4R was found to be significantly linked with elevated TG levels in cord blood of term infants and preterm infants . These observations suggest that at some CpG sites, MC4R methylation patterns in early life may be influenced by environmental factors rather than by heritable traits. Longitudinal studies and larger sample sizes are required to determine whether MC4R methylation status could impair metabolic profiles and possibly increase the risk of MetS. This study demonstrated BMI-independent associations between the methylation status of HNF4A-CpG4 in P1 and of HNF4A-CpG3 in P2, and TC levels. In humans, HNF4A methylation levels have been reported to differ in neonates with intrauterine growth restriction . These findings concur with those of Ribel–Madsen et al. who found HNF4A was hypermethylated in individuals with type 2 diabetes .
This study had several limitations. First, we did not examine methylation levels in target organ tissues, such as, hypothalamus or liver. In addition, we did not take into account cellular heterogeneity in blood samples which may influence methylation levels. Because the study was conducted on children, it was difficult to obtain tissues and samples sizes were limited, and thus, we evaluated DNA methylation levels in whole blood, which is a readily accessible for the assessment of biomarkers. Second, we did not evaluate the environment or social-economic factors, which may affect methylation. Nonetheless, it was identified the methylation status of POMC, MC4R and HNF4A were related to metabolic profiles in Korean children. We also confirmed the associations between the methylation status of three genes and metabolic profiles in 1) cord blood of term and preterm infants in a case-control study ; 2) cord blood and blood in 90 children aged 7–9 years in a cohort study . Additionally, these methylation status could influence metabolic profiles. Consequently, we suggest that changes in DNA methylation patterns in childhood may increase the risk of MetS, and longitudinal studies with larger sample sizes would be needed to investigate further.
This study demonstrates the associations between the methylation status of POMC, MC4R, and HNF4A and metabolic profiles in Korean children. Furthermore, we found BMI–independent associations between DNA methylations of three genes and TC and HDL–c levels. These findings indicate the clinical utility as potential epigenetic biomarkers of MetS, and larger cohort studies are required for further evaluation and verification.
Body mass index
High-density lipoprotein cholesterol
Hepatocyte nuclear factor 4 alpha
Melanocortin 4 receptor
Transcription start site
Ebbeling CB, Db P, Ludwig DS. Childhood obesity: public-health crisis, common sense cure. Lancet. 2002;360:473–782.
Jang M, Berry D. Overweight, obesity, and metabolic syndrome in adults and children in South Korea: a review of the literature. Clin Nurs Res. 2011;20:276–91.
Moon JS. Secular trends of body sizes in Korean children and adolescents: from 1965 to 2010. Korean J Pediatr. 2011;54:436–42.
Lifshitz F. Obesity in children. J Clin Res Pediatr Endocrinol. 2008;1:53–60.
Freedman DS, Khan LK, Serdula MK, Dietz WH, Srinivasan SR, Berenson GS. The relation of childhood BMI to adult adiposity: the Bogalusa heart study. Pediatrics. 2005;115:22–7.
Ford ES, Giles WH, Dietz WH. Prevalence of the metabolic syndrome among US adults: findings from the third National Health and nutrition examination survey. JAMA. 2002;287:356–9.
Yang W, Kelly T, He J. Genetic epidemiology of obesity. Epidemiol Rev. 2007;29:49–61.
Waalen J. The genetics of human obesity. Transl Res. 2014;164:293–301.
Henikoff S, Matzke MA. Exploring and explaining epigenetic effects. Trends Genet. 1997;13:293–5.
Heard E, Martienssen RA. Transgenerational epigenetic inheritance: myths and mechanisms. Cell. 2014;157:95–109.
Lillycrop KA, Slater-Jefferies JL, Hanson MA, Godfrey KM, Jackson AA, Burdge GC. Induction of altered epigenetic regulation of the hepatic glucocorticoid receptor in the offspring of rats fed a protein-restricted diet during pregnancy suggests that reduced DNA methyltransferase-1 expression is involved in impaired DNA methylation and changes in histone modifications. Br J Nutr. 2007;97:1064–73.
Tobi EW, Goeman JJ, Monajemi R, Gu H, Putter H, Zhang Y, Slieker RC, Stok AP, Thijssen PE, Müller F, van Zwet EW, Bock C, Meissner A, Lumey LH, Eline Slagboom P, Heijmans BT. DNA methylation signatures link prenatal famine exposure to growth and metabolism. Nat Commun. 2014;5:5592.
Yoo JY, Lee S, Lee HA, Park H, Park YJ, Ha EH, Kim YJ. Can proopiomelanocortin methylation be used as an early predictor of metabolic syndrome? Diabetes Care. 2014;37:734–9.
Kwon EJ, Lee HA, You YA, Park H, Cho SJ, Ha EH, Kim YJ. DNA methylations of MC4R and HNF4α are associated with increased triglyceride levels in cord blood of preterm infants. Medicine (Baltimore). 2016;95:e4590.
Berthoud HR, Morrison C. The brain, appetite, and obesity. Annu Rev Psychol. 2008;59:55–92.
Farooqi IS, Drop S, Clements A, Keogh JM, Biernacka J, Lowenbein S, Challis BG, O'Rahilly S. Heterozygosity for a POMC-null mutation and increased obesity risk in humans. Diabetes. 2006;55:2549–53.
Martinelli CE, Keogh JM, Greenfield JR, Henning E, van der Klaauw AA, Blackwood A, O'Rahilly S, Roelfsema F, Camacho-Hübner C, Pijl H, Farooqi IS. Obesity due to melanocortin 4 receptor (MC4R) deficiency is associated with increased linear growth and final height, fasting hyperinsulinemia, and incompletely suppressed growth hormone secretion. J Clin Endocrinol Metab. 2011;96:E181–8.
Stoffel M, Duncan SA. The maturity-onset diabetes of the young (MODY1) transcription factor HNF4alpha regulates expression of genes required for glucose transport and metabolism. Proc Natl Acad Sci U S A. 1997;94:13209–14.
Love-Gregory LD, Wasson J, Ma J, Jin CH, Glaser B, Suarez BK, Permutt MA. A common polymorphism in the upstream promoter region of the hepatocyte nuclear factor-4 alpha gene on chromosome 20q is associated with type 2 diabetes and appears to contribute to the evidence for linkage in an ashkenazi jewish population. Diabetes. 2004;53:1134–40.
Lehto M, Bitzén PO, Isomaa B, Wipemo C, Wessman Y, Forsblom C, Tuomi T, Taskinen MR, Groop L. Mutation in the HNF-4alpha gene affects insulin secretion and triglyceride metabolism. Diabetes. 1999;48:423–5.
Pritchard LE, Turnbull AV, White A. Pro-opiomelanocortin processing in the hypothalamus: impact on melanocortin signaling and obesity. J Endocrinol. 2002;172:411–21.
Farooqi IS, Keogh JM, Yeo GS, Lank EJ, Cheetham T, O'Rahilly S. Clinical spectrum of obesity and mutations in the melanocortin 4 receptor gene. N Engl J Med. 2003;348:1085–95.
Min JW, Kong KA, Park BH, Hong JH, Park EA, Cho SJ, Ha EH, Park H. Effect of postnatal catch-up growth on blood pressure in children at 3 years of age. J Hum Hypertens. 2007;21:868–74.
Korea Pediatric Society: Korean infants and adolescents growth standards. Available from: http://www.cdc.go.kr/CDC/notice/CdcKrInfo0301.jsp?menuIds=HOME001-MNU0004-MNU0036-MNU0037&cid=12103. Accessed Novermber 2007.
Ba HJ, Chen HS, Su Z, Du ML, Chen QL, Li YH, Ma HM. Associations between serum apelin-12 levels and obesity-related markers in Chinese children. PLoS One. 2014;9:e86577.
Azab SF, Saleh SH, Elsaeed WF, Elshafie MA, Sherief LM, Esh AM. Serum trace elements in obese Egyptian children: a case-control study. Ital J Pediatr. 2014;40:20.
Woo KS, Chook P, Yu CW, Sung RY, Qiao M, Leung SS, Lam CW, Metreweli C, Celermajer DS. Overweight in children is associated with arterial endothelial dysfunction and intima-media thickening. Int J Obes Relat Metab Disord. 2004;28:852–7.
Fan W, Dinulescu DM, Butler AA, Zhou J, Marks DL, Cone RD. The central melanocortin system can directly regulate serum insulin levels. Endocrinology. 2000;141(9):3072.
Hayhurst GP, Lee YH, Lambert G, Ward JM, Gonzalez FJ. Hepatocyte nuclear factor alpha (nuclear receptor 2A1) is essential for maintenance of hepatic gene expression and lipid homeostasis. Mol Cell Biol. 2001;21:1393–403.
Crujeiras AB, Campion J, Díaz-Lagares A, Milagro FI, Goyenechea E, Abete I, Casanueva FF, Martínez JA. Association of weight regain with specific methylation levels in the NPY and POMC promoters in leukocytes of obese men: a translational study. Regul Pept. 2013;186:1–6.
Kuehnen P, Mischke M, Wiegand S, Sers C, Horsthemke B, Lau S, Keil T, Lee YA, Grueters A, Krude H. An Alu element-associated hypermethylation variant of the POMC gene is associated with childhood obesity. PLoS Genet. 2012;8:e1002543.
Perez-Martinez P, Garcia-Rios A, Delgado-Lista J, Delgado-Casado N, Malagon MM, Marin C, Gomez-Luna P, Caballero J, Perez-Jimenez F, Lopez-Miranda J. A variant near the melanocortin-4 receptor gene regulates postprandial lipid metabolism in a healthy Caucasian population. Br J Nutr. 2011;106:468–71.
Widiker S, Karst S, Wagener A, Brockmann GA. High-fat diet leads to a decreased methylation of the Mc4r gene in the obese BFMI and the lean B6 mouse lines. J Appl Genet. 2010;51:193–7.
Einstein F, Thompson RF, Bhagat TD, Fazzari MJ, Verma A, Barzilai N, Greally JM. Cytosine methylation dysregulation in neonates following intrauterine growth restriction. PLoS One. 2010;5:e8887.
Ribel-Madsen R, Fraga MF, Jacobsen S, Bork-Jensen J, Lara E, Calvanese V, Fernandez AF, Friedrichsen M, Vind BF, Højlund K, Beck-Nielsen H, Esteller M, Vaag A, Poulsen P. Genome-wide analysis of DNA methylation differences in muscle and fat from monozygotic twins discordant for type 2 diabetes. PLoS One. 2012;7:e51302.
The authors are thankful to all of the consultants from the Cohort group.
This study was supported by the Basic Science Research Program through the National Research Foundation of Korea funded by the Ministry of Education (2016 R1D1A1A09918620) and by the Ministry of Health & Welfare of the Republic of Korea (HI15C2059) through the Korea Health Industry Development Institute.
Availability of data and materials
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Ethics approval and consent to participate
We obtained written informed consent from parents or guardians before enrollment of the participant in this study. This study was approved by the Institutional Review Board of the Ewha Womans University Hospital (ECT 13-01A-13).
Consent for publication
All authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.