Next Article in Journal
Maternal Consumption of Non-Staple Food in the First Trimester and Risk of Neural Tube Defects in Offspring
Next Article in Special Issue
Patterns of Food Parenting Practices and Children’s Intake of Energy-Dense Snack Foods
Previous Article in Journal
Ginseng Berry Extract Supplementation Improves Age-Related Decline of Insulin Signaling in Mice
Previous Article in Special Issue
Use of Food Practices by Childcare Staff and the Association with Dietary Intake of Children at Childcare
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 7
Issue 4
10.3390/nu7043054
nutrients-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interactions between Obesity-Related Copy Number Variants and Dietary Behaviors in Childhood Obesity

by
Dandan Zhang
1,2,†,
Zhenli Li
1,2,†,
Hao Wang
1,2,
Min Yang
3,
Li Liang
4,
Junfen Fu
5,
Chunling Wang
4,
Jie Ling
6,
Yan Zhang
6,
Shuai Zhang
1,2,
Yuyang Xu
6,
Yimin Zhu
6,* and
Maode Lai
1,2,*
1
Department of Pathology, Zhejiang University School of Medicine, 866 Yu-hang-tang Road, Hangzhou 310058, Zhejiang, China
2
Key Laboratory of Disease Proteomics of Zhejiang Province, 866 Yu-hang-tang Road, Hangzhou 310058, Zhejiang, China
3
Department of Nutrition, Zhejiang University School of Public Health, 866 Yu-hang-tang Road, Hangzhou 310058, Zhejiang, China
4
Department of Pediatrics, the First Affiliated Hospital of College of Medicine, Zhejiang University, 79 Qing-chun Road, Hangzhou 310003, Zhejiang, China
5
Department of Endocrinology, Children's Hospital of College of Medicine, Zhejiang University, 25 Guang-fu Road, Hangzhou 310003, Zhejiang, China
6
Department of Epidemiology & Biostatistics, Zhejiang University School of Public Health, 866 Yu-hang-tang Road, Hangzhou 310058, Zhejiang, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Nutrients 2015, 7(4), 3054-3066; https://doi.org/10.3390/nu7043054
Submission received: 28 January 2015 / Revised: 3 April 2015 / Accepted: 14 April 2015 / Published: 22 April 2015
(This article belongs to the Special Issue Nutritional Approaches to the Prevention and Management of Childhood Obesity)
Browse Figures
Versions Notes

Abstract

:
Copy number variants (CNVs) have been implicated as an important genetic marker of obesity, and gene-environment interaction has been found to modulate risk of obesity. To evaluate the associations between CNVs and childhood obesity, as well as the interactions between CNVs and dietary behaviors, we recruited 534 obese children and 508 controls from six cities in China and six candidate CNVs were screened through published genome-wide studies (GWAS) on childhood obesity. We found three loci (10q11.22, 4q25 and 11q11) to be significantly associated with obesity after false discovery rate (FDR) correction (all the p ≤ 0.05). Cumulative effect of the three positive loci was measured by the genetic risk score (GRS), showing a significant relationship with the risk of obesity (Ptrend < 0.001). The OR of obesity increased to 21.38 (95% CI = 21.19–21.55) among the 10q11.22 deletion carriers who had meat-based diets, indicating prominent multiplicative interaction (MI) between deletions of 10q11.22 and preference for a meat-based diet. Simultaneous deletions of 5q13.2 and duplications of 6q14.1 had significant MI with a preference for salty foods. Our results suggested that CNVs may contribute to the genetic susceptibility of childhood obesity, and the CNV-diet interactions modulate the risk of obesity.

1. Introduction

Obesity is a concerning public health issue within modern society, and childhood obesity has become a global epidemic [1]. The prevalence of childhood obesity has increased greatly in the past 25 years in China [2,3]. Childhood obesity is considered to be a risk factor for many chronic diseases such as hypertension [4], coronary artery disease [5], insulin resistance [6] and diabetes mellitus [7]. It is also associated with an increased risk of behavioral traits such as attention deficit hyperactivity disorder [8], lower intelligence [9] and various anxiety disorders [10].
Inheritance or genetic variations are crucial contributors to obesity, with high body mass index (BMI) heritability (0.47 to 0.90) estimated from twin studies [11]. Copy number variants (CNVs) are deletions or duplications of DNA segments ranging from one kilobase to several megabases, and are important constituents of genomic structural variants. A set of CNVs has been reported to be associated with early onset obesity by genome-wide analysis [12,13,14,15]. Deletions of 16p11.2 were first reported to be associated with severe early-onset obesity and developmental delay by Bochukova et al. [14]. Later, the results were replicated successfully by a larger case-control study utilizing array comparative hybridization [16]. A meta-analysis of 15 genome-wide analyses discovered that an upstream deletion polymorphism of NEGR1 (neuronal growth regulator 1) was a candidate variant of BMI [17]. Two subsequent genome-wide studies (GWAS) on childhood obesity detected the variant as well, which confirmed the importance of CNVs upstream of NEGR1 in childhood obesity [12,15]. Another important CNV is located at 10q11.22, which encompasses a BMI associated gene, PPYR1 (pancreatic polypeptide receptor 1), which was reported in a Chinese GWAS study [18]. Jarick et al. also detected copy number changes at 10q11.22 in European children [12].
Besides genetic factors, dietary behaviors play a vital role in obesity [19,20]. Recent studies suggested that genetic factors are involved in modulating the response to diet, thus affect susceptibility to obesity [21,22,23]. Hiroi M et al. reported that the PRDX3 gene was associated with obesity in subjects with high dietary fat intake, but not in subjects with low dietary fat intake [21]. An intervention study found that children with MC4R mutations had more difficulty regulating weight after lifestyle interventions compared to those without such mutations [24]. Such evidence provides a hint of gene-diet interactions in obesity, but the interactions between the main genetic factor—CNV—and diet have not been examined before.
The aims of the present study were to examine the associations of candidate obesity-related common CNVs in obese Han Chinese children and further explore the biological interactions of candidate CNVs and dietary behaviors.

