Next Article in Journal
The Role of Nutrition in Cardiovascular Disease: Current Concepts and Trends
Previous Article in Journal
Does Better Diet Quality Offset the Association between Depression and Metabolic Syndrome?
Previous Article in Special Issue
Association of Polygenic Variants Involved in Immunity and Inflammation with Duodenal Ulcer Risk and Their Interaction with Irregular Eating Habits
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association of Plant-Based and High-Protein Diets with a Lower Obesity Risk Defined by Fat Mass in Middle-Aged and Elderly Persons with a High Genetic Risk of Obesity

1
Department of R & D, Daily Manufacturing Inc., Rockwell, NC 28138, USA
2
Department of Food & Nutrition, Obesity/Diabetes Center, Hoseo University, Asan 31499, Republic of Korea
*
Author to whom correspondence should be addressed.
Nutrients 2023, 15(4), 1063; https://doi.org/10.3390/nu15041063
Submission received: 6 December 2022 / Revised: 24 January 2023 / Accepted: 17 February 2023 / Published: 20 February 2023
(This article belongs to the Special Issue Nutrigenomic and Metabolism)

Abstract

:
Obesity has become a severe public health challenge globally. The present study aimed to identify separate and interactive dietary, genetic, and other factors that increase the risk of obesity as measured by body fat (BF) mass. We utilized a genome-wide association study to identify genetic variants associated with high fat mass (obesity; n = 10,502) and combined them to generate polygenic risk scores (PRS) of genetic variants interacting with each other in adults aged over 40 while excluding body-fat-related diseases in a city-hospital-based cohort (n = 53,828). It was validated in Ansan/Ansung plus rural cohorts (n = 13,007). We then evaluated dietary and lifestyle factors in subjects to assess what factors might help overcome a genetic propensity for higher BF. The three-SNP model included brain-derived neurotrophic factor (BDNF)_rs6265, fat-mass- and obesity-associated protein (FTO)_rs1421085, and SEC16B_rs509325. The genes with the minor alleles of ADCY3_rs6545790 and BAIAP2_rs35867081 increased their gene expression in the visceral and subcutaneous adipocytes, but their gene expression decreased in the hypothalamus in eQTL analysis. In the three-SNP model, the PRS was associated with BF mass by 1.408 and 1.396 times after adjusting covariates 1 (age, gender, survey year, residence area, education, and income) and 2 (covariates in model 1 plus energy intake, alcohol intake, regular exercise, and smoking status), respectively. However, when separating subjects by PRS of the three-SNP model, a plant-based diet was the most significant factor associated with low BF, followed by high-protein diets and lower energy intakes. They could offset the effects of high genetic risk for high BF. In conclusion, modulating nutrient intakes might overcome a high genetic risk for obesity. Dietary choices favoring more plant-based and higher-protein foods might help prevent increased BF in Asians and potentially people of other ethnicities with high polygenetic risk scores.

1. Introduction

Obesity is a multi-factored pathological condition involving genetic, environmental, and behavioral factors. Despite intense research efforts over many years, no practical solutions for the prevention and treatment of obesity have been realized. It has been estimated that 40–75% of the variability in body mass index (BMI) is genetic [1]. However, known combinations of single-nucleotide polymorphisms (SNPs) can only predict about 6% of the BMI variability in people of European descent [2]. Although genetic susceptibility to obesity is incompletely described, genome-wide association studies (GWAS) have identified patterns of genetic variations that increase the risk of obesity [3]. Applicability of these genetic factors among all populations is impossible since there are ethnic differences in the risk alleles [4], and many study results may not apply to people of different races. It is also unknown whether genetic predispositions to obesity play a causative or permissive role in most obesity. There have always been examples of people with severe obesity due to inherited traits, but that is uncommon. The recent rapid global increases in obesity rates suggest a more permissive role since the genomes of populations are slow to evolve [5]. It is well known that many people who carry obesity-related genes do not become obese, and for most people, genetics do not equate to destiny [5]. Genetic traits also interact with lifestyles, including dietary patterns. Dietary patterns include a Korean balanced diet (KBD), plant-based diet (PBD), Western-style diet (WSD), and rice-main diet (RMD) in Korea. Their intakes influence obesity differently when interacting with a genetic predisposition [6].
Many studies use BMI as a measure of obesity, but BMI may not represent obesity equally in people, especially in those with a high muscle mass or different fat distributions. It has been demonstrated that people of different ethnicities with similar BMI may exhibit rather significant differences in fat mass [7]. BMI reflects lean body mass and fat mass, but obesity is best defined by fat mass, which is difficult to measure accurately. Furthermore, separating subcutaneous fat from visceral fat is difficult. However, body fat mass can be determined by dual-energy X-ray absorptiometry (DEXA), bioelectrical impedance analysis (BIA), and hydrodensitometry. DEXA and hydrodensitometry can estimate body fat more accurately than BIA [8], and DEXA can predict visceral and subcutaneous fat mass separately. However, DEXA and hydrodensitometry may not be applicable to epidemiological studies [9]. The BIA method is substantially accurate for measuring body fat mass. Therefore, we conducted a GWAS study in a large Korean cohort and evaluated dietary and lifestyle factors that interact with genetic predispositions to obesity defined by body fat mass measured by the BIA method.
The association of fat mass with genetics needs to be studied in large cohorts since numerous dietary and lifestyle factors make small but statistically significant contributions to fat mass. In the Korean Genome and Epidemiology Study (KoGES), fat mass was measured in the Ansan/Ansung cohort (n = 8842) by BIA and estimated in a city-hospital-based cohort study (n = 58,701) by a prediction model generated by a machine learning approach. In an attempt to isolate more significant effects in real-life interventions, we took the analysis a step further and developed polygenic risk scores (PRS) for obesity and separated the population according to those into high and low PRSs and then divided each into those with high and low body fat composition. We hypothesized that investigating people with a high PRS that do not develop obesity would provide information on the most important dietary factors required to maintain a lean body type even in those predisposed to obesity. This study aimed to determine the fundamental approaches to avoiding obesity in the Korean population.

2. Methods

2.1. Participants and Setting

The institutional review board (IRB) of the Korean National Institute of Health and the IRB of Hoseo University approved the KoGES and the present study (KBP-2015-055 and 1041231-150811-HR-034-01, respectively). All participants signed a written informed consent form. In this study, 58,701 participants aged 40–74 years were recruited from a large city-hospital-based cohort that formed part of the Korean Genome and Epidemiology Study (KoGES), which was conducted during 2010–2014 [10]. Adults aged over 40 (n = 13598) that participated in the Ansan/Ansung plus rural cohorts during 2001–2002 were used as the replicate study for exploring fat-mass-related genetic variants. People with a disease history that might influence energy metabolisms, such as cancers, thyroid diseases, chronic kidney disease, and brain-related diseases, were excluded from the study. The 4873 and 591 participants were excluded from the urban-hospital-based cohort and Ansan/Ansung plus rural cohorts, respectively. The participants included in the further study were 53,828 for the city-hospital-based cohort and 13,007 for the Ansan/Ansung plus rural cohorts for the replicate study.

2.2. Demographic, Anthropometric, and Biochemical Parameters of the Participants

The information on age, gender, education, income, alcohol and coffee consumption, physical activity, and smoking history was received during a health interview. Alcohol intake was calculated as the daily amount (g/day) by multiplying the frequency of drinking and the amount at an event [11]. The smoking status was categorized into current smokers (smoked at least 20 cigarettes in the past six months), past smokers (not smoking for at least six months), and non-smokers [11]. The coffee intake was assessed as the weekly drinking frequency and was categorized into three groups by the tertiles of daily coffee intake. Regular physical activity was defined as more than 30 min of moderate physical activity for three or more days per week.
The anthropometric characteristics (height and weight) were measured at the initial visit, as described previously [11]. Body fat and skeletal muscle masses were determined using the Inbody 3.0 (Cheonan, Republic of Korea) based on the BIA method in Ansan/Ansung cohort. The prediction models for body fat and skeletal muscle mass were developed and validated in Ansan/Ansung cohort using a machine learning prediction model [12]. The prediction models generated from the XGBoost and deep neural network (DNN) algorithms were applied to the city-hospital-based cohort, and body fat and skeletal muscle masses were estimated in the participants of a city-hospital-based cohort. The prediction models by XGBoost and DNN showed good accuracy for body fat mass and skeletal muscle (R2 = 0.82 and 0.89, respectively) [12]. The predicted values of body fat mass were used for the present study.
A doctor measured blood pressure three times via a sphygmomanometer under resting conditions, and the average systolic blood pressure (SBP) and diastolic blood pressure (DBP) were reported. After fasting for more than 12 h, blood was collected in heparin-treated and non-treated tubes, and the plasma glucose concentrations were measured using a Hitachi 7600 Automatic Analyzer (Hitachi, Tokyo, Japan). HbA1c in heparin-treated blood was measured using an automatic analyzer (ZEUS 9.9; Takeda, Tokyo, Japan). Serum total cholesterol, high-density lipoprotein cholesterol (HDL-C), triglyceride and creatinine concentrations, and alanine aminotransferase (ALT) and aspartate aminotransferase (AST) activities were assessed using a Hitachi 7600 Automatic Analyzer. The serum high-sensitive C-reactive protein (hs-CRP) concentrations were measured using a high-sensitivity ELISA kit (Thermofisher, Waltham, MA, USA). The white blood cells (WBCs) were counted from the EDTA-treated blood. The estimated glomerular filtration rate (eGFR) was calculated using the equation by Modification of the Diet in Renal Disease (MDRD) study [13].

2.3. Definition of Obesity and Metabolic Syndrome

