Serum levels of the novel adipokine isthmin-1 are associated with obesity in pubertal boys

Implications and Contribution

Study findings reveal differences in the circulating levels of the novel insulin-like adipokine ISM1 between pubertal children with obesity and normal-weight, and different associations with metabolic traits, which may be further investigated to be used as a potential biomarker of obesity in young populations.

Childhood obesity is increasing globally and total adiposity is the key driver of metabolic risk in children and adolescent, which represent a strong risk factor for insulin resistance (IR) and future type 2 diabetes (T2D) [1]. Puberty is a time of metabolic and hormonal changes, and it is associated with a reduced insulin sensitivity that recovers at puberty completion in only some children [2-4]. In growing children, adipocyte hypertrophy is associated with inflammation and local and systemic IR, independently of BMI and fat mass, being adipose tissue essential to maintain a functional metabolism [5]. However, the molecular mechanism of IR is still unknown, particularly in children. This is in part because processes such as growth and puberty affect insulin secretion and sensitivity [6]. Hence, understanding the molecular and biological processes underlying metabolic changes (glucose and lipid regulation) during puberty, and identifying pathways and biomarkers that might help to increase peripheral glucose uptake would be beneficial to reduce the impact of obesity and to prevent T2D. A recent paper reported that the protein isthmin-1 (ISM1) is secreted by mature adipocytes and triggers a signaling cascade similar to that of insulin. The novel adipokine acts through a unidentified receptor tyrosine kinase and, at pharmacological doses in mice, ISM1 ameliorates metabolic disturbances associated with T2D, including hyperglycemia and liver steatosis [7, 8]. Here, we report that serum levels of the novel insulin-like adipokine ISM1 is, indeed, associated with obesity in pubertal boys, but not in girls in a well-characterized population of Spanish children under a cross-sectional design. Furthermore, we identified DNA methylation in two-enhancer-related CpG sites of the ISM1 region (cg14269097 and cg14269097) associated with serum levels of ISM1 in children with obesity.

Methods

This study was conducted within the context of the multicenter PUBMEP study in Spain “Puberty and metabolic risk in children with obesity”, previously published [9, 10]. Here, a sub-population of 119 pubertal children (54 boys and 65 girls) from the whole PUBMEP cohort was selected for analysis. A total of 38 were normal-weight (17 boys), 36 overweight (17 boys) and 45 children with obesity (20 boys). The following characteristics were considered as exclusion criteria: birth weight <2500 g; intake of any drug that could alter blood glucose, blood pressure or lipid metabolism; not being able to comply with the study procedures and being participating or having participated in the last three months in an investigation project. This study was conducted according to the guidelines set out in the Declaration of Helsinki (Edinburgh 2000 revised), and all procedures were approved by the Ethics and Research Committee of Galicia Autonomous Community (2011/198 and 2016/522). Written consent was obtained from the parents of all the children. Anthropometric measurements such as body weight (kg), height (cm), hip circumference (cm) and waist circumference (WC) (cm) were measured using standardized procedures. BMI z-score was calculated based on the Spanish reference standards published. Blood pressure was measured three times for each individual by the same examiner using a mercury sphygmomanometer and following international recommendations [11]. Measures of lipid and glucose metabolism, hormones and classical biochemical parameters were performed at the laboratories of each participating hospital following internationally accepted quality control protocols. Blood samples were collected in overnight fasting conditions, centrifuged, and plasma and serum were stored at -80°C. The presence of IR in children was defined according to the HOMA insulin resistance (HOMA-IR) index. The cut-off points were obtained from a previously well-described Spanish cohort composed of children and adolescents [12, 13]. The cut-off points for IR were based on the 95^th HOMA-IR percentile, considering sex (HOMA-IR ≥ 3.38 in boys and HOMA-IR ≥ 3.90 in girls). These cut-off points have already been tested and validated as good metabolic risk classifiers in our population according to the results from a previous PUBMEP report [10]. Plasma adipokines, inflammation, and cardiovascular risk biomarkers (adiponectin, leptin, resistin, tumor necrosis factor alpha (TNF-α), high-sensitivity CRP (hsCRP), interleukin (IL)-6, IL-8, total plasminogen activator inhibitor-1 (PAI-1), P-selectin, myeloperoxidase (MPO), monocyte chemoattractant protein 1 (MCP-1), soluble intercellular cell adhesion molecule-1 sICAM-1, and soluble vascular cell adhesion molecule-1 (sVCAM)) were analyzed in all samples and time points using XMap technology (Luminex Corporation, Austin, TX) and human monoclonal antibodies (Milliplex Map Kit; Millipore, Billerica, MA) as previously reported [10, 14]. Descriptive data are expressed as mean (standard deviation) or median [min-max] if not normally distributed. One-way ANOVA, Kruskal-Wallis and the Welch test were employed to assess group differences. DNA methylation analysis was carried out by using the Infinium MethylationEPIC microarray using bead chip technology (Illumina, San Diego, CA, USA) as previously described [10]. ISM1 protein levels were determined in serum using the Human SEQ515Hu for ISM 1 (Cloud-Clone Corp., USA), an enzyme-linked immune-absorbent assay kit according to the manufacturers’ instructions. The coefficient of variance was 4%. Two-way ANOVA and Tukey’s multiple comparisons test were employed to assess groups’ differences in ISM1 levels between boys and girls and normal-weight, overweight and obese children. Multiple linear regression (MLR) analyses were applied for all continuous variables to study their association with ISM1 levels. In these analyses, origin, age, gender, pubertal stage, height, BMI Z-score, and insulin were adjusted as confounders when necessary. A P-value <0.05 was considered as significant. Given the number of analyzed outcomes, we considered false discovery rate (FDR) as in Benjamini and Hochberg to correct for multiple hypothesis testing. MLRs were also applied for all calculated deltas to study their correlation with the change in ISM1 levels. All described analyses were performed in R environment version 3.6.0.

