Next Article in Journal
Regulatory Effects of Quercetin on M1/M2 Macrophage Polarization and Oxidative/Antioxidative Balance
Next Article in Special Issue
Influence of a Nutrigenetic Intervention on Self-Efficacy, Emotions, and Rewarding Behaviors in Unhealthy Eating among Mexicans: An Exploratory Pilot Study
Previous Article in Journal
Dietary Supplementation for Attenuating Exercise-Induced Muscle Damage and Delayed-Onset Muscle Soreness in Humans
Previous Article in Special Issue
High-Salt Diet Impairs the Neurons Plasticity and the Neurotransmitters-Related Biological Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Association between Fasting Glucose and Sugar Sweetened Beverages Intake Is Greater in Latin Americans with a High Polygenic Risk Score for Type 2 Diabetes Mellitus

by
María Lourdes López-Portillo
1,
Andrea Huidobro
2,3,
Eduardo Tobar-Calfucoy
1,
Cristian Yáñez
1,
Rocío Retamales-Ortega
1,
Macarena Garrido-Tapia
4,
Johanna Acevedo
4,
Fabio Paredes
4,
Vicente Cid-Ossandon
4,
Catterina Ferreccio
3,4 and
Ricardo A. Verdugo
1,5,*
1
Programa de Genética Humana, ICBM, Facultad de Medicina, Universidad de Chile, Santiago 8380000, Chile
2
Departamento de Ciencias Preclínicas, Facultad de Medicina, Universidad Católica del Maule, Talca 3460000, Chile
3
Advanced Center for Chronic Diseases, ACCDiS, Santiago 8380000, Chile
4
Facultad de Medicina, Pontificia Universidad Católica de Chile, Santiago 8328888, Chile
5
Departamento de Oncología Básico-Clínica, Facultad de Medicina, Universidad de Chile, Santiago 8380000, Chile
*
Author to whom correspondence should be addressed.
Nutrients 2022, 14(1), 69; https://doi.org/10.3390/nu14010069
Submission received: 21 October 2021 / Revised: 16 December 2021 / Accepted: 17 December 2021 / Published: 24 December 2021
(This article belongs to the Special Issue The Effect of Gene-Diet Interactions in Human Health)

Abstract

:
Chile is one of the largest consumers of sugar-sweetened beverages (SSB) world-wide. However, it is unknown whether the effects from this highly industrialized food will mimic those reported in industrialized countries or whether they will be modified by local lifestyle or population genetics. Our goal is to evaluate the interaction effect between SSB intake and T2D susceptibility on fasting glucose. We calculated a weighted genetic risk score (GRSw) based on 16 T2D risk SNPs in 2828 non-diabetic participants of the MAUCO cohort. SSB intake was categorized in four levels using a food frequency questionnaire. Log-fasting glucose was regressed on SSB and GRSw tertiles while accounting for socio-demography, lifestyle, obesity, and Amerindian ancestry. Fasting glucose increased systematically per unit of GRSw (β = 0.02 ± 0.006, p = 0.00002) and by SSB intake (β[cat4] = 0.04 ± 0.01, p = 0.0001), showing a significant interaction, where the strongest effect was observed in the highest GRSw-tertile and in the highest SSB consumption category (β = 0.05 ± 0.02, p = 0.02). SNP-wise, SSB interacted with additive effects of rs7903146 (TCF7L2) (β = 0.05 ± 0.01, p = 0.002) and with the G/G genotype of rs10830963 (MTNRB1B) (β = 0.19 ± 0.05, p = 0.001). Conclusions: The association between SSB intake and fasting glucose in the Chilean population without diabetes is modified by T2D genetic susceptibility.

1. Introduction

Sugar-sweetened beverages (SSB) have high quantity of rapidly absorbed sugars and poor nutritional value [1]. Chile is one of the highest per capita consumers of sugar-sweetened beverages (SSB) at international level [2,3]. However, the biological impacts from such frequent consumption in this population and how these may be modified by its particular genetic composition is unknown. Chileans are largely of admixed ancestry from European (55%), Amerindian (42%), and African (3%) origin [4]. Amerindian ancestry is related to the Aymara and Mapuche ethnic groups and proportions depends on latitude [5]. Further, the population is structured in a way that has preserved a socio-economic cline associated with ancestry [6].
SSB consumption has been related to increased blood insulin concentration, impaired fasting glucose (IFG) [7,8], glucose intolerance (GI), insulin resistance (IR) [9,10,11] and high risk of type 2 diabetes (T2D) [12,13,14]. IFG is considered as a physical state of evolution towards T2D [15], so its control is essential to prevent the T2D onset and cardiovascular diseases [16]. There is also evidence that SSB promote weight gain and increased adiposity due to their high energy content and because they promote higher caloric intake [17,18,19].
Genome-wide studies (GWAS) have found more than 140 T2D variants in independent loci but with a modest association that even in combination explain a low proportion (~20%) of disease liability [20]. Presence pervasive genotype-by-environment interactions (GxE) and rare alleles with large effects, which would be missed from GWAS focusing in common variants, have been presented as possible explanations for the missing heritability [21,22]. However, large whole-genome and whole-exome efforts have found only a handful of rare variants that explain even a smaller fraction of phenotypic variation [23,24]. The relevance of GxE in predisposition to T2D could be used to detect groups at higher risk and indicate that interventions are best conceived in a population-specific or even individual-specific manner [25].
There is interest in identifying genetic loci implicated in fasting glucose and lifestyle-gene interactions that modify the incidence of T2D because they may suggest relevant biological pathways [26,27]. Interaction between fruit intake and polygenic genetic risk for type 2 diabetes have shown an effect on both diabetes risk and on fasting glucose [28]. Similar studies have explored interaction with SSB, without conclusive results [29]. A study on subjects of European origin evaluated the interaction of a T2D-GRS and foods related to T2D risk (including sugar-sweetened beverages) but did not have significant results [30]. A meta-analysis of 11 cohorts including 34,748 participants without type 2 diabetes that included loci related to the CHREB pathway of fructose metabolism found a significant association of SSB consumption with insulin concentration and fasting glucose, but no interaction on was found [8].
However, consumption of carbohydrate-rich foods has shown interaction with individual T2D-risk variants, among which Transcription Factor 7 Like-2 (TCF7L2) has been consistently replicated [25]. For instance, a prospective study including 5477 participants of European origin found that a higher intake of whole grain was associated with a 34% lower risk to deteriorate in glucose tolerance to prediabetes or T2D. They show that the protective effect of whole-grain intake was lost in subjects who carried the T allele of rs7903146 in TCF7L2 [31]. Recently, in a Swedish case-control study conducted in a European population investigate whether the association between sweetened beverages and T2D is modified by TCF7L2 but there were no indications of synergistic effects [32].
The present study aims to evaluate the association between sugar-sweetened beverage consumption and fasting glucose, and to establish whether this relationship is affected by Chilean genetic variation. In a sample of subjects participating in the MAUCO cohort without T2D self-reporting, we investigated the presence of interaction effects between sugar-sweetened beverage consumption and T2D genetic susceptibility considering ancestry and obesity. We provide the first genotype-SSB assessment of interaction effect on fasting glucose in a Latin American population.

2. Materials and Methods

2.1. Study Population and Data Collection

MAUCO is a prospective cohort study initiated in 2014 in the city of Molina, Maule Region, 200 km south of Chile’s capital Santiago in which 8970 eligible participants were enrolled. The aim of this project is to study the natural history of some chronic diseases and Cancer and associated risk factors. In brief, all adult inhabitants of Molina were invited to participate. A baseline health interview was applied at home, where a validated food frequency questionnaire of Mediterranean diet was applied [33].
The survey also assessed demography, socioeconomic status, educational level, lifestyle factors such as frequency of 30 min of physical activity (three or more times per week, once or twice per week, less than 4 times per month), smoking status, and alcohol use. Subsequently, all participants were called to a health station to collect anthropometric, biochemical measurements and collection of biological samples (urine, saliva, and blood). Waist circumference (WC) was measured with a non-distensible fiberglass tape and was classified as normal or high according to the Adult Treatment Panel III (ATPIII) criteria. Weight and height were measured with a SECA 700 scale with height rod and a up to 200 cm. Body mass index (BMI) > 25 and <30 points were classified as overweight and those with 30 or more as obese. The biochemical measurements obtained for this analysis were HDL cholesterol (mg/dL) and triglycerides (mg/dL), which were classified as normal or low/high according to ATPIII criteria. Fasting glucose was measured in venous blood after 8 to 12 h of fasting using the glucose oxidase method in a CT 600i spectrophotometer (Weiner Lab, Rosario, Argentina). The detailed methodology is described elsewhere [33].
We conducted a cross-sectional study with a sample of 3321 subjects aged 36–77 years participating in the MAUCO cohort without self-reported T2D. subjects with a genotyping rate of less than 80% of selected SNPs and without SSB consumption register were eliminated. The final sample was 2828 subjects with complete genomic and environmental information. A workflow representation of sample sizes and data filtering is shown in Figure 1. All participants who met the criteria information for this investigation signed the informed consent form. This research was approved by the Ethics Committee of the Pontificia Universidad Católica de Chile and the Ministry of Health of the Maule Regional Health Services.

2.2. Evaluation of Sugar-Sweetened Beverage Consumption

Dietary information was obtained with a semi-quantitative food frequency questionnaire of Mediterranean diet applied by previously trained personnel. The average consumption per day during the last year was specified according to 4 levels as follows: category 1:no consumption (0 = 0), category 2: >0 and <1 servings/day, category 3: 1≤ and <2 servings/day, and category 4: ≥2 servings/day. One serving is equivalent to 350 mL. Sugar-sweetened beverages were defined by containing caloric sweeteners such as sucrose or high fructose corn syrup, providing the following examples: regular soda, industrialized fruit juices with natural or artificial ingredients, energy drinks, and sports drinks. “Light”, “zero”, “low”, or “reduced calorie” beverages as well as natural juices were excluded. High-risk alcohol consumption was defined by weekly servings of ≥30 for men ≥20 for women.

2.3. Selection of T2D Risk SNPs and Ancestry

Single nucleotide polymorphisms (SNPs) for T2D risk were obtained from a Genome Wide Association Studies (GWAS) for T2D and from meta-analyses. SNPs were selected by statistical significance of ≤3 × 10−8. From these, the set of not redundant (not in LD) were selected, while prioritizing those with the largest OR (odds ratio) and those that were assessed in populations of Latin American ancestry. A total of 21 SNPs were selected, of which 17 resulted in successful PCR assays (Supplementary Table S1).

2.4. Genotyping

The 17 T2D SNPs were genotyped separately from 150 ancestry informative markers (AIMs) in two multiplexed reactions using a modified version of GT-Seq protocol [34]. Briefly, genomic DNA was distributed in 96-well plates and amplified using the PCR (polymerase chain reaction) technique (PCR1) with locus-specific primers. A second reaction (PCR2) added Illumina sequencing adapters and index pairs to uniquely identify wells and plates. Next, amplicon concentration across wells was homogenized using SequalPrep™ Normalization Plate Kit (Applied Biosystems, Waltham, MA, USA) before mixing samples into one pool per plate. Finally, an equimolar mixture of pools was created and sequenced in a NextSeq 550 (Illumina Inc., San Diego, CA, USA) with Mid Output v2 kits in 1 × 125 configuration. Genotypes for 17 T2D SNPs and 150 AIMS were successfully obtained but rs188827514 was excluded from the analysis because it had a minor allele frequency less than 5%.

2.5. Bioinformatic Analysis

The bcl files generated were converted to a fastq format with bcl2fastq v2.20.0.422 and quality-controlled with the fastqc v11.5. Primer homodimers and heterodimers were removed and one fastq file per sample was generated according to the sequence of the i7 (plate) and the i5 (well) indices by using scripts available at https://github.com/GTseq/GTseq-Pipeline (accessed on 26 January 2018). On-target reads, were defined as having both the sequence of locus-specific forward primer and a known sequence of 17 base pairs around the selected SNP. The genotype assignment was performed by counting reads matching to each allele of the selected SNPs in each sample.

2.6. Ancestry Inference

A set reference dataset with 224 samples was created by joining whole-genome sequencing (WGS) data for 30 individuals of African (YRI) ancestry and 30 of European (CEU) ancestry from the 1000 Genomes Project and Axiom LAT1 microarray data from 110 Aymaras from Puno, Peru (Andrés Moreno-Estrada, personal communication), 54 high Amerindian ancestry from central and southern Chile from the ChileGenomico and PatagoniaDNA projects [5,35]. Global ancestry proportion per individual was calculated with the ADMIXTURE program (version 1.23) counting the reference group allele frequencies and those of the mestizo sample under study. A script was implemented to calculate the ancestry proportion in thousands of individuals by creating bins of 50 MAUCO participants in order to minimize bias due to the large sample size of this dataset.

2.7. Polygenic Risk Score Calculation

All T2D SNPs that met the quality criteria were included to construct a weighted Genetic Risk Score (GRSw) as follows: first, a count of T2D-associated risk alleles of each polymorphism was performed (0 = protective homozygous, 1 = heterozygous, 2 = risk homozygous). The GRSw calculation was performed adding the Beta effect (β) obtained from the natural logarithm of the odds ratio (OR) reported in genome-wide studies (GWAS) performed in Latinos or multi-ethnic studies that considered this ancestry, multiplied by the number of risk alleles carried by each subject as follows: GRSwi = j = 1 k β j x i , where i is a subject index, k is the number independent SNPs genotyped in each subject, β i is the estimated effect for the risk allele of each SNP j and x i is the number of risk alleles in subject i. β i was obtained from publications. The odds ratio (OR) was reported, β i wasestimated as ln O R . If the OR for the protective alleles was reported, ln 1 / O R was used.

2.8. Statistical Analysis

