Introduction

Although obesity and insulin resistance (IR) are closely linked, BMI explains only 22% of the variance of IR in the general population [1]. Individuals with the metabolic syndrome who do not meet the criteria for central obesity can still be insulin resistant, and exhibit a similarly increased risk of coronary artery disease [1, 2]. Several studies have shown strong correlations between dysregulated lipid metabolism and IR, with some reports proposing impaired mitochondrial handling of lipids and consequent accumulation of intracellular lipid-derived metabolites, such as diacylglycerols and ceramides [35], whereas others implicate increased fatty acid flux and accumulation of incompletely oxidised lipid metabolites in the mitochondria [68]. Protein and amino acid metabolism are also dysregulated in obese and insulin-resistant states or type 2 diabetes [914], but the extent to which this contributes to IR in normal-weight individuals is incompletely understood. In addition to changes in circulating lipids and amino acids, obese and insulin-resistant individuals exhibit changes in a wide array of circulating hormones and inflammatory mediators [14], some of which may make a direct contribution to IR [7, 8].

We have recently used comprehensive metabolic profiling (‘metabolomics’) to identify multiple metabolites from different biochemical classes that differ in obese compared with lean humans from the south-eastern USA. Multivariate statistical analysis showed that the principal component (cluster of metabolites) most strongly associated with IR (as measured by HOMA index) was one comprising branched-chain amino acids (BCAAs) and related metabolites [14]. A subsequent cross-sectional study in sedentary and dyslipidaemic individuals of varying body weight confirmed a strong relationship between a similar BCAA-related metabolite cluster and the insulin sensitivity index [15]. The source of increased BCAAs in the blood of insulin-resistant individuals is unknown; contributing factors could include increased protein intake in obese individuals or changes in rates of protein turnover or amino acid catabolism. A mechanistic link between elevated BCAAs and IR was suggested in studies in rats, in which supplementation of a high-fat diet with BCAAs caused an increase in circulating BCAAs and their metabolites, and contributed to development of IR [14].

The foregoing summary highlights the need for deeper understanding of metabolic, hormonal and inflammatory factors contributing to the development of IR. A particular opportunity may reside in studies of development of IR in Asian populations. The epidemic of type 2 diabetes mellitus is occurring most rapidly in Asia [16], and Asians often develop IR at relatively low BMI [17, 18]. Singapore comprises three major ethnic groups (Chinese, Malaysian and Asian-Indian) living in a small geographical area. The country has become completely urbanised over the last 50 years, accompanied by a rapid increase in the prevalence of diabetes mellitus, from 4.7% in 1984 [19] to >8% in 1998 [20]. Compared with Chinese, Asian-Indians in Singapore have a much higher prevalence of type 2 diabetes mellitus, and exhibit greater IR, pathologies not accounted for by differences in BMI [20, 21]. Asian-Indians also have greater abdominal adiposity (as measured by waist circumference) [22], lower HDL-cholesterol and a higher risk of coronary artery disease, but lower BP than Chinese [23]. It is thus possible that differences in the phenotype associated with IR are related to differences in the causes of IR in the different ethnic groups, but this idea has never been explored in a systematic fashion.

In the current study, we have measured a wide array of metabolic, hormonal and inflammatory mediators in Asian-Indian and Chinese individuals from Singapore with relatively normal body weight (mean BMI approximately 24 kg/m2), separated into upper and lower tertiles of insulin sensitivity. By studying Singaporean residents, we minimised demographic and societal variables that might influence metabolic control and insulin sensitivity. With this approach, we sought to address the following key questions: (1) which of the elements of the metabolic profile associated with obesity is independently associated with IR and (2) when controlled for BMI, which metabolic features differ in the two ethnic groups living in Singapore? In this study population, we find a lack of association between IR and traditional IR biomarkers such as inflammatory mediators and fatty acids. IR is instead associated with changes in amino acid metabolism, visceral adiposity, the levels of two adipokines, leptin and adiponectin, and IGF-I-binding protein (IGFBP)-1.

Methods

Participant recruitment and characteristics

