| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
National Heart, Lung, and Blood Institutes Framingham Heart Study (C.S.F., N.L.H.-C., R.S.V., J.M.M., R.B.D., L.D.A.), Framingham, Massachusetts 01702; National Heart, Lung, and Blood Institute (C.S.F.), Bethesda, Maryland 20892; Department of Endocrinology, Hypertension, and Metabolism, Brigham and Womens Hospital and Harvard Medical School (C.S.F.), Boston, Massachusetts 02115; and Department of Neurology (N.L.H.-C., L.D.A.), Division of General Internal Medicine (J.M.M.), Division of Cardiology (R.S.V.), and Department of Mathematics and Statistics (R.B.D.), Boston University School of Medicine, Boston, Massachusetts 02118
Address all correspondence and requests for reprints to: Dr. Caroline S. Fox, M.D., MPH, 73 Mount Wayte Avenue, Suite 2, Framingham, Massachusetts 01702. E-mail: foxca{at}nhlbi.nih.gov.
| Abstract |
|---|
|
|
|---|
Methods: Variance components linkage analysis using microsatellites was performed on 336 families from the Framingham Heart Study offspring cohort, using a 10-cM genome-wide linkage analysis. We evaluated linkage to two traits: short-term (8-yr) weight change and long-term (up to 24-yr) weight change. Models were adjusted for age, age squared, baseline weight, smoking status, and menopausal status.
Results: Mean short-term weight change ranged from 1.43.8 kg, and mean long-term weight change was 7.7 kg. The heritability of long-term weight change was 0.24; weight change was minimally heritable among younger individuals and over shorter follow-up intervals. We found significant evidence for linkage for long-term weight change, with a peak LOD score of 3.10 on chromosome 20 at 63.7 cM (nearest marker, D20S481). We also found suggestive evidence for linkage on chromosome 1 at 239.7 cM (LOD score, 2.28; nearest marker, D1S1644).
Conclusion: Long-term weight change is heritable, and evidence for linkage exists on chromosomes 1 and 20. Potential candidate genes include MC3R, ASIP, AGT, and HSD11B1. Additional research is necessary to uncover the genetic underpinnings of weight change that might contribute to associated adverse metabolic profiles.
| Introduction |
|---|
|
|
|---|
Risk factor profiles can be modified by minimal gains and losses in weight (8, 9), suggesting that overall measures of adiposity may not accurately reflect the dynamic interaction of weight change on metabolic profiles. Reductions in weight of 510% can result in significant improvements in lipoprotein profiles (10), insulin sensitivity (11), and blood pressure (12), whereas weight gain of more than 10 kg is associated with an increased risk for hypertension, elevated low density lipoprotein cholesterol, and higher triglycerides (13). Although we have previously reported quantitative trait loci (QTLs) for body mass index (BMI) and waist circumference to 6q2324 and 11q14 in the Framingham Heart Study (14, 15), it is possible that different loci underlie the genetics involved with weight change over time.
In fact, changes in body mass and adiposity have been shown to be heritable, with estimates ranging from 0.230.86 for changes in BMI (16, 17), 0.70 for weight gain (18), and 0.45 for change in waist circumference (16). However, limited data exist looking at linkage to weight change. Thus, we determined heritability and linkage to weight change in the Framingham Heart Study offspring study.
| Subjects and Methods |
|---|
|
|
|---|
Details regarding the methods of risk factor measurement and laboratory analysis have been described previously (21). Each examination included anthropometry, an extensive cardiovascular disease assessment, 12-lead electrocardiogram, and blood testing. Weight was measured to the nearest pound with the participant wearing only a gown without shoes or slippers and standing in the middle of the scale with weight equally distributed on both feet. For this analysis, we only used information regarding age, sex, weight, smoking status (defined as smoker at baseline or follow-up), and menopausal status (pre- or postmenopausal). BMI was calculated by dividing weight (in kilograms) by height (in meters squared). Adjustments were made for age and age squared to be consistent with our previously published work (14). We also adjusted for menopausal status, because the menopausal transition is associated with increases in body weight (22), increased visceral (23) and central fat depots (24), and higher BMI (24) compared with premenopausal women. We adjusted for smoking status, because smoking cessation has been associated with body weight gain (25). We adjusted for baseline weight to investigate genetic influences on weight change independent of baseline weight, because the magnitude of absolute weight change depends on baseline weight.
Leukocyte DNA was extracted from 510 ml whole blood or buffy coat specimens using a standard protocol (26). Aliquots of DNA from members of the largest Framingham Heart Study families were sent in four batches to the Mammalian Genotyping Service Laboratory at Marshfield Clinic (Marshfield, WI), consisting of 336 pedigrees. A 10-cM density genome-wide scan was performed (average heterozygosity, 0.77). Genotype data cleaning, including verification of family relationships and Mendelian inconsistencies, have been previously described (27).
For the purpose of this analysis, we created three datasets using information from participant examinations that occurred at roughly 8-yr intervals (19711975 to 19791982, 19791982 to 19871992, and 19871992 to 19951998), enabling the creation of a weight change variable over 8 yr (short-term weight change), and a fourth dataset that examined long-term weight change (19711975 to 19951998) over approximately 24 yr. Participants were excluded from these analyses if they were under age 20 yr or had a BMI less than 18 kg/m2 at any of the baseline examinations. Severe weight cyclers, defined as participants who repeatedly lost and gained 9 kg in at least three examination cycles, were also excluded from analyses. The rationale for this particular exclusion is derived from data from prior observational studies that suggest that individuals with weight fluctuation are at increased risk for all-cause mortality (28), indicating that it may be a different phenotype from more typical weight change over time. Because linkage signals are often the strongest when phenotypes are parsed into their most basic components, we removed the weight cyclers from our analysis. We also performed a secondary analysis in which we limited our study sample to participants under the age of 65 yr. Although obesity increases with age in young and middle-aged adults, declines are seen after age 60 yr (29), which may confound the assessment of the genetic component.
Statistical analysis
Using SAS (30), sex-specific residuals were computed using weight change, as defined by follow-up weight minus baseline weight. Sex-specific residuals were used to account for gender differences in adiposity. Residuals included adjustment for age, age squared, baseline weight, smoking status, and menopausal status (in women). Residuals were obtained from a regression-based linear model and were used to predict values of weight change after correction for variation attributable to covariates. The residuals were computed by calculating the difference between the actual and predicted values obtained by linear regression. In this way the correction for covariates was fixed across the dataset. The use of residuals was necessary to standardize values between different examination cycles and genders as well as to account for the effects of covariates. A residual value greater than 4 SD from the mean was reduced to the 4 SD value (i.e. Winsorized) to reduce the skewness and kurtosis of the data to an acceptable level (31). Specifically, the sex-specific residuals were derived for each of the four datasets separately (resulting in eight regression equations). Residuals were derived separately in men and women, then merged. The regression equations in women for weight change residual were: exams 1 to 6, 35.1 + (0.19 x age) (0.01 x age2) (0.03 x wt) + (1.94 x smok_baseline) (7.15 x smok_follow) + (0.06 x menopausal status); exams 1 to 2, 0.23 + (age x 0.53) (0.01 x age2) (0.01 x wt) + (0.54 x smok_baseline) (2.89 x smok_follow) + (0.33 x menopausal); exams 2 to 4, 9.6 + (age x 0.91) (0.01 x age 2) + (0.02 x wt) + (4.84 x smok_baseline) (5.48 x smok_follow) (2.8 x menopausal); exams 4 to 6, 40.5 (age x 0.53) (0.0003 x age2) (0.05 x wt) + (4.28 x smok_baseline) (6.06 x smok_follow) + (2.01 x menopausal), where wt stands for weight, smok_baseline is smoking status at baseline, smok_follow is smoking status at follow-up, and menopausal is menopausal status. The regression equations among men for weight change residual were: exams 1 to 6, 45.6 (0.52 x age) (0.004 x age2) (0.06 x wt) + (4.70 x smok_baseline) (4.84 x smok_follow); exams 1 to 2, 25.5 (0.37 x age) + (0.002 x age2) (0.08 x wt) + (1.50 x smok_baseline) (1.34 x smok_follow). exams 2 to 4, 23.9-(0.21 x age) (0.001 x age2) (0.05 x wt) + (5.03 x smok_baseline) (4.22 x smok_follow); and exams 4 to 6, 43.5 (1.23 x age) + (0.01 x age2) + (0.001 x wt) + (8.22 x smok_baseline) (9.44 x smok_follow). R2 values ranged from 0.020.13. For women, we included menopausal status as a covariate. Otherwise, models were identical for all datasets. Residuals were computed and then merged together for use in Genehunter to compute the heritability and linkage estimates.
Variance components linkage analysis (32, 33), as implemented in Genehunter (34), was used to perform the genomewide linkage analysis using residuals for weight change as the phenotype of interest. This approach uses the genotype information at a locus to decompose the phenotypic variance into a component attributable to the locus [known as a quantitative trait locus (QTL)], a polygenic component and an environmental component. Genotype information at a locus is characterized by the probability that two related individuals share zero, one, or two alleles identical by descent (IBD). IBD computation grows exponentially with pedigree size in Genehunter. To make the IBD computation feasible, 19 of the largest families were split into smaller pedigrees. The resulting small loss of genetic information should be conservative with respect to linkage results.
All variance components were estimated by maximum likelihood. Linkage was tested by a likelihood ratio test in which the hypothesis that the QTL variance component is equal to zero was compared with it being greater than zero. The resulting
2 statistic was converted to a traditional LOD score by dividing by 2 x ln (10). It is possible to obtain estimates of the proportion of the total phenotypic variation due to the QTL. However, two recent studies (35, 36) have shown that the estimate of this effect is strongly correlated with the LOD score estimate and is poorly correlated with the true proportion. Therefore, we will not present effect estimates here. The Framingham Heart Study is population based and is not selected for any particular trait; therefore, no ascertainment correction is necessary.
| Results |
|---|
|
|
|---|
|
|
|
|
| Discussion |
|---|
|
|
|---|
Previous heritability studies have demonstrated that changes in measures of adiposity over time, including BMI, weight, and waist circumference, are heritable (16, 17, 18). These previous estimates were higher than ours, with estimates of heritability as high as 0.70 for weight gain.
Limited data exist on linkage to weight gain in the published literature. As part of the Genetic Analysis Workshop proceedings, linkage analysis was performed on the slope of weight change in the Framingham Heart Study, which demonstrated a peak multipoint LOD score of 1.6 on chromosome 8 at 152 cM for body weight change and 1.9 on chromosome 8 at 102 cM for body weight gain and body weight change up to 50 yr of age (37). Different methodologies, including use of the regression slope instead of actual weight change, may contribute to differences from our findings.
The QTL that we identified on chromosome 20q13 has been replicated in several other linkage analyses that have examined BMI, percent body fat, and eating behavior, with LOD scores ranging from 2.23.2 (38). This region harbors several interesting candidate genes, including the melanocortin 3 receptor, MC3R, which has been shown to be associated with BMI (39) and severe obesity (40); the agouti-signaling protein, ASIP; GHRH, which has been associated with central adiposity (41); and the adenosine deaminase gene, ADA, which has been associated with BMI (42).
We also found suggestive linkage on chromosome 1 to weight change, with a LOD score of 2.28. Interestingly, this QTL is the same location where long-term fasting glucose values were mapped from the Framingham Heart Study (43), raising the possibility that weight change may be an important genetic component in the development of diabetes. Furthermore, the angiotensinogen (AGT) gene lies in this region, and polymorphisms in this gene have been associated with weight gain among treated hypertensive patients (44). This region also contains the 11ß-hydroxysteroid dehydrogenase 1 gene (HSD11B1), which has been shown to be associated with obesity in expression studies (45).
Candidate gene analyses have been performed looking at weight change. Polymorphisms in the leptin receptor and the ß2-adrenergic receptor were shown to be associated with weight gain in a case-control study comparing participants who gained an average of 12.8 kg over 7 yr with those who remained weight stable (46). Polymorphisms in the serotonin 5-HT2c receptor gene have been associated with clozapine-induced weight gain among patients with schizophrenia (47). Candidate gene studies of overfeeding also shed light on genetic susceptibility to weight gain. In 12 pairs of identical male twins, polymorphisms in the UCP gene were involved in recovery from overfeeding (48), and polymorphisms in the ß2-adrenergic receptor (49) and the resistin gene (50) were shown to be associated with depot-specific fat accumulation after overfeeding. These data suggest that weight gain is an important phenotype that warrants additional study in the general population.
Strengths of our study include repeated measures of weight collected over a 27-yr period using standardized methods and well documented covariate data. Some limitations do exist. The predominantly Caucasian sample that comprises the majority of the Framingham offspring cohort may limit the generality of our findings. However, cardiovascular disease risk factor relationships from Framingham have been validated in six ethnically and geographically diverse cohorts and were found to be applicable in other populations, reinforcing the generalizability of our data (51). Our genetic data were not collected until 1995; therefore, there may have been a survival bias, because subjects with a greater genetic predisposition to more severe weight change may have died prematurely. Our peak LOD score of 3.1 is significant for genome-wide significance, but there is a possibility that our results may represent a false positive finding. Our sample size and the marker spacing of approximately 10 cM in the genome-wide linkage analysis may limit our power to localize disease-influencing QTL in linkage to weight change. Although there may be additional QTLs not detected in our study, simulation studies using finer marker maps have demonstrated that the location error is low in maps using a marker spacing of 10 cM compared with 0.5 cM (52), suggesting that this should not have affected our results significantly. Lastly, we were unable to adjust for physical activity and diet in our multivariable regressions.
In summary, we have found evidence for linkage to long-term weight gain to two genomic regions. Both of these regions have been implicated in independent studies of obesity-related phenotypes. Additional studies are warranted to understand the genetic underpinnings to weight change, which affects blood pressure, lipids, and glycemia.
| Footnotes |
|---|
First Published Online March 15, 2005
Abbreviations: BMI, Body mass index; IBD, identical by descent; QTL, quantitative trait locus.
Received September 7, 2004.
Accepted March 7, 2005.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |