Hepatocyte Nuclear Factor-4 P2 Promoter Haplotypes Are Associated With Type 2 Diabetes in the Japanese Population
http://www.100md.com
糖尿病学杂志 2006年第5期
1 Department of Metabolic Diseases, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
2 Department of Clinical Bioinformatics, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
3 CREST, Japan Science and Technology Corporation (JST), Tokyo, Japan
4 National Institute of Health and Nutrition, Tokyo, Japan
5 Institute for Diabetes Care and Research, Asahi Life Foundation, Tokyo, Japan
6 Hiroshima Atomic Bomb Casualty Council Health Management Center, Hiroshima, Japan
7 Department of Endocrinology and Metabolism, International Medical Center of Japan, Tokyo, Japan
8 Department of Human Genetics, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
9 CNRS 8090, Institut de Biologie de Lille, Institute Pasteur and CHU Lille, Lille, France
10 Department of Cardiovascular Medicine, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
HNF, hepatocyte nuclear factor; LD, linkage disequilibrium; MODY, maturity-onset diabetes of the young; SNP, single nucleotide polymorphism
ABSTRACT
Hepatocyte nuclear factor (HNF)-4 is a transcription factor known as a key molecule in the development and functions of the -cells. In a previously performed genome-wide scan of Japanese type 2 diabetic sibpairs, we observed linkage of type 2 diabetes to chromosome 20q12-q13, a region in which the HNF4A gene is located. Recent studies have reported associations between type 2 diabetes and polymorphisms in the P2 promoter region specific to -cells. In this study, we attempted to assess whether the HNF4A gene plays a role in the genetic susceptibility to type 2 diabetes in the Japanese population by analyzing polymorphisms and haplotypes of the HNF4A gene. Linkage disequilibrium across the P2 promoter region was preserved in the Japanese population, consistent with previous reports. Although none of the individual polymorphisms examined showed any significant association with type 2 diabetes, we found very strong evidence of the association between type 2 diabetes and the haplotype consisting of two polymorphisms in the P2 promoter region of the HNF4A gene (P = 3.82 x 10eC4). In contrast, there was no association between type 2 diabetes and haplotypes consisting of polymorphisms not located in the P2 promoter region, suggesting that the type 2 diabetes susceptibility loci are localized in the P2 promoter region of the HNF4A gene. The association was replicated using two additional cohorts (P = 1.51 x 10eC4 and 0.019, respectively). The results of the present analysis revealed that the HNF4A gene might be a type 2 diabetes susceptibility gene common to different ethnic groups. The study also suggested the possible existence of an as-yet-unidentified but functional polymorphism in the P2 promoter region of the HNF4A gene that directly influences susceptibility to type 2 diabetes.
We previously conducted a 10-cM genome-wide scan for regions linked to type 2 diabetes in 224 affected Japanese sibpairs and found one suggestive linked region and seven potentially linked regions, including 20q12-q13 (the logarithm of odds [LOD] in a multipoint analysis was 2.32 [P = 0.00102] at D20S119 in a population subset whose maximum BMI was <30 kg/m2), near the gene for hepatocyte nuclear factor (HNF)-4 (1). Evidence for a type 2 diabetes locus in chromosome 20q12-q13 (OMIM no. 603694) has also been reported from studies in several Caucasian and Chinese populations (2eC5), suggesting that this region may harbor a susceptibility gene for type 2 diabetes common to these ethnic groups. The transcription factor HNF-4 (HNF4A) gene, the gene for maturity-onset diabetes of the young (MODY) type 1 (6), a dominantly inherited, early-onset type 2 diabetes characterized by defective glucose-dependent insulin secretion (7), is located in this region. The HNF4A has a complex expression pattern, in part due to alternative splicing, and is expressed in many tissues, including the liver and pancreas. Three of the isoforms are transcribed by an alternative P2 promoter, located about 46 kb upstream of the P1 promoter and the coding exons. Transcripts from both the P1 and P2 promoters have been detected in pancreatic -cells, but the P2 promoter is suggested to be the major transcription start site in these cells (8eC10). Mutations of HNF4A have been identified in families of MODY type 1 diabetic families, in both the coding and the regulatory regions of the gene, including the P2 promoter region (11). In the -cells, HNF4A regulates the expression of genes involved in glucose metabolism and insulin secretion (12,13). Therefore, HNF4A is currently one of the most attractive candidates as the type 2 diabetes susceptibility gene. Indeed, significant associations between single nucleotide polymorphisms (SNPs) in the P2 promoter region of HNF4A and type 2 diabetes have been shown in three Caucasian populations (14eC16). However, the precise SNPs directly influencing the susceptibility to type 2 diabetes remain to be identified.
In this study, we conducted a haplotype analysis to investigate whether SNPs in the HNF4A gene influence the susceptibility to type 2 diabetes in the Japanese population.
RESEARCH DESIGN AND METHODS
The inclusion criteria used for diabetic and nondiabetic subjects in this study have been previously described (17). Diabetes was diagnosed according to the criteria of the World Health Organization (WHO) (18). All subjects enrolled in this study were of full Japanese ancestry. SNPs were genotyped in 192 nondiabetic subjects (mean age 69.1 ± 6.93 years and mean BMI 23.9 ± 2.77 kg/m2) and type 2 diabetic subjects (mean age 63.1 ± 4.16 years and mean BMI 24.6 ± 2.78 kg/m2) (first case-control cohort). We confirmed the results derived from our first case-control cohort in a second cohort (384 type 2 diabetic subjects with a mean age of 62.5 ± 9.78 years and a mean BMI of 24.4 ± 3.92 kg/m2 and 384 nondiabetic subjects with a mean age of 68.4 ± 9.80 years and a mean BMI of 23.4 ± 5.87 kg/m2) and a third cohort (192 type 2 diabetic subjects with a mean age of 63.7 ± 8.31 years and mean BMI of 23.8 ± 1.39 kg/m2 and 192 nondiabetic subjects with a mean age of 69.1 ± 8.32 years and a mean BMI of 22.9 ± 2.7 kg/m2). The three cohorts of subjects (first, second, and third cohorts) were recruited during different times (in the order of first, second, and third cohorts). The type 2 diabetic subjects of the first and second cohorts were recruited from the outpatient clinic of University of Tokyo Hospital; they were all residents of the Tokyo metropolitan area. The type 2 diabetic subjects of the third cohort and the nondiabetic subjects of all the three cohorts were recruited from the Hiroshima Atomic Bomb Casualty Council Health Management Center; they were all residents of Hiroshima, a city 500 miles west of Tokyo (note that all of the nondiabetic subjects were recruited from Hiroshima). Informed consent was obtained from all of the subjects enrolled in the study, and the study was approved by the Ethics Committee of the University of Tokyo.
Screening and selection of SNPs in HNF4A.
To establish an SNP map encompassing HNF4A, SNPs were ascertained by both direct sequencing and a search of databases. We conducted a screening for SNPs by direct sequencing at all the exons of HNF4A and their exon-intron boundaries and along a 1,500-bp region 5' upstream of exon D, which has been reported to be a -celleCspecific promoter region. The conditions and sequences of the primers used in the PCR shall be made available upon request. The SNPs were identified based on the sequences reported in GenBank that contain the HNF4A gene (accession no. NT_086910). From the public database, rs6065723, rs1884612, rs4812822, rs4810424, rs1884613, rs1884614, rs2144908, rs2273618, and rs3818247 were selected and validated in 30 type 2 diabetic subjects. The SNPs were genotyped in type 2 diabetic subjects and nondiabetic subjects by direct sequencing. PCR was performed under standard conditions, and the sequencing reactions were performed using the BigDye terminator kit (Applied Biosystems, Foster City, CA) and resolved using an ABI 3700 automated DNA sequencer (Applied Biosystems). The results were integrated using a Sequencher (Gene Codes Corporation, Ann Arbor, MI) and individual SNPs were manually genotyped. Ambiguous base callings were eliminated from further analysis. There was no genotyping error based on blind replicates for two SNPs in 192 samples.
Statistical analysis.
The characteristics of the populations were described as means ± SD. The proportions of specific genotypes or alleles in subjects with and without type 2 diabetes were compared using the 2 test. The differences among subjects with different SNP genotypes were statistically tested using ANOVA. The statistical analyses, except for the haplotype estimation, were performed using JMP for Windows, Version 4.00 (SAS Institute, Cary, NC).
Haplotype analysis.
The frequencies of each haplotype were estimated, and the differences in the haplotype frequencies between nondiabetic and diabetic subjects were assessed using a software based on the expectation-maximization (EM) algorithm, SNPAlyze (Dynacom, Tokyo, Japan), and PHASE, Version 2 (http://www.stat.washington.edu/stephens/software.html). The differences in the haplotype frequencies were then analyzed using the 2 and permutation tests. The linkage disequilibrium (LD) structure in the HNF4A region was plotted using GOLD software.
RESULTS AND DISCUSSION
Genomic DNAs from 192 unrelated type 2 diabetic and nondiabetic Japanese individuals were genotyped at an average of one SNP every 4.1 kb across an 78-kb region harboring the HNF4A gene and its alternative upstream promoter, P2. All the genotypic distributions of the SNPs that were identified were in Hardy-Weinberg equilibrium in both the nondiabetic and type 2 diabetic subjects (P > 0.05). Among those identified, SNPs with a minor allele frequency >0.10 were investigated for LD in HNF4A and for their association with type 2 diabetes. As shown in Fig. 1A, SNPs were used to determine the pattern of LD in the chromosomal region of HNF4A. The LD plot illustrates that the chromosomal region beginning with rs6065723, eC27 kb upstream of the P2 promoter, and rs6073418, 17kb downstream of the P2 promoter, represents one block of strong LD, which is consistent with previous reports (15,16,20). No difference in the distributions of either the genotypes or the SNPs in HNF4A were observed between the nondiabetic and type 2 diabetic subjects (Table 1). We then tested the association of haplotypes with susceptibility to type 2 diabetes and found very strong evidence of an association between type 2 diabetes and a certain haplotype in the promoter region of the HNF4A gene. This highly significant association was seen with the haplotype consisting of two SNPs, namely rs1884614 and rs2144908, in the P2 promoter region of the HFN4A gene (overall P value = 3.82 x 10eC4) (Table 2). In contrast, there was no significant association between type 2 diabetes and haplotypes consisting of polymorphisms not located in the P2 promoter region. Indeed, the extent of association between haplotypes and type 2 diabetes decreased with increasing distance from the P2 promoter region (Fig. 1B). This result supports the idea that a type 2 diabetes susceptibility polymorphism might be located in the P2 promoter region but not in any other region of the HNF4A gene. Martin et al. (19) reported that haplotype analysis allowed more accurate mapping of the Alzheimer’s disease susceptibility locus to the APOE gene than analysis of individual SNPs, similar to our findings. The association of type 2 diabetes with the rs1884614-rs2144908 haplotype was replicated in two other case-control cohorts. The second cohort (overall P value = 1.51 x 10eC4) consisted of 384 diabetic and 384 nondiabetic subjects. The third cohort (overall P value = 0.019) (Table 2) included both type 2 diabetic and nondiabetic subjects recruited from the same area of Japan to exclude the possibility of false haplotype associations as a result of population stratification between samples enrolled from different areas of Japan (Tokyo and Hiroshima). The T-G rs1884614-rs2144908 haplotype was found consistently more frequently in the type 2 diabetic subjects than in the nondiabetic subjects of the first (0.052 vs. 0.024, haplotype-specific P value = 0.04, odds ratio [OR] 2.23 [95% CI 1.06eC5.35]), second (0.042 vs. 0.012, haplotype-specific P value = 3.6 x 10eC4, OR 3.50 [95% CI 1.67eC7.11]), and third cohorts (0.039 vs. 0.007, haplotype-specific P value = 0.003, OR 5.83 [95% CI 1.69eC22.4]) (Table 2).
We found no association between individual SNPs in the P2 promoter region and type 2 diabetes, which differs from the results reported by Silander et al. (15) in a Finnish population and by Love-Gregory et al. (16) in Ashkenazi Jews. However, our result was consistent with that reported in two recent articles published by Bagwell et al. (20) and Winckler et al. (21), who reported no evidence of any association of single SNPs with type 2 diabetes in Scandinavians and American Caucasians. These discrepant results among the studies may be attributable to the different genetic backgrounds of the subjects. In contrast, we found significant and consistent association between certain haplotypes and type 2 diabetes, which is consistent with the report of Bagwell et al. in which haplotypes consisting of SNPs in the P2 promoter region, but not any individual SNPs, were significantly associated with type 2 diabetes. The haplotype found to be associated with type 2 diabetes in the present study is relatively uncommon. However, it is noteworthy that the association between this haplotype and type 2 diabetes was observed consistently in all the three cohorts, suggesting the existence of an as-yet-unidentified polymorphism that might be relatively uncommon but associated with susceptibility to type 2 diabetes.
The results of the present study revealed that the P2 region of the HNF4A gene might be a susceptibility region for type 2 diabetes common to different ethnic groups, including the Japanese. Given that specific haplotypes, but not individual SNPs, in the P2 promoter region were significantly associated with type 2 diabetes, there might exist an as-yet-unidentified but functional polymorphism in the P2 promoter region of the HNF4A gene directly influencing the susceptibility to type 2 diabetes. It is also possible that SNPs constituting the haplotype in the P2 promoter region of the HNF4A gene associated with type 2 diabetes may coordinately affect the susceptibility to type 2 diabetes.
ACKNOWLEDGMENTS
This study was supported by the Program for Promotion of Fundamental Studies in Health Sciences of the National Institute of Biomedical Innovation (to T.K.); Grant-in-Aid for Scientific Research on Priority Areas "Applied Genomics" from the Ministry of Education, Culture, Sports, Science and Technology of Japan; and Grant-in-Aid from the 21st Century Program (to R.N.).
We thank Naoko Miyama and Yuko Okada for technical assistance.
FOOTNOTES
K.H. and M.H. contributed equally to this study.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
Mori Y, Otabe S, Dina C, Yasuda K, Populaire C, Lecoeur C, Vatin V, Durand E, Hara K, Okada T, Tobe K, Boutin P, Kadowaki T, Froguel P: Genome-wide search for type 2 diabetes in Japanese affected sib-pairs confirms susceptibility genes on 3q, 15q, and 20q and identifies two new candidate loci on 7p and 11p. Diabetes 51:1247eC1255, 2002
Ghosh S, Watanabe RM, Hauser ER, Valle T, Magnuson VL, Erdos MR, Langefeld CD, Balow J Jr, Ally DS, Kohtamaki K, Chines P, Birznieks G, Kaleta HS, Musick A, Te C, Tannenbaum J, Eldridge W, Shapiro S, Martin C, Witt A, So A, Chang J, Shurtleff B, Porter R, Boehnke M.: Type 2 diabetes: evidence for linkage on chromosome 20 in 716 Finnish affected sib pairs. Proc Natl Acad Sci U S A 96:2198eC2203, 1999
Bowden DW, Sale M, Howard TD, Qadri A, Spray BJ, Rothschild CB, Akots G, Rich SS, Freedman BI: Linkage of genetic markers on human chromosomes 20 and 12 to NIDDM in Caucasian sib pairs with a history of diabetic nephropathy. Diabetes 46:882eC886, 1997
Ji L, Malecki M, Warram JH, Yang Y, Rich SS, Krolewski AS: New susceptibility locus for NIDDM is localized to human chromosome 20q. Diabetes 46:876eC881, 1997
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 20q, proximal to the phosphoenolpyruvate carboxykinase gene. Hum Mol Genet 6:1401eC1408, 1997
Fajans SS, Bell GI, Polonsky KS: Molecular mechanisms and clinical pathophysiology of maturity-onset diabetes of the young. N Engl J Med 345:971eC980, 2001
Winter WE: Newly defined genetic diabetes syndromes: maturity onset diabetes of the young. Rev Endocr Metab Disord 4:43eC51, 2003
Nakhei H, Lingott A, Lemm I, Ryffel GU: An alternative splice variant of the tissue specific transcription factor HNF4alpha predominates in undifferentiated murine cell types. Nucleic Acid Res 26:497eC504, 1998
Boj SF, Parrizas M, Maestro MA, Ferrer J: A transcription factor regulatory circuit in differentiated pancreatic cells. Proc Natl Acad Sci U S A 98:14481eC14486, 2001
Eeckhoute J, Moerman E, Bouckenooghe T, Lukoviak B, Pattou F, Formstecher P, Kerr-Conte J, Vandewalle B, Laine B: Hepatocyte nuclear factor 4 alpha isoforms originated from the P1 promoter are expressed in human pancreatic beta-cells and exhibit stronger transcriptional potentials than P2 promoter-driven isoforms. Endocrinology 144:1686eC1694, 2003
Hansen SK, Parrizas M, Jensen ML, Pruhova S, Ek J, Boj SF, Johansen A, Maestro MA, Rivera F, Eiberg H, Andel M, Lebl J, Pedersen O, Ferrer J, Hansen T: Genetic evidence that HNF-1alpha-dependent transcriptional control of HNF-4alpha is essential for human pancreatic beta cell function. J Clin Invest 110:827eC833, 2002
Stoffel M, Duncan SA: The maturity-onset diabetes of the young (MODY1) transcription factor HNF4alpha regulates expression of genes required for glucose transport and metabolism. Proc Natl Acad Sci U S A 94:13209eC13214, 1997
Wang H, Maechler P, Antinozzi PA, Hagenfeldt KA, Wollheim CB: Hepatocyte nuclear factor 4alpha regulates the expression of pancreatic beta-cell genes implicated in glucose metabolism and nutrient-induced insulin secretion. J Biol Chem 275:35953eC35959, 2000
Damcott CM, Hoppman N, Ott SH, Reinhart LJ, Wang J, Pollin TI, O’connell JR, Mitchell BD, Shuldiner AR: Polymorphisms in both promoters of hepatocyte nuclear factor 4a are associated with type 2 diabetes in the Amish. Diabetes 53:3337eC3341, 2004
Silander K, Mohlke KL, Scott LJ, Peck EC, Hollstein P, Skol AD, Jackson AU, Deloukas P, Hunt S, Stavrides G, Chines PS, Erdos MR, Narisu N, Conneely KN, Li C, Fingerlin TE, Dhanjal SK, Valle TT, Bergman RN, Tuomilehto J, Watanabe RM, Boehnke M, Collins FS: Genetic variation near the hepatocyte nuclear factor-4 gene predicts susceptibility to type 2 diabetes. Diabetes 53:1141eC1149, 2004
Love-Gregory LD, Wasson J, Ma J, Jin CH, Glaser B, Suarez BK, Permutt MA: A common polymorphism in the upstream promoter region of the hepatocyte nuclear factor-4 gene on chromosome 20q is associated with type 2 diabetes and appears to contribute to the evidence for linkage in an Ashkenazi Jewish population. Diabetes 53:1134eC1140, 2004
Hara K, Boutin P, Mori Y, Tobe K, Dina C, Yasuda K, Yamauchi T, Otabe S, Okada T, Eto K, Kadowaki H, Hagura R, Akanuma Y, Yazaki Y, Nagai R, Taniyama M, Matsubara K, Yoda M, Nakano Y, Tomita M, Kimura S, Ito C, Froguel P, Kadowaki T: Genetic variation in the gene encoding adiponectin is associated with increased risk of type 2 diabetes in the Japanese population. Diabetes 51:536eC550, 2002
Alberti KG, Zimmet PZ: Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet Med 15:539eC553, 1998
Martin ER, Lai EH, Gilbert JR, Rogala AR, Afshari AJ, Riley J, Finch KL, Stevens JF, Livak KJ, Slotterbeck BD, Slifer SH, Warren LL, Conneally PM, Schmechel DE, Purvis I, Pericak-Vance MA, Roses AD, Vance JM: SNPing away at complex diseases: analysis of single-nucleotide polymorphisms around APOE in Alzheimer disease. Am J Hum Genet 67:383eC394, 2000
Bagwell AM, Bento JL, Mychaleckyj JC, Freedman BI, Langefeld CD, Bowden DW: Genetic analysis of HNF4A polymorphisms in Caucasian-American type 2 diabetes. Diabetes 54:1185eC1190, 2005
Winckler W, Graham RR, de Bakker PI, Sun M, Almgren P, Tuomi T, Gaudet D, Hudson TJ, Ardlie KG, Daly MJ, Hirschhorn JN, Groop L, Altshuler D: Association testing of variants in the hepatocyte nuclear factor 4 gene with risk of type 2 diabetes in 7,883 people. Diabetes 54:886eC892, 2005(Kazuo Hara, Momoko Horiko)
2 Department of Clinical Bioinformatics, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
3 CREST, Japan Science and Technology Corporation (JST), Tokyo, Japan
4 National Institute of Health and Nutrition, Tokyo, Japan
5 Institute for Diabetes Care and Research, Asahi Life Foundation, Tokyo, Japan
6 Hiroshima Atomic Bomb Casualty Council Health Management Center, Hiroshima, Japan
7 Department of Endocrinology and Metabolism, International Medical Center of Japan, Tokyo, Japan
8 Department of Human Genetics, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
9 CNRS 8090, Institut de Biologie de Lille, Institute Pasteur and CHU Lille, Lille, France
10 Department of Cardiovascular Medicine, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
HNF, hepatocyte nuclear factor; LD, linkage disequilibrium; MODY, maturity-onset diabetes of the young; SNP, single nucleotide polymorphism
ABSTRACT
Hepatocyte nuclear factor (HNF)-4 is a transcription factor known as a key molecule in the development and functions of the -cells. In a previously performed genome-wide scan of Japanese type 2 diabetic sibpairs, we observed linkage of type 2 diabetes to chromosome 20q12-q13, a region in which the HNF4A gene is located. Recent studies have reported associations between type 2 diabetes and polymorphisms in the P2 promoter region specific to -cells. In this study, we attempted to assess whether the HNF4A gene plays a role in the genetic susceptibility to type 2 diabetes in the Japanese population by analyzing polymorphisms and haplotypes of the HNF4A gene. Linkage disequilibrium across the P2 promoter region was preserved in the Japanese population, consistent with previous reports. Although none of the individual polymorphisms examined showed any significant association with type 2 diabetes, we found very strong evidence of the association between type 2 diabetes and the haplotype consisting of two polymorphisms in the P2 promoter region of the HNF4A gene (P = 3.82 x 10eC4). In contrast, there was no association between type 2 diabetes and haplotypes consisting of polymorphisms not located in the P2 promoter region, suggesting that the type 2 diabetes susceptibility loci are localized in the P2 promoter region of the HNF4A gene. The association was replicated using two additional cohorts (P = 1.51 x 10eC4 and 0.019, respectively). The results of the present analysis revealed that the HNF4A gene might be a type 2 diabetes susceptibility gene common to different ethnic groups. The study also suggested the possible existence of an as-yet-unidentified but functional polymorphism in the P2 promoter region of the HNF4A gene that directly influences susceptibility to type 2 diabetes.
We previously conducted a 10-cM genome-wide scan for regions linked to type 2 diabetes in 224 affected Japanese sibpairs and found one suggestive linked region and seven potentially linked regions, including 20q12-q13 (the logarithm of odds [LOD] in a multipoint analysis was 2.32 [P = 0.00102] at D20S119 in a population subset whose maximum BMI was <30 kg/m2), near the gene for hepatocyte nuclear factor (HNF)-4 (1). Evidence for a type 2 diabetes locus in chromosome 20q12-q13 (OMIM no. 603694) has also been reported from studies in several Caucasian and Chinese populations (2eC5), suggesting that this region may harbor a susceptibility gene for type 2 diabetes common to these ethnic groups. The transcription factor HNF-4 (HNF4A) gene, the gene for maturity-onset diabetes of the young (MODY) type 1 (6), a dominantly inherited, early-onset type 2 diabetes characterized by defective glucose-dependent insulin secretion (7), is located in this region. The HNF4A has a complex expression pattern, in part due to alternative splicing, and is expressed in many tissues, including the liver and pancreas. Three of the isoforms are transcribed by an alternative P2 promoter, located about 46 kb upstream of the P1 promoter and the coding exons. Transcripts from both the P1 and P2 promoters have been detected in pancreatic -cells, but the P2 promoter is suggested to be the major transcription start site in these cells (8eC10). Mutations of HNF4A have been identified in families of MODY type 1 diabetic families, in both the coding and the regulatory regions of the gene, including the P2 promoter region (11). In the -cells, HNF4A regulates the expression of genes involved in glucose metabolism and insulin secretion (12,13). Therefore, HNF4A is currently one of the most attractive candidates as the type 2 diabetes susceptibility gene. Indeed, significant associations between single nucleotide polymorphisms (SNPs) in the P2 promoter region of HNF4A and type 2 diabetes have been shown in three Caucasian populations (14eC16). However, the precise SNPs directly influencing the susceptibility to type 2 diabetes remain to be identified.
In this study, we conducted a haplotype analysis to investigate whether SNPs in the HNF4A gene influence the susceptibility to type 2 diabetes in the Japanese population.
RESEARCH DESIGN AND METHODS
The inclusion criteria used for diabetic and nondiabetic subjects in this study have been previously described (17). Diabetes was diagnosed according to the criteria of the World Health Organization (WHO) (18). All subjects enrolled in this study were of full Japanese ancestry. SNPs were genotyped in 192 nondiabetic subjects (mean age 69.1 ± 6.93 years and mean BMI 23.9 ± 2.77 kg/m2) and type 2 diabetic subjects (mean age 63.1 ± 4.16 years and mean BMI 24.6 ± 2.78 kg/m2) (first case-control cohort). We confirmed the results derived from our first case-control cohort in a second cohort (384 type 2 diabetic subjects with a mean age of 62.5 ± 9.78 years and a mean BMI of 24.4 ± 3.92 kg/m2 and 384 nondiabetic subjects with a mean age of 68.4 ± 9.80 years and a mean BMI of 23.4 ± 5.87 kg/m2) and a third cohort (192 type 2 diabetic subjects with a mean age of 63.7 ± 8.31 years and mean BMI of 23.8 ± 1.39 kg/m2 and 192 nondiabetic subjects with a mean age of 69.1 ± 8.32 years and a mean BMI of 22.9 ± 2.7 kg/m2). The three cohorts of subjects (first, second, and third cohorts) were recruited during different times (in the order of first, second, and third cohorts). The type 2 diabetic subjects of the first and second cohorts were recruited from the outpatient clinic of University of Tokyo Hospital; they were all residents of the Tokyo metropolitan area. The type 2 diabetic subjects of the third cohort and the nondiabetic subjects of all the three cohorts were recruited from the Hiroshima Atomic Bomb Casualty Council Health Management Center; they were all residents of Hiroshima, a city 500 miles west of Tokyo (note that all of the nondiabetic subjects were recruited from Hiroshima). Informed consent was obtained from all of the subjects enrolled in the study, and the study was approved by the Ethics Committee of the University of Tokyo.
Screening and selection of SNPs in HNF4A.
To establish an SNP map encompassing HNF4A, SNPs were ascertained by both direct sequencing and a search of databases. We conducted a screening for SNPs by direct sequencing at all the exons of HNF4A and their exon-intron boundaries and along a 1,500-bp region 5' upstream of exon D, which has been reported to be a -celleCspecific promoter region. The conditions and sequences of the primers used in the PCR shall be made available upon request. The SNPs were identified based on the sequences reported in GenBank that contain the HNF4A gene (accession no. NT_086910). From the public database, rs6065723, rs1884612, rs4812822, rs4810424, rs1884613, rs1884614, rs2144908, rs2273618, and rs3818247 were selected and validated in 30 type 2 diabetic subjects. The SNPs were genotyped in type 2 diabetic subjects and nondiabetic subjects by direct sequencing. PCR was performed under standard conditions, and the sequencing reactions were performed using the BigDye terminator kit (Applied Biosystems, Foster City, CA) and resolved using an ABI 3700 automated DNA sequencer (Applied Biosystems). The results were integrated using a Sequencher (Gene Codes Corporation, Ann Arbor, MI) and individual SNPs were manually genotyped. Ambiguous base callings were eliminated from further analysis. There was no genotyping error based on blind replicates for two SNPs in 192 samples.
Statistical analysis.
The characteristics of the populations were described as means ± SD. The proportions of specific genotypes or alleles in subjects with and without type 2 diabetes were compared using the 2 test. The differences among subjects with different SNP genotypes were statistically tested using ANOVA. The statistical analyses, except for the haplotype estimation, were performed using JMP for Windows, Version 4.00 (SAS Institute, Cary, NC).
Haplotype analysis.
The frequencies of each haplotype were estimated, and the differences in the haplotype frequencies between nondiabetic and diabetic subjects were assessed using a software based on the expectation-maximization (EM) algorithm, SNPAlyze (Dynacom, Tokyo, Japan), and PHASE, Version 2 (http://www.stat.washington.edu/stephens/software.html). The differences in the haplotype frequencies were then analyzed using the 2 and permutation tests. The linkage disequilibrium (LD) structure in the HNF4A region was plotted using GOLD software.
RESULTS AND DISCUSSION
Genomic DNAs from 192 unrelated type 2 diabetic and nondiabetic Japanese individuals were genotyped at an average of one SNP every 4.1 kb across an 78-kb region harboring the HNF4A gene and its alternative upstream promoter, P2. All the genotypic distributions of the SNPs that were identified were in Hardy-Weinberg equilibrium in both the nondiabetic and type 2 diabetic subjects (P > 0.05). Among those identified, SNPs with a minor allele frequency >0.10 were investigated for LD in HNF4A and for their association with type 2 diabetes. As shown in Fig. 1A, SNPs were used to determine the pattern of LD in the chromosomal region of HNF4A. The LD plot illustrates that the chromosomal region beginning with rs6065723, eC27 kb upstream of the P2 promoter, and rs6073418, 17kb downstream of the P2 promoter, represents one block of strong LD, which is consistent with previous reports (15,16,20). No difference in the distributions of either the genotypes or the SNPs in HNF4A were observed between the nondiabetic and type 2 diabetic subjects (Table 1). We then tested the association of haplotypes with susceptibility to type 2 diabetes and found very strong evidence of an association between type 2 diabetes and a certain haplotype in the promoter region of the HNF4A gene. This highly significant association was seen with the haplotype consisting of two SNPs, namely rs1884614 and rs2144908, in the P2 promoter region of the HFN4A gene (overall P value = 3.82 x 10eC4) (Table 2). In contrast, there was no significant association between type 2 diabetes and haplotypes consisting of polymorphisms not located in the P2 promoter region. Indeed, the extent of association between haplotypes and type 2 diabetes decreased with increasing distance from the P2 promoter region (Fig. 1B). This result supports the idea that a type 2 diabetes susceptibility polymorphism might be located in the P2 promoter region but not in any other region of the HNF4A gene. Martin et al. (19) reported that haplotype analysis allowed more accurate mapping of the Alzheimer’s disease susceptibility locus to the APOE gene than analysis of individual SNPs, similar to our findings. The association of type 2 diabetes with the rs1884614-rs2144908 haplotype was replicated in two other case-control cohorts. The second cohort (overall P value = 1.51 x 10eC4) consisted of 384 diabetic and 384 nondiabetic subjects. The third cohort (overall P value = 0.019) (Table 2) included both type 2 diabetic and nondiabetic subjects recruited from the same area of Japan to exclude the possibility of false haplotype associations as a result of population stratification between samples enrolled from different areas of Japan (Tokyo and Hiroshima). The T-G rs1884614-rs2144908 haplotype was found consistently more frequently in the type 2 diabetic subjects than in the nondiabetic subjects of the first (0.052 vs. 0.024, haplotype-specific P value = 0.04, odds ratio [OR] 2.23 [95% CI 1.06eC5.35]), second (0.042 vs. 0.012, haplotype-specific P value = 3.6 x 10eC4, OR 3.50 [95% CI 1.67eC7.11]), and third cohorts (0.039 vs. 0.007, haplotype-specific P value = 0.003, OR 5.83 [95% CI 1.69eC22.4]) (Table 2).
We found no association between individual SNPs in the P2 promoter region and type 2 diabetes, which differs from the results reported by Silander et al. (15) in a Finnish population and by Love-Gregory et al. (16) in Ashkenazi Jews. However, our result was consistent with that reported in two recent articles published by Bagwell et al. (20) and Winckler et al. (21), who reported no evidence of any association of single SNPs with type 2 diabetes in Scandinavians and American Caucasians. These discrepant results among the studies may be attributable to the different genetic backgrounds of the subjects. In contrast, we found significant and consistent association between certain haplotypes and type 2 diabetes, which is consistent with the report of Bagwell et al. in which haplotypes consisting of SNPs in the P2 promoter region, but not any individual SNPs, were significantly associated with type 2 diabetes. The haplotype found to be associated with type 2 diabetes in the present study is relatively uncommon. However, it is noteworthy that the association between this haplotype and type 2 diabetes was observed consistently in all the three cohorts, suggesting the existence of an as-yet-unidentified polymorphism that might be relatively uncommon but associated with susceptibility to type 2 diabetes.
The results of the present study revealed that the P2 region of the HNF4A gene might be a susceptibility region for type 2 diabetes common to different ethnic groups, including the Japanese. Given that specific haplotypes, but not individual SNPs, in the P2 promoter region were significantly associated with type 2 diabetes, there might exist an as-yet-unidentified but functional polymorphism in the P2 promoter region of the HNF4A gene directly influencing the susceptibility to type 2 diabetes. It is also possible that SNPs constituting the haplotype in the P2 promoter region of the HNF4A gene associated with type 2 diabetes may coordinately affect the susceptibility to type 2 diabetes.
ACKNOWLEDGMENTS
This study was supported by the Program for Promotion of Fundamental Studies in Health Sciences of the National Institute of Biomedical Innovation (to T.K.); Grant-in-Aid for Scientific Research on Priority Areas "Applied Genomics" from the Ministry of Education, Culture, Sports, Science and Technology of Japan; and Grant-in-Aid from the 21st Century Program (to R.N.).
We thank Naoko Miyama and Yuko Okada for technical assistance.
FOOTNOTES
K.H. and M.H. contributed equally to this study.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
REFERENCES
Mori Y, Otabe S, Dina C, Yasuda K, Populaire C, Lecoeur C, Vatin V, Durand E, Hara K, Okada T, Tobe K, Boutin P, Kadowaki T, Froguel P: Genome-wide search for type 2 diabetes in Japanese affected sib-pairs confirms susceptibility genes on 3q, 15q, and 20q and identifies two new candidate loci on 7p and 11p. Diabetes 51:1247eC1255, 2002
Ghosh S, Watanabe RM, Hauser ER, Valle T, Magnuson VL, Erdos MR, Langefeld CD, Balow J Jr, Ally DS, Kohtamaki K, Chines P, Birznieks G, Kaleta HS, Musick A, Te C, Tannenbaum J, Eldridge W, Shapiro S, Martin C, Witt A, So A, Chang J, Shurtleff B, Porter R, Boehnke M.: Type 2 diabetes: evidence for linkage on chromosome 20 in 716 Finnish affected sib pairs. Proc Natl Acad Sci U S A 96:2198eC2203, 1999
Bowden DW, Sale M, Howard TD, Qadri A, Spray BJ, Rothschild CB, Akots G, Rich SS, Freedman BI: Linkage of genetic markers on human chromosomes 20 and 12 to NIDDM in Caucasian sib pairs with a history of diabetic nephropathy. Diabetes 46:882eC886, 1997
Ji L, Malecki M, Warram JH, Yang Y, Rich SS, Krolewski AS: New susceptibility locus for NIDDM is localized to human chromosome 20q. Diabetes 46:876eC881, 1997
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 20q, proximal to the phosphoenolpyruvate carboxykinase gene. Hum Mol Genet 6:1401eC1408, 1997
Fajans SS, Bell GI, Polonsky KS: Molecular mechanisms and clinical pathophysiology of maturity-onset diabetes of the young. N Engl J Med 345:971eC980, 2001
Winter WE: Newly defined genetic diabetes syndromes: maturity onset diabetes of the young. Rev Endocr Metab Disord 4:43eC51, 2003
Nakhei H, Lingott A, Lemm I, Ryffel GU: An alternative splice variant of the tissue specific transcription factor HNF4alpha predominates in undifferentiated murine cell types. Nucleic Acid Res 26:497eC504, 1998
Boj SF, Parrizas M, Maestro MA, Ferrer J: A transcription factor regulatory circuit in differentiated pancreatic cells. Proc Natl Acad Sci U S A 98:14481eC14486, 2001
Eeckhoute J, Moerman E, Bouckenooghe T, Lukoviak B, Pattou F, Formstecher P, Kerr-Conte J, Vandewalle B, Laine B: Hepatocyte nuclear factor 4 alpha isoforms originated from the P1 promoter are expressed in human pancreatic beta-cells and exhibit stronger transcriptional potentials than P2 promoter-driven isoforms. Endocrinology 144:1686eC1694, 2003
Hansen SK, Parrizas M, Jensen ML, Pruhova S, Ek J, Boj SF, Johansen A, Maestro MA, Rivera F, Eiberg H, Andel M, Lebl J, Pedersen O, Ferrer J, Hansen T: Genetic evidence that HNF-1alpha-dependent transcriptional control of HNF-4alpha is essential for human pancreatic beta cell function. J Clin Invest 110:827eC833, 2002
Stoffel M, Duncan SA: The maturity-onset diabetes of the young (MODY1) transcription factor HNF4alpha regulates expression of genes required for glucose transport and metabolism. Proc Natl Acad Sci U S A 94:13209eC13214, 1997
Wang H, Maechler P, Antinozzi PA, Hagenfeldt KA, Wollheim CB: Hepatocyte nuclear factor 4alpha regulates the expression of pancreatic beta-cell genes implicated in glucose metabolism and nutrient-induced insulin secretion. J Biol Chem 275:35953eC35959, 2000
Damcott CM, Hoppman N, Ott SH, Reinhart LJ, Wang J, Pollin TI, O’connell JR, Mitchell BD, Shuldiner AR: Polymorphisms in both promoters of hepatocyte nuclear factor 4a are associated with type 2 diabetes in the Amish. Diabetes 53:3337eC3341, 2004
Silander K, Mohlke KL, Scott LJ, Peck EC, Hollstein P, Skol AD, Jackson AU, Deloukas P, Hunt S, Stavrides G, Chines PS, Erdos MR, Narisu N, Conneely KN, Li C, Fingerlin TE, Dhanjal SK, Valle TT, Bergman RN, Tuomilehto J, Watanabe RM, Boehnke M, Collins FS: Genetic variation near the hepatocyte nuclear factor-4 gene predicts susceptibility to type 2 diabetes. Diabetes 53:1141eC1149, 2004
Love-Gregory LD, Wasson J, Ma J, Jin CH, Glaser B, Suarez BK, Permutt MA: A common polymorphism in the upstream promoter region of the hepatocyte nuclear factor-4 gene on chromosome 20q is associated with type 2 diabetes and appears to contribute to the evidence for linkage in an Ashkenazi Jewish population. Diabetes 53:1134eC1140, 2004
Hara K, Boutin P, Mori Y, Tobe K, Dina C, Yasuda K, Yamauchi T, Otabe S, Okada T, Eto K, Kadowaki H, Hagura R, Akanuma Y, Yazaki Y, Nagai R, Taniyama M, Matsubara K, Yoda M, Nakano Y, Tomita M, Kimura S, Ito C, Froguel P, Kadowaki T: Genetic variation in the gene encoding adiponectin is associated with increased risk of type 2 diabetes in the Japanese population. Diabetes 51:536eC550, 2002
Alberti KG, Zimmet PZ: Definition, diagnosis and classification of diabetes mellitus and its complications. Part 1: diagnosis and classification of diabetes mellitus provisional report of a WHO consultation. Diabet Med 15:539eC553, 1998
Martin ER, Lai EH, Gilbert JR, Rogala AR, Afshari AJ, Riley J, Finch KL, Stevens JF, Livak KJ, Slotterbeck BD, Slifer SH, Warren LL, Conneally PM, Schmechel DE, Purvis I, Pericak-Vance MA, Roses AD, Vance JM: SNPing away at complex diseases: analysis of single-nucleotide polymorphisms around APOE in Alzheimer disease. Am J Hum Genet 67:383eC394, 2000
Bagwell AM, Bento JL, Mychaleckyj JC, Freedman BI, Langefeld CD, Bowden DW: Genetic analysis of HNF4A polymorphisms in Caucasian-American type 2 diabetes. Diabetes 54:1185eC1190, 2005
Winckler W, Graham RR, de Bakker PI, Sun M, Almgren P, Tuomi T, Gaudet D, Hudson TJ, Ardlie KG, Daly MJ, Hirschhorn JN, Groop L, Altshuler D: Association testing of variants in the hepatocyte nuclear factor 4 gene with risk of type 2 diabetes in 7,883 people. Diabetes 54:886eC892, 2005(Kazuo Hara, Momoko Horiko)