Interaction between TCF7L2 polymorphism and dietary fat intake on high density lipoprotein cholesterol

Abstract

Recent evidence suggests that lifestyle factors influence the association between the Melanocortin 4 receptor (MC4R) and Transcription Factor 7-Like 2 (TCF7L2) gene variants and cardio-metabolic traits in several populations; however, the available research is limited among the Asian Indian population. Hence, the present study examined whether the association between the MC4R single nucleotide polymorphism (SNP) (rs17782313) and two SNPs of the TCF7L2 gene (rs12255372 and rs7903146) and cardio-metabolic traits is modified by dietary factors and physical activity. This cross sectional study included a random sample of normal glucose tolerant (NGT) (n = 821) and participants with type 2 diabetes (T2D) (n = 861) recruited from the urban part of the Chennai Urban Rural Epidemiology Study (CURES). A validated food frequency questionnaire (FFQ) was used for dietary assessment and self-reported physical activity measures were collected. The threshold for significance was set at P = 0.00023 based on Bonferroni correction for multiple testing [(0.05/210 (3 SNPs x 14 outcomes x 5 lifestyle factors)]. After Bonferroni correction, there was a significant interaction between the TCF7L2 rs12255372 SNP and fat intake (g/day) (Pinteraction = 0.0001) on high-density lipoprotein cholesterol (HDL-C), where the ‘T’ allele carriers in the lowest tertile of total fat intake had higher HDL-C (P = 0.008) and those in the highest tertile (P = 0.017) had lower HDL-C compared to the GG homozygotes. In a secondary analysis of SNPs with the subtypes of fat, there was also a significant interaction between the SNP rs12255372 and polyunsaturated fatty acids (PUFA, g/day) (Pinteraction<0.0001) on HDL-C, where the minor allele carriers had higher HDL-C in the lowest PUFA tertile (P = 0.024) and those in the highest PUFA tertile had lower HDL-C (P = 0.028) than GG homozygotes. In addition, a significant interaction was also seen between TCF7L2 SNP rs12255372 and fibre intake (g/day) on HDL-C (Pinteraction<0.0001). None of the other interactions between the SNPs and lifestyle factors were statistically significant after correction for multiple testing. Our findings indicate that the association between TCF7L2 SNP rs12255372 and HDL-C may be modified by dietary fat intake in this Asian Indian population.

About the Authors:

Dhanasekaran Bodhini

Roles Formal analysis, Investigation, Writing – review & editing

Affiliation: Department of Molecular Genetics, Madras Diabetes Research Foundation, Chennai, India

Szilvia Gaal

Roles Writing – original draft

Affiliation: Hugh Sinclair Unit of Human Nutrition and Institute for Cardiovascular and Metabolic Research (ICMR), Department of Food and Nutritional Sciences, University of Reading, Reading, United Kingdom

Israa Shatwan

Roles Formal analysis

Affiliation: Hugh Sinclair Unit of Human Nutrition and Institute for Cardiovascular and Metabolic Research (ICMR), Department of Food and Nutritional Sciences, University of Reading, Reading, United Kingdom

Kandaswamy Ramya

Roles Investigation

Affiliation: Department of Molecular Genetics, Madras Diabetes Research Foundation, Chennai, India

Basma Ellahi

Roles Writing – review & editing

Affiliation: Faculty of Health and Social Care, University of Chester, Chester, United Kingdom

Shelini Surendran

Roles Software

Affiliation: Hugh Sinclair Unit of Human Nutrition and Institute for Cardiovascular and Metabolic Research (ICMR), Department of Food and Nutritional Sciences, University of Reading, Reading, United Kingdom

Vasudevan Sudha

Roles Investigation, Supervision, Writing – review & editing

Affiliation: Department of Foods, Nutrition and Dietetics Research, Madras Diabetes Research Foundation, Chennai, India

Mohan R. Anjana

Roles Project administration, Writing – review & editing

Affiliations Department of Diabetology, Madras Diabetes Research Foundation, Chennai, India, Dr. Mohan’s Diabetes Specialties Centre, WHO Collaborating Centre for Non-communicable Diseases Prevention and Control, Chennai, India

Viswanathan Mohan

Roles Conceptualization, Supervision, Writing – review & editing

Affiliations Department of Diabetology, Madras Diabetes Research Foundation, Chennai, India, Dr. Mohan’s Diabetes Specialties Centre, WHO Collaborating Centre for Non-communicable Diseases Prevention and Control, Chennai, India

Julie A. Lovegrove

Roles Supervision, Writing – review & editing

Affiliation: Hugh Sinclair Unit of Human Nutrition and Institute for Cardiovascular and Metabolic Research (ICMR), Department of Food and Nutritional Sciences, University of Reading, Reading, United Kingdom

Venkatesan Radha

Roles Methodology, Supervision

Affiliation: Department of Molecular Genetics, Madras Diabetes Research Foundation, Chennai, India

Karani Santhanakrishnan Vimaleswaran

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Hugh Sinclair Unit of Human Nutrition and Institute for Cardiovascular and Metabolic Research (ICMR), Department of Food and Nutritional Sciences, University of Reading, Reading, United Kingdom

ORCID http://orcid.org/0000-0002-8485-8930

Introduction

Genetic variants, unhealthy dietary intake, physical inactivity and their multiple interactions are considered to be contributory factors to the development of obesity and type 2 diabetes (T2D) [1–3]. After China, India has the highest number of people with T2D in the world and according to the Indian Council of Medical Research–INdia DIABetes (ICMR–INDIAB) study, T2D cases have reached 62.4 million and 77.2 million people are pre-diabetic [4]. Furthermore, obesity and T2D are risk factors for non-communicable diseases (NCDs) such as cardiovascular disease (CVD) and it is estimated that India will have the highest rate of CVD mortality in the world [5]. Asian Indians have different biochemical characteristics from other populations from birth, often referred to as the ‘Asian Indian Phenotype’ which consists of increased visceral fat and waist circumference, hyperinsulinemia and insulin resistance [6].

The association between several genes and metabolic diseases has been identified by the candidate gene approach and genome-wide scans; to date, besides the FTO (Fat mass and obesity associated) gene, the strongest obesity risk loci known so far [1, 3, 7, 8], two commonly studied candidates for obesity and T2D have been the Melanocortin 4 Receptor (MC4R) and Transcription Factor 7-Like 2 (TCF7L2) genes. Strong association between the MC4R gene and risk of obesity was identified by a genome-wide association (GWAS) study [9] whereas the association between the Transcription Factor 7-Like 2 (TCF7L2) gene and risk of T2D was identified by a genome wide linkage study [10]. MC4R is expressed in the hypothalamus within the brain and therefore it is suggested that it contributes to body weight regulation by its effect on food intake and energy homeostasis [11]. A strong association was identified between MC4R rs17782313 genetic variant and risk of obesity in a European population [12] which was then replicated in other populations [13–15], including Asian Indians [16]. The TCF7L2 gene is involved in the Wnt signalling pathway where it affects the expression of pro-glucagon and consequently blood glucose regulation [10]. In addition, the effect of pro-glucagon on Glucagon-like peptide 1 (GLP-1) also influences blood glucose regulation with insulin [17]. Decreased insulin secretion and increased glucose production in the liver is suggested to be the result of over expression of the TCF7L2 gene [18]. A couple of studies [19, 20] have shown strong associations between the two TCF7L2 single nucleotide polymorphisms (SNPs) (rs7903146, rs12255372) and risk of T2D among Asian Indians living in India, in addition to a meta-analysis [17].

In recent years, several studies have examined whether the association between the genetic variants of MC4R and TCF7L2 genes and cardio-metabolic traits is modified by lifestyle factors such as diet and physical activity in various populations [21–25], however there are no studies to date among Asian Indians [26]. Whilst most studies in European populations found no significant interactions between the MC4R SNP rs17782313 and dietary factors on obesity traits [11, 22, 23], a prospective cohort study reported significant interactions between the SNP and fat and protein intake on body mass index (BMI) and risk of T2D [27]. Similarly some studies identified significant interactions between TCF7L2 SNPs and fibre and fat intake on T2D[28–30]; however there were discrepancies between the studies which could be due to differences in sample size, study design, dietary assessment and genetic heterogeneity. The objectives of this paper were to determine whether the MC4R SNP rs17782313 and TCF7L2 SNPs (rs7903146, rs12255372) were associated with cardio-metabolic traits and whether the association was modified by diet and physical activity in a sample size of up to 1,682 adults from the Chennai Urban Rural Epidemiology Study (CURES).

Materials and methods

Study participants

A random, unrelated sample of normal glucose tolerant (NGT) (n = 821) and T2D (n = 861) participants were recruited from the urban part of the cross sectional Chennai Urban Rural Epidemiology Study (CURES) which included 26,001 individuals in total as a representative sample of Chennai city. Full details of the methodology have been explained previously [31]. Briefly 26,001 adult subjects (>20 years of age) were recruited in Phase 1 of CURES using a systematic random sampling method covering the whole Chennai city. This included 1,529 ‘self-reported’ diabetic subjects. In Phase 2, all self-reported or ‘known diabetic’ subjects were invited to our centre for detailed studies of whom 1382 responded (response rate: 90.4%). In Phase 3, every tenth subject from Phase 1, excluding those with self-reported diabetes, underwent an oral glucose tolerance test (OGTT). Subjects who had fasting plasma glucose <5.6 mmol/l (100mg/dl) and 2 hr plasma glucose value 7.8 mmol/l (140mg/dl) were categorized as having NGT [32]. Those who were confirmed by OGTT to have 2 hr plasma glucose value 11.1 mmol/l (200 mg/dl) were classed as ‘newly detected diabetic subjects’ (n = 222).For the present study, the NGT subjects were selected from Phase 3 and subjects with T2D which included ‘known diabetic’ and ‘newly detected diabetic’ subjects were selected from Phase 2 and Phase 3 of the Chennai Urban Rural Epidemiology Study. The Madras Diabetes Research Foundation Institutional Ethics Committee granted the ethical approval and informed consent was obtained from the study participants.

Anthropometric and biochemical measurements

Weight, height and waist circumference were measured by standard methods. The BMI calculation was based on the body weight (kg) divided by the square of body height (m). Roche Diagnostics (Mannheim) provided the equipments in order to be able to carry out the biochemical analyses on a Hitachi-912 Auto Analyzer (Hitachi, Mannheim, Germany)(8). Fasting plasma glucose, serum cholesterol, serum triglycerides (TG) and high-density lipoprotein cholesterol were measured by glucose oxidase-peroxidase, cholesterol oxidase-phenol-4-amino-antipyrene peroxidase, glycerol phosphatase oxidase-phenol-4-amino-antipyrene peroxidise and polyethylene glycol-pretreated enzyme methods respectively(8). The Friedewald formula was used to estimate low-density lipoprotein cholesterol concentrations [33]. Glycated haemoglobin (HbA1c) and serum insulin were determined by high-performance liquid chromatography (HPLC) on a Variant instrument (Bio-Rad, Hercules, CA, USA) and an enzyme-linked immunosorbent assay (Dako, Glostrup, Denmark) respectively [8].

Assessment of dietary intake and physical activity

A validated, interviewer administered semi-quantitative food frequency questionnaire (FFQ) [34] consisting of a list of 222 different foods was used in order to evaluate dietary intake for the previous year including macronutrient and total energy intake. Frequencies (per day, week, month, year, never) and portion sizes were estimated by the participants with the help of visual aids of measurement equipments and food sizes. Daily average food and nutrient intake was calculated by the EpiNu database system.

A validated self-report questionnaire was used to measure physical activity [5]. Individuals were divided into the vigorously active group when they both exercised and engaged in demanding work activities whereas within the moderately active group the participants either exercised or carried out heavy physical work. The remainder of the study participants were separated into the sedentary group.

SNP selection and genotyping

A strong association was identified between MC4R rs17782313 genetic variant and risk of obesity in a European population [12] which was then replicated in other populations (13–15), including Asian Indians [16]. Strong associations have also been found between TCF7L2 SNPs (rs12255372 and rs7903146) and risk of T2D in a meta-analysis [17] in the Japanese population [35] as well as in Asian Indians [19, 20]. Although the TCF7L2 SNPs (rs12255372 and rs7903146) have showed significant linkage disequilibrium (r2 = 0.746), the rs7903146 variant has been shown to have the strongest effect in Caucasian populations [36, 37]. Based on the previous studies, the above mentioned three SNPs were selected for the present study.

Phenol-chloroform method was used to extract DNA from whole blood. The methodology for genotyping TCF7L2 rs12255372 (G/T) and rs7903146 (C/T) SNPs has been previously published [20]. Direct sequencing by an ABI 310 genetic analyzer (Applied Biosystems, Foster City, CA) helped to confirm the efficiency of the genotyping which was in 99% concordance based on random duplicates of 20% of the samples [20].

MC4R rs17782313 (T/C) SNP: PCR volume (10 μl) consisted of 1X reaction buffer, 200 mmol of dNTP, 1.5 mmol of MgCl2, 1 U taq DNA polymerase and 100 ng of genomic DNA. The concentrations for primers included 15 pmol of common primer, 15 pmol of allele 1 primer and 1 pmol of allele 2 primer, equal to a 15:1 ratio of short primers to long primer. Cycles for PCR were carried out at 96°C for 12 min, then 35 cycles at 94°C for 30 sec, followed by 30 sec at 57°C, 30 sec at 72°C and finally for 10 min at 72°C. Electrophoresis was carried out with a 3% agarose gel.

Inner primers

F: AAGTTTAAAGCAGGAGAGATTGTATACC (C allele 222bp)

R: GCTTTTCTTGTCATTTCCAGCA (T allele 149bp)

Outer primer

F: TTACTGATTTTAAGGGCATAAGCAA

R: TATCATGCTGAGACAGGTTCATAAA (321bp)

Statistical analyses

Statistical analyses were carried out by using SPSS software (version 21). BMI ≥ 25 kg/m2 was categorised as obese and BMI < 25 kg/m2 as non-obese. Descriptive statistics for continuous variables are shown as means and standard deviation (SD). Fasting serum insulin and triglyceride values were log transformed to obtain normal distribution. Genotype frequencies between cases and controls were compared by Chi Square test. The difference in the means of continuous variables between the genotypes was analysed by independent sample t test. Based on an additive model of analysis, dominant models were used for all 3 single nucleotide polymorphisms (SNPs) where the common homozygous allele was compared to the combined heterozygous and rare homozygous alleles due to the low allele frequency of rare homozygotes. Association analyses between SNPs and continuous and categorical variables were carried out by linear and logistic regression models, respectively, adjusting for age, gender with the addition of BMI, when T2D was the outcome, and, adjusting for T2D, when obesity was the outcome. Linear and logistic regression models were also used for interaction analyses between SNPs and dietary factors (continuous variables) / physical activity (categorical variable) on continuous and categorical outcomes respectively, where the interaction terms were included into the models and were adjusted for age, gender, BMI, T2D and total energy intake when appropriate. Bonferroni correction for multiple testing was calculated by multiplying 3 SNPs with 14 outcomes (T2D, obesity, BMI, waist circumference, fasting blood glucose, HbA1c, fasting insulin, systolic and diastolic blood pressure, HDL, LDL, VLDL, TG and total cholesterol) and 4 dietary factors (carbohydrate, protein, fat, fibre) and physical activity level. The P value of 0.05 was then divided by 210 (3 SNPs x 14 outcomes x 5 lifestyle factors) which set the significant p value for all results at P = 0.00023.

As a secondary analysis, given the significant SNP-fat intake interaction, individuals were grouped into tertiles based on the fatty acid subtypes [monounsaturated (MUFA) and polyunsaturated (PUFA) fatty acids] for testing the interaction between SNPs and these fatty acid subtypes on lipids. In addition, PUFA was further stratified into ALA (alpha linolenic acid) and LA (linoleic acid) for the interaction analysis.

Given that there are no previously reported effect sizes for the gene-diet interaction, we were unable to perform a prospective power calculation. However, based on the most significant interaction observed in the present study, we performed retrospective power calculations using QUANTO software, Version 1.2.4 (May 2009). We performed power calculations in the form of least detectable effects based on the assumption of significance levels and powers of 5 and 80%, respectively. At 80% power, the minimum detectable effects ranged from beta 0.02 mg/dl (HDL-C) for a SNP with MAF of 5% to beta 1.0 mg/dl for a SNP with MAF 50% in the case-control analysis. For the TCF7L2 SNP–fat intake interaction on HDL-C (most significant interaction), the beta was 0.067, which is within the range of effect sizes for which the power was calculated.

Results

Phenotypic associations

Based on the clinical and biochemical characteristics of the individuals from the CURES study as illustrated in Table 1, individuals with T2D were older (P<0.0001), had higher BMI (P<0.0001), waist circumference (P<0.0001), fasting plasma glucose (P<0.0001), HbA1c (P<0.0001), fasting plasma insulin (P<0.0001), systolic and diastolic blood pressure (P<0.0001), low density and very low density lipoproteins (P<0.0001), total cholesterol (P<0.0001), TG (P<0.0001) and lower high density lipoprotein (P = 0.001) than NGT individuals.

[Figure omitted. See PDF.]

Table 1. Anthropometric and biochemical characteristics of T2D and NGT participants.

https://doi.org/10.1371/journal.pone.0188382.t001

Genetic associations

The minor allele (‘T’) of SNPs (rs12255372 and rs7903146) of TCF7L2 gene showed significant susceptibility to T2D; however, in the present study after correction for multiple testing only the association between SNP rs7903146 and T2D remained statistically significant (P = 0.0001) (S1 Table).After Bonferroni correction, statistically significant association between the ‘C’ allele of theMC4R SNP rs17782313 and T2D was also observed (P = 0.00022) (S2 Table). None of the associations between the three SNPs and continuous variables remained significant after correction for multiple testing (P>0.00023) (S1 and S2 Tables).

TCF7L2 –dietary fat intake interactions on HDL-C

Individuals were grouped into tertiles based on their fat and subtypes of fat intake (g/day). The means for fat intake in the 1st tertile: 41 g/d; 2nd tertile: 62 g/d; 3rd tertile: 95 g/d. Means for PUFA intake in the 1st tertile: 9g/d; 2nd tertile: 17 g/d; 3rd tertile: 29 g/d. Means for MUFA intake in the 1st tertile: 12 g/d; 2nd tertile: 18 g/d; 3rd tertile: 29 g/d. Means for ALA intake in the 1st tertile: 0.38 g/d; 2nd tertile: 0.58 g/d; 3rd tertile: 0.89 g/d. Means for LA intake in the 1st tertile: 8 g/d; 2nd tertile: 17 g/d; 3rd tertile: 29 g/d.

The interaction between the TCF7L2 SNP rs12255372 and fat intake (g/day) on HDL-C was statistically significant after correction for multiple testing (Pinteraction = 0.0001). The ‘T’ allele carriers of the TCF7L2 SNP rs12255372 had 2.26 mg/dl higher HDL-C level in the lowest tertile of fat intake (mean: 41 g/day) than the ‘GG’ homozygotes (P = 0.008) and in the highest tertile of fat intake (mean: 95 g/day), HDL-C was 1.87 mg/dl lower in the risk ‘T’ allele carriers in comparison to the ‘GG’ homozygotes (p = 0.017) (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Interaction of the TCF7L2 gene polymorphism (rs12255372) with fat (g) intake, PUFA intake and Alpha Linolenic Acid (g) intake on HDL-C.

Individuals carrying the ‘XT’ genotype had 2.26 mg/dl higher HDL-C in the lowest fat tertile (P = 0.008), while those in the highest tertile had 1.87 mg/dl lower HDL-C (P = 0.017) than those who carry the ‘GG’ allele. Carriers of the ‘XT’ genotype had 1.96 mg/dl higher HDL-C in the 1st tertile of PUFA intake (g) (P = 0.024), while those in the 3rd tertile had 1.64 mg/dl lower HDL-C in comparison to the carriers of the ‘GG’ genotype (P = 0.028). In the 1st tertile of Alpha Linolenic acid intake (g), individuals with the ‘XT’ genotype had 2.42 mg/dl higher HDL-C than the ‘GG’ homozygotes (P = 0.004).

https://doi.org/10.1371/journal.pone.0188382.g001

Stratification to fat subgroups showed significant interactions between the SNP rs12255372 and PUFA (g/day) on HDL-C (Pinteraction<0.0001), where the ‘T’ allele carriers had 1.96 mg/dl higher HDL-C (P = 0.024) in the low PUFA tertile (mean: 9 g/day) in comparison to the ‘GG’ homozygotes and in the 3rd tertile (mean: 29 g/day), the HDL-C level of the ‘T’ allele carriers was 1.64 mg/dl lower than the ‘GG’ homozygotes (P = 0.028) (Fig 1). A similar interaction was also found between the SNP rs12255372 and MUFA (g/day) on HDL-C (Pinteraction = 0.0003), where the ‘T’ allele carriers had 1.77 (mg/dl) higher HDL-C in the lowest MUFA tertile (mean: 12 g/day) (P = 0.03) and had 1.61 (mg/dl) higher HDL-C in the 2nd tertile (mean: 18 g/day) (P = 0.045) than the ‘GG’ carriers, however in the highest MUFA tertile (mean: 29 g/day) the ‘T’ allele carriers had 1.59 (mg/dl) decreased HDL-C (P = 0.041) than individuals with the ‘GG’ genotype.

PUFA was further stratified to linoleic acid (LA) and alpha linolenic acid (ALA) to investigate whether omega-3 and omega-6 fatty acids modified the association between the TCF7L2 SNP rs12255372 and HDL-C. Significant interaction was found between the SNP and ALA on HDL-C (Pinteraction = 0.012), where the ‘T’ allele carriers had 2.42 (mg/dl) higher HDL-C than the ‘GG’ homozygotes (P = 0.004) in the lowest tertile (mean: 0.38 g/day) (Fig 1).A similar interaction was also found between the SNP rs12255372 and LA (g/day) on HDL-C (Pinteraction<0.0001). S1 Table shows the interactions of TCF7L2 SNP rs12255372 with fat, PUFA and ALA intakes on HDL-C under an additive and dominant model.

Additional gene-diet interactions

In addition to the main significant findings, there was a significant interaction between TCF7L2 SNP rs12255372 and fibre intake (g/day) on HDL-C (Pinteraction< 0.0001), where in the lowest tertile (mean: 20 g/day) individuals carrying the ‘T’ allele had 1.95 (mg/dl) higher HDL-C (P = 0.02), in the 2nd tertile (mean: 30 g/day), the ‘T’ allele carriers had 2.39 (mg/dl) higher HDL-C (P = 0.003) and in the highest tertile (mean: 44 g/day), the HDL-C level of the ‘T’ allele carriers was 2.37 mg/dl lower in comparison to the ‘GG’ homozygotes (P = 0.002). There were several other interactions which did not reach statistical significance after correction for multiple testing and these interactions are shown in the Table 2 and S2 and S3 Tables).

[Figure omitted. See PDF.]

Table 2. Interactions of TCF7L2 SNPs rs12255372 and rs7903146 with carbohydrate, fibre and protein intake on HDL-C, fasting blood glucose and diastolic blood pressure.

https://doi.org/10.1371/journal.pone.0188382.t002

Gene-physical activity interactions on cardio-metabolic traits

No statistically significant interactions were observed after correction for multiple testing between the three SNPs and physical activity on obesity- and T2D- related traits (Pinteraction>0.0002). Interactions observed between TCF7L2 SNP rs7903146 and physical activity on VLDL (Pinteraction = 0.012) and TG (Pinteraction = 0.014), where the risk ‘T’ allele carriers had higher TG and VLDL levels, did not reach statistical significance after correction for multiple testing (S4 Table).

Discussion

This is the first study to investigate interactions between TCF7L2 and MC4R SNPs and lifestyle factors on cardio-metabolic traits among Asian Indians. The main findings suggest that total fat and PUFA intakes may modify the association between the TCF7L2 SNP rs12255372 and HDL-C. In the lowest total fat and PUFA tertiles, participants carrying the ‘XT’ genotype (GT+TT) had significantly higher HDL-C, whereas in the highest tertile participants carrying the ‘XT’ genotype had lower HDL-C, as compared to those with the GG genotype. This finding is of public health significance given that Asian Indians tend to have low HDL-C which puts them at markedly increased risk for CVD [38, 39].

The fat intake consists of invisible, hidden and visible fat (vegetable oils, ghee, butter) in India where in the urban areas the minimum average daily intake of visible fat is approximately 22–45 g/day [40]. The recommended omega 3 PUFA intake is 0.1 gram/day, however the median of combined intake of LA and ALA is 13 g/day and the recommended ratio between LA and ALA is 5:1 to 10:1 [40]. In comparison, the results of our study suggest that total fat intake below 53 g/day, PUFA intake below 13 g/day, and ALA intake below 0.49 g/day may help maintain high levels of HDL-C in the ‘T’ allele carriers of the TCF7L2 SNP rs12255372. The findings that MUFA intake above 22 g/day reduces HDL-C and only ALA intake below 0.5 g/day maintains high HDL-C level in risk carriers are unexpected results; however, the effect of polygenic traits cannot be ruled out.

HDL-C is generally considered to be protective against CVD due to its role in the reverse cholesterol transport; however, the recent Mendelian randomization (MR) studies [41, 42] have failed to show a causal effect of low HDL-C on cardiovascular disease risk. But, none of these MR studies have taken dietary factors into account. Furthermore, fatty acids have shown to have different modulating effect on HDL-C for which the mechanism is not fully understood [43]. On the other hand, it is argued that the total fat intake is more influential on postprandial lipoprotein abnormalities, which is a characteristic of T2D, than the type of fatty acids [44]. However, studies have also shown that different fatty acids can have differential effects on postprandial TG [45–47]. In our study, we found that those in the lowest fat intake subgroups had higher HDL-C while those in the highest tertile groups had lower HDL-C among the risk allele carriers. A previous intervention study in the US population [30] also found significant interaction between SNP rs12255372 and fat intake where positive changes in body composition were observed in the ‘T’ risk allele carriers only on the low-fat diet (20% from total energy). Other studies only reported significant interactions between SNP rs7903146 and high saturated fat intake on increased metabolic syndrome risk [29] and high (n-6) PUFA intake (≥6.62% of energy intake) on increased VLDL and TG [24] among the risk ‘T’ allele carriers.

High omega 6 to omega 3 ratio has been shown to increase HDL in mice [43], whereas the results from this study in humans contradict those findings. It has been shown in cell and animal studies that high omega 6 PUFA intake is pro-inflammatory leading to an increased risk for CVD and diabetes [48] and in the present study reduced HDL was observed in the minor allele carriers in response to high intake of omega 6 PUFA (>21.6 g/d). According to our findings, a low omega 6 (<12.8 g/d) and a low ALA intake (<0.5 g/d) may help maintain HDL above 44 (mg/dl) in the risk allele carriers.

Though there was no significant difference in mean HDL levels in the different tertiles of total fat and PUFA intakes in the study participants (data not shown), in the presence of TCF7L2 rs12255372 genotype, a clear interaction between the genotype and fat/ PUFA intake on HDL-C was observed. While the carriers of ‘XT’ genotype had increased HDL-C in the presence of a low fat/ low PUFA diet, there was a decrease in HDL-C levels in the carriers of ‘XT’ genotype in the presence of a high fat / high PUFA diet. It is to be noted that among those who carry the GG genotype, irrespective of the dietary fat/PUFA intake, there is no effect on HDL-C levels.

Another interesting interaction was the one between TCF7L2 SNP rs12255372 and fibre intake on HDL-C. The results suggest that low and medium fibre intake (means: 20, 30 g/day), respectively, increase HDL-C whereas high fibre intake (mean: 44 g/day) may reduce HDL-C among risk allele carriers of the TCF7L2 SNP rs12255372. The average fibre intake in India is 30–40 g/day [40] which is consistent with the mean intake in the medium and high tertiles and also higher compared to the mean fibre intake in the UK (~18g/day) and the US (~16g/day). A previous study investigated interaction between SNP rs7903146 and fibre intake on T2D risk [28],where the minor allele carriers had increased T2D risk with higher fibre intake (mean intake: 13.1 ± 2.2 g/4,184 kJ). High fibre intake is generally recommended to improve glycemic control in T2D individuals [49] and to reduce total and LDL-C in order to reduce CVD risk. A meta-analysis of data from clinical studies (n = 2,990) [50] indicated that high fibre diets (20–30 g/day) reduced HDL-C. Similarly, another meta-analysis of data from 24 clinical studies also suggested that medium and high carbohydrate, high fibre (≥20 g/day) diets also decreased HDL by 4% [51].However, it was argued, that the decrease in LDL-C and TG values would reduce CVD risk by 16.4% which would outweigh the increased risk of CVD by 11.9% due to decreased HDL [51]. But it is of note that these meta-analyses did not consider the genetic component and, furthermore, we also observed gene-diet interaction on other lipid outcomes such as VLDL and TG, which did not remain significant after correction for multiple testing.

The risk of coronary artery disease is increased in individuals with circulating HDL-C concentration of <40 mg/dl [52]. In our study, within the high fibre intake tertile (≥35g/day), the mean HDL-C level of the minor risk allele carriers was 39.46 mg/dL in comparison to 42.27 mg/dl of the ‘GG’ carriers of the TCF7L2 SNP rs12255372. The results of our study suggest that high fibre intake (35–102 g/day) may help maintain HDL-C level above 40 mg/dl whereas below 35 g/day may lower HDL-C in the ‘GG’ genotype carriers of SNP rs12255372. Mechanisms on how low/high fibre intake decreases/increases HDL-C are not well understood and more research is needed to clarify the effects of dietary fibre on HDL-C metabolism [52, 53].

No studies, to date, have investigated interactions between the TCF7L2 SNPs (rs12255372 and rs7903146) and MC4R SNP rs17782313 and physical activity on cardio-metabolic traits in Asian Indians. Despite the fact that the majority of people are physically inactive in India [54], no significant interactions were found after correction for multiple testing between the three polymorphisms and physical activity on cardio-metabolic traits, which could be due to a small sample size and measurement bias associated with self-reported physical activity questionnaire. However, our finding is in support of the previous study in a Spanish population [21] with a much larger sample size (n = 7,052) which also did not find a significant interaction between MC4R SNP rs17782313 and physical activity on obesity traits. Though the inclusion of a representative sample of Chennai for analysis and the use of a comprehensive, validated, interviewer administered semi-quantitative FFQ for dietary assessment could be considered the strengths of this study, there are some underlying limitations. The data used to calculate the measures of dietary intake and physical activity came from self-report and hence, measurement bias associated with self-reported questionnaire cannot be ruled out. Given that obesity and diabetes are multifactorial traits, several genetic and lifestyle factors are likely to contribute to the disease. While 97 loci have been shown to be associated with body weight [55], the present study examined only three common variants, given their consistent associations with obesity and diabetes, respectively, in Europeans and Asian Indians. The cross sectional study design gives only a snapshot of the prevalence and cause—effect cannot be established due to lack of follow up which is another limitation of this study. However, several outcomes and risk factors were assessed and there was no loss to follow up.

In conclusion, this study has found significant interactions between the TCF7L2 SNP rs12255372 and dietary factors on HDL-C in this Asian Indian population. The results of this study indicate that high total fat and PUFA intakes may be associated with lower HDL-C whereas low intake is associated with higher HDL among the risk allele carriers. More research is required to better understand the interactions between the TCF7L2 gene variant and lifestyle factors on cardio-metabolic traits. Exact mechanisms identifying the effect of different fatty acids on HDL-C and whether/how high fat- and PUFA- intake may reduce HDL-C should also be established before public health recommendations and personalised nutrition advice can be developed for this Asian Indian population in order to reduce the burden of cardiometabolic diseases.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Interactions of TCF7L2 SNPs rs12255372 with fat, PUFA and ALA intakes on HDL-C under an additive and a dominant model of analysis.

https://doi.org/10.1371/journal.pone.0188382.s001

(DOCX)

S2 Table. Associations between TCF7L2 SNPs rs12255372, rs7903146 and obesity, T2D and related traits.

https://doi.org/10.1371/journal.pone.0188382.s002

(DOCX)

S3 Table. Association and interaction between MC4R SNP rs17782313 and obesity, T2D and related traits.

https://doi.org/10.1371/journal.pone.0188382.s003

(DOCX)

S4 Table. Interactions between TCF7L2 SNP rs7903146 and physical activity on VLDL and TG.

https://doi.org/10.1371/journal.pone.0188382.s004

(DOCX)

Citation: Bodhini D, Gaal S, Shatwan I, Ramya K, Ellahi B, Surendran S, et al. (2017) Interaction between TCF7L2 polymorphism and dietary fat intake on high density lipoprotein cholesterol. PLoS ONE 12(11): e0188382. https://doi.org/10.1371/journal.pone.0188382

References

1. Vimaleswaran KS, Loos RJ. Progress in the genetics of common obesity and type 2 diabetes. Expert Reviews in Molecular Medicine. 2010;12(0):e7.

2. Vimaleswaran KS, Li S, Zhao JH, Luan J, Bingham SA, Khaw KT, et al. Physical activity attenuates the body mass index-increasing influence of genetic variation in the FTO gene. The American journal of clinical nutrition. 2009;90(2):425–8. pmid:19553294

3. Vimaleswaran KS, Angquist L, Hansen RD, van der AD, Bouatia-Naji N, Holst C, et al. Association between FTO variant and change in body weight and its interaction with dietary factors: the DiOGenes study. Obesity (Silver Spring). 2012;20(8):1669–74.

4. Anjana RM, Pradeepa R, Deepa M, Datta M, Sudha V, Unnikrishnan R, et al. Prevalence of diabetes and prediabetes (impaired fasting glucose and/or impaired glucose tolerance) in urban and rural India: phase I results of the Indian Council of Medical Research-INdia DIABetes (ICMR-INDIAB) study. Diabetologia. 2011;54(12):3022–7. pmid:21959957

5. Mohan V, Sandeep S, Deepa M, Gokulakrishnan K, Datta M, Deepa R. A diabetes risk score helps identify metabolic syndrome and cardiovascular risk in Indians—the Chennai Urban Rural Epidemiology Study (CURES-38). Diabetes Obesity and Metabolism. 2007;9(3):337–43.

6. Mohan V, Deepa R. Adipocytokines and The Expanding ‘Asian Indian Phenotype’. Journal of the Association of Physicians of India. 2006;54(0):685–6.

7. Vimaleswaran KS, Tachmazidou I, Zhao JH, Hirschhorn JN, Dudbridge F, Loos RJ. Candidate genes for obesity-susceptibility show enriched association within a large genome-wide association study for BMI. Hum Mol Genet. 2012;21(20):4537–42. pmid:22791748

8. Vimaleswaran KS, Bodhini D, Lakshmipriya N, Ramya K, Anjana RM, Sudha V, et al. Interaction between FTO gene variants and lifestyle factors on metabolic traits in an Asian Indian population. Nutrition and Metabolism (London). 2016;13(0):39.

9. Farooqi IS, Keogh JM, Yeo GSH, Lank EJ, Cheetham T, O'Rahilly S. Clinical Spectrum of Obesity and Mutations in the Melanocortin 4 Receptor Gene. The New England Journal of Medicine. 2003;348:1085–95. pmid:12646665

10. Grant SF, Thorleifsson G, Reynisdottir I, Benediktsson R, Manolescu A, Sainz J, et al. Variant of transcription factor 7-like 2 (TCF7L2) gene confers risk of type 2 diabetes. Nature Genetics. 2006;38(3):320–3. pmid:16415884

11. Bauer F, Elbers CC, Adan RA, Loos RJ, Onland-Moret NC, Grobbee DE, et al. Obesity genes identified in genome-wide association studies are associated with adiposity measures and potentially with nutrient-specific food preference. The American journal of clinical nutrition. 2009;90(4):951–9. pmid:19692490

12. Loos RJ, Lindgren CM, Li S, Wheeler E, Zhao JH, Prokopenko I, et al. Common variants near MC4R are associated with fat mass, weight and risk of obesity. Nature Genetics. 2008;40(6):768–75. pmid:18454148

13. Sull JW, Lee M, Jee SH. Replication of genetic effects of MC4R polymorphisms on body mass index in a Korean population. Endocrine. 2013;44(3):675–9. pmid:23460509

14. Renstrom F, Payne F, Nordstrom A, Brito EC, Rolandsson O, Hallmans G, et al. Replication and extension of genome-wide association study results for obesity in 4923 adults from northern Sweden. Human Molecular Genetics. 2009;18(8):1489–96. pmid:19164386

15. Zobel DP, Andreasen CH, Grarup N, Eiberg H, Sorensen TI, Sandbaek A, et al. Variants near MC4R are associated with obesity and influence obesity-related quantitative traits in a population of middle-aged people: studies of 14,940 Danes. Diabetes. 2009;58(3):757–64. pmid:19073769

16. Vasan SK, Fall T, Neville MJ, Antonisamy B, Fall CH, Geethanjali FS, et al. Associations of variants in FTO and near MC4R with obesity traits in South Asian Indians. Obesity (Silver Spring). 2012;20(11):2268–77.

17. Tong Y, Lin Y, Zhang Y, Yang J, Zhang Y, Liu H, et al. Association between TCF7L2 gene polymorphisms and susceptibility to type 2 diabetes mellitus: a large Human Genome Epidemiology (HuGE) review and meta-analysis. BMC Medical Genetics. 2009;10(0):15.

18. Lyssenko V, Lupi R, Marchetti P, Del Guerra S, Orho-Melander M, Almgren P, et al. Mechanisms by which common variants in the TCF7L2 gene increase risk of type 2 diabetes. Journal of Clinical Investigation. 2007;117(8):2155–63. pmid:17671651

19. Chandak GR, Janipalli CS, Bhaskar S, Kulkarni SR, Mohankrishna P, Hattersley AT, et al. Common variants in the TCF7L2 gene are strongly associated with type 2 diabetes mellitus in the Indian population. Diabetologia. 2007;50(1):63–7. pmid:17093941

20. Bodhini D, Radha V, Dhar M, Narayani N, Mohan V. The rs12255372(G/T) and rs7903146(C/T) polymorphisms of the TCF7L2 gene are associated with type 2 diabetes mellitus in Asian Indians. Metabolism. 2007;56(9):1174–8. pmid:17697858

21. Corella D, Ortega-Azorin C, Sorli JV, Covas MI, Carrasco P, Salas-Salvado J, et al. Statistical and biological gene-lifestyle interactions of MC4R and FTO with diet and physical activity on obesity: new effects on alcohol consumption. PloS one. 2012;7(12):e52344. pmid:23284998

22. Hasselbalch AL, Angquist L, Christiansen L, Heitmann BL, Kyvik KO, Sorensen TI. A variant in the fat mass and obesity-associated gene (FTO) and variants near the melanocortin-4 receptor gene (MC4R) do not influence dietary intake. Journal of Nutrition. 2010;140(4):831–4. pmid:20181787

23. Holzapfel C, Grallert H, Huth C, Wahl S, Fischer B, Doring A, et al. Genes and lifestyle factors in obesity: results from 12,462 subjects from MONICA/KORA. Int J Obes (Lond). 2010;34(10):1538–45.

24. Warodomwichit D, Arnett DK, Kabagambe EK, Tsai MY, Hixson JE, Straka RJ, et al. Polyunsaturated fatty acids modulate the effect of TCF7L2 gene variants on postprandial lipemia. Journal of Nutrition. 2009;139(3):439–46. pmid:19141698

25. Wirstrom T, Hilding A, Gu HF, Ostenson CG, Bjorklund A. Consumption of whole grain reduces risk of deteriorating glucose tolerance, including progression to prediabetes. The American journal of clinical nutrition. 2013;97(1):179–87. pmid:23235198

26. Taylor AE, Sandeep MN, Janipalli CS, Giambartolomei C, Evans DM, Kranthi Kumar MV, et al. Associations of FTO and MC4R Variants with Obesity Traits in Indians and the Role of Rural/Urban Environment as a Possible Effect Modifier. Journal of Obesity. 2011;2011(0):307542.

27. Qi L, Kraft P, Hunter DJ, Hu FB. The common obesity variant near MC4R gene is associated with higher intakes of total energy and dietary fat, weight change and diabetes risk in women. Human Molecular Genetics. 2008;17(22):3502–8. pmid:18697794

28. Hindy G, Sonestedt E, Ericson U, Jing XJ, Zhou Y, Hansson O, et al. Role of TCF7L2 risk variant and dietary fibre intake on incident type 2 diabetes. Diabetologia. 2012;55(10):2646–54. pmid:22782288

29. Phillips CM, Goumidi L, Bertrais S, Field MR, McManus R, Hercberg S, et al. Dietary saturated fat, gender and genetic variation at the TCF7L2 locus predict the development of metabolic syndrome. Journal of Nutritional Biochemistry. 2012;23(3):239–44. pmid:21543200

30. Mattei J, Qi Q, Hu FB, Sacks FM, Qi L. TCF7L2 genetic variants modulate the effect of dietary fat intake on changes in body composition during a weight-loss intervention. The American Journal of Clinical Nutrition. 2012;96(5):1129–36. pmid:23034957

31. Deepa M, Pradeepa R, Rema M, Mohan A, Deepa R, Shanthirani S, et al. The Chennai Urban Rural Epidemiology Study (CURES)—Study Design And Methodology (Urban Component) (CURES—1). Journal of the Association of Physicians of India. 2003;51(0):863–70.

32. Alberti KG, Zimmet PZ. Definition, Diagnosis and Classification of Diabetes Mellitus and its Complications Part 1: Diagnosis and Classification of Diabetes Mellitus Provisional Report of a WHO Consultation. Diabetic Medicine. 1998;15(7):539–53. pmid:9686693

33. Friedewald WT, Levy RI, Fredrickson DS. Estimationof the Concentration of Low-Density Lipoprotein Cholesterol in Plasma, Without Use of the Preparative Ultracentrifuge. Chlinical Chemistry. 1972;18(0):499–502.

34. Sudha V, Radhika G, Sathya RM, Ganesan A, Mohan V. Reproducibility and validity of an interviewer-administered semi-quantitative food frequency questionnaire to assess dietary intake of urban adults in southern India. International Journal of Food Science and Nutrition. 2006;57(7–8):481–93.

35. Miyake K, Horikawa Y, Hara K, Yasuda K, Osawa H, Furuta H, et al. Association of TCF7L2 polymorphisms with susceptibility to type 2 diabetes in 4,087 Japanese subjects. Journal of Human Genetics. 2008;53(2):174–80. pmid:18097733

36. Florez JC. The new type 2 diabetes gene TCF7L2. Current opinion in clinical nutrition and metabolic care. 2007;10(4):391–6. pmid:17563454

37. Ip W, Chiang YT, Jin T. The involvement of the wnt signaling pathway and TCF7L2 in diabetes mellitus: The current understanding, dispute, and perspective. Cell & bioscience. 2012;2(1):28.

38. Misra R, Patel T, Kotha P, Raji A, Ganda O, Banerji M, et al. Prevalence of diabetes, metabolic syndrome, and cardiovascular risk factors in US Asian Indians: results from a national study. Journal of Diabetes Complications. 2010;24(3):145–53.

39. Joshi SR, Anjana RM, Deepa M, Pradeepa R, Bhansali A, Dhandania VK, et al. Prevalence of dyslipidemia in urban and rural India: the ICMR-INDIAB study. PloS one. 2014;9(5):e96808. pmid:24817067

40. Research ICoM. NUTRIENT REQUIREMENTS AND RECOMMENDED DIETARY ALLOWANCES FOR INDIANS. Hyderabad: NATIONAL INSTITUTE OF NUTRITION, 2010.

41. Norby FL, Eryd SA, Niemeijer MN, Rose LM, Smith AV, Yin X, et al. Association of Lipid-Related Genetic Variants with the Incidence of Atrial Fibrillation: The AFGen Consortium. PloS one. 2016;11(3):e0151932. pmid:26999784

42. Voight BF, Peloso GM, Orho-Melander M, Frikke-Schmidt R, Barbalic M, Jensen MK, et al. Plasma HDL cholesterol and risk of myocardial infarction: a mendelian randomisation study. Lancet. 2012;380(9841):572–80. pmid:22607825

43. Zhang L, Geng Y, Xiao N, Yin M, Mao L, Ren G, et al. High Dietary n-6/n-3 PUFA Ratio Promotes HDL Cholesterol Level, but does not Suppress Atherogenesis in Apolipoprotein E-Null Mice 1. Journal of Atherosclerosis and Thrombosis. 2009;16(0):463–71. pmid:19729862

44. Riccardi G, Bozzetto L, Annuzzi G. Postprandial lipid metabolism. Scandinavian Journal of Food & Nutrition. 2006;50(1):99–106.

45. Thomsen C, Storm H, Holst JJ, Hermansen K. Differential effects of saturated and monounsaturated fats on postprandial lipemia and glucagon-like peptide 1 responses in patients with type 2 diabetes. The American journal of clinical nutrition. 2003;77(3):605–11. pmid:12600850

46. Thomsen C, Rasmussen O, Lousen T, Holst JJ, Fenselau S, Schrezenmeir J, et al. Differential effects of saturated and monounsaturated fatty acids on postprandial lipemia and incretin responses in healthy subjects. The American journal of clinical nutrition. 1999;69(6):1135–43. pmid:10357731

47. Rasmussen O, Lauszus FF, Christiansen C, Thomsen C, Hermansen K. Differential effects of saturated and monounsaturated fat on blood glucose and insulin responses in subjects with non-insulin-dependent diabetes mellitus. The American journal of clinical nutrition. 1996;63(2):249–53. pmid:8561067

48. Simopoulos AP. An Increase in the Omega-6/Omega-3 Fatty Acid Ratio Increases the Risk for Obesity. Nutrients. 2016;8(3):128. pmid:26950145

49. Silva FM, Kramer CK, de Almeida JC, Steemburgo T, Gross JL, Azevedo MJ. Fiber intake and glycemic control in patients with type 2 diabetes mellitus: a systematic review with meta-analysis of randomized controlled trials. Nutrition Reviews. 2013;71(12):790–801. pmid:24180564

50. Brown L, Rosner B, Willett WW, Sacks FM. Cholesterol-lowering effects of dietary fiber: a meta-analysis. The American journal of clinical nutrition. 1999;69:30–42. pmid:9925120

51. Anderson JW, Randles KM, Kendall CWC, Jenkins DJA. Carbohydrate and Fiber Recommendations for Individuals with Diabetes: A Quantitative Assessment and Meta-Analysis of the Evidence. Journal of the American College of Nutrition. 2004;23(1):5–17. pmid:14963049

52. Yanai H, Katsuyama H, Hamasaki H, Abe S, Tada N, Sako A. Effects of Carbohydrate and Dietary Fiber Intake, Glycemic Index and Glycemic Load on HDL Metabolism in Asian Populations. Journal of Clinical Medicine Research. 2014;6(5):321–6. pmid:25110535

53. Choi H, Song S, Kim J, Chung J, Yoon J, Paik HY, et al. High carbohydrate intake was inversely associated with high-density lipoprotein cholesterol among Korean adults. Nutrition Research. 2012;32(2):100–6. pmid:22348458

54. Anjana RM, Pradeepa R, Das AK, Deepa M, Bhansali A, Joshi SR, et al. Physical activity and inactivity patterns in India–results from the ICMR-INDIAB study (Phase-1) [ICMR-INDIAB-5]. International Journal of Behavioral Nutrition and Physical Activity. 2014;11(1):26. pmid:24571915

55. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518(7538):197–206. pmid:25673413