Association of ADRB2 rs1042713 with Obesity and Obesity-Related Phenotypes and Its Interaction with Dietary Fat in Modulating Glycaemic Indices in Malaysian Adults

Gene-diet interaction studies have reported that individual variations in phenotypic traits may be due to variations in individual diet. Our study aimed to evaluate (i) the association of ADRB2 rs1042713 with obesity and obesity-related metabolic parameters and (ii) the effect of dietary nutrients on these associations in Malaysian adults. ADRB2 genotyping, dietary, physical activity, anthropometric, and biochemical data were collected from 79 obese and 99 nonobese individuals. Logistic regression revealed no association between ADRB2 rs1042713 and obesity (p=0.725). However, the carriers of G allele (AG + GG genotypes) of rs1042713 were associated with increased odds of insulin resistance, 2.83 (CI = 1.04–7.70, adjusted p=0.042), in the dominant model, even after adjusting for potential confounders. Obese individuals carrying the G allele were associated with higher total cholesterol (p=0.011), LDL cholesterol levels (p=0.008), and total cholesterol/HDL cholesterol ratio (p=0.048), compared to the noncarriers (AA), even after adjusting for potential confounders. Irrespective of obesity, the carriers of GG genotype had significantly lower fasting glucose levels with low saturated fatty acid intake (<7.3% of TE/day) (4.92 ± 0.1 mmol/L vs 5.80 ± 0.3 mmol/L, p=0.011) and high intake of polyunsaturated fatty acid:saturated fatty acid ratio (≥0.8/day) (4.83 ± 0.1 mmol/L vs 5.93 ± 0.4 mmol/L, p=0.006). Moreover, the carriers of GG genotype with high polyunsaturated fatty acid intake (≥6% of TE/day) had significantly lower HOMA-IR (1.5 ± 0.3 vs 3.0 ± 0.7, p=0.026) and fasting insulin levels (6.8 ± 1.6 µU/mL vs 11.4 ± 2.1 µU/mL, p=0.036). These effects were not found in the noncarriers (AA). In conclusion, G allele carriers of ADRB2 rs1042713 were associated with increased odds of insulin resistance. Obese individuals carrying G allele were compromised with higher blood lipid levels. Although it is premature to report gene-diet interaction on the regulation of glucose and insulin levels in Malaysians, we suggest that higher quantity of PUFA-rich food sources in regular diet may benefit overweight and obese Malaysian adults metabolically. Large-scale studies are required to replicate and confirm the current findings in the Malaysian population.


Background
Obesity and related chronic diseases have become the leading causes of morbidity and mortality worldwide. Longterm consumption of energy dense food coupled with sedentary lifestyle are the main contributors to the development of obesity and related complications [1]. Multiple genetic loci determined by genome-wide association studies (GWAS) have been found to be associated with increased susceptibility to obesity, diabetes, and dyslipidemia [2]. Earlier studies have reported on the adrenergic receptor system for its role in the stimulation of thermogenesis and in activating lipid mobilization from the fat stores. With respect to this, the beta-2 adrenoceptor gene is a candidate gene since it is the dominating lipolytic receptor in the human white adipose tissue [3]. It stimulates lipid mobilization through lipolysis in adipocytes and regulates body fat accumulation and energy expenditure [4].
A strong association between obesity and a single nucleotide polymorphism located at codon 16 substituting arginine for glycine (rs1042713/Arg16Gly) has been reported by Large et al. [5]. Masuo et al. reported that insulinresistant subjects had higher frequencies of the G allele of rs1042713 [6]. Total body fat mass and blood pressure levels were higher in nonobese and nonhypertensive men with G allele in the Japanese population. e authors speculated that insulin resistance could, in part, be determined by the genetic variant of the beta-2 adrenoceptor gene and that polymorphism and higher plasma adrenaline could increase insulin resistance, adiposity, and high blood pressure in their subjects. us, the G allele could lead to heightened sympathetic nerve activity, insulin resistance, and higher blood pressure and adiposity in nonobese and nonhypertensive individuals. Other studies have adequately reported that insulin resistance was strongly associated with heightened sympathetic nerve activity [7]. In other words, adrenergic receptor defects lead to the sympathetic nervous system over activity that may play a role in the development of insulin resistance, hypertension, and obesity [8].
A study on the Swedish population reported that ADRB2 rs1042713 was significantly associated with elevated central body fat, systolic blood pressure, serum leptin, and triglyceride levels but not with obesity [9]. Studies from Saudi Arabia reported significant association between ADRB2 rs1042713 polymorphism and the development of insulin resistance, dyslipidemia, overweight, and obesity [10,11]. However, findings from Asian populations (Japanese and Korean) reported negative association between obesity and ADRB2 gene polymorphisms. Moreover, these studies did not investigate interaction of gene variants with dietary nutrients [12,13].
Individual or population differences in the development of obesity-related metabolic diseases may result not only from genetic variation but may also be the modulatory effect of dietary nutrients on gene and gene variants [14]. In the Malaysian population, relatively little is known with respect to the interaction between dietary nutrients and ADRB2 gene variations on obesity, insulin resistance, and glucose homeostasis. It is believed that early identification of the candidate gene variants and their interaction with diet may allow for the provision of good quality personalised dietary recommendation to achieve effective weight loss and reduction in metabolic risk factors [15]. Although a couple of intervention studies have reported lipid outcomes associated with ADRB2 rs1042713, there is no such study done in Malaysian adults. Since there has been conflicting results with respect to ADRB2 rs1042713 between Asian, Caucasian, and Arabic populations, the current study on the Malaysian population is valuable and will shed light and add to the existing evidence on gene variants and phenotypic outcomes and influence of diet on the latter. To the best of our knowledge, this is the first study in Malaysian adults that investigates the interaction between ADRB2 rs1042713 and dietary nutrients on obesity-related metabolic traits. is study is nested in a broader study investigating the association of single nucleotide polymorphisms in genes that have widely been reported to influence obesity and obesityrelated metabolic disorders in human individuals. In an earlier publication, we have reported that FTO rs9930506 may interact with dietary protein and Vitamin E intake and modulate hsCRP levels in our Malaysian participants [16].
e aim of the current study was to evaluate (i) the effect of ADRB2 rs1042713 on obesity and obesity-related anthropometric and blood biochemical parameters and (ii) the influence of diet on the association between ADRB2 rs1042713 and obesity phenotypes, in Malaysian adults.

Study Design.
is cross-sectional study was conducted from 2014-2017 on Malaysian adults aged between 18 and 74 years. e study investigated (i) the association between ADRB2 rs1042713 with obesity and insulin resistance; (ii) the association between ADRB2 rs1042713 and phenotypes in obese and nonobese individuals; and (iii) the interaction between ADRB2 rs1042713 and dietary nutrients on phenotypic traits.

Participant Selection.
Detailed information on the study design and methods can be found in our earlier publication [16]. erefore, with respect to assessment of anthropometric parameters, dietary nutrients analysis, and physical activity level assessment, we refer the readers to our earlier publication [16].

Blood Collection and Biochemical
Analysis. 10-12 hour fasting blood was collected into vacutainer tubes containing fluoride oxalate for plasma glucose analysis and vacutainer tubes with clot activator and gel (Becton Dickinson, Oxford, United Kingdom) for serum lipid profile (total cholesterol, triglyceride, and HDL cholesterol levels), insulin, and highsensitivity C-reactive protein (hsCRP) analysis. e analysis for above biochemical parameters were assessed using Abbott Architect CI8200 Automatic System following the manufacturer's instructions. Homeostatic model assessment to estimate insulin resistance (HOMA-IR) was calculated by multiplying fasting plasma glucose (mmol/L) by fasting serum insulin (µU/ml) and divided by 22.5 [17]. Low-density lipoprotein (LDL) cholesterol was calculated using the Friedewald formula: LDL cholesterol � total cholesterol − ((triglyceride/ 5) + HDL cholesterol) [18]. is information has been reported in our earlier publication [16]. is study investigates ADRB2 rs1042713 gene polymorphism in the same population with the aim to evaluate (i) the effect of ADRB2 rs1042713 on obesity and obesity-related anthropometric and blood biochemical parameters and (ii) the influence of diet on the association between ADRB2 rs1042713 and obesity phenotypes.

Genotyping of ADRB2 rs1042713 Gene Polymorphism.
Five milliliters of whole blood was drawn from an antecubital vein into vacutainer tubes (Becton, Dickinson and Co., Franklin Lakes, NJ) containing EDTA. e genomic DNA was extracted from leukocytes using the MasterPure DNA Purification kit (Lucigen Corporation Middleton, WI, USA) according to the manufacturer's instructions. DNA samples were stored at −20°C until use. A DNA fragment of 310 bp containing rs1042713 was amplified by using polymerase chain reaction (PCR) for the identification of ADRB2 rs1042713 gene polymorphism with specific primers (forward primer: 5′-CCGCCGTGGGTCCGCC-3′ and reverse primer: 5′-CCATGACCAGATCAGCAC-3′) derived from an earlier study [19]. PCR was performed using 5 µl of genomic DNA (∼1 ng/µl), 0.2 µM forward primer, and reverse primer with 5 µl of Taq 5X Master Mix. ermal cycling was performed as follows: initial denaturing at 95°C for 5 min; 35 cycles of denaturation at 94°C for 45 s, annealing at 64°C for 40 s and extension at 72°C for 45 sec, and then a final extension at 72°C for 5 min. e amplicons were verified by using electrophoresis on 2% agarose gel and visualized under ultraviolet illumination after staining by ethidium bromide. e verified amplicons were then sequenced by using BigDye ® terminator v3.1 cycle sequencing kit chemistry.

Power and Sample Size Calculation.
We performed power calculation using software QUANTO, Version 1.2.4, to find the minimum detectable effect for a given sample size.
is calculation takes into account the type 1 error rate of 0.05, and the population prevalence of insulin resistance (using a cut off of HOMA-IR ≥ 1.7) of 45%, as reported in the present study. Given that the minor allele frequency (A) of ADRB2 rs1042713 in our study was 0.49 with 57 insulinresistant and 69 non-insulin-resistant participants, we had 68.2% power to detect an effect of 2.83 (odds ratio) (dominant model).
With regard to the gene-diet interaction, given that the mean of fasting glucose levels in our population was 5.2 mmol/L, with a SD of 2.0, environmental effect (differences of fasting glucose levels between high and low PUFA : SFA ratio) of −0.4 (5.0 mmol/L-5.4 mmol/L), genetic effect (differences of fasting glucose levels between the carriers of G allele of ADRB2 rs1042713 and the noncarriers (A)) of 0.3 (5.4 mmol/L − 5.1 mmol/L), and interaction effect of −1.1 (4.83 mmol/L-5.93 mmol/L), a power of 61% for the genediet interaction was computed.