The output variable (fasting glucose) was log transformed in base e in order to approximate a normal distribution. Univariate association with categorical variables was tested by the Chi-square tests. Association with continuous variables was assessed with the Student’s t-test, ANOVA or Kruskal–Wallis for non-parametric variables.
Association between fasting glucose-log and SSB, GRSw and individual SNPs was evaluated using multiple linear regression models, while adjusting for age, sex, waist circumference, physical activity, educational and socioeconomic level, smoking status, high consumption of fruits, vegetables, sugar, processed meats, and alcoholic beverages (as defined in Table 1), Amerindian ancestry tertile, and GRSw. GRSw was included both as a continuous variable and as a categorial factor by classifying individuals into tertiles of GRSw. All two-way interaction terms between sex, SSB, Ancestry, and GRSw where considered. When testing interaction between SSB and individual SNPs, no GRSw term was included and all two-way interaction terms between sex, SSB, Ancestry, and SNP where considered. One SNP at a time was included. p values were calculated from t-tests for individual each regression coefficient (H0:β = 0) and from F-tests for each factor by comparing the full model with a reduced model that dropped a single term using the drop1 function in R. Significance was set at 0.05 and Bonferroni correction was used when testing association and interactions for individual SNPs only. The proportion of variance explained by each SNP was calculated by R2 = SSSNP/SSTot, were SSSNP is the sum of squares due to the SNP’s marginal and interaction effects and SSSNP is the total sum of squares, i.e., i = 1 n y i y ¯ 2 .
Post-hoc power estimation to detect (replicate) genotype-to-phenotype associations were performed using the online web tool (https://clincalc.com/stats/power.aspx, accessed on 14 December 2021) based on the marginal effects of association between SNPs (categorical) and fasting glucose (log) (continuous). All statistical analyses were performed with RStudio version 1.31093.

3. Results

3.1. General Characteristic of the Study Population

The studied sample consisted of 3866 subjects aged 36 to 77 years, 32.2% (n = 1247) are male and 67.8% (n = 2619) are female with an average age of 54(9.1) years of which 42.7% (n = 1652) with overweight (SP) and 40.9% (n = 1580) with obesity (OB) so the combined prevalence of SP and OB is 83.6% (n = 3247). The prevalence of central obesity measured by waist circumference was 71.2% (n = 2754) and 14.1% (n = 545) reported being diabetic or taking glucose-lowering medication.
The final sample size including complete genomic, anthropometric, socio-demographic and lifestyle information was 2828 subjects. General characteristics and their association with fasting glucose in the whole sample and stratified by sex are shown in Table 1 and Supplementary Tables S2 and S3.

3.2. General Characteristics by SSB Consumption

Results show significant differences in general characteristics of participants by intake level. High SSB consumption was more frequent in men (p < 0.0001) and in younger individuals (p < 0.03), active smokers (p < 0.001), and risky alcohol drinkers (p < 0.0001). In addition, consumers of sugar-sweetened beverages reported lower consumption of fruits and vegetables (p < 0.0001) and higher intake of sugar and processed meats in their diet (p < 0.0001). Progressive increase on fasting blood glucose-log and WC are observed for each increase in SSB category of consumption (p < 0.0001) (Supplementary Table S4).
SSB consumption was associated with higher log-fasting glucose when comparing the highest (≥2 servings/day) versus the lowest category (0 servings/day) of intake in the whole sample (β = 0.04 ± 0.01, p < 0.0001) and in women (β = 0.05 ± 0.01, p < 0.001). When BMI was excluded from models, the association remained significant in the whole sample (Supplementary Table S5). No significant association was observed in males. There was also no evidence of a sex-by-SSB interaction (p > 0.51) (Supplementary Table S2).

3.3. Fasting Glucose and SNPs: Association and Interaction Effect

Among the 17 selected SNPs, rs188827514 (CCND2) was excluded because it had a MAF (Minor Allele Frequency) lower than 5%. The association between each SNPs that had significant association with T2D in GWAS studies and fasting glucose-log. SNPs were located in 12 different chromosomes, and all are distant enough so that they are nonredundant (r2 < 0.80).
Our results showed a positive association between log-fasting glucose and variants rs516946 in Ankyrin 1 (ANK1), rs4402960 in Insuline-Like Growth Factor 2 mRNA Binding Protein (IGF2BP2), and rs7903146 in Transcription Factor 7 Like 2 (TCF7L2) genes, while adjusting for age, sex, BMI and Amerindian Ancestry, after Bonferroni correction (p < 0.0004, 0.0007 and 0.0009, respectively) (Supplementary Table S6). The association remain significative when BMI is excluded from the model, indicating they are not mediated by adiposity.
An interaction effect was found between rs7903146 (TCFL2) (as continuous) and the highest SSB category (β= 0.05 ± 0.01, p < 0.002) (Supplementary Table S7). When regressed by number of risk alleles, significance was lost after Bonferroni correction (p < 0.005) (Table 2 and Supplementary Table S8). When genotypes were coded as a categorical variable, rs10830963 Melatonin Receptor 1B (MTNR1B) showed significant positive interaction between the G/G genotype of rs10830963 (MTNR1B) variant and the highest category of SSB consumption (cat4) on log-fasting glucose in the overall sample after Bonferroni correction (β = 0.19 ± 0.05, p < 0.001). Post-hoc power estimations showed that our sample size provided ≥80% power to detect effects of OR ≥ 1.01 for genotypes with high frequency genotyping rate, but it was limited for other SNPs (Supplementary Table S9). Results were unaffected when BMI was excluded from the models (data not shown). The interaction results for these two SNPs are found in Table 2 and for all SNPs in Supplementary Table S8.

3.4. Fasting Glucose and GRSw: Association and Interaction Effect

We constructed an effect-weighted risk score from the 16 selected SNPs of T2D risk. β estimates for association (used in the score) came from studies including Latin American participants whenever possible. The GRSw had an average value of 4.16 (SD = 0.48, min 2.02, max 5.85) and a gaussian distribution (not shown). Linear regression accounting for demographics, lifestyles and ancestry showed that fasting glucose-log increased 0.02 per unit of GRSw (β = 0.02 ± 0.006, p < 0.00002) in the whole sample, i.e., increased 1.02 mg/dL for each GRSw increment (Table 3). When GRSw was discretized in tertiles, participants in the highest tertile incremented log-glycemia also in 0.2 versus the lowest tertile (β = 0.02 ± 0.007, p < 0.0003) (Table 3). The effect was marginally larger in females (β = 0.03 ± 0.008, p < 0.0024) than in males (β = 0.02 ± 0.01, p < 0.08) but no significant interaction was observed between sex and GRSw (Supplementary Table S2). Results remained significant when BMI removed from the models (data not shown).
We observed a tendency of increasing effects from SSB consumption as the tertile of GRSw increased (Figure 2). There was a significant and positive interaction between continuous SSB and continuous GRSw (β = 0.02 ± 0.006, p < 0.004) (Table 4). That is, the increment in log-glycemia by each level of consumption is 2% higher with every unit of increment in GRSw (e0.02 = 1.02). When SSB was considered as categorical variable, the trend by GRSw was higher in the highest category of SSB consumption in the whole sample (β = 0.06 ± 0.02, p < 0.01) and in men (β = 0.08 ± 0.04, p < 0.03) (Table 4). Although the interaction was significant only in the global sample and in men, we did not detect a significant interaction of sex with either SSB or GRSw (Supplementary Table S2). When both SSB and GRSw where rendered categorical, the interaction effect was only detected in the highest levels of SSB consumption (cat4) and GRSw (tertil3) in the whole sample (β = 0.05 ± 0.02, p = 0.02). The interaction was also observed in men (β = 0.07 ± 0.05, p = 0.009) (Table 4). The effect modification remained when including or excluding BMI as a covariate (not shown).

3.5. Sensitivity Analysis

After excluding subjects older than 65 years old, who were at higher risk because age is a risk factor, the interaction effects were attenuated. However, the interaction between sugary drinks consumption and GRSw (as continuous) was preserved (p < 0.04). Similarly, the P trend of the interaction between SSB consumption and rs7903146 (TCF7L2) was significant in double risk allele carriers (p < 0.003), but significance is lost after Bonferroni correction (p < 0.03). In a second analysis, including two levels of consumption were considered (less than 1 serving/day = 0, and 1 or more servings/day = 1) in order to balance sample sizes. Interaction remains significant, and was accentuated in the subjects belonging to the highest category of consumption (Supplementary Tables S10–S13).

4. Discussion

Here we preset results from a cross-sectional study that included 2828 participants from the MAUCO cohort without self-reported T2D diagnosis. The positive association of sugar-sweetened beverages and fasting glucose has been inconsistent in the international literature. However, majority of studies and a recent meta-analysis tended to support the unfavorable effect of these beverages on glycemic traits, including fasting glucose [7,8,36].
In agreement with previous studies our results indicate a positive association between SSB consumption and fasting glucose when adjusting for all possible confounders. Furthermore, significance was preserved when adjusting for GRSw and proportion of Amerindian ancestry. The general characteristics clearly showed that consumption of fruits and vegetables decrease while that of processed meats increase when SSB intake is higher. The same trend of increase was observed on fasting glucose levels, obesity, waist circumference and triglycerides. HDL-c cholesterol decrease when SSB consumption increases.
We used a weighted polygenic risk score (GRSw) based on 16 T2D SNPs previously observed in Latinos. Most of them have been related to HOMA-β which is an indirect measure of the first phase of insulin secretion [37,38,39]. Results showed linear association between GRSw and fasting glucose while accounting for clinical, socio-demographic, and lifestyle factors, as well as Amerindian ancestry. The highest glucose levels correspond to subjects with the highest DM2 genetic risk. That is, the score increases as blood glucose rises, so is therefore considered to be a useful tool for measuring genetic susceptibility.
SSB-FG association was intensified in subjects with the highest consumption and the highest GRSw. In other words, the association between SSB consumption and fasting glucose was amplified in subjects with the highest genetic susceptibility that consume two or more servings of SSB per day. All our results were not change when BMI was excluded from models, suggesting that this genetic effect is not mediated by obesity.
When analyzing each selected SNP, we found association between rs7903146 (TCF7L2), rs4402960 (IGF2BP2), rs516946 (ANK1) and fasting glucose. TCF7L2 is a transcription factor located in the non-coding region of Chromosome 10 and the association with fasting glucose was previously observed in Chileans by Peterman and coworkers [40]. This variant have the greatest effect on T2D risk and is recognized as the most important regulator of pro-insulin expression and processing [41]. It has been observed that TCF7L2 influences regulation of glucose metabolism through the Wnt signaling pathway [42] and that carriers of the double risk allele (T/T) were more likely to have progression from impaired glucose tolerance to diabetes than were CC homozygotes and decreased insuline secretion [43]. IGF2BP2 (insulin like growth factor 2 MRNA binding protein 2) gene encodes a protein that binds the 5′ UTR of insulin-like growth factor 2 (IGF2) mRNA and regulates insulin translation. This variant increase T2D risk through an insulin secretory mechanism, including the lower first-phase insulin response [44]. ANK1 is a protein coding gene called ankyrin-1 expressed mainly in red blood cells and has been associated with HbA1c (glycated hemoglobin) [45]. ANK1 has been found in skeletal muscle but the mechanism by which sAnk1 in skeletal muscle might be involved in T2D remains uncertain [46].
When considering diet-dependent genetic effects, we found positive interaction between SSB consumption and T2D-SNPs on fasting glucose (FG). In two of the sixteen selected SNPs, rs7903146 (TCF7L2) and rs10830963 (MTNR1B), a positive interaction with SSB was observed. The largest effect was observed among homozygous risk carriers belonging to the highest category of SSB consumption. This research reports the first SSB-genotype interaction analysis in Latin American subjects. Therefore, this study represents an approximation to the knowledge of the genetic risk on glycemic balance and the genes involved in Latinos.
The first scientific evidence on diet-gene interaction has been provided by studies that have employed BMI risk scores to measure the genetic effect on adiposity level [47,48,49]. There are few studies that have evaluated the food-gene interaction when the outcome is T2D or its quantitative traits. In a prospective nested case-control study including 1196 diabetic and 1337 non-diabetic men a significant interaction (p < 0.02) was observed between T2D-GRS and Western dietary pattern on T2D risk. The effect were more evident among men with high genetic score and the highest quartile of Western dietary pattern [50]. Ericson and coworkers examined the interaction between a GRS constructed from T2D-48 SNPs, and a diet risk score of four foods associated with T2D (processed meat, sugar-sweetened beverages, whole grain and coffee) in 25,069 individuals of European origin, but no significant interaction between the GRS or food components was observed [30].
Previous research showed that variants in MTNR1B were consistently associated with fasting glucose and the strongest signal was observed at rs10830963, where each G allele was associated with an increase in fasting glucose and reduced beta-cell function [51]. A meta-analysis evaluating whether diet modifies the association between circadian-related variants and cardio-metabolic risk factors subjects of European origin showed the protective association of a higher carbohydrate intake in rs1387153 (MTNR1B) was weaker for each additional T allele. That is, the T allele attenuates the inverse association between carbohydrate intake and FG [52]. Finally, in a randomized cross-over trial evaluating glucose-tolerance in late-dinner and early-dinner of homozygous carriers and non-carriers of the MTNR1B risk allele suggest that moving the dinner to an earlier time may result in better glucose-tolerance specially in MTNR1B carriers [53].
Other researchers have also been interested in testing SNP-by-SSB interaction. A recent meta-analysis conducted by the CHARGE Consortium in 11 cohorts evaluated the effect of SSB intake and 18 genetic variants in 11 genes related to fructose metabolism (ChREB pathway) on glycemic traits in individuals of European descent. They only found a suggestive interaction between KLB (rs1542423) and SSB consumption on fasting insulin, but not other significant interactions on fasting glucose were found [8]. In parallel, a Swedish case–control study including 1253 cases and 1545 controls evaluated the association and interaction effects of rs7903146 and SSB consumption on T2D risk [32]. Although they found a positive association across genotypes, no evidence of interaction effects was found.
Strengths of our study lie in the large number of socio-demographic and lifestyle variables collected by highly trained personnel producing dietary and environmental information with few missing data. The high consumption of SSB (exposure factor) in the Chilean population facilitated assessing its association with the outcome variable. Limitations lie in the low number of genetic loci that were evaluated and in the potential measurement errors of SSB consumption since it was indirectly estimated though a food frequency questionnaire and the relatively limited sample size. Post-hoc power equal showed that our sample size was sufficient to detect associations with an OR ≥ 1.01, depending on allele frequency and genotyping success rate. However, the genetic risk score was constructed with a small proportion of polymorphisms explaining only a small proportion of the genetic risk which may represent a limitation in the estimation of genetic susceptibility.
The fact that we did not identify an association between Ancestry and fasting glucose must not be overinterpreted. This may result from the little ancestry variation present in the MAUCO cohort, reducing its power to detect an association. In addition, Amerindian ancestry was positively and significantly associated with SSB consumption, probably mediated by differences in socio-economic level, indicating that any test for association with ancestry must also consider these other variables to avoid confounding effects. It is important to consider a possible bias and overestimation of the effect size from SSB since it may be associated with other food components of habitual diet as well as detrimental lifestyle factors not included in the survey or models. Furthermore, because the exposure and the outcome were measured at a specific time it is not possible establish a causal relationship.
Despite the described limitations, our study provides the first evidence of genotype-by-diet interaction in a Latin American population that modifies the risk for T2D before clinical sings of the disease. Further studies will be necessary to confirm our findings and comprehensive genomic profiles will be needed to identify genetic variation that may affect T2D risk factors in this population. These types of studies will contribute to devising new hypothesis for potential mechanisms of genetic effects on increasing susceptibility to diabetes and may propose targeted public-health strategies to control the diabetes pandemic in Latin America.

5. Conclusions

This cross-sectional study provides evidence of interaction effects between sugar-sweetened beverages consumption and genetic susceptibility to T2D on fasting glucose in a sample of Chileans who did not report to be diabetic. The adverse effects of SSBs on fasting glucose is observed in those consuming at least two servings per day and is magnified in individuals at the highest genetic susceptibility level in a non-adiposity-mediated manner. The biological pathways triggered by SSB consumption may involve the TCF7L2 and MTNR1B genes.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nu14010069/s1, Supplementary Table S1: SNPs selected from the literature. Supplementary Table S2: analysis of variance for log-glycemia. Supplementary Table S3: Association between selected variables and log-glycemia. Supplementary Table S4: General characteristics of all participants by category of sugar sweetened beverage intake. Supplementary Table S5: Association between SSB intake and fasting glucose (log) in whole sample, men and women. Supplementary Table S6: Association between D2M SNPs and fasting glucose (log) in all participants. Supplementary Table S7: Interaction effects between categorical SSB intake and additive effects from D2M SNPs on fasting glucose. Supplementary Table S8: Interaction between SSB intake categories and genotypes of selected variants on fasting glucose (log). Supplementary Table S9: Association between SNPs (categorical) and fasting glucose (log) (continuous) and pos-hoc power estimation in 2895 subjects. Supplementary Table S10: Sensitivity analysis: interaction between rs7903146 and rs10830963 (as continuous and categorical variable) and SSB categories on log-fasting glucose in subjects younger than 65 years. Supplementary Table S11: Sensitivity analysis: interaction effects between GRSw and sugar-sweetened beverages intake on fasting blood glucose in subjects younger than 65 years. Supplementary Table S12: Sensitivity analysis: Interaction effect between SNPs (rs7903146 and rs10830963) and 2 levels of SSB consumption (continuous and categorical) on log-fasting glucose. Supplementary Table S13: Sensitivity analysis, interaction effect between GRSw (continuous and categorical) and 2 levels of SSB consumption on log-fasting blood glucose.

Author Contributions

Conceptualization, M.L.L.-P., R.A.V. and C.F.; methodology, M.L.L.-P., R.A.V., A.H., E.T.-C., R.R.-O. and C.F.; formal analysis, M.L.L.-P., R.A.V., C.Y., F.P. and V.C.-O.; investigation, A.H., J.A., E.T.-C. and M.L.L.-P.; resources C.F.; data curation, R.R.-O., F.P., V.C.-O., C.Y. and M.L.L.-P.; writing—original draft, M.L.L.-P. and R.A.V.; writing—review and editing A.H., C.F. and R.A.V.; supervision, R.A.V., A.H. and C.F.; project administration, J.A. and M.G.-T.; funding acquisition, R.A.V. and C.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Commission for Scientific and Technological Research (CONICYT) through the funds FONDEF (Grant Number D10E1007), FONDAP (Grant Number 15130011), FONDEQUIP (Grant Number EQM140157), and Doctorate Fellowship 21150938, all from the Government of Chile.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committees of the University of Chile (Protocol Code 055-2016, approved on 23 May 2017) and of the Pontifical University of Chile (Protocol Code 171011004, approved on 30 November 2017).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data available on request due to privacy restrictions.

Acknowledgments

We acknowledged the comments and advice provided by Francisco Pérez, Pilar Durruty, Francisco Cubillos, and Fernando Carrasco as members of the Ph.D. committee for Lourdes López-Portillo. We also acknowledge the board of principal investigators of the Advanced Center for Chronic Diseases (ACCDiS) for reviewing the proposal of this study and facilitating access to the materials and data used.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Malik, V.S.; Hu, F.B. Sugar-Sweetened Beverages and Cardiometabolic Health: An Update of the Evidence. Nutrients 2019, 11, 1840. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Popkin, B.M.; Hawkes, C. Sweetening of the Global Diet, Particularly Beverages: Patterns, Trends, and Policy Responses. Lancet Diabetes Endocrinol. 2016, 4, 174–186. [Google Scholar] [CrossRef] [Green Version]
  3. Taillie, L.S.; Reyes, M.; Colchero, M.A.; Popkin, B.; Corvalán, C. An Evaluation of Chile’s Law of Food Labeling and Advertising on Sugar-Sweetened Beverage Purchases from 2015 to 2017: A before-and-after Study. PLoS Med. 2020, 17, e1003015. [Google Scholar] [CrossRef] [PubMed]
  4. Eyheramendy, S.; Martinez, F.I.; Manevy, F.; Vial, C.; Repetto, G.M. Genetic Structure Characterization of Chileans Reflects Historical Immigration Patterns. Nat. Commun. 2015, 6, 6472. [Google Scholar] [CrossRef]
  5. Verdugo, R.A.; Di Genova, A.; Herrera, L.; Moraga, M.; Acuña, M.; Berríos, S.; Llop, E.; Valenzuela, C.Y.; Bustamante, M.L.; Digman, D.; et al. Development of a Small Panel of SNPs to Infer Ancestry in Chileans That Distinguishes Aymara and Mapuche Components. Biol. Res. 2020, 53, 15. [Google Scholar] [CrossRef] [Green Version]
  6. Barozet, E.; Valenzuela, C.Y.; Cifuentes, L.; Verdugo, R.A.; Herrera, L.; Acuña, M.; Llop, E.; Moraga, M.; Berríos, S.; Di Genova, A.; et al. The Chilean Socio-Ethno-Genomic Cline. Biodemogr. Soc. Biol. 2021, 66, 156–171. [Google Scholar] [CrossRef]
  7. Dhingra, R.; Sullivan, L.; Jacques, P.F.; Wang, T.J.; Fox, C.S.; Meigs, J.B.; D’Agostino, R.B.; Gaziano, J.M.; Vasan, R.S. Soft Drink Consumption and Risk of Developing Cardiometabolic Risk Factors and the Metabolic Syndrome in Middle-Aged Adults in the Community. Circulation 2007, 116, 480–488. [Google Scholar] [CrossRef] [Green Version]
  8. McKeown, N.M.; Dashti, H.S.; Ma, J.; Haslam, D.E.; Kiefte-de Jong, J.C.; Smith, C.E.; Tanaka, T.; Graff, M.; Lemaitre, R.N.; Rybin, D.; et al. Sugar-Sweetened Beverage Intake Associations with Fasting Glucose and Insulin Concentrations Are Not Modified by Selected Genetic Variants in a ChREBP-FGF21 Pathway: A Meta-Analysis. Diabetologia 2018, 61, 317–330. [Google Scholar] [CrossRef] [Green Version]
  9. Ma, J.; Jacques, P.F.; Meigs, J.B.; Fox, C.S.; Rogers, G.T.; Smith, C.E.; Hruby, A.; Saltzman, E.; McKeown, N.M. Sugar-Sweetened Beverage but Not Diet Soda Consumption Is Positively Associated with Progression of Insulin Resistance and Prediabetes. J. Nutr. 2016, 146, 2544–2550. [Google Scholar] [CrossRef] [Green Version]
  10. Kondaki, K.; Grammatikaki, E.; Jiménez-Pavón, D.; De Henauw, S.; González-Gross, M.; Sjöstrom, M.; Gottrand, F.; Molnar, D.; Moreno, L.A.; Kafatos, A.; et al. Daily Sugar-Sweetened Beverage Consumption and Insulin Resistance in European Adolescents: The HELENA (Healthy Lifestyle in Europe by Nutrition in Adolescence) Study. Public Health Nutr. 2013, 16, 479–486. [Google Scholar] [CrossRef] [Green Version]
  11. Lin, W.-T.; Chan, T.-F.; Huang, H.-L.; Lee, C.-Y.; Tsai, S.; Wu, P.-W.; Yang, Y.-C.; Wang, T.-N.; Lee, C.-H. Fructose-Rich Beverage Intake and Central Adiposity, Uric Acid, and Pediatric Insulin Resistance. J. Pediatr. 2016, 171, 90–96. [Google Scholar] [CrossRef] [Green Version]
  12. Imamura, F.; O’Connor, L.; Ye, Z.; Mursu, J.; Hayashino, Y.; Bhupathiraju, S.N.; Forouhi, N.G. Consumption of Sugar Sweetened Beverages, Artificially Sweetened Beverages, and Fruit Juice and Incidence of Type 2 Diabetes: Systematic Review, Meta-Analysis, and Estimation of Population Attributable Fraction. BMJ 2015, 351, h3576. [Google Scholar] [CrossRef] [Green Version]
  13. Stern, D.; Mazariegos, M.; Ortiz-Panozo, E.; Campos, H.; Malik, V.S.; Lajous, M.; López-Ridaura, R. Sugar-Sweetened Soda Consumption Increases Diabetes Risk Among Mexican Women. J. Nutr. 2019, 149, 795–803. [Google Scholar] [CrossRef]
  14. Hu, F.B.; Malik, V.S. Sugar-Sweetened Beverages and Risk of Obesity and Type 2 Diabetes: Epidemiologic Evidence. Physiol. Behav. 2010, 100, 47–54. [Google Scholar] [CrossRef] [Green Version]
  15. Meigs, J.B.; Muller, D.C.; Nathan, D.M.; Blake, D.R.; Andres, R. The Natural History of Progression From Normal Glucose Tolerance to Type 2 Diabetes in the Baltimore Longitudinal Study of Aging. Diabetes 2003, 52, 1475–1484. [Google Scholar] [CrossRef]
  16. Meigs, J.B.; Nathan, D.M.; D’Agostino, R.B.; Wilson, P.W.F. Fasting and Postchallenge Glycemia and Cardiovascular Disease Risk: The Framingham Offspring Study. Diabetes Care 2002, 25, 1845–1850. [Google Scholar] [CrossRef] [Green Version]
  17. DellaValle, D.M.; Roe, L.S.; Rolls, B.J. Does the Consumption of Caloric and Non-Caloric Beverages with a Meal Affect Energy Intake? Appetite 2005, 44, 187–193. [Google Scholar] [CrossRef]
  18. Malik, V.S.; Pan, A.; Willett, W.C.; Hu, F.B. Sugar-Sweetened Beverages and Weight Gain in Children and Adults: A Systematic Review and Meta-Analysis. Am. J. Clin. Nutr. 2013, 98, 1084–1102. [Google Scholar] [CrossRef] [Green Version]
  19. DiMeglio, D.; Mattes, R. Liquid versus Solid Carbohydrate: Effects on Food Intake and Body Weight. Int. J. Obes. 2000, 24, 794–800. [Google Scholar] [CrossRef] [Green Version]
  20. eQTLGen Consortium; Xue, A.; Wu, Y.; Zhu, Z.; Zhang, F.; Kemper, K.E.; Zheng, Z.; Yengo, L.; Lloyd-Jones, L.R.; Sidorenko, J.; et al. Genome-Wide Association Analyses Identify 143 Risk Variants and Putative Regulatory Mechanisms for Type 2 Diabetes. Nat. Commun. 2018, 9, 2941. [Google Scholar] [CrossRef] [Green Version]
  21. Manolio, T.A.; Collins, F.S.; Cox, N.J.; Goldstein, D.B.; Hindorff, L.A.; Hunter, D.J.; McCarthy, M.I.; Ramos, E.M.; Cardon, L.R.; Chakravarti, A.; et al. Finding the Missing Heritability of Complex Diseases. Nature 2009, 461, 747–753. [Google Scholar] [CrossRef] [Green Version]
  22. Flannick, J.; Florez, J.C. Type 2 Diabetes: Genetic Data Sharing to Advance Complex Disease Research. Nat. Rev. Genet. 2016, 17, 535–549. [Google Scholar] [CrossRef]
  23. Fuchsberger, C.; Flannick, J.; Teslovich, T.M.; Mahajan, A.; Agarwala, V.; Gaulton, K.J.; Ma, C.; Fontanillas, P.; Moutsianas, L.; McCarthy, D.J.; et al. The Genetic Architecture of Type 2 Diabetes. Nature 2016, 536, 41–47. [Google Scholar] [CrossRef] [Green Version]
  24. Broad Genomics Platform; DiscovEHR Collaboration; CHARGE; LuCamp; ProDiGY; GoT2D; ESP; SIGMA-T2D; T2D-GENES; AMP-T2D-GENES; et al. Exome Sequencing of 20,791 Cases of Type 2 Diabetes and 24,440 Controls. Nature 2019, 570, 71–76. [Google Scholar] [CrossRef] [Green Version]
  25. Geng, T.; Huang, T. Gene-Environment Interactions and Type 2 Diabetes. Asia Pac. J. Clin. Nutr. 2020, 29, 220–226. [Google Scholar] [CrossRef]
  26. DIAGRAM Consortium; GIANT Consortium; Global BPgen Consortium; Anders Hamsten on behalf of Procardis Consortium; the MAGIC Investigators; Dupuis, J.; Langenberg, C.; Prokopenko, I.; Saxena, R.; Soranzo, N.; et al. New Genetic Loci Implicated in Fasting Glucose Homeostasis and Their Impact on Type 2 Diabetes Risk. Nat. Genet. 2010, 42, 105–116. [Google Scholar] [CrossRef]
  27. Dietrich, S.; Jacobs, S.; Zheng, J.; Meidtner, K.; Schwingshackl, L.; Schulze, M.B. Gene-lifestyle Interaction on Risk of Type 2 Diabetes: A Systematic Review. Obes. Rev. 2019, 20, 1557–1571. [Google Scholar] [CrossRef]
  28. Jia, X.; Xuan, L.; Dai, H.; Zhu, W.; Deng, C.; Wang, T.; Li, M.; Zhao, Z.; Xu, Y.; Lu, J.; et al. Fruit Intake, Genetic Risk and Type 2 Diabetes: A Population-Based Gene–Diet Interaction Analysis. Eur. J. Nutr. 2021, 60, 2769–2779. [Google Scholar] [CrossRef]
  29. Haslam, D.E.; McKeown, N.M.; Herman, M.A.; Lichtenstein, A.H.; Dashti, H.S. Interactions between Genetics and Sugar-Sweetened Beverage Consumption on Health Outcomes: A Review of Gene–Diet Interaction Studies. Front. Endocrinol. 2018, 8, 368. [Google Scholar] [CrossRef] [Green Version]
  30. Ericson, U.; Hindy, G.; Drake, I.; Schulz, C.-A.; Brunkwall, L.; Hellstrand, S.; Almgren, P.; Orho-Melander, M. Dietary and Genetic Risk Scores and Incidence of Type 2 Diabetes. Genes Nutr. 2018, 13, 13. [Google Scholar] [CrossRef] [Green Version]
  31. Wirström, T.; Hilding, A.; Gu, H.F.; Östenson, C.-G.; Björklund, A. Consumption of Whole Grain Reduces Risk of Deteriorating Glucose Tolerance, Including Progression to Prediabetes. Am. J. Clin. Nutr. 2013, 97, 179–187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Löfvenborg, J.E.; Ahlqvist, E.; Alfredsson, L.; Andersson, T.; Dorkhan, M.; Groop, L.; Tuomi, T.; Wolk, A.; Carlsson, S. Genotypes of HLA, TCF7L2, and FTO as Potential Modifiers of the Association between Sweetened Beverage Consumption and Risk of LADA and Type 2 Diabetes. Eur. J. Nutr. 2020, 59, 127–135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ferreccio, C.; Roa, J.C.; Bambs, C.; Vives, A.; Corvalán, A.H.; Cortés, S.; Foerster, C.; Acevedo, J.; Huidobro, A.; Passi, A.; et al. Study Protocol for the Maule Cohort (MAUCO) of Chronic Diseases, Chile 2014–2024. BMC Public Health 2015, 16, 122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Campbell, N.R.; Harmon, S.A.; Narum, S.R. Genotyping-in-Thousands by Sequencing (GT-seq): A Cost Effective SNP Genotyping Method Based on Custom Amplicon Sequencing. Mol. Ecol. Resour. 2015, 15, 855–867. [Google Scholar] [CrossRef]
  35. De la Fuente, C.; Ávila-Arcos, M.C.; Galimany, J.; Carpenter, M.L.; Homburger, J.R.; Blanco, A.; Contreras, P.; Cruz Dávalos, D.; Reyes, O.; San Roman, M.; et al. Genomic Insights into the Origin and Diversification of Late Maritime Hunter-Gatherers from the Chilean Patagonia. Proc. Natl. Acad. Sci. USA 2018, 115, E4006–E4012. [Google Scholar] [CrossRef] [Green Version]
  36. Schwingshackl, L.; Hoffmann, G.; Iqbal, K.; Schwedhelm, C.; Boeing, H. Food Groups and Intermediate Disease Markers: A Systematic Review and Network Meta-Analysis of Randomized Trials. Am. J. Clin. Nutr. 2018, 108, 576–586. [Google Scholar] [CrossRef]
  37. Qi, Q.; Stilp, A.M.; Sofer, T.; Moon, J.-Y.; Hidalgo, B.; Szpiro, A.A.; Wang, T.; Ng, M.C.Y.; Guo, X.; MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium; et al. Genetics of Type 2 Diabetes in U.S. Hispanic/Latino Individuals: Results From the Hispanic Community Health Study/Study of Latinos (HCHS/SOL). Diabetes 2017, 66, 1419–1425. [Google Scholar] [CrossRef] [Green Version]
  38. DIAbetes Genetics Replication and Meta-analysis (DIAGRAM) Consortium; Scott, R.A.; Lagou, V.; Welch, R.P.; Wheeler, E.; Montasser, M.E.; Luan, J.; Mägi, R.; Strawbridge, R.J.; Rehnberg, E.; et al. Large-Scale Association Analyses Identify New Loci Influencing Glycemic Traits and Provide Insight into the Underlying Biological Pathways. Nat. Genet. 2012, 44, 991–1005. [Google Scholar] [CrossRef] [Green Version]
  39. DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium; Asian Genetic Epidemiology Network Type 2 Diabetes (AGEN-T2D) Consortium; South Asian Type 2 Diabetes (SAT2D) Consortium; Mexican American Type 2 Diabetes (MAT2D) Consortium; Type 2 Diabetes Genetic Exploration by Nex-Generation Sequencing in Muylti-Ethnic Samples (T2D-GENES) Consortium; Mahajan, A.; Go, M.J.; Zhang, W.; Below, J.E.; Gaulton, K.J.; et al. Genome-Wide Trans-Ancestry Meta-Analysis Provides Insight into the Genetic Architecture of Type 2 Diabetes Susceptibility. Nat. Genet. 2014, 46, 234–244. [Google Scholar] [CrossRef]
  40. Petermann-Rocha, F.; Lasserre-Laso, N.; Villagrán, M.; Mardones, L.; Martínez, M.A.; Leiva, A.M.; Ulloa, N.; Celis-Morales, C. Asociación Del Polimorfismo Rs7903146, Del Gen TCF7L2, Con Marcadores de Adiposidad y Metabólicos En Población Chilena—Resultados Del Estudio GENADIO. Rev. Méd. Chile 2019, 147, 965–976. [Google Scholar] [CrossRef] [Green Version]
  41. Zhou, Y.; Park, S.-Y.; Su, J.; Bailey, K.; Ottosson-Laakso, E.; Shcherbina, L.; Oskolkov, N.; Zhang, E.; Thevenin, T.; Fadista, J.; et al. TCF7L2 Is a Master Regulator of Insulin Production and Processing. Human Mol. Genet. 2014, 23, 6419–6431. [Google Scholar] [CrossRef] [Green Version]
  42. Jin, T. Current Understanding on Role of the Wnt Signaling Pathway Effector TCF7L2 in Glucose Homeostasis. Endocr. Rev. 2016, 37, 254–277. [Google Scholar] [CrossRef] [Green Version]
  43. Florez, J.C.; Jablonski, K.A.; Bayley, N.; Pollin, T.I.; De Bakker, P.I.W.; Shuldiner, A.R.; Knowler, W.C.; Nathan, D.M.; Altshuler, D. TCF7L2 Polymorphisms and Progression to Diabetes in the Diabetes Prevention Program. N. Engl. J. Med. 2006, 355, 241–250. [Google Scholar] [CrossRef] [Green Version]
  44. Wood, A.R.; Jonsson, A.; Jackson, A.U.; Wang, N.; Van Leewen, N.; Palmer, N.D.; Kobes, S.; Deelen, J.; Boquete-Vilarino, L.; Paananen, J.; et al. A Genome-Wide Association Study of IVGTT-Based Measures of First-Phase Insulin Secretion Refines the Underlying Physiology of Type 2 Diabetes Variants. Diabetes 2017, 66, 2296–2309. [Google Scholar] [CrossRef] [Green Version]
  45. Grimsby, J.L.; Porneala, B.C.; Vassy, J.L.; Yang, Q.; Florez, J.C.; Dupuis, J.; Liu, T.; Yesupriya, A.; Chang, M.-H.; Ned, R.M.; et al. Race-Ethnic Differences in the Association of Genetic Loci with HbA1c Levels and Mortality in U.S. Adults: The Third National Health and Nutrition Examination Survey (NHANES III). BMC Med. Genet. 2012, 13, 30. [Google Scholar] [CrossRef] [Green Version]
  46. Yan, R.; Lai, S.; Yang, Y.; Shi, H.; Cai, Z.; Sorrentino, V.; Du, H.; Chen, H. A Novel Type 2 Diabetes Risk Allele Increases the Promoter Activity of the Muscle-Specific Small Ankyrin 1 Gene. Sci. Rep. 2016, 6, 25105. [Google Scholar] [CrossRef] [Green Version]
  47. Olsen, N.J.; Ängquist, L.; Larsen, S.C.; Linneberg, A.; Skaaby, T.; Husemoen, L.L.N.; Toft, U.; Tjønneland, A.; Halkjær, J.; Hansen, T.; et al. Interactions between Genetic Variants Associated with Adiposity Traits and Soft Drinks in Relation to Longitudinal Changes in Body Weight and Waist Circumference. Am. J. Clin. Nutr. 2016, 104, 816–826. [Google Scholar] [CrossRef]
  48. Brunkwall, L.; Chen, Y.; Hindy, G.; Rukh, G.; Ericson, U.; Barroso, I.; Johansson, I.; Franks, P.W.; Orho-Melander, M.; Renström, F. Sugar-Sweetened Beverage Consumption and Genetic Predisposition to Obesity in 2 Swedish Cohorts. Am. J. Clin. Nutr. 2016, 104, 809–815. [Google Scholar] [CrossRef]
  49. Qi, Q.; Chu, A.Y.; Kang, J.H.; Jensen, M.K.; Curhan, G.C.; Pasquale, L.R.; Ridker, P.M.; Hunter, D.J.; Willett, W.C.; Rimm, E.B.; et al. Sugar-Sweetened Beverages and Genetic Risk of Obesity. N. Engl. J. Med. 2012, 367, 1387–1396. [Google Scholar] [CrossRef] [Green Version]
  50. Qi, L.; Cornelis, M.C.; Zhang, C.; Van Dam, R.M.; Hu, F.B. Genetic Predisposition, Western Dietary Pattern, and the Risk of Type 2 Diabetes in Men. Am. J. Clin. Nutr. 2009, 89, 1453–1458. [Google Scholar] [CrossRef] [Green Version]
  51. Prokopenko, I.; Langenberg, C.; Florez, J.C.; Saxena, R.; Soranzo, N.; Thorleifsson, G.; Loos, R.J.F.; Manning, A.K.; Jackson, A.U.; Aulchenko, Y.; et al. Variants in MTNR1B Influence Fasting Glucose Levels. Nat. Genet. 2009, 41, 77–81. [Google Scholar] [CrossRef]
  52. Dashti, H.S.; Follis, J.L.; Smith, C.E.; Tanaka, T.; Garaulet, M.; Gottlieb, D.J.; Hruby, A.; Jacques, P.F.; Kiefte-de Jong, J.C.; Lamon-Fava, S.; et al. Gene-Environment Interactions of Circadian-Related Genes for Cardiometabolic Traits. Dia Care 2015, 38, 1456–1466. [Google Scholar] [CrossRef] [Green Version]
  53. Lopez-Minguez, J.; Saxena, R.; Bandín, C.; Scheer, F.A.; Garaulet, M. Late Dinner Impairs Glucose Tolerance in MTNR1B Risk Allele Carriers: A Randomized, Cross-over Study. Clin. Nutr. 2018, 37, 1133–1140. [Google Scholar] [CrossRef]
Figure 1. Number of participants and data filtering steps.
Figure 1. Number of participants and data filtering steps.
Nutrients 14 00069 g001
Figure 2. Fasting glucose association with SSB is dependent on genetic risk for T2D. Average log-fasting glycemia is plotted SSB category against tertiles of GRSw (A) and genotypes of rs10830963 (MTNR1B) (B) and rs7903146 (TCF7L2) (C), n = 2828.
Figure 2. Fasting glucose association with SSB is dependent on genetic risk for T2D. Average log-fasting glycemia is plotted SSB category against tertiles of GRSw (A) and genotypes of rs10830963 (MTNR1B) (B) and rs7903146 (TCF7L2) (C), n = 2828.
Nutrients 14 00069 g002
Table 1. General characteristics of participants.
Table 1. General characteristics of participants.
Overall
(n = 2828)
pMen
(n = 971)
pWomen
(n = 1857)
p
Sex (%)- 34.3 65.7
Age (years)53.4 (9.5)0.000154.0 (9.6)0.00353.1 (9.4)0.0001
Education (% low/medium/high)49.0/39.2/11.80.1447.4/40.8/11.80.3449.8/38.4/11.80.33
Socioeconomic (% low/med/high)24.3/58.9/16.80.1316.2/60.8/22.10.7528.0/57.9/14.10.81
Physical Activity (% active/inactive)19.96/80.040.3420.0/80.00.9019.97/80.030.24
Smoking (% never/current/former)42.8/33.1/24.10.000132.9/35.9/31.20.4847.9/31.7/20.40.003
Amerindian ancestry 0.35 (0.08)0.330.34 (0.07)0.190.35 (0.08)0.46
High risk alcohol consumption (%)15.30.3527.60.579.00.37
Fruit intake ≥ 1 portions/day (%)48.30.4538.40.6353.30.03
Vegetable intake ≥ 1 portions/day (%)68.70.3360.00.0773.10.25
Processed meat ≥ 4 portions/week (%)10.10.7212.30.358.60.48
Sugar ≥ 4 teaspoons/day (%)38.80.9447.30.0434.40.0001
SSB category (%) (n)
0 servings/day24.3 (686)0.000113.0 (126)0.1630.2 (560)0.02
>0 and <1 servings/day37.6 (1064) 38.0 (369) 37.4 (695)
≤1 and <2 servings/day25.1 (710) 27.8 (270) 23.7 (440)
≥2 servings/day13.0 (368) 21.2 (206) 8.7 (162)
Glycemia ≥ 100 mg/dL (%) (n)21.9 (618)-30.0 (291)-17.6 (327)-
BMI (kg/m2)29.3 (4.8)0.000128.9 (4.1)0.000129.6 (5.0)0.0001
BMI (normal/overweight/obesity) (%)16.7/44.2/39.10.000116.6/47.0/36.40.00116.7/42.9/40.40.0001
WC (cm)98.1 (10.8)0.0001101.0 (9.4)0.000196.5 (11.0)0.0001
High WC (%)68.9 (1949)0.000145.8 (443)0.000181.0 (1506)0.0001
Triglycerides (mg/dL)162.8 (122.4)0.0001181.6 (153.0)0.0001153.0 (101.5)0.0001
Hypertriglyceridemia (%)44.2 (1251)0.000150.5 (490)0.000141.0 (761)0.0001
HDL-c (mg/dL)45.8 (11.1)0.000142.3 (10.4)0.1947.7 (11.0)0.06
Low HDL-c (%)44.5 (1258)0.2256.2 (544)0.2238.4 (714)0.01
Data are expressed as mean and standard deviation (SD) for continuous variables and as percentage (%) and sample size (n) for categorical variables. p-value for association of each variable with log-fasting glucose are shown. p-values in bold mean significance (<0.05).
Table 2. Interaction between rs7903146 and rs10830963 (as continuous and categorical variable) and levels of SSB categories intake on log-fasting glucose.
Table 2. Interaction between rs7903146 and rs10830963 (as continuous and categorical variable) and levels of SSB categories intake on log-fasting glucose.
rsIDGenGenotypeSSB Category 2SSB Category 3SSB Category 4Pi
Β (SE)PtΒ (SE)PtΒ (SE)Pt
rs7903146TCF7L20,1,20.0004 (0.01) 0.970.020 (0.01) 0.110.05 (0.01)0.0010.002
C/C4.37 (0.053)0.324.20 (0.06)0.444.33 (0.072)0.800.005
C/T−0.006 (0.02)0.100.03 (0.02)0.100.03 (0.02)0.09
T/T0.009 (0.03)0.770.02 (0.03)0.530.14 (0.04)0.0006
rs10830963MTNR1B0,1,20.005 (0.01) 0.670.01 (0.01) 0.430.013 (0.01) 0.470.83
C/C4.21 (0.053)0.264.30 (0.058)0.614.28 (0.072)0.880.001
C/G0.01 (0.02)0.540.01 (0.02)0.41−0.03 (0.02)0.16
G/G−0.004 (0.04)0.910.01 (0.05)0.750.19 (0.05)0.0008
β: intercept of the reference (first) genotype or regression coefficient for number of risk alleles and for non-reference genotypes when treated as a categorical variable, SE: standard error, Pt: p value for trend, Pi: p value for interaction. Model adjusted for: age, sex, BMI, waist circumference, physical activity, education, socio-economic level, smoking, consumptions: fruit, vegetables, sugar, processed meats, alcoholic beverages and % Amerindian ancestry, p values in bold are significant after Bonferroni correction (p < 0.0031), n = 2828.
Table 3. Association between GRSw and fasting blood glucose levels (log) as continuous and categorical variable.
Table 3. Association between GRSw and fasting blood glucose levels (log) as continuous and categorical variable.
GlobalMenWomen
β (SE)pβ (SE)pβ (SE)p
GRSw (continuous)
GRSw0.02 (0.006)0.000020.03 (0.01)0.0050.02 (0.007)0.003
GRSw (categorical)
GRSw tertile 20.001 (0.006)0.85−0.007 (0.01)0.620.004 (0.008)0.60
GRSw tertile 30.02 (0.007)0.00030.02 (0.01)0.080.03 (0.008)0.002
Data are expressed in β (standard error), 95% confidence interval (95% CI). Model adjusted for: age, sex, BMI, waist circumference, physical activity, schooling, socioeconomic level, smoking, consumptions: fruit, vegetables, sugar, processed meats, alcoholic beverages and % Amerindian ancestry, n = 2828. p values in bold are significant (p < 0.05).
Table 4. Interaction effects between GRSw and sugar-sweetened beverages intake on fasting blood glucose levels as continuous/categorical variables.
Table 4. Interaction effects between GRSw and sugar-sweetened beverages intake on fasting blood glucose levels as continuous/categorical variables.
SSB CategoryGlobalMenWomen
β (SE)PtPiβ (SE)PtPiβ (SE)PtPi
SSB and GRSw as continuous
0.02 (0.006)0.0040.0040.03 (0.01)0.010.010.008 (0.006)0.250.25
Categorical SSB and continuous GRSw
20.01 (0.01)0.490.010.003 (0.04)0.920.030.02 (0.02)0.260.46
30.01 (0.02)0.41 0.02 (0.04)0.54 0.007 (0.02)0.71
40.06 (0.02)0.001 0.08 (0.04)0.04 0.03 (0.02)0.16
SSB and GRSw as categorical
GRSw Tertile 2
20.03 (0.02)0.090.020.02 (0.04)0.620.0090.04 (0.02)0.0530.57
30.03 (0.02)0.15 0.06 (0.05)0.22 0.01 (0.02)0.50
40.0005 (0.02)0.98 0.03 (0.05)0.60 0.03 (0.03)0.26
GRSw Tertile 3
20.007 (0.02)0.70 −0.01 (0.04)0.52 0.02 (0.02)0.30
30.02 (0.02)0.41 0.04 (0.05)0.38 0.006 (0.02)0.76
40.05 (0.02)0.02 0.07 (0.05)0.12 0.04 (0.03)0.26
β estimated interaction effect, SE: standard error, Pt: p value for trend, Pi: p value for interaction. Models adjusted for: age, sex, BMI, waist circumference, physical activity, schooling, socio-economic level, smoking, consumptions of fruit, vegetables, sugar, processed meat, and alcoholic beverages, and Amerindian ancestry, n = 2828. p values in bold are significant (p < 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

López-Portillo, M.L.; Huidobro, A.; Tobar-Calfucoy, E.; Yáñez, C.; Retamales-Ortega, R.; Garrido-Tapia, M.; Acevedo, J.; Paredes, F.; Cid-Ossandon, V.; Ferreccio, C.; et al. The Association between Fasting Glucose and Sugar Sweetened Beverages Intake Is Greater in Latin Americans with a High Polygenic Risk Score for Type 2 Diabetes Mellitus. Nutrients 2022, 14, 69. https://doi.org/10.3390/nu14010069

AMA Style

López-Portillo ML, Huidobro A, Tobar-Calfucoy E, Yáñez C, Retamales-Ortega R, Garrido-Tapia M, Acevedo J, Paredes F, Cid-Ossandon V, Ferreccio C, et al. The Association between Fasting Glucose and Sugar Sweetened Beverages Intake Is Greater in Latin Americans with a High Polygenic Risk Score for Type 2 Diabetes Mellitus. Nutrients. 2022; 14(1):69. https://doi.org/10.3390/nu14010069

Chicago/Turabian Style

López-Portillo, María Lourdes, Andrea Huidobro, Eduardo Tobar-Calfucoy, Cristian Yáñez, Rocío Retamales-Ortega, Macarena Garrido-Tapia, Johanna Acevedo, Fabio Paredes, Vicente Cid-Ossandon, Catterina Ferreccio, and et al. 2022. "The Association between Fasting Glucose and Sugar Sweetened Beverages Intake Is Greater in Latin Americans with a High Polygenic Risk Score for Type 2 Diabetes Mellitus" Nutrients 14, no. 1: 69. https://doi.org/10.3390/nu14010069

APA Style

López-Portillo, M. L., Huidobro, A., Tobar-Calfucoy, E., Yáñez, C., Retamales-Ortega, R., Garrido-Tapia, M., Acevedo, J., Paredes, F., Cid-Ossandon, V., Ferreccio, C., & Verdugo, R. A. (2022). The Association between Fasting Glucose and Sugar Sweetened Beverages Intake Is Greater in Latin Americans with a High Polygenic Risk Score for Type 2 Diabetes Mellitus. Nutrients, 14(1), 69. https://doi.org/10.3390/nu14010069

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop