|
|
||||||||
ORIGINAL RESEARCH COMMUNICATION |
1 From the University of Maryland School of Medicine, Baltimore (NIS, SS, TIP, BDM, and ARS); the Veterans Administration Geriatric Research and Education Clinical Center, Baltimore (ARS); the Southwest Foundation for Biomedical Research, San Antonio, TX (W-CH); Axys Pharmaceuticals, La Jolla, CA (HS and CJB); and Glaxo Wellcome Inc, Research Triangle Park, NC (PLSJ).
2 Supported by grants from Glaxo Wellcome Inc, Sequana Therapeutics, the NIH (R01 DK54261 and K24 DK02673), the American Diabetes Association, and the Hopkins Bayview General Clinical Research Center.
3 Reprints not available. Address correspondence to AR Shuldiner, Division of Endocrinology, Diabetes and Nutrition, University of Maryland School of Medicine, 660 West Redmond Street, Room 494, Baltimore, MD 21201. E-mail: ashuldin{at}medicine.umaryland.edu.
| ABSTRACT |
|---|
|
|
|---|
Objective: We investigated the genetic underpinnings of eating behavior.
Design: We administered a standardized eating behavior inventory to 624 adults from 28 families participating in the Amish Family Diabetes Study. Three quantifiable components of eating behavior were measured: restraint, disinhibition, and hunger. Associations between eating behavior scores and physical characteristics were evaluated. Heritability analysis and a genome-wide multipoint linkage analysis were performed.
Results: Eating behavior scores were associated with obesity and obesity-related phenotypes. Heritability estimates were 0.28 ± 0.09 for restraint, 0.40 ± 0.10 for disinhibition, and 0.23 ± 0.09 for hunger (P < 0.001). The linkage analysis showed 4 regions of suggestive linkage. We observed suggestive evidence for linkage of restraint scores to 2 chromosomal regions, near markers D3S1304 [LOD (log of odds) = 2.5, P = 0.0003] and D6S276 (LOD = 2.3, P = 0.0006). We previously reported that D3S1304 is linked to a locus influencing percentage body fat in this same population (LOD = 1.6), suggesting that this behavioral phenotype may be secondary to obesity. The maximum LOD scores for disinhibition were 1.6 (P = 0.003) near marker D7S657 and 1.4 (P = 0.005) near marker D16S752. The maximum LOD score for hunger was 1.4 (P = 0.005) near marker D3S1278.
Conclusion: Significant familial effects on eating behavior and suggestive genetic linkage were found in Amish adults.
Key Words: Restraint obesity behavioral genetics Amish heritability linkage analysis
| INTRODUCTION |
|---|
|
|
|---|
A few rare syndromes involving obesity and genetic mutations have been identified in humans (5). For example, Prader-Willi syndrome, an autosomal dominant disorder caused by a deletion in chromosome region 15q11.2q12, is characterized by excessive eating and morbid obesity. Very rare monogenic forms of human obesity are due to mutations in the genes for leptin (6, 7), proconvertase 1 (8), leptin receptor (9), proopiomelanocortin (10), or melanocortin 4 receptor (11, 12), all of which result in a phenotype of excessive energy intake relative to expenditure. However, gene variants that increase the risk of typical obesity in humans remain largely unknown. Understanding the genetic basis for the control of eating behavior has implications for developing models of the pathophysiology of nutrition-related diseases and for prevention, diagnosis, and treatment of these disorders.
Eating behavior is difficult to define and quantify. The Three-Factor Eating Questionnaire (13, 14) measures 3 behavioral factors: restraint, disinhibition, and hunger. Restraint is a cognitive avoidance of eating to control body weight. An example of a true-or-false question that addresses restraint is "I do not eat some foods because they make me fat." Disinhibition is loss of restraint resulting in overeating. An example of a disinhibition question is "When I am with someone who overeats I usually overeat too." Hunger measures the perceived need for food, eg, "I get so hungry that my stomach often feels like a bottomless pit."
To examine the possible relations between eating behavior, obesity and related traits, and genetics, we administered the Three-Factor Eating Questionnaire to participants in the Amish Family Diabetes Study, a large family study investigating the genetics of type 2 diabetes and obesity. First, we evaluated the association between eating behavior and obesity phenotypes. Then we estimated the heritability of eating behavior. Finally, we conducted a genome-wide linkage analysis to identify chromosomal regions that are likely to harbor genes that regulate eating behavior traits.
| SUBJECTS AND METHODS |
|---|
|
|
|---|
Subject recruitment for the Amish Family Diabetes Study began in early 1995. Details of the study design, recruitment, procedures, and pedigree construction were previously described (17). Subjects were enrolled through the Amish Research Clinic in Strasburg, PA. Individuals with type 2 diabetes were identified through door-to-door interviews and by word of mouth. Probands were defined as individuals with previously diagnosed diabetes with an age at diagnosis of 3565 y. All first- and second-degree family members aged
18 y were recruited around the probands. If another individual with diabetes was identified in the family (eg, aunt or uncle), the family was expanded further to include the first- and second-degree relatives aged
18 y of that individual.
Because dietary education for management of diabetes may alter eating behavior scores, subjects with prevalent diabetes (n = 30) were excluded from these analyses. Thus, this report is based on data from 624 relatives of the probands who were not known to have diabetes at the time of their evaluation. Although all of our subjects can be related by tracing their ancestors back multiple generations (18), to reduce computational difficulties, we divided the sample into 28 discrete families, ranging in size from 3 to 69 individuals. The sample included a large number of relative pairs, including 436 parent-offspring pairs, 1326 sibling pairs, 1342 avuncular pairs (aunt or uncle with a niece or nephew), and 1311 first cousin pairs. The study was conducted in accordance with the standards involving human subjects of the Institutional Review Board of the University of Maryland, and informed consent was obtained from all subjects.
Phenotypic characterization
Data regarding medical and family history were obtained during one clinic visit. Medication history was reviewed. No subject was taking any medication known to affect appetite. Height and weight were measured with the use of a stadiometer and a calibrated scale. The subjects removed their shoes and wore light clothing for these measurements. Waist circumference was measured at the level of the umbilicus, and hip circumference was measured at the widest protuberance across the pelvis. Body mass index (BMI) was calculated as weight (kg) divided by height squared (m2) to estimate total adiposity. The waist-to-hip ratio (WHR) and subscapular-to-triceps skinfold thickness ratio were calculated as indexes of abdominal (central) adipose distribution. Percentage body fat was determined by bioimpedance with the use of standard formulas (17).
Fasting lipid profile (total cholesterol, HDL cholesterol, and triacylglycerols) was assayed by Quest diagnostics (Baltimore), and the interassay CVs were 1.6% for total cholesterol, 5.0% for HDL cholesterol, and 1.6% for triacylglycerols. Fasting leptin concentrations were assayed at the Johns Hopkins Bayview General Clinical Research Center core laboratory by radioimmunoassay (Linco, Inc, St Louis, MO; interassay CV: 4.25%). Oral glucose tolerance tests were administered to subjects with no history of diabetes.
The Three-Factor Eating Questionnaire was administered at the same visit. This questionnaire has 51 items, 36 true-false and the remainder asking for a rating on a scale of 14 or 15 (13, 14). Each question is assigned a weight of one point, and scores are categorized by eating behavior trait. The range of scores for restraint is 021, for disinhibition 016, and for hunger 014.
To our knowledge, the Three-Factor Eating Questionnaire has not been used previously to study the Amish. Therefore, we investigated the means and distributions of answers to each question as well as the overall scores for restraint, disinhibition, and hunger. We found that the responses to each question were distributed across the range of possible responses in the Amish. Furthermore, eating behavior scores spanned all possible values for each of the 3 behaviors, and the configuration of frequency of scores was compatible with a Gaussian distribution. Finally, we compared the mean scores from subjects in other studies with our results and found that the mean scores of Amish subjects were similar to those reported in other populations (14).
Analysis of polymorphic short tandem repeat markers
Whole blood was collected, and genomic DNA was isolated with the use of Qiagen Maxipreps and the manufacturer's instructions (Qiagen, Inc, Santa Clarita, CA). For our genome scan we used fluorescent semiautomated genotype analysis with the Prism Linkage Mapping Set (Applied Biosystems Division, Perkin-Elmer, Foster City, CA). Polymerase chain reactions were performed on 25 ng genomic DNA with the use of 96-well plates. Polymerase chain reaction mixtures were pooled according to the manufacturer's recommendations and loaded onto a 5% denaturing polyacrylamide gel (Long Ranger; FMC Bioproducts, Rockland, ME) positioned within the ABI 377 automated DNA sequencer (Applied Biosystems Division, Perkin-Elmer). Polymerase chain reaction failure rates were low (generally <1%). The gel runs were visually inspected and were analyzed with the use of ABI GENESCAN software (Applied Biosystems). The inspected data were imported into ABI GENOTYPER software (Applied Biosystems), and alleles were assigned and binned. To assess genotyping accuracy, several pairs of replicate DNA samples were genotyped. On the basis of blind replicates of the 373 markers (357 autosomal and 16 X chromosome) that were used in the linkage analysis, the genotyping error rate was 0.16%. Mendel error checks were performed with the use of MENDEL (19); we detected no cases of nonpaternity. The small number of genotypes exhibiting Mendel errors that could not be resolved by visual inspection of the gel were archived and deleted from the database before linkage analysis. On average, we obtained usable genotypes from 96.3 ± 2.4% (
± SD; range: 81.999.8%) for each of the 373 markers used in the linkage analysis. The marker order and distances between markers were estimated from our data with the use of CRI-MAP (20). The average interval between markers was 10 centimorgans (cM), and the largest gap between markers was 25.4 cM, occurring on chromosome 7. The mean marker heterozygosity was 0.75 (range: 0.330.91).
Statistical methods
Association of eating questionnaire scores with clinical characteristics
The clinical characteristics of the study subjects were analyzed according to sex. The statistical significance of differences between men and women was tested with the use of the analysis of variance procedure in SYSTAT 8.0 (SPSS Inc, Chicago), adjusting for the effects of age and BMI. Two-factor analysis of variance was used to determine whether eating behavior scores were associated with sex or BMI.
We next conducted more complex analyses that took into account the complicated pedigree structure to estimate the effect of each eating behavior score on each measure of obesity. A series of models were considered, each including a single obesity measure as the dependent variable and a single eating behavior score as the independent variable. We then used maximum likelihood methods to estimate the association between eating behavior score and obesity measure while simultaneously adjusting for the effects of age and sex. Analyses were conditioned on the pedigree structures to account for the nonindependence between related individuals. We tested the statistical significance of the association by likelihood ratio test; that is, we compared the likelihood of the pedigree data under a full model (allowing the effect of eating behavior on obesity to be estimated) with that under a nested model (constraining the effect of eating behavior to be zero). Two times the difference between the log-likelihoods of the 2 models is distributed asymptotically as a chi-square distribution with degrees of freedom equal to the difference in the number of dependent factors in the models being compared. In our analysis, the chi-square tests had one degree of freedom. These tests were performed with the use of the SOLAR program (21).
Heritability and linkage analyses
Heritabilities of the eating behavior scores were assessed as the proportion of total phenotypic variance that could be accounted for by the covariance among relatives, while simultaneously adjusting for the effects of age, sex-specific age, and age2. We used maximum likelihood procedures as implemented in the SOLAR software package (21) for these analyses. We also considered the possibility that the familial effects on eating behavior may be due to shared household rather than shared genes. We therefore estimated the household effects on eating behavior by assigning a unique household number to individuals sharing the same street address at the time of their examination. The household effect was parameterized as a random effect by assigning a 1 if the relative pair shared the same household and a 0 if they did not.
We carried out quantitative trait locus (QTL) linkage analysis with the use of a variance components method in which variation in eating behavior scores was partitioned into components attributable to environmental covariates, the additive effects of genes (ie, residual heritability), a specific QTL (ie, the linkage component), and random errors. For these analyses, also, we used maximum likelihood procedures as implemented in the SOLAR software package (21). The additive genetic effect was modeled as a function of the expected genetic covariances between relatives, and the QTL effect was modeled as a function of the identity-by-descent relations at the marker locus. The hypothesis of linkage was evaluated by the likelihood ratio test, in which we evaluated whether the locus-specific effect was significantly greater than zero (ie, Ho:
2QTL = 0 compared with HA:
2QTL > 0). Both multipoint and two-point linkage analyses were carried out. Scores for restraint, disinhibition, and hunger were analyzed separately, and in each analysis, we simultaneously adjusted for the effects of age, sex-specific age, and age2 on the outcome score. Because nonnormality of the data distribution can artificially inflate nominal LOD (log of odds) scores, we derived the distribution of LOD scores under the null hypothesis of no linkage by Monte Carlo simulation (22). Briefly, we simulated an unlinked marker locus with 5 equifrequent alleles, assigned genotypes to each founder, and then dropped genotypes down through the pedigree on the basis of Mendelian expectations and the founder genotypes. We then conducted linkage analysis of each eating behavior score with the simulated genotypes. We conducted 5000 replicates and defined the probability of obtaining a false positive result as the proportion of replicates for which we obtained a specified LOD score or higher. These probabilities were converted back into LOD scores, and these LOD scores can be interpreted as the probabilities of detecting linkage to an unlinked marker.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
19 cM from the p terminus near marker D3S1304. The 2-point LOD scores for markers in this region were 2.29 (D3S1304 at position 16.7 cM), 1.04 (D3S1263 at position 26.7 cM), and 2.64 (D3S3608 at position 31.0 cM). We also detected suggestive evidence for linkage of restraint scores to a region on chromosome 6,
40 cM from the p terminus near marker D6S276 (LOD = 2.3). The 2-point LOD scores for markers in this region were 0.91 (D6S470 at position 22.6 cM), 1.12 (D6S289 at position 32.9 cM), and 1.98 (D6S276 at position 46.5 cM). Evidence for linkage to the other eating behavior traits was less impressive. For disinhibition scores, the maximum multipoint LOD score occurred on chromosome 7 near marker D7S657 (LOD = 1.6) and on chromosome 16 near marker D16S752 (LOD = 1.4). For hunger scores, the maximum multipoint LOD score was 1.4, occurring on chromosome 3 near marker D3S1278. | DISCUSSION |
|---|
|
|
|---|
As in a Swedish study by Bjorvell et al (28), disinhibition scores in the present study were higher in obese persons than in normal-weight subjects. Because the brain influences food intake, it is probable that certain inherent behavioral traits, such as a tendency for disinhibition, may promote obesity in the appropriate environment. However, some studies showed changes in the scores for restraint, disinhibition, and hunger in response to treatment for obesity (14, 24), suggesting that eating behavior is pliable to some extent. Some studies also found a negative correlation between restraint and mean daily energy intake (14, 24), suggesting that restraint may serve to counteract weight gain. However, personality traits that favor dietary restraint may unintentionally promote weight gain by precipitating disinhibition and overeating (14). In the present study, we found that obesity was associated with high restraint scores. Design features may have improved our ability to detect a relation that has eluded other investigators. It is likely that the influence of body composition on eating behavior is culturally determined and that a reduced stigmatization of obese individuals among the Old Order Amish may diminish the influence of obesity on behavior and thus make the influence of inherent behavioral traits on body composition more apparent. A large number of subjects is another advantage of the present study. To avoid the confounding effect of professional dietary advice, subjects with prevalent diabetes (n = 30) were excluded. In contrast with other studies, ours did not select subjects on the basis of their interest in a weight-loss program.
The evidence for a genetic component to eating behaviors comes from our heritability and linkage analyses. Heritability was highest for disinhibition, with the magnitude of the additive genetic effects rivaling that of many nonbehavioral measures of obesity (29). Disinhibition was also the factor that showed the strongest association with obesity phenotypes. To the extent that related persons may also live together, it could be argued that a component of the heritability between related individuals is attributable to their shared household. With the exception of disinhibition, we found that the contribution of household effect to our heritability estimates of eating behavior traits was negligible, supporting the theory that biological rather than social inheritance explains the observations. Although subjects with diabetes were specifically excluded because of potential effects of nutrition counseling on eating behaviors, it is possible that households with diabetic subjects may exhibit a familial disinhibition that results from a shared effect of the nutrition counseling provided to the diabetic family member.
Linkage of eating behavior traits to specific chromosomal regions further suggests that specific genetic variants may account for some of these genetic effects. However, these results must be interpreted cautiously. Evidence for linkage was expressed in terms of LOD scores, defined as the logarithm of the ratio of the likelihood for linkage to the likelihood for no linkage. Some experts have proposed that when many genetic markers are analyzed simultaneously (as in genome-wide searches), LOD-score thresholds of 3.3 and 1.9 be used to assert evidence for significant linkage and suggestive linkage, respectively (30). It is becoming increasingly clear that for complex traits, these criteria are likely to be too rigorous. Importantly, the LOD scores that we report were derived by simulation studies; thus, they correspond to specific empirically derived P values that are not likely to be inflated. The relevance of these linkages must await replication by others or identification of the actual gene variants accounting for phenotypic variation in eating behavior traits.
The burgeoning Human Genome Project database has provided the opportunity to identify many of the genes that reside within the putative regions of linkage. Indeed, several candidate genes in the regions of these linkages may influence eating behavior. However, because we do not yet know all of the genes expressed within these regions, nor do we know the functions of many of the genes that reside in these regions, examination of positional candidate genes must be regarded with caution. A functional candidate gene on chromosome 3, which is located within our region of peak linkage to disinhibition at 10.3810.5 cM, encodes peroxisome proliferator activated receptor
(PPAR-
). PPAR-
is expressed in many tissues, but the highest concentrations are found in adipose tissue and large intestine (31). PPAR-
is an important mediator of adipocyte differentiation and energy storage. PPAR-
agonists stimulate eating by decreasing leptin concentrations (32). Furthermore, it was shown that mutations in PPAR-
are associated with obesity (33, 34). Indeed, we found modest evidence for linkage of this same region on chromosome 3 to obesity-related traits in the Amish (29). We also found a linkage for restraint on chromosome 6. The area of linkage lies in proximity to the genes encoding glucagon-like peptide 1 (GLP-1) receptor, tumor necrosis factor
(TNF-
), and lymphotoxin. GLP-1 was shown to be involved in eating behavior in human and animal models (35, 36). GLP-1 receptors are found in the hypothalamus and are the target for hypothalamic neurons containing glucose transporter 2 and glucokinase. GLP-1 secretion after a meal stimulates insulin release, lowers blood glucose, and reduces food intake (36). TNF-
and lymphotoxin are cytokines known to be expressed in adipocytes. TNF-
may have a role in the development of insulin resistance and obesity (3739). Evidence of association and linkage of the TNF-
gene with percentage body fat and BMI was reported by others (5). Linkage analysis of disinhibition scores yielded suggestive evidence for linkage on chromosome 7. Interestingly, the region of linkage lies in proximity to the leptin and plasminogen activator inhibitor 1 (PAI-1) genes. Other researchers found evidence of association of leptin with weight loss, body weight, leptin concentrations, and leptin concentration in response to diet (5, 3943). PAI-1 plays a role in the development of coronary artery disease. Plasma concentrations of PAI-1 may be higher in persons with obesity, insulin resistance, and type 2 diabetes and in their family members (45, 46).
To our knowledge this is the largest family study of eating behavior published to date. We found that eating behavior is familial in the Amish and for the first time show linkage of specific chromosomal regions to eating behavior traits. Replication of these studies in other populations and ultimately the positional cloning of the genes involved will be required to confirm and extend our findings. Insights gained by identifying genes that influence eating behavior may provide new molecular targets for the early detection of individuals susceptible to obesity, which may advance the development of preventive interventions and provide new drug targets for the treatment of obesity and related disorders.
| ACKNOWLEDGMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. V. E. Kral and M. S. Faith Influences on Child Eating and Weight Development from a Behavioral Genetics Perspective J. Pediatr. Psychol., April 10, 2008; (2008) jsn037v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. de Krom, Y. T. van der Schouw, J. Hendriks, R. A. Ophoff, C. H. van Gils, R. P. Stolk, D. E. Grobbee, and R. Adan Common Genetic Variations in CCK, Leptin, and Leptin Receptor Genes Are Associated With Specific Human Eating Patterns Diabetes, January 1, 2007; 56(1): 276 - 280. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. HAINER, K. KABRNOVA, B. ALDHOON, M. KUNESOVA, and M. WAGENKNECHT Serotonin and Norepinephrine Reuptake Inhibition and Eating Behavior Ann. N.Y. Acad. Sci., November 1, 2006; 1083(1): 252 - 269. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. I. Steinle, T. I. Pollin, J. R. O'Connell, B. D. Mitchell, and A. R. Shuldiner Variants in the Ghrelin Gene Are Associated with Metabolic Syndrome in the Old Order Amish J. Clin. Endocrinol. Metab., December 1, 2005; 90(12): 6672 - 6677. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Tholin, F. Rasmussen, P. Tynelius, and J. Karlsson Genetic and environmental influences on eating behavior: the Swedish Young Male Twins Study Am. J. Clinical Nutrition, March 1, 2005; 81(3): 564 - 569. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. I. Pollin, K. Tanner, J. R. O'Connell, S. H. Ott, C. M. Damcott, A. R. Shuldiner, J. C. McLenithan, and B. D. Mitchell Linkage of Plasma Adiponectin Levels to 3q27 Explained by Association With Variation in the APM1 Gene Diabetes, January 1, 2005; 54(1): 268 - 274. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Bouchard, V. Drapeau, V. Provencher, S. Lemieux, Y. Chagnon, T. Rice, D. Rao, M.-C. Vohl, A. Tremblay, C. Bouchard, et al. Neuromedin {beta}: a strong candidate gene linking eating behaviors and susceptibility to obesity Am. J. Clinical Nutrition, December 1, 2004; 80(6): 1478 - 1486. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Cai, S. A Cole, R. A Bastarrachea-Sosa, J. W MacCluer, J. Blangero, and A. G Comuzzie Quantitative trait locus determining dietary macronutrient intakes is located on human chromosome 2p22 Am. J. Clinical Nutrition, November 1, 2004; 80(5): 1410 - 1414. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S Faith, K. L Keller, S. L Johnson, A. Pietrobelli, P. E Matz, S. Must, M. A. Jorge, J. Cooperberg, S. B Heymsfield, and D. B Allison Familial aggregation of energy intake in children Am. J. Clinical Nutrition, May 1, 2004; 79(5): 844 - 850. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Collaku, T. Rankinen, T. Rice, A. S Leon, D. Rao, J. S Skinner, J. H Wilmore, and C. Bouchard A genome-wide linkage scan for dietary energy and nutrient intakes: the Health, Risk Factors, Exercise Training, and Genetics (HERITAGE) Family Study Am. J. Clinical Nutrition, May 1, 2004; 79(5): 881 - 886. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |