Skip to main content
  • Research article
  • Open access
  • Published:

Association between TCF7L2gene polymorphisms and susceptibility to Type 2 Diabetes Mellitus: a large Human Genome Epidemiology (HuGE) review and meta-analysis

Abstract

Background

Transcription factor 7-like 2 (TCF7L2) has been shown to be associated with type 2 diabetes mellitus (T2MD) in multiple ethnic groups in the past two years, but, contradictory results were reported for Chinese and Pima Indian populations. The authors then performed a large meta-analysis of 36 studies examining the association of type 2 diabetes mellitus (T2DM) with polymorphisms in the TCF7L2 gene in various ethnicities, containing rs7903146 C-to-T (IVS3C>T), rs7901695 T-to-C (IVS3T>C), a rs12255372 G-to-T (IVS4G>T), and rs11196205 G-to-C (IVS4G>C) polymorphisms and to evaluate the size of gene effect and the possible genetic mode of action.

Methods

Literature-based searching was conducted to collect data and three methods, that is, fixed-effects, random-effects and Bayesian multivariate mete-analysis, were performed to pool the odds ratio (OR). Publication bias and study-between heterogeneity were also examined.

Results

The studies included 35,843 cases of T2DM and 39,123 controls, using mainly primary data. For T2DM and IVS3C>T polymorphism, the Bayesian OR for TT homozygotes and TC heterozygotes versus CC homozygote was 1.968 (95% credible interval (CrI): 1.790, 2.157), 1.406 (95% CrI: 1.341, 1.476), respectively, and the population attributable risk (PAR) for the TT/TC genotypes of this variant is 16.9% for overall. For T2DM and IVS4G>T polymorphism, TT homozygotes and TG heterozygotes versus GG homozygote was 1.885 (95%CrI: 1.698, 2.088), 1.360 (95% CrI: 1.291, 1.433), respectively. Four ORs among these two polymorphisms all yielded significant between-study heterogeneity (P < 0.05) and the main source of heterogeneity was ethnic differences. Data also showed significant associations between T2DM and the other two polymorphisms, but with low heterogeneity (P > 0.10). Pooled ORs fit a codominant, multiplicative genetic model for all the four polymorphisms of TCF7L2 gene, and this model was also confirmed in different ethnic populations when stratification of IVS3C>T and IVS4G>T polymorphisms except for Africans, where a dominant, additive genetic mode is suggested for IVS3C>T polymorphism.

Conclusion

This meta-analysis demonstrates that four variants of TCF7L2 gene are all associated with T2DM, and indicates a multiplicative genetic model for all the four polymorphisms, as well as suggests the TCF7L2 gene involved in near 1/5 of all T2MD. Potential gene-gene and gene-environmental interactions by which common variants in the TCF7L2 gene influence the risk of T2MD need further exploration.

Peer Review reports

Background

Type 2 diabetes mellitus (T2DM) is characterized by hyperglycemia that can occur through mechanisms such as impaired insulin secretion, insulin resistance in peripheral tissues and increased glucose output by liver [1]. Most individuals with T2DM suffer serious complications of chronic hyperglycemia, involved in nephropathy, neuropathy, retinopathy, and accelerated development of cardiovascular disease [2]. The prevalence of T2DM has increased sharply during the past two decades in the world and is close to 6% [3]. But estimations suggest that this trend will continue to rise over the next decade owing to increasing age of population and surge of obesity [4]. Accounting for 90~95% of populations with diabetes, T2DM is projected to 300 million worldwide by year 2025 [2, 4], and its incidence exhibits significant geographic and racial difference, ranging from more than 0.3 to 17.9% in Africa [5–11], 1.2 to 14.6% in Asia [12–19], about 2.5% in Australia [20], 0.7 to 11.6% in Europe [21–32], 4.6 to 40% in the Middle East [33–36], 6.69 to 28.2% in North America [3, 37–39], and 2.01 to 17.4% in South America [40–43]. Recent years, new diagnosis of T2DM subjects in developing nations with a mass of populations such as China and India has increased rapidly [19, 44, 45], and it has also found remarkable augmented in rural that being with lower prevalence before [14, 17, 19, 45]. Incidence rates are diverse in females and males in Europe, North America, Africa, Latin America and East Asia, and there are scarcely reports with accordant result.

Transcription factor 7-like 2 (TCF7L2) gene spans a 215,863 bases region on chromosome 10q25.3 [46] (114700201–114916063, NCBI build 36.2), and its product is a high-mobility box-containing transcription factor that has a role in activating many genes downstream of the Wnt signaling pathway and in T2DM [47–50]. The mechanisms through which TCF7L2 affects the susceptibility to T2DM remain to be elucidated. The bipartite transcription factor, cat/TCF7L2, activates many genes downstream of the Wnt signaling cascade [51]. One of the genes transcriptionally regulated by cat/TCF7L2 is proglucagon, which encodes the insulinotropic hormone glucagon-like peptide 1 (GLP-1) [49, 50]. Because GLP1, in concert with insulin, plays a critical role in blood glucose homeostasis, it has been postulated that TCF7L2 gene variants may affect the susceptibility to type2 diabetes by indirectly altering GLP-1 levels [50].

The human TCF7L2 gene consists of 14 exons and 13 introns (NCBI build 36.2). However, a previous study [52] has shown that TCF7L2 has 17 exons, of which five are alternative. There are at least four well-studied polymorphic markers that have brought the most attention, that is, (1) a C-to-T (genomic position: 114748339) substitution at snp rs7903146 of the intron 3 (IVS3C>T) [2, 47, 53, 54]; (2) a T-to-C (genomic position: 114744078) substitution at snp rs7901695 of the intron 3(IVS3T>C) [2, 53, 55]; (3) a G-to-T (genomic position: 14798892) substitution at snp rs12255372 of the intron 4 (IVS4G>T) [2, 56]; and (4) a G-to-C (genomic position: 114797037)substitution at snp rs11196205 of the intron 4 (IVS4G>C) [57, 58]. Most published epidemiologic studies have emphasized on the IVS3C>T and IVS4G>T polymorphisms [2, 47, 59–65]. There are plenty of evidence on strong linkage disequilibrium between these four polymorphisms [2, 47, 57, 58]. The T allele of the IVS3C>T, C allele of the IVS3T>C, and T allele of IVS4G>T, as well as C allele of IVS4G>C polymorphisms have higher frequency in Caucasian [59, 66], North European [2, 47], African [67, 68], and Indian [69, 70] populations than those in East Asian [57, 61, 62, 71] populations.

The polygenic basis of T2DM has been intensely investigated by large consortia worldwide, and consensus on a range of valid susceptibility genes has emerged. Such genes, for example PPARG, KCNJ11 and CPN10, appear to contain common variation that confers only 10–20% additional risk of T2DM, and only studies with large numbers of case and controls have had sufficient statistical power to examine their role robustly [72]. Recent years, some novel loci in genes such as TCF7L2, SLC30A8, HHEX, ETX2, CDKN2A/CDKN2B, IGF2BP2, CDKAL1, FTO, PPARG, and KCNJ11 have been found with susceptibility to T2DM [2, 73–76], and these has promoted the development of genetics on T2DM. From 2006 to 2008, study on TCF7L2 gene polymorphisms and its association with T2DM has become the focus in diabetes mellitus, therefore, some meaningful results have reported one by one. In March, 2006, Grant and colleagues [2] first proposed an association with TCF7L2 gene variants with T2DM. The study [2] consisted of three samples from Iceland, Denmark and USA with strikingly similar results. Since then, this potential association has been examined in over 50 samples from around the world. Given the amount of accumulated data, it is important and urgent to perform a quantitative synthesis of these evidences. We therefore undertake a large meta-analysis of studies of the association between TCF7L2 gene polymorphisms and T2DM. To exclude heterogeneity of these studies, we stratified the data by Hardy-Weinberg equilibrium (HWE), publishing time, ethnicity, body mass index (BMI), diagnosis and age.

Although one meta-analysis [77] and several narrative reviews about TCF7L2 gene and T2DM had been reported [78–80], Succedent studies have now made more data available. Compared with the only one meta-analysis written by Cauchi and his colleagues [77] and published online in May, 2007, 28 more studies from 14 new papers [53, 55, 57, 63, 68, 69, 71, 73–76, 81–83]were available to us. And that, this meta-analysis of the association between TCF7L2 gene and T2DM, only summarized IVS3C>T polymorphism by just analyze allele frequency [77], the authors did not generalize the association between other important variations such as IVS4G>T, IVS3T>C and IVS4G>C polymorphisms and T2DM. Further, the only meta-analysis was based on published data, namely, without stratification, and, was thus limited in the interpretations they allowed. However, using the original data, we are able to explore more deeply, searching for sources of heterogeneity. Thus, we carried out a large meta-analysis to conduct a broad comprehensive and quantitative assessment on the association between variants of TCF7L2 gene and T2DM, and then, to perform a detailed examination of sources of heterogeneity. We sought to estimate magnitude of the genetic association as well as the genetic mode of action, and to determine the extent of heterogeneity in the strength of associations between different studies.

Methods

Search strategy and inclusion criteria

We conducted a comprehensive search of the US National Library of Medicine's PubMed, OMIM, ISI Web of Science, and Embase databases for all genetic association studies on at least one of the above four polymorphisms in the TCF7L2 gene and T2DM published from March, 2006, when the association between TCF7L2 gene and T2DM was first reported [2], through March, 2008. The search strategy was based on combination of the key words "transcription factor 7-like 2" or "TCF7L2" or "TCF7L2 gene polymorphism", and "diabetes" or "diabetes 2" or "type 2 diabetes" or "type 2 diabetes mellitus" or "T2D" or "T2DM". The references of all computer-identified publications were searched for additional studies, and the PubMed option "Related Articles" was used to search for potentially relevant papers. Reference lists in retrieved articles were also screened.

Searching was performed in duplicate by two independent reviewers (YZ and JYY). Without any language restriction (we will translate articles published in other languages into English by professional interpreters), we only selected published manuscripts (including their online supporting materials and the original data kindly sent by the authors) if they met the following criteria: 1) the publication was a population-based association study (family-based study design with linkage considerations was excluded), regardless of sample size, and 2) there were enough results for extraction of data, that is, the number of subjects with each allele and genotype in the T2DM and control groups. Where eligible papers had insufficient information, we contacted authors by e-mail or mail for additional information. In the case of sequential or multiple publications of analyses of the same data or overlapping data sets, the publication that reported data from the largest or most recent study was included.

Data extraction

Data were extracted independently and in duplicate by two investigators (YZ and JJY) who used the recommended guidelines for reporting on meta-analyses of observational studies [84]. The following data were extracted from the eligible studies: authors, journal and year of publication, country of origin, selection and characteristics of cases and controls, demographic data, ethnicity of the study population (Caucasian, North European, East Asian, African, Middle Eastern, American, or Pacific Asian), numbers of eligible and genotyped cases and controls, and genotype distributions in cases and controls and available subgroups thereof. Furthermore, we examined whether matching had been used; whether there was specific mention of blinding of the genotyping personnel to the clinical status of subjects; whether the genotyping method used had been validated; and whether genotype frequencies in control groups were in HWE. Any disagreement was adjudicated by a third author (YWZ).

Statistical Analysis

We used the odds ratio as the metric of choice and it was estimated for each study. In order to explore the possible association between TCF7L2 gene polymorphisms and T2MD and to avoid excessive comparisons, we calculated the odds ratio by two methods, namely, allele comparison (i.e., T allele vs. C allele in IVS3C>T polymorphism), and genotype comparison, that is, comparing the risk variant homozygotes and heterozygotes with wild homozygotes, respectively (i.e., TT vs. CC (OR 1) and TC vs. CC (OR 2) in IVS3C>T polymorphism; TT vs. GG (OR 3) and TG vs. GG (OR 4) in IVS4G>T polymorphism; CC vs. TT (OR 5) and CT vs. TT (OR 6) in IVS3T>C polymorphism; CC vs. TT (OR 7) and CT vs. TT (OR 8) in IVS4G>C polymorphism). We estimated and characterized the summary prevalence of risk allele by using only the data on controls. When we analyze genotype data for meta-analysis, cells with zero count had 0.5 been added. In addition, we calculated the population attributable risk (PAR) of risk allele according to method of Chang et al [71].

We first compared alleles for cases and controls in order to detect overall differences and genetic association. Allele frequencies were computed for studies reporting only genotypic data. Pooled odds ratio were computed three times, by the fixed-effects method of Mantel and Haenszel [85], by the random-effects method of DerSimonian and Laird [86], and by the Bayesian random-effects method of Warn et al [87]. Unless stated otherwise, random-effects estimates (DerSimonian and Laird) are reported here.

Our primary analysis for TCF7L2 gene polymorphisms and T2DM is based on genotype comparisons in order that magnitude of the genetic association and the genetic mode of action can be identified exactly. Once an overall gene effect was confirmed by the allele comparison, the genotype effects and genetic model were estimated by using the genetic model-free approach suggested by Minelli et al [88], where no assumptions about genetic models were required. The OR 1 to OR 8 were calculated by using multivariate meta-analysis with Bayesian method [88]. The logarithm (log) odds ratios were modeled on the base of both between and within study variations. Four separate stochastic lambda values (reflecting the genetic model), consisting of the ratios of logOR 2 versus logOR 1 for λ1, logOR 4 versus logOR 3 for λ2, logOR 6 versus logOR 5 for λ3, and logOR 8 versus logOR 7 for λ4, were calculated [88]. The parameter λ capture information about the genetic mode of action, briefly, the model is a recessive model if λ = 0, a codominant model if λ = 0.5, a dominant model if λ = 1, and homozygous or heterosis model if λ < 0 or λ > 1. Pooled odds ratio were also computed by the fixed-effects method of Mantel and Haenszel [85] and by the random-effects method of DerSimonian and Laird [86]. Unless stated otherwise, Bayesian estimates are reported here.

We tested for deviations from Hard-Weinberg equilibrium (HWE) in both case and control populations by using the exact method reported by Emigh et al [89] whether the author provide the P values of HWE or not. For above allele-based and genotype-base analyses, we took two methods to handle Hard-Weinberg disequilibrium. Firstly, we performed sensitivity analyses by including by including and excluding studies in Hard-Weinberg disequilibrium (case and control group were all considered). Secondly, we included all studies regardless of HWE and instead adjusted for the degree of disequilibrium by using the inbreeding coefficient (F) suggested by Trikalinos et al [90]. Briefly, the F was estimated for each study by using data in the control group. Predicted genotype frequencies were estimated [91] and used instead of the observed frequencies in the summary analysis of magnitude and genetic model.

In sensitivity analysis, we estimated between-study heterogeneity across all eligible comparisons using Cochran's Q statistic [92]. Heterogeneity was considered significant at P value less than 0.10 [92]. We also report I 2 statistic, which describes the percentage of variability in point estimates due to sample heterogeneity rather than sampling error [93, 94], and can quantify heterogeneity irrespective of the number of studies [94, 95]. As a guide, I 2 values less than 25 percent may be considered "low", values of 31 percent to 56 percent have been defined as indicating "low" to "moderate" heterogeneity, whereas values greater than 56 percent are considered indicative of "notable" heterogeneity, and "large" heterogeneity is stipulated for I 2 values of larger than or equal to 75 percent [94, 95]. The examination of heterogeneity was performed separately for two odds ratios (OR 1 and OR 2 for IVS3C>T polymorphism, OR 3 and OR 4 for IVS4G>T polymorphism, OR 6 and OR 7 for IVS3T>C polymorphism, OR 7 and OR 8 for IVS4G>C polymorphism) in genotype-based analysis. If there was heterogeneity on at least one of these two odds ratios, the source of the heterogeneity was explored by fitting a covariable (such as age, gender, ethnicity, BMI) in a meta-regression model [96–98]. Possible subgroup analyses were involved in ethnicity (Caucasians, North Europeans, East Asians, Indians and other racial groups), age (< 65 years vs. ≥ 65 years), diagnosis (normal diagnosis vs. no diagnosis or false diagnosis), publishing time (2006 vs. 2007 and 2008) and BMI (normal (< 30.0 kg/m2) vs. abnormal (≥ 30.0 kg/m2) according to the meta-regression analysis. Logistic regression analysis was used to compare odds ratios among these possible groups. In addition, genetic association studies may be especially susceptibility to the selective publication of positive findings on gene association [99]. In order to investigate the publication bias among studies, funnel plots [100] and cumulative meta-analysis [101] were used.

All analyses were conducted in Stata software, version 9.0 [102], using meta, metan, metabias, metacum, and metareg commands, except for the Bayesian method of genotype-based analysis. We fit the Bayesian models by using Markov chain Monte Carlo (MCMC) methods with Bayesian framework and perform our inferences using WinBUGS 1.4.3 [103], taking advantage of its flexibility as well as its ability to incorporate full uncertainty in all unknown parameters. Bayesian analyses yield credible intervals rather than confidence intervals; a 95% credible interval (CrI) describes a range in which the probability that the unknown quantity lies within this interval, after seeing the data, is 95%. A 'burn-in' of 10,000 iterations is carried out for models, followed by 50,000 iterations for parameter estimates. A P value less than 0.05 was considered statistically significant, except for tests of heterogeneity, where a level of 0.10 was used.

Results

Eligible studies

Totally eighty-nine papers were identified across Medline and Embase on the base of our search strategies. Two of the studies [57, 71]did not have genotypic data, but the authors friendly sent the supplementary information to us. Parra et al [68], reported four samples, that is, two Mexican American populations, one Spanish and one Nigerian population. With no controls of last three populations, we combined the two Mexican samples and then, discarded the Spanish and Nigerian samples in meta-analysis of polymorphism of IVS3C>T and IVS4G>T. Therefore, 25 eligible papers [2, 53, 55–71, 74, 75, 77, 81–83], including 36 studies (multiple studies performed in one paper in particular: 3 studies in paper of Grant et al [2], Horikoshi et al [62], Humphries et al [67], and Miyake et al [82]; 2 studies in paper of Scott et al [74], Sladek et al [75], and Cauchi et al [77]), examining the association between TCF7L2 gene polymorphisms and T2DM were identified by the inclusion criteria and all of them were written in English (Table 1, Additional file 1). The eligible studies for analysis involved in a total of 35,843 cases with T2DM and 39,123 controls. 35 studies from 24 papers [2, 53, 55, 57–71, 74, 75, 77, 81–83] contained data for the IVS3C>T polymorphism (Table 2, Additional file 1), 10 studies from 8 papers contained data for the IVS3T>C polymorphism (Table 3, Additional file 1), 21 papers involved in 29 studies contained data for the IVS4G>T polymorphism (Table 4, Additional file 1), and 9 papers including 15 samples contained data for IVS4G>C polymorphism (Table 5, Additional file 1).

There was considerable diversity of ethnic groups. Eligibility criteria for T2DM patients were shown in Table 1 (Additional file 1). Controls were mainly healthy populations who were described as nondiabetic and/or, Normoglycemic, and/or normal glucose tolerance (NGT), although varying details were presented regarding the extent of testing that had been done to excluded controls with impaired glucose tolerance and diabetes (table 1, Additional file 1). Remaining 3 investigations did not report the details of controls [60, 66, 67]. 1 study matched for geographic region [53]; 1 study matched for age, race, and BMI [56]; 5 studies matched for race [55, 67, 70, 83]; and one study 1 study matched for age, sex, and geographic region [58, 65]. The other investigations did not provide for the details of the matching information. Detailed characteristics of each investigation, along with P values for testing HWE, are listed in the Table 2–5 (Addional file 1) for IVS3C>T polymorphism, IVS3T>C polymorphism, IVS4G>T polymorphism and IVS4G>C polymorphism, respectively.

IVS3C>T polymorphism

The eligible studies for analysis of IVS3C>T polymorphism included a total of 33,135 cases with T2DM and 36,316 controls (Table 2, Additional file 1). Allele and genotype data were available for all the cases and controls in the eligible studies. There was no evidence of Hard-Weinberg disequilibrium in both cases and controls.

Bias diagnostics

We did not find any evidence of publication bias of the eligible studies. Funnel plots for the comparisons most commonly made, that is, for T allele vs. C allele, T homozygotes (TT) vs. C homozygotes (CC), and heterozygotes (TC) vs. C homozygotes (CC) are indicated in figure 1, 2 and 3. These three plots give corrected P = 0.410 (corrected z = 0.82), corrected P = 0.379 (corrected z = 0.88), and corrected P = 0.426 (corrected z = 0.80), respectively. Figure 4, 5, and 6 indicate the cumulative meta-analysis after each study from March 2006 to February 2008, for T allele vs. C allele, T homozygotes vs. C homozygotes, and heterozygotes vs. C homozygotes. The random-effects odds ratio of T homozygotes vs. C homozygotes waved slightly among studies performed in 2006 and has changed very little from around 2.0 after the eleventh study (figure 5). Cumulative meta-analysis of the other two comparisons by year of publication showed that the association between T2DM and TCF7L2 IVS3C>T polymorphism has remained significant and been consistent over time (figure 4 and figure 6).

Figure 1
figure 1

Funnel plots of rs7903146 polymorphism studies for T allele vs. C allele (corrected z = 0.82, corrected P = 0.410).

Figure 2
figure 2

Funnel plots of rs7903146 polymorphism studies for T homozygotes vs. C homozygotes (corrected z = 0.02, corrected P = 0.987).

Figure 3
figure 3

Funnel plots of rs7903146 polymorphism studies for heterozygotes vs. C homozygotes (corrected z = 0.80, corrected P = 0.426).

Figure 4
figure 4

Cumulative meta-analysis of the odds ratios of association between rs7903146 polymorphism of TCF7L2 gene and T2DM for T allele vs. C allele.

Figure 5
figure 5

Cumulative meta-analysis of the odds ratios of association between rs7903146 polymorphism of TCF7L2 gene and T2DM for T homozygotes vs. C homozygotes.

Figure 6
figure 6

Cumulative meta-analysis of the odds ratios of association between rs7903146 polymorphism of TCF7L2 gene and T2DM for heterozygotes vs. C homozygotes.

Allele comparisons

Data from control groups were used for calculating the summary allele frequency. The frequency of the risk allele T at TCF7L2 gene among controls was 20.2 percent (95% confidence interval (CI): 16.3, 24.1) and was considerable variation in different ethnic populations (for Caucasians, 29.0% (95% CI: 28.3, 29.7); for North Europeans, 21.9% (95% CI: 16.2, 27.6); for East Asians, 3.2 (95% CI: 2.6, 3.9); for Indians, 28.1% (95% CI: 25.9, 30.3); for others, 23.0% (95% CI: 9.4, 36.6)). We first contained all the 35 studies to evaluate the overall association between TCF7L2 IVS3C>T polymorphism and T2DM. With our allele-based analysis, the pooled odds ratio suggested a 1.416-fold increase in susceptibility to T2DM among persons with the T allele (random-effects OR TvC = 1.416, 95% CI: 1.346, 1.491), with a finding that was highly statistically significant (P = 0.000), but there was significant between-study heterogeneity (Q = 104.966, P = 0.000; I 2 = 67.6%) (Table 6 [Additional file 1], figure 7).

Figure 7
figure 7

Odds ratios of association between rs7903146 polymorphism of TCF7L2 gene and T2DM for T allele vs. C allele.

Genotype comparisons

The genotype frequency of TCF7L2 IVS3C>T polymorphism between case and control groups is characterized in web table. The genotype effects for TT versus CC (OR 1) and TC versus CC (OR 2) were calculated for each study (data are not shown and available from the first author on request). In our primary concern, multivariate meta-analysis was conducted to estimate the pooled effects and there was an significant increased risk of T2DM among populations with both homozygous variant TT genotype (Bayesian random-effect OR 1 = 1.968, 95 percent credible interval (CrI): 1.790, 2.157) and heterozygous variant TC genotype (Bayesian random-effect OR 2 = 1.406, 95 percent CrI: 1.341, 1.476), with a "moderate" between-study heterogeneity (Q = 66.102, P = 0.001, I 2 = 48.6%, figure 8) in the former and a "notable" between-study heterogeneity (Q = 84.028, P = 0.000, I 2 = 59.5%, figure 9) in the latter. The estimated stochastic parameter λ1 was 0.504 (95 percent CrI: 0.461, 0.551), which suggested a codominant genetic mode of action.

Figure 8
figure 8

Odds ratios of association between rs7903146 polymorphism of TCF7L2 gene and T2DM for T homozygotes vs. C homozygotes.

Figure 9
figure 9

Odds ratios of association between rs7903146 polymorphism of TCF7L2 gene and T2DM for heterozygotes vs. C homozygotes.

Examination of heterogeneity

Meta-regression was used to explore the source of heterogeneity and it was found that ethnicity was the only co-variable associated with both lnOR 1 (meta-regression beta coefficient (β) = -0.0993, P = 0.002) and lnOR 2 (β = -0.0360, P = 0.002). We then did a subgroup analysis by stratifying study samples into five ethnic groups, that is, Caucasians (11 samples), North European (6 samples), East Asian (9 samples), Indians (4 samples) and other racial descent group (5 samples), and found that this stratification can remove all the between-study heterogeneity among North Europeans, East Asians and Indians for both OR 1 and OR 2, rather than that among Caucasians and other racial group (Table 6, Additional file 1). We then did a further subgroup analysis of Caucasians by stratifying them into French Caucasians versus other Caucasians and extracted the studies with African populations for an independent analysis, and found that such analyses can only remove heterogeneity among these three populations for OR 1 rather than OR 2 (Table 6, Additional file 1). The pooled OR 1 for six ethnic groups by descending ranking, that is, French Caucasians, other Caucasians, North Europeans, Africans, Indians, and East Asians, was 3.051(95%CrI: 2.039, 4.906), 1.989 (95%CrI: 1.807, 2.216), 1.944 (95%CrI: 1.639, 2.284), 1.802 (95%CrI: 1.292, 2.480), 1.669 (95%CrI: 1.373. 2.104), and 1.595 (95%CrI: 0.537, 2.653), respectively, and the pooled OR 2 by descending ranking was 1.766 (95%CrI: 1.427, 2.278) for French Caucasians, 1.696 (95%CrI:1.258, 2.242) for Africans, 1.406 (95%CrI: 1.230, 1.617) for Indians, 1.357 (95%CrI: 1.255, 1.474) for North Europeans, 1.352 (95%CrI: 1.160, 1.591) for East Asians, and 1.349 (95%CrI: 1.278, 1.431) for other Caucasians (Table 2, Additional file 1). The estimated parameter λ1 by descending ranking was 0.981 (95 percent CrI: 0.759, 0.999) for Africans, 0.676 (95 percent CrI: 0.452, 0.957) for Indians, 0.670 (95 percent CrI: 0.278, 0.995) for East Asians, 0.511 (95 percent CrI: 0.449, 0.578) for French Caucasians, 0.472 (95 percent CrI: 0.345, 0.624) for North Europeans, and 0.436 (95 percent CrI: 0.366, 0.511) for other Caucasians, which suggested codominant genetic mode of action among these different ethnic populations except for Africans where dominant mode of action was found.