We invited 10,445 individuals from four previously conducted population-based cross-sectional surveys carried out in Singapore (1982–1998) [20, 22, 2426] to participate in a repeat examination from 2004 to 2007. All studies were a random sample of Chinese, Malay and Asian-Indian individuals from the Singaporean population. Individuals who had died by the time of follow-up were excluded (n = 517). Also excluded were six individuals who had emigrated and 85 individuals who had errors in the records for the Identity Card number. All contactable individuals were invited to attend a health examination for additional tests and collection of biological specimens shortly after the home visit. A total of 7,742 participants completed the questionnaire, of whom 5,157 also attended the health examination.

Data on demographics and lifestyle factors, as well as medical history (including physician-diagnosed hypertension, diabetes mellitus and hyperlipidaemia) were collected using interviewer-administered questionnaires. Dietary intake was assessed using a semi-quantitative food frequency questionnaire, as described before [2730].

For the health examination, participants were examined in the morning following a 10 h overnight fast. Weight was measured in light clothing using the same digital scale (SECA, model 782 2321009; Vogel and Halke, Hamburg, Germany) for all participants. Waist circumference was measured with a non-elastic tape at the mid-point between the costal margin and the iliac crest in the mid-axillary line. Hip circumference was measured at the level of the greater trochanter. Two readings of BP were taken from participants after a 5 min rest, using an automated BP monitor (Dinamap Pro100V2; Criticon, Norderstedt, Germany). Hypertension was defined as systolic BP >140 mmHg or diastolic BP >90 mmHg or history of hypertension or currently taking anti-hypertensive medications.

Participant selection for metabolic profiling

We studied only men in order to reduce the potential confounding effects of sex as well as metabolic alterations occurring at different phases of the menstrual cycle. Furthermore, we chose to study only Chinese and Asian-Indian men, as these were the ethnic groups that show the greatest differences in IR in our population [21, 22]. In order to examine the metabolic and hormonal correlates of IR as opposed to obesity, we first separated Asian-Indian and Chinese men into the upper and lower tertiles of IR, as measured by HOMA. These opposing tertiles corresponded to HOMA indices of ≤1.06 and ≥1.93. We next excluded all those individuals who were taking medications for diabetes or hyperlipidaemia, those who did not provide a urine specimen, and those with missing height and weight information. Of the 56 Asian-Indians with low HOMA index, 30 had a BMI between 21 and 29 kg/m2 (non-obese). For each of these individuals, we selected up to three individuals from the other three groups, with BMI within 0.3 kg/m2. Of all the participants available we were able to identify 97 Chinese with low HOMA, 83 Chinese with high HOMA and 53 Asian-Indians with high HOMA with BMI within 0.3 kg/m2 of an Asian-Indian individual with low HOMA, giving a total sample size of 263 men.

Measurement of insulin sensitivity

Insulin sensitivity was measured by HOMA [31]. In what follows, we will refer to Asian-Indian or Chinese individuals in the upper tertile for HOMA interchangeably as ‘high HOMA’ or ‘insulin resistant’.

Biochemical measurements

Venous blood was drawn and collected in plain and fluoride oxalate tubes in the morning following a 10 h overnight fast. Some biochemical analyses on blood were carried out at the National University Hospital Referral Laboratory in Singapore. Serum total cholesterol, triacylglycerol and HDL-cholesterol were measured using an automated autoanalyser (ADVIA 2400; Bayer Diagnostics, Leverkusen, Germany). LDL-cholesterol levels were calculated using the Friedewald formula. Plasma glucose was also assayed using enzymatic methods (ADVIA 2400, Bayer Diagnostics) using blood collected in fluoride oxalate tubes.

All other metabolic, hormonal and inflammatory marker assays were performed at the Sarah W. Stedman Nutrition and Metabolism Center, Duke University. Commercial RIAs were used to measure leptin (Millipore, St Charles, MO, USA) and IGF-I (Beckman-Coulter, Webster, TX, USA). Searchlight chemiluminescent technology (Thermo Pierce-Endogen, Woburn, MA, USA) was used for measuring three IGFBPs (IGFBP1, IGFBP2, IGFBP3) using a custom multiplex assay. Adiponectin was measured via technology from Meso Scale Discovery (Gaithersburg, MD, USA), as was a nine plex cytokine panel for granulocyte–macrophage colony-stimulating factor, IFN-γ, IL-10, IL-12p70, IL-1β, IL-2, IL-6, IL-8 and TNF-α. Total NEFAs, ketones and β-hydroxybutyrate (reagents from Wako, Richmond, VA, USA) and lactate (Roche, Indianapolis, IN, USA) were measured on a Hitachi 911 autoanalyser, as described previously [14].

Acylcarnitines and amino acids were analysed in plasma by tandem MS (MS/MS), and urine organic acids by GC/MS, as described previously [14, 3235]. All MS analyses employed stable-isotope dilution with internal standards from Isotec (St Louis, MO, USA), Cambridge Isotope Laboratories (Andover, MA, USA) and CDN Isotopes (Pointe-Claire, QC, Canada) [14, 33].

Ethics approval was obtained from two Institutional Review Boards (National University of Singapore and Singapore General Hospital). Informed consent was obtained before the conduct of the study.

Statistical analyses

Statistical analysis was carried out using STATA 8.0 (StataCorp LP; College Station, TX, USA and SAS v9.1.3, Cary, USA). Categorical variables were compared using the χ 2 test. Continuous variables were compared using ANOVA. We first used an F test to identify those variables for which at least one of the four groups showed a statistically significant difference from the others. Subsequently, the levels of metabolites in high- and low-HOMA groups were compared using one-way ANOVA with a Bonferroni multiple-comparison test, with and without adjustment for age and BMI. p < 0.05 was taken as statistically significant.

We also employed principal components analysis (PCA) to reduce the large number of correlated variables into uncorrelated factors. PCA was constructed on the full set of study individuals using all measured metabolites and hormones, without prior knowledge of ethnicity or HOMA. Factors with an eigen value ≥1.0 were identified and varimax rotation performed to produce interpretable factors. Metabolites with a factor load ≥0.4 were reported as composing a factor. Scoring coefficients were constructed from the initial group and used to calculate factor scores for each individual (weighted sum of the standardised metabolites within that factor, weighted on the factor loading for each metabolite). Statistical analyses of these metabolite factors was similar as that for the individual metabolites, namely an F test was first used to identify factors that were significantly different between groups, and subsequently factor levels in high- and low-HOMA groups were compared using ANOVA, with and without adjustment for age and BMI.

Results

Demographic and clinical characteristics of participants

Because we measured a total of 100 separate analytes in this study, Tables 1, 2, 3 and 4 summarise only those analytes for which at least one of the four groups showed a statistically significant difference from the others with p < 0.05. Data for analytes that did not differ significantly among groups can be found in Electronic supplementary material (ESM) Tables 1 and 2. Table 1 presents demographic and clinical characteristics of the Chinese and Asian-Indian study participants, stratified by IR. Individuals from the top tertile of IR in both ethnic groups had HOMA indices three to four times higher than those from the lowest tertile. These differences were driven almost exclusively by elevations in fasting insulin levels in the high-HOMA groups. Despite our attempts to match for BMI (mean BMI values among the four groups were within a tight range of 23.89–24.95 kg/m2), BMI was slightly higher in individuals in the high-HOMA groups. Hip circumference was identical between the groups, whereas waist circumference, and consequently WHR, varied significantly between high- and low-HOMA participants in both ethnic groups. The levels of HDL-cholesterol were significantly lower in both groups of Asian-Indian individuals than in Chinese individuals, and HDL-cholesterol levels were lower in the high- vs low-HOMA individuals within both ethnic groups, whereas LDL-cholesterol levels were not associated with HOMA tertile or ethnic origin. Triacylglycerol levels were clearly elevated in the high-HOMA individuals in both the Chinese and Asian-Indian groups (Table 1), whereas NEFA and ketone levels did not differ (Table 3 and ESM Table 1). Finally, BP (both systolic and diastolic) was modestly higher in the high-HOMA vs low-HOMA individuals within both ethnic groups (Table 1).

