The threshold value for identifying insulin resistance (HOMA-IR) in an admixed adolescent population: A hyperglycemic clamp validated study

ABSTRACT Objectives: To validate the homeostasis model assessment (HOMA) of insulin resistance (IR) as a surrogate to the hyperglycemic clamp to measure IR in both pubertal and postpubertal adolescents, and determine the HOMA-IR cutoff values for detecting IR in both pubertal stages. Subjects and methods: The study sample comprised 80 adolescents of both sexes (aged 10-18 years; 37 pubertal), in which IR was assessed with the HOMA-IR and the hyperglycemic clamp. Results: In the multivariable linear regression analysis, adjusted for sex, age, and waist circumference, the HOMA-IR was independently and negatively associated with the clamp-derived insulin sensitivity index in both pubertal (unstandardized coefficient – B = −0.087, 95% confidence interval [CI] = −0.135 to −0.040) and postpubertal (B = −0.101, 95% CI, −0.145 to −0.058) adolescents. Bland-Altman plots showed agreement between the predicted insulin sensitivity index and measured clamp-derived insulin sensitivity index in both pubertal stages (mean = −0.00 for pubertal and postpubertal); all P > 0.05. The HOMA-IR showed a good discriminatory power for detecting IR with an area under the receiver operator characteristic curve of 0.870 (95% CI, 0.718-0.957) in pubertal and 0.861 (95% CI, 0.721-0.947) in postpubertal adolescents; all P < 0.001. The optimal cutoff values of the HOMA-IR for detecting IR were > 3.22 (sensitivity, 85.7; 95% CI, 57.2-98.2; specificity, 82.6; 95% CI, 61.2-95.0) for pubertal and > 2.91 (sensitivity, 63.6; 95% CI, 30.8-89.1, specificity, 93.7; 95%CI, 79.2-99.2) for postpubertal adolescents. Conclusion: The threshold value of the HOMA-IR for identifying insulin resistance was > 3.22 for pubertal and > 2.91 for postpubertal adolescents.


INTRODUCTION
T he increase in the prevalence of type 2 diabetes mellitus (1,2) and metabolic syndrome (3) in children and adolescents is an important public health concern. Insulin resistance has a key role in the pathogenesis of type 2 diabetes mellitus (1,4) and metabolic syndrome (3,5), which can lead to the development of coronary artery disease (3,6). Hence, a valid, practical, and accessible method of assessing insulin resistance in this age group must be developed to monitor its progression over time, to identify adolescents at risk of developing associated factors, and to establish strategies for preventing and mitigating the transition from normal glucose tolerance to impaired fasting glucose and type 2 diabetes mellitus.
Insulin resistance can be assessed in vivo by several methods. The euglycemic-hyperinsulinemic clamp technique is considered the gold standard for assessing insulin sensitivity/resistance (7,8). However, it is not applicable to large-scale epidemiological studies or clinical practice due to being a complex, invasive, expensive, and time-consuming method (8). The homeostasis model assessment (HOMA) of insulin resistance (IR) is a surrogate marker that estimates insulin resistance based on basal measurements of plasma insulin and glucose (9). It has been widely validated and used in clinical and epidemiological studies of adult populations (10)(11)(12). The HOMA-IR has also been validated as a surrogate to the clamp technique as a measure of insulin resistance in adolescents (13)(14)(15)(16)(17). However, to our knowledge, studies have not separately validated the HOMA-IR in pubertal and postpubertal adolescents. Cross-sectional and longitudinal studies have shown that a significant physiological change in insulin sensitivity occurs during the transition from late childhood throughout adolescence, with increased insulin resistance at the onset of puberty and subsequent normalization towards the end of pubertal development (18,19).
This study aimed to validate the HOMA-IR as a surrogate to the hyperglycemic clamp technique to measure insulin resistance in both pubertal and postpubertal adolescents; and determine the HOMA-IR cutoff values for detecting insulin resistance in both pubertal stages.

