|
|
||||||||
1 Osteoporosis Research Center, Creighton University Medical Center, Omaha, Nebraska
2 Department of Biomedical Sciences, Creighton University Medical Center, Omaha, Nebraska
3 Laboratory of Molecular and Statistical Genetics, College of Life Sciences, Hunan Normal University, ChangSha, Hunan, China
| ABSTRACT |
|---|
|
|
|---|
D'
0.940, P < 0.0001) were observed between pairs of SNPs. In single-locus analyses, we found statistically significant association (P = 0.034) and linkage (P = 0.031) between 55 C/T and BMI. This polymorphism explains 2.29% of BMI variation, and subjects carrying the T allele had an average of 3.5% lower BMI than those without it (P = 0.003). In haplotype analyses, we also observed evidence of linkage (P = 0.002) and association (P = 0.035) with BMI. In summary, our results suggest that UCP3 gene polymorphisms may contribute to BMI variation in this Caucasian population. uncoupling protein 3; quantitative transmission disequalibrium test
| INTRODUCTION |
|---|
|
|
|---|
Uncoupling proteins (UCPs) are a family of inner mitochondrial membrane transporters that dissipate the proton gradient, releasing stored energy as heat, and are considered potentially important determinant of defense against obesity (9). Three distinct UCPs (UCP1, UCP2, and UCP3) have been identified. UCP1, existing only in brown adipose tissue, is unlikely to be a major gene involved in body weight regulation in humans (24). UCP2 has a wide tissue distribution at varying levels (15). UCP3 is expressed predominantly in skeletal muscle, a major site of thermogenesis in humans, making UCP3 an attractive target for studies toward body weight regulation (42). In rodents, UCP3 expression was regulated by thyroid hormone, ß3-adrenergic agonists, leptin, and fat feeding (17, 28, 33). Transgenic mice overexpressing UCP3 in muscle were hyperphagic but lean and resistant to diet-induced obesity (8). In Pima Indians, UCP3 mRNA expression was positively correlated with the sleeping metabolic rate (39). UCP3 has also been shown to be involved in the handling of fatty acids to maintain mitochondrial oxidative capacity (5, 22, 47).
The human UCP3 gene contains seven exons and spans
8.7 kb (42). UCP3 was mapped within 7 kb of UCP2 on chromosome 11q13 (42). In humans, linkage was reported on chromosome 11q13 in the vicinity of the UCP2-UCP3 genes for resting metabolic rate, percent body fat, and fat mass (6). In mice, a quantitative trait locus (QTL) for obesity was mapped to a region that contains the UCP2-UCP3 gene locus on chromosome 7 (44). A number of polymorphisms in the UCP3 gene have been identified, and the correlation with obesity-related phenotypes was evaluated (9). However, the results have been inconsistent.
The purpose of this study was to investigate whether the UCP3 gene contributes to obesity risk in the general population. Using the tests implemented in the statistical package quantitative transmission disequilibrium test (QTDT), we examined association and linkage between UCP3 gene polymorphisms and obesity-related phenotypes in a large sample of 405 Caucasian nuclear families comprising 1,873 subjects.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
1%.
Measurements
BMI was calculated as the ratio of total body weight to the square of height (in units of kg/m2). Body weight was measured in light indoor clothing, using a calibrated balance beam scale, and height was measured using a calibrated stadiometer. Fat mass and lean mass were measured by dual-energy X-ray absorptiometry using a Hologic 2000+ or 4500 scanner (Hologic; Bedford, MA). Both machines were calibrated daily. The body composition bar was used on every whole body scan on the Hologic 2000+. On the Hologic 4500, the bar was not needed for the body scans; instead, it was scanned every week. Percent fat mass (PFM) is the ratio of fat mass to the sum of fat mass, lean mass, and bone mass.
The measurement precision of BMI as reflected by the coefficient of variation was 0.2%. The coefficients of variation for fat mass, PFM, and lean mass were 2.2%, 2.2%, and 1.0%, respectively, for measurements obtained on the Hologic 2000+ and were 1.2%, 1.1%, and 0.7%, respectively, for measurements on the Hologic 4500. Members of the same nuclear family were usually measured on the same type of machine. Except for lean mass and fat mass, the phenotypes are correlated significantly. The phenotypic correlations were, respectively, 0.88 (between BMI and fat mass), 0.49 (between BMI and lean mass), 0.54 (between BMI and PFM), 0.12 (between fat mass and lean mass), 0.81 (between fat mass and PFM), and 0.34 (between lean mass and PFM).
Statistical Analyses
Single-locus analyses.
The allele frequencies of the tested SNPs were estimated in all the subjects using a maximum-likelihood method implemented in the program SOLAR (4). Using a variance component analysis for quantitative traits (which is implemented in SOLAR), we estimated heritability of the four obesity phenotypes. Age and sex were adjusted as covariates in the analyses.
For quantitative traits, the essence of the transmission disequilibrium test is to examine the difference between average phenotypic values of children with different alleles transmitted from a heterozygous parent (3). Under a flexible variance-component framework, tests of population stratification, linkage, total association, and within-family association between each of the SNPs and obesity phenotypes are implemented in the statistical software package QTDT (1). The permutation procedure built in the QTDT may yield significance levels (P values) of the tests that are not biased. The orthogonal model of Abecasis et al. (1) was adopted in the analyses, where the genotype score is decomposed into orthogonal between-family (ßb) and within-family (ßw) components. Population stratification is examined by testing whether ßb = ßw. Linkage tests are based on the standard variance-component methods and the identity by descent among relatives. Total association tests use all information including ßb and ßw components and may yield false positive/negative results due to population stratification. Within-family association (via the transmission disequilibrium test), however, is significant only if there is linkage disequilibrium (LD), and it is robust to population stratification. When significant association is observed, the approximate phenotypic variation due to the detected marker is calculated as 2p(1 p)a2/VP, where VP is the total phenotypic variance (VP = Ve + Vg + Va), p is the allele frequency of the marker, and a is the estimate of additive effect, i.e., E(ßw) = a (1). Ve is the residual component of variance, Vg is the polygenic component of variance, and Va is the additive component of variance.
A trade off for conducting the within-family association test is the reduction in power to detect allelic association (45). In the absence of population stratification, total association confers more power and is more sensitive than within-family association to detect correlation between a marker and an underlying trait (45). Thus, in light of such power consideration, we first tested population stratification. If no evidence of population stratification was observed, the more powerful test of total association would be utilized.
Variance-component methods implemented in the QTDT make a critical assumption that the underlying trait follows a multivariate normal distribution. In the present study, all four obesity phenotypes showed a marked departure from normal distribution (P < 0.01) using the Anderson-Darling test and, therefore, were transformed to approach normality using the Box-Cox procedure implemented in the statistical software MINITAB (Minitab; State College, PA). Before the Box-Cox transformation, the values of kurtosis for BMI, fat mass, PFM, and lean mass were 1.59, 1.38, 0.43, and 0.82, respectively, and the values of skewness were 0.96, 0.98, 0.04, and 0.55, respectively. After Box-Cox transformation, the values of kurtosis for BMI, fat mass, PFM, and lean mass were 0.30, 0.01, 0.42, and 0.43, respectively, and the values of skewness were 0.02, 0.05, 0.02, and 0.16, respectively. All the statistical analyses were performed on the transformed data with age and sex as covariates, both significantly affecting the obesity phenotypes.
Haplotype analyses.
Pair-wise LD between the SNPs was calculated by the normalized measure D' (29). The statistical significance of the observed LD was examined by Monte-Carlo approximation of Fisher's exact test (38). Haplotypes were constructed using the program Genehunter version 2.1 (26). Genehunter extracts complete multipoint inheritance information to infer maximum likelihood haplotypes for all individuals in nuclear families. Haplotype frequencies were estimated in unrelated subjects (parents from each nuclear family). For QTDT haplotype analyses, we first utilized the QTDT multiallelic tests option, which include alleles in one test as a categorical variable with allele frequencies of <5% being pooled to give a global significance. We then performed the allele-wise test on each allele separately, with all other alleles being pooled into another category. Employing the aforementioned sequential tests within QTDT,
2 and P values were computed for those haplotypes present in at least 30 informative offspring.
| RESULTS |
|---|
|
|
|---|
Because D' can be inflated in samples of small sizes, especially for SNPs with rare alleles (48), to avoid potential bias, we characterized the pair-wise LD and constructed haplotypes for the three highly polymorphic SNPs (i.e., 55 C/T, Tyr99Tyr, and Tyr210Tyr), which were genotyped in the total sample of 1,873 subjects. The three SNPs were arranged in the order of 55 C/T Tyr99Tyr Tyr210Tyr in haplotype construction. All eight possible haplotypes were observed. The three most common haplotypes were CTC, CTT, and TCT, with allele frequencies of 40.6%, 31.2%, and 23.5%, respectively. These three haplotypes accounted for >95% of the total allele frequencies. The allele frequencies of the selected five SNPs and haplotypes defined by 55 C/T Tyr99Tyr Tyr210Tyr are shown in Table 2.
|
D'
0.940, P < 0.0001) were found between pairs of SNPs. We did not find significant deviation from the Hardy-Weinberg equilibrium for any SNP/phenotype.
Linkage and Association Between Individual SNPs and Obesity Phenotypes
Table 3 summarizes the results of population stratification, association, and linkage tests. In our sample, the heritability of BMI, fat mass, PFM, and lean mass was 42%, 33%, 58%, and 51%, respectively. We did not find significant evidence of population stratification for any of the phenotypes. Because no population stratification was observed, we utilized the more powerful test of total association instead of within-family association. In single locus analyses, evidence of association was observed between 55 C/T and BMI (P = 0.034). This polymorphism contributes to 2.29% of BMI variation. In the children's group, BMI values of the subjects with CC, CT, and TT genotypes were 27.32 ± 5.33, 26.74 ± 5.21, and 25.96 ± 5.18, respectively. Subjects carrying the T allele had an average of 3.5% lower BMI than those without it (P = 0.003). We also found evidence of linkage between 55 C/T and BMI (P = 0.031). Evidence of linkage was also observed for both Tyr99Tyr (P = 0.032) and Tyr210Tyr (P = 0.027) with BMI, suggesting that the UCP3 gene locus is linked to BMI in this study population. However, we did not find a significant association between these two SNPs and BMI. For three other phenotypes, fat mass, PFM, and lean mass, we did not find evidence of linkage or association for any of the SNPs.
|
| DISCUSSION |
|---|
|
|
|---|
The UCP3 55C/T transition is placed 55 bp upstream of the most commonly used transcription initiation site of skeletal muscle (14). This polymorphism is potentially interesting because it is located only 6 bp from the TATA box and 4 bp downstream of a putative peroxisome proliferator-activated receptor (PPAR) responsive element and could modify the PPAR responsiveness of the UCP3 gene (2, 41). Thus the UCP3 gene could be one of the PPAR-
targets involved in the modulation of lipid metabolism and insulin sensitivity. In French Caucasians, the subjects bearing the TT genotype of the UCP3 55C/T polymorphism had a lower risk for developing Type 2 diabetes than others did (34). In Pima Indians, subjects who carried the 55 T allele had significantly higher UCP3 mRNA concentration than 55 C homozygotes (40), suggesting that the 55 T allele could increase UCP3 mRNA expression compared with the C allele. In addition, the mRNA expression of UCP3 was positively correlated with the sleeping metabolic rate in Pima Indians, and a decrease in the full-length UCP3 message was associated with a increase in BMI, which is conceived to be due to an observed reduction in fat oxidation rates (39). Therefore, a possible explanation for our results is that the 55 T allele is associated with an increase in the UCP3 message, which, if similar to the situation in Pimas, is associated with increased fat oxidation and reduced BMI. Further studies are needed to assess the functional consequence of the 55 C/T polymorphism by examining its direct effects on expression of UCP3.
We also tested the UCP3 Tyr99Tyr and Tyr210Tyr polymorphisms for linkage and association with obesity phenotypes. These two SNPs were selected because they were highly polymorphic in our sample and may render relatively high statistical power, as the statistical power of a family-based association study depends directly on the degree of allelic heterozygosity. Although we did not find any significant association for Tyr99Tyr and Tyr210Tyr with obesity phenotypes, there was some evidence of linkage, suggesting that they might be close to a functional variant. Haplotype-based analysis, which considers the variants segregating at multiple loci, can provide additional information in studies of complex traits. This has been reflected in our linkage tests for haplotypes, in which the evidence of significance has increased markedly. Haplotypes can also provide more information for the transmission disequilibrium test analyses. However, one trade off is the increase in the degrees of freedom of the test. This may partially explain why in allele-wise haplotype tests the significance of association between haplotype CTC and BMI (P = 0.063) has decreased a bit compared with that between 55 C/T and BMI (P = 0.034) in single locus analyses.
Recently, there are some debates on whether the primary physiological function of UCP3 is the regulation of energy expenditure (21). For instance, fasting, an energy expenditure attenuating condition, unregulated UCP3 expression (35), and UCP3 knockout mice had a normal metabolic rate and normal body weight, even though their mitochondria showed improved coupling (18, 46). Although the exact function of UCP3 remains to be elucidated, there has been accumulating evidence showing that UCP3 is involved in the regulation of the production of reactive oxygen species, mitochondrial fatty acid transport, and the regulation of glucose metabolism in skeletal muscle (21). Thus an important function of UCP3 might be influencing, if not necessarily regulating, energy metabolism and body weight (21). As such, UCP3 remains an interesting target for the intervention of obesity.
Earlier studies on associations between the UCP3 gene and obesity have yielded inconsistent results. For example, association was reported in a French cohort of morbidly obese subjects, in which the 55 T/T genotype was associated with increased BMI (37). However, in Danish Caucasian subjects, no association was observed between 55 C/T and BMI or long-term body weight change (10). In African American women, no significant variation in resting energy expenditure (REE) was seen for 55 C/T; but for Tyr210Tyr, REE was significantly lower in subjects with the CC genotype than in those with the TT genotype (23).
There are several possible reasons for this inconsistency. First, there might be true variability in association among different populations, particularly of different ethnic groups (31). Even when a causal variant is under investigation, it might be more or less important in different populations, especially if the variant has low genetic effects, variable penetrance, and variable allele frequencies in different populations. In terms of allele frequencies, differences exist for the UCP3 55 C/T polymorphism between our United States Caucasian nuclear families and French cohort (37) and African American women (23), with the minor allele frequencies being 26.5%, 22%, and 21%, respectively. For Tyr210Tyr, the minor allele frequency was 44.7% in our United States Caucasian families but was 22.7% in African American women (23). Second, diversity in study designs may also contribute to the situation. For example, population stratification, if not checked and controlled, could lead to spurious association outcomes (12). Recent empirical studies demonstrated that modest amount of stratification can exist even in well-designed studies (16, 32). In this regard, our results should be robust because we adopted the transmission disequalibrium test (a family-based disequilibrium test), which obviates concerns about population stratification by examining the transmission of alleles from parents to offspring.
False positive associations may exist in the literature, but there are also many real associations lurking among the data (31). These true associations probably confer modest but real effects on common disease risk. Studies using large samples may convincingly identify such variants, whereas studies of small samples may miss the chance of detection. Our large sample containing 1,873 subjects from 405 nuclear families confers a relatively high statistical power. Assuming a marker is in strong LD with a functional mutation (D' = 0.9) that explains 2% of phenotypic variation, this sample has more than 80% power to detect the association via the transmission disequalibrium test.
In this study, the effect of UCP3 55 C/T on BMI was not reflected on variation in three other phenotypes (fat mass, PFM, and lean mass), although BMI is phenotypically correlated with them. The lack of association may be partially due to insufficient power to detect changes in these phenotypes. On the other hand, a significant and high phenotypic correlation does not necessarily imply a significant and high genetic correlation (13). A limitation with this study is that the quantitative traits studied here were focused on BMI, fat mass, PFM, and lean mass. Some other phenotypes may also be valuable in studying obesity risk. For example, the waist-to-hip ratio is a better index than BMI in reflecting central obesity, which is associated with increased cardiovascular risk. Earlier studies (e.g., Refs. 7 and 20) have found that the UCP3 55 C/T polymorphism was associated with the waist-to-hip ratio in Caucasian and Asian populations.
In summary, our results suggest that the UCP3 55 C/T polymorphism may contribute to BMI variation in the general population. Confirmation of the association observed herein awaits further studies such as meta-analysis and/or independent replication studies. Further in vitro functional analyses and clinical studies are also necessary to define the precise contribution of this variant to human obesity.
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: H.-W. Deng, Osteoporosis Research Center, Creighton Univ., 601 N. 30th St., Suite 6787, Omaha, NE 68131 (E-mail: deng{at}creighton.edu).
10.1152/physiolgenomics.00031.2005.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. F M van Abeelen, M. de Krom, J. Hendriks, D. E Grobbee, R. A H Adan, and Y. T van der Schouw Variations in the uncoupling protein-3 gene are associated with specific obesity phenotypes Eur. J. Endocrinol., May 1, 2008; 158(5): 669 - 676. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Bezaire, E. L. Seifert, and M.-E. Harper Uncoupling protein-3: clues in an ongoing mitochondrial mystery FASEB J, February 1, 2007; 21(2): 312 - 324. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Liang and B. Ventura Physiological genomics in PG and beyond: July to September 2005 Physiol Genomics, October 17, 2005; 23(2): 119 - 124. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |