DOI: 10.2337/diabetes.55.02.06.db05-0997 © 2006 by the American Diabetes Association
Identification of TUB as a Novel Candidate Gene Influencing Body Weight in Humans
1 Department of Molecular Genetics, Maastricht University, Maastricht, the Netherlands Address correspondence and reprint requests to Marten H. Hofker, Moleculaire Genetica (UNS50/11), Universiteit Maastricht, Universiteitssingel 50, Postbus 616, Maastricht 6200 MD, Netherlands. E-mail: m.hofker{at}gen.unimaas.nl
Key Words: SNP, single nucleotide polymorphism
Previously, we identified a locus on 11p influencing obesity in families with type 2 diabetes. Based on mouse studies, we selected TUB as a functional candidate gene and performed association studies to determine whether this controls obesity. We analyzed the genotypes of 13 single nucleotide polymorphisms (SNPs) around TUB in 492 unrelated type 2 diabetic patients with known BMI values. One SNP (rs1528133) was found to have a significant effect on BMI (1.54 kg/m2, P = 0.006). This association was confirmed in a population enriched for type 2 diabetes, using 750 individuals who were not selected for type 2 diabetes. Two SNPs in linkage disequilibrium with rs1528133 and mapping to the 3' end of TUB, rs2272382, and rs2272383 also affected BMI by 1.3 kg/m2 (P = 0.016 and P = 0.010, respectively). Combined analysis confirmed this association (P = 0.005 and P = 0.002, respectively). Moreover, comparing 349 obese subjects (BMI >30 kg/m2) from the combined cohort with 289 normal subjects (BMI <25 kg/m2) revealed that the protective alleles have a lower frequency in obese subjects (odds ratio 1.32 [95% CI 1.04–1.67], P = 0.022). Altogether, data from the tubby mouse as well as these data suggest that TUB could be an important factor in controlling the central regulation of body weight in humans. Until recently, type 2 diabetes was considered a disease of the elderly. Currently, worldwide, a large proportion of newly diagnosed patients are adolescents. The increase in the prevalence of type 2 diabetes, as well as the rapid spreading to a younger age at onset, is largely due to environmental factors, including modern eating habits and reduced physical activity. Obesity is a major risk factor for type 2 diabetes, and the obesity epidemic coincides with the type 2 diabetes epidemic. Genetic studies have suggested a common basis for both diseases. Recently, we identified a locus on 11p15 (95% CI 4,868,745–10,676,565 bps) influencing obesity in families with type 2 diabetes (1). The TUB gene was found to be the most relevant candidate gene within this locus, and we hypothesized that TUB was associated with obesity susceptibility.
TUB is the founding-member of the tubby-like proteins and is conserved among vertebrate genomes (2). Interestingly, a loss of function mutation of tub results in the tubby mouse syndrome, which is characterized by late-onset obesity and neurosensory deficits. The tubby mouse begins to diverge in weight at To investigate the role of TUB in human obesity, 13 single nucleotide polymorphisms (SNPs) around TUB (Fig. 1A and online appendix Table 1 [available from http://diabetes.diabetesjournals.org]) were genotyped in 492 unrelated type 2 diabetic patients from the Breda cohort, for whom BMI values were available (Table 1). Linear regression analysis, adjusted for age and sex, revealed that the minor allele of SNP rs1528133 had a marked effect on BMI (+1.54 kg/m2, P = 0.006) (Fig. 1B, Table 2). SNP rs1528133 is located 22 kb distal to TUB in the flanking gene RIC3, which has an unknown function. SNP rs1528133 is in strong linkage disequilibrium with SNPs in the 3' end of TUB (Fig. 1C), indicating that the COOH-terminus of TUB may be associated with obesity. A replication study of rs1528133 and four SNPs located in the 3' of TUB and in linkage disequilibrium with rs1528133 (rs2272382, rs2272383, rs3750955, and rs1406095) was performed in a population enriched for type 2 diabetes Dutch cohort (RIVM) of 750 individuals (Table 1). We were able to validate our initial result and found a significant effect on BMI for the major alleles of both rs2272382 (–1.28 kg/m2, P = 0.016) and rs2272383 (–1.26 kg/m2, P = 0.01) (Fig. 1B, Table 2). These two SNPs are in a strong linkage disequilibrium with rs1528133 (D' = 0.8 and D' = 0.9, respectively). Combining both cohorts (n = 1,242) strengthened the effect, as rs2272382 and rs2272383 showed even stronger association with BMI (–1.12 kg/m2, P = 0.005, –1.11 kg/m2, P = 0.002) (Fig. 1B, Table 2). These two SNPs remained significant in the combined cohort after stringent Bonferroni correction for five SNPs (Pc < 0.05/5 = 0.01).
Comparing the 349 obese subjects from the combined cohort (BMI >30 kg/m2) with 289 normal subjects (BMI <25 kg/m2) revealed that the minor allele of SNP rs2272382 was significantly more common in the obese subjects (37.3%) than normal subjects (31.1%) (odds ratio 1.32 [95% CI 1.04–1.67], P = 0.022), while the minor allele of rs2272383 showed borderline significance (1.23 [0.98–1.55], P = 0.078) (Table 3). In addition, the genotype homozygous for the minor allele of SNP rs2272382 was also significantly more common in the obese subjects (14.2%) compared with normal subjects (7.1%) (P = 0.016) (online appendix Table 2). Notably, the RIVM cohort includes 55% type 2 diabetic subjects. In a case-control study comparing patients with type 2 diabetes with healthy subjects from the Breda cohort for all 13 SNPs and RIVM cohort for all five SNPs, we found no associations (data not shown), clearly confirming the genes direct influence on body weight independent of type 2 diabetes. Therefore, the data suggests that TUB may influence BMI in general and may not be involved in type 2 diabetes directly.
The two associated SNPs, rs2272382 and rs2272383, are located in the highly conserved COOH-terminal domain of tubby, which contains the DNA-binding domain and the phosphatidylinositol-binding region. Strikingly, mutations in the murine TUB allele and tubby-like proteins were found to be in the same region (6,7). Comparing the linkage disequilibrium plot (D') between these three SNPs in the Breda case subjects and the case and control subjects from RIVM revealed no major differences (data not shown). Moreover, the allele frequencies were similar. It is therefore possible that the differences between the Breda and RIVM cohorts can be attributed to different patterns of association between the marker allele and the true, not yet identified, disease-associated allele. Alternatively, it is possible that different alleles at the same locus are responsible for increased risk in different populations. It is tempting to speculate that a variation in this region may have affected a domain important to the biological function of tubby. However, it should be noted that the region showing association also encompasses the RIC3 gene, which is known to be involved in neurotransmission (8). Although we cannot formally exclude that RIC3 is contributing to body weight, the TUB gene is a more attractive candidate gene for body weight due to its function. To our knowledge, the potential contribution of TUB to human obesity was assessed only once (9). The authors have examined five polymorphic microsatellite markers around TUB in 716 Pima Indians comprising 217 nuclear families for sibpair linkage to BMI. No significant linkages were found in an analysis of all sibships or in an analysis restricted to discordant sibpairs. However, this linkage study was done in a different population and not with the same SNPs that were used by us. Thus, the role of tubby as an important factor in metabolism and obesity has not been shown in humans. Our study is the first to provide an initial estimate of the association of the TUB gene with a quantitative measure of obesity. Recently, the effect of MC4R, a major gene for obesity, on weight regulation was estimated to be –0.52 kg/m2 (10), and since the effect of TUB on weight regulation is more then twice as high (–1.3 kg/m2), this indicates the importance of TUBs contribution to polygenetically regulated body weight. Altogether, data from the tubby mouse as well as these data suggest that tubby has a protective role in the central regulation of body weight in humans. Mutations in TUB can lead to increased body weight and contribute to obesity. This discovery will reveal new molecular pathways in weight regulation that should lead to innovative therapies, preventive measures, and insights into the pharmacogenetics of such intervention strategies.
Subjects DNA, isolated from whole blood, was available from two independent cohorts of Dutch subjects: the Breda and the RIVM cohorts (Table 1). The Breda cohort comprised 492 individuals (267 women and 224 men). The level of obesity in each individual was given by the BMI, defined as weight in kilograms divided by the square of height in meters. All the patients were diagnosed according to World Health Organization criteria. There was also a cohort of control subjects available (n = 356) comprised of general random subjects. The RIVM cohort consisted of participants in a health monitoring project from two Dutch towns (Doetinchem and Maastricht). Subjects were diagnosed with diabetes according to venous plasma glucose concentrations ( 11.1 mmol/l for nonfasting and 7.0 mmol/l for fasting or were being treated with tablets or insulin at the time of the health monitoring project [1993–1997] or of a follow-up questionnaire [1998]). Control subjects were normoglycemic (venous plasma glucose <5.5 mmol/l for nonfasting and <6.1 mmol/l for fasting participants). We excluded subjects if either of their parents was not born in a European country. We calculated the BMI as weight (minus 1 kg for measured weight to correct for clothing) divided by the square of height in meters (kg/m2). The Breda study was approved by the University Medical Center Utrecht. The RIVM study was approved by the TNO Leiden medical ethical committee. Written informed consent was obtained from all subjects. Experiments were conducted according to the principles expressed in the Declaration of Helsinki.
SNP selection and genotyping.
Genotype analysis.
Statistical analysis.
We thank the Dutch Diabetes Foundation (DFN 97.144), the Dutch Ministry of Welfare and Sport, the European Vascular Genomics Network, and the Netherlands Heart Foundation for financial support. A.C. is supported by a Kootstra Fellowship, Maastricht University. The authors thank Jackie Senior for editing the manuscript.
J.H.O.V.T. is currently affiliated with the Department of Human Biology, Maastricht University, Maastricht, the Netherlands. Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org. Received for publication August 4, 2005 and accepted in revised form November 8, 2005
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||