Results

General characteristics of the 119 children in this cross-sectional study are shown in Supplementary Table 1. ISM1 serum levels according to obesity by sex are shown in Fig. 1A. Higher ISM1 serum levels were observed in boys with obesity when compared with normal-weight (P = 0.0037) and overweight (P=0.0071), non-adjusted. However, no changes were observed in girls. When all subjects of the sample were compared together, we found higher levels in children with obesity compared with those with normal weight (P = 0.0413) and with those overweight (P = 0.0102). A trend to increase the ISM1 levels in boys with IR was observed, but no statistical differences were shown between normal-weight non-IR, non-IR and IR children (Fig. 1B). To elucidate the relationship between ISM1 and obesity, MLRs were further conducted in a wide range of metabolic outcomes separated by sex (Tables 1-2) and properly adjusted by confounders such as age, origin, height, BMI and IR when applicable. In boys, the strongest association was found for BMI Z-score (P=0.005), but not in girls (Fig. 2A and S1A). MPO (P=0.043) was negatively associated with ISM1 serum levels in boys. All of them were properly adjusted for confounders (please, see Table 1 and 2 legends), both in boys and girls. Surprisingly, though both ISM1 and leptin levels are significantly higher in children with obesity, no significant associations were identified with leptin neither in boys nor in girls (Fig. 2B). Fifty-one methylation sites were selected from the Infinium Methylation EPIC microarray of which two were annotated as promoter-associated CpGs (Supplementary Table 2). All the CpGs were annotated as open sea. We found a positive significant association between the methylation status of the probe cg03304641 and ISM1 serum levels in pubertal children (P=0.0062) (Fig. 3A, S2A and Supplementary Table 2). However, for the probe cg14269097 there was a negative association with ISM1 serum levels (P=0.038) (Fig. 3B, S2B and Supplementary Table 2).

• Download figure
• Open in new tab

Figure 1.

Group comparisons for ISM1 serum levels (ng/mL) in the pubertal population of 119 children. A) Comparison between normal-weight, overweight and obese. B) Comparison between normal-weight non-IR, non-IR and IR in children. The two-way ANOVA, Tukey’s multiple comparisons test was employed to assess group differences in ISM1 levels according to standard statistical assumptions. *P<0.05

• Download figure
• Open in new tab

Figure 2.