Statistical Analysis.
Statistical analysis was performed using the statistical package for social sciences (IBM SPSS statistic, Chicago, IL, USA, version 22). Data were expressed as mean ± standard error (SE) or number (percentage). Log transformation was performed to transform nonnormally distributed data into normally distributed data. Independent t-test and chi-squared test were performed to assess the differences between the two genotype groups (AA vs AG + GG) on baseline continuous variables and categorical variables, respectively. Allele frequency was estimated by gene counting and chi-squared test was used to assess deviation from Hardy-Weinberg equilibrium (HWE) [20]. To study the effect of ADRB2 rs1042713 on obesity, data were dichotomised into obese and nonobese groups (obesity was defined by BMI ≥ 27.5 kg/m 2 ) [21]. e AA genotype of ADRB2 rs1042713 was used as the reference group in both codominant and dominant models, whereas the combination of AG and GG genotypes was used as the reference group in recessive model. Odds ratios (ORs) with 95% confidence intervals (95% CIs) were estimated for each genotype by logistic regression to determine the odds of obesity associated with gene variants, after adjusting for covariates age, gender, physical activity status, smoking status, and alcohol consumption. Same analysis was performed to study the odds of ADRB2 gene variants on insulin resistance. Data were dichotomised into (i) non-insulinresistant and insulin-resistant groups (using a cutoff of HOMA-IR ≥ 1.7 [22]).
Differences in means between gene variants in anthropometric, blood biochemical, and dietary parameters in obese and nonobese groups were assessed by using one-way analysis of covariance (ANCOVA). Adjustment for covariates such as age, gender, physical activity status, smoking status, alcohol consumption, BMI, WC, fat mass, body fat percent, and total energy intake were applied where appropriate. e intake of macronutrients (energy-adjusted) was dichotomised into two groups based on the median intake of the population. A multivariate general linear model (GLM) was used to investigate the effect of the interaction between dietary macronutrients and ADRB2 rs1042713 on obesity-related metabolic traits, after adjusting for potential confounders (age, gender, physical activity status, smoking status, alcohol consumption, BMI, and total energy intake). A statistical probability level of p < 0.05 (two-sided) was considered significant. No significant association was found between rs1042713 and protein and carbohydrate intake on obesity-related traits. erefore, the results were not reported.

Baseline Characteristics.
In total, 178 Malaysian adults (female � 154; male � 24) were recruited for anthropometric measurement and genetic analysis. For biochemical and dietary analysis, 126 participants (female � 106, male � 20) were available. General characteristics of the study participants are reported in Table 1. e ages between the two groups (AA vs AG + GG) did not differ significantly (43.4 ± 1.8 y vs 41.1 ± 1.0 y; p � 0.286). e gender distribution was not significantly different between the gene variants in the genotype groups (AA vs AG + GG), females (84.1% vs 87.3%), and males (15.9% vs 12.7%) (p � 0.587). Ethnicity distribution was not significantly different between the two genotype groups (p � 0.556). No significant difference was found in physical activity (p � 0.256), smoking (p � 0.606), and alcohol consumption (p � 0.717) between the two genotype groups. No significant difference was found in height (p � 0.907), body weight (p � 0.505), and BMI (p � 0.427) between the two genotype groups (p > 0.05). Due to the higher number of female participants, all data analyses were adjusted for by gender.

Allele Frequencies of the Gene Variants of ADRB2
rs1042713. G allele was the most frequent variant in our study population (51.1%). e allele frequency for the minor allele of rs1042713 (A allele) was 0.49, which did not deviate from Hardy-Weinberg equilibrium as tested by the chisquared test, χ 2 � 0.36 (χ 2 < 3.841) ( Table 2).

Association between ADRB2 rs1042713 and Obesity and
Insulin Resistance. Logistic regression was performed to examine the independent effect of ADRB2 rs1042713 on the odds of obesity and insulin resistance. We found no significant association between ADRB2 rs1042713 and obesity (obesity as defined by BMI ≥ 27.5 kg/m 2 ) under codominant (AG p � 0.548 and GG, p � 0.884), dominant (p � 0.725), and recessive (p � 0.538) models, after adjusting for covariates age, gender, physical activity status, smoking status, and alcohol consumption (Table 2). However, our results revealed significant association between ADRB2 rs1042713 and insulin resistance (using a cutoff of HOMA-IR ≥ 1.7) (Table 3). e carriers of GG genotype of rs1042713 had increased odds of insulin resistance, compared to AA genotype in both codominant and dominant models, 4.43 (CI � 1.31-15.0, adjusted p � 0.016) and 2.83 (CI � 1.04-7.70, adjusted p � 0.042), respectively, even after adjusting for covariates age, gender, BMI, physical activity status, smoking status, and alcohol consumption. No significant association was found in the recessive model (adjusted p � 0.060).