Study design and participants
This study used data from the Brazilian Metabolic Syndrome Study (BRAMS), a cross-sectional study conducted in the state of São Paulo, Brazil. The BRAMS studied the insulin resistance in an intentional non-probabilistic sample composed of adolescents aged from 10 to 19 years and 11 months. Out of 1,033 enrolled participants in the BRAMS study, data from 80 adolescents (aged 10-18 years, 40 females) who underwent the hyperglycemic clamp technique were analyzed. The adolescents were recruited from public schools and the University of Campinas Teaching Hospital. The following exclusion criteria were applied: prepubertal children due to the small sample size, pregnancy, use of either systemic corticosteroids or drugs with hypoglycemic properties, malnutrition, hepatopathy, nephropathy, metabolic disorders (e.g., hypothyroidism, hyperthyroidism, and type 1 and 2 diabetes), genetic syndrome diagnosis, and delayed neuropsychomotor development.
The study was approved by the Research Ethics

Clinical evaluation
Pubertal development was assessed by self-assessments (20) according to Tanner's criteria (21). The self-assessment method in the BRAMS study has been reported in detail elsewhere (22). Participants were divided into two groups: pubertal (Tanner II-IV) and postpubertal (Tanner V).

Anthropometric and body composition assessments
The body mass index (BMI) was calculated as weight (kilograms) divided by height in meters squared. The BMI-for-age z-score was calculated using the Epi Info version 3.5.2 software (Centers for Disease Control and Prevention, Atlanta, Georgia, USA). The nutritional status was defined using the Centers for Disease Control and Prevention criteria (23). Waist circumference was measured at the midpoint between the lower margin of the last palpable rib and the top of the iliac crest (24). The amount of lean body mass was determined using tetrapolar bioimpedance (Biodynamics, model 310, Shoreline, Washington, USA) (25).

Insulin resistance assessment
Participants underwent a 2-hour hyperglycemic clamp (with blood glucose acutely raised and maintained at approximately 225 mg/dL; to convert to millimoles per liter, multiply by 0.0555) according to the protocol previously described by Arslanian (8). The insulin sensitivity index (ISI) from the hyperglycemic clamp technique was calculated as the mean exogenous glucose infusion rate from 60 to 120 minutes of the clamp technique minus the urinary glucose excretion, divided by the mean insulin concentration of five determinations during the same time period, and it was then corrected for lean body mass (LBM) (ISI LBM ; Threshold value for HOMA-IR in adolescent Arch Endocrinol Metab. 2023;67/1 milligrams of glucose infused per kilogram of lean body mass per minute, multiplied by 100) (8,26). The HOMA-IR index was calculated as the product of the fasting plasma insulin level (in milliunits per liter) and the fasting plasma glucose level (in millimoles per liter), divided by 22.5 (9).

Statistical analysis
The Shapiro-Wilk test was used to check the distribution of variables. Data are reported as the mean ± standard deviation or median (interquartile range). Relationship between two variables was evaluated with the Spearman's correlation coefficient. We used multivariable linear regression models to evaluate the associations between the HOMA-IR and the clamp-derived ISI (independent variables, sex, age, and HOMA-IR in Model 1 and the variables in Model 1 plus waist circumference in Model 2). Residuals were evaluated for normal distribution by the Shapiro-Wilk test. Preliminary prediction models demonstrated non-normality of the residuals and clamp-derived ISI was therefore transformed to the logarithmic scale. The agreement between the predicted ISI and measured clamp-derived ISI LBM was evaluated by Bland-Altman plots (27). For each plot, average bias and 95% limits of agreement were estimated, with results analyzed further by a One-sample T-test to assess the significance of any bias between measured clamp-derived ISI LBM and predicted ISI (28). The bias represents the mean difference between the two methods (29). For a method to be considered of good agreement, the mean differences should not different from zero (29). The area under the receiver-operating characteristic (ROC) curve (AUC), sensitivity, specificity, positive predictive value, and negative predictive value were calculated to evaluate the accuracy of the HOMA-IR for detecting insulin resistance in both pubertal and postpubertal adolescents. The optimal cutoff points of the HOMA-IR were obtained from the Youden index, which is defined as the maximum sensitivity + specificity − 1 (30). To classify insulin resistance for the ROC curve analysis, the cutoff of the insulin sensitivity index corrected for lean body mass (lower 10th percentile) derived from the normal-weight group was used (cutoff, < 0.08 mg/kg LBM /min per mU/L for the pubertal and < 0.07 mg/kg LBM /min per mU/L for the postpubertal adolescents). Statistical analyses were performed using IBM SPSS Statistics for Windows, version 20.0 (IBM Corporation, Armonk, NY, USA) and MedCalc for Windows, version 18.5 (MedCalc Software, Ostend, Belgium). In all statistical tests, P values < 0.05 were considered significant. Table 1 shows clinical, anthropometric, and biochemical data. Data are presented as the mean ± standard deviation or as the median (interquartile ranges).

