Next Article in Journal
Nitrogen Metabolism and Antioxidant Capacity of Selected Vegetables from Organic and Conventional Crops
Next Article in Special Issue
Phase Angle and Bioelectrical Impedance Vector Analysis (BIVA) in Amyotrophic Lateral Sclerosis (ALS) Patients
Previous Article in Journal
A Study on the Load-Bearing Characteristics and Load Transfer Mechanism of Bag Grouting Pile in Soft Soil Areas
Previous Article in Special Issue
Relationship between Body Composition and Gross Motor Coordination in Six-Year-Old Boys and Girls
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dietary Behavior of Spanish Schoolchildren in Relation to the Polygenic Risk of Obesity

by
Andrea Calderón García
1,2,3,
Roberto Pedrero Tomé
2,4,
Ana Alaminos-Torres
2,4,
Consuelo Prado Martínez
2,5,
Jesús Román Martínez Álvarez
1,2,
Noemí López Ejeda
1,2,4,
María Dolores Cabañas Armesilla
2 and
María Dolores Marrodán Serrano
1,2,4,*
1
Sociedad Española de Dietética y Ciencias de la Alimentación (SEDCA), 28080 Madrid, Spain
2
Research Group EPINUT (Nutritional Epidemiology), Faculty of Medicine, Complutense University of Madrid, 28040 Madrid, Spain
3
Department of Nursing and Nutrition, Faculty of Biomedical Sciences, European University of Madrid, UEM, Villaviciosa de Odón, 28670 Madrid, Spain
4
Department of Biodiversity, Ecology and Evolution, Faculty of Biology, Complutense University of Madrid, 28040 Madrid, Spain
5
Department of Biology, Faculty of Sciences, Universidad Autónoma de Madrid, 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(20), 11169; https://doi.org/10.3390/app132011169
Submission received: 13 September 2023 / Revised: 1 October 2023 / Accepted: 7 October 2023 / Published: 11 October 2023

Abstract

:
Several precedents support an association between single nucleotide genetic polymorphisms (SNPs), the obese phenotype, and eating behavior in the infant-juvenile population. This study aims to study this aspect in depth, analyzing the eating behavior of a sample of schoolchildren from Madrid in regard to their genetic predisposition to obesity. A total of 258 schoolchildren, aged 6 to 16 years, were evaluated through the Children’s Eating Behaviour Questionnaire (CEBQ) and the genotyping of 32 SNPs. Associations were observed between the total genetic risk score and eating behaviors related to emotional eating and food responsiveness. Individually, different SNPs were associated with eating behaviors, primarily those related to pro-eating behaviors or increased risk of developing obesity. However, diverse results are obtained, depending on the SNP. These results highlighted the strongest associations between the rs1801725 SNP risk allele (CASR) and increased enjoyment of food by 1846-fold. Likewise, the satiety response was associated with SNP rs11676272 (ADCY3) 2.39 and SNP rs7566605 (INSIG2), increasing this response by 2.39 and 1.63 times, respectively. Emotional anti-ingesting behaviors were inversely associated with SNP rs1421085 (FTO) and SNP rs987237 (TFAP2B). In contrast, SNP rs55915917 (CRHR1) increased the risk of these behaviors. SNPs rs4788099 (SH2B1), rs6857 (NECTIN2), and rs180172 (ADCY3) were associated with slow feeding. In conclusion, associations were found between most of the analyzed SNPs and the CEBQ items. This suggests that feeding behavior exists as a mediator between genotype and obesity phenotype, beginning in childhood.

1. Introduction

Obesity in children and adolescents is a public health problem with a multifactorial origin in which eating behavior plays a role [1]. This behavior, although influenced by the environment, can also present as genetic conditioning [2], which can be analyzed through the detection of single nucleotide polymorphisms (SNPs) that are associated with the obese phenotype, mediated by eating behavior [3].
Numerous studies support the association between SNPs and the risk of developing overweight or obesity. This is explained by the fact that specific genotypes seem to condition attitudes towards food, such as appetite or satiety, speed of intake, frequency of food consumption, or the tendency to experience emotional hunger. With this approach, research groups such as Jääskelainen et al. [4], Llewellyn et al. [5], Steinbeck et al. [6], and Monnereau et al. [7] have associated eating behavior with polygenic risk, quantified in a range between 8 and 32 SNPs.
Other investigations have related obesity with lower responsiveness to internal satiety signals and greater receptivity to olfactory or gustatory stimuli, which would lead to a higher caloric intake [8]. It has also been shown that an infant’s feeding self-regulation arises from genetic–environmental interactions responsible for behaviors such as rapid eating and delayed satiety [8,9]. Hardle et al. [10] investigated the role of eating behavior in genetic susceptibility to obesity, and their thorough literature analysis over the decade of 2010–2020 concluded that emotional and uncontrolled eating behaviors lead to being overweight.
A pioneering study from the United Kingdom by Wardle et al. [11] investigated the association of the SNP rs9939609 genotype of the FTO gene with adiposity and satiety control indicators in children. It proved the usefulness of psychometric scales (Satiety et al. of Food and the Child Eating Behavior Questionnaire) as objective measures of eating behavior. This is also emphasized by Hardle’s literature review [10], in which the authors highlight the value of psychometric tests that objectively measure attitudes toward food.
Undoubtedly, the development of eating habits is a process with a complex etiology, including psychological aspects, particularly the influence of family, social, and affective factors. Some research has highlighted a particular circumstance, reporting that in overweight subjects, there is a positive relationship between an emotional eating pattern and body mass index [12]. Hence, food is eaten to satisfy physiological hunger, to fill an inner emptiness, or to comfort the mood.
Determining the interactions between genotype and eating behavior in greater depth is particularly interesting in childhood and adolescence because it is at these stages when most future habits are shaped. This is a window of opportunity for prevention, and eventually, for an individualized nutritional approach. Therefore, this study aims to analyze the eating behavior of a sample of Madrid schoolchildren concerning their genetic predisposition to obesity through the conduction of the Children’s Eating Behavior Questionnaire (CEBQ), in combination with a genotyping assessment from a panel of 32 SNPs.

2. Methodology

2.1. Participants