Differences in Means between ADRB2 rs1042713 Gene Variants in Anthropometric, Biochemical, and Dietary Parameters in Obese and Nonobese Groups.
e age, anthropometric, biochemical, and dietary parameters of the study participants between gene variants in obese and nonobese groups have been reported in Table 4. In obese participants, we found that the carriers of G allele of rs1042713 had significantly higher total cholesterol (p � 0.011), LDL cholesterol levels (p � 0.008), and total cholesterol per HDL cholesterol ratio (p � 0.048), compared to the noncarriers (AA), even after adjusting for covariates age, gender, BMI, WC, fat mass, body fat percent, physical activity status, smoking status, and alcohol consumption. Interestingly, such differences in blood biochemical parameters within genotypes were not observed in the nonobese group. With respect to dietary parameters, we found that the carriers of the G allele had significantly lower consumption of PUFA compared to the noncarriers (AA), even after adjusting for covariates (p � 0.036). No significance association was found between ADRB2 rs1042713 and others dietary parameters.

Discussion
Single-nucleotide polymorphism (SNP) in ADRB2 rs1042713 causes alterations in the structural conformation of the receptor which eventually affect the function of β-adrenergic receptors (ADRB) [27]. is may influence the binding of catecholamines to the beta-2 adrenoceptors and hence alter lipolysis. A meta-analysis involving 18 published articles revealed that there was no association between rs1042713 and obesity [28]. Contrary to the latter, studies on the Saudi population reported significant association between ADRB2 rs1042713 polymorphism and the development of obesity, as also with insulin resistance and dyslipidemia [10,11]. However, findings from Asian populations (Japanese and Korean) reported negative Odds ratios (ORs) with 95% confidence intervals (95% CIs) were estimated for each genotype by logistic regression to determine the risk of insulin resistance associated with ADRB2 rs1042713. p ≤ 0.05 was considered as significant. * Adjusted for age, BMI, gender, physical activity status, smoking status, and alcohol consumption. Non-insulin-resistant and insulin-resistant groups were dichotomised using a cutoff of HOMA-IR ≥ 1.7 [22]. e agreement of genotype frequencies with Hardy-Weinberg equilibrium was tested by using the chi-squared test, with χ 2 < 3.841 considered as no deviation from Hardy-Weinberg equilibrium. Logistic regression was conducted to determine the risk of obesity associated with gene variants. Odds ratios (ORs) with 95% confidence intervals (95% CIs) were estimated for each genotype. * Adjusted for age, gender, physical activity status, smoking status, and alcohol consumption. p < 0.05 was considered significant. One-way analysis of covariance was performed to determine the differences between means between gene variants in anthropometric, biochemical, and dietary parameters in obese and nonobese participants, after adjusting for covariates in different model 1 age, gender, physical activity status, smoking status, and alcohol consumption; 2 model 1 + BMI, WC, body fat mass, and body fat percent; 3 model 1 + BMI; and 4 model 3 + total energy intake. * p < 0.05 was considered significant. HOMA-IR: homeostatic model assessment-insulin resistance; HDL: high-density lipoprotein; LDL: low-density lipoprotein; hsCRP: high-sensitivity C-reactive protein; SFA: saturated fatty acid; MUFA: monounsaturated fatty acid; PUFA: polyunsaturated fatty acid; % TE: percentage of total energy intake.  [22] 0.003 * Total cholesterol (mmol/L) 5.7 ± 0.1 5.2 [26] 0.001 * Triglyceride (mmol/L) 1.4 ± 0.1 1.7 [24] <0.001 * HDL cholesterol (mmol/L) 1.5 ± 0.1 1.0 [24] <0.001 * LDL cholesterol (mmol/L) 3.6 ± 0.1 2.6 [26] <0.001 * Total cholesterol/HDL cholesterol 3.9 ± 0.1 -e one-sample t-test was performed to assess the difference between blood biochemical parameters and respective diagnostic cutoffs * p < 0.05 was considered significant. HOMA-IR: homeostatic model assessment-insulin resistance; HDL: high-density lipoprotein; LDL, low-density lipoprotein. 6 Journal of Nutrition and Metabolism association between obesity and ADRB2 gene polymorphisms [12,13]. Moreover, these studies did not investigate interaction of gene variants, phenotypes, and dietary nutrients. To date, there are no data reported on ADRB2 rs1042713 in the Malaysian population. To the best of our knowledge, this is the first study that investigates the effect of ADRB2 rs1042713 on obesity and obesity-related metabolic parameters and its interaction with dietary nutrients in Malaysian adults.

Association between ADRB2 rs1042713 and Obesity and
Insulin Resistance. In our study, we found no association between ADRB2 rs1042713 and odds of obesity (obesity was defined as BMI ≥ 27.5 kg/m 2 ). However, ADRB2 rs1042713 was associated with insulin resistance (using a cutoff of HOMA-IR ≥ 1.7). We found that the carriers of G allele of rs1042713 had increased odds of insulin resistance compared to the noncarriers (AA), in the dominant model, even after adjusting for potential confounders. ese findings  Figure 2: e effect of the interaction between ADRB2 rs1042713 and PUFA intake on (a) HOMA-IR and (b) fasting insulin levels. PUFA intake was dichotomised into two groups by using the median value of 6% of TE/day (for all participants) for analysis. Gene-diet interaction was evaluated by using the multivariate general linear model after adjusting for age, gender, BMI, physical activity status, smoking status, alcohol consumption, and total energy intake. * p < 0.05 was considered statistically significant. PUFA: polyunsaturated fatty acid; HOMA-IR: homeostatic model assessment-insulin resistance.
suggest that variations in ADRB2 rs1042713 may interfere with glucose homeostasis and cause insulin resistance. Prior et al. reported that ADRB2 Arg16Gly-Gln27Glu haplotype was associated with glucose intolerance and insulin resistance in obese postmenopausal women [29]. A possible explanation for this observation could be the alteration in the structural conformation of the receptor, which may have enhanced sympathetic stimulation leading to increased lipolysis [8]. is overstimulation of ADRB2 is found to be associated with the pathogenesis of insulin resistance as it inhibits the insulin-induced translocation of GLUT4 and reduces glucose uptake via the cAMP-dependent protein kinase A-dependent pathways [30]. We found that HOMA-IR in the obese individuals carrying G allele of ADRB2 rs1042713 was above the diagnostic cutoff (1.7). It is now an established fact that higher levels of nonesterified fatty acids (NEFA) in the blood can induce preferential use of free fatty acids over glucose to generate ATP even in the presence of insulin in muscle and adipose tissue resulting in hyperglycaemia [31,32]. Stimulation by noradrenaline of adipose tissue ADRB2 increases the release of NEFA. In addition, free fatty acids (FFAs) can stimulate hepatic gluconeogenesis and alter pancreatic insulin release and subsequent metabolism in individuals with impaired glucose metabolism [33].

Association between ADRB2 rs1042713 and Blood Lipid
Levels. In the present study, we report that obese individuals carrying the G allele of rs1042713 had significantly higher total cholesterol, LDL cholesterol levels, and total cholesterol per HDL cholesterol ratio compared to the noncarriers (AA). Total cholesterol (5.7 ± 0.1 mmol/L) and LDL cholesterol levels (3.6 ± 0.1 mmol/L) in obese individuals carrying the G allele were above the diagnostic cutoffs (5.2 mmol/L and 2.6 mmol/L, respectively [26]). However, these differences in biochemical parameters were not observed in nonobese individuals. In our participants, excess FFAs in circulation in obese individuals carrying the gene variants of ADRB2 may have driven dyslipidemia. Increased lipolysis due to the polymorphisms of ADRB2 gene may have caused increased levels of NEFA induced hepatic production of VLDL and hence higher LDL levels in our participants [34].