Table 1 Demographic and clinical characteristics
Table 2 Dietary intake and physical activity
Table 3 Metabolite concentrations for analytes that differed among study groups
Table 4 Hormone profiles

Dietary intake and physical activity of participants

Table 2 shows the dietary intake and physical activity in these individuals. As a proportion of total energy intake, Asian-Indians consumed less protein and more saturated fat than Chinese. However, no dietary differences were noted between high- and low-HOMA individuals in either ethnic group, including no difference in protein intake. There was a trend towards less physical activity in individuals with high HOMA that was not statistically significant (p = 1.0 between high- and low-HOMA groups in both Chinese and Asian-Indians).

Blood and urinary metabolites in study participants

Table 3 shows the concentrations of individual metabolites in the various groups. Chinese individuals with high HOMA had higher levels of several amino acids than Chinese individuals with low HOMA. These included alanine, proline, valine, leucine/isoleucine, phenylalanine, tyrosine, glutamate/glutamine and ornithine. With the exception of phenylalanine, all these differences remained statistically significant when adjusted for age and BMI. Associations in the same direction were observed in Asian-Indians but reached statistical significance only for alanine, tyrosine and glutamate/glutamine.

We also measured a group of 20 urinary organic acids in all of the participant groups (Table 3). Lactate and pyruvate levels were increased in urine in the high-HOMA individuals of both ethnic groups; these changes were accompanied by significant increases in blood lactate levels as a function of HOMA, again in both ethnic groups. Interestingly, urinary isobutyrylglycine was significantly lower in the high-HOMA vs low-HOMA group in Asian-Indians but not in Chinese, and isovalerylglycine showed a similar trend. Methylmalonate, methylsuccinate and hydroxymethylglutarate levels were all higher in urine of Asian-Indian compared with Chinese individuals, but did not change significantly with HOMA in either ethnic group. Finally, among 45 acylcarnitine species measured by MS/MS, C10:1 and C8:1-dicarboxyl were low in insulin-resistant individuals, but only in Asian-Indians and not in Chinese. C3 and C5 acylcarnitines, previously shown to be elevated in association with the increase in BCAAs in obese vs lean individuals from the south-eastern USA [14], tended to be higher in the high-HOMA group in both the Asian-Indian and Chinese men, but this change did not achieve statistical significance (ESM Table 1).

Hormone and cytokine profiling

Table 4 shows the hormonal profile of the various groups. In both ethnic groups, IR was associated with elevated plasma leptin and reduced adiponectin, but in addition leptin levels were clearly elevated in the Asian-Indians compared with the Chinese, regardless of HOMA group. IGF-I levels were negatively associated with HOMA in the Chinese, but this was no longer statistically significant after adjustment for age and BMI. More strikingly, IGF-I levels were higher in Chinese than in Asian-Indian individuals, irrespective of the degree of IR. IGFBP1 levels were sharply decreased as a function of HOMA in both ethnic groups, but overall levels were not different between Chinese and Asian-Indian individuals. Among cytokines, TNF-α levels were elevated in high-HOMA vs low-HOMA quartiles, but these changes only reached statistical significance for Chinese individuals and were attenuated after adjustment for age and BMI.

PCA

PCA was applied to all analytes measured and all study participants, and identified 26 factors grouping in biologically plausible clusters (ESM Table 3). Factor 4 was composed of the BCAAs leucine/isoleucine, as well as phenylalanine, tyrosine and methionine. Although not meeting our a priori criteria for inclusion based on factor loading (factor load ≥ 0.4), factor 4 also included glutamate/glutamine (0.29), alanine (0.30), proline (0.22) and valine (0.22). Consistent with the results for the individual metabolites, factor 4 scores differed significantly across groups (p value for overall F test 0.003), and were significantly higher in individuals with high-HOMA compared with individuals with low HOMA, in both Chinese (p = 0.005) and Asian-Indians (p = 0.02), even after adjusting for age and BMI (p = 0.04 and 0.02, respectively) (Table 5). In addition, factor 4 was significantly correlated with HOMA across all individuals (r 2 = 0.18, p = 0.003).

