| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Sackler Faculty of Medicine (G.L.), Tel Aviv University, Tel Aviv 69978, Israel; Twin Research and Genetic Epidemiology (G.L., B.S.K., T.D.S.), Unit St. Thomas Hospital, Kings College London, London SE1 7EH, United Kingdom; and Department of Endocrinology and Diabetes (S.G.W.), Sir Charles Gairdner Hospital, Nedlands, Western Australia 6009
Address all correspondence and requests for reprints to: Tim D. Spector, Twin Research and Genetic Epidemiology Unit, St. Thomas Hospital, Kings College London, London SE1 7EH, United Kingdom. E-mail: tim.spector{at}kcl.ac.uk.
| Abstract |
|---|
|
|
|---|
Objective: Our objective was to examine the contribution of genetic factors for LEAN-tot variation, including a genome-wide search for the genes.
Research Methods: Dual-energy x-ray absorptiometry measurements of LEAN-tot were obtained from each of the 3180 United Kingdom females (509 monozygotic and 1081 dizygotic twin pairs). Contribution of genetic factors was assessed using variance component analysis. A genome-wide linkage analysis was performed on the dizygotic twins using a modified version of the Haseman-Elston method.
Results: Age, body height, total fat, and bone mass were correlated with LEAN-tot, and commonly explained 52% of the LEAN-tot variation. The crude heritability estimate was 74.0 ± 4.0%, after adjustment for the aforementioned factors; 65.2 ± 4.6% was attributable to independent genetic effects. Significant (P < 0.001) genetic correlations were found between LEAN-tot and bone mass, and LEAN-tot and total fat. Adjusted only for age, LEAN-tot showed no significant linkage. After adjustment for all covariates, significant linkage (LOD = 4.49 and 3.62) was observed at chromosome 12q24.3 and 14q22.3, respectively. Additional peaks of interest were on 7p15.3-15.1 (LOD = 2.86) and 8p22 (LOD = 2.83).
Conclusions: LEAN-tot measured by dual-energy x-ray absorptiometry is highly heritable, independent of other body measures. This first genomic search for genes associated with the lean component of body mass suggests significant linkage to quantitative trait loci on chromosomes 12 and 14.
| Introduction |
|---|
|
|
|---|
Heritability, candidate genes and genome-wide linkage analyses for bone mass and fat mass have been widely reported in numerous publications (11, 12). Surprisingly, however, there are only very few studies that examined familial resemblance of lean mass (13, 14, 15). Moreover, we are not aware of any genome-wide linkage or association study of lean mass. This may have added importance because substantial common genetic effects are shared with bone mass and structure (8, 15).
Dual-energy x-ray absorptiometry (DXA) provides a definitive method of accurately determining these three body composition components, and allows investigation of their independent effects on various aspects of human physiology in normalcy and in pathology (16). The present study used DXA body composition data in a very large population-based sample of female twins and provides the first quantitative trait loci (QTL) from a genome-wide linkage scan of total lean body mass (LEAN-tot).
| Subjects and Methods |
|---|
|
|
|---|
The data examined in the present study were from the TwinsUK Adult Twin Registry, described in detail elsewhere (17). All participants gave written informed consent before entering the study, and the St. Thomas Hospital research ethics committee approved the project. The present study comprised 1018 monozygotic (MZ) and 2162 dizygotic (DZ) twins having available whole body densitometry data for body components. All participants were female volunteers ascertained from the general population through national media campaigns in the United Kingdom, and unselected for any disease or trait.
Phenotype
Each individual in the sample was assessed for body weight and height in addition to body composition measurements. Body composition was measured using the whole body DXA method, following manufacturers recommendations (QDR 4500W system; Hologic, Inc., Bedford, MA). All scan printouts were reviewed by an expert reader to ensure proper positioning and analysis. The subject lay supine on a bed and was scanned from head to toe for determination of total bone mineral density (BMD-tot), total lean tissue (LEAN-tot) mass, and total fat (FAT-tot) mass, as described by us elsewhere (13). Note that although LEAN-tot includes mostly muscle mass, it also includes tendons and soft tissue, such as viscera.
Genotyping
DNA was extracted using the BACC2 DNA extraction kit (Nucleon Biosciences, Coatbridge, Lanarkshire, UK) from 10 ml venous blood collected into 1.6 mg/ml EDTA. Genome-wide linkage analysis included the genotyping of 737 highly polymorphic markers, using standard fluorescence-based genotyping methodologies (18) from the ABI Prism linkage mapping set (Applied Biosystems, Foster City, CA), and ordered according to the combined linkage-physical map of the human genome (19). Consistencies among genotypes, family relationships, and twin zygosity were routinely investigated, and this included analysis of patterns of Mendelian inheritance and identity-by-state relationships. Random duplicate genotyping was routinely undertaken throughout the study and indicated a mean genotyping error rate of less than 1% for the microsatellite genotyping. Discrepant genotypes and MZ twins identified by this routine were excluded from further analyses.
Statistical analysis
First, all MZ and DZ twin study participants were included in analysis of distribution and adjustment, using the multiple linear regression technique (StataCorp, version 9; StataCorp LP, College Station, TX). Three types of adjustment were undertaken: 1) age and age-derivatives were used for adjustment of variation of LEAN-tot to obtain the crude heritability and linkage signal estimates of the trait; 2) to examine the correlation between the three body composition components, free of body size effect, each of them was adjusted for age and body height; and 3) LEAN-tot was adjusted for age, height, FAT-tot, and BMD-tot to test the residuals created at this stage for lean-specific heritability and linkage.
To estimate heritability, we used quantitative-model fitting to twin data. Briefly, this approach is based on comparisons of the correlations between MZ and DZ twin pairs, and allows for the partitioning of observed phenotypical variance into additive genetic component, and common and/or unique environmental components (20). In this study the variance component estimates were obtained by the maximum-likelihood based method implemented in the MAN-7 software package (21).
Because our data showed significant correlation between the three body composition components, we also conducted a bivariate genetic analysis between each pair of these variables (21). This analysis represents an extension from univariate to multivariate models, and allows exploration of the question as to whether the origin of the covariance between the different variables is genetically and/or environmentally determined. The general model estimates contribution of the common genetic factors, genetic correlation (rG) and shared environment, environmental correlation (rE) between two traits. The statistical significance of rG and rE was studied using nested models, and examining the change in
2 values between the models. A more detailed description is given elsewhere (20). We used this analysis to ascertain the extent to which encountered phenotypical correlations are attributable to common (pleiotropic) genetic and shared environmental effects. To achieve this, we estimated rG and rE between the BMD-tot and FAT-tot, both adjusted for age and body height, as well as between the BMD-tot and LEAN-tot adjusted for the same covariates.
Linkage analysis
Only DZ twins were examined at this stage. A total of 2162 individuals with LEAN-tot were genotyped and available for the analysis. Twin pairs formed independent families. The estimation of the proportion of alleles shared identical-by-descent (IBD) by a twin pair was obtained using GENEHUNTER 2 software (Whitehead Institute, MIT, Cambridge, MA). Multipoint genome-wide linkage analysis was conducted using generalized linear modeling in Stata (StataCorp LP). This regression technique is based on optimal Haseman-Elston methods (22), in which the square of the sibling difference in LEAN-tot phenotype (adjusted for age and other covariates) is regressed on the estimated proportion of alleles IBD. This method is algebraically equivalent to other likelihood techniques (23, 24). The prioritized regions were selected if the LOD score obtained was suggestive, i.e. more than 1.9 (23). The confirmation of the positive results was made by computing the genome-wide significance level (empirical P value) using a permutation approach (25). To obtain empirical P value estimates of the established LOD scores of interest, 1000 permutations of the data sets were performed for each LOD score, keeping both IBD structure and family structure intact. To keep twin pairs reserved, we computed the squared difference for each twin pair, which implies for 1081 DZ twin pairs we obtain 1081 "squared differences," i.e. we used the same phenotype as our phenotype for the genome scan. Phenotype values were permuted among individuals between the twin pairs. To obtain additional confirmatory evidence for the presence of the potential QTL in the chromosomal region(s) showing significant LOD scores, we used a bootstrap approach (26). Basically, the bootstrap is a procedure that involves choosing random samples with replacement from a data set and analyzing each sample the same way. Similarly to permutations, the bootstrap samples were sampled from the dependent variable consisting of 1081 squared differences. In this study the random samples of the same size as the original data set were drawn 2000 times with replacement. In each bootstrap replication b (b = 1, ... , 2000), we conducted a genome scan and each time estimated corresponding LOD scores.
It should be mentioned that concerning the X-chromosome, female twins share at least one paternal allele; the oversharing on X-chromosome may lead to a different distribution of the test statistic compared with autosomes. Because "Stata" package does not correct for the aforementioned bias and because of the poor coverage on X-chromosome, we will not consider X-chromosome associated linkage peaks in this study.
| Results |
|---|
|
|
|---|
Table 1
provides the basic descriptive statistics of the studied traits for MZ and DZ twins separately. All measurements in both samples were very similar. LEAN-tot significantly correlated (P < 0.001) with age, anthropometrics, and FAT-tot and BMD. Table 2
summarizes the results of multiple regression analyses of LEAN-tot on all the potential covariates in the total sample. As seen, body height, FAT-tot, BMD-tot, and age were retained in the regression equation, collectively explaining some 51.5% of the LEAN-tot total variation. The distribution of lean data adjusted for age, age and height, and for all significant covariates closely followed the normality assumptions by the Kolmogorov-Smirnov test (d = 0.018; P > 0.10).
|
|
|
Our bivariate genetic analysis included simultaneous estimates of variance components and adjustment for significant covariates (adjustment No. 2). Both genetic and environmental correlations between LEAN-tot and BMD-tot were statistically highly significant by likelihood ratio tests (rG = 0.286 ± 0.019 and rE = 0.336 ± 0.031) and all statistically different from zero (P < 0.001). Similarly, high and significant estimates were obtained in bivariate analysis of LEAN-tot and FAT-tot (rG = 0.479 ± 0.017 and rE = 0.501 ± 0.022). Our data showed no detectable common genetic factors between BMD-tot and FAT-tot (rG = 0.000), however, shared environmental effects for these two variables were significant (rE = 0.209 ± 0.032).
Linkage analysis
First, LEAN-tot adjusted for only age (adjustment No. 1) was subjected to the multipoint linkage analyses with all available markers. The analysis revealed no significant linkage signals, and there was only one suggestive peak on X-chromosome at 160cM: LOD = 1.94. However, when LEAN-tot was adjusted for all significant covariates (adjustment No. 3), significant and suggestive linkage peaks were observed on chromosomal regions 1p36.13, 2q31.1, 7p15.3–15.1, 8p22, 12q24.3, and 14q22.3 (Table 3
). Table 3
shows only empirical LOD scores higher than 1.9. The maximum LOD score of 4.49 at the marker D12S1675 was observed on chromosome 12q24.3. Significant LOD = 3.62 at marker D14S276 was also observed on chromosome 14q22.3. Using permutation technique, we computed empirical P values for all six aforementioned LOD scores (Table 3
). Of these, four linkage peaks appeared reliable: chromosomes 7, 8, 12, and 14, with P = 0.0004 for chromosome 12 being the most impressive result. Finally, as an additional test of the reliability of the LOD scores observed for these four chromosomes, we conducted the bootstrap analysis of the data. It should be mentioned that 95% confidence intervals in all instances included nonsignificant and not indicative LOD scores (<1.9). However, as seen in Table 4
for chromosome 12q and 14q, the mean and 50% LOD scores were significant (>3.0), and marginally significant (>2.8) for the remaining two chromosomes. In case of the 12q, the proportion of the bootstrapped LOD scores greater than three was observed in 72%.
|
|
| Discussion |
|---|
|
|
|---|
The major goal of this study was to evaluate the contribution of the putative genetic factors to interindividual differences in LEAN-tot and to identify specific chromosomal regions that can harbor corresponding genes. Present results suggest that lean mass in females is strongly determined by genes. About 65.2% (±4.6%) of the variation of LEAN-tot, adjusted for significant covariates, was attributable to genetic factors effect in this study. The impact of heritability on body lean mass has been evaluated in several previous studies (14, 15, 16). However, heritabilities are population specific, and some caution is needed when comparing the results obtained in ethnically different populations. Nevertheless, comparison of the heritability estimates for LEAN-tot obtained in various studies showed a relatively narrow range of the interpopulation variation, ranging from 0.56 ± 0.05 (9) to 0.60 ± 0.11 (14), perhaps because of the weak independent age effect and common environmental influence.
To our knowledge, our genome scan of 2162 DZ twins represents the first large-scale genome scan for lean mass and suggests reliable linkage of this phenotype to QTL on two chromosomal segments. The strongest signal was observed on chromosomes 12q24.32 (associated with marker D12S1675) and 14q22.3 (D14S276). The multipoint LOD scores were 4.49 and 3.62, respectively, with accompanied empirical P values = 0.0004 and 0.0013. These LOD scores correspond to nominal P = 0.000005 and 0.000042 (Table 3
), and therefore achieved the P value 0.000049, significant for a genome-wide significance level of 0.05, according to guidelines proposed by Kruglyak and Lander (23), especially for the 12q segment. Both peaks were reasonably well supported by the bootstrapping of the data, which is important in light of the low reproducibility of many linkage results. They show that even for such a high LOD score as 4.49, observed in a large sample (n = 2162), a substantial proportion of the bootstrapped LOD scores (
28%) may be less than 3.0, and 12%, maybe even smaller than 1.9.
The additional marginally significant linkage peaks (LOD = 2.83 and 2.86) were observed on chromosome 7 (40cM from p terminus, D7S629 and D7S516, with an empirical P value = 0.0042) and on chromosome 8 (at 35cM, flanked by marker D8S258, with an empirical P value = 0.0044).
The best evidence for linkage in this study was detected at 12q24. Although, there are no prior linkage studies on total body lean mass, this region may certainly be of interest. First, recently our team found suggestive linkage for central fat mass to this chromosomal segment (29). In that study additional analysis of a number of available single nucleotide polymorphisms provided evidence for an association between central fat mass and two genes located on chromosome 12q24, PLA2G1B, and P2RX4 (12q23-q24.1; 119.24-.25Mb and 12q24.32; 120.13-.16Mb, respectively), with P values of 0.0067 and 0.017, respectively. However, of further interest in regard to LEAN-tot may be the fact that this region also harbors the acetyl-coenzyme A carboxylase-ß (ACACB alias ACCB, 12q24.11; 108.06-.19Mb) gene. The product of this gene is primarily expressed in heart and skeletal muscles. In these tissues it is thought to be deeply involved in control of fatty acid oxidation and the regulation of energy homeostasis. Moreover, Schrauwen et al. (30) found that even a minimal amount of physical training tends to increase fat oxidation and leads to marked changes in the expression of ACACB. In their study ACACB mRNA expression was significantly decreased after training (P = 0.005), whereas lipoprotein lipase mRNA expression tended to increase (P = 0.07). It has been shown in animal models that mutations in this gene may strongly be associated with adipose tissue metabolism and can accumulate a lesser amount of fat than wild-type mice on the same diet. The knockout mice at this gene do not survive even embryonic development (31). Relatively recently, Oh et al. (32) showed that the enzymatic product of this gene in humans is differently expressed in skeletal muscles and liver. They found that there are two promoters, P-I and P-II, controlling the transcription of the ACACB gene, of which P-I is active in the skeletal muscles and heart, but not in the liver. It would certainly be of importance for a future study to test whether the linkage we observed is associated with P-I.
As mentioned, currently there are very scarce published data on the genetics of total body lean mass. Two similar linkage studies to ours using a slightly different phenotype are the Quebec Family Study (33), which used 748 subjects genotyped for 292 microsatellite markers to find genes linked to fat-free mass (FFM), and the HERITAGE Study (34), which performed genome scan (with the use of 344 markers) for FFM at baseline and after 20 wk of endurance exercise in 364 sib pairs. However, their method implements underwater weighing, and it does not differentiate between bone and muscle mass. Therefore, their results are not directly comparable. Nevertheless, they reported among their significant linkage results chromosome 7p15.3 (P = 0.0002), observed also in the present study. The authors suggested that this linkage is probably related more to the skeletal muscle component of FFM than to BMD.
Because linkage peaks have wide confidence intervals (35, 36), which harbor hundreds of genes, the observed linkage signals are difficult to translate into obvious candidate genes with any confidence. However, within chromosome 7p15.1–15.3 mapped at least two genes of interest: the GSBS (7p15; 31.69-.71Mb) gene causing susceptibility to hypercholesterolemia (37), and DSMAV (7p15; 30.60-.64Mb) gene, mutations in which were associated with spinal muscular atrophy type V (38). However, there are no available data suggesting the association of these genes with lean mass, and ACACB gene looks, therefore, to be a much more promising candidate for follow-up.
There are some limitations to this study. The measurement method implemented in this project does not differentiate between the muscles and tendons from the one side and other soft tissues (viscera) on the other side. It is assumed that this method is inferior to computed tomography, although the recent comparison demonstrated good concordance between DXA and computed tomography for abdominal total tissue mass and fat mass (39). Moreover, in this study DXA also showed excellent reliability among three different operators to determine fat and lean body mass in the L1–L4 region of interest. Another limitation is that the obtained results may not be directly applicable to males having in general much more developed musculature and more massive bones. On the other side, using only females may not necessarily be a limitation but rather a possible advantage because it is not prone to potential genotype x sex interactions, which are likely the rule rather than exception (40). This is particularly important for body composition components, which are substantially different between the sexes due to fundamental differences in sex-related hormonal status.
There are additional major strengths of this study. These families were selected randomly from the relatively homogeneous ethnic background, and, therefore, are informative for genetic-mapping studies and representative of the general population. If there were substantial genetic heterogeneity, our study would have limited power to detect linkage signals. Age matching of twins and shared household decreased substantially random effects of environment, and made heritability estimates and linkage results more accurate than might be expected for the heterogeneous composition of the complex pedigree.
In summary, we have shown that LEAN-tot variation is a highly heritable trait. We have observed for the first time reliable linkage between LEAN-tot and DNA markers located on chromosome 12q and 14q, and good evidence of linkage on chromosome 7p and 8p. The linkage peaks on chromosomes 12q and 7p were confirmed in previous studies of related phenotypes. Of special interest is the linkage on 12q that may be caused by the physiologically important and relevant ACACB gene. Replication studies from other populations followed by fine mapping should be the next step in identifying the individual genes involved that could have important physiological effects on bone and fat, as well as muscle.
| Footnotes |
|---|
Disclosure Statement: The authors have nothing to declare.
First Published Online June 5, 2007
Abbreviations: ACACB, Acetyl-coenzyme A carboxylase-ß; BMD-tot, total bone mineral density; DXA, dual-energy x-ray absorptiometry; DZ, dizygotic; FAT-tot, total fat; FFM, fat-free mass; IBD, identical-by-descent; LEAN-tot, total lean body mass; MZ, monozygotic; QTL, quantitative trait loci; rE, environmental correlation; rG, genetic correlation.
Received February 23, 2007.
Accepted May 30, 2007.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
F. E. Dewey, D. Rosenthal, D. J. Murphy Jr, V. F. Froelicher, and E. A. Ashley Does Size Matter?: Clinical Applications of Scaling Cardiac Size and Function for Body Size Circulation, April 29, 2008; 117(17): 2279 - 2287. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 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 |