2. Experimental Section

2.1. Subjects

The subjects were recruited from a cross-sectional study on metabolic syndrome of children and adolescents throughout six Chinese cities (Beijing, Shanghai, Tianjin, Hangzhou, Chongqing, and Nanjing) in 2010. A total of 1042 unrelated individuals (534 obese children and 508 non-obese samples) aged 7–17 were finally included in the study. Subject recruitment has been described in detail elsewhere [25]. The obese children were defined those with a BMI of greater than the 95th percentile in the same age and gender group, and individuals with a BMI between the 85th and 15th percentiles were defined as controls [26]. Samples were excluded according to the following criteria: (1) currently with or having a history of sever chronic disease or cancer; (2) possessing incomplete clinical information or a history of poor compliance; (3) having singular values in one or more areas of clinical or other study data; and (4) two or more individuals sharing a blood lineage (as evaluated via random selection). Obese cases and non-obese samples were frequency-matched by gender, age and living area. Informed parental consent of enrolled children was obtained for both blood collection and subsequent genetic analysis prior to the procedures. All experimental procedures were approved by the ethics committee of Children's Hospital, School of Medicine, Zhejiang University (Approval Number 2,009,020).

2.2. Anthropometric Measurements and Epidemiologic Investigation

Anthropometric measurements and epidemiologic investigation has been described previously [25]. In brief, anthropometric indices, including height, weight, waist circumference (WC) and hip circumference (HC) were measured. BMI, a person’s weight divided by their height squared, and waist to height ratio (WHtR), WC divided by height, were calculated. Data for dietary behaviors including dietary patterns (meat-based, vegetable-based or balanced diets), preference for salty foods (categorized as like, dislike or no strong preference) and preference for sweet foods (categorized as like, dislike or no strong preference) were collected via face-to-face interview by trained investigators using a standardized questionnaire.

2.3. CNVs Selection

Since rare variants are generally population specific and hard to replicate, we focused on the common CNVs in this study. We screened the published GWAS and aCGH (array-based comparative genomic hybridization studies) on childhood obesity compiled prior to December 2011; three GWASs and two aCGH studies were published [12,13,14,16,27]. However, only Jarick et al. [12] focused on the common variants. Therefore, common CNVs discussed in Jarick et al. were selected for this study [12]. Two CNV regions (11q11 and 1p31.1) suggesting significance in both discovery and replication datasets were included. In addition, we selected loci with adjusted P values of under 0.005 from case-control discovery samples. Thirteen loci were sent to be designed via Multiplex AccuCopy™ Kit (Genesky Biotehnologies Inc., Shanghai, China); they were 11q11, 1p31.1, 4q32.1, 3q21.3, 4q25, 5q13.2, 15q13.2, 10q11.22(1), 2p22.1, 10q11.22(2), 17q21.31, 1q21.1 and 6q14.1. Six loci were successfully designed for the CNV assay in one panel (see Table 1).
Table
Table 1. Basic information of selected candidate CNV loci.
Table 1. Basic information of selected candidate CNV loci.
CNV_IDLocusPosition *Size (bp)Overlap Genes
CNV_111q11chr11: 55,130,596–55,210,16579,569OR4P1P OR4V1P OR4P4 OR4S2 OR4C6
CNV_25q13.2chr5: 68,903,038–70,343,3131,440,275SMN1 SMN2 SERF1A SMA4 SMA5 GUSBP1 LOC643367 LOC643784 LOC653080 LOC653188 GTF2H2B LOC653391 LOC728488 SERF1B LOC728499 LOC728506 LOC728519 LOC728526 LOC728535 LOC728555 LOC100093625 LOC100132218 LOC100133280
CNV_410q11.22(1)chr10: 46,338,178–46,812,351474,173GLUDP2 PPYR1 GPRIN2 SYT15 BMS1P2 LOC642826 LOC643650 ANXA8L1 CTGLF7 LOC728643 LOC728657 LOC100132646 FAM25B LOC100133189
CNV_1010q11.22(2)chr10: 47,011,183–47,145,122133,939LOC340844 LOC728684
CNV_126q14.1chr6: 81,340,896–81,346,2665370None
CNV_134q25chr4: 108,285,188–108,293,2708082None
* Chromosome positions were based on 2006 (NCBI36/hg18); CNV, copy number variation.

2.4. CNVs Genotyping by AccuCopy™ Technique

One and a half milliliters fasting whole blood samples were collected and frozen at −40 °C for subsequent DNA extraction. The TOYOBO MagExtractor Genomic DNA Purification Kit (Toyobo, Osaka, Japan) was used to extract DNA from whole blood. DNA concentration and purity were measured by Nano drop 2000 (Thermo Scientific) and DNA was stored at −80 °C. Qualified DNA was used to detect the copy number changes utilizing AccuCopy® assay, a technology of multiple competitive real-time Polymerase Chain Reaction (PCR) developed by Genesky BioTech (Shanghai, China) [28]. In brief, the CNV detection assay for each sample was carried out in a 20 μL reaction containing 1 × Multiplex PCR Master Mix, 1 × Fluorescence Primer Mix, 1 × Competitive DNA mix, and ~10 ng genomic DNA. The primers for the target genes are listed in Supporting Information Table S1. The PCR program was as follows: 95 °C 10 min; 11 cycles of 94 °C 20 s, 65 °C–0.5 °C/cycle 40 s, 72 °C 1.5 min; 24 cycles of 94 °C 20 s, 59 °C 30 s, 72 °C 1.5 min; 60 min at 60 °C and keep at 4 °C to completion. PCR products were diluted 20-fold and loaded on ABI 3130XL sequencer (Life technologies, Carlsbad, CA, USA) for quantification detection.

2.5. Statistical Analysis

All statistical analyses were performed using SPSS version 19.0 (SPSS Inc., Chicago, IL, USA). Participants with missing CNV genotype data were excluded from the analysis. Normality of distribution of continuous variables was evaluated by Kolmogorov-Smirnov test. The quantitative data were described as means ± standard deviations (SD) for normally distributed variables, and median (lower quartile–upper quartile) for non-normally distributed data. Independent t-test for normally distributed variables, Mann-U Whitney test for non-normally distributed variables and χ2 test for categorical variables were utilized to examine the differences between obese and non-obese subjects. The standardized data of Chinese children was not available, so the WHO 2007 growth reference was used to calculate the age- and sex-specific BMI z-scores in order to exclude the effects of age and sex [29]. Odds ratio, 95% confidence interval and P values of deletions or duplications for the risk of obesity were computed via logistic regression adjusted for gender and age, using a normal CNV state (two copies) as reference. The same analysis was done when using the suggested age and gender specific cut-off of BMI by the International Obesity Task Force (IOTF) to assess obesity status [30]. The p-values were further adjusted by false discovery rate (FDR) using R software (version 2.15.3). Covariance analysis was used to evaluate the BMI z-score and WHtR difference between different copy number states of the positive loci using sex and age as covariates. For multiple comparisons, the least significant difference (LSD) test was used. Genetic risk score (GRS) was calculated as the number of the CNV risk alleles of the three positive loci on obesity, and chi-square trend test was performed to evaluate the dose response association between the GRS and the risk of obesity.
A multiplicative model of interaction was used to assess interactions between candidate CNV loci and dietary behaviors. As the effects ranging from deletions to duplications of the same loci were not necessarily linear, the interactions of different CNV states (deletion versus normal, duplication vs. normal) between dietary behaviors were calculated separately. Next, interaction terms were added to a multivariate logistic regression model to examine effects of MI, adjusted for sex and age. The bootstrap method with 1000 re-sampling was used to estimate the odds ratio and a 95% confidence interval of the logistic regression to reduce small sample bias. AP (attributable proportion due to interaction) and its 95% CI were calculated to evaluate the proportion of the risk of obesity due to interactions between candidate CNV loci and dietary behaviors [31].

3. Results