Table 5 PCA-derived factor scores

Factor 8, composed of pyruvate, lactate and arginine, was also significantly different between groups, with higher levels in high-HOMA individuals in both ethnic groups. This finding mirrored results for the individual analytes of this factor. Factor 13, composed of isobutyrylglycine and isovalerylglycine, was significantly lower in the high-HOMA vs low-HOMA group in Asian-Indians but not in Chinese, again similar to the results for individual metabolites in the factor. Factor 21, composed of IGF-I and IGFBP3 was significantly lower in high-HOMA Chinese individuals vs low-HOMA Chinese individuals, and was not significant in Asian-Indians. No other factors were significantly different between high- and low-HOMA individuals.

Discussion

Using a comprehensive metabolic profiling approach, we have found a clear association between a subset of amino acids and IR in Chinese and Asian-Indian individuals living in Singapore. In a previous study that applied the same tools to a group of obese (BMI 37 kg/m2) vs lean (BMI 23 kg/m2) individuals residing in the southeastern USA, striking increases in amino acid levels were noted in the obese individuals, including BCAAs and several of their metabolites [14]. The current findings strongly support the idea that these changes are related to IR rather than obesity per se, since we observed changes in many of the same individual amino acids in high vs low tertiles of IR in a group of Chinese individuals with a mean BMI of approximately 24 kg/m2. Although the number of Asian-Indian individuals available limited the power of this study to detect associations between individual metabolites and IR in this ethnic group, BCAAs and related amino acids trended in the same direction in Asian-Indians as observed in the Chinese, and the PCA-derived factor 4 composed of BCAAs and related amino acids was associated significantly with IR in both ethnic groups.

As essential amino acids, BCAA levels in blood are controlled by dietary intake, their rate of use for the anabolic processes of protein biosynthesis and cell growth, their rate of catabolism through transamination and the branched-chain keto acid dehydrogenase complex, and the rate of protein turnover/hydrolysis. In our prior study, we reported that obese individuals consumed more protein than lean individuals [14]. In the present study, whereas dietary protein intake differed between ethnic groups, it did not differ between individuals with high or low HOMA. Possible explanations for the increased levels of BCAAs in individuals with high HOMA in the current study could include increased protein turnover and/or reduced rates of BCAA catabolism. Support for the former idea comes from studies showing higher rates of leucine appearance in obese compared with lean individuals, and even higher rates of protein turnover in obese diabetic individuals compared with obese non-diabetic individuals [11]. On the other hand, a separate group reported no change in overall rate of leucine appearance in obese compared with lean individuals, and an actual decrease in appearance of leucine from forearm and adipose tissue in the obese individuals [36]. A major difference in the two studies was that individuals were studied in the post-absorptive state or after short-term fasting in the former case, and after 22 h of fasting in the latter case. Obesity has been reported to be associated with changes in the abundance of BCAA-catabolising enzymes in adipose tissue [37], but the contribution of these changes to plasma BCAA levels has not been demonstrated. Interestingly, in the current study, we saw increases in some metabolites reflective of BCAA catabolism (e.g. glutamate/glutamine, alanine, C3 and C5 acylcarnitines), but decreases in others (urinary isobutyrylglycine and isovalerylglycine). We are not aware of any data on protein turnover or expression of BCAA-catabolic enzymes in individuals matched for a relatively low body weight but differing in insulin sensitivity. Thus, further work will be required to fully address the mechanisms underlying the increases in amino acid levels in insulin-resistant individuals of the type studied here.

Despite similar BMI, insulin-resistant and insulin-sensitive individuals in the current study had markedly different waist circumferences. This suggests that despite similar BMI, individuals who were insulin resistant in this study may have had increased central adiposity, which is often associated with reduced lean body mass, although neither of these variables was directly measured in the current study. Central adiposity accompanied by reduced lean body mass is a condition that has been termed ‘sarcopenic obesity’ [38], and its presence is consistent with the strongly elevated plasma leptin levels observed in the individuals with high HOMA from both ethnic groups in the current study. Assuming that the increased levels of select amino acids observed in insulin-resistant individuals in the current study are indicative of increased protein turnover, this profile could be consistent with sarcopenia. Also consistent with this idea, IR is independently correlated with poor muscle strength and with accelerated loss of muscle strength and quality in older individuals with type 2 diabetes [3942], although such relationships are not apparent in all studies [43, 44].