Interaction between Dietary Fats and ADRB2 rs1042713
on Glycaemic Indices. We report that the level of fasting glucose was modulated by the types of dietary fatty acids such as SFA and PUFA in the carriers of GG genotype of ADRB2 rs1042713. e carriers of GG genotype of rs1042713 had significantly lower fasting glucose levels with intake of relatively higher PUFA:SFA ratio (≥0.8/day) (Figure 1(b)) and lower SFA intake (<7.3% of TE/day) (Figure 1(a)). Moreover, the carriers of GG genotype of rs1042713 had significantly lower fasting insulin levels and HOMA-IR when consuming higher PUFA intake (≥6% of TE/day). ese findings suggest that carriers of GG genotype of ADRB2 rs1042713 consuming higher percentage of PUFA demonstrated better homeostatic control of fasting blood glucose and insulin sensitivity. A meta-analysis of 102 randomised controlled feeding trials with 4,200 subjects has reported that PUFA had the most beneficial effects in improving glycaemia, insulin resistance, and insulin secretion in comparison to dietary carbohydrate, SFA, and MUFA [35]. e anti-inflammatory properties of PUFA increase the production of adiponectin via PPARα activation and alleviate adipose tissue inflammation via GPR120 and resolvins/protectins, which favour insulin sensitivity. It also suppresses oxidative stress and pancreatic lipotoxicity, reduces toxicity of tissue free fatty acids, and increases membrane fluidity [36].
is body of evidence indicates that the composition of dietary fatty acid intake plays an important role in affecting glucose metabolism and insulin sensitivity.
As per Malaysian recommendations, 10% of total energy should come from SFA and 3 to 8% of total energy should come from PUFA [23]. On an average, dietary intake of our participants with respect to SFA and PUFA were within these ranges. However, individuals with G allele of rs1042713 may need to consume higher quantity of PUFA to combat diet-related noncommunicable diseases. Replacing isocaloric quantity of foods rich in SFA with PUFA may improve homeostatic control of blood glucose and enhance insulin sensitivity (HOMA-IR) in such individuals.

Limitations
A major limitation of our study is the small sample size. We acknowledge that the current study is exploratory and is underpowered to detect the gene-diet interaction between ADRB2 rs1042713 and dietary nutrients on phenotypic and metabolic alterations in our population. However, with association and interaction analysis in the current study, we have generated a hypothesis. In future, large-scale studies are required to confirm such findings in the Malaysian population. In this study, we did not stratify our participants by gender for analysis due to the small sample size of male participants. To account for this, we have adjusted for gender statistically in all our data analysis to eliminate Type 1 error.

Conclusion
In conclusion, our study revealed that there was no association between ADRB2 rs1042713 and obesity. However, ADRB2 rs1042713 was associated with insulin resistance in Malaysian adults. e carriers of G allele of rs1042713 had increased odds of insulin resistance compared to noncarriers (AA). Obese individuals carrying the G allele of rs1042713 had significantly higher total cholesterol, LDL cholesterol levels, and total cholesterol per HDL cholesterol ratio compared to the noncarriers (AA). ese differences were not observed in nonobese individuals. ere is evidence from earlier studies that high PUFA intake is associated with favourable effects on glycaemia and insulin resistance. Over and above the latter, we found that higher PUFA intake was beneficial in individuals carrying the G allele with respect to glycaemic indices compared to the noncarriers. Although it is premature to report gene-diet interaction on the