Quantitative trait linkage analysis of longitudinal change in body weight
© Golla et al; licensee BioMed Central Ltd 2003
Published: 31 December 2003
One of the great strengths of the Framingham Heart Study data, provided for the Genetic Analysis Workshop 13, is the long-term survey of phenotypic data. We used this unique data to create new phenotypes representing the pattern of longitudinal change of the provided phenotypes, especially systolic blood pressure and body weight. We performed a linear regression of body weight and systolic blood pressure on age and took the slopes as new phenotypes for quantitative trait linkage analysis using the SOLAR package. There was no evidence for heritability of systolic blood pressure change. Heritability was estimated as 0.15 for adult life "body weight change", measured as the regression slope, and "body weight gain" (including only individuals with a positive regression slope), and as 0.22 for body weight "change up to 50" (regression slope of weight on age up to an age of 50). With multipoint analysis, two regions on the long arm of chromosome 8 showed the highest LOD scores of 1.6 at 152 cM for "body weight change" and of >1.9 around location 102 cM for "body weight gain" and "change up to 50". The latter two LOD scores almost reach the threshold for suggestive linkage. We conclude that the chromosome 8 region may harbor a gene acting on long-term body weight regulation, thereby contributing to the development of the metabolic syndrome.
One of the special qualities of the Framingham Heart Study data, provided by the Genetic Analysis Workshop 13 (GAW13), is the long-term survey of phenotypic data. The quantitative phenotypes systolic blood pressure, body height, body weight, fasting serum glucose, total plasma cholesterol, high density lipoprotein cholesterol, and triglycerides were measured repeatedly over a period of 40 years in the first generation of study participants and over a period of 20 years in their offspring. We used this unique data to create new phenotypes representing the pattern of longitudinal change of the provided phenotypes, assuming that not only the measured traits, but also their pattern of longitudinal change, could be genetically determined. We chose the phenotypes systolic blood pressure and body weight because of their known tendency to increase during adult lifespan, and because of their most complete longitudinal survey by repeated measurements at each examination during the study. Our own previous analyses of the population sample of Framingham Cohort 1 and 2 showed that body weight and systolic blood pressure rise linearly with age, but the gradient is steeper when the regression is performed only up to an age of 50 compared with the regression on all available measurements. Another group participating in Genetic Analysis Workshop 13 (GAW13) analyzed the changes of body mass index with age (BMI = weight/height2, which corresponds to weight change, because there is almost no change of height during adult life). They found a linear gain-phase up to a maximum value (at about age 53 on average) and a decline-phase later in life . According to this we analyzed "body weight change" during whole life, "body weight gain" (using phenotypes of individuals with a positive gain of weight only), and body weight "change up to 50", the regression slope of weight on age up to an age of 50.
Materials and Methods
The Framingham Heart Study data were provided as Problem 1 of GAW13. The data consisted of family structure information, genotypic information (genome scan), and long-term phenotypic information for 1702 subjects from 330 families. Phenotypic information was provided over a period of 40 years for the parent generation and 20 years for the offspring generation. Phenotypic information consisted of values of repeated measurements of systolic blood pressure, height, weight, lipids (total plasma cholesterol, high density lipoprotein cholesterol, triglycerides), and fasting plasma glucose. We chose the phenotypes systolic blood pressure and body weight for our analysis.
Linear regression was perfomed using SAS  for the values of systolic blood pressure and body weight, at each time of examination, on age, for each individual separately. To obtain the values for adult life only, individuals younger than 18 years at the first examination were excluded. We omitted the values of these individuals instead of using the examinations with age 18 as a starting point because of lower reliability of the regression slopes (based on four time points at most and often less due to younger age and missing values) and because of their relatively low contribution due to small frequency (only 6% of individuals of Cohort 2 and none of Cohort 1 were younger than 18 years at their first examination). For systolic blood pressure, measurements were omitted if high blood pressure treatment was stated at this examination or any previous examination. The slope estimates were obtained for all individuals and used as a new phenotype providing a measure of the adult life change of systolic blood pressure and body weight. Further analyses of body weight change were conducted with regression slope up to an age of 50 as the phenotype "change up to 50" and with only those individuals who showed a positive weight slope (≥ 0.01), representing the phenotype adult life "body weight gain". Heritability estimates for the new phenotypes were obtained using variance-component methodology as implemented in the SOLAR package . Heritability estimates were done analyzing the covariates sex, age at first examination and weight at first examination in addition. Only age at first examination contributed to the variance of the three phenotypes, explaining about 12% of the variance of "body weight change", 6% of the variance of "body weight gain", and 9% of the variance of "change up to 50". Because further examinations were at the same regular intervals for all individuals of the two cohorts, age at first examination represents age-dependence during the whole study. Two-point and multipoint quantitative trait linkage analyses were conducted on the regression slope estimates for adult life "body weight change", "body weight gain", and "change up to 50", using the SOLAR package. In the case of "body weight gain", we included age at first examination as a covariate. The information content of chromosomal regions with positive LOD scores was obtained using GENEHUNTER .
Adult life increase of systolic blood pressure and of body weight are known phenomena, at least in western lifestyle populations. Environmental factors like lack of exercise and high food intake are major factors influencing both of them. Yet little is known as to whether heritable factors may contribute to the pattern of change of blood pressure and body weight during lifetime. Our analysis using regression slopes as a measure for lifetime changes showed no evidence of a heritable basis of blood pressure change, but we found heritabilities of 0.15 for "body weight change" and "body weight gain" and of 0.22 for body weight "change up to 50". We identified two regions 50 cM apart from one another on the long arm of chromosome 8 with linkage to "body weight change" (distal region) and to "body weight gain" and "change up to 50" (proximal region). Although not significant, the LOD scores for "body weight gain" of 1.97 and "change up to 50" of 1.93 at locations 102 and 103 cM on chromosome 8 almost reach the threshold of 2.2 for suggestive evidence of linkage . Very large sample sizes are required for detection of a quantitative trait locus (QTL) attributable to 10% or 15% heritability of a trait with significant LOD scores. For instance, if 10% heritability is due to a single QTL, 3300–33,000 individuals (depending on family structure) are required to obtain a LOD score of 3.0 with 80% power . The Framingham Heart Study sample with 1702 subjects from 330 families is very big, yet it may well be not big enough to obtain statistically significant linkage results for traits with total heritabilities of 0.15 to 0.22, even if a single locus would attribute to 0.67 or 0.50 of the heritability. We therefore think, that "almost suggestive linkage" is probably the best obtainable result for the body change traits we analyzed, and that confirmation would require analysis (or meta-analysis) of still larger sample sizes.
Interestingly, two other groups participating in the Derived Phenotypes Group of GAW13 reported related results with different (but pathophysiologically related) phenotypes. Jun et al. found statistically significant evidence of linkage between the region on chromosome 8 around marker GATA21C12 and fasting glucose levels progression (exactly the same region which we identified for "body weight change") . Wilcox et al. derived four phenotypically distinct groups using correspondence analysis and clustering. One of these groups resembled atherogenic dyslipidemia and the metabolic syndrome. Besides other chromosomal regions, they found the same chromosome 8 location (two-point LOD scores above two for markers GATA7G07, UT721, GATA21C12) for their metabolic-syndrome-like phenotype as we found for "body weight change" . The region on chromosome 3, where we obtained a relatively high two-point LOD score for "body weight change", contains AOMS1 (abdominal obesity-metabolic syndrome QTL1), the first known metabolic syndrome QTL . These findings underscore the pathophysiological relationship between body weight change (body weight gain), blood glucose level progression, and the metabolic syndrome, and point to a long-term body weight regulating gene on chromosome 8 that may play a role in the development of the metabolic syndrome.
Longitudinal data like that collected in the Framingham Heart Study provided for GAW13 enable the search not only for genes influencing the phenotypes themselves, but also for genes contributing to the pattern of change during lifetime. Using regression slope estimates as a parameter for longitudinal change of body weight, we found evidence for a long-term body weight regulating gene located on the long arm of chromosome 8. Combining our results with the findings of other GAW13 participants [7, 8], we think that this gene might be involved in the development of the metabolic syndrome by influencing long-term body weight and consequently blood glucose level. Location of genes involved in long-term regulation of physiological processes can be found using phenotypes derived from longitudinal data, and exactly those long-term regulatory genes may play key roles in frequent diseases of late onset.
- Strug L, Sun L, Corey M: Genetics of longitudinal and cross-sectional BMI. BMC Genetics. 2003, 4 (suppl 1): S14-10.1186/1471-2156-4-S1-S14.PubMed CentralView ArticlePubMedGoogle Scholar
- SAS Institute Inc.: SAS/STAT User's Guide, Version 8. Cary, NC, SAS Institute Inc. 1999Google Scholar
- Almasy L, Blangero J: Multipoint quantitative-trait linkage analysis in general pedigrees. Am J Hum Genet. 1998, 62: 1198-1211. 10.1086/301844.PubMed CentralView ArticlePubMedGoogle Scholar
- Krugylak L, Daly MJ, Reeve-Daly MP, Lander ES: Parametric and non-parametric linkage analysis: A unified multipoint approach. Am J Hum Genet. 1996, 58: 1347-1363.Google Scholar
- Lander ES, Krugylak L: Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nat Genet. 1995, 11: 241-247. 10.1038/ng1195-241.View ArticlePubMedGoogle Scholar
- Williams JT, Blangero J: Power of variance component linkage analysis to detect quantitative trait loci. Ann Hum Genet. 1999, 63: 545-563. 10.1046/j.1469-1809.1999.6360545.x.View ArticlePubMedGoogle Scholar
- Jun G, Song Y, Stein CM, Iyengar SK: An autosome-wide search using longitudinal data for loci linked to type 2 diabetes progression. BMC Genetics. 2003, 4 (suppl 1): S8-10.1186/1471-2156-4-S1-S8.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilcox MA, Wyszynski D, Panhuysen CI, Ma Q, Yip A, Farrell J, Farrer LA: Empirically derived phenotypic subgroups: qualitative and quantitative trait analyses. BMC Genetics. 2003, 4 (suppl 1): S15-10.1186/1471-2156-4-S1-S15.PubMed CentralView ArticlePubMedGoogle Scholar
- Kissebah AH, Sonnenberg GE, Myklebust J, Goldstein M, Broman K, James RG, Marks JA, Krakower GR, Jacob HJ, Weber J, Martin L, Blangero J, Comuzzie AG: Quantitative trait loci on chromosome 3 and 17 influence phenotypes of the metabolic syndrome. Proc Natl Acad Sci USA. 2000, 19: 14478-14483. 10.1073/pnas.97.26.14478.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.