Relationship of the HOMA-IR index and clampderived insulin sensitivity index
The HOMA-IR index showed a strong correlation with the clamp-derived ISI in both the pubertal (r = −0.77; P < 0.001) and postpubertal (r = −0.83; P < 0.001) adolescents.

Association between the clamp-derived insulin sensitivity index and HOMA-IR index
In the multivariable linear regression analysis, adjusted for sex and age, the HOMA-IR index was independently and negatively associated with the clamp-derived ISI in both the pubertal and postpubertal adolescents ( Table 2). The HOMA-IR index remained negatively associated with the clamp-derived ISI, in both the pubertal and postpubertal adolescents, even after further adjustment for waist circumference ( Table 2).

Agreement between the predicted ISI and measured clamp-derived ISI LBM
The Bland-Altman plot shows a bias equal to zero (Figure 1). The mean differences between predicted ISI and measured clamp-derived ISI LBM were not significantly different from zero in both pubertal (t(36) = 0.00, P > 0.05) and postpubertal (t(42) = 0.00, P > 0.05) adolescents.

DISCUSSION
This study of pubertal and postpubertal adolescents indicates that the HOMA-IR index is strongly related with the clamp-derived ISI in both pubertal stages. To our knowledge, this is the first study to explore the association between the HOMA-IR index and the clamp-derived ISI in these two pubertal stages separately. We found that the HOMA-IR index was negatively associated with clamp-derived insulin sensitivity, even after adjustment for waist circumference, in both pubertal and postpubertal adolescents. Additionally, Bland-Altman plots showed agreement between the predicted ISI and measured clamp-derived ISI LBM in both pubertal stages. Finally, we found that the HOMA-IR index was capable of accurately detecting insulin resistance in both pubertal and postpubertal adolescents. The cutoff points for detecting insulin resistance using the HOMA model were different between the pubertal (>3.22) and postpubertal (>2.91) adolescents.
Studies have compared the HOMA-IR index with clamp-derived measures in pediatric populations (13)(14)(15)(16)(17), and as well as in this study, they found a significant correlation between the two methods. However, previous studies have not demonstrated whether the HOMA-IR is capable of separately estimating insulin resistance in pubertal and postpubertal adolescents. Although Gungor and cols. (14) reported correlations for pubertal adolescents, they defined pubertal adolescents as Tanner stages II to V. However, in our study, adolescents were considered as pubertal if they presented Tanner stages  (31,32), whereas the cutoff for postpubertal adolescents (>2.91) was higher than previously reported value (2.7) (32). These differences may be related to BMI differences, population age, although they are primarily related to the accuracy of the methodology used for determining the cutoff point (hyperglycemic clamp versus oral glucose tolerance test (31) or 95th percentile of the HOMA-IR (32)).
A potential limitation of the current study is its cross-sectional design, which does not allow for inferences of causality. Another limitation is the use of the hyperglycemic clamp technique to evaluate insulin resistance. Although the hyperglycemic clamp technique is not the gold standard for estimating insulin resistance, studies comparing this technique with the euglycemic-hyperinsulinemic clamp technique (gold standard for quantifying insulin resistance) reported an excellent correlation between both clamp techniques in children and adolescents (14,26,33). The evaluation of sexual maturation was performed via self-assessments (20,34) to increase the participation rate, due to privacy concerns, cultural, and emotional factors. Studies that evaluated the agreement between self-assessment sexual maturation and physical examination performed by a physician (20,34) suggest that the self-assessment can be used in epidemiologic studies for evaluating sexual maturation when the physician exam is impossible. Also, this study did not report race/ethnicity-stratified results, because Brazil has one of the most admixed populations and this distinction is unfeasible.
In summary, the threshold value of the HOMA-IR for identifying insulin resistance was > 3.22 for pubertal and > 2.91 for postpubertal adolescents. The HOMA-IR is a low-cost approach with potential clinical and epidemiological applications and these cutoff points can improve the detection and control of metabolic diseases in pubertal and postpubertal adolescents.