Downregulation of the growth hormone axis has also been implicated in the link between IR and sarcopenia. Levels of IGF-I were reduced in insulin-resistant individuals in the present study, similar to our prior findings in obese, insulin-resistant North American individuals [14]. This could potentially result in decreased anabolic use of amino acids, including BCAAs, for protein synthesis and cell growth, contributing to their rise in the circulation. Interestingly, Chinese individuals had higher IGF-I levels than Asian-Indian individuals, possibly contributing to the particular susceptibility of the latter group to develop metabolic diseases. We also considered the possibility that increased inflammation associated with obesity and IR may result in accelerated muscle catabolism [45, 46]. However, the levels of pro-inflammatory cytokines were only minimally elevated in our insulin-resistant individuals.

When BMI is controlled as in the current study, no correlation is evident between HOMA and either NEFA levels or levels of lipid-derived metabolites such as acylcarnitines or ketones. Increases in NEFAs have been reported in obese and insulin-resistant individuals in many prior studies, leading to extensive studies on the potential ‘lipotoxic’ effects of these agents in development of IR [37]. It should be noted that the absence of a rise in circulating NEFA in the insulin-resistant individuals in the current study does not mean that rates of fatty acid turnover were unchanged [47, 48], and flux measurements would be needed to fully address this point. Another striking finding of our study is the dramatic lowering of IGFBP1 in the high-HOMA individuals in both ethnic groups. Insulin is known to suppress production of IGFBP1 in the liver [49]. These findings suggest that the elevated plasma insulin in insulin-resistant individuals was largely effective for suppression of lipolysis and NEFA levels, as well as hepatic IGFBP1 production. However, these same elevated insulin levels apparently failed to control the levels of BCAAs and related amino acids. Whether this is indicative of selective IR at the level of suppression of protein turnover and amino acid catabolism or whether the reduced IGF-I levels play a more prominent and primary role is a matter that will require further investigation.

Several limitations of the study should be noted. First, we have only studied men, and the relevance of these findings to women is not known at this time. Second, we have only studied these individuals in the fasting or post-absorptive state. Third, HOMA is an imprecise method for assessing IR. However, individuals from opposing tertiles of HOMA index had clear differences in the levels of blood lipids (triacylglycerol and HDL-cholesterol) and adiponectin, supporting the idea that these individuals had true differences in insulin sensitivity.

In conclusion, significant evidence supports the role of NEFAs and inflammation in the pathogenesis of obesity-associated IR. Our study complements this knowledge, by showing that when IR develops at relatively low body weight it is no longer associated with NEFA and inflammation but instead with a metabolic signature of elevated BCAAs and related metabolites. It remains to be determined if the rise in BCAAs and related metabolites is an early manifestation of IR caused by other factors, or is actually one of the primary causes of the insulin-resistant state. Support for a primary role of amino acids includes recent studies in animals showing that feeding of BCAAs causes IR, particularly in the context of a high-fat diet [14]. On the other hand, hepatic steatosis is an example of a factor that could trigger an early state of IR [50], leading to a secondary rise in BCAAs and their metabolites. It will also be important to determine if insulin-resistant individuals such as those studied here with relatively low body weight actually do have enhanced protein catabolism and other markers of reduction of lean body mass/sarcopenia. If so, and given that increased physical activity is known to attenuate muscle loss during ageing or weight loss [51], future work might be directed at investigation of the impact of physical activity on amino acid profiles and IR in non-obese individuals. Regardless of the outcome of these additional studies, the current work has shown that BCAAs and related metabolites are associated with IR even in the absence of obesity, and suggest that dysregulation of amino acid metabolism may be an early event in the progression to more severe IR and type 2 diabetes.