Although BMI was not associated with both lnOR 1 and lnOR 2 by separate meta-regression analysis, the significant association was found when we evaluate it together with ethnicity (for lnOR 1, β = -0.0512, P = 0.038; for lnOR 2, β = -0.0311, P = 0.004). In addition, some references reported that BMI may interact with TCF7L2 gene to increase the risk of T2MD. Therefore, we did a subgroup analysis for BMI (normal (< 30.0 kg/m2) vs. abnormal (≥ 30.0 kg/m2). However, such stratification cannot remove heterogeneity for both lnOR 1 (Q = 35.017, P = 0.006, I 2 = 52.1%) and lnOR 2 (Q = 28.046, P = 0.013, I 2 = 50.0%). The pooled OR 1, OR 2, and λ1 for populations with normal BMI were 2.017 (95%CrI: 1.711, 2.347), 1.470 (95%CrI: 1.350, 1.609), and 0.554 (95%CrI: 0.481, 0.637), respectively. The pooled OR 1, OR 2, and λ1 for populations with abnormal BMI were 1.931 (95%CrI: 1.670, 2.213), 1.337 (95%CrI: 1.250, 1.431), and 0.441 (95%CrI: 0.387, 0.500), respectively. Logistic regression analysis indicated that there was no significant difference between populations with normal and abnormal BMI for both lnOR 1 (P = 0.855) and lnOR 2 (P = 0.716).

IVS4G>T

Our searches identified 29 studies (totally 57,235 participants; 28,188 cases of T2DM, 29,047 controls) meeting our inclusion criteria. All of these studies examined the IVS4G>T polymorphism (Table 4, Additional file 1). Allele and genotype data were available for all the cases and controls in the eligible studies. There was no evidence of Hard-Weinberg disequilibrium in both cases and controls.

Bias diagnostics

There was no evidence of small study bias and publication bias for any comparisons. Funnel plots for the comparisons made for T allele vs. G allele, T homozygotes (TT) vs. G homozygotes (GG), and heterozygotes (TG) vs. G homozygotes are shown in figure 10, 11 and 12. These three funnel plots give corrected P = 0.339 (corrected z = 0.96), corrected P = 0.536 (corrected z = 0.62), and corrected P = 0.442 (corrected z = 0.77), respectively. Figure 13, 14 and 15 indicate the cumulative meta-analysis after each study from March 2006 to February 2008, for allele-based analysis and two genotype-based models of IVS4G>T polymorphism. The random-effects odds ratio of allele-based analysis was significant more than 1 and changed very little from around 1.379 after the first study (figure 13). The random-effects odds ratio of T homozygotes vs. G homozygotes waved slightly among studies performed in 2006 and has changed very little from around 1.884 after the eleventh study (figure 14). Cumulative meta-analysis of heterozygotes vs. G homozygotes by year of publication showed that the association between T2DM and TCF7L2 IVS4G>T polymorphism has remained significant and been consistent over time (figure 15).

Figure 10
figure 10

Funnel plots of rs12255372 polymorphism studies for T allele vs. G allele (corrected z = 0.96, corrected P = 0.339).

Figure 11
figure 11

Funnel plots of rs12255372 polymorphism studies for T homozygotes vs. G homozygotes (corrected z = 0.62, corrected P = 0.536).

Figure 12
figure 12

Funnel plots of rs12255372 polymorphism studies for heterozygotes vs. G homozygotes (corrected z = 0.77, corrected P = 0.442).

Figure 13
figure 13

Cumulative meta-analysis of the odds ratios of association between rs12255372 polymorphism of TCF7L2 gene and T2DM for T allele vs. H allele.

Figure 14
figure 14

Cumulative meta-analysis of the odds ratios of association between rs12255372 polymorphism of TCF7L2 gene and T2DM for T homozygotes vs. G homozygotes.

Figure 15
figure 15

Cumulative meta-analysis of the odds ratios of association between rs12255372 polymorphism of TCF7L2 gene and T2DM for heterozygotes vs. G homozygotes.

Allele comparisons

The frequency of the risk allele T at TCF7L2 gene among controls was 18.0 percent (95% CI: 13.5, 22.5) and was considerable variation in different ethnic populations (for Caucasians, 27.9% (95% CI: 27.0, 28.9); for North Europeans, 22.3% (95% CI: 12.6, 32.1); for East Asians, 2.1 (95% CI: 1.4, 2.7); and for Indians, 23.1% (95% CI: 17.8, 28.3)). We first contained all the 29 studies to evaluate the overall association between TCF7L2 IVS4G>T polymorphism and T2DM. With our allele-based analysis, the summary odds ratio suggested a 1.379-fold increase in susceptibility to T2DM among persons with the T allele (random-effects OR TvG = 1.379, 95% CI: 1.307, 1.454), with a finding that was highly statistically significant (P = 0.000), but there was significant between-study heterogeneity (Q = 66.191, P = 0.000; I 2 = 57.7%) (Table 6 [Additional file 1], figure 16).

Figure 16
figure 16

Odds ratios of association between rs12255372 polymorphism of TCF7L2 gene and T2DM for T allele vs. H allele.

Genotype comparisons

The genotype frequency of TCF7L2 IVS4G>T polymorphism between case and control groups is characterized in Table 4, Additional file 1. The genotype effects for TT versus GG (OR 3) and TG versus GG (OR 4) were calculated for each study (data are not shown and available from the first author on request). In our primary concern, multivariate meta-analysis was conducted to estimate the pooled effects and there was an significant increased risk of T2DM among populations with both homozygous variant TT genotype (Bayesian random-effect OR 2 = 1.885, 95 percent credible interval (CrI): 1.698, 2.088) and heterozygous variant TG genotype (Bayesian random-effect OR 4 = 1.360, 95 percent CrI: 1.291, 1.433), with a "moderate" between-study heterogeneity both in the former (Q = 46.515, P = 0.015, I 2 = 39.8%, figure 17) and in the latter (Q = 47.980, P = 0.011, I 2 = 41.6%, figure 18). The estimated stochastic parameter λ2 was 0.486 (95 percent CrI: 0.431, 0.545), which suggested a codominant genetic mode of action.

Figure 17
figure 17

Odds ratios of association between rs12255372 polymorphism of TCF7L2 gene and T2DM for T homozygotes vs. G homozygotes.

Figure 18
figure 18

Odds ratios of association between rs12255372 polymorphism of TCF7L2 gene and T2DM for heterozygotes vs. G homozygotes.

Examination of heterogeneity

Meta-regression was used to explore the source of heterogeneity and it was found that ethnicity was the only co-variable associated with both lnOR 3 (β = -0.1083, P = 0.005) and lnOR 4 (β = -0.0418, P = 0.047). We then did a subgroup analysis by stratifying study samples into five racial descent groups, that is, Caucasians (9 samples), North European (4 samples), East Asian (8 samples), Indians (4 samples) and other racial descent group (4 samples), and found that this stratification can remove all the between-study heterogeneity among North Europeans, East Asians and Indians for both OR 3 and OR 4, rather than that among Caucasians and other ethnic group (Table 6, Additional file 1). We then extracted the studies with other Caucasians (7 studies) for an independent analysis, and found that such analysis can well-done remove the heterogeneity among these studies for both OR 3 (Q = 8.407, P = 0.210, I 2 = 28.6%) and OR 4 (Q = 4.327, P = 0.633, I 2 = 0.0%) (Table 6, Additional file 1). We did not do a further analysis of other ethnic group because of different origin of populations in these studies. The pooled OR 3 for four ethnic groups (data of other ethnic group not shown for the heterogeneity of these studies) by descending ranking, that is, Caucasians, North Europeans, East Asians and Indians, was 2.091 (95%CrI: 1.755, 2.528; for other Caucasians, OR 3 = 1.860 (95%CrI: 1.672, 2.012)), 1.889 (95%CrI: 1.485, 2.412), 1.875 (95%CrI: 0.464~3.439), and 1.593 (95%CrI: 1.286, 2.086), respectively, and the pooled OR 4 by descending ranking was 1.569 (95%CrI: 1.304, 1.887) for East Asians, 1.413 (95%CrI: 1.299, 1.555; for other Caucasians, OR 4 = 1.326 (95%CrI: 1.251, 1.412)) for Caucasians, 1.380 (95%CrI: 1.222, 1.569) for North Europeans, and 1.358 (95%CrI: 1.176, 1.579) for Indians (Table 6, Additional file 1). The estimated parameter λ2 by descending ranking was 0.744 (95 percent CrI: 0.341, 1.996) for East Asians, 0.657 (95 percent CrI: 0.403, 0.982) for Indians, 0.507 (95 percent CrI: 0.347, 0.748) for North Europeans, and 0.470 (95 percent CrI: 0.410, 0.534) for Caucasians (for other Caucasians, λ2 = 0.455, 95 percent CrI: 0.376, 0.541), which suggested codominant genetic mode of action among all these four ethnic populations.

Although BMI was not associated with both lnOR 3 and lnOR 4 by separate meta-regression analysis, the significant association was found when we evaluate it together with racial descent (for lnOR 3, β = -0.0600, P = 0.029; for lnOR 4, β = -0.0467, P = 0.000). In addition, some references reported that BMI may interact with TCF7L2 gene to increase the risk of T2MD. Therefore, we also did a subgroup analysis for BMI in this polymorphism. However, such stratification can remove few heterogeneity for both lnOR 3 (Q = 22.517, P = 0.084, I 2 = 37.2%) and lnOR 4 (Q = 19.733, P = 0.042, I 2 = 47.0%). The pooled OR 3, OR 4, and λ2 for populations with normal BMI were 1.795 (95%CrI: 1.480, 2.211), 1.397 (95%CrI: 1.254, 1.567), and 0.572 (95%CrI: 0.448, 0.731), respectively. The pooled OR 3, OR 4, and λ2 for populations with abnormal BMI were 1.957 (95%CrI: 1.710, 2.210), 1.343 (95%CrI: 1.261, 1.429), and 0.442 (95%CrI: 0.377, 0.511), respectively. Logistic regression analysis indicated that there was no significant difference between populations with normal and abnormal BMI for both lnOR 1 (P = 0.069) and lnOR 2 (P = 0.152).

Other polymorphisms

The C allele of IVS3T>C polymorphism (10 studies with 15,718 participants (8,175 cases of T2DM, 7,543 controls); random-effects OR CvT = 1.323, 95% CI: 1.232, 1.419, P = 0.000; Q = 12.615, P = 0.181, I 2 = 28.7 percent) and the C allele of IVS4G>C polymorphism (15 studies with 21,161 participants (11,260 cases of T2DM, 9,901 controls); random-effects OR CvG = 1.238, 95 percent CI: 1.153, 1.330, P = 0.000; Q = 20.968, P = 0.102, I 2 = 33.2 percent) were significantly associated with T2DM (Table 6, Additional file 1). The pooled OR 5, OR 6, and λ3 for IVS3T>C polymorphism were 1.771 (95%CrI: 1.525, 2.062), 1.292 (95%CrI: 1.199, 1.402), and 0.450 (95%CrI: 0.330, 0.598), respectively. The pooled OR 7, OR 8, and λ4 for IVS4G>C polymorphism were 1.471 (95%CrI: 1.300, 1.706), 1.242 (95%CrI: 1.152, 1.355), and 0.561(95%CrI: 0.393, 0.777), respectively. These all suggested a codominant genetic mode of action among such two polymorphisms. We did not find any evidence of small sample bias and publication bias for any of the studies (IVS3T>C polymorphism: for OR 5, corrected z = 0.54, corrected P = 0.592; for OR 6, corrected z = 1.43, corrected P = 0.152; IVS4G>C polymorphism: for OR 7, corrected z = 0.40, corrected P = 0.692, for OR 8, corrected z = 1.39, corrected P = 0.166), and cumulative meta-analysis results for both allele-based and genotype-based analysis of the above two polymorphisms were stable (data and figure are not shown and available from the first author on request).

Discussion

Main findings

This is a large meta-analysis, including data from 25 papers involved in 36 genetic association studies with exactly 35,843 cases of T2DM and 39,123 controls, carefully avoiding the double-counting of participants in the study. The HuGE systematic review provides the most recent and comprehensive evaluation of the association between four TCF7L2 gene polymorphisms and susceptibility to T2DM. We find that near 70,000 subjects (33,135 cases of T2DM and 36,316 controls) were from 35 studies concerning the IVS3C>T polymorphism, and over 55,000 (28,188 cases of T2DM and 29,047 controls) subjects were from 29 studies investigating the IVS4G>T polymorphism. These two variants were the main study focus on the association between TCF7L2 gene polymorphisms and susceptibility to T2DM.

On one hand, the results indicates notable associations between two main TCF7L2 gene polymorphisms, namely, IVS3C>T as well as IVS4G>T, and T2DM. The magnitudes of this association were moderate, however, statistically significant. Our primary analysis finds that among IVS3C>T polymorphism, TC heterozygotes carry just over a 1.4-fold increased risk of T2MD, and TT homozygous variants carry near a 2.0-fold increase in T2MD risk when compared with CC homozygotes, and that among IVS4G>T polymorphism, TG heterozygotes carry near a 1.4-fold increased risk of T2MD, and TT homozygous variants carry approximate a 1.9-fold increase in T2MD risk when compared with GG homozygotes. Significant magnitudes of genetic effect for heterozygotes of these two variants were confirmed by the analysis on the pooled odds ratios among different ethnic groups with slight differences except for East Asians where higher pooled odds ratio was observed in IVS4G>T polymorphism when compared with that of the other ethnic populations, and genetic effect sizes for homozygotes of these two variants were also strengthened by the analysis on the pooled odds ratios among different ethnic groups including Caucasians, North Europeans, Africans and Indians by showing a small variation. However, we found that there was no significant association between such two homozygotes and disease among East Asians, and therefore, conflicting results have appeared, in which studies on Chinese provided negative results [57, 71], however, positive results were observed among Japanese populations [61, 62, 82]. Our meta-analysis on 9 studies of IVS3C>T polymorphism and 7 studies of IVS4G>T polymorphism among East Asians indicates a moderate statistical association and shows that the heterozygous variants contributed all the increased risk of T2MD among this population. This distinct disparity may suggest different mechanism of gene-disease between East Asians and other ethnic population. Moreover, although less samples and subjects were provided, the IVS3T>C and IVS4G>C polymorphisms are also significantly associated with T2DM. We then suggest that TCF7L2 is the most common susceptible gene for T2DM among various ethnic groups in the world.

On the other hand, we made no assumptions about genetic models and finds that the pooled odds ratios can clearly fit a multiplicative model and the lambda (λ) parameter can well-done indicate a codominant genetic mode of action with tight CIs among all these four TCF7L2 gene polymorphisms. We also explored the genetic mode among different ethnic populations and similar results were observed in Caucasians, North Europeans, East Asians and Indians except for Africans, where dominant genetic mode is suggested for IVS3C>T polymorphism. We cannot verify the genetic mode of action for Africans in the other three polymorphisms due to few studies. In addition, genetic mode of IVS3C>T variant among Africans also need a further verification by more studies. Even so, we nevertheless suggest a potential multiplicative genetic model for the four TCF7L2 gene polymorphisms.

The effect size of overall association between IVS3C>T variant and T2MD is near to 1.42, which is lower than that in a former meta-analysis about TCF7L2 and T2MD performed by Cauchi et al [77], and larger than that in a pooled analysis of three genome-wide association (GWA) study [73, 74, 76]. we supplement some studies about East Asians and Pima Indians [55, 57, 71, 88] in which risk homozygous variant TT was scanty, and excluded several studies without genotype frequency data, which may decrease the overall effect but provide a more comprehensive understanding of the association. Even so, we cannot find any significant difference of effect size between our study with the former meta-analysis as well as the pooled GWA studies. In addition, we computed the PAR of TCF7L2 according to IVS3C>T variant and the PAR for the combined genotypes TT and TC were 16.9, 23.2, 14.1, 2.5, 17.9, 27.0 for overall, Caucasians, North Europeans, East Asians, Indians, and Africans, respectively, suggesting this gene polymorphism may contribute near 1/5 of all T2MD in the globe except for East Asians.

Our findings were based on some gene-association studies and tens of thousands participants and were robust to each of the planned sensitivity analyses used. We cannot find any evidence of publication bias and small study bias by funnel plots and cumulative meta-analysis, but, considerable between-study heterogeneity was found. Between-study heterogeneity may be due to differences in sample selection (e.g., in age, sex, diagnosis, sample content, etc), or in methods (e.g., genotyping method), or it may be due to real differences in populations (e.g., in race), or in interactions with other risk factors (genetic or/and environmental factors). The results of primary analysis (genotype-based) on the four polymorphisms showed statistically significant between-study heterogeneity for both IVS3C>T polymorphism (OR 1 and OR 2), and IVS4G>C polymorphism (OR 3 and OR 4). In this Human Genome Epidemiology review, we examined five potential sources of between-study heterogeneity for genotype-based model by meta-regression and it was suggested that ethnicity was the only co-variable associated with the four ORs. In fact, more detailed manner of ethnic stratification, namely, French Caucasians, other Caucasians, North Europeans, East Asians, Indians, and Africans, can remove near all the heterogeneity of studies about both IVS3C>T and IVS4G>C polymorphisms for OR 1 to OR 4. Study results reported by several articles showed that BMI may cooperate with the TCF7L2 gene to increase the risk of the T2DM. However, disagreement of the results was published by different researchers. Cauchi et al[59], Horikoshi et al [62] and Humphries et al [67] separated to analyze the nonobese type 2 diabetic subjects (BMI<30 Kg/m2) and more higher allelic ORs or RRs were obtained in IVS3C>T and IVS4G>C polymorphisms; whereas, Miyake et al [82], Chandak et al[70], and Dahlgren et al [104] found that there was a slight increase or no change of odds ratio by adjustment of BMI. So we stratified the studied by BMI (normal vs. abnormal) and found that significant between-study heterogeneity nevertheless appeared in all the two subgroups. Both studies on IVS3C>T and IVS4G>C polymorphisms provided us for the similar results. But, we found a significant association when it was evaluated together with ethnicity. These results, on the one hand, suggested that BMI may be a possible factor to impact on the effect of TCF7L2 gene on T2DM by a indirect action; on the other hand, were difficult for us to explain the deep-seated reasons that how the BMI influence the effect of TCF7L2 gene on T2DM because we cannot gain the original data of BMI in each article. This is a general limitation of meta-analysis. So we hope the reasonable interpretation will be presented in the subsequent studies. Although sex [56, 63, 68, 71, 105], age of cases [63, 68, 71, 105], drugs[106], constitution[64, 69], and lifestyle [64, 69, 105] may regulate the effects of TCF7L2 gene, we cannot test the effects of age, sex, drugs and lifestyle on heterogeneity due to without sufficient related data being used. In addition, the age of the control NGT groups may be associated with the heterogeneity and the young controls may develop diabetes at a later date. However, there are also insufficient data.

Potential gene-gene and gene-environment interaction

T2DM seems to result from a complicated interplay of genetic and environmental factors influencing a number of intermediate traits of relevance to the diabetic phenotype (e.g., insulin secretion, insulin action, fat distribution, obesity). As a matter of fact, T2DM appears to be composed of subtypes where genetic susceptibility is strongly associated with environmental factors at one end of the spectrum, and highly genetic forms at the other end. Thus, there are several possible interactions between gene and gene or between gene and environmental factors.

Recent years, a series of new loci in some genes have been identified to contribute about 10~30% population attributable risk (PAR) of T2DM [73–76, 107]. But, the overall PAR seems to far less than the cumulative effect. So there must be some potential interactions among genes. The strongest known T2DM association (random-effects OR TvC≈1.42) was recently mapped to the transcription factor TCF7L2, a gene that is a target of the Wnt signaling pathway [50]. Whereafter, some other genes such as HHEX, IDE, DKK3 and KIF11, which are also in the Wnt signaling pathway [108–113] may have some interactions with TCF7L2 on the risk of T2DM. Cauchi et al [59]found that the expression of TCF7L2 gene was downregulated in obese subjects developing T2DM and they then brought forward a hypothesis of a potential interaction between TCF7L2 and Calpain-10 by a supporting material [114] that suggested that the Calpain system was involved in the constitutive regulation of β-catenin signaling functions. Wnt signaling has recently been shown to regulate pancreatic β-cell proliferation and the author suggested a possible interation between TCF7L2 gene and β-catenin gene [115]. However, there is not any reported publication about the interaction between TCF7L2 other related genes. These hypothetical interactions and their specific effects on T2DM, as well as functional analyses, will be required to further elucidate the role of variation in TCF7L2 in the pathogenesis of T2MD and very large samples are needed.

A variety of environmental factors can be implicated in the phenotype of T2DM, such as obesity, hypertension, bad lifestyle (smoking, drinking, high-energy diet), short of exercise, malnutrition, and some drugs [116]. All the above environmental factors may interact with the TCF7L2 gene to influence T2DM. So far, few specific gene-environment interactions have been described for TCF7L2 gene polymorphisms. Because obesity is a major determinant of development of T2DM, most patients with T2DM are obese when they develop diabetes, and obesity aggravates the insulin resistance. BMI and waist circumstance may become the main focus on attention of interactions, particularly in studies of IVS3C>T and IVS4G>C polymorphism among T2DM. Cauchi and colleagues [117] reported that IVS3C>T risk allele is more prevalent in non-obese than in obese diabetics. Humphries et al [67] found a statistically significant interaction between IVS3C>T as well as IVS4G>C genotype and BMI. Wang et al [118] replicated the result of Humphries et al. In succession, Helgason et al [47] discovered that two haplotypes, that is, HapA and HapBT2D had a interaction with BMI, and that the former was associated with increased BMI, whereas the latter was associated decreased BMI. Florez et al [106], observed that another index of obesity, waist circumstance, showed a nominally significant interation with both IVS3C>T and IVS4G>T in spite of no interaction between BMI and genotype. In addition, plasma triglyceride (TG) and C-reactive protein (CRP) have also been found to interact with genotype of both IVS3C>T and IVS4G>T [67]. These early findings suggestive of gene-environment interactions with the different TCF7L2 gene polymorphisms should be interpreted extremely cautiously, however, and much larger and more detailed studies are looked forward to substantiating such putative interactions with appropriate power and rigor. Complementarily, Smoking and drinking were important risk factors to influence on T2DM and they can also interact with genetic factor to increase the risk of this disease [116]. We found a higher proportion of current or/and past smokers and drinkers in the study populations. But, it was regretted for us that there was no any report about the interactions between TCF7L2 gene polymorphisms and smoking or/and drinking in all the open publications. This is a meaningful and valuable domain to explore.

Limitation of this meta-analysis

The lack of information especially the genotype data from some articles was the main limitation, we have done our possible to contact exhaustibly to the authors of publication that did not provide the original data and obtained full data on about 75 percent of subjects involved in various ethnic populations. Next, the quality of diagnosis, match, and genotyping varied a lot among studies, but the summary results were not changed when these questions were discussed carefully. There was considerable heterogeneity in IVS3C>T and IVS4G>T polymorphisms in our initial results. However, we discovered the main sources of the heterogeneity and were able to remove it finally. A question we cannot resolve is that how the four snps discussed here in strong LD will lead to differences in T2DM association among diverse populations. Besides, since environment factors and genes, such as obesity, smoking, drinking, and hypertension as well as plasma TG, and gene related to the Wnt signaling pathway are strong candidates for a potential interaction with TCF7L2 gene polymorphisms, the lack of data available to us on most of those environmental factors (except for obesity) and genes was a limitation, which we hope will be demonstrated by the following studies.

Conclusions and recommendation for future researches

This large meta-analysis summarizes the strong evidence for an association between TCF7L2 gene and T2DM both overall and in Caucasians, North Europeans, East Asians, Indians, and Africans, and suggested a potential multiplicative genetic model for all the four polymorphisms of TCF7L2 gene among different ethnic populations except for Africans, where additive genetic mode is suggested for IVS3C>T polymorphism. Our study results also suggest that IVS3C>T and IVS4G>T variants of TCF7L2 gene can be taken as reference loci for exploring T2DM susceptibility. We found and removed the main source of heterogeneity; however, we cannot find any evidence of bias. Furthermore, we estimated the potential gene-gene and gene-environmental interactions by which common variants in the TCF7L2 gene influence risk of T2DM. These two domains are just the recommendation for future researches.

Abbreviations

HuGE:

Human Genome Epidemiology

T2DM:

type 2 diabetes mellitus

TCF7L2:

Transcription factor 7-like 2

CI:

confidence interval

CrI:

credible interval

OR:

odds ratio

NGT:

normal glucose tolerant

BMI:

body mass index

HWE:

Hardy-Weinberg equilibrium.

References

  1. American Diabetes Association: Diagnosis and classification of diabetes mellitus. Diabetes Care. 2004, 27 (Suppl 1): S5-S10.

    Google Scholar 

  2. Grant SF, Thorleifsson G, Reynisdottir I, Benediktsson R, Manolescu A, Sainz J, Helgason A, Stefansson H, Emilsson V, Helgadottir A, Styrkarsdottir U, Magnusson KP, Walters GB, Palsdottir E, Jonsdottir T, Gudmundsdottir T, Gylfason A, Saemundsdottir J, Wilensky RL, Reilly MP, Rader DJ, Bagger Y, Christiansen C, Gudnason V, Sigurdsson G, Thorsteinsdottir U, Gulcher JR, Kong A, Stefansson K: Variant of transcription factor 7-like 2 (TCF7L2) gene confers risk of type 2 diabetes. Nat Genet. 2006, 38 (3): 320-323.

    CAS  PubMed  Google Scholar 

  3. King H, Aubert RE, Herman WH: Global burden of diabetes, 1995–2025: prevalence, numerical estimates, and projections. Diabetes Care. 1998, 21 (9): 1414-1431.

    CAS  PubMed  Google Scholar 

  4. Zimmet P, Alberti KG, Shaw J: Global and societal implications of the diabetes epidemic. Nature. 2001, 414 (6865): 782-787.

    CAS  PubMed  Google Scholar 

  5. Ceesay MM, Morgan MW, Kamanda MO, Willoughby VR, Lisk DR: Prevalence of diabetes in rural and urban populations in southern Sierra Leone: a preliminary survey. Trop Med Int Health. 1997, 2 (3): 272-277.

    CAS  PubMed  Google Scholar 

  6. Ducorps M, Baleynaud S, Mayaudon H, Castagne C, Bauduceau B: A prevalence survey of diabetes in Mauritania. Diabetes Care. 1996, 19 (7): 761-763.

    CAS  PubMed  Google Scholar 

  7. Elbagir MN, Eltom MA, Elmahadi EM, Kadam IM, Berne C: A population-based study of the prevalence of diabetes and impaired glucose tolerance in adults in northern Sudan. Diabetes Care. 1996, 19 (10): 1126-1128.

    CAS  PubMed  Google Scholar 

  8. Owoaje EE, Rotimi CN, Kaufman JS, Tracy J, Cooper RS: Prevalence of adult diabetes in Ibadan, Nigeria. East Afr Med J. 1997, 74 (5): 299-302.

    CAS  PubMed  Google Scholar 

  9. Soderberg S, Zimmet P, Tuomilehto J, de Courten M, Dowse GK, Chitson P, Gareeboo H, Alberti KG, Shaw JE: Increasing prevalence of Type 2 diabetes mellitus in all ethnic groups in Mauritius. Diabet Med. 2005, 22 (1): 61-68.

    CAS  PubMed  Google Scholar 

  10. Sande van der MA, Bailey R, Faal H, Banya WA, Dolin P, Nyan OA, Ceesay SM, Walraven GE, Johnson GJ, McAdam KP: Nationwide prevalence study of hypertension and related non-communicable diseases in The Gambia. Trop Med Int Health. 1997, 2 (11): 1039-1048.

    PubMed  Google Scholar 

  11. Malek R, Belateche F, Laouamri S, Hamdi-Cherif M, Touabti A, Bendib W, Nechadi A, Mekideche FZ, Hanat S: [Prevalence of type 2 diabetes mellitus and glucose intolerance in the Setif area (Algeria)]. Diabetes Metab. 2001, 27 (2 Pt 1): 164-171.

    CAS  PubMed  Google Scholar 

  12. Abu Sayeed M, Mahtab H, Akter Khanam P, Abul Ahsan K, Banu A, Rashid AN, Azad Khan AK: Diabetes and impaired fasting glycemia in the tribes of Khagrachari hill tracts of Bangladesh. Diabetes Care. 2004, 27 (5): 1054-1059.

    CAS  PubMed  Google Scholar 

  13. Chang C, Lu F, Yang YC, Wu JS, Wu TJ, Chen MS, Chuang LM, Tai TY: Epidemiologic study of type 2 diabetes in Taiwan. Diabetes Res Clin Pract. 2000, 50 (Suppl 2): S49-59.

    PubMed  Google Scholar 

  14. Hussain A, Rahim MA, Azad Khan AK, Ali SM, Vaaler S: Type 2 diabetes in rural and urban population: diverse prevalence and associated risk factors in Bangladesh. Diabet Med. 2005, 22 (7): 931-936.

    CAS  PubMed  Google Scholar 

  15. Karki P, Baral N, Lamsal M, Rijal S, Koner BC, Dhungel S, Koirala S: Prevalence of non-insulin dependent diabetes mellitus in urban areas of eastern Nepal: a hospital based study. Southeast Asian J Trop Med Public Health. 2000, 31 (1): 163-166.

    CAS  PubMed  Google Scholar 

  16. Kim SM, Lee JS, Lee J, Na JK, Han JH, Yoon DK, Baik SH, Choi DS, Choi KM: Prevalence of diabetes and impaired fasting glucose in Korea: Korean National Health and Nutrition Survey 2001. Diabetes Care. 2006, 29 (2): 226-231.

    CAS  PubMed  Google Scholar 

  17. Kuzuya T: Prevalence of diabetes mellitus in Japan compiled from literature. Diabetes Res Clin Pract. 1994, 24 (Suppl): S15-21.

    PubMed  Google Scholar 

  18. Lee WR: The changing demography of diabetes mellitus in Singapore. Diabetes Res Clin Pract. 2000, 50 (Suppl 2): S35-39.

    PubMed  Google Scholar 

  19. Ramachandran A, Snehalatha C, Latha E, Vijay V, Viswanathan M: Rising prevalence of NIDDM in an urban population in India. Diabetologia. 1997, 40 (2): 232-237.

    CAS  PubMed  Google Scholar 

  20. Glatthaar C, Welborn TA, Stenhouse NS, Garcia-Webb P: Diabetes and impaired glucose tolerance. A prevalence estimate based on the Busselton 1981 survey. Med J Aust. 1985, 143 (10): 436-440.

    CAS  PubMed  Google Scholar 

  21. Andersson DK, Svardsudd K, Tibblin G: Prevalence and incidence of diabetes in a Swedish community 1972–1987. Diabet Med. 1991, 8 (5): 428-434.

    CAS  PubMed  Google Scholar 

  22. Baan CA, Feskens EJ: [Prevention of diabetes mellitus type 2]. Ned Tijdschr Geneeskd. 2001, 145 (35): 1677-1680.

    CAS  PubMed  Google Scholar 

  23. Botas Cervero P, Delgado Alvarez E, Castano Fernandez G, Diaz De Grenu C, Prieto Santiago J, Diaz Cadorniga FJ: [Prevalence of diabetes mellitus and glucose intolerance in the population aged 30 to 75 years in Asturias, Spain]. Rev Clin Esp. 2002, 202 (8): 421-429.

    CAS  PubMed  Google Scholar 

  24. Brohall G, Behre CJ, Hulthe J, Wikstrand J, Fagerberg B: Prevalence of diabetes and impaired glucose tolerance in 64-year-old Swedish women: experiences of using repeated oral glucose tolerance tests. Diabetes Care. 2006, 29 (2): 363-367.

    PubMed  Google Scholar 

  25. Detournay B, Vauzelle-Kervroedan F, Charles MA, Forhan A, Fagnani F, Fender P, Eschwege E: [Epidemiology, management and costs of type 2 diabetes in France in 1998]. Diabetes Metab. 1999, 25 (4): 356-365.

    CAS  PubMed  Google Scholar 

  26. Fabian W, Majkowska L, Stefanski A, Moleda P: [Prevalence of diabetes, antidiabetic treatment and chronic diabetic complications reported by general practitioners]. Przegl Lek. 2005, 62 (4): 201-205.

    PubMed  Google Scholar 

  27. Garancini MP, Calori G, Ruotolo G, Manara E, Izzo A, Ebbli E, Bozzetti AM, Boari L, Lazzari P, Gallus G: Prevalence of NIDDM and impaired glucose tolerance in Italy: an OGTT-based population study. Diabetologia. 1995, 38 (3): 306-313.

    CAS  PubMed  Google Scholar 

  28. Gatling W, Budd S, Walters D, Mullee MA, Goddard JR, Hill RD: Evidence of an increasing prevalence of diagnosed diabetes mellitus in the Poole area from 1983 to 1996. Diabet Med. 1998, 15 (12): 1015-1021.

    CAS  PubMed  Google Scholar 

  29. Glumer C, Jorgensen T, Borch-Johnsen K: Prevalences of diabetes and impaired glucose regulation in a Danish population: the Inter99 study. Diabetes Care. 2003, 26 (8): 2335-2340.

    PubMed  Google Scholar 

  30. Gourdy P, Ruidavets JB, Ferrieres J, Ducimetiere P, Amouyel P, Arveiler D, Cottel D, Lamamy N, Bingham A, Hanaire-Broutin H: Prevalence of type 2 diabetes and impaired fasting glucose in the middle-aged population of three French regions – The MONICA study 1995–97. Diabetes Metab. 2001, 27 (3): 347-358.

    CAS  PubMed  Google Scholar 

  31. Kuraeva TL, Sergeev AS, Lebedev NB, Babadzhanova G, Kerimi NB: [Diabetes mellitus incidence and its prevalence in Moscow]. Probl Endokrinol (Mosk). 1993, 39 (6): 4-7.

    CAS  Google Scholar 

  32. Pilotto L, Gaggioli A, Lo Noce C, Dima F, Palmieri L, Uguccioni M, Pede S, Giampaoli S, Vanuzzo D: [Diabetes in Italy: a public health problem]. Ital Heart J Suppl. 2004, 5 (6): 480-486.

    PubMed  Google Scholar 

  33. Abdella N, Al Arouj M, Al Nakhi A, Al Assoussi A, Moussa M: Non-insulin-dependent diabetes in Kuwait: prevalence rates and associated risk factors. Diabetes Res Clin Pract. 1998, 42 (3): 187-196.

    CAS  PubMed  Google Scholar 

  34. Al-Habori M, Al-Mamari M, Al-Meeri A: Type II Diabetes Mellitus and impaired glucose tolerance in Yemen: prevalence, associated metabolic changes and risk factors. Diabetes Res Clin Pract. 2004, 65 (3): 275-281.

    CAS  PubMed  Google Scholar 

  35. Al-Nuaim AR: Prevalence of glucose intolerance in urban and rural communities in Saudi Arabia. Diabet Med. 1997, 14 (7): 595-602.

    CAS  PubMed  Google Scholar 

  36. Husseini A, Abdul-Rahim H, Awartani F, Giacaman R, Jervell J, Bjertness E: Type 2 diabetes mellitus, impaired glucose tolerance and associated factors in a rural Palestinian village. Diabet Med. 2000, 17 (10): 746-748.

    CAS  PubMed  Google Scholar 

  37. Blanchard JF, Dean H, Anderson K, Wajda A, Ludwig S, Depew N: Incidence and prevalence of diabetes in children aged 0–14 years in Manitoba, Canada, 1985–1993. Diabetes Care. 1997, 20 (4): 512-515.

    CAS  PubMed  Google Scholar 

  38. Carter J, Horowitz R, Wilson R, Sava S, Sinnock P, Gohdes D: Tribal differences in diabetes: prevalence among American Indians in New Mexico. Public Health Rep. 1989, 104 (6): 665-669.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Castro-Sanchez H, Escobedo-de la Pena J: [Prevalence of non insulin dependent diabetes mellitus and associated risk factors in the Mazatec population of the State of Oaxaca, Mexico]. Gac Med Mex. 1997, 133 (6): 527-534.

    CAS  PubMed  Google Scholar 

  40. Malerbi DA, Franco LJ: Multicenter study of the prevalence of diabetes mellitus and impaired glucose tolerance in the urban Brazilian population aged 30–69 yr. The Brazilian Cooperative Group on the Study of Diabetes Prevalence. Diabetes Care. 1992, 15 (11): 1509-1516.

    CAS  PubMed  Google Scholar 

  41. Perez-Bravo F, Carrasco E, Santos JL, Calvillan M, Larenas G, Albala C: Prevalence of type 2 diabetes and obesity in rural Mapuche population from Chile. Nutrition. 2001, 17 (3): 236-238.

    CAS  PubMed  Google Scholar 

  42. Torquato MT, Montenegro Junior RM, Viana LA, de Souza RA, Lanna CM, Lucas JC, Bidurin C, Foss MC: Prevalence of diabetes mellitus and impaired glucose tolerance in the urban population aged 30–69 years in Ribeirao Preto (Sao Paulo), Brazil. Sao Paulo Med J. 2003, 121 (6): 224-230.

    PubMed  Google Scholar 

  43. Carrasco EP, Perez FB, Angel BB, Albala CB, Santos JL, Larenas GY, Montalvo DV: [Prevalence of type 2 diabetes and obesity in two Chilean aboriginal populations living in urban zones]. Rev Med Chil. 2004, 132 (10): 1189-1197.

    PubMed  Google Scholar 

  44. Ramachandran A, Snehalatha C, Kapur A, Vijay V, Mohan V, Das AK, Rao PV, Yajnik CS, Prasanna Kumar KM, Nair JD: High prevalence of diabetes and impaired glucose tolerance in India: National Urban Diabetes Survey. Diabetologia. 2001, 44 (9): 1094-1101.

    CAS  PubMed  Google Scholar 

  45. Wang K, Li T, Xiang H: [Study on the epidemiological characteristics of diabetes mellitus and IGT in China]. Zhonghua Liu Xing Bing Xue Za Zhi. 1998, 19 (5): 282-285.

    CAS  PubMed  Google Scholar 

  46. Duval A, Busson-Leconiat M, Berger R, Hamelin R: Assignment of the TCF-4 gene (TCF7L2) to human chromosome band 10q25.3. Cytogenet Cell Genet. 2000, 88 (3–4): 264-265.

    CAS  PubMed  Google Scholar 

  47. Helgason A, Palsson S, Thorleifsson G, Grant SF, Emilsson V, Gunnarsdottir S, Adeyemo A, Chen Y, Chen Y, Reynisdottir I, Benediktsson R, Hinney A, Hansen T, Andersen G, Borch-Johnsen K, Jorgensen T, Schafer H, Faruque M, Doumatey A, Zhou J, Wilensky RL, Reilly MP, Rader DJ, Bagger Y, Christiansen C, Sigurdsson G, Hebebrand J, Pedersen O, Thorsteinsdottir U, Gulcher JR, Kong A, Rotimi C, Stefansson K: Refining the impact of TCF7L2 gene variants on type 2 diabetes and adaptive evolution. Nat Genet. 2007, 39 (2): 218-225.

    CAS  PubMed  Google Scholar 

  48. Papadopoulou S, Edlund H: Attenuated Wnt signaling perturbs pancreatic growth but not pancreatic function. Diabetes. 2005, 54 (10): 2844-2851.

    CAS  PubMed  Google Scholar 

  49. Prunier C, Hocevar BA, Howe PH: Wnt signaling: physiology and pathology. Growth Factors. 2004, 22 (3): 141-150.

    CAS  PubMed  Google Scholar 

  50. Yi F, Brubaker PL, Jin T: TCF-4 mediates cell type-specific regulation of proglucagon gene expression by beta-catenin and glycogen synthase kinase-3beta. J Biol Chem. 2005, 280 (2): 1457-1464.

    CAS  PubMed  Google Scholar 

  51. Nelson WJ, Nusse R: Convergence of Wnt, beta-catenin, and cadherin pathways. Science. 2004, 303 (5663): 1483-1487.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Duval A, Rolland S, Tubacher E, Bui H, Thomas G, Hamelin R: The human T-cell transcription factor-4 gene: structure, extensive characterization of alternative splicings, and mutational analysis in colorectal cancer cell lines. Cancer Res. 2000, 60 (14): 3872-3879.

    CAS  PubMed  Google Scholar 

  53. Sale MM, Smith SG, Mychaleckyj JC, Keene KL, Langefeld CD, Leak TS, Hicks PJ, Bowden DW, Rich SS, Freedman BI: Variants of the transcription factor 7-like 2 (TCF7L2) gene are associated with type 2 diabetes in an African-American population enriched for nephropathy. Diabetes. 2007, 56 (10): 2638-2642.

    CAS  PubMed  Google Scholar 

  54. Saxena R, Gianniny L, Burtt NP, Lyssenko V, Giuducci C, Sjogren M, Florez JC, Almgren P, Isomaa B, Orho-Melander M, Lindblad U, Daly MJ, Tuomi T, Hirschhorn JN, Ardlie KG, Groop LC, Altshuler D: Common single nucleotide polymorphisms in TCF7L2 are reproducibly associated with type 2 diabetes and reduce the insulin response to glucose in nondiabetic individuals. Diabetes. 2006, 55 (10): 2890-2895.

    CAS  PubMed  Google Scholar 

  55. Guo T, Hanson RL, Traurig M, Muller YL, Ma L, Mack J, Kobes S, Knowler WC, Bogardus C, Baier LJ: TCF7L2 is not a major susceptibility gene for type 2 diabetes in Pima Indians: analysis of 3,501 individuals. Diabetes. 2007, 56 (12): 3082-3088.

    CAS  PubMed  Google Scholar 

  56. Zhang C, Qi L, Hunter DJ, Meigs JB, Manson JE, van Dam RM, Hu FB: Variant of transcription factor 7-like 2 (TCF7L2) gene and the risk of type 2 diabetes in large cohorts of U.S. women and men. Diabetes. 2006, 55 (9): 2645-2648.

    CAS  PubMed  Google Scholar 

  57. Ng MC, Tam CH, Lam VK, So WY, Ma RC, Chan JC: Replication and identification of novel variants at TCF7L2 associated with type 2 diabetes in Hong Kong Chinese. J Clin Endocrinol Metab. 2007, 92 (9): 3733-3737.

    CAS  PubMed  Google Scholar 

  58. Scott LJ, Bonnycastle LL, Willer CJ, Sprau AG, Jackson AU, Narisu N, Duren WL, Chines PS, Stringham HM, Erdos MR, Valle TT, Tuomilehto J, Bergman RN, Mohlke KL, Collins FS, Boehnke M: Association of transcription factor 7-like 2 (TCF7L2) variants with type 2 diabetes in a Finnish sample. Diabetes. 2006, 55 (9): 2649-2653.

    CAS  PubMed  Google Scholar 

  59. Cauchi S, Meyre D, Dina C, Choquet H, Samson C, Gallina S, Balkau B, Charpentier G, Pattou F, Stetsyuk V, Scharfmann R, Staels B, Fruhbeck G, Froguel P: Transcription factor TCF7L2 genetic study in the French population: expression in human beta-cells and adipose tissue and strong association with type 2 diabetes. Diabetes. 2006, 55 (10): 2903-2908.

    CAS  PubMed  Google Scholar 

  60. Groves CJ, Zeggini E, Minton J, Frayling TM, Weedon MN, Rayner NW, Hitman GA, Walker M, Wiltshire S, Hattersley AT, McCarthy MI: Association analysis of 6,736 U.K. subjects provides replication and confirms TCF7L2 as a type 2 diabetes susceptibility gene with a substantial effect on individual risk. Diabetes. 2006, 55 (9): 2640-2644.

    CAS  PubMed  Google Scholar 

  61. Hayashi T, Iwamoto Y, Kaku K, Hirose H, Maeda S: Replication study for the association of TCF7L2 with susceptibility to type 2 diabetes in a Japanese population. Diabetologia. 2007, 50 (5): 980-984.

    CAS  PubMed  Google Scholar 

  62. Horikoshi M, Hara K, Ito C, Nagai R, Froguel P, Kadowaki T: A genetic variation of the transcription factor 7-like 2 gene is associated with risk of type 2 diabetes in the Japanese population. Diabetologia. 2007, 50 (4): 747-751.

    CAS  PubMed  Google Scholar 

  63. Kimber CH, Doney AS, Pearson ER, McCarthy MI, Hattersley AT, Leese GP, Morris AD, Palmer CN: TCF7L2 in the Go-DARTS study: evidence for a gene dose effect on both diabetes susceptibility and control of glucose levels. Diabetologia. 2007, 50 (6): 1186-1191.

    CAS  PubMed  Google Scholar 

  64. Marzi C, Huth C, Kolz M, Grallert H, Meisinger C, Wichmann HE, Rathmann W, Herder C, Illig T: Variants of the transcription factor 7-like 2 gene (TCF7L2) are strongly associated with type 2 diabetes but not with the metabolic syndrome in the MONICA/KORA surveys. Horm Metab Res. 2007, 39 (1): 46-52.

    CAS  PubMed  Google Scholar 

  65. Mayans S, Lackovic K, Lindgren P, Ruikka K, Agren A, Eliasson M, Holmberg D: TCF7L2 polymorphisms are associated with type 2 diabetes in northern Sweden. Eur J Hum Genet. 2007, 15 (3): 342-346.

    CAS  PubMed  Google Scholar 

  66. van Vliet-Ostaptchouk JV, Shiri-Sverdlov R, Zhernakova A, Strengman E, van Haeften TW, Hofker MH, Wijmenga C: Association of variants of transcription factor 7-like 2 (TCF7L2) with susceptibility to type 2 diabetes in the Dutch Breda cohort. Diabetologia. 2007, 50 (1): 59-62.

    CAS  PubMed  Google Scholar 

  67. Humphries SE, Gable D, Cooper JA, Ireland H, Stephens JW, Hurel SJ, Li KW, Palmen J, Miller MA, Cappuccio FP, Elkeles R, Godsland I, Miller GJ, Talmud PJ: Common variants in the TCF7L2 gene and predisposition to type 2 diabetes in UK European Whites, Indian Asians and Afro-Caribbean men and women. J Mol Med. 2006, 84 (12): 1005-1014.

    CAS  PubMed  Google Scholar 

  68. Parra EJ, Cameron E, Simmonds L, Valladares A, McKeigue P, Shriver M, Wacher N, Kumate J, Kittles R, Cruz M: Association of TCF7L2 polymorphisms with type 2 diabetes in Mexico City. Clin Genet. 2007, 71 (4): 359-366.

    CAS  PubMed  Google Scholar 

  69. Bodhini D, Radha V, Dhar M, Narayani N, Mohan V: The rs12255372(G/T) and rs7903146(C/T) polymorphisms of the TCF7L2 gene are associated with type 2 diabetes mellitus in Asian Indians. Metabolism. 2007, 56 (9): 1174-1178.

    CAS  PubMed  Google Scholar 

  70. Chandak GR, Janipalli CS, Bhaskar S, Kulkarni SR, Mohankrishna P, Hattersley AT, Frayling TM, Yajnik CS: Common variants in the TCF7L2 gene are strongly associated with type 2 diabetes mellitus in the Indian population. Diabetologia. 2007, 50 (1): 63-67.

    CAS  PubMed  Google Scholar 

  71. Chang YC, Chang TJ, Jiang YD, Kuo SS, Lee KC, Chiu KC, Chuang LM: Association study of the genetic polymorphisms of the transcription factor 7-like 2 (TCF7L2) gene and type 2 diabetes in the Chinese population. Diabetes. 2007, 56 (10): 2631-2637.

    CAS  PubMed  Google Scholar 

  72. McCarthy MI, Zeggini E: Genetics of type 2 diabetes. Curr Diab Rep. 2006, 6 (2): 147-154.

    CAS  PubMed  Google Scholar 

  73. Saxena RVB, Lyssenko V, Burtt NP, de Bakker PI, Chen H, Roix JJ, Kathiresan S, Hirschhorn JN, Daly MJ, Hughes TE, Groop L, Altshuler D, Almgren P, Florez JC, Meyer J, Ardlie K, Bengtsson Bostrom KIB, Lettre G, Lindblad U, Lyon HN, Melander O, Newton-Cheh C, Nilsson P, Orho-Melander M, Rastam L, Speliotes EK, Taskinen MR, Tuomi T, Guiducci C, Berglund A, Carlson J, Gianniny L, Hackett R, Hall L, Holmkvist J, Laurila E, et al: Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science. 2007, 316 (5829): 1331-1336.

    CAS  PubMed  Google Scholar 

  74. Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, Erdos MR, Stringham HM, Chines PS, Jackson AU, Prokunina-Olsson L, Ding CJ, Swift AJ, Narisu N, Hu T, Pruim R, Xiao R, Li XY, Conneely KN, Riebow NL, Sprau AG, Tong M, White PP, Hetrick KN, Barnhart MW, Bark CW, Goldstein JL, Watkins L, Xiang F, Saramies J, Buchanan TA, Watanabe RM, Valle TT, Kinnunen L, Abecasis GR, Pugh EW, Doheny KF, Bergman RN, Tuomilehto J, Collins FS, Boehnke M: A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science. 2007, 316 (5829): 1341-1345.

    CAS  PubMed  PubMed Central  Google Scholar 

  75. Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, Boutin P, Vincent D, Belisle A, Hadjadj S, Balkau B, Heude B, Charpentier G, Hudson TJ, Montpetit A, Pshezhetsky AV, Prentki M, Posner BI, Balding DJ, Meyre D, Polychronakos C, Froguel P: A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 2007, 445 (7130): 881-885.

    CAS  PubMed  Google Scholar 

  76. Zeggini E, Weedon MN, Lindgren CM, Frayling TM, Elliott KS, Lango H, Timpson NJ, Perry JR, Rayner NW, Freathy RM, Barrett JC, Shields B, Morris AP, Ellard S, Groves CJ, Harries LW, Marchini JL, Owen KR, Knight B, Cardon LR, Walker M, Hitman GA, Morris AD, Doney AS, McCarthy MI, Hattersley AT: Replication of genome-wide association signals in UK samples reveals risk loci for type 2 diabetes. Science. 2007, 316 (5829): 1336-1341.

    CAS  PubMed  PubMed Central  Google Scholar 

  77. Cauchi S, El Achhab Y, Choquet H, Dina C, Krempler F, Weitgasser R, Nejjari C, Patsch W, Chikri M, Meyre D, Froguel P: TCF7L2 is reproducibly associated with type 2 diabetes in various ethnic groups: a global meta-analysis. J Mol Med. 2007, 85 (7): 777-782.

    CAS  PubMed  Google Scholar 

  78. Florez JC: The new type 2 diabetes gene TCF7L2. Curr Opin Clin Nutr Metab Care. 2007, 10 (4): 391-396.

    CAS  PubMed  Google Scholar 

  79. Frayling TM: A new era in finding Type 2 diabetes genes-the unusual suspects. Diabet Med. 2007, 24 (7): 696-701.

    CAS  PubMed  Google Scholar 

  80. Weedon MN: The importance of TCF7L2. Diabet Med. 2007, 24 (10): 1062-1066.

    CAS  PubMed  Google Scholar 

  81. De Silva NM, Steele A, Shields B, Knight B, Parnell K, Weedon MN, Hattersley AT, Frayling TM: The transcription factor 7-like 2 (TCF7L2) gene is associated with Type 2 diabetes in UK community-based cases, but the risk allele frequency is reduced compared with UK cases selected for genetic studies. Diabet Med. 2007, 24 (10): 1067-1072.

    CAS  PubMed  Google Scholar 

  82. Miyake K, Horikawa Y, Hara K, Yasuda K, Osawa H, Furuta H, Hirota Y, Yamagata K, Hinokio Y, Oka Y, Iwasaki N, Iwamoto Y, Yamada Y, Seino Y, Maegawa H, Kashiwagi A, Yamamoto K, Tokunaga K, Takeda J, Makino H, Nanjo K, Kadowaki T, Kasuga M: Association of TCF7L2 polymorphisms with susceptibility to type 2 diabetes in 4,087 Japanese subjects. J Hum Genet. 2007

    Google Scholar 

  83. Rees SD, Bellary S, Britten AC, O'Hare JP, Kumar S, Barnett AH, Kelly MA: Common variants of the TCF7L2 gene are associated with increased risk of type 2 diabetes mellitus in a UK-resident South Asian population. BMC Med Genet. 2008, 9 (1): 8-

    PubMed  PubMed Central  Google Scholar 

  84. Stroup DF, Berlin JA, Morton SC, Olkin I, Williamson GD, Rennie D, Moher D, Becker BJ, Sipe TA, Thacker SB: Meta-analysis of observational studies in epidemiology: a proposal for reporting. Meta-analysis Of Observational Studies in Epidemiology (MOOSE) group. Jama. 2000, 283 (15): 2008-2012.

    CAS  PubMed  Google Scholar 

  85. Mantel N, Haenszel W: Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst. 1959, 22 (4): 719-748.

    CAS  PubMed  Google Scholar 

  86. DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7 (3): 177-188.

    CAS  PubMed  Google Scholar 

  87. Warn DE, Thompson SG, Spiegelhalter DJ: Bayesian random effects meta-analysis of trials with binary outcomes: methods for the absolute risk difference and relative risk scales. Stat Med. 2002, 21 (11): 1601-1623.

    CAS  PubMed  Google Scholar 

  88. Minelli C, Thompson JR, Abrams KR, Lambert PC: Bayesian implementation of a genetic model-free approach to the meta-analysis of genetic association studies. Stat Med. 2005, 24 (24): 3845-3861.

    PubMed  Google Scholar 

  89. Emigh T: A comparison of tests for Hardy-Weinberg equilibrium. Biometrics. 1980, 36: 627-642.

    CAS  PubMed  Google Scholar 

  90. Trikalinos TA, Salanti G, Khoury MJ, Ioannidis JP: Impact of violations and deviations in Hardy-Weinberg equilibrium on postulated gene-disease associations. Am J Epidemiol. 2006, 163 (4): 300-309.

    PubMed  Google Scholar 

  91. Hernandez JL, Weir BS: A disequilibrium coefficient approach to Hardy-Weinberg testing. Biometrics. 1989, 45 (1): 53-70.

    CAS  PubMed  Google Scholar 

  92. Lau J, Ioannidis JP, Schmid CH: Quantitative synthesis in systematic reviews. Ann Intern Med. 1997, 127 (9): 820-826.

    CAS  PubMed  Google Scholar 

  93. Sterne J, Bradburn M, Egger M: Meta-analysis in Stata. Systematic reviews in health care. Edited by: Egger M, Davey Smith G, Altman D. 2001, Boston, MA: Blackwell BMJ Books, 347-369. 2

    Google Scholar 

  94. Higgins JP, Thompson SG, Deeks JJ, Altman DG: Measuring inconsistency in meta-analyses. Bmj. 2003, 327 (7414): 557-560.

    PubMed  PubMed Central  Google Scholar 

  95. Higgins JP, Thompson SG: Quantifying heterogeneity in a meta-analysis. Stat Med. 2002, 21 (11): 1539-1558.

    PubMed  Google Scholar 

  96. Thompson SG: Why sources of heterogeneity in meta-analysis should be investigated. Bmj. 1994, 309 (6965): 1351-1355.

    CAS  PubMed  PubMed Central  Google Scholar 

  97. Thompson SG, Sharp SJ: Explaining heterogeneity in meta-analysis: a comparison of methods. Stat Med. 1999, 18 (20): 2693-2708.

    CAS  PubMed  Google Scholar 

  98. Thompson SG, Smith TC, Sharp SJ: Investigating underlying risk as a source of heterogeneity in meta-analysis. Stat Med. 1997, 16 (23): 2741-2758.

    CAS  PubMed  Google Scholar 

  99. Salanti G, Sanderson S, Higgins JP: Obstacles and opportunities in meta-analysis of genetic association studies. Genet Med. 2005, 7 (1): 13-20.

    PubMed  Google Scholar 

  100. Egger M, Davey Smith G, Schneider M, Minder C: Bias in meta-analysis detected by a simple, graphical test. Bmj. 1997, 315 (7109): 629-634.

    CAS  PubMed  PubMed Central  Google Scholar 

  101. Lau J, Antman EM, Jimenez-Silva J, Kupelnick B, Mosteller F, Chalmers TC: Cumulative meta-analysis of therapeutic trials for myocardial infarction. N Engl J Med. 1992, 327 (4): 248-254.

    CAS  PubMed  Google Scholar 

  102. Stata Corporation. Stata statistical software, release 9.0. 2005, College Station, TX: Stata Corporation

  103. Spiegelhalter D, Thomas A, Best N, et al: WinBUGS user manual. Version 1.4, January 2003. 2003, Cambridge, United Kingdom: MRC Biostatistics Unit, Institute of Public Health, [http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/manual14.pdf]

    Google Scholar 

  104. Dahlgren A, Zethelius B, Jensevik K, Syvanen AC, Berne C: Variants of the TCF7L2 gene are associated with beta cell dysfunction and confer an increased risk of type 2 diabetes mellitus in the ULSAM cohort of Swedish elderly men. Diabetologia. 2007, 50 (9): 1852-1857.

    CAS  PubMed  Google Scholar 

  105. Elbein SC, Chu WS, Das SK, Yao-Borengasser A, Hasstedt SJ, Wang H, Rasouli N, Kern PA: Transcription factor 7-like 2 polymorphisms and type 2 diabetes, glucose homeostasis traits and gene expression in US participants of European and African descent. Diabetologia. 2007, 50 (8): 1621-1630.

    CAS  PubMed  Google Scholar 

  106. Florez JC, Jablonski KA, Bayley N, Pollin TI, de Bakker PI, Shuldiner AR, Knowler WC, Nathan DM, Altshuler D: TCF7L2 polymorphisms and progression to diabetes in the Diabetes Prevention Program. N Engl J Med. 2006, 355 (3): 241-250.

    CAS  PubMed  PubMed Central  Google Scholar 

  107. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447 (7145): 661-678.

  108. Bort R, Martinez-Barbera JP, Beddington RS, Zaret KS: Hex homeobox gene-dependent tissue positioning is required for organogenesis of the ventral pancreas. Development. 2004, 131 (4): 797-806.

    CAS  PubMed  Google Scholar 

  109. Bort R, Signore M, Tremblay K, Martinez Barbera JP, Zaret KS: Hex homeobox gene controls the transition of the endoderm to a pseudostratified, cell emergent epithelium for liver bud development. Dev Biol. 2006, 290 (1): 44-56.

    CAS  PubMed  Google Scholar 

  110. Foley AC, Mercola M: Heart induction by Wnt antagonists depends on the homeodomain transcription factor Hex. Genes Dev. 2005, 19 (3): 387-396.

    CAS  PubMed  PubMed Central  Google Scholar 

  111. Groves CJ, Wiltshire S, Smedley D, Owen KR, Frayling TM, Walker M, Hitman GA, Levy JC, O'Rahilly S, Menzel S, et al: Association and haplotype analysis of the insulin-degrading enzyme (IDE) gene, a strong positional and biological candidate for type 2 diabetes susceptibility. Diabetes. 2003, 52 (5): 1300-1305.

    CAS  PubMed  Google Scholar 

  112. Karamohamed S, Demissie S, Volcjak J, Liu C, Heard-Costa N, Liu J, Shoemaker CM, Panhuysen CI, Meigs JB, Wilson P, et al: Polymorphisms in the insulin-degrading enzyme gene are associated with type 2 diabetes in men from the NHLBI Framingham Heart Study. Diabetes. 2003, 52 (6): 1562-1567.

    CAS  PubMed  Google Scholar 

  113. Florez JC, Wiltshire S, Agapakis CM, Burtt NP, de Bakker PI, Almgren P, Bengtsson Bostrom K, Tuomi T, Gaudet D, Daly MJ, Hirschhorn JN, McCarthy MI, Altshuler D, Groop L: High-density haplotype structure and association testing of the insulin-degrading enzyme (IDE) gene with type 2 diabetes in 4,206 people. Diabetes. 2006, 55 (1): 128-135.

    CAS  PubMed  Google Scholar 

  114. Benetti R, Copetti T, Dell'Orso S, Melloni E, Brancolini C, Monte M, Schneider C: The calpain system is involved in the constitutive regulation of beta-catenin signaling functions. J Biol Chem. 2005, 280 (23): 22070-22080.

    CAS  PubMed  Google Scholar 

  115. Rulifson IC, Karnik SK, Heiser PW, ten Berge D, Chen H, Gu X, Taketo MM, Nusse R, Hebrok M, Kim SK: Wnt signaling regulates pancreatic beta cell proliferation. Proc Natl Acad Sci USA. 2007, 104 (15): 6247-6252.

    CAS  PubMed  PubMed Central  Google Scholar 

  116. Zhan SY: Epidemiology of Diabetes Mellitus. Epidemiology. Edited by: Li LM. 2003, People's Publishing House of Health, Beijing, Chin, 418-436. 5

    Google Scholar 

  117. Cauchi S, Choquet H, Gutierrez-Aguilar R, Capel F, Grau K, Proenca C, Dina C, Duval A, Balkau B, Marre M, Potoczna N, Langin D, Horber F, Sorensen TI, Charpentier G, Meyre D, Froguel P: Effects of TCF7L2 polymorphisms on obesity in European populations. Obesity (Silver Spring). 2008, 16 (2): 476-482.

    CAS  Google Scholar 

  118. Wang J, Kuusisto J, Vanttinen M, Kuulasmaa T, Lindstrom J, Tuomilehto J, Uusitupa M, Laakso M: Variants of transcription factor 7-like 2 (TCF7L2) gene predict conversion to type 2 diabetes in the Finnish Diabetes Prevention Study and are associated with impaired glucose regulation and impaired insulin secretion. Diabetologia. 2007, 50 (6): 1192-1200.

    CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

This work was partially supported by grants from China's National Natural Science Foundation (contact No. 30771219) and the Science Research Foundation of Chinese Ministry of Health (contact No. WKJ2007-3-001).

The authors thank Dr. Huaigong Chen for checking coding of some of the data. They also thank the authors who kindly provided genotype and allele frequency data for the meta-analyses.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ben Zhang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

BZ, YT, YL, and YZ carried out the design of this meta-analysis, analyzed the data and drafted the manuscript. YZ and JYY conducted a searching, data extraction, and statistical analysis on the manuscript. YWZ performed a searching and data extraction on the manuscript. HCL participated in designing and writing for the manuscript. All authors read and approved the final manuscript.

Yu Tong, Ying Lin, Yuan Zhang contributed equally to this work.

Electronic supplementary material

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Authors’ original file for figure 8

Authors’ original file for figure 9

Authors’ original file for figure 10

Authors’ original file for figure 11

Authors’ original file for figure 12

Authors’ original file for figure 13

Authors’ original file for figure 14

Authors’ original file for figure 15

Authors’ original file for figure 16

Authors’ original file for figure 17

Authors’ original file for figure 18

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Tong, Y., Lin, Y., Zhang, Y. et al. Association between TCF7L2gene polymorphisms and susceptibility to Type 2 Diabetes Mellitus: a large Human Genome Epidemiology (HuGE) review and meta-analysis. BMC Med Genet 10, 15 (2009). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2350-10-15

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2350-10-15

Keywords