Data were collected between 2019 and 2021 during a cross-sectional descriptive study in which 258 schoolchildren, aged 6–16 (67.83% male; 32.17% female), from different communities of Madrid, Spain, participated. The data collected were anonymized and disaggregated to remove information that could identify the subject. In order to be included, informed consent was obtained from the parents or guardians of each subject, respecting the bioethical principles of the Declaration of Helsinki in its most updated version [13]. The Ethics Committee of the Universidad Autónoma de Madrid approved the study.

2.2. Instruments

Each participant was assessed anthropometrically, and their polygenic risk of obesity was calculated by genotyping 32 SNPs. Parents or guardians completed the Child Eating Behaviour Questionnaire (CEBQ) to assess eating the behavior of the subjects [14].

2.2.1. CEBQ Questionnaire

The CEBQ [14] is a validated questionnaire that indicates the participant’s responses to satiety, taste for food, speed of intake, and emotional food consumption, among other parameters. It consists of 35 items that evaluate eight subscales of eating behavior and whose questions are answered using a Likert-type scale, with an option to score each item from 1 to 5, according to the intensity of the behavior: never (1), seldom (2), sometimes (3), often (4), always (5).
The items are classified into eight subscales: food responsiveness (FR; 5 items), enjoyment of food (EF; 4 items), emotional overeating (EOE; 4 items), desire to drink (DD; 3 items), slowness in eating (SE; 4 items), satiety responsiveness (SR; 5 items), food fussiness (FF; 6 items), and emotional undereating (EUE; 4 items). The first four subscales have a positive focus or are related to increased food intake (pro-eating dimension). In contrast, the last four subscales are related to food avoidance or negative food-related responses (anti-eating dimension). The present study used a Spanish-language translated and adapted version of the CEBQ that was previously validated among Spanish and Chilean populations [15].

2.2.2. Genetic Analysis

A total of 32 SNPs were analyzed (Table A1).
The extraction protocol employed was that of BioTools B&M S.A. Laboratories for saliva [16], and the Speedtools Tissue DNA Extraction kit was used. Subsequently, genotyping was carried out in the laboratories of the University of Santiago de Compostela Node of the National Genotyping Center (CeGen), which is part of the Network Platform of Biomolecular and Bioinformatics Resources of the Carlos III Institute of Health. The iPLEX® Gold technology for the MassArray platform of Agena Bioscience Inc. [17] was used. Once genotyping was completed, all samples were burned for the data protection of the participants.
A total of 32 SNPs were analyzed: all risk alleles of the selected SNPs have been associated with the childhood and adolescent obese phenotype or eating behaviors in previous research [16,17,18]. The SNPs selected are reported in Table A1 (Supplementary file: Table A1).
The extraction protocol employed was that of BioTools B&M S.A. Laboratories for saliva [19], and the Speedtools Tissue DNA Extraction kit was used. Subsequently, genotyping was carried out in the laboratories of the University of Santiago de Compostela Node of the National Genotyping Center (CeGen), which is part of the Network Platform of Bio-molecular and Bioinformatics Resources of the Carlos III Institute of Health. The iPLEX® Gold technology for the MassArray platform of Agena Bioscience Inc. (Hamburg, Germany) [20] was used. Once genotyping was completed, all samples were burned for the data protection of the participants.

Calculation of the Genetic Risk for Obesity

The genotype of each SNP was ranked by awarding 2 points if the individual was homozygous for the risk allele, 1 point if heterozygous, and 0 points if homozygous for the alternate allele. The total genetic risk score (GRS) was obtained from the sum of the scores of the 32 SNPs. The simplest possible model uses unweighted scores, which is appropriate in this type of study, especially when the sample is insignificant [16]. The total score can be obtained from 0 to 64 points (0: not a single risk allele; 64: homozygous for risk at all 32 SNPs). Quartiles were calculated for GRS by establishing four categories, from lowest to highest risk (GRS ≤ Q1; GRS > Q1–≤Q2; GRS > Q2–≤Q3; PRG > Q3).

2.2.3. Anthropometric Study

To evaluate nutritional status, height (cm), weight (kg), waist circumference (cm), and bicipital, tricipital, subscapular, and suprailiac adipose skinfolds (mm) were measured. For the prevalence analysis, the sample was stratified by sex, and nutritional categories were established based on the body mass index [BMI= weight (kg)/height (m2)] using the cutoff points of Cole et al. [17] and the waist-to-hip ratio [WHR = waist circumference (cm)/height (cm)], using the cutoff points established by Marrodán et al. [18]. The percentage of body fat (%BF) was calculated using the formulas of Siri 1961 [17] after calculating the density (Brook, 1971; Durnin and Rahaman, 1967) [21,22]. The degree of adiposity and its distribution were classified following the patterns of Marrodán et al. [23].

2.3. Statistical Procedures

Cronbach’s alpha was applied to evaluate the internal consistency of the CEBQ subscales, which was found to be adequate (greater than 0.7). Kruskal–Wallis tests were performed to determine whether GRS and individual SNPs were associated with these subscales, followed by linear logistic regression tests. A CEBQ score ≥ 50th percentile was the dependent variable, and at least one reference allele at each SNP was used as a reference point to measure genetic risk. All models were adjusted for sex and age. Statistical processing was performed using the IBM SPSS 24 statistical package.

3. Results

3.1. Anthropometric and Genetic Description of the Sample

The anthropometric description of the sample is shown in Table 1. The genotypic frequencies and risk allele frequency of each selected SNP are presented in Table A2, and it can be observed that no allele showed a frequency lower than 0.05, which corresponds with the minimal data quality criterion for the association study (Hardy–Weinberg equilibrium).

The Genetic Risk Score and Prevalence of Overweight and Adiposity

The resulting genetic risk categories were Q1: ≤19 points; Q2: 20–22 points; Q3: 23–25 points; and Q4: >25 points.
As shown in Table 2, the proportion of overweight and obesity categorized by BMI is significantly higher among subjects with a higher genetic predisposition to obesity (p = 0.040). The WHR and %BF trend also reflects a higher prevalence of overweight or adiposity in genetically predisposed individuals, without significance.

3.2. Genetic Risk Score for Obesity and Eating Behavior