Obesity was defined according to total body fat mass estimated by BIA: men with ≥25% (n = 3298) and women with ≥30% (n = 7204) were considered to be obese and assigned to the high-body-fat (high-BF) group (the case group) by Korean definitions of obesity. Metabolic syndrome is a cluster of energy, glucose, and lipid disorders categorized according to the 2005 revised National Cholesterol Education Program—Adult Treatment Panel III criteria for Asia [14,15]. The criteria for MetS were as follows: (1) abdominal obesity (waist circumference ≥ 90 cm for men and ≥85 cm for women), (2) elevated fasting blood glucose level (≥100 mmol/L) or current use of anti-diabetic medication, (3) elevated blood pressure (average systolic blood pressure ≥ 130 mmHg or diastolic blood pressure ≥ 85 mmHg) or current blood pressure medication use, (4) low HDL-C level (<40 mg/dL for men and <50 mg/dL for women), or (5) elevated serum triglyceride level (≥150 mmol/L) or current use of anti-dyslipidemic medication. Participants meeting three or more criteria were considered to have MetS.

2.4. Usual Food Intake Using a Semi-Quantitative Food Frequency Questionnaire (SQFFQ)

During the last 12 months of the interview, the food intake of each participant was determined using an SQFFQ designed for the Korean diet. The SQFFQ accuracy and reproducibility were validated with the three-day food records in four seasons in Koreans [16,17]. The SQFFQ includes 106 food items that Koreans commonly consume, and the food frequencies were categorized into never or seldom, once per month, two to three times monthly, once or twice weekly, three or four times weekly, five or six times weekly, daily, twice daily, and ≥3 times daily. The food amount at a meal was scored as more than, equal to, or less than the regular portion size visualized by photographs of 106 foods. The participants checked the frequencies and the portion size of 106 food items in the SQFFQ. The daily food intake was calculated by multiplying the median of the weekly consumed frequencies by portion sizes. The food intake was given in grams/day. The daily energy, carbohydrates, fat, protein, vitamin, and mineral intakes were calculated from the SQFFQ results using Can-Pro 2.0 nutrient assessment software designed by the Korean Nutrition Society.

2.5. Dietary Patterns by Principal Components Analysis

For dietary pattern analysis, 106 food items in the SQFFQ were classified into 30 predefined food groups, as reported previously [18]. The dietary patterns were designed using principal component analysis of the 30 food groups based on eigenvalues > 1.5, and four dietary patterns explained the criteria [18]. The orthogonal rotation procedure (varimax) yielded four dietary patterns uncorrelated with each other, and foods with ≥0.40 factor-loading values were considered the predominant contributors to the assigned dietary pattern [18]. Table 1A lists the predominant foods found as >0.4 of factor-loading values in each dietary pattern. These patterns indicate the participant’s diet types, and they were divided into the KBD, PBD, WSD, or RMD. Nutrient intake in each dietary pattern is shown in Table 1B. Energy intake with high carbohydrates and low fat was lower in RMD, but energy intake with high fat and low carbohydrate was higher in the WSD group.

2.6. Dietary Inflammatory Index (DII)

DII, an index of the pro-inflammatory potential of diets, was calculated from the equation with assigned food and nutrient intakes using their dietary inflammatory weights including energy, 32 nutrients, four food products, four spices, and caffeine, as described previously [19]. Because the SQFFQ did not include garlic, ginger, saffron, and turmeric, their intakes were excluded from the DII calculation. DII was calculated by multiplying the dietary inflammatory scores of the 38 food and nutrient components by the daily intakes, and the sums of the scores of 38 items were divided by 100 [19,20].

2.7. Genotyping Using a Korean CHIP and Quality Control

The Center for Genome Science at the Korea National Institute of Health determined the participants’ genotypes in the Ansan/Ansung and city-hospital-based cohorts. The genomic DNA was isolated from whole blood, and genotypes were measured using a Korean Chip (Affymetrix, Santa Clara, CA, USA) designed to examine the disease-related single-nucleotide polymorphisms (SNPs) in Koreans [10]. The genotyping accuracy was estimated using Bayesian robust linear modeling in the Mahalanobis distance genotyping algorithm. The inclusion criteria of the genotyping accuracy, missing genotype call rate, and heterozygosity were ≥98%, <4%, and <30%, respectively, and the data showed no gender bias. The genetic variants that were included satisfied the Hardy–Weinberg equilibrium (HWE) at p > 0.05 and minor allele frequency (MAF) at >1% [21]. Manhattan and quantile–quantile (Q-Q) plots indicated the accuracy of GWAS data using the Fastman library in the R program [21]. A Manhattan plot of genetic variants was displayed with the negative logarithms of the association of p-values for high body fat mass. A Q-Q plot is a probability plot to show the goodness of fit of the actual data distribution to the theoretical data distribution. The Q-Q plot of genotype data displayed the quantile distribution of observed p-values (on the y-axis) versus the quantile distribution of expected p-values (on the x-axis). The Q-Q plot was constructed to ensure that the lambda value of the Q-Q plot was close to 1 and confirmed that the GWAS genotypes were ideal. The pathways linked to the genetic variants associated with high body fat mass having p-values < 0.05 for Bonferroni correction were selected using the MAGMA gene-set analysis in SNP2GENE of FUMA web application, available through the git repository (https://github.com/Kyoko-wtnb/FUMA-webapp/, accessed on 18 May 2022).

2.8. Selection of the Genetic Variants That Influence Obesity Defined by Fat Mass and the Best Model with SNP–SNP Interactions

Figure 1 presents the procedure of selecting genetic variants for high body fat risk and investigating the best model for SNP–SNP interactions. The GWAS was conducted to explore genetic variants associated with obesity risk in the urban-hospital-based cohort (p < 5 × 10−5). From the GWAS associated with obesity risk, 1992 genetic variants were selected at p < 5 × 10−5. We eliminated 387 genetic variants that did not meet MAF (<1%) and HWE (p < 0.05). In the gene name search using g:Profiler (https://biit.cs.ut.ee/gprofiler/snpense, accessed on 2 June 2022), 178 SNPs were not identified with gene names, and 126 gene names of the 1427 genetic variants were identified. The linkage disequilibrium (LD) analyses were performed on the SNPs of the 1427 genetic variants using Haploview 4.2 in PLINK. The potential genetic variants in the same chromosome were not strongly correlated (D’ < 0.2). The SNPs with high D’ values were not included in the generalized multifactor dimensionality reduction (GMDR) because they provided the same information on the genetic impact. There were 19 SNPs in 18 genes selected, and genes associated with fat mass were selected using HuGE Navigator (https://phgkb.cdc.gov/PHGKB/hNHome.action, accessed on 20 June 2022).
Among the 14 genetic variants selected for obesity risk, ten SNPs with an SNP–SNP interaction were selected automatically by GMDR. The best SNP–SNP interaction model was selected in a sign-rank test of trained balanced accuracy (TRBA) and testing balanced accuracy (TEBA) while adjusting for the covariates using a GMDR program and a p-value threshold of 0.05 [10]. The covariates used were age, gender, survey year, residence area, education, and income for models 1 and 2 plus energy intake, alcohol intake, regular exercise, and smoking status when carried out GMDR and logistic regression between BF groups and genetic variants. Models 1 and 2 indicated the results with different covariates. Ten-fold cross-validation was also checked for cross-validation consistency (CVC) because the sample size was larger than 1000 [22]. In total, 10 out of 10 in the CVC met the perfect cross-validation criteria. The association between body fat mass and ten SNPs selected from a city-hospital-based cohort was validated in the Ansan/Ansung plus rural cohorts (n = 13,007). Their significance level was considered as 0.05 since the validation was conducted in the selected genetic variants.
The risk allele number of each SNP was counted to generate the PRS of the best model. For example, the genetic score for the SNP was 2, 1, and 0 when the participants had AA, AG, and GG of one SNP, and the A allele was the risk allele, respectively. The polygenic risk score (PRS) of the best model was assessed by summing the number of the risk alleles from each selected SNP in the best gene–gene interaction model [3,23]. The criteria of the best mode were the p-value of the sign test for TEBA, which was 0.001, and CVC, which was 10. The two models with the smallest genetic variants were selected among the models meeting the criteria. From the GMDR analysis, the PRS in the three- and six-SNP models were divided into three categories according to the number of risk alleles. They were classified as low PRS, middle PRS, and high PRS when the number of risk alleles in the PRS was 0–2 (n = 19,686), 3–4 (n = 30,513), and ≥5 (n = 3629) in the three-SNP model and 0–5 (n = 27,212), 6–7 (n = 20375), and ≥8 (n = 1822) in the six-SNP model, respectively. Among the best models to meet the p-value of the sign test and CVC, the model with the lowest SNP number (three-SNP model) was used to interact with the lifestyle parameters.

2.9. Expression Quantitative Trait Locus (eQTL) Analysis

The eQTL analysis is a direct approach to estimating the candidate gene expression of the genetic variants at risk loci. The allele variants are involved in the corresponding gene expression, and the expression of candidate susceptible genes with risk alleles is estimated to influence various diseases [23]. Gene expressions corresponding to the genetic variants related to the abdominal obesity risk were identified by eQTL analysis in the Genotype-Tissue Expression (GTE) × eQTL calculator (https://gtexportal.org/home/testyourown, accessed on 19 July 2021). The gene expressions in the subcutaneous and visceral adipose tissues, skeletal muscles, liver, and brain were calculated using the GTE × eQTL calculator.

2.10. Statistical Analysis

The statistical analysis was performed using SAS (version 9.3; SAS Institute, Cary, NC, USA). A sample size of 53,828 was sufficient to achieve significance at α = 0.05 and β = 0.99 at an odds ratio of 1.05 in the logistic analysis using a G-power calculator. The descriptive statistics for categorical variables, such as gender and dietary habits, were obtained by determining the frequency distributions, which were analyzed statistically according to the low-BF and high-BF groups of the classification variables using a chi-square test. Descriptive statistics of the continuous variables were analyzed as the adjusted means with standard deviations after adjusting for the covariates. The statistical differences among the gender and insulin-resistance groups were compared using a two-way analysis of covariance (ANCOVA) [24]. Multiple comparisons of the groups were performed using Tukey’s test.
The association of insulin resistance on metabolic parameters was examined by logistic regression analysis with the low-BF group as the reference after adjustment for covariates. The results are presented as the odds ratios (ORs) and 95% confidence intervals (CI) of each biochemical parameter for the high-BF and low-BF groups. The first model was generated after adjusting for age, residence area, survey year, lean body mass, education, and income. The second model was produced with the adjustments for covariates in model 1 and the energy intake, physical activity, smoking status, and alcohol consumption.
The lifestyle-related parameters were categorized into the high or low groups using the criteria defined by the predesignated cutoffs, such as the dietary reference intake or 30th percentiles of each variable, to determine the interactions between the fat mass and lifestyle parameters. Two-way ANCOVA was used to analyze the interactions between the fat mass groups and the lifestyle parameters, including dietary intake, smoking, and physical activity. The main effects were insulin resistance and lifestyle-related parameters with their interaction terms after adjusting for covariates. The ORs and 95% CI of fat mass with lifestyle-related parameters were also calculated by logistic regression analysis in the high-BF and low-BF groups of the lifestyle-related parameters. The significant difference in the high-BF percentage was analyzed using the PRS groups in the χ2 test in the low and high groups of lifestyle-related parameters.

3. Results

3.1. Demographic Characteristics and Lifestyles According to Genders and Obesity

Age was significantly higher in the high-BF group than in the low-BF group for both genders. Education and income also significantly affected body fat mass for both genders: The participants with less than a high school education had a significantly higher incidence of high body fat for both genders. Men in the lowest income group had a significantly lower incidence of high body fat, whereas women in the lowest income group had a significantly higher incidence of high body fat mass (Table 2). Energy intake was not significantly different between the high- and low-BF groups in women; however, men in the low-BF group had slightly but significantly higher energy intakes than those in the high-BF group. Carbohydrate and fat proportions were not significantly different according to body fat. In dietary fat composition, MUFA intake was higher, and PUFA intake was lower in the high-BF group than in the low-BF group. Vitamin C, fiber, and Ca intakes were lower in the high-BF group than in the low-BF group only in women, but vitamin D intake showed the same trend in both genders (Table 2). DII was the same between fat mass groups in both genders although DII was higher in men than women. Among dietary patterns, the number of men consuming a KBD was higher in the high-fat mass group, but it was the opposite in women. The number of participants with high PBD was lower in the high-BF group only in women, but those with a high WSD were higher in the high-BF group only in men (Table 2). The participants who were smokers were lower in number in the high-BF group than the low-BF group only in men. Alcohol intake did not differ between low- and high-BF groups for either gender. However, the number of participants who participated in regular exercise was lower in the high-BF group than in the low-BF group (Table 2).

3.2. Anthropometric and Biochemical Parameters According to Genders and Obesity

The participants with high body fat mass were shorter and had higher BMIs and waist circumferences but lower SMI than those with low body fat mass (Table 3). Height and SMI were inversely associated with body fat mass, and BMI and waist circumferences were positively linked to body fat mass. MetS was higher in the high-BF group than in the low-BF group and positively associated with body fat mass by 6.29 times (Table 3). Serum glucose and blood HbA1c concentrations were significantly lower in men in the low-BF group than in men in the high-BF group. Serum glucose but not HbA1c was significantly lower in women in the low-BF group as compared to the high-BF group (Table 3). Furthermore, dyslipidemia, including hypercholesterolemia, hyper-LDL-cholesterolemia, hypo-HDL-cholesterolemia, and hypertriglyceridemia, was also positively linked to fat mass (Table 3). However, hypertension was less associated with fat mass. eGFR was lower in high-BF than low-BF in men but not women. Serum AST and ALT concentrations, known as liver damage indexes, were higher in the high-BF group than in the low-BF group, whereas they were positively associated with fat mass (Table 3). Serum hs-CRP concentrations were higher in the high-BF group than in the low-BF group and were positively related to fat mass (Table 3).

3.3. Polygenetic Variants and Their Interactions Related to Obesity Defined by Fat Mass

The statistical association of genetic variants with body fat mass is shown in a Manhattan plot (Figure 2A). The Q-Q plot (Figure 2B) shows the quantile distribution of the log of observed p-values versus that of the log of expected p-values. The lambda value was 1.116. These results suggested that genetic variants for fat mass were in an acceptable range although they might be somewhat inflated.
The gene-set analyses used the full distribution of SNP p-values, and we found which GO pathways were related to the prioritized genes. The genes of selected genetic variants found in GWAS between low-BF and high-BF groups were clustered with the GO pathways. Table 4 shows the GO pathways related to the genetic variants linked to body fat mass (p < 0.0002). Genes of the genetic variants associated with fat mass risk were highly involved in the sodium ion transmembrane transport pathway in the biological process of gene ontology (GO) (Table 4). Interestingly, the genes were associated with myoblast proliferation and the apoptotic process of neurons in the biological process of GO. In the curated gene sets, the biocarta flumazenil pathway and Reactome signaling to p38 via RIT and RIN were associated with genes related to fat mass (Table 4). These pathways were related to toxic compound metabolism by MAPK and RAS.
Ten genetic variants with the genetic variant–genetic variant interaction are presented in Table 5. The ten genetic variants were rs509325_SEC16B, rs6545790_ADCY3, rs7560575_PSME4, rs2196476_SLIT2, rs6265_BDNF, rs587056_FARP1, rs1421085_FTO, rs35867081_BAIAP2, rs60259426_SYMPK, and rs6089240_CDH4. The genes are known to be involved in obesity. They met the MAF > 1% and p-value for HWE > 0.05, whereas the adjusted ORs were 0.9278–1.35, and their p-values were <5 × 10−6 and <0.0069 in the city-hospital-based cohort and Ansan/Ansung cohort plus rural cohort, respectively.

3.4. The Gene Expression According to the Alleles of the Selected SNPs in Different Tissues from GTEx v8

The expression of ADCY3 with the minor allele of rs6545790 was elevated with 0.21 of slope in visceral fat (p = 4.1 × 10−14) and 0.23 in subcutaneous fat (p = 2.2 × 10−19) (Figure 3). It indicated that ADCY3 having a risk allele had an increased expression compared to that having a non-risk allele. By contrast, it significantly decreased with a −0.18 and −0.28 slope in the amygdala (p = 0.0064) and brain cortex (p = 5.4 × 10−10), respectively. The SEC16B expression with the minor allele of rs509325 increased with a slope of 0.24 in the adrenal gland. The BAIAP2 expression with the minor allele of rs35867081 was elevated with a 0.066 (p = 0.025) and −0.14 (p = 0.000035) slope in the subcutaneous adipose tissue and the hippocampus, respectively (Figure 3). However, the minor allele of BDNF_rs6265 was not significantly associated with gene expression (p = 0.55) although it was a missense mutation.

3.5. The Best Model of Genetic Variants with SNP–SNP Interaction for Obesity

The best models with genetic variant–genetic variant interactions influencing obesity determined by fat mass were selected when satisfying p-value < 0.05 for the sign test of TEBA and CVC 10/10. The models to meet the criteria included 3-, 6-, 7-, 8-, and 10-SNPs (Supplementary Table S1). The three-SNP model included BDNF_rs6265, FTO_rs1421085, and SEC16B_rs509325. In the three-SNP model, the PRS was associated with body fat mass by 1.408 and 1.396 times after adjusting for covariates 1 and 2, respectively (Figure 4). The six-SNP model contained ADCY3_rs6545790, BDNF_rs6265, SEC16B_rs509325, BAIAP2_rs35867081, SYMPK_rs60259426, and CDH4_rs6089240. The PRS of the six-SNP model was related to body fat mass by 1.36 and 1.37 times, respectively. However, the adjusted ORs in the six-SNP model were decreased compared to the three-SNP model (Figure 4). Therefore, the three-SNP model was used as the best model for body fat mass risk.

3.6. Interaction of PRS with Lifestyles to Influence Obesity

Energy intake interacted with PRS of the three-SNP model to affect body fat mass, and positive associations of PRS with body fat mass risk were 1.228 and 1.157 times, respectively, in low and high energy intakes (Table 6). The proportion of the participants was higher in the high-PRS group than in the low-PRS group in both low and high energy intake (Figure 5A). The results suggested that the participants with high PRS increased fat mass even in low energy intake. Protein intake also interacted with PRS to influence fat mass (Table 6). A high protein intake could prevent the increase of fat mass in the participants in the high-PRS group compared to those in the low-PRS group (Figure 5B). Other nutrient intakes did not interact with PRS affecting fat mass. Among four dietary patterns, only PBD intake interacted with PRS for fat mass (Table 6). The high body fat mass proportion was higher in the participants with high PRS than those with low PRS in both low and high PBD intake. However, the increase was much more significant in the low PBD intake than in the high PRB intake (Figure 5C). Exercise and smoking status did not interact with PRS to affect fat mass. Therefore, these results suggest that Koreans with high PRS can decrease their risk of high body fat by consuming a PBD diet with less energy and high protein intake.

4. Discussion

The present study investigated genetic, lifestyle, and dietary factors that differentiate people with healthy fat mass levels from those with excessive body fat. The study utilized a genome-wide association study (GWAS) followed by an analysis of PRS for high body fat mass (obesity) and comparisons of body fat mass status in subjects according to high and low risk scores. This approach was used to highlight what dietary and lifestyle patterns might lead to an individual being lean despite a high genetic risk of obesity and vice versa.
This study also looked at four dietary patterns: KBD, PBD, WSD, and RMD. The diets varied in energy and nutrient content, varying from 179% of the estimated energy requirements for the WSD to 87% for the RMD. Nevertheless, there was no clear distinction between the dietary patterns except that only the PBD was clearly associated with a lower fat mass (110% of the estimated energy requirements). This may not be too surprising. In a recent perspective article, Mozaffarian published a graph that clearly showed a large increase in the percentage of obesity in the USA population from the year 2000 (30%) to 2018 (42%) [25]. Despite the rapid increase in the obesity rate, there was no change in energy intake. It highlights the need to look for causes of obesity beyond energy intake. This study attempts to contribute to understanding the dietary causes of obesity while realizing that there are unlikely to be straightforward explanations.
The three genetic variants constituted the best model for high body fat mass contained BDNF, FTO, and SEC16B in the present study. Genetic variants that result in decreased production of BDNF have been shown to increase the obesity risk significantly in humans [26]. The BDNF_rs6265 variant causes a methionine substitution for valine, which results in multiple physiological and psychological changes, including impaired fat oxidation in muscle, type 2 diabetes, depression, and disruption of hunger and satiation signals [24,27]. Interestingly, the consequences of the BDNF_rs6265 are associated with obesity risk, and it interacts with lifestyles, including nutrient intake [28]. Furthermore, it interacts with other obesity-related genetic variants to exacerbate obesity risk, and PRS, including BDNF_rs6265, interacts with lifestyles to modulate the obesity risk [3]. These results suggest that neuronal circuit integrity, survival, and synaptic plasticity may be associated with obesity risk.
The FTO gene, located at chromosome 16q12.2, encodes for the fat-mass- and obesity-associated protein. The FTO protein functions as a nucleic acid demethylase, and thus, it is involved in controlling the methylation of DNA and RNA. However, the FTO_rs1421085 non-coding variant does not affect the expression or the functionality of the FTO protein. Mutated FTO (rs1421085) interacts with Iroquois-class homeodomain proteins IRX-3 and IRX-5, which are also linked to obesity in humans and animals [28,29]. Laber et al. [30] have demonstrated that the FTO_rs 1421085 interaction with IRS-3 and IRX-5 modulates numerous organismal phenotypes related to diet-induced weight gain and observes that the FTO_rs1421085 genetic variant mediates cross-species conserved changes in steroid patterns following nutritional challenges.
The third member of the best model group, SEC16B_rs509325, is a variant of the SEC16B gene, also known as RGPR-p117, a gene that encodes a protein known as regucalcin that is characterized by its leucine zipper motif. Regucalcin is known to be a transporter of proteins from the endoplasmic reticulum, but evidence to date has not demonstrated the role of the transporter in lipid metabolism [31]. Nevertheless, variants of the SEC16B gene have consistently been shown to be linked to human obesity, especially in Asian populations [3,32]. Although there has not been any link to fat metabolism established for SEC16B, it is linked to transport peptides involved in regulating satiety, including neuropeptide Y and proopiomelanocortin [3,33]. All three genes composing the best model affect different aspects of metabolism.
One of the objectives of this study was to determine which dietary patterns and/or dietary factors might help overcome a genetic propensity to obesity. Three factors interacted with the PRS to lower obesity risk. A plant-based diet was the most highly significant factor (p = 0.0026) in ameliorating the risk of obesity in participants with a high genetic risk for obesity. KBD, WSD, and RBD had no statistically significant interaction with a high PRS. Energy intake and protein intake also significantly interacted with a high PRS, and a high-protein diet seemed to mostly eradicate the effect of a high PRS on fat mass.
The most important finding of the current study was that the primary dietary factors linked to low body fat content other than lower energy intake were high dietary protein consumption and, even more potent, a more plant-based diet. Among the four dietary patterns studied, most had no significant effect on the risk for high body fat content; however, diets considered low in plant-based foods and low in protein resulted in a significantly higher percentage of obese subjects among those with medium and high PRS. This result seems counterintuitive since plant-based diets are typically lower in protein than diets high in animal foods [34]. However, these results concur with studies in other ethnic groups that have found that vegetarian and/or plant-based diets support a leaner body type [35,36].
There are various hypotheses about how a plant-based diet may help maintain a healthy body weight. First, a plant-based diet has a lower energy density, resulting in greater satiety and lower energy consumption [22,36]. Furthermore, plant-based diets contain lower fat and higher fiber, which may impact the microbiome, which is known to have substantive impacts on energy utilization in humans [37,38,39]. However, microbiome research has not progressed to the point that changes in the gut microbiome can be considered a cause or cure for excess body weight. There is also evidence that various phytochemicals from plant-based foods and dietary supplements may have anti-obesity effects by various mechanisms [23]. It has been shown that the anti-inflammatory properties of plant-based diets can prevent obesity and even result in a healthier obese phenotype in people who become obese [40]. The present study revealed a weak inverse association between high body fat mass and dietary inflammation index (OR = 0.927; 95% CI = 0884–0.973) despite a high inverse association between PBD and DII (OR = 0.368; 95% CI = 0.348–0.339). Even though it has been difficult to demonstrate specific phytochemicals that profoundly affect body fat mass, it has been previously demonstrated that a phytochemical index score was associated with a lower prevalence of obesity in the Korean population [41]. Therefore, the accumulation of current evidence supports the finding in this study that consuming a plant-based diet provides significant benefits for managing body fat.
Previous studies have also found an association between higher protein intake and lower body fat mass. It has also been observed in a recent meta-analysis and systematic review that diets high in protein (18–59 energy %) were associated with an average of 1.6 kg decrease in body weight [42]. The present study did not determine a more precise dietary protein level needed for weight reduction. However, the efficacy of increased protein intake for weight loss has support from a human trial with biscuits enriched with plant-based protein. The cutoff of protein intake was 15 En%, so our results suggest that higher than 15 En% from protein may decrease fat mass. The previous study found that the group with the protein-added biscuits lost significantly more body weight and fat mass than the control group [43].
In this study, having a plant-based diet had the most significant association with lower body fat, followed by a high-protein diet. It suggests that a plant-based diet high in protein might be optimal for the Korean population and possibly people of other ethnic groups. The findings support the possibility that substituting plant protein for animal protein can decrease the incidence of type 2 diabetes, probably by decreasing obesity-induced inflammation [44]. It was also demonstrated that elderly Korean men with concurrent sarcopenia, osteoporosis, and obesity had the same total and animal protein intake but a significantly lower intake of plant protein than controls [45].
This research revealed that modest differences in dietary choices could have a significant impact on fat accumulation in the Korean population. Other studies have shown that people who adhere to vegetarian diets are less obese than their omnivorous counterparts [46] and that vegetarians are healthier even at the same BMI [47]. However, the subjects in the present study with a plant-based diet were not necessarily vegetarians. A plant-based diet in this study indicated that plant foods contributed relatively more to the total energy intake than animal foods. However, the subjects did not adhere to any particular dietary restrictions. The lack of a clearly defined diet could be both a strength and a limitation of the study. There was more variability in the dietary type and potentially less of an effect than if all subjects with a plant-based diet were vegetarians. On the other hand, this represents a more “real world” scenario for many people and shows that modest changes can have profound effects.

5. Conclusions

In conclusion, this study suggested that BDNF_rs6265, FTO_rs1421085, and SEC16B_rs509325 interact with high body fat mass, that plant-based and high-protein diets are beneficial for maintaining a lower body fat in adults with a high PRS of the three genetic variants. The ADCY3_rs6545790 and BAIAP2_rs35867081 minor alleles included in the three-SNP model increased their gene expression in the visceral and subcutaneous adipocytes. However, their gene expression decreased in the hypothalamus in eQTL analysis, suggesting the changes in their expression could be involved in increased body fat mass. This study also demonstrated that the discovery of PRS–environment interactions in body fat was applicable in the clinical setting. Therapies targeting the PRS could help prevent obesity and manage weight by applying personalized nutrition programs emphasizing a plant-based diet high in protein. A large clinical randomized trial is needed to show the cause-and-effect relationship between the three SNP combinations and PBD in the development of obesity. If confirmed, it can be applied to personalized nutrition.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu15041063/s1, Table S1: The characteristics of the ten genetic variants of genes related to fat mass in adults.

Author Contributions

S.P. conceived, designed, and analyzed the data; J.W.D. wrote the draft of the article. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Research Foundation of Korea.

Institutional Review Board Statement

The institutional review board (IRB) of the Korean National Institute of Health and the IRB of Hoseo University approved the KoGES and the present study (KBP-2015-055 and 1041231-150811-HR-034-01, respectively). All participants signed a written informed consent form.

Informed Consent Statement

All participants signed a written informed consent form.

Data Availability Statement

The data can be retrieved from the corresponding author upon reasonable request.

Acknowledgments

We appreciate the participants in the KoGES to provide time and efforts.

Conflicts of Interest

The authors have no potential conflict of interest.

References

  1. Jacob, R.; Bertrand, C.; Llewellyn, C.; Couture, C.; Labonté, M.; Tremblay, A.; Bouchard, C.; Drapeau, V.; Pérusse, L. Dietary Mediators of the Genetic Susceptibility to Obesity-Results from the Quebec Family Study. J. Nutr. 2022, 152, 49–58. [Google Scholar] [CrossRef] [PubMed]
  2. Yengo, L.; Sidorenko, J.; Kemper, K.E.; Zheng, Z.; Wood, A.R.; Weedon, M.N.; Frayling, T.M.; Hirschhorn, J.; Yang, J.; Visscher, P.M. Meta-analysis of genome-wide association studies for height and body mass index in ~700,000 individuals of European ancestry. Hum. Mol. Genet. 2018, 27, 3641–3649. [Google Scholar] [CrossRef] [PubMed]
  3. Park, S.; Yang, H.J.; Kim, M.J.; Hur, H.J.; Kim, S.H.; Kim, M.S. Interactions between Polygenic Risk Scores, Dietary Pattern, and Menarche Age with the Obesity Risk in a Large Hospital-Based Cohort. Nutrients 2021, 13, 3772. [Google Scholar] [CrossRef] [PubMed]
  4. Konuma, T.; Okada, Y. Statistical genetics and polygenic risk score for precision medicine. Inflamm. Regen. 2021, 41, 18. [Google Scholar] [CrossRef]
  5. Hu, F.B. 437Genetic Predictors of Obesity. In Obesity Epidemiology; Hu, F.B., Ed.; Oxford University Press: Oxford, UK, 2008. [Google Scholar]
  6. Lim, H.; Kim, S.Y.; Wang, Y.; Lee, S.J.; Oh, K.; Sohn, C.Y.; Moon, Y.M.; Jee, S.H. Preservation of a traditional Korean dietary pattern and emergence of a fruit and dairy dietary pattern among adults in South Korea: Secular transitions in dietary patterns of a prospective study from 1998 to 2010. Nutr. Res. 2014, 34, 760–770. [Google Scholar] [CrossRef]
  7. Heymsfield, S.B.; Peterson, C.M.; Thomas, D.M.; Heo, M.; Schuna, J.M., Jr. Why are there race/ethnic differences in adult body mass index-adiposity relationships? A quantitative critical review. Obes. Rev. 2016, 17, 262–275. [Google Scholar] [CrossRef] [Green Version]
  8. Hamilton-James, K.; Collet, T.-H.; Pichard, C.; Genton, L.; Dupertuis, Y.M. Precision and accuracy of bioelectrical impedance analysis devices in supine versus standing position with or without retractable handle in Caucasian subjects. Clin. Nutr. ESPEN 2021, 45, 267–274. [Google Scholar] [CrossRef]
  9. Fang, W.-H.; Yang, J.-R.; Lin, C.-Y.; Hsiao, P.-J.; Tu, M.-Y.; Chen, C.-F.; Tsai, D.-J.; Su, W.; Huang, G.-S.; Chang, H.; et al. Accuracy augmentation of body composition measurement by bioelectrical impedance analyzer in elderly population. Medicine 2020, 99, e19103. [Google Scholar] [CrossRef]
  10. Kim, Y.; Han, B.G. Cohort Profile: The Korean Genome and Epidemiology Study (KoGES) Consortium. Int. J. Epidemiol. 2017, 46, e20. [Google Scholar] [CrossRef]
  11. Park, S. Association between polygenetic risk scores related to sarcopenia risk and their interactions with regular exercise in a large cohort of Korean adults. Clin. Nutr. 2021, 40, 5355–5364. [Google Scholar] [CrossRef]
  12. Wu, X.; Park, S. An Inverse Relation between Hyperglycemia and Skeletal Muscle Mass Predicted by Using a Machine Learning Approach in Middle-Aged and Older Adults in Large Cohorts. J. Clin. Med. 2021, 10, 2133. [Google Scholar] [CrossRef] [PubMed]
  13. Park, S.; Kim, K.; Lee, B.K.; Ahn, J. A Healthy Diet Rich in Calcium and Vitamin C Is Inversely Associated with Metabolic Syndrome Risk in Korean Adults from the KNHANES 2013–2017. Nutrients 2021, 13, 1312. [Google Scholar] [CrossRef] [PubMed]
  14. Yang, H.J.; Song, M.Y.; Kim, M.J.; Park, S. Associations between metabolic syndrome and urinary Na-to-K ratio and glomerular filtration rate in middle-aged adults regardless of Na and K intakes. Clin. Exp. Nephrol. 2020, 24, 1015–1024. [Google Scholar] [CrossRef] [PubMed]
  15. Park, S.; Daily, J.W.; Song, M.Y.; Kwon, H.K. Gene-gene and gene-lifestyle interactions of AKAP11, KCNMA1, PUM1, SPTBN1, and EPDR1 on osteoporosis risk in middle-aged adults. Nutrition 2020, 79–80, 110859. [Google Scholar] [CrossRef]
  16. Park, S.; Kang, S. A Western-style diet interacts with genetic variants of the LDL receptor to hyper-LDL cholesterolemia in Korean adults. Public Health Nutr. 2021, 24, 2964–2974. [Google Scholar] [CrossRef]
  17. Ahn, Y.; Kwon, E.; Shim, J.E.; Park, M.K.; Joo, Y.; Kimm, K.; Park, C.; Kim, D.H. Validation and reproducibility of food frequency questionnaire for Korean genome epidemiologic study. Eur. J. Clin. Nutr. 2007, 61, 1435–1441. [Google Scholar] [CrossRef]
  18. Shivappa, N.; Steck, S.E.; Hurley, T.G.; Hussey, J.R.; Hébert, J.R. Designing and developing a literature-derived, population-based dietary inflammatory index. Public Health Nutr. 2014, 17, 1689–1696. [Google Scholar] [CrossRef] [Green Version]
  19. Hébert, J.R.; Shivappa, N.; Wirth, M.D.; Hussey, J.R.; Hurley, T.G. Perspective: The Dietary Inflammatory Index (DII)-Lessons Learned, Improvements Made, and Future Directions. Adv. Nutr. 2019, 10, 185–195. [Google Scholar] [CrossRef]
  20. Liu, M.; Park, S. A Causal Relationship between Vitamin C Intake with Hyperglycemia and Metabolic Syndrome Risk: A Two-Sample Mendelian Randomization Study. Antioxidants 2022, 11, 857. [Google Scholar] [CrossRef]
  21. Burns, K.F.; Robson, S.M.; Raynor, H.A. A low-energy-dense dietary pattern: A narrative review. Nutr. Rev. 2022, 80, 2301–2311. [Google Scholar] [CrossRef]
  22. Li, X.; Zheng, L.; Zhang, B.; Deng, Z.Y.; Luo, T. The Structure Basis of Phytochemicals as Metabolic Signals for Combating Obesity. Front. Nutr. 2022, 9, 913883. [Google Scholar] [CrossRef] [PubMed]
  23. Mou, Z.; Hyde, T.M.; Lipska, B.K.; Martinowich, K.; Wei, P.; Ong, C.J.; Hunter, L.A.; Palaguachi, G.I.; Morgun, E.; Teng, R.; et al. Human Obesity Associated with an Intronic SNP in the Brain-Derived Neurotrophic Factor Locus. Cell Rep. 2015, 13, 1073–1080. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Nguyen, V.T.; Hill, B.; Sims, N.; Heck, A.; Negron, M.; Lusk, C.; Galindo, C.L. Brain-derived neurotrophic factor rs6265 (Val66Met) single nucleotide polymorphism as a master modifier of human pathophysiology. Neural Regen. Res. 2023, 18, 102–106. [Google Scholar] [PubMed]
  25. Mozaffarian, D. Perspective: Obesity—And unexplained epidemic. Am. J. Clin. Nutr. 2022, 115, 1445–1450. [Google Scholar] [CrossRef] [PubMed]
  26. Daily, J.W.; Park, S. Interaction of BDNF rs6265 variants and energy and protein intake in the risk for glucose intolerance and type 2 diabetes in middle-aged adults. Nutrition 2017, 33, 187–194. [Google Scholar] [CrossRef]
  27. Ma, X.-Y.; Qiu, W.Q.; Smith, C.E.; Parnell, L.D.; Jiang, Z.-Y.; Ordovas, J.M.; Tucker, K.L.; Lai, C.-Q. Association between BDNF rs6265 and Obesity in the Boston Puerto Rican Health Study. J. Obes. 2012, 2012, 102942. [Google Scholar] [CrossRef] [Green Version]
  28. Smemo, S.; Tena, J.J.; Kim, K.H.; Gamazon, E.R.; Sakabe, N.J.; Gómez-Marín, C.; Aneas, I.; Credidio, F.L.; Sobreira, D.R.; Wasserman, N.F.; et al. Obesity-associated variants within FTO form long-range functional connections with IRX3. Nature 2014, 507, 371–375. [Google Scholar] [CrossRef] [Green Version]
  29. Claussnitzer, M.; Dankel, S.N.; Kim, K.-H.; Quon, G.; Meuleman, W.; Haugen, C.; Glunk, V.; Sousa, I.S.; Beaudry, J.L.; Puviindran, V.; et al. FTO Obesity Variant Circuitry and Adipocyte Browning in Humans. N. Engl. J. Med. 2015, 373, 895–907. [Google Scholar] [CrossRef] [Green Version]
  30. Laber, S.; Forcisi, S.; Bentley, L.; Petzold, J.; Moritz, F.; Smirnov, K.S.; Al Sadat, L.; Williamson, I.; Strobel, S.; Agnew, T.; et al. Linking the FTO obesity rs1421085 variant circuitry to cellular, metabolic, and organismal phenotypes in vivo. Sci. Adv. 2021, 7, eabg0108. [Google Scholar] [CrossRef]
  31. Yamaguchi, M.; Murata, T. Involvement of regucalcin gene promoter region-related protein-p117, a transcription factor, in human obesity. Biomed. Rep. 2017, 6, 374–378. [Google Scholar] [CrossRef] [Green Version]
  32. Ng, M.C.; Tam, C.H.; So, W.Y.; Ho, J.S.; Chan, A.W.; Lee, H.M.; Wang, Y.; Lam, V.K.; Chan, J.C.; Ma, R.C. Implication of genetic variants near NEGR1, SEC16B, TMEM18, ETV5/DGKG, GNPDA2, LIN7C/BDNF, MTCH2, BCDIN3D/FAIM2, SH2B1, FTO, MC4R, and KCTD15 with obesity and type 2 diabetes in 7705 Chinese. J. Clin. Endocrinol. Metab. 2010, 95, 2418–2425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Bhattacharyya, D.; Glick, B.S. Two mammalian Sec16 homologues have nonredundant functions in endoplasmic reticulum (ER) export and transitional ER organization. Mol. Biol. Cell 2007, 18, 839–849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Domić, J.; Grootswagers, P.; van Loon, L.J.C.; de Groot, L. Perspective: Vegan Diets for Older Adults? A Perspective On the Potential Impact On Muscle Mass and Strength. Adv. Nutr. 2022, 13, 712–725. [Google Scholar] [CrossRef]
  35. Vang, A.; Singh, P.N.; Lee, J.W.; Haddad, E.H.; Brinegar, C.H. Meats, processed meats, obesity, weight gain and occurrence of diabetes among adults: Findings from Adventist Health Studies. Ann. Nutr. Metab. 2008, 52, 96–104. [Google Scholar] [CrossRef] [PubMed]
  36. Fuente González, C.E.; Chávez-Servín, J.L.; de la Torre-Carbot, K.; Ronquillo González, D.; Aguilera Barreiro, M.; Ojeda Navarro, L.R. Relationship between Emotional Eating, Consumption of Hyperpalatable Energy-Dense Foods, and Indicators of Nutritional Status: A Systematic Review. J. Obes. 2022, 2022, 4243868. [Google Scholar] [CrossRef] [PubMed]
  37. Riedl, R.A.; Burnett, C.M.L.; Pearson, N.A.; Reho, J.J.; Mokadem, M.; Edwards, R.A.; Kindel, T.L.; Kirby, J.R.; Grobe, J.L. Gut Microbiota Represent a Major Thermogenic Biomass. Function 2021, 2, zqab019. [Google Scholar] [CrossRef]
  38. Armstrong, L.E.; Casa, D.J.; Belval, L.N. Metabolism, bioenergetics and thermal physiology: Influences of the human intestinal microbiota. Nutr. Res. Rev. 2019, 32, 205–217. [Google Scholar] [CrossRef]
  39. Li, D.; Li, Y.; Yang, S.; Lu, J.; Jin, X.; Wu, M. Diet-gut microbiota-epigenetics in metabolic diseases: From mechanisms to therapeutics. Biomed. Pharmacother. 2022, 153, 113290. [Google Scholar] [CrossRef]
  40. Mohamadi, A.; Shiraseb, F.; Mirzababaei, A.; Hosseininasab, D.; Rasaei, N.; Clark, C.C.T.; Mirzaei, K. Circulating Inflammatory Markers May Mediate the Relationship between Healthy Plant-Based Diet and Metabolic Phenotype Obesity in Women: A Cross-Sectional Study. Int. J. Clin. Pract. 2022, 2022, 8099382. [Google Scholar] [CrossRef]
  41. Im, J.; Kim, M.; Park, K. Association between the Phytochemical Index and Lower Prevalence of Obesity/Abdominal Obesity in Korean Adults. Nutrients 2020, 12, 2312. [Google Scholar] [CrossRef]
  42. Hansen, T.T.; Astrup, A.; Sjödin, A. Are Dietary Proteins the Key to Successful Body Weight Management? A Systematic Review and Meta-Analysis of Studies Assessing Body Weight Outcomes after Interventions with Increased Dietary Protein. Nutrients 2021, 13, 3193. [Google Scholar] [CrossRef] [PubMed]
  43. Binou, P.; Yanni, A.E.; Kartsioti, K.; Barmpagianni, A.; Konstantopoulos, P.; Karathanos, V.T.; Kokkinos, A. Wheat Biscuits Enriched with Plant-Based Protein Contribute to Weight Loss and Beneficial Metabolic Effects in Subjects with Overweight/Obesity. Nutrients 2022, 14, 2516. [Google Scholar] [CrossRef] [PubMed]
  44. Li, J.; Glenn, A.J.; Yang, Q.; Ding, D.; Zheng, L.; Bao, W.; Beasley, J.; LeBlanc, E.; Lo, K.; Manson, J.E.; et al. Dietary Protein Sources, Mediating Biomarkers, and Incidence of Type 2 Diabetes: Findings from the Women’s Health Initiative and the UK Biobank. Diabetes Care 2022, 45, 1742–1753. [Google Scholar] [CrossRef]
  45. Choi, M.K.; Bae, Y.J. Protein intake and osteosarcopenic adiposity in Korean adults aged 50 years and older. Osteoporos. Int. 2020, 31, 2363–2372. [Google Scholar] [CrossRef] [PubMed]
  46. Orlich, M.J.; Fraser, G.E. Vegetarian diets in the Adventist Health Study 2: A review of initial published findings. Am. J. Clin. Nutr. 2014, 100 (Suppl. S1), 353s–358s. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Slywitch, E.; Savalli, C.; Duarte, A.C.; Escrivão, M. Obese Vegetarians and Omnivores Show Different Metabolic Changes: Analysis of 1340 Individuals. Nutrients 2022, 14, 2204. [Google Scholar] [CrossRef]
Figure 1. The scheme of searching the genetic variants related to body fat mass and polygenic risk factor (PRS) of the selected genetic variants’ interaction with lifestyles. 1 High body fat mass (obesity) was defined as body fat higher than 25% and 30% of body weight for men and women, respectively. 2 Covariates included age, gender, residence area, survey year, daily energy intake, and education and income levels. 3 Covariates for model 1 were age, gender, residence area, survey year, education, and income; those for model 2 were covariates for model 1 plus energy intake, alcohol intake, regular exercise, and smoking status.
Figure 1. The scheme of searching the genetic variants related to body fat mass and polygenic risk factor (PRS) of the selected genetic variants’ interaction with lifestyles. 1 High body fat mass (obesity) was defined as body fat higher than 25% and 30% of body weight for men and women, respectively. 2 Covariates included age, gender, residence area, survey year, daily energy intake, and education and income levels. 3 Covariates for model 1 were age, gender, residence area, survey year, education, and income; those for model 2 were covariates for model 1 plus energy intake, alcohol intake, regular exercise, and smoking status.
Nutrients 15 01063 g001
Figure 2. Distribution of genetic variants for high body fat mass risk by genome-wide association study. (A) Manhattan plot of the p-value of genetic variants for high body fat risk. Red-dot and blue-dot lines indicate the p-value of 5 × 10−8 and 5 × 10−5. (B) Q–Q plot of observed and expected p-values for body fat risk. Red and black lines represent theoretical and actual distributions, respectively.
Figure 2. Distribution of genetic variants for high body fat mass risk by genome-wide association study. (A) Manhattan plot of the p-value of genetic variants for high body fat risk. Red-dot and blue-dot lines indicate the p-value of 5 × 10−8 and 5 × 10−5. (B) Q–Q plot of observed and expected p-values for body fat risk. Red and black lines represent theoretical and actual distributions, respectively.
Nutrients 15 01063 g002
Figure 3. Gene expression according to the alleles of the selected SNPs for high body fat risk in different tissues. (A) ADCY3_rs6545790 in visceral adipose tissues (slope of the SNP alleles = 0.21, p = 4.1 × 10−14). (B) ADCY3_rs6545790 in subcutaneous adipose tissues (slope of the SNP alleles = 0.23, p = 2.2 × 10−19). (C) ADCY3_rs6545790 in the amygdala (slope of the SNP alleles = −0.18, p = 0.0064). (D) ADCY3_rs6545790 in the brain cortex (slope of the SNP alleles = −0.28, p = 5.4 × 10−10). (E) SEC16B_rs509325 in the adrenal gland (slope of the SNP alleles = 0.24, p = 0.0026). (F) BAIAP2_rs35867081 in subcutaneous adipose tissues (slope of the SNP alleles = 0.066, p = 0.025). (G) BAIAP2_rs35867081 in the hippocampus (slope of the SNP alleles = 0.066, p = 0.025). (H) FARP1_rs587056 in the hypothalamus (slope of the SNP alleles = −0.19, p = 0.00072).
Figure 3. Gene expression according to the alleles of the selected SNPs for high body fat risk in different tissues. (A) ADCY3_rs6545790 in visceral adipose tissues (slope of the SNP alleles = 0.21, p = 4.1 × 10−14). (B) ADCY3_rs6545790 in subcutaneous adipose tissues (slope of the SNP alleles = 0.23, p = 2.2 × 10−19). (C) ADCY3_rs6545790 in the amygdala (slope of the SNP alleles = −0.18, p = 0.0064). (D) ADCY3_rs6545790 in the brain cortex (slope of the SNP alleles = −0.28, p = 5.4 × 10−10). (E) SEC16B_rs509325 in the adrenal gland (slope of the SNP alleles = 0.24, p = 0.0026). (F) BAIAP2_rs35867081 in subcutaneous adipose tissues (slope of the SNP alleles = 0.066, p = 0.025). (G) BAIAP2_rs35867081 in the hippocampus (slope of the SNP alleles = 0.066, p = 0.025). (H) FARP1_rs587056 in the hypothalamus (slope of the SNP alleles = −0.19, p = 0.00072).
Nutrients 15 01063 g003
Figure 4. Adjusted odds ratio (ORs) and 95% confidence intervals (CI) of three-SNP PRS and six-SNP PRS for obesity defined with body fat mass. PRS was generated with the sum of the number of risk alleles in each SNP, and it was classified as low PRS, medium PRS, and high PRS according to the PRS was 0–2, 3–4, and ≥5 in the three-SNP model and 0–5, 6–7, and ≥8 in the six-SNP model, respectively. Models 1 and 2 were conducted with different covariates. Covariates were age, gender, residence area, education, and income for models 1 and 2 plus energy intake, alcohol intake, regular exercise, and smoking status.
Figure 4. Adjusted odds ratio (ORs) and 95% confidence intervals (CI) of three-SNP PRS and six-SNP PRS for obesity defined with body fat mass. PRS was generated with the sum of the number of risk alleles in each SNP, and it was classified as low PRS, medium PRS, and high PRS according to the PRS was 0–2, 3–4, and ≥5 in the three-SNP model and 0–5, 6–7, and ≥8 in the six-SNP model, respectively. Models 1 and 2 were conducted with different covariates. Covariates were age, gender, residence area, education, and income for models 1 and 2 plus energy intake, alcohol intake, regular exercise, and smoking status.
Nutrients 15 01063 g004
Figure 5. Percentage of the participants with high body fat mass. (A) According to energy intake. (B) According to protein intake. (C) According to a plant-based diet.
Figure 5. Percentage of the participants with high body fat mass. (A) According to energy intake. (B) According to protein intake. (C) According to a plant-based diet.
Nutrients 15 01063 g005aNutrients 15 01063 g005b
Table 1. (A) The factor-loading values of the predefined 29 food groups in each dietary pattern. (B) Energy intake in each dietary pattern.
Table 1. (A) The factor-loading values of the predefined 29 food groups in each dietary pattern. (B) Energy intake in each dietary pattern.
(A)
Korean Balanced DietPlant-Based DietWestern-Style DietRice-Main Diet
Rice−3 −7 6 93*
Whole grain8 −4 −3 −93*
Bread−1 5 63*1
Cookie−8 37 50 −5
Noodles−8 33 29*6
Bean32 47*3 2
Potato26 49*4 −3
Kimchi50*0 0 −2
Egg8 44*15 4
Fast food−5 18 73*−4
Green vegetables68*40*−2 −1
White vegetables71*27 3 2
Mushroom50*34 −6 −3
Fatty fish53*21 12 0
Whitefish66*15 15 0
Crab48*3 21 1
Processed meats18 15 7 −1
Red meat43*−6 44*7
Chicken16 4 63 −5
Soups32 −5 40*3
Seaweeds45*39*−2 −4
Milk12 49*1 0
Beverage20 31 6 2
Coffee9 0 20 14
Tea12 −7 26 13
Fruits21 46*−6 −5
Pickles50*−2 6 2
Alcohol16 −27 18 5
Nuts−1 50*5 −5
Variance explained by each dietary pattern3.57 2.44 2.30 1.79
(B)
KBDPBDWSDRMD
Energy intake (EER%)125.6 ± 0.434 b109.8 ± 0.254 c178.9 ± 1.193 a87.1 ± 0.143 c
Carbohydrates (En%)64.8_0.103 c69.4 ± 0.061 b63.7 ± 0.287 d73.2 ± 0.04 a
Protein (EN%)17.1 ± 0.037 a14.0 ± 0.02 b14.8 ± 0.103 b12.8 ± 0.012 c
Fat (En%)18.0 ± 0.08 b16.3 ± 0.047 c20.8 ± 0.221 a12.6 ± 0.02 d
SFA (En%)4.15 + 0.13 d5.67 + 0.42 b4.89 + 0.24 c0.64 + 0.12 a
MUFA (En%)5.18 + 0.15 d6.89 + 0.45 b6.13 + 0.257 c8.87 + 0.12 a
PUFA (En%)2.93 + 0.16 d4.06 + 0.38 b3.51 + 0.22 c4.75 + 0.11 a
Fiber (g)26.0 ± 0.12315.1 ± 0.0714.1 ± 0.3413.3 ± 0.04
Vitamin C (mg)169.7 ± 0.885 a125.5 ± 0507 b54.5 ± 2.44 d91.6 ± 0.289 c
Values were factor-loading values. When their absolute values were >0.4, they were flagged by *, and the name of each dietary pattern was assigned from the predominant food groups. EER, estimated energy requirement; En%, energy %; SAF, saturated fatty acid; MUFA, monounsaturated fatty acid; PUFA, polyunsaturated fatty acid. a,b,c,d different letters on the bar indicate significant differences among the groups in Tukey’s test at p < 0.05.
Table 2. Demographic parameters and nutrient intake according to gender and body fat mass (BF).
Table 2. Demographic parameters and nutrient intake according to gender and body fat mass (BF).
Men (n = 19,444)Women (n = 34,384)
Low-BF (n = 16,146)High-BF (n = 3298)Low-BF (n = 27,180)High-BF (n = 7204)
Age (year)55.7 ± 0.06 b56.9 ± 0.13 a51.9 ± 0.05 d53.8 ± 0.09 c***+++###
Education ≤ Middle school 1318 (13.2)343 (16.8) ‡‡‡4290 (20.3)1727 (28.4) ‡‡‡
High school7540 (76.2)1491 (72.8)15,479 (73.4)4074 (67.0)
≥College 1040 (10.5)214 (10.5)1322 (6.27)276 (4.54)
Income
≤USD 20001258 (8.24)251 (7.76) ‡2848 (11.2)868 (12.4) ‡‡‡
USD 2000–40006565 (42.8)1311 (40.5)10,903 (43.0)3369 (48.0)
>USD 40007479 (49.0)1672 (51.7)11,579 (45.7)2777 (39.6)
Energy (EER%) 86.0 ± 0.06 b85.4 ± 0.13 c104 ± 0.05 a104 ± 0.09 a***+++
Carbohydrates (En%) 71.6 ± 0.0871.3 ± 0.1771.7 ± 0.0671.6 ± 0.12 ++
Proteins (En%)13.3 ± 0.03 b13.4 ± 0.05 b13.6 ± 0.02 a13.6 ± 0.04 a***#
Fat (En%) 13.9 ± 0.0614.2 ± 0.1213.9 ± 0.0414.1 ± 0.09
SFA (En%)4.46 ± 0.024.55 ± 0.064.45 ± 0.024.43 + 0.04
MUFA (En%)5.62 ± 0.03 a5.80 ± 0.06 a5.45 ± 0.02 b5.48 ± 0.05 b***+
PUFA (En%)3.26 ± 0.03 a3.21 ± 0.05 a3.12 ± 0.02 b3.07 ± 0.04 b**
Cholesterol (mg/d)169 ± 1.07173 ± 2.25170 ± 0.81171 ± 1.52
Vitamin C (mg/d)93.8 ± 0.68 c91.6 ± 1.50 c121 ± 0.53 a115 ± 1.05 b***++
Vitamin D (ug/d)5.68 ± 0.05 c5.39 ± 0.10 d6.94 ± 0.04 a6.57 ± 0.07 b***+++
Fiber (g/d) 14.6 ± 0.09 b14.1 ± 0.20 b15.3 ± 0.07 a14.4 ± 0.14 b***+++
DII (scores)−19.96 ± 0.02 b−21.07 ± 0.04 b−21.26 ± 0.01 a−21.34 ± 0.03 a**+#
Ca (mg/d)417 ± 3.42 b414 ± 5.56 b491 ± 2.31 a484 ± 3.81 a***
KBD (Yes, %)6169 (38.2)1324(40.2) 7379 (31.1)3057 (28.7) ‡‡‡
PBD (Yes, %)3399 (21.1)693 (21.0)9675 (40.7)4136(38.9) ‡‡
WSD (Yes, %)8218 (50.9)1877 (50.9) ‡‡‡7871 (33.2)3684 (34.6) ‡‡
RMD (Yes, %)5089 (31.5)1016 (30.8)8180 (34.5)3628 (34.1)
Non-Smokers (Yes, %)4712 (29.3)834 (25.3) ‡‡‡22,864 (96.8)10,296 (96.9)
Former smokers (Yes, %)6770 (42.1)1546 (46.9)284 (1.2)132(1.24)
Smokers (Yes, %)4612 (28.7)917 (27.8)532 (1.97)203 (1.91)
Alcohol drinking (g/day)35.1 ± 0.40 b38.9 ± 0.84 a5.33 ± 0.29 c5.51 ± 0.56 c***+++###
Regular exercise (Yes, %)10155 (60.5)1797 (52.4) ‡‡‡16350 (54.0)3674 (45.9) ‡‡‡
Values represent adjusted means ± standard error for continuous variables and the number and percentages of participants for categorical variables. Adjusted for the covariates including age, BMI, education, income, energy intake, alcohol intake, smoking, and total activity. High body fat (case) was defined as 25% and 30% of body weight for men and women, respectively. EER%, the percentage of energy intake based on estimated energy requirement; En%, energy intake percentage. ** Significantly different by genders at p < 0.01; *** at p < 0.001. + Significantly different by fat mass at p < 0.05; ++ at p < 0.01; +++ at p < 0.001. # Significant interactions at p < 0.05; ### at p < 0.001. Significantly different from the low-fat mass group at p < 0.05; ‡‡ at p < 0.01; ‡‡‡ at p < 0.001. a,b,c,d, different letters on the bar indicate significant differences among the groups in Tukey’s test at p < 0.05.
Table 3. Adjusted means and association of anthropometric and biochemical variables in the urban-hospital-based cohort according to genders and body fat mass (BF).
Table 3. Adjusted means and association of anthropometric and biochemical variables in the urban-hospital-based cohort according to genders and body fat mass (BF).
Men (n = 19,444)Women (n = 34,384)Adjusted ORs and 95% CI
Low-BF (n = 16,146)High-BF (n = 3298)Low-BF (n = 27,180)High-BF (n = 7204)
Height (cm) 1 169.0 ± 0.04 b167.1 ± 0.09 a157.0 ± 0.03 d154.9 ± 0.07 c***+++#3.093 (2.881–3.320)
BMI (mg/kg2) 2 23.9 ± 0.02 c26.9 ± 0.05 a22.9 ± 0.02 d26.1 ± 0.03 b***###+++18.05 (17.03–19.14)
Waist (cm) 384.3 ± 0.06 c91.7 ± 0.13 a76.7 ± 0.05 d83.4 ± 0.09 b***###+++5.038 (4.794 5.293)
SMI (%) 4 7.84 ± 0.01 a7.50 ± 0.01 b7.05 ± 0.003 c6.68 ± 0.006 d***###+++0.858 (0.745–0.954)
Fat mass (%) 522.4 ± 0.05 d26.3 ± 0.09 c29.8 ± 0.04 b33.9 ± 0.06 a***+++
MetS (%) 62393 (14.2)1205 (35.1) ‡‡‡2933 (9.65)1769 (22.1) ‡‡‡6.289 (5.833–6.780)
glucose (mg/dL) 798.2 ± 0.17 b99.6 ± 0.37 a93.2 ± 0.13 d94.1 ± 0.25 c***+++1.178 (1.103–1.257)
HbA1c (%) 85.72 ± 0.01 b5.77 ± 0.01 a5.70 ± 0.01 b5.69 ± 0.01 b***##1.454 (1.321–1.601)
Total cholesterol 9190.4 ± 0.31 d194.9 ± 0.68 c199.4 ± 0.24 b207.3 ± 0.46 a***++1.073 (1.001–1.150)
HDL (mg/dL) 1049.2 ± 0.11 c49.8 ± 0.24 c55.8 ± 0.08 a57.6 ± 0.16 b***+++###1.210 (1.135–1.290)
LDL (mg/dL) 11113 ± 0.28 b114 ± 0.62 c121 ± 0.22 a126 ± 0.42 a***+++###1.089 (1.007–1.177)
TG (mg/dL) 12142 ± 0.72 b158 ± 1.59 a115 ± 0.57 c117 ± 1.07 c***+++###1.229 (1.156–1.308)
Hs-CRP (mg/dL) 13 0.153 ± 0.004 b0.187 ± 0.007 a0.118 ± 0.003 c0.154 ± 0.005 b***+++1.266 (1.008–1.589)
SBP (mmHg) 14125.0 ± 0.14 a124.9 ± 0.24 a121.0 ± 0.10 b121.4 ± 0.20 b***1.105 (1.040–1.175)
DBP (mmHg) 1578.1 ± 0.10 a77.9 ± 0.16 a74.4 ± 0.07 b74.6 ± 0.13 b***#1.040 (0.948–1.140)
eGFR (mL/min) 1684.5 ± 0.16 b83.0 ± 0.26 c87.2 ± 0.11 a87.7 ± 0.23 a***+###0.774 (0.687–0.873)
AST (U/L) 1724.6 ± 0.24 b25.7 ± 0.40 a23.1 ± 0.16 c23.2 ± 0.34 c**+++1.425 (1.250–1.624)
ALT(U/L) 1825.1 ± 0.24 b27.9 ± 0.39 a20.5 ± 0.16 c20.9 ± 0.33 b***+++###1.333 (1.227–1.448)
Values represent adjusted means ± standard error for continuous variables and the number and percentages of participants for categorical variables. Adjusted for age, BMI, education, income, energy intake, alcohol intake, smoking, and total activity. High body fat (case) was defined as 25% and 30% of body weight for men and women, respectively. MetS, metabolic syndrome. The cutoff points of the reference for logistic regression with adjusted for age, BMI, education, income, energy intake, alcohol intake, smoking, and total activity. They were as follows: 1 <172.5 cm for men and <160 cm for women; 2 <25 kg/m2 for BMI; 3 <90 cm for men and 85 cm for women waist circumferences; 4 <29.0% for men and 22.8% for women in skeletal muscle index (SMI; defined as appendicular skeletal muscle mass/weight); 5 <25% for men and 30% for women for fat mass; 6 metabolic syndrome (MetS); 7 <126 mL/dL fasting serum glucose plus diabetic drug intake; 8 <6.5% HbA1c plus diabetic drug intake; 9 <230 mg/dL plasma total cholesterol concentrations; 10 >40 mg/dL for men and 50 mg/dL for women plasma HDL cholesterol; 11 <160 mg/dL plasma LDL cholesterol concentrations; 12 <150 mg/dL plasma triglyceride (TG) concentrations; 13 <0.5 mg/dL serum high-sensitive C-reactive protein (hs-CRP) concentrations; 14 <140 mmHg SBP; 15 <90 mmHg DBP plus hypertension medication; 16 estimated glomerular filtration rate (eGFR) < 70; 17 aspartate aminotransferase < 40 U/L; 18 alanine aminotransferase <35 U/L. ** Significantly different by genders at p < 0.01; *** at p < 0.001. ++ Significantly different by fat mass at p < 0.01; +++ at p < 0.001. # Significant interactions at p < 0.05; ## at p < 0.01; ### at p < 0.001. ‡‡‡ Significantly different from the low-fat mass group at p < 0.001. a,b,c,d, different letters on the bar indicate significant differences among the groups in Tukey’s test at p < 0.05.
Table 4. Pathways related to genetic variants for body fat mass.
Table 4. Pathways related to genetic variants for body fat mass.
PathwaysNo. of GenesBetaStdp-Valuep-Value for Bonferroni Correction
GO BP: GO sodium ion transmembrane transport1330.3230.0270.0813.18 × 10−5
GO BP: GO myoblast proliferation170.6930.0210.1795.41 × 10−5
Curated gene sets: Biocarta flumazenil pathway81.2080.0250.3260.000104
GO BP: GO Neuron apoptotic process2240.2110.0230.0580.000129
Curated gene sets: Reactome RMTS methylate histone arginine600.3870.0220.1080.000172
GO MF: GO Peptide hormone binding470.4880.0240.1370.000183
GO BP: GO Positive regulation of vascular endothelial cell proliferation130.7780.0200.2190.000193
No., number; GO, gene ontology; BP, biological process; MF, molecular function; Std, standard deviation; Bon, Bonferroni correction.
Table 5. The characteristics of the ten genetic variants of genes related to body fat mass in adults using the generalized multifactor dimensionality reduction analysis.
Table 5. The characteristics of the ten genetic variants of genes related to body fat mass in adults using the generalized multifactor dimensionality reduction analysis.
CHR 1SNP 2Location Mi 3Ma 4OR 5SE 6p-Value for OR 7p-Value for OR 8GenesFeature9 MAF10 HWE
1rs509325177894591GT1.123 0.01451.37 × 10−158.28 × 10−4SEC16BIntron0.28450.5035
2rs6545790 25109302AG1.0630.01334.87 × 10−64.75 × 10−4ADCY3Intron0.43740.4064
2rs756057554142030CT0.98890.06702.33 × 10−71.94 × 10−2PSME4Transcript0.01120.2631
4rs219647620270600GA1.0970.01911.26 × 10−60.0069SLIT2Intron0.1350.4579
11rs6265 27679916CT0.9230.01332.45 × 10−100.0048BDNFMissense0.45880.1481
13rs587056 98976374TC1.350.06362.33 × 10−100.0067FARP1Intron0.01060.5496
16rs1421085 53800954TC1.1730.01976.24 × 10−162.82 × 10−6FTOTranscript0.12450.4604
17rs35867081 79047278GA1.0650.01373.99 × 10−60.0052BAIAP2Transcript0.36510.2094
19rs60259426 46340832GA0.9380.01341.54 × 10−70.00072SYMPKTranscript0.42020.3135
20rs6089240 60152260AG0.92780.01261.60 × 10−85.81 × 10−4CDH4Intron0.46411.0
1 Chromosome; 2 single-nucleotide polymorphism; 3 minor allele; 4 major allele; 5 odds ratio for body fat mass by GWAS in a city-hospital-based cohort; 6 standard error for body fat mass by GWAS in a city-hospital-based cohort; 7 p-value for OR for body fat mass after adjusting for age, gender, residence area, survey year, daily energy intake, levels of education, and income in a city-hospital-based cohort; 8 p-value for OR for body fat mass after adjusting for covariates stated above in an Ansan/Ansung cohort; 9 minor allele frequency; 10 Hardy–Weinberg equilibrium.
Table 6. Adjusted odds ratios for fat mass risk by polygenetic risk scores of the three-SNP model (PRS) for gene–gene interaction after covariate adjustments according to the patterns of lifestyles.
Table 6. Adjusted odds ratios for fat mass risk by polygenetic risk scores of the three-SNP model (PRS) for gene–gene interaction after covariate adjustments according to the patterns of lifestyles.
Low-PRS
(n = 19,686)
Medium-PRS
(n = 30,513)
High-PRS
(n = 3629)
Gene–Nutrient Interaction
p-Value
Low energy 1
High energy
11.150 (1.081–1.223)
1.102 (1.021–1.189)
1.228 (1.088–1.387)
1.157 (0.997–1.348)
0.0231
Low CHO 2 70
High CHO
11.045 (0.893–1.223)
0.991 (0.927–1.060)
1.006 (0.866–1.169)
1.079 (1.014–1.149)
0.3493
Low protein 3 13
High protein
11.123 (1.052–1.198)
0.961 (0.881–1.049)
1.255 (1.103–1.427)
1.019 (0.938–1.107)
0.0272
Low fat 4 15
Moderate fat
11.000 (0.926–1.080)
0.994 (0.896–1.102)
1.088 (1.013–1.169)
1.028 (0.932–1.133)
0.0664
Low alcohol 5 20
High alcohol
11.139 (1.075–1.208)
1.129 (1.050–1.213)
1.218 (1.086–1.366)
1.204 (1.046–1.386)
0.8083
Low KBD 6
High KBD
10.997 (0.938–1.061)
1.002 (0.929–1.080)
1.067 (1.007–1.130)
1.083 (1.009–1.162)
0.5368
Low PBD 6
High PBD
11.132 (1.081–1.185)
1.023 (0.961–1.088)
1.258 (1.126–1.406)
1.110 (1.006–1.224)
0.0026
Low WSD 6
High WSD
10.997 (0.938–1.061)
0.998 (0.927–1.075)
1.067 (1.007–1.130)
1.074 (1.002–1.151)
0.1356
Low RMD 6
High RMD
11.112 (1.057–1.169)
1.159 (1.096–1.225)
1.197 (1.134–1.263)
1.254 (1.128–1.395)
0.8419
Low exercise 7
High exercise
11.034 (0.949–1.128)
0.997 (0.938–1.061)
1.103 (1.018–1.196)
1.067 (1.007–1.130)
0.1778
Non-smoke 8
Smoke
10.997 (0.938–1.061)
1.071 (0.948–1.210)
1.067 (1.007–1.130)
1.041 (0.926–1.169)
0.1328
Values were expressed as odds ratio and 95% confidence intervals. PRS of the three-SNP model (BDNF_rs6265, FTO_rs1421085, and SEC16B_rs509325) was divided into three categories (0–2, 3–4, and ≥5) by three groups as the low, medium, and high groups of the best model of GMDR. The cutoff points of the parameters were defined as follows: 1 <estimated energy intake, 2 <70% carbohydrate (CHO), 3 <13% protein, 4 <15% fat, 5 <20 g/day alcohol, 6 <75th percentile, 7 30 min moderate exercise for 3 times a week, and 8 smoking. Multivariate regression models include the main effects, interaction terms of gene and main effects (energy and nutrient intake), and potential confounders, including sex, age, BMI, the status of smoking and drinking, levels of income and education, physical activity, energy intake, and percent intake for carbohydrate and fat. The reference was low PRS. KBD, Korean-style balanced diet; PBD, plant-based diet; WSD, Western-style diet; RMD, rice-main diet.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Daily, J.W.; Park, S. Association of Plant-Based and High-Protein Diets with a Lower Obesity Risk Defined by Fat Mass in Middle-Aged and Elderly Persons with a High Genetic Risk of Obesity. Nutrients 2023, 15, 1063. https://doi.org/10.3390/nu15041063

AMA Style

Daily JW, Park S. Association of Plant-Based and High-Protein Diets with a Lower Obesity Risk Defined by Fat Mass in Middle-Aged and Elderly Persons with a High Genetic Risk of Obesity. Nutrients. 2023; 15(4):1063. https://doi.org/10.3390/nu15041063

Chicago/Turabian Style

Daily, James W., and Sunmin Park. 2023. "Association of Plant-Based and High-Protein Diets with a Lower Obesity Risk Defined by Fat Mass in Middle-Aged and Elderly Persons with a High Genetic Risk of Obesity" Nutrients 15, no. 4: 1063. https://doi.org/10.3390/nu15041063

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