Multiple linear regressions analyses between the change in ISM1 serum levels (ng/mL) and the changes in BMI Z-score and leptin (µg/L) in the cross-sectional cohort. A reports the linear model with delta BMI Z-score as dependent variable, B refers to the model for delta leptin as dependent variable.

• Download figure
• Open in new tab

Figure 3.

Multiple linear regressions analyses between the change in ISM1DNAmethylation status and ISM1 serum levels (ng/mL). A reports the linear model with delta ISM1 serum levels as dependent variable and cg03304641 as independent variable. B refers the linear model with delta ISM1 serum levels as dependent variable and cg14269097 as independent variable.

Discussion

In the present work, we show higher serum levels of the novel insulin-like adipokine ISM1 in pubertal children with obesity, with a strong association with BMI Z-score in boys. Moreover, ISM1 was also associated with MPO, an enzyme associated with IR and inflammation in individuals with overweight [15]. These findings illustrate how this protein might be of special relevance as a new biomarker of obesity in children.

In obesity, hypertrophic adipocytes and adipose tissue-resident immune cells accelerate a low-grade and chronic, proinflammatory profile with altered secretion of adipokines and lipokines, thereby exacerbating cardiometabolic disease [16]. In this context, the production and secretion of adipokines and lipokines, which contribute to systemic energy metabolism by different mechanisms, are dependent on the energy status of adipose tissue. Hence, further preclinical and clinical studies exhibiting the activation or inhibiting the signaling of specific adipokines or lipokines (e.g. by using adipokine-neutralizing antibodies) may be an approach suitable to treat or prevent the development of metabolic disease. Nevertheless, efficacy and safety in humans need to be confirmed [17]. ISM1 is a secreted protein, originally discovered in fetal brain development and expressed in the brain, lung, vasculature, skin and immune cells [18]. ISM1 has been recently identified in mouse and human adipocytes, regulating glucose uptake while suppressing hepatic lipid synthesis, thus improving hyperglycemia and reducing lipid accumulation in mouse models. Besides, circulating plasma levels of human ISM1 has been detected at an average of 50 pg/mL and trend to positively correlate with BMI, but not with glucose in female individuals [7]. Now, we reveal a relationship between circulating levels of ISM1 and obesity is investigated in a pubertal population. The motivation for focusing on the identification of new biomarkers in puberty lies in the fact that sexual maturation has been presented as a significant metabolic risk period for children with obesity [4], and indeed, we previously described the role of S100A4 in IR through a multi-omics approach in children, providing an interesting knowledge into the plausible molecular mechanism underlying that association [10].

Here, in a pubertal sample of 119 children, we show how higher levels of ISM1 in children with obesity compared to normal-weight and overweight, and a robust positive association between ISM1 serum levels and BMI Z-score in boys, indicating a putative sexual dimorphism. However, MPO was negatively associated with circulating levels of ISM1. Indeed, MPO is the most abundant protein in human neutrophils, playing a major role in inflammation, oxidative stress, lipoprotein oxidation, and atherosclerosis [15, 19]. Moreover, MPO deficient mice are resistant to diet-induced obesity and IR, and inhibition of MPO activity in the neutrophils decreases diet-induced IR in obese mice, and activation of MPO is associated with the development of obesity and obesity-associated IR [20]. Group comparisons for ISM1 levels also revealed significant results in boys, although no significant changes were observed when comparing extreme experimental conditions in relation to IR (normal-weight vs. obese with IR). Jiang et al. (2021) [7] determined that ISM1 signaling is dependent on PI3K and shares downstream phosphorylation targets with insulin signaling, such as -AKT, p-AKT, p-ERK1/2 and p-S6. Outstandingly, ISM1 activates a PI3K-AKT signaling pathway independently of the insulin and insulin-like growth factor 1 receptors, being most likely to signal through another, yet to be identified, receptor tyrosine kinase. We found a correlation with obesity but not with HOMA-IR in children, which would point to the adipokine ISM1 with a direct role in obesity but not in the metabolic status derived from IR. Remarkably, while the glucoregulatory function of the novel adipokine is shared with insulin, ISM1 also neutralizes lipid accumulation in the liver by inhibiting de novo lipogenesis, promoting protein synthesis and preventing hepatic steatosis in a diet-induced fatty liver mouse model. Nevertheless, in addition to a disturbed hepatic and postprandial lipoprotein metabolism, enhanced triacylglycerol lipolysis in adipocytes and subsequent fatty acid flux to the liver are major determinants of hepatic steatosis [21]. As insulin is a major anti-lipolytic hormone in adipocytes, therefore, it is plausible that ISM1 signaling indirectly modulates hepatic lipid accumulation by inhibiting fatty acid release from adipose tissue in mice (Heeren 2021a). Additionally, in the paper of Jiang et al. (2021) [7], the therapeutic dosing of recombinant ISM1 improves glucose tolerance to the same degree as metformin, enhances diabetes in diet-induced obese mice and ameliorates hepatic steatosis in a diet-induced fatty liver mouse model, establishing that the recombinant Ism1 and its derivatives may be explored for therapeutic purposes and may offer certain advantages over current monotherapies. Nonetheless, the observed higher circulating ISM1 levels in pubertal children suggest that maybe an ISM1 resistance is present, as administration of ISM1 into mice with established disease improves glucose and lipid dysfunction in the diet-induced obesity. Additionally, we should consider that pubertal children did not show hyperglycemia, which it might mask the observed effects of ISM1 in mice.

The strengths of our findings are the relatively high number of recruited children from different centers in the country (Andalucia, Galicia and Aragón); the novelty of the recently described new adipokine, which has not been reported elsewhere in humans, and the possibility to correlate with several plasma adipokines, inflammation, and cardiovascular risk biomarkers in the subcohort. Furthermore, the availability of DNA methylation analysis by using the Infinium MethylationEPIC microarray using bead chip technology in all the children population allowed us to observe two-enhancer related CpG sites of ISM1 (cg03304641 and cg14269097) associated with serum levels of ISM1 in children. As a limitation, though ISM1 serum levels are associated with BMI Z-score, we cannot distinguish the origin of the circulating levels simply because ISM1 is secreted by other non-adipose cells, and a secondary validation population would be needed to confirm these findings.

In conclusion, the circulating levels of the novel insulin-like adipokine ISM1 are significantly higher in pubertal children with obesity, and strongly associated with BMI Z-score in boys. Furthermore, we reveal two DNA methylations in two-enhancer-related CpG sites of the ISM1 region associated with serum levels of the protein in children with obesity.

Conflicts of interest

The authors have no conflict of interests to disclose. No author received any form of payment to produce this paper.

Funding source

This work was supported by the Plan Nacional de Investigacio’n Científica, Desarrollo e Innovación Tecnológica (I + D + I), Instituto de Salud Carlos III-Health Research Funding (FONDOS FEDER) (PI051968, PI1102042 and PI1600871); Redes temáticas de Investigación cooperativa RETIC (Red SAMID RD12/0026/0015) and the Mapfre Foundation.

Author contributions

CAG is the guarantor of this work. AG and CAG contributed to the study concept of the cross-sectional study. FJRO, AAR and CAG design this study. RL, GB, MG, and LM participated in the child recruitment and anthropometric measures (data acquisition). CAG and AAR revised DNA extraction and methylation analyses. FJRO and MCR planned and performed the ISM1 analysis in serum of the children. AAR and FJRO performed all data analyses. AAR, FJRO and CAG took part in the interpretation of data. FJRO and CAG wrote the manuscript. All authors took park in the critical revision of the manuscript. CAG, RL and GB obtained funding. All authors approved the final version of the manuscript.

Supplementary Figures

• Download figure
• Open in new tab

Figure S1 related to Figure 2.

Multiple linear regressions were employed with M values as independent variables and each outcome as the dependent variable. Models were adjusted for confounders when necessary. Because percentage methylation is easily interpretable, beta values were employed for the graphical representation of results. A refers to BMI Z-score as dependent variable.

• Download figure
• Open in new tab

Figure S2 related to Figure 3.

Cross-sectional associations between ISM1DNAmethylation status and ISM1 serum levels (ng/mL). Multiple linear regressions were employed with M values as independent variables and each outcome as the dependent variable. Models were adjusted for confounders when necessary. Because percentage methylation is easily interpretable, beta values were employed for the graphical representation of results.