Five obese and ten non-obese children failed in the genotyping assay and thus were excluded in the analysis. The clinical characteristics of the remaining 529 obese children and 498 non-obese children were summarized in Table 2. The obese group had higher a BMI, BMI z-score, WHtR, WC, and HC compared with the control group (all p < 0.001). Approximately 23.5% of the obese children preferred a meat-based diet compared to 11.0% of the normal weight children (p <0.001). A significantly higher percentage of cases than controls preferred salty foods (25.8% vs. 10.8%, p < 0.001). No significant difference in preference for sweet food was observed. Subjects who preferred salty food and a meat-based diet had a significantly increased risk of obesity.
The association analysis results of the six CNV loci are shown in Table 3. Deletions of three candidate loci (11q11, 10q11.22 and 4q25) demonstrated significant difference in frequencies between cases and controls (p < 0.05). The significance persisted in all the three loci after multiple corrections by FDR. The deletion of 10q11.22 was closely associated with an increased risk of childhood obesity (OR = 2.40, 95% CI = 1.43–4.00; p = 9.1 × 10−4). Similar results were observed for 4q25 (deletion in 4q25: OR = 3.77, 95% CI = 1.52–9.36, p = 0.004). In addition, the three CNV loci showed consistent results when the obese definition of the IOTF was applied: the deletions of 11q11 and 10q11.2 were significantly associated with obesity after FDR correction (both p < 0.05) and 4q25 showed marginal significance, with a p-value of 0.06 (Supporting Information Table S2).
Table
Table 2. Basic characteristics of the study subjects.
Table 2. Basic characteristics of the study subjects.
CharacteristicsObese Children (n = 529)Non-Obese Children (n = 498)p-Value
Age (years)11.7 ± 2.911.7 ± 2.90.85
Gender (M/F)345/184326/1720.93
BMI (Kg m−2)27.3 ± 4.217.6 ± 2.3<0.001
BMI z-score2.71 ± 0.71−0.17 ± 0.88<0.001
WHtR0.54 ± 0.060.41 ± 0.03<0.001
WHR0.89 ± 0.080.81 ± 0.06<0.001
WC (cm)82.8 ± 13.460.4 ± 8.1<0.001
HC (cm)96.0 (82.0–102.0)75.0 (66.0–84.0)<0.001
Food preference (n, %) <0.001
Meat-based diet120 (23.5%)54 (11.0%)
Balanced diet358 (70.1%)392 (79.7%)
Vegetable-based diet33 (6.5%)46 (9.3%)
Salt preference (n, %) <0.001
Like130 (25.8%)53 (10.8%)
No strong preference/dislike374 (74.2%)439 (89.2%)
Sweet taste (n, %) 0.75
Like319 (66.6%)327 (67.6%)
No strong preference/dislike160 (33.4%)157 (32.4%)
HC was expressed as median (lower quartile–upper quartile), and other variables were expressed as means ± SD; p-values were calculated by independent t-test for normally distributed continuous variables; Mann-U Whitney test for non-normally distributed variables and χ2 test for categorical variables comparing obese and non-obese subjects; BMI, body mass index; BMI z-score, BMI for sex and age z-score; WHtR, waist to height ratio; WHR, waist to hip ratio; WC, waist circumference; HC, hip circumference.
Furthermore, the association analysis of the three positive loci showed consistent results on obesity-related index, BMI z-score and WHtR (Table 4). Deletions of 10q11.22 and 4q25 were significantly associated with higher BMI z-score and WHtR (all p < 0.01), and 11q11 showed marginal significance with the two indexes (p = 0.063 and p = 0.014 for BMI z-score and WHtR, respectively).
Table
Table 3. Associations between candidate CNV loci and the risk of obesity.
Table 3. Associations between candidate CNV loci and the risk of obesity.
LocusTypeObese Children (n, %)Non-Obese Children (n, %)OR (95%CI) 1p-Value 2p (FDR) 3
11q11normal218 (41.2)232 (46.6)1
deletion123 (23.3)92 (18.5)1.43 (1.03–1.98)0.030.05
duplication188 (35.5)174 (34.9)1.15 (0.87–1.52)0.320.75
5q13.2normal386 (73.0)369 (74.1)1
deletion68 (12.9)59 (11.8)1.10 (0.76–1.61)0.610.61
duplication75 (14.2)70 (14.1)1.02 (0.72–1.46)0.900.90
10q11.22normal459 (86.8)464 (93.2)1
deletion52 (9.8)22 (4.4)2.40 (1.43–4.00)9.1 × 10−40.006
duplication18 (3.4)12 (2.4)1.51 (0.72–3.18)0.280.75
10q11.22(2)normal490 (92.6)472 (94.8)1
deletion25 (4.7)15 (3.0)1.61 (0.84–3.10)0.150.21
duplication14 (2.6)11 (2.2)1.23 (0.55–2.73)0.620.87
6q14.1normal270 (51.0)257 (51.6)1
deletion202 (38.2)180 (36.1)1.07 (0.82–1.39)0.610.61
duplication56 (10.6)61 (12.2)0.87 (0.58–1.30)0.500.87
4q25normal500 (94.5)487 (97.8)1
deletion23 (4.3)6 (1.2)3.77 (1.52–9.36)0.0040.01
duplication6 (1.1)5 (1.0)1.16 (0.35–3.83)0.810.90
1,2 Odds ratio, 95% confidence interval and p-values of deletions and duplications were calculated by logistic regression adjusted for gender and age, using normal CNV state as reference; 3 p (FDR) values represented the adjusted p-values for false discovery rate. CNV, copy number variation.
Table
Table 4. Associations between the three positive CNV loci and BMI z-score, WHtR.
Table 4. Associations between the three positive CNV loci and BMI z-score, WHtR.
LocusTypeBMI z-Score 1p-Value 3WHtR 2p-Value 4
11q11normal1.25 ± 1.55 0.47 ± 0.08
deletion1.49 ± 1.690.0630.49 ± 0.090.014
duplication1.29 ± 1.720.6490.48 ± 0.080.258
10q11.22normal1.26 ± 1.61 0.47 ± 0.08
deletion1.87 ± 1.680.0020.50 ± 0.080.006
duplication1.75 ± 2.120.0690.51 ± 0.120.022
4q25normal1.29 ± 1.64 0.48 ± 0.08
deletion2.18 ± 1.450.0070.52 ± 0.090.008
duplication1.31 ± 1.740.8520.51 ± 0.080.187
1,2 BMI z-score and WHtR of each CNV state was presented as means ± SD; 3,4 p-Values were calculated by covariance analysis with adjustment for gender and age, comparing the BMI z-score and WHtR of deletion or duplication group with the copy number = 2 group; WHtR, waist to height ratio; CNV, copy number variation; BMI, body mass index.
To analyze the cumulative effect of associated CNVs, all subjects were divided into four groups according to the GRS (GRS = 0, 1, 2, 3). The ORs of each group were calculated using the group with GRS equal to 0 as reference (OR = 1). The ORs of the subjects with GRS equal to 1, 2 and 3 increased to 1.16 (95% CI = 0.86–1.57), 4.13 (95% CI = 1.67–10.23) and 10.45 (95% CI = 1.33–82.25), respectively. Significant dose-response relationship was observed between GRS and the risk of obesity (p for trend = 2.92 × 10−4) (Figure 1). The risk of obesity increased as the GRS increased.
Nutrients 07 03054 g001 1024
Figure 1. Cumulative effects of risk CNV loci on obesity. GRS (X-axis) was the sum of the CNV risk alleles of the three positive loci (11q11, 10q11.22, 4q25). Odds ratio (Y-axis) was calculated by logistic regression adjusted for gender and age, using the group with no risk type of the three loci as reference. The Ptrend was calculated by Chi-square trend test.
Figure 1. Cumulative effects of risk CNV loci on obesity. GRS (X-axis) was the sum of the CNV risk alleles of the three positive loci (11q11, 10q11.22, 4q25). Odds ratio (Y-axis) was calculated by logistic regression adjusted for gender and age, using the group with no risk type of the three loci as reference. The Ptrend was calculated by Chi-square trend test.
Nutrients 07 03054 g001
The duplications of 6q14.1 showed a marginal association with preference for sweet foods. No significant association was observed between other CNV loci and dietary behaviors (Supporting Information Table S3).
Interestingly, three CNV loci (10q11.22, 5q13.2, 6q13.1) showed significant MI with dietary behaviors (Details are shown in Supporting Information Table S4). Strong MI association between deletion of 10q11.22 and preference for meat-based diet was observed. All deleted carriers with meat-based preference had an increased risk of developing obesity, and the estimated risk (OR) for obesity was 21.38-fold (95% CI = 21.19–21.55) compared to reference (11q11.22 diploid case preferred vegetable-based or balanced diets) (Figure 2a). The OR of MI between deletions at 10q11.22 and preference for a meat-based diet was 19.64 (95% CI = 18.88–20.28), estimated by bootstrap.
Meanwhile, the deletion of 10q11.22 also revealed MI with sweet taste preference, with the OR of the interaction being 3.29 (95% CI = 1.10–9.83), and the AP being 0.76 (95% CI = 0.46–1.06) (Figure 2b). MI between 5q13.2 deletions and salty taste preference were also found to have an OR of 9.89 (95% CI =1.23–79.80), and the estimated AP was 0.90 (95% CI = 0.68–1.11) (Figure 2c). Similar results were also found between duplications of 6q14.1 and salty taste (OR for MI was 3.84, 95% CI = 0.97–15.28; AP = 0.69, 95% CI = 0.27–1.10) (Figure 2d).

4. Discussion

In the present study, we comprehensively evaluated 6 obesity-related CNV loci in childhood obesity of Han Chinese. We confirmed that CNVs at 10q11.22, 4q25 and 11q11 were associated with childhood obesity. Furthermore, significant interactions between 10q11.22, 5q13.2 and 6q14.1 and diet preferences were observed.
The deletions of 10q11.22 (CNV_4) demonstrated a highly significant association with the risk of childhood obesity in our study (p = 9.1 × 10−4). One study on Chinese population also detected that deletions at PPYR1 were significantly associated with higher BMI [18], whereas another study found no copy number changes at PPYR1 in 536 Chinese obese individuals and 263 healthy controls [32]. 10q11.22 contains 14 genes, and PPYR1 is the most interesting gene in the region. PPYR1, also termed neuropeptide Y receptor Y4 (NPY4R), is one of the receptors of pancreatic polypeptide (PP). A series of studies have found that PP was associated with reduced food intake and lowered body weight [33,34,35]. Furthermore, it has been confirmed by several studies that PP regulates the intake of food via PPYR1; PPYR1 knockout mice had no response to PP injections [36,37,38].
Nutrients 07 03054 g002 1024
Figure 2. Positive CNV-diet interactions on the risk of obesity. Odds ratio (OR) (95% CI) for the risk of developing obesity were calculated by logistic regression adjusted for age and sex.
Figure 2. Positive CNV-diet interactions on the risk of obesity. Odds ratio (OR) (95% CI) for the risk of developing obesity were calculated by logistic regression adjusted for age and sex.
Nutrients 07 03054 g002
Interestingly, the interaction between deletion of 10q11.22 and dietary behaviors had prominent effects on the risk of obesity, especially the interaction with a preference for meat-based diets. The calculated AP for dietary pattern was almost 1, meaning that approximately 100% of the risk was due to the interaction. Although this might, to some extent, be due to the limited sample size of our study, it indicated strong interactivity between dietary patterns and 10q11.22 in obesity. Such strong interactions might be due to the regulation PP production. It has been well studied that protein and dietary fat can stimulate PP secretion, and glucose can also lead to a minor increase in PP levels [39]. Regularly, the intake of meat may stimulate PP secretion and PP regulating the intake of food via PPYR1. However, PP cannot function normally in subjects with 10q11.22 deletions. Thus, the 10q11.22 deletion carriers with meat-based diets might be more susceptible to developing obesity due to pathologic regulation of food intake. Interestingly, the 10q11.22 deletion carriers with meat dominant diets in the present study were all obese cases. In addition, the deletion of 10q11.22 also had significant multiplicative interaction (MI) with preference for sweet foods. This may be partly due to PP release stimulated by glucose consumption. These results may help explain how genetic factors responded to diet and are intricately involved in the pathophysiologic development of obesity, i.e., how gene-diet interactions operate on a cellular level.
4q25 is a gene desert region, and deletions of CNV 4q25 conferred a highest risk of childhood obesity in our analysis (OR = 3.77, 95% CI = 1.52–9.36, p = 0.004). Little is known about the biological function of this region. There may reside some long-range regulators such as enhancers or repressors. CNVs may exert their effects on the causative gene that the distance can be as far as 1 Mb away [40]. Additionally, a marginal MI was shown between the duplication of 6q14.1 and salty preference. The duplication of 6q14.1 demonstrated no effect on obesity and the risk of obesity increased 7.09 (95% CI = 2.05–24.58) times when interacted with salty preference. Interestingly, Jarick et al. found that the duplication of the region was associated with obesity, whereas little is known about the region itself [12]. These data suggest that duplications at 6q14.1 play a role in the development of children obesity.
Not surprisingly, a strong cumulative association between risk alleles and childhood obesity was observed, with the risk of obesity increasing as individuals harbor greater numbers of CNV alleles at the three positive CNV loci (10q11.22, 4q25 and 11q11). The cumulative effect is highly significant with a p-value for trend of 1.98 × 10−5. We estimate that children with three risk alleles have an odds ratio of 10.45 for developing obesity. It may be possible to use the GRS to assess the risk of childhood obesity, but this strategy will have to be tested in a prospective study before being utilized for any such risk estimation.

5. Conclusions

In conclusion, our results provide support for the role of CNVs in the risk of childhood obesity, which is modified by lifestyle factors. The present study successfully replicated three obesity-related loci 10q11.22, 4q25 and 11q11 in Han Chinese children. A strong cumulative effect of these loci on the risk of obesity was noted. Furthermore, we observed strong MI between CNVs and dietary behaviors. Our findings highlight the importance of the interplay between genotypes and lifestyle factors. The full elucidation of complex interactions between genes and lifestyle factors may help establish public guidelines for the general population and have great implications on public health systems worldwide.

Supplementary Files

Supplementary File 1

Acknowledgments

This study has been supported by National Key Technology R & D Program of China (2012BAI02B03), National Natural Science Foundation of China (81172755). Zhejiang Provincial Program for the Cultivation of High-Level Innovative Health Talents, Fundamental Research Funds for the Central Universities and Program for Zhejiang Leading Team of Science and Technology Innovation (2010R50050) and the 111 Project (B13026).

Author Contributions

ML, YZ, DZ, and ZL conceived and designed the experiments. HW, MY, LL, JF, CW, JL, YX, YZ and SZ performed the experiments. DZ and ZL analyzed data. DZ, ZL, YZ and ML wrote the paper. The authors appreciated the help of Dr. Erik Matro at Sir Run Run Shaw Hospital of Zhejiang University for his editing on this manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Onis, M.; Blossner, M.; Borghi, E. Global prevalence and trends of overweight and obesity among preschool children. Am. J. Clin. Nutr. 2010, 92, 1257–1264. [Google Scholar] [CrossRef] [PubMed]
  2. Ji, C.Y.; Chen, T.J. Empirical changes in the prevalence of overweight and obesity among chinese students from 1985 to 2010 and corresponding preventive strategies. Biomed. Environ. Sci. 2013, 26, 1–12. [Google Scholar] [PubMed]
  3. Dong, B.; Wang, H.J.; Wang, Z.; Liu, J.S.; Ma, J. Trends in blood pressure and body mass index among chinese children and adolescents from 2005 to 2010. Am. J. Hypertens. 2013, 26, 997–1004. [Google Scholar] [CrossRef] [PubMed]
  4. Lu, X.; Shi, P.; Luo, C.Y.; Zhou, Y.F.; Yu, H.T.; Guo, C.Y.; Wu, F. Prevalence of hypertension in overweight and obese children from a large school-based population in shanghai, China. BMC Public Health 2013, 13, 24. [Google Scholar] [CrossRef] [PubMed]
  5. Rocchini, A.P. Childhood obesity and coronary heart disease. N. Engl. J. Med. 2011, 365, 1927–1929. [Google Scholar] [CrossRef] [PubMed]
  6. Chiarelli, F.; Marcovecchio, M.L. Insulin resistance and obesity in childhood. Eur. J. Endocrinol. 2008, 159, S67–S74. [Google Scholar] [CrossRef] [PubMed]
  7. Verbeeten, K.C.; Elks, C.E.; Daneman, D.; Ong, K.K. Association between childhood obesity and subsequent type 1 diabetes: A systematic review and meta-analysis. Diabet. Med. 2011, 28, 10–18. [Google Scholar] [CrossRef] [PubMed]
  8. Yang, R.; Mao, S.; Zhang, S.; Li, R.; Zhao, Z. Prevalence of obesity and overweight among chinese children with attention deficit hyperactivity disorder: A survey in zhejiang province, China. BMC Psychiatry 2013, 13, 133. [Google Scholar] [CrossRef] [PubMed]
  9. Cserjesi, R.; Molnar, D.; Luminet, O.; Lenard, L. Is there any relationship between obesity and mental flexibility in children? Appetite 2007, 49, 675–678. [Google Scholar] [CrossRef] [PubMed]
  10. Anderson, S.E.; Cohen, P.; Naumova, E.N.; Jacques, P.F.; Must, A. Adolescent obesity and risk for subsequent major depressive disorder and anxiety disorder: Prospective evidence. Psychosom. Med. 2007, 69, 740–747. [Google Scholar] [CrossRef] [PubMed]
  11. Elks, C.E.; den Hoed, M.; Zhao, J.H.; Sharp, S.J.; Wareham, N.J.; Loos, R.J.; Ong, K.K. Variability in the heritability of body mass index: A systematic review and meta-regression. Front. Endocrinol. 2012, 3, 29. [Google Scholar] [CrossRef]
  12. Jarick, I.; Vogel, C.I.; Scherag, S.; Schafer, H.; Hebebrand, J.; Hinney, A.; Scherag, A. Novel common copy number variation for early onset extreme obesity on chromosome 11q11 identified by a genome-wide analysis. Hum. Mol. Genet. 2011, 20, 840–852. [Google Scholar] [CrossRef] [PubMed]
  13. Glessner, J.T.; Bradfield, J.P.; Wang, K.; Takahashi, N.; Zhang, H.; Sleiman, P.M.; Mentch, F.D.; Kim, C.E.; Hou, C.; Thomas, K.A.; et al. A genome-wide study reveals copy number variants exclusive to childhood obesity cases. Am. J. Hum. Genet. 2010, 87, 661–666. [Google Scholar] [CrossRef] [PubMed]
  14. Bochukova, E.G.; Huang, N.; Keogh, J.; Henning, E.; Purmann, C.; Blaszczyk, K.; Saeed, S.; Hamilton-Shield, J.; Clayton-Smith, J.; O’Rahilly, S.; et al. Large, rare chromosomal deletions associated with severe early-onset obesity. Nature 2010, 463, 666–670. [Google Scholar] [CrossRef] [PubMed]
  15. Wheeler, E.; Huang, N.; Bochukova, E.G.; Keogh, J.M.; Lindsay, S.; Garg, S.; Henning, E.; Blackburn, H.; Loos, R.J.; Wareham, N.J.; et al. Genome-wide snp and cnv analysis identifies common and low-frequency variants associated with severe early-onset obesity. Nat. Genet. 2013, 45, 513–517. [Google Scholar] [CrossRef] [PubMed]
  16. Bachmann-Gagescu, R.; Mefford, H.C.; Cowan, C.; Glew, G.M.; Hing, A.V.; Wallace, S.; Bader, P.I.; Hamati, A.; Reitnauer, P.J.; Smith, R.; et al. Recurrent 200-kb deletions of 16p11.2 that include the sh2b1 gene are associated with developmental delay and obesity. Genet. Med. 2010, 12, 641–647. [Google Scholar] [CrossRef] [PubMed]
  17. Willer, C.J.; Speliotes, E.K.; Loos, R.J.; Li, S.; Lindgren, C.M.; Heid, I.M.; Berndt, S.I.; Elliott, A.L.; Jackson, A.U.; Lamina, C.; et al. Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat. Genet. 2009, 41, 25–34. [Google Scholar] [CrossRef] [PubMed]
  18. Sha, B.Y.; Yang, T.L.; Zhao, L.J.; Chen, X.D.; Guo, Y.; Chen, Y.; Pan, F.; Zhang, Z.X.; Dong, S.S.; Xu, X.H.; et al. Genome-wide association study suggested copy number variation may be associated with body mass index in the chinese population. J. Hum. Genet. 2009, 54, 199–202. [Google Scholar] [CrossRef] [PubMed]
  19. Kimm, S.Y. The role of dietary fiber in the development and treatment of childhood obesity. Pediatrics 1995, 96, 1010–1014. [Google Scholar] [PubMed]
  20. Gallos, L.K.; Barttfeld, P.; Havlin, S.; Sigman, M.; Makse, H.A. Collective behavior in the spatial spreading of obesity. Sci. Rep. 2012, 2, 454. [Google Scholar] [CrossRef] [PubMed]
  21. Hiroi, M.; Nagahara, Y.; Miyauchi, R.; Misaki, Y.; Goda, T.; Kasezawa, N.; Sasaki, S.; Yamakawa-Kobayashi, K. The combination of genetic variations in the prdx3 gene and dietary fat intake contribute to obesity risk. Obesity 2011, 19, 882–887. [Google Scholar] [CrossRef] [PubMed]
  22. Marti, A.; Martinez-Gonzalez, M.A.; Martinez, J.A. Interaction between genes and lifestyle factors on obesity. Proc. Nutr. Soc. 2008, 67, 1–8. [Google Scholar] [CrossRef] [PubMed]
  23. Basic, M.; Butorac, A.; Landeka Jurcevic, I.; Bacun-Druzina, V. Obesity: Genome and environment interactions. Arh. Hig. Rada Toksikol. 2012, 63, 395–405. [Google Scholar] [CrossRef] [PubMed]
  24. Reinehr, T.; Hebebrand, J.; Friedel, S.; Toschke, A.M.; Brumm, H.; Biebermann, H.; Hinney, A. Lifestyle intervention in obese children with variations in the melanocortin 4 receptor gene. Obesity 2009, 17, 382–389. [Google Scholar] [CrossRef] [PubMed]
  25. Yang, M.; Xu, Y.; Liang, L.; Fu, J.; Xiong, F.; Liu, G.; Gong, C.; Luo, F.; Chen, S.; Xu, C.; et al. The effects of genetic variation in fto rs9939609 on obesity and dietary preferences in Chinese han children and adolescents. PLoS ONE 2014, 9, e104574. [Google Scholar] [CrossRef] [PubMed]
  26. Group of China Obesity Task Force. Body mass index reference norm for screening overweight and obesity in chinese children and adolescents. Zhonghua Liu Xing Bing Xue Za Zhi 2004, 25, 97–102. [Google Scholar]
  27. Yu, Y.; Zhu, H.; Miller, D.T.; Gusella, J.F.; Platt, O.S.; Wu, B.L.; Shen, Y. Age- and gender-dependent obesity in individuals with 16p11.2 deletion. J. Genet. Genomics 2011, 38, 403–409. [Google Scholar] [CrossRef] [PubMed]
  28. Du, R.; Lu, C.; Jiang, Z.; Li, S.; Ma, R.; An, H.; Xu, M.; An, Y.; Xia, Y.; Jin, L.; et al. Efficient typing of copy number variations in a segmental duplication-mediated rearrangement hotspot using multiplex competitive amplification. J. Hum. Genet. 2012, 57, 545–551. [Google Scholar] [CrossRef] [PubMed]
  29. De Onis, M.; Onyango, A.W.; Borghi, E.; Siyam, A.; Nishida, C.; Siekmann, J. Development of a who growth reference for school-aged children and adolescents. Bull. World Health Organ. 2007, 85, 660–667. [Google Scholar] [CrossRef] [PubMed]
  30. 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] [PubMed]
  31. Andersson, T.; Alfredsson, L.; Kallberg, H.; Zdravkovic, S.; Ahlbom, A. Calculating measures of biological interaction. Eur. J. Epidemiol. 2005, 20, 575–579. [Google Scholar] [CrossRef] [PubMed]
  32. Sun, C.; Cao, M.; Shi, J.; Li, L.; Miao, L.; Hong, J.; Cui, B.; Ning, G. Copy number variations of obesity relevant loci associated with body mass index in young chinese. Gene 2013, 516, 198–203. [Google Scholar] [CrossRef] [PubMed]
  33. Asakawa, A.; Inui, A.; Yuzuriha, H.; Ueno, N.; Katsuura, G.; Fujimiya, M.; Fujino, M.A.; Niijima, A.; Meguid, M.M.; Kasuga, M.; et al. Characterization of the effects of pancreatic polypeptide in the regulation of energy balance. Gastroenterology 2003, 124, 1325–1336. [Google Scholar] [CrossRef] [PubMed]
  34. Batterham, R.L.; le Roux, C.W.; Cohen, M.A.; Park, A.J.; Ellis, S.M.; Patterson, M.; Frost, G.S.; Ghatei, M.A.; Bloom, S.R. Pancreatic polypeptide reduces appetite and food intake in humans. J. Clin. Endocrinol. Metab. 2003, 88, 3989–3992. [Google Scholar] [CrossRef] [PubMed]
  35. Ueno, N.; Inui, A.; Iwamoto, M.; Kaga, T.; Asakawa, A.; Okita, M.; Fujimiya, M.; Nakajima, Y.; Ohmoto, Y.; Ohnaka, M.; et al. Decreased food intake and body weight in pancreatic polypeptide-overexpressing mice. Gastroenterology 1999, 117, 1427–1432. [Google Scholar] [CrossRef] [PubMed]
  36. Balasubramaniam, A.; Mullins, D.E.; Lin, S.; Zhai, W.; Tao, Z.; Dhawan, V.C.; Guzzi, M.; Knittel, J.J.; Slack, K.; Herzog, H.; et al. Neuropeptide y (npy) y4 receptor selective agonists based on npy(32–36): Development of an anorectic y4 receptor selective agonist with picomolar affinity. J. Med. Chem. 2006, 49, 2661–2665. [Google Scholar] [CrossRef] [PubMed]
  37. Lin, S.; Shi, Y.C.; Yulyaningsih, E.; Aljanova, A.; Zhang, L.; Macia, L.; Nguyen, A.D.; Lin, E.J.; During, M.J.; Herzog, H.; et al. Critical role of arcuate y4 receptors and the melanocortin system in pancreatic polypeptide-induced reduction in food intake in mice. PLoS ONE 2009, 4, e8488. [Google Scholar] [CrossRef] [PubMed]
  38. Sainsbury, A.; Shi, Y.C.; Zhang, L.; Aljanova, A.; Lin, Z.; Nguyen, A.D.; Herzog, H.; Lin, S. Y4 receptors and pancreatic polypeptide regulate food intake via hypothalamic orexin and brain-derived neurotropic factor dependent pathways. Neuropeptides 2010, 44, 261–268. [Google Scholar] [CrossRef] [PubMed]
  39. Karhunen, L.J.; Juvonen, K.R.; Huotari, A.; Purhonen, A.K.; Herzig, K.H. Effect of protein, fat, carbohydrate and fibre on gastrointestinal peptide release in humans. Regul. Pept. 2008, 149, 70–78. [Google Scholar] [CrossRef] [PubMed]
  40. Klopocki, E.; Mundlos, S. Copy-number variations, noncoding sequences, and human phenotypes. Annu. Rev. Genomics Hum. Genet. 2011, 12, 53–72. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Zhang, D.; Li, Z.; Wang, H.; Yang, M.; Liang, L.; Fu, J.; Wang, C.; Ling, J.; Zhang, Y.; Zhang, S.; et al. Interactions between Obesity-Related Copy Number Variants and Dietary Behaviors in Childhood Obesity. Nutrients 2015, 7, 3054-3066. https://doi.org/10.3390/nu7043054

AMA Style

Zhang D, Li Z, Wang H, Yang M, Liang L, Fu J, Wang C, Ling J, Zhang Y, Zhang S, et al. Interactions between Obesity-Related Copy Number Variants and Dietary Behaviors in Childhood Obesity. Nutrients. 2015; 7(4):3054-3066. https://doi.org/10.3390/nu7043054

Chicago/Turabian Style

Zhang, Dandan, Zhenli Li, Hao Wang, Min Yang, Li Liang, Junfen Fu, Chunling Wang, Jie Ling, Yan Zhang, Shuai Zhang, and et al. 2015. "Interactions between Obesity-Related Copy Number Variants and Dietary Behaviors in Childhood Obesity" Nutrients 7, no. 4: 3054-3066. https://doi.org/10.3390/nu7043054

Article Metrics

Back to TopTop