Based on the sample participant’s anthropometric and genotypic description, the resulting genetic risk score (GRS) categories were Q1: ≤19 points; Q2: 20–22 points; Q3: 23–25 points; and Q4: >25 points. Table 3 presents the mean scores of the eight subscales of the CEBQ, stratifying the results according to GRS quartile categories. The trend is increasing for the pro-eating scales and decreasing for the anti-eating scales.
The pro-eating subscales of food responsiveness (FR) and emotional overeating (EOE), as well as the anti-eating subscale of emotional undereating (EUE), were the eating behavioral traits significantly associated with GRS. As genetic predisposition to obesity increases, items related to food cravings and emotional hunger score higher.
The multinomial logistic regression model highlighted associations between all pro-intake subscales and selected SNPs, analyzed independently (Figure 1).
The SNPs rs1801725 (CASR) and rs1137101 (LEPR) were independently associated with enjoyment of food (EF), showing, in the first case, an increase in the score of this scale with an OR of 1.846 (95% CI: 1.023–3.329; p = 0.024). In the second case, the same scale (EF) was reduced, with an OR of 0.586 (95% CI: 0.329–0.980; p = 0.024).
Food responsiveness (FR) was increased by 1.827-fold (95% IC: 1.090–3.062; p = 0.041) for each risk allele in SNP rs10887741 (PAPSS2) and by 2.224-fold in SNP rs1421085 (95% IC: 1.117–4.427 p = 0.042). In contrast, for each risk allele in SNP rs693839 (SPRY2), FR was reduced, with an OR of 0.519 (95% IC: 0.309–0.872; p = 0.019).
Emotional overeating (EOE) subscale scores were increased by 1.628-fold for each risk allele of SNP rs10887741 (PAPSS2) (95% IC: 1.016–2.698; p = 0.049). Finally, the desire to drink (DD) subscale score was reduced for each risk allele of SNP rs3761445 (PLA2G6), with an OR of 0.501 (95% IC: 0.291–0.864; p = 0.021).
To investigate the degree to which the risk allele of each SNP increases the possible obesogenic behavior, the inverse OR values were plotted, corresponding to the anti-intake subscales, when constructing Figure 2. Thus, the obesogenic behaviors are positioned on the right and the non-obesogenic behaviors on the left of Figure 2.
Associations were found with several SNPs except for the food fussiness (FF) subscale. Each risk allele of SNP rs11676272 (ADCY3) reduced satiety responsiveness (SR) by 2.390-fold (95% CI: 1.367–4.184; p < 0.001). Similarly, each risk allele of SNP rs7566605 (INSIG2) reduced it by 1.627-fold (95% CI: 1.060–2.775; p = 0.033).
The SNPs rs4788099 (SH2B1), rs6857 (NECTIN2), and rs180172 (ADCY3) showed a decrease in slowness in eating (SE), or in other words, an increase in eating speed of 1.65 (95% CI: 1.258–2.945; p = 0.010), 2.06 (95% CI: 1.179–3.709; p = 0.008), and 2.37 times (95% CI: 1. 124–5.000; p = 0.025) respectively. In contrast, for each copy of the risk allele of rs987237 (TFAP2B), intake speed was reduced, with an OR of 0.492 (95% CI: 0.263–0.921 p = 0.042).
Emotional undereating (EUE) was nominally associated with 4 SNPs. The risk alleles of SNPs rs55915917 (CRHR1) and rs7132908 (FAI2M) reduced the probability of emotional undereating (EUE) by 1.779 (95% CI: 1.044–3.030; p = 0.014) and 1.751 (95% CI: 1.007–3.719; p = 0.040) times, respectively. In contrast, the risk alleles corresponding to SNPs rs1421085 (FTO) and rs987237 (TFAP2B), indicated on the left side of the graph, were associated with emotional undereating anti-intake behavior, with an OR of 0.385 (95% CI: 0.189–0.781; p = 0.009) and 0.391 (95% CI: 0.214–0.713; p = 0.011), respectively. Finally, for each risk allele of SNP rs6857 (INSIG2), the lower food fussiness (FF) is increased by 2.303 times (1.045–5.142; p = 0.023).

4. Discussion

The scientific literature supports that childhood eating behaviors may be modulated by genetic factors [25]. The behavioral susceptibility theory hypothesizes that genetics influences the obesity phenotype by regulating appetite, expressed in eating behavior. The review by Herle et al. [10] is illuminating, as it highlights the value of studies associating the genetic risk scores and the psychometric measures of eating behaviors. The same authors highlight the need for further research, particularly longitudinal studies in children.
In the present study, an association was found between GRS and obesity and the food responsiveness (FR) subscale, in addition to the subscales related to emotional eating, both by excess (EOE) and deficiency (EUE). The results using other subscales are generally in line with the expectations that the average scores increase with higher genetic risk scores in the pro-intake scales and decrease in the anti-intake scales, following the trend of previous studies, which highlight the more significant evidence concerning the ability to satiate and the demand to food [10].
Some studies have observed that schoolchildren genetically predisposed to obesity show a more significant response to or enjoyment of food, greater appetite, and less responsiveness [24,26]. In the present study, despite not finding such an association with the total risk score, the trend of the averages is along the same lines.
In turn, in the literature, a more significant response to food, understood as a greater desire to consume food, is associated with a greater risk of being overweight. The literature concludes that children who enjoy or express a greater desire for food consume more daily energy on average [27]. These schoolchildren also eat more processed products of low nutritional quality, which are high in calories, refined and saturated fats, and free or added sugars, a dietary pattern associated with an increased risk of obesity [28].
Genetic predisposition to emotional eating is controversial. Further evidence is still required in school populations, and previous studies reached contradictory results. Some authors suggest that the genetic component predisposes a greater or lesser tendency to emotional eating [29], an observation in line with the findings of the present study. However, other studies find that the heritability of the emotional eating behavior is low and, in any case, is less deterministic than learned behavior [30,31]. Certain studies report associations between a lower genetic risk for obesity and slower eating, leading to earlier satiety and lower food intake, possibly because of less responsiveness to internal satiety cues and greater appetite [32]. Our study finds both trends, although without significant association with GRS.
Children who are pickier and who reject new foods might have a less varied and obesogenic diet [10]. However, the present study found no relationship between food pickiness and genetic predisposition to obesity, as previous studies have described [33].
On the other hand, 14 of 32 SNPs analyzed could be independently associated with eating behaviors. Thus, five SNPs from different loci were associated with the anti-eating subscale of emotional undereating (EUE) and two other SNPs with the pro-intake subscale of increased risk of emotional overeating (EOE). In particular, the SNPs of the FTO gene, rs1421085 (significant in the present study) and rs9939609, have been associated in the literature with uncontrolled eating (loss of control eating), which is usually accompanied by a preference for more palatable foods, usually higher in calories and rich in fat or sugar [34].
In the anti-intake subscale of slowness in eating (SE), lower scores have also been found in children with the risk genotype in three SNPs, implying a possible higher eating speed. Scientific evidence shows that eating slowly reduces food intake, a protective factor against weight gain. Eating slowly promotes earlier satiety, while eating fast means higher total intake [6]. A recent study by Jackson et al. [35] concludes that children with lower scores on the slowness to eat subscale of the CEBQ have a higher daily energy intake, around 169 kcal on average, compared to those with higher scores.
In our study, schoolchildren with the risk allele in the SNP rs11676272 (ADCY3) exhibited lower average scores in the anti-eating subscale of satiety responsiveness (SR). The SNP rs11676272 (ADCY3) has been previously associated with lower satiety capacity, possibly due to alterations in the regulation of intake and energy expenditure [36]. The study by Jackson et al., mentioned previously [37], concludes that children with lower scores in this subscale of the CEBQ have a higher daily energy intake of 83 kcal on average. In the literature, other SNPs associated with lower satiety capacity have not reached significant association in our study, highlighting those in or near the FTO, MC4R, GNPDA2, LMX1B, or SH2B1 genes [27,38].
The limitations of this study include a relatively small sample size. It would be interesting to replicate the study in a larger population because this may reveal more associations between other eating behaviors and SNPs that have not obtained significance in this study. Similarly, due to the young age of the participants, there may sometimes be limited control over their intake and eating behaviors. However, the continuity of such behaviors in early life has been demonstrated. Likewise, genetic predisposition may explain some variations in eating behaviors during childhood and adolescence. However, environmental influences undoubtedly play a fundamental role, especially in adolescence, when environmental factors are of significant importance [39].

