Genome-wide meta-analysis of genetic susceptible genes for Type 2 Diabetes
© Hale et al; licensee BioMed Central Ltd. 2012
Published: 17 December 2012
Many genetic studies, including single gene studies and Genome-wide association studies (GWAS), aim to identify risk alleles for genetic diseases such as Type II Diabetes (T2D). However, in T2D studies, there is a significant amount of the hereditary risk that cannot be simply explained by individual risk genes. There is a need for developing systems biology approaches to integrate comprehensive genetic information and provide new insight on T2D biology.
We performed comprehensive integrative analysis of Single Nucleotide Polymorphisms (SNP's) individually curated from T2D GWAS results and mapped them to T2D candidate risk genes. Using protein-protein interaction data, we constructed a T2D-specific molecular interaction network consisting of T2D genetic risk genes and their interacting gene partners. We then studied the relationship between these T2D genes and curated gene sets.
We determined that T2D candidate risk genes are concentrated in certain parts of the genome, specifically in chromosome 20. Using the T2D genetic network, we identified highly-interconnected network "hub" genes. By incorporating T2D GWAS results, T2D pathways, and T2D genes' functional category information, we further ranked T2D risk genes, T2D-related pathways, and T2D-related functional categories. We found that highly-interconnected T2D disease network “hub” genes most highly associated to T2D genetic risks to be PI3KR1, ESR1, and ENPP1. The well-characterized TCF7L2, contractor to our expectation, was not among the highest-ranked T2D gene list. Many interacted pathways play a role in T2D genetic risks, which includes insulin signalling pathway, type II diabetes pathway, maturity onset diabetes of the young, adipocytokine signalling pathway, and pathways in cancer. We also observed significant crosstalk among T2D gene subnetworks which include insulin secretion, regulation of insulin secretion, response to peptide hormone stimulus, response to insulin stimulus, peptide secretion, glucose homeostasis, and hormone transport. Overview maps involving T2D genes, gene sets, pathways, and their interactions are all reported.
Large-scale systems biology meta-analyses of GWAS results can improve interpretations of genetic variations and genetic risk factors. T2D genetic risks can be attributable to the summative genetic effects of many genes involved in a broad range of signalling pathways and functional networks. The framework developed for T2D studies may serve as a guide for studying other complex diseases.
Type 2 Diabetes (T2D) is a complex metabolic disease that affects 25.8 million Americans in 2011, according to statistics reported by Centers for Disease Control and Prevention (CDC). T2D occurs when the body develops resistance to insulin due to the malfunction of insulin producing β-cells. The developmental process of T2D involves a complex interplay between genetic and environmental factors. However, it is not clear how the underlying genetic defects give rise to T2D pathogenesis over time. Recent T2D genetic study results, particularly those from genome-wide association studies (GWAS), have yielded insights to the molecular mechanisms and underlying genetic risk factors of T2D . Among the many risk genes identified are: transcription factor 7-like 2 (TCF7L2)[2–4], peroxisome proliferator-activated receptor gamma (PPARG)[5–7], and potassium inwardly-rectifying channel, subfamily J, member 11 (KCNJ11)[5, 6].
These GWAS results were challenging to interpret. Many single nucleotide polymorphisms (SNPs) identified from GWAS tend to show strong sample biases and may not extrapolate from one population to another. In T2D, only approximately 28% of the disease heritability may be explained by identified individual SNPs that showed statistical significance in these samples/population--a problem known as missing heritability . The combined effects of multiple risk SNP's can increase the overall odds ratio of T2D by 1.24 per allele for up to 8.68 among 18 risk alleles in one study  and by 1.265 per allele in another study . The additive effect suggests the presence of molecular system structures that are essential to T2D pathogenesis.
To confirm the presence of molecular systems structures that may better explain missing heritability problems for T2D, we adopted a Systems Biology approach to studying T2D genetic risk gene networks as a whole rather than the risk genes individually. Prior to this study, several reports [10, 11] examined genes implicated T2D differential expressions in affected tissues. In this study, we used T2D-associated SNP information curated from the Type 2 Diabetes Genetic Association Database (T2DGADB), which integrated comprehensively reported SNPs, their odds ratios, population description, and all related metadata from various T2D GWAS performed worldwide . We further annotated individual SNPs collected from T2DGADB with information from the DbSNP database , including information such as nearby genes, Chromosomal location, gene functional class, and base changes. To create a model for T2D genetic risk gene molecular systems structure, we built a gene interaction network seeded by T2D risk genes collected from T2DGADB and expanded with high-confidence protein interaction data collected from the Human Annotated and Predicted Protein Interaction database (HAPPI) . We also ranked risk genes in the network according to these high confidence interactions.
T2D risk SNPs and risk genes data collection and curation
Data from both the ftp site and web pages of T2DGADB were downloaded. On the ftp site, only a gene list and a SNP list without annotation were available for download. Therefore, the complete information from individual web pages of T2DGADB was extracted into a single Excel file manually. Data entries with dbSNP SNP cross-references were kept and entries without dbSNP SNP cross-reference information were removed from this study. Gene annotation information is derived from the VEGA  database. The Excel file was imported into the ORACLE 11 g database for subsequent efficient database querying.
Statistical significance testing of T2D risk SNPs and risk genes
Once we collected information integrated from T2DGADB and dbSNP, we applied standard hyper-geometric tests (using an R software package called phyper) to the data set to determine which chromosomes were over-represented/under-represented. We determined significance on three data sets. First, the distribution of dbSNP of all human SNP's in current build of known origins were compared against that of risk SNP's across all chromosomes. Second, the distribution of genes where the risk SNPs can be mapped to were compared against that of all the genes across all chromosomes. Third, the distribution of protein-coding genes where the risk SNPs can be mapped to were compared against that of all the protein-coding genes across all chromosomes.
Construction and analysis of the T2D risk genes network
In this equation, r p score measures the relative significance of a gene/protein in the subnetwork. p and q are proteins in the subnetwork. k is a constant set at two for our purposes. conf (p, q) is the confidence score of the interaction between the two proteins provided by the HAPPI database and is 0 if p and q do not interact. N (p, q) is 1 if the proteins have an interaction and is 0 otherwise. With all of this we are able to generate a network map of all involved genes.
r mod adjusts r p score of any risk gene with both the count of populations, P count , in which significant risk SNPs were identified, and the average odds ratio (OR avj ) of reported risk SNPs found in these studies. Genes containing only one significant study can still be adjusted using the formula provided here. In constructing the final T2D risk gene network, we modified the original network to exclude studies in which the risk genes were determined to be insignificant (r p score < 2) before we calculated r mod scores for risk genes. Cytoscape software was used to visualize network relationships among genes and gene sets (to be described next).
Construction and analysis of the T2D risk genes network
To further sift the results and explore functional connections, we also mapped genes onto known gene sets. For this purpose, we used DAVID [22, 23] to search for enriched KEGG  pathways. We also used GARNET  to identify enriched Gene Ontology categories and their relationships.
SNPs identified from various T2D GWAS
Based on T2DGADB, we collected 4358 T2D SNP entries that cover 518 PubMed articles reporting T2D Genome-wide association studies (GWAS) results worldwide. Since not all study reported statistics on all SNPs, there are only 3720 SNP entries with P-values, 2715 SNP entries with complete odds ratios, 2406 SNP entries with sample size and minor allele frequency values. All together, there are 1269 SNP entries with the above-mentioned complete set of statistics. After comparing collected information against dbSNP entries manually, we validated 333 SNP gene annotations, re-annotated 11 SNP gene annotations, and flagged 140 additional genes that do not appear to be consistent with dbSNP curated locus information. The "mis-annotations" in T2DGADB were partly due to the presence of two genes, e.g., ABCC8 and KCNJ11, for some reported SNPs, therefore confounding manual curations. In other cases, gene symbols that are similar to each other, e.g., AGER and RAGE, which do not even appear on the same chromosome, also seemed to have been mixed up. The flagging of putative genes, e.g., LOC387761, according to dbSNP was performed, because we wanted to prioritize genes with known gene functions. In the end, we collected 4085 distinct SNP entries that cover 1539 SNPs in 370 different genes. Among these SNP entries, 598 SNPs from 255 different genes passed a P-Value significance cutoff of 0.05.
Observation of chromosomal specificity for T2D risk genes
P-value of Chromosomal Specificity Significant Test.
P-value (all genes)
P-value (coding genes)
Contributing to this result are genes including Hepatocyte Nuclear Factor 4 Alpha (HNF4A) and Protein Tyrosine Phosphotase non-receptor type 1 (PTPN1), both of which have been extensively studied for their roles to T2D genetics. Linkage studies dating back to 1997 showed a modest association with this region of chromosome 20 [26, 27], when whole genome data was not available. Particularly interesting is that the interplay between genetics and environment also may act on this chromosomal region, e.g., the epigenetic effect of diet on the promoter regions of HNF4A .
Strong network centrality in T2D high-risk genes
Top-ranking T2D risk genes ordered by their r p scores in the T2D risk gene protein interaction network.
r p Score
Top-ranking T2D risk genes ordered by their r mod scores in the T2D risk gene protein interaction network.
Rank by r mod Score
r mod Score
r p Score
Original Rank by r p Score
Enriched pathways identified in the T2D risk genes.
Adipocytokine signalling pathway
Type II diabetes mellitus
Insulin signalling pathway
Maturity onset diabetes of the young
PPAR signalling pathway
Calcium signalling pathway
Renal cell carcinoma
Hypertrophic cardiomyopathy (HCM)
Aldosterone-regulated sodium reabsorption
VEGF signalling pathway
Pathways in cancer
Our study suggests PI3KR1 may be the gene with highest T2D genetic risk associated. PIK3R1 is a regulatory subunit of the phosphoinositide-3-kinase, which is a protein known to be involved in insulin actions, cancer signaling, and cytokine signaling. While the coverage of the gene's functional relationship to T2D risks in T2DGADB is very limited (with one article only) , there is increasing evidence, including a recent SNP UTR study  and a mixed methods meta-analysis , that supports our finding.
The second highest-ranking gene in the T2D risk gene network by r p scores is ESR1, the estrogen receptor 1 gene. The gene encodes a transcription factor that responds to estrogen action and cancer, and will also form a heterodimer with ESR2. There are two articles cited in the T2DGADB database [32, 33].
The third highest-ranking gene in the network is ENPP1, Ectonucleotide Pyrophosphatase/Phosphodiesterase 1, a trans-membrane glycoprotein involved in metabolism and has been shown to have an effect of insulin signaling and glucose metabolism . ENPP1 was well studied among 20 GWAS-related publications collected through T2DGADB and 10 of those studies returned positive results in the population examined.
To demonstrate that the network hub genes are indeed functionally associated with T2D risks, we performed a t-test on the distribution of risk SNP per gene between the top ranked 25% of risk genes and the bottom 75% of risk genes. When ranks are given by the original r p scores, the results showed significant difference (P = 0.01) between the two groups for the reported risk SNP per gene raition (3.06 SNP/gene for the top 25% "hub genes" vs. 1.67 SNP/gene for the bottom 75% "non-hub genes"). When ranks are given by the modified r p scores, the results showed even higher significant difference (P = 3.55E-4) between the same two groups (3.74 SNP/gene for the top 25% "hub genes" vs. 1.45 SNP/gene for the bottom 75% "non-hub genes").
Functional heterogeneity and cohesion of T2D high-risk genes
To gain an overview of functional categories represented by the T2D high-risk genes, we mapped these genes onto curated pathways. In Table 4, we list top significantly over-represented pathways identified by the DAVID software. Two highly ranked "novel" pathways identified (excluding known T2D pathways) are the Adipocytokine signaling and the PPAR signaling pathway.
Enriched gene ontology categories identified in the T2D risk genes.
Gene Ontology Term
regulation of glucose transport
regulation of insulin secretion
response to insulin stimulus
response to peptide hormone stimulus
cellular response to hormone stimulus
positive regulation of glucose import
positive regulation of fatty acid metabolic process
cellular response to insulin stimulus
negative regulation of macrophage derived foam cell differentiation
positive regulation of glucose metabolic process
regulation of lipid metabolic process
In this study, we showed our findings of T2D genetic risk gens. Genes from the chromosome 20 collectively accounted for the highest T2D genetic risks of all the chromosomes. However, the individual contribution of these chromosome 20 risk genes is relatively small. HNF4A as the most significant gene on chromosome 20 has a relatively small r p score of 4.72, which is far lower than all r p scores shown for the top-ranking T2D risk genes in Table 2. Nonetheless, when all other information derived from GWAS results are integrated into the r mod score of 118, the significant contribution of HNF4A to T2D risks becomes clear. The "missing inheritability" problem of T2D genetic risks are therefore partially explained with the calculated integration of T2D molecular interaction network information and T2D genotype-phenotype association study results.
GWAS results show population-specific effectiveness in using TCF7L2 for T2D genetic risk profiling.
Number of Studies
Average Odds Ratio
Maximum Odds Ratio
Minimum Odds Ratio
TCF7L2 was not among the highest ranked T2D genes, primarily due to the emphasis on the quality of protein interaction data that we bring in. The HAPPI database reported 238 protein interactions for TCF7L2 but only 11 of those were above the confidence threshold of 0.8. This is in contrast to ENPP1, in which we identified 743 protein interactions from the HAPPI database and 87 of them passed our confidence threshold of 0.8. The relatively low network centrality explained why TCF7L2 is not ranked at the top overall, although it is a population target for many T2D GWAS.
Future analysis that is built upon this work could benefit by integrating additional genomics and functional genomics information, e.g., available miRNA or mRNA expression information, available copy number variations results, and whole genome sequencing data.
Large-scale systems biology meta-analyses of GWAS results can improve interpretations of genetic variations and genetic risk factors. In this work, we determined that T2D candidate risk genes are located in higher concentration in certain parts of the genome, specifically in chromosome 20. Using the T2D genetic network, we identified highly interconnected network "hub" genes. By incorporating T2D GWAS results, T2D pathways, and T2D genes' functional category information, we further ranked T2D risk genes, T2D-related pathways, and T2D-related functional categories. Overview maps involving T2D genes, gene sets, pathways, and their interactions are all reported. Moreover, we demonstrate a computational framework built upon disease-specific data integration, model construction, and data analysis. The framework developed for T2D studies may serve as a guide for studying other complex diseases.
This article has been published as part of BMC Systems Biology Volume 6 Supplement 3, 2012: Proceedings of The International Conference on Intelligent Biology and Medicine (ICIBM) - Systems Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcsystbiol/supplements/6/S3.
- Wheeler E, Barroso I: Genome-wide association studies and type 2 diabetes. Brief Funct Genomics. 2011, 10 (2): 52-60. 10.1093/bfgp/elr008.View ArticlePubMedGoogle Scholar
- Cruz M, Valladares-Salgado A, Garcia-Mena J, Ross K, Edwards M, Angeles-Martinez J, Ortega-Camarillo C, de la Pena JE, Burguete-Garcia AI, Wacher-Rodarte N, et al: Candidate gene association study conditioning on individual ancestry in patients with type 2 diabetes and metabolic syndrome from Mexico City. Diabetes Metab Res Rev. 2010, 26 (4): 261-270. 10.1002/dmrr.1082.View ArticlePubMedGoogle Scholar
- Lin Y, Li P, Cai L, Zhang B, Tang X, Zhang X, Li Y, Xian Y, Yang Y, Wang L, et al: Association study of genetic variants in eight genes/loci with type 2 diabetes in a Han Chinese population. BMC Med Genet. 2010, 11: 97-PubMed CentralView ArticlePubMedGoogle Scholar
- Sladek R, Rocheleau G, Rung J, Dina C, Shen L, Serre D, Boutin P, Vincent D, Belisle A, Hadjadj S, et al: A genome-wide association study identifies novel risk loci for type 2 diabetes. Nature. 2007, 445 (7130): 881-885. 10.1038/nature05616.View ArticlePubMedGoogle Scholar
- Gaulton KJ, Willer CJ, Li Y, Scott LJ, Conneely KN, Jackson AU, Duren WL, Chines PS, Narisu N, Bonnycastle LL, et al: Comprehensive association study of type 2 diabetes and related quantitative traits with 222 candidate genes. Diabetes. 2008, 57 (11): 3136-3144. 10.2337/db07-1731.PubMed CentralView ArticlePubMedGoogle Scholar
- Hu C, Zhang R, Wang C, Wang J, Ma X, Lu J, Qin W, Hou X, Bao Y, Xiang K, et al: PPARG, KCNJ11, CDKAL1, CDKN2A-CDKN2B, IDE-KIF11-HHEX, IGF2BP2 and SLC30A8 are associated with type 2 diabetes in a Chinese population. PLoS One. 2009, 4 (10): e7643-10.1371/journal.pone.0007643.PubMed CentralView ArticlePubMedGoogle Scholar
- Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, Duren WL, Erdos MR, Stringham HM, Chines PS, Jackson AU, et al: A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science. 2007, 316 (5829): 1341-1345. 10.1126/science.1142382.PubMed CentralView ArticlePubMedGoogle Scholar
- So HC, Gui AH, Cherny SS, Sham PC: Evaluating the heritability explained by known susceptibility variants: a survey of ten complex diseases. Genet Epidemiol. 2011, 35 (5): 310-317. 10.1002/gepi.20579.View ArticlePubMedGoogle Scholar
- Cauchi S, Meyre D, Durand E, Proenca C, Marre M, Hadjadj S, Choquet H, De Graeve F, Gaget S, Allegaert F, et al: Post genome-wide association studies of novel genes associated with type 2 diabetes show gene-gene interaction and high predictive value. PLoS One. 2008, 3 (5): e2031-10.1371/journal.pone.0002031.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharma A, Chavali S, Tabassum R, Tandon N, Bharadwaj D: Gene prioritization in Type 2 Diabetes using domain interactions and network analysis. BMC Genomics. 2010, 11: 84-10.1186/1471-2164-11-84.PubMed CentralView ArticlePubMedGoogle Scholar
- Jesmin J, Rashid MS, Jamil H, Hontecillas R, Bassaganya-Riera J: Gene regulatory network reveals oxidative stress as the underlying molecular mechanism of type 2 diabetes and hypertension. BMC Med Genomics. 2010, 3 (1): 45-10.1186/1755-8794-3-45.PubMed CentralView ArticlePubMedGoogle Scholar
- Lim JE, Hong KW, Jin HS, Kim YS, Park HK, Oh B: Type 2 diabetes genetic association database manually curated for the study design and odds ratio. BMC Med Inform Decis Mak. 2010, 10: 76-10.1186/1472-6947-10-76.PubMed CentralView ArticlePubMedGoogle Scholar
- Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K: dbSNP: the NCBI database of genetic variation. Nucleic acids research. 2001, 29 (1): 308-311. 10.1093/nar/29.1.308.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen JY, Mamidipalli S, Huan T: HAPPI: an online database of comprehensive human annotated and predicted protein interactions. BMC genomics. 2009, 10 (Suppl 1): S16-10.1186/1471-2164-10-S1-S16.PubMed CentralView ArticlePubMedGoogle Scholar
- Loveland J: VEGA, the genome browser with a difference. Brief Bioinform. 2005, 6 (2): 189-193. 10.1093/bib/6.2.189.View ArticlePubMedGoogle Scholar
- von Mering C, Jensen LJ, Kuhn M, Chaffron S, Doerks T, Kruger B, Snel B, Bork P: STRING 7--recent developments in the integration and prediction of protein interactions. Nucleic acids research. 2007, 35 (Database): D358-362. 10.1093/nar/gkl825.PubMed CentralView ArticlePubMedGoogle Scholar
- Brown KR, Jurisica I: Online predicted human interaction database. Bioinformatics. 2005, 21 (9): 2076-2082. 10.1093/bioinformatics/bti273.View ArticlePubMedGoogle Scholar
- Alfarano C, Andrade CE, Anthony K, Bahroos N, Bajec M, Bantoft K, Betel D, Bobechko B, Boutilier K, Burgess E, et al: The Biomolecular Interaction Network Database and related tools 2005 update. Nucleic acids research. 2005, 33 (Database): D418-424.PubMed CentralPubMedGoogle Scholar
- Mishra GR, Suresh M, Kumaran K, Kannabiran N, Suresh S, Bala P, Shivakumar K, Anuradha N, Reddy R, Raghavan TM, et al: Human protein reference database--2006 update. Nucleic acids research. 2006, 34 (Database): D411-414.PubMed CentralView ArticlePubMedGoogle Scholar
- Chatr-aryamontri A, Ceol A, Palazzi LM, Nardelli G, Schneider MV, Castagnoli L, Cesareni G: MINT: the Molecular INTeraction database. Nucleic acids research. 2007, 35 (Database): D572-574. 10.1093/nar/gkl950.PubMed CentralView ArticlePubMedGoogle Scholar
- Li J, Zhu X, Chen JY: Building disease-specific drug-protein connectivity maps from molecular interaction networks and PubMed abstracts. PLoS Comput Biol. 2009, 5 (7): e1000450-10.1371/journal.pcbi.1000450.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang da W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic acids research. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic acids research. 2012, 40 (Database): D109-114.PubMed CentralView ArticlePubMedGoogle Scholar
- Rho K, Kim B, Jang Y, Lee S, Bae T, Seo J, Seo C, Lee J, Kang H, Yu U, et al: GARNET--gene set analysis with exploration of annotation relations. BMC Bioinformatics. 2011, 12 (Suppl 1): S25-10.1186/1471-2105-12-S1-S25.PubMed CentralView ArticlePubMedGoogle Scholar
- Ji L, Malecki M, Warram JH, Yang Y, Rich SS, Krolewski AS: New susceptibility locus for NIDDM is localized to human chromosome 20 q. Diabetes. 1997, 46 (5): 876-881. 10.2337/diabetes.46.5.876.View ArticlePubMedGoogle Scholar
- Zouali H, Hani EH, Philippi A, Vionnet N, Beckmann JS, Demenais F, Froguel P: A susceptibility locus for early-onset non-insulin dependent (type 2) diabetes mellitus maps to chromosome 20 q, proximal to the phosphoenolpyruvate carboxykinase gene. Hum Mol Genet. 1997, 6 (9): 1401-1408. 10.1093/hmg/6.9.1401.View ArticlePubMedGoogle Scholar
- Sandovici I, Smith NH, Nitert MD, Ackers-Johnson M, Uribe-Lewis S, Ito Y, Jones RH, Marquez VE, Cairns W, Tadayyon M, et al: Maternal diet and aging alter the epigenetic control of a promoter-enhancer interaction at the Hnf4a gene in rat pancreatic islets. Proc Natl Acad Sci USA. 2011, 108 (13): 5449-5454. 10.1073/pnas.1019007108.PubMed CentralView ArticlePubMedGoogle Scholar
- Barroso I, Luan J, Middelberg RP, Harding AH, Franks PW, Jakes RW, Clayton D, Schafer AJ, O'Rahilly S, Wareham NJ: Candidate gene association study in type 2 diabetes indicates a role for genes involved in beta-cell function as well as insulin action. PLoS Biol. 2003, 1 (1): E20-10.1371/journal.pbio.0000020.PubMed CentralView ArticlePubMedGoogle Scholar
- Malodobra M, Pilecka A, Gworys B, Adamiec R: Single nucleotide polymorphisms within functional regions of genes implicated in insulin action and association with the insulin resistant phenotype. Mol Cell Biochem. 2011, 349 (1-2): 187-193. 10.1007/s11010-010-0673-5.View ArticlePubMedGoogle Scholar
- Rasche A, Al-Hasani H, Herwig R: Meta-analysis approach identifies candidate genes and associated molecular networks for type-2 diabetes mellitus. BMC Genomics. 2008, 9: 310-10.1186/1471-2164-9-310.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang Q, Wang TH, Lu WS, Mu PW, Yang YF, Liang WW, Li CX, Lin GP: Estrogen receptor alpha gene polymorphism associated with type 2 diabetes mellitus and the serum lipid concentration in Chinese women in Guangzhou. Chin Med J (Engl). 2006, 119 (21): 1794-1801.Google Scholar
- Dahlman I, Vaxillaire M, Nilsson M, Lecoeur C, Gu HF, Cavalcanti-Proenca C, Efendic S, Ostenson CG, Brismar K, Charpentier G, et al: Estrogen receptor alpha gene variants associate with type 2 diabetes and fasting plasma glucose. Pharmacogenet Genomics. 2008, 18 (11): 967-975. 10.1097/FPC.0b013e32831101ef.View ArticlePubMedGoogle Scholar
- Di Paola R, Caporarello N, Marucci A, Dimatteo C, Iadicicco C, Del Guerra S, Prudente S, Sudano D, Miele C, Parrino C, et al: ENPP1 affects insulin action and secretion: evidences from in vitro studies. PLoS One. 2011, 6 (5): e19462-10.1371/journal.pone.0019462.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. 10.1007/s00125-007-0656-6.View ArticlePubMedGoogle Scholar
- 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. 10.1007/s00125-006-0588-6.View ArticlePubMedGoogle Scholar
- Ren Q, Han XY, Wang F, Zhang XY, Han LC, Luo YY, Zhou XH, Ji LN: Exon sequencing and association analysis of polymorphisms in TCF7L2 with type 2 diabetes in a Chinese population. Diabetologia. 2008, 51 (7): 1146-1152. 10.1007/s00125-008-1039-3.View ArticlePubMedGoogle Scholar
- 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. 10.2337/db07-0621.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.