5. Conclusions

No strong associations were found between GRS and obesity and most eating behaviors. However, the present study highlights the role in subscales related to emotional eating, either by over or under eating. In addition, a potential relationship was found between the emotional eating score and different SNPs in regards to eating behaviors, again highlighting its influence on emotional underfeeding, the speed of intake and enjoyment, and response to food. The relationship between genetic phenotypes and eating behavior could be consolidated as an approach of interest in the future for the comprehensive approach to evaluating obesity from an early age.

Author Contributions

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

Funding

This research was funded by Project PR41/17_21008 BANCO DE SANTANDER.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of the Autonomous University of Madrid (CEI-91-1699-2018).

Informed Consent Statement

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

Data Availability Statement

The original data are not published in any repository.

Acknowledgments

We thank the participating children, their families, and their teachers for their contributions to this study.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Description of the genotypic and risk allele frequencies of the 32 SNPs selected for inclusion in the study.
Table A1. Description of the genotypic and risk allele frequencies of the 32 SNPs selected for inclusion in the study.
ChromosomeAssociated GeneSNPPositionAlelles *
1LEPRrs113710165.592.830A/G
1SEC16Brs543874177.920.345A/G
1GPR61rs7550711110.082.886C/T
2TMEM18rs6548238634.905C/T
2TMEM18rs4854349647.861T/C
2INSIG2rs7566605118.836.025C/G
2ADCY3rs1167627224.918.669A/G
2COBLL1rs6738627164.687.940G/A
3CASRrs1801725122.284.910G/T
4GNPDA2rs1093839745.182.530A/G
4CLOCKrs180126056.301.369A/G
6TFAP2Brs98723750.835.337A/G
7EXOC4rs7804463133.762.898T/C
8ELP3rs1325311128.204.457A/G
9FAM120AOSrs94499096.191.004C/T
9LMX1Brs3829849126.628.521C/T
10PAPSS2rs1088774187.683.553T/C
12FAI2Mrs713290849.869.365G/A
12FAIM2rs713880350.242.468A/G
13OLFM4rs1242954554.102.206A/G
13SPRY2zrs69383980.384.153T/C
16FTOrs155890253.803.547A/T
16FTOrs1781744953.813.367G/T
16FTOrs993960953.820.527A/T
16FTOrs142108553.767.042T/C
16IRX3rs375172354.286.285G/T
16SH2B1rs478809928.844.406A/G
18MC4Rrs656716060.161.902T/C
18RAB27Brs809250354.812.256A/G
19NECTIN2rs685744.888.997T/C
19CRTC1rs75731818.709.498C/A
22PLA2G6rs376144538.199.404G/A
* The risk allele is highlighted in bold.
Table A2. Description of the genotypic and risk allele frequencies of the 32 obesity-associated polymorphisms included in the study sample.
Table A2. Description of the genotypic and risk allele frequencies of the 32 obesity-associated polymorphisms included in the study sample.
SNPGenotype Frequency (%)FAREuropean FAR *H-W Equilibrium
rs1137101 (LEPR)
Risk allele: G
AA: 33.20AG: 48.00GG: 18.800.430.45p = 0.844
rs543874 (SEC16B)
Risk allele: G
AA: 69.50AG: 27.70GG: 2.800.170.18p = 0.984
rs7550711 (GPR61)
Risk allele: T
CC: 94.90CT: 0.00TT: 5.100.070.04p = 0.070
rs6548238(TMEM 18)
Risk allele: T
CC: 70.80CT: 27.40TT: 1.800.160.15p = 0.646
rs4854349 (TMEM18)
Risk allele: C
TT: 3.10TC: 33.20CC: 63.700.820.83p = 0.621
rs7566605 (INSIG2)
Risk allele: C
GG: 43.80GC: 43.60CC: 12.600.340.32p = 0.734
rs11676272 (ADCY3)
Risk allele: G
AA: 27.70AG: 53.50GG: 18.800.460.48p = 0.432
rs6738627 (COBLL1)
Risk allele: A
GG: 39.80GA: 50.40AA: 9.800.350.36p = 0.281
rs1801725 (CASR)
Risk allele: T
GG: 73.40GT: 25.00TT: 1.600.140.15p = 0.748
rs10938397(GNPDA2)
Risk allele: G
AA: 32.20AG: 50.00GG: 17.800.430.37p = 0.838
rs1801260 (CLOCK)
Risk allele: G
AA: 53.40AG: 37.00GG: 9.600.280.29p = 0.399
rs987237 (TFAP2B)
Risk allele: G
AA: 70.30AG: 27.30GG: 2.300.160.19p = 0.854
rs7804463 (EXOC4)
Risk allele: C
TT: 25.80TC: 53.10CC: 21.100.480.46p = 0.519
rs13253111 (ELP3)
Risk allele: G
AA: 23.80AG: 50.00GG: 26.200.510.44p = 0.995
rs944990 (FAM120AOS)
Risk allele: T
CC: 50.50CT: 42.70TT: 6.800.280.28p = 0.578
rs3829849 (LMX1B)
Risk allele: T
CC: 39.50CT: 48.40TT: 12.100.360.36p = 0.217
rs10887741 (PAPSS2)
Risk allele: C
TT: 46.10TC: 40.20CC: 13.700.340.39p = 0.984
rs7132908 (FAI2M)
Risk allele: A
GG: 35.90GA: 49.60AA: 14.500.340.39p = 0.309
rs7138803 (FAIM2)
Risk allele: A
GG: 42.20GA: 42.90AA: 14.080.360.38p = 0.555
rs12429545 (OLFM4)
Risk allele: A
GG: 74.00GA: 23.50AA: 2.500.140.07p = 0.701
rs693839 (SPRY2z)
Risk allele: C
TT: 45.70TC: 42.60CC: 11.700.330.31p = 0.714
rs1558902 (FTO)
Risk allele: A
TT: 34.50TA: 44.70AA: 20.800.430.42p = 0.374
rs17817449 (FTO)
Risk allele: G
TT: 36.10GT: 44.30GG: 19.600.420.41p = 0.374
rs9939609 (FTO)
Risk allele: A
TT: 36.10TA: 44.10AA: 19.900.420.41p = 0.984
SNPGenotype frequency (%)FAREuropean FAR *Balance H-W
rs1421085 (FTO)
Risk allele: C
TT: 35.20TC: 46.10CC: 18.800.420.42p = 0.372
rs3751723 (IRX3)
Risk allele: T
GG: 41.40GT: 45.70TT: 12.900.360.39p = 0.959
rs4788099 (SH2B1)
Risk allele: G
AA: 43.00AG: 43.70GG: 13.300.350.38p = 0.678
rs6567160 (MC4R)
Risk allele: C
TT: 59.80TC: 35.90CC: 4.300.220.22p = 0.706
rs8092503 (RAB27B)
Risk allele: G
AA: 66.80AG: 28.90GG: 4.300.190.23p = 0.607
rs6857 (NECTIN2)
Risk allele: T
CC: 85.20CT: 14.10TT: 0.800.080.16p = 0.798
rs757318 (CRTC1)
Risk allele: A
CC: 27.30CA: 52.00AA: 20.700.470.48p = 0.656
rs3761445 (PLA2G6)
Risk allele: G
AA: 34.00AG: 52.70GG: 13.300.400.42p = 0.311
SNP: single nucleotide polymorphism; FAR: frequency of risk allele; * FAR Spain: risk allele frequency according to ALFA Project (NCBI). Iberian population reference in Spain [40].

References

  1. Liberali, R.; Kupek, E.; Assis, M.A.A. Dietary Patterns and Childhood Obesity Risk: A Systematic Review. Child Obes. 2020, 16, 70–85. [Google Scholar] [CrossRef]
  2. Vesnina, A.; Prosekov, A.; Kozlova, O.; Atuchin, V. Genes and Eating Preferences, Their Roles in Personalized Nutrition. Genes 2020, 11, 357. [Google Scholar] [CrossRef]
  3. Comuzzie, A.G.; Cole, S.A.; Laston, S.L.; Voruganti, V.S.; Haack, K.; Gibbs, R.A.; Butte, N.F. Novel genetic loci identified for the pathophysiology of childhood obesity in the Hispanic population. PLoS ONE 2012, 12, e51954. [Google Scholar] [CrossRef]
  4. Jääskeläinen, A.; Schwab, U.; Kolehmainen, M.; Kaakinen, M.; Savolainen, M.J.; Froguel, P.; Cauchi, S.; Järvelin, M.-R.; Laitinen, J. Meal frequencies modify the effect of common genetic variants on body mass index in adolescents of the northern Finland birth cohort 1986. PLoS ONE 2013, 8, e73802. [Google Scholar] [CrossRef]
  5. Llewellyn, C.H.; Fildes, A. Behavioral Susceptibility Theory: Professor Jane Wardle and the role of appetite in the genetic risk of obesity. Curr. Obes. Rep. 2017, 6, 38–45. [Google Scholar] [CrossRef]
  6. Steinsbekk, S.; Belsky, D.; Guzey, I.C.; Wardle, J.; Wichstrøm, L. Polygenic Risk, Appetite Traits, and Weight Gain in Middle Childhood: A Longitudinal Study. JAMA Pediatr. 2016, 170, e154472. [Google Scholar] [CrossRef]
  7. Monnereau, C.; Jansen, P.W.; Tiemeier, H.; Jaddoe, V.W.V.; Felix, J.F. Influence of genetic variants associated with body mass index on eating behavior in childhood. Obesity 2017, 25, 765–772. [Google Scholar] [CrossRef]
  8. Faith, M.S.; Carnell, S.; Kral, T.V. Genetics of food intake self-regulation in childhood: Literature review and research opportunities. Hum. Hered. 2013, 75, 80–89. [Google Scholar] [CrossRef]
  9. Hughes, S.O.; Frazier-Wood, A.C. Satiety and the Self-Regulation of Food Take in Children: A Potential Role for Gene-Environment Interplay. Curr. Obes. Rep. 2016, 5, 81–87. [Google Scholar] [CrossRef]
  10. Herle, M.; Smith, A.D.; Kininmonth, A.; Llewellyn, C. The Role of Eating Behaviours in Genetic Susceptibility to Obesity. Curr. Obes. Rep. 2020, 9, 512–521. [Google Scholar] [CrossRef]
  11. Wardle, J.; Carnell, S.; Haworth, C.M.; Farooqi, I.S.; O’Rahilly, S.; Plomin, R. Obesity Associated Genetic Variation in FTO Is Associated with Diminished Satiety. J. Clin. Endocrinol. Metab. 2008, 93, 3640–3643. [Google Scholar] [CrossRef] [PubMed]
  12. Brytek-Matera, A.; Czepczor-Bernat, K.; Olejniczak, D. Food-related behaviours among individuals with overweight/obesity and normal body weight. Nutr. J. 2018, 17, 93. [Google Scholar] [CrossRef] [PubMed]
  13. WMA—World Medical Association. Helsinki Declaration—Ethical Principles for Medical Research Involving Human Subjects. 64ª Asamblea General, Fortaleza, Brasil. 2013. Available online: https://www.wma.net/policies-post/wma-declaration-ofhelsinki-ethical-principles-for-medical-research-involving-human-subjects/ (accessed on 12 June 2023).
  14. Wardle, J.; Guthrie, C.A.; Sanderson, S.; Rapoport, L. Development of the Children’s Eating Behaviour Questionnaire. J. Child Psychol. Psychiatry 2001, 42, 963–970. [Google Scholar] [CrossRef] [PubMed]
  15. Jimeno-Martínez, A.; Maneschy, I.; Moreno, L.A.; Bueno-Lozano, G.; De Miguel-Etayo, P.; Flores-Rojas, K.; Jurado-Castro, J.M.; de Lamas, C.; Vázquez-Cobela, R.; Martinez-Lacruz, R.; et al. Reliability and Validation of the Child Eating Behavior Questionnaire in 3- to 6-Year-Old Spanish Children. Front. Psychol. 2022, 13, 705912. [Google Scholar] [CrossRef]
  16. Zhu, J.; Loos, R.J.F.; Lu, L.; Zong, G.; Gan, W.; Ye, X.; Sun, L.; Li, H.; Lin, X. Associations of genetic risk score with obesity and related traits and the modifying effect of physical activity in a Chinese Han population. PLoS ONE 2014, 9, e91442. [Google Scholar] [CrossRef]
  17. Warrington, N.M.; Howe, L.D.; Paternoster, L.; Kaakinen, M.; Herrala, S.; Huikari, V.; Wu, Y.Y.; Kemp, J.P.; Timpson, N.J.; Pourcain, B.S.; et al. A genome-wide association study of body mass index across early life and childhood. Int. J. Epidemiol. 2015, 44, 700–712. [Google Scholar] [CrossRef]
  18. Martos-Moreno, G.A.; Serra-Juhé, C.; Pérez- Juraod, L.A.; Argente, J. Aspectos genéticos de la obesidad. Rev. Esp. Endocrinol. Pediatr. 2017, 8 (Suppl. S1), 21–32. [Google Scholar] [CrossRef]
  19. BioTools, B.M.; Labs, S.A. Speedtools Tissue DNA Extraction Kit. Kit Para la Extracción y Purificación de ADN Genómico a Partir de Tejidos, Células, Bacterias, Hongos y Fluidos Orgánicos. Instrucciones de Uso. (Ref. 21.136/7). 7ª Ed. 2015. Available online: http://www.biotools.eu/documentospdf/Speedtools_Tissue_DNA_Extraction_Kit.esp.ed07.Mayo15.pdf (accessed on 12 June 2023).
  20. Oeth, P.; Beaulieu, M.; Park, C.; Kosman, D.; del Mistro, G.; van deen Boom, D.; Jurinke, C. Sequenom Application Note. iPLEX® Assay: Increased Pleaxing Efficiency and Flexibility for MassARRAY® System through Single Base Primer Extensión with Mass-Modified Terminators. 2005. Available online: http://www.usc.es/cegen/wp-content/uploads/2016/03/Sequenom-iplex-assay.pdf (accessed on 14 June 2023).
  21. Belsky, D.W.; Moffitt, T.E.; Sugden, K.; Williams, B.; Houts, R.; McCarthy, J.; Caspi, A. Development and evaluation of a genetic risk score for obesity. Biodemogr. Soc. Biol. 2013, 59, 85–100. [Google Scholar] [CrossRef]
  22. Cole, T.J.; Bellizzi, M.C.; Flegal, K.M.; Dietz, W.H. Establishing a standard definition for child overweight and obesity worldwide international survey. BMJ 2000, 320, 1240–1243. [Google Scholar] [CrossRef]
  23. Marrodán, M.D.; Martínez, J.R.; González-Montero, M.; López-Ejeda, N.; Cabañas, M.D.; Prado, C. Precisión diagnóstica del índice cintura-talla para la identificación del sobrepeso y de la obesidad infantil. Med. Clin. 2012, 140, 296–301. [Google Scholar] [CrossRef]
  24. Brook, C.G. Determination of body composition of children from skinfold measurements. Arch. Dis. Child. 1971, 46, 182–184. [Google Scholar] [CrossRef] [PubMed]
  25. Siri, W.E. Body composition from fluid spaces and density; analysis of methods. In Techniques for Measuring Body Composition; Brozek, J., Henschel, A., Eds.; National Academy of Sciences: Washington, DC, USA, 1961; pp. 223–244. [Google Scholar]
  26. Durnin, J.V.; Rahaman, M.M. The assessment of the amount of fat in the human body from measurements of skinfold thickness. Br. J. Nutr. 1967, 21, 681–689. [Google Scholar] [CrossRef] [PubMed]
  27. Serrano, M.M.; Santurino, M.M.; Díaz, J.A.; Soblechero, B.A.; Caballero, P.B.; Hernández, L.D.; Yepes, M.G.; Castelar, J.L.; Rosa, J.R.; de Espinosa, M.G.-M. Diagnosis de la obesidad: Actualización de criterios y su validez clínica y poblacional. An. Pediatría 2006, 65, 5–14. [Google Scholar] [CrossRef] [PubMed]
  28. Cecil, J.; Dalton, M.; Finlayson, G.; Blundell, J.; Hetherington, M.; Palmer, C. Obesity and eating behaviour in children and adolescents: Contribution of common gene polymorphisms. Int. Rev. Psychiatry 2012, 24, 200–210. [Google Scholar] [CrossRef]
  29. Llewellyn, C.H.; Trzaskowski, M.; van Jaarsveld, C.H.M.; Plomin, R.; Wardle, J. Satiety mechanisms in genetic risk of obesity. JAMA Pediatr. 2014, 168, 338–344. [Google Scholar] [CrossRef]
  30. Steinsbekk, S.; Wichstrøm, L. Predictors of Change in BMI from the Age of 4 to 8. J. Pediatr. Psychol. 2015, 40, 1056–1064. [Google Scholar] [CrossRef]
  31. Fildes, A.; Mallan, K.M.; Cooke, L.; van Jaarsveld, C.H.; Llewellyn, C.H.; Fisher, A.; Daniels, L. The relationship between appetite and food preferences in British and Australian children. Int. J. Behav. Nutr. Phys. Act. 2015, 12, 116. [Google Scholar] [CrossRef]
  32. Rudy, E.; Bauer, K.W.; Hughes, S.O.; O’Connor, T.M.; Vollrath, K.; Davey, A.; Correa, N.E.M.; Chen, T.; Fisher, J.O. Interrelationships of child appetite, weight and snacking among Hispanic preschoolers. Pediatr. Obes. 2018, 13, 38–45. [Google Scholar] [CrossRef]
  33. Panichsillaphakit, E.; Chongpison, Y.; Saengpanit, P.; Kwanbunbumpen, T.; Uaariyapanichkul, J.; Chomtho, S.; Pancharoen, C.; Visuthranukul, C. Cuestionario de comportamiento alimentario infantil correlacionado con composiciones corporales de niños y adolescentes tailandeses con obesidad: Un estudio piloto. J. Nutr. Metab. 2021, 2021, 6496134. [Google Scholar] [CrossRef]
  34. Herle, M.; Fildes, A.; Llewellyn, C.H. Emotional eating is learned not inherited in children, regardless of obesity risk. Pediatr. Obes. 2018, 13, 628–631. [Google Scholar] [CrossRef]
  35. Herle, M.; Fildes, A.; Steinsbekk, S.; Rijsdijk, F.; Llewellyn, C.H. Emotional over- and under-eating in early childhood are learned not inherited. Sci. Rep. 2017, 7, 9092. [Google Scholar] [CrossRef] [PubMed]
  36. Gibson, E.L.; Cooke, L. Understanding Food Fussiness and Its Implications for Food Choice, Health, Weight and Interventions in Young Children: The Impact of Professor Jane Wardle. Curr. Obes. Rep. 2017, 6, 46–56. [Google Scholar] [CrossRef] [PubMed]
  37. Jackson, R.; Haszard, J.J.; Morrison, S.; Galland, B.C.; McIntosh, D.; Ward, A.L.; Meredith-Jones, K.A.; Taylor, R.W. Measuring short-term eating behaviour and desire to eat: Validation of the child eating behaviour questionnaire and a computerized ‘desire to eat’ computerized questionnaire. Appetite 2021, 167, 105661. [Google Scholar] [CrossRef]
  38. Tanofsky-Kraff, M.; Han, J.C.; Anandalingam, K.; Shomaker, L.B.; Columbo, K.M.; Wolkoff, L.E.; Kozlosky, M.; Elliott, C.; Ranzenhofer, L.M.; Roza, C.A.; et al. The FTO gene rs9939609 obesity-risk allele and loss of control over eating. Am. J. Clin. Nutr. 2009, 90, 1483–1488. [Google Scholar] [CrossRef] [PubMed]
  39. Housni, F.E.; Sáenz-Pardo-Reyes, E.; López-Espinoza, A.; Moreno, A.G.M.; Galindo, M.d.R.P.; Saucedo, G.V. Efecto de las técnicas y estrategias de modificación de la velocidad al comer sobre la ingesta de alimentos o energía: Revisión sistemática y metaanálisis. Nutr. Hosp. 2021, 38, 631–644. [Google Scholar] [CrossRef]
  40. Phan, L.; Jin, Y.; Zhang, H.; Qiang, W.; Shekhtman, E.; Shao, D.; Revoe, D.; Villamarin, R.; Ivanchenko, E.; Kimura, M.; et al. ALFA: Allele Frequency Aggregator. National Center for Biotechnology Information, U.S. National Library of Medicine, 10 March 2020. Available online: https://www.ncbi.nlm.nih.gov/snp/docs/gsr/alfa/ (accessed on 23 March 2023).
Figure 1. Associations between some independent SNPs and the pro-ingestion subscales of the CEBQ (N= 258). CEBQ-EF: food enjoyment subscale; CEBQ-FR: food response subscale; CEBQ-EOE: emotional overeating subscale; CEBQ-DD: desire to drink subscale. Binary logistic regression analysis was applied. Values are OR (95% CI) for sex- and age-adjusted models. Significant association * p < 0.05. Independent variable: containing or not at least one risk allele of each SNP. Reference group = category having no risk allele of each SNP. Dependent variable: CEBQ score classified from p50.
Figure 1. Associations between some independent SNPs and the pro-ingestion subscales of the CEBQ (N= 258). CEBQ-EF: food enjoyment subscale; CEBQ-FR: food response subscale; CEBQ-EOE: emotional overeating subscale; CEBQ-DD: desire to drink subscale. Binary logistic regression analysis was applied. Values are OR (95% CI) for sex- and age-adjusted models. Significant association * p < 0.05. Independent variable: containing or not at least one risk allele of each SNP. Reference group = category having no risk allele of each SNP. Dependent variable: CEBQ score classified from p50.
Applsci 13 11169 g001
Figure 2. Associations between some independent SNPs and the anti-ingestion subscales of the CEBQ (N= 258). CEBQ-SR: response to satiety subscale; CEBQ-SE: slowness to eat subscale; CEBQ-EUE: emotional undereating subscale; CEBQ-FF: food demand subscale. A binary logistic regression analysis was applied. Values are OR (95% CI) for sex- and age-adjusted models. Significant association * p < 0.05, ** p < 0.01, *** p < 0.0001. Independent variable: containing or not at least one risk allele of each SNP. Reference group = category having no risk allele of each SNP. Dependent variable: CEBQ score classified from p50.
Figure 2. Associations between some independent SNPs and the anti-ingestion subscales of the CEBQ (N= 258). CEBQ-SR: response to satiety subscale; CEBQ-SE: slowness to eat subscale; CEBQ-EUE: emotional undereating subscale; CEBQ-FF: food demand subscale. A binary logistic regression analysis was applied. Values are OR (95% CI) for sex- and age-adjusted models. Significant association * p < 0.05, ** p < 0.01, *** p < 0.0001. Independent variable: containing or not at least one risk allele of each SNP. Reference group = category having no risk allele of each SNP. Dependent variable: CEBQ score classified from p50.
Applsci 13 11169 g002
Table 1. Direct anthropometric measurements of the sample as a function of sexual dimorphism by age group.
Table 1. Direct anthropometric measurements of the sample as a function of sexual dimorphism by age group.
6–10 Years
(Mean ± SD)
11–16 Years
(Mean ± SD)
MaleFemalepMaleFemalep
Weight (kg) 33.54 ± 10.5034.46 ± 8.180.44051.28 ± 13.3050.53 ± 12.520.939
Height (cm)134.42 ± 10.15136.23 ± 9.420.161155.48 ± 11.46155.92 ± 8.300.179
Waist circumference (cm)65.10 ± 10.5463.91 ± 8.120.32375.45 ± 11.3170.61 ± 8.70<0.001 *
Tricipital skinfold (mm)12.64 ± 6.1713.82 ± 4.940.10014.70 ± 9.6015.76 ± 6.120.179
Bicipital skinfold (mm)8.06 ± 4.838.76 ± 3.750.2029.66 ± 5.349.61 ± 4.250.838
Subescapular skinfold (mm)9.14 ± 6.3710.00 ± 4.870.23512.01 ± 7.5012.52 ± 6.500.548
Suprailiac skinfold (mm) 10.57 ± 7.8411.24 ± 6.570.47715.16 ± 8.7714.30 ± 6.920.372
* Significant differences between groups.
Table 2. Prevalence of overweight and adiposity compared between established obesity genetic risk groups based on median genetic risk score (p50 = 22 points).
Table 2. Prevalence of overweight and adiposity compared between established obesity genetic risk groups based on median genetic risk score (p50 = 22 points).
GRS ≤ Q2
(≤22 Points)
GRS ≥ Q3
(≥23 Points)
Comparison
BMI 1Normal weight69.60%74.30%X2 = 6.245
p = 0.040 *
Overweight24.10%14.60%
Obesity6.30%11.10%
WHR 2Normal weight58.90%52.80%X2 = 1.235
p = 0.556
Abdominal overweight14.60%16.70%
Abdominal
obesity
26.60%30.60%
%BF 3Medium (p < 90)54.90%52.10%X2 = 2.679
p = 0.220
High (p90–p97)18.20%20.80%
Very High
(>p97)
26.90%27.10%
Q: quartile; BMI: body mass index; WHI: waist-height index; NS: not significant, according to a Chi-square test (all variables show non-normal distribution): (*) significant differences. 1 Applying the cutoff points for BMI proposed by Cole et al. [21]. 2 Applying the cutoff points for WHR proposed by Marrodán et al. [22]. 3 Applying the cutoff points for fat percentage proposed by Marrodán et al. [24].
Table 3. Comparison of mean scores of CEBQ subscales according to categories of polygenic risk of obesity (GRS), calculated with 32 SNPs.
Table 3. Comparison of mean scores of CEBQ subscales according to categories of polygenic risk of obesity (GRS), calculated with 32 SNPs.
DimensionSubscaleGRS (Q1: ≤19 Points)GRS (Q2: 20–22 Points)GRS (Q3: 23–25 Points)GRS (Q4: >25 Points)Comparison
Pro-intake1. Enjoyment of food (EF) (X ± SD)3.84 ± 0.683.92 ± 0.673.90 ± 0.643.92 ± 0.63F = 0.513
p = 0.916
2. Food responsiveness (FR)
(X ± SD)
2.29 ± 0.912.13 ± 0.922.53 ± 0.902.51 ± 1.21F = 7.118
p = 0.048 *
3. Emotional overeating (EOE)
(X ± SD)
2.19 ± 0.791.91 ± 0.922.30 ± 0.952.18 ± 0.92F = 7.334
p = 0.042 *
4. Desire to drink (DD) (X ± SD)2.44 ± 0.932.18 ± 0.872.45 ± 0.922.40 ± 0.92F = 4.372
p = 0.224
Anti-intake5. Satiety responsiveness (SR) (X ± SD)2.46 ± 0.652.46 ± 0.482.45 ± 0.712.38 ± 0.70F = 0.736
p = 0.865
6. Slowness in eating (SE) (X ± SD)2.53 ± 1.052.49 ± 0.832.55 ± 0.902.28 ± 1.02F = 2.810
p = 0.122
7. Emotional undereating (EUE) (X ± SD)2.50 ± 0.832.65 ± 0.972.19 ± 0.862.41 ± 0.90F = 8.561
p = 0.036 *
8. Food fussiness (FF) (X ± SD)2.63 ± 0.862.76 ± 0.982.87 ± 0.972.55 ± 0.79F = 3.452
p = 0.327
Q: quartile; X: mean; SD: standard deviation. Analysis by Kruskal–Wallis test; (*) significant differences.
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

Calderón García, A.; Pedrero Tomé, R.; Alaminos-Torres, A.; Prado Martínez, C.; Martínez Álvarez, J.R.; López Ejeda, N.; Cabañas Armesilla, M.D.; Marrodán Serrano, M.D. Dietary Behavior of Spanish Schoolchildren in Relation to the Polygenic Risk of Obesity. Appl. Sci. 2023, 13, 11169. https://doi.org/10.3390/app132011169

AMA Style

Calderón García A, Pedrero Tomé R, Alaminos-Torres A, Prado Martínez C, Martínez Álvarez JR, López Ejeda N, Cabañas Armesilla MD, Marrodán Serrano MD. Dietary Behavior of Spanish Schoolchildren in Relation to the Polygenic Risk of Obesity. Applied Sciences. 2023; 13(20):11169. https://doi.org/10.3390/app132011169

Chicago/Turabian Style

Calderón García, Andrea, Roberto Pedrero Tomé, Ana Alaminos-Torres, Consuelo Prado Martínez, Jesús Román Martínez Álvarez, Noemí López Ejeda, María Dolores Cabañas Armesilla, and María Dolores Marrodán Serrano. 2023. "Dietary Behavior of Spanish Schoolchildren in Relation to the Polygenic Risk of Obesity" Applied Sciences 13, no. 20: 11169. https://doi.org/10.3390/app132011169

APA Style

Calderón García, A., Pedrero Tomé, R., Alaminos-Torres, A., Prado Martínez, C., Martínez Álvarez, J. R., López Ejeda, N., Cabañas Armesilla, M. D., & Marrodán Serrano, M. D. (2023). Dietary Behavior of Spanish Schoolchildren in Relation to the Polygenic Risk of Obesity. Applied Sciences, 13(20), 11169. https://doi.org/10.3390/app132011169

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