Diabetes 54:3035-3039, 2005 © 2005 by the American Diabetes Association, Inc.
Variants in Hepatocyte Nuclear Factor 4
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
(HNF4
) gene are associated with type 2 diabetes in Finns and Ashkenazi Jews. Previous studies in both populations have reported linkage to type 2 diabetes near the HNF4
locus on chromosome 20q12-13. To investigate whether HNF4
is a diabetes susceptibility gene in Pima Indians, a population with the highest reported prevalence of type 2 diabetes but with no evidence for linkage of the disease on chromosome 20q, 19 SNPs across the promoter and coding region of HNF4
were genotyped for association analysis. In a group of 1,037 Pima Indians (573 diabetic and 464 nondiabetic subjects), three SNPs in HNF4
(rs3212183 and rs2071197 located in introns 3 and 1, respectively, and rs6031558, an extremely rare SNP located in the P2 promoter region) were modestly associated with type 2 diabetes (rs3212183 odds ratio [OR] 1.34 [95% CI 1.07–1.67], P = 0.009; rs2071197 1.34 [1.07–1.66], P = 0.008; and rs6031558 3.18 [1.03–9.84], P = 0.04, adjusted for age, sex, birth year, heritage, and family membership). We conclude that variants in HNF4
do not appear to be major determinants for type 2 diabetes in Pima Indians; however, HNF4
may have a minor role in type 2 diabetes susceptibility within this Native American population.
Abbreviations:
FUSION, Finnish-U.S. Investigation of Type 2 Diabetes; HNF4
, hepatocyte nuclear factor 4
; SNP, single nucleotide polymorphism
Hepatocyte nuclear factor 4
(HNF4
) is a transcription factor expressed in many tissues including liver and pancreas (1). Its expression is controlled by two alternative promoters, P1 and P2, which produce at least nine transcripts. Transcripts from the P2 promoter are the predominant isoforms in pancreatic ß-cells (2). HNF4
has been suggested to have a role in metabolic pathways involved in glucose metabolism and insulin secretion (3,4). Mutations in both the coding and regulatory regions of HNF4
have been associated with maturity-onset diabetes of young (MODY type I), a dominantly inherited, early-onset form of type 2 diabetes (5).
Several genome-wide scans for type 2 diabetes susceptibility loci have identified linkage on chromosome 20q12-13 in a region that encompasses the HNF4
locus (6–10). Recent fine-mapping follow-up studies independently identified multiple variants within the HNF4
gene that are associated with type 2 diabetes in the Finnish-U.S. Investigation of Type 2 Diabetes (FUSION) study population and in the Ashkenazi Jewish population (11,12). In both studies, the strongest association was observed with four SNPs (rs4810424, rs1884613, rs1882614, and rs2144908) from the P2 promoter region (11,12). However, several SNPs in the P1 promoter region are also modestly associated with type 2 diabetes in the FUSION study population but not in the Ashkenazi Jewish cohort. To investigate whether variants in HNF4
are associated with type 2 diabetes in Pima Indians who have a high prevalence of type 2 diabetes but show no evidence for linkage to type 2 diabetes on chromosome 20q13, a total of 19 SNPs across a
79-kb region encompassing both the P2 and P1 promoters and the coding region of HNF4
were genotyped in this population.
| RESEARCH DESIGN AND METHODS |
|---|
|
|
|---|
25 years of age (mean age of onset 16.0 ± 0.5 years and mean BMI 35.0 ± 0.5 kg/m2) were selected for sequencing of the HNF4
coding region. Variants were genotyped for a family-based association study in 1,037 Pima Indians (from 332 nuclear families) who are participants of our ongoing longitudinal study of the etiology of type 2 diabetes among the Gila River Indian Community in Arizona (13,14). Among these 1,037 subjects, 573 had type 2 diabetes (64.4% women, mean BMI 34.4 ± 8.5 kg/m2, mean study age 44.7 ± 13.0 years, and mean age at diagnosis 33.8 ± 10.8 years) and 464 were nondiabetic (43.5% women, mean BMI 33.5 ± 8.4 kg/m2, and mean study age 34.0 ± 11.0 years). Diabetes status was determined according to the criteria of the World Health Organization (15).
Two variants were additionally genotyped in an independent case/control group of 170 full-heritage, non–first-degree–related Pima Indians for analysis of early-onset diabetes. The case group consisted of the same 96 subjects used for the initial sequencing of the HNF4
gene, as described above. The control group consisted of 74 subjects with normal glucose tolerance (mean age 54.0 ± 0.4 years at last medical exam and mean BMI 36.0 ± 1 kg/m2).
Metabolic measurements of risk factors for type 2 diabetes in normal glucose-tolerant subjects.
Among the 1,037 subjects who were genotyped for the type 2 diabetes association analysis, 179 subjects were normal glucose-tolerant, full-heritage Pima Indian subjects who had undergone detailed metabolic testing as inpatients in our clinical research center. Glucose tolerance was determined by a 75-g oral glucose tolerance test with measurements of fasting 30-, 60-, 120-, and 180-min plasma glucose and insulin concentrations (15). To measure the acute insulin response, blood samples were collected before a 25-g glucose bolus infusion and at 3, 4, 5, 6, 8, and 10 min after infusion. The acute insulin response was calculated as the mean increment in plasma insulin concentrations from 3 to 5 min (16). Insulin sensitivity was assessed using the hyperinsulinemic-euglycemic clamp technique as described previously (16,17). Body composition was estimated by underwater weighing until January 1996 and by dual-energy X-ray absorptiometry (DPX-1; Lunar Radiation) thereafter (18).
Sequencing of HNF4
and genotyping of SNPs.
Sequencing was done using Big Dye terminator (Applied Biosystems) on an automated DNA capillary sequencer (model 3730; Applied Biosystems). Variants identified by either sequencing or databases were genotyped using the TaqMan Allelic Discrimination (AD) Assay (Applied BioSystems) on an ABI Prism 7700 (Applied BioSystems).
Statistical analysis.
Associations were assessed using the statistical analysis system of the SAS institute (Cary, NC). The relationship between trait and marker was generally assessed via an ordered categorical variable representing genotype (i.e., an additive model). For continuous variables, the general estimating equation procedure was used to adjust for covariates. These analyses account for the correlation among family members (i.e., siblings). In the family-based study, the association with diabetes was assessed by logistic regression with control for age, sex, heritage, and year of birth. These analyses were also performed with the general estimating equation to account for the correlation among siblings. In the case-control study of early-onset diabetes, the association was assessed by logistic regression with control for sex and heritage. Haplotype frequencies for pairs of SNPs were calculated in all 1,037 individuals with the Estimating Haplotypes (EH) program (19). D' was calculated as a measure of allelic association and
2 (r2) as a measure of concordance. Association between traits and individual haplotypes were examined with a modification of the zero-recombinant haplotyping procedure, as described previously (20).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
in Pima Indians (Fig. 1). Two of these variants were novel (Met86Arg and IVS9 + 378T-C) and two were present in databases (Thr130Ile = rs1800961 and C/T in the 3' untranslated region = rs6130615). In addition, 10 SNPs (rs1884614, rs2144908, rs2425637, rs2425640, rs6031552, rs6031558, rs3212183, rs1885088, rs1028583, and rs3818247) were selected for genotypic analysis in Pima Indians based on results from the FUSION and Ashkenazi Jewish studies of HNF4
(11,12). Two of these variants, rs1884614 and rs2144908, were associated with type 2 diabetes in both the FUSION and Ashkenazi Jewish studies, whereas the remaining eight variants were associated with type 2 diabetes in only one of the studies. These SNPs span a
79-kb region of HNF4
including the P2 and P1 promoters and the entire coding region of HNF4
. To approximate an evenly spaced dense SNP map across the HNF4
locus, with an average inter-SNP distance of 3–4 kb, seven additional database SNPs (hcv11183312, rs6017335, rs4364072, rs4812831, rs3212198, rs6103731, and rs1028583) were also selected, based on their relative positions, for genotyping in Pima Indians (Fig. 1).
|
variants in 1,037 Pima subjects are given in Table 1. All SNPs were in Hardy-Weinberg equilibrium. Three SNPs, rs3212183, rs2071197, and rs6031558, were modestly associated with type 2 diabetes. SNP rs3212183, positioned in intron 3, had a minor allele (C) frequency of 0.43 in diabetic subjects versus 0.37 in nondiabetic subjects (OR 1.34 [95% CI 1.07–1.67], P = 0.009). This variant was also modestly associated with type 2 diabetes in the FUSION study, where the C allele was also the risk allele. SNP rs2071197, which was not reported in either the FUSION or Ashkenazi Jewish study, was also associated with type 2 diabetes, with an allele (G) frequency of 0.51 in diabetic subjects versus 0.45 in nondiabetic subjects (1.34 [1.07–1.66], P = 0.008). These two SNPs (rs3212183 and rs2071197) were in very high linkage disequalibrium (D' = 0.99) in Pima Indians and were also in high, but not complete, genotypic concordance (
2 = 0.69). SNP rs6031558, positioned in the P2 promoter region, had a minor allele (C) frequency of 0.02 in diabetic subjects versus 0.01 in nondiabetic subjects. Although this SNP had a significant association with type 2 diabetes (3.18 [1.03–9.84], P = 0.04), the rare frequency of this variant makes this analysis unreliable.
|
25 years). Neither SNP was associated with early-onset diabetes in this case/control study (P = 0.24, OR 0.75 [95% CI 0.47–0.121], and P = 0.56, 0.87 [0.56–1.37] for rs3212183 and rs2071197, respectively).
We further investigated whether rs3212183 and rs2071197, which were in extremely high linkage disequilibrium and associated with type 2 diabetes, were associated with metabolic risk factors that predict type 2 diabetes in 179 normal glucose-tolerant, full-heritage Pima Indians who had undergone detailed metabolic testing. The various tests assessed insulin action, insulin secretion, and rates of glucose turnover. Neither SNP was associated with various measures of insulin action (data only shown for rs3212183 in Table 2). Although HNF4
is known to have a primary role in insulin secretion, there was no association between these SNPs and insulin secretion, as measured by the acute insulin response following an intravenous bolus of glucose. It is possible that we lacked the power to detect a subtle association due to the relatively small sample number (n = 179). In support of this, subjects homozygous for the type 2 diabetes risk allele did have a trend toward lower insulin secretion compared with subjects heterozygous or homozygous for the nonrisk allele (Table 2).
|
The pattern of linkage disequilibrium across the HNF4
locus was evaluated in Pima Indians by a pair-wise analysis of genotypes from the 19 variants in the 1,037 subjects (Fig. 2). Similar to that observed in the Finnish and Ashkenazi Jewish populations, SNPs in the P2 promoter region were in extremely high linkage disequilibrium. Linkage disequilibrium did not decay until
33 kb from the beginning of the P2 promoter, with the exception of a rare SNP rs6031558 that is positioned in the middle of this block and breaks the block due to low linkage disequilibrium (<0.2) with rs6017335 and rs4364072. SNPs in exons 4–11 were also in high linkage disequilibrium. However, SNPs across the P1 promoter and exons 1–3 tended to have lower linkage disequilibrium among themselves but still had D' values of >0.5, with some SNPs in the P2 region. When SNPs containing common variants (minor allele frequency >0.05) were analyzed, there were five haplotype "blocks" among these SNPs, as defined by contiguous SNPs with D' >0.9. Analyses of the individual haplotypes for SNPs within these blocks showed that the haplotype consisting of "G" at rs2071197 and "C" at rs3212183 had a higher prevalence of diabetes compared with the other combined haplotypes (OR 1.37 [95% CI 1.10–1.72], P = 0.02), a finding that largely reflects the single marker associations. There were no statistically significant associations with any of the other haplotypes (data not shown).
|
are consistent between different ethnic groups, while others are quite different. For example, the allele (C) for rs3212183 had a frequency of 0.43 in both diabetic Pima Indians and diabetic Finns and a frequency of 0.37 and 0.38 in nondiabetic Pima Indians and Finns, respectively. In contrast, for rs1884614 the risk allele (T) in the Finns and Ashkenazi Jews had a frequency of 0.20 and 0.27, respectively, but the T allele in Pima Indians had a frequency of 0.83. This suggests that different populations have different HNF4
haplotypes. A difference in allele frequencies of HNF4
variants between ethnic groups may also reflect the various population-attributed risks to type 2 diabetes due to these variants.
The strongest associations with type 2 diabetes in both the FUSION and the Ashkenazi Jews studies were observed with SNPs in the P2 promoter of HNF4
. Recent studies (21,22) in the Amish and U.K. populations have also replicated associations of variants in HNF4
with type 2 diabetes. In the Amish study (21), the strongest association was with rs2425640 located in the P1 promoter. The U.K. study reported significant associations with two SNPs (rs4810424 and rs2144908) that tag the P2 promoter (22). In Pima Indians, the two common SNPs associated with type 2 diabetes were in introns 1 and 3. The fact that different SNPs in the HNF4
region are associated with diabetes in different populations suggests that none of these alleles themselves are causative functional variants but that they may be in linkage disequilibrium with a nearby functional allele. Alternatively, some of these alleles may be causative, but allelic heterogeneity across populations may make their identification difficult. Given the previous evidence implicating HNF4
variants in diabetes, we considered that it is desirable to preserve statistical power. Therefore no corrections of multiple comparisons were performed. Nonetheless, one cannot exclude chance as a potential explanation for the associations observed in the present study.
Although the sample size of the present study is fairly large in comparison with many other studies, it may nonetheless have low statistical power to detect variants with modest effects. We estimate that the power of the study is >70% to detect a significant (P < 0.05) association with a causal variant with minor allele frequency >0.05 accounting for 1% of the variance in liability to diabetes. Given an allele frequency of 0.2, this corresponds to an OR of 1.32, approximately the effect observed in the FUSION study. It is noteworthy that the ORs in the FUSION and Ashkenazi studies are likely an overestimation of the true effect due to testing of multiple variants within this region, and the power to detect a true, smaller, underlying effect will be lower. Examination of the 95% CIs for the ORs in Table 1 provides an evaluation of the extent of effect sizes most compatible with the data.
Although our genome-wide linkage scan in Pima Indians did not detect evidence for linkage to type 2 diabetes on chromosome 20q12-13, this association study showed that variants in introns 1 and 3 of HNF4
are modestly associated with type 2 diabetes. These results suggest that the HNF4
locus, or another gene near this locus, may be a common genetic determinant of type 2 diabetes in multiple populations, but its relative contribution may differ between populations.
| ACKNOWLEDGMENTS |
|---|
Received for publication January 3, 2005 and accepted in revised form July 12, 2005
| REFERENCES |
|---|
|
|
|---|
gene predicts susceptibility to type 2 diabetes. Diabetes 53:1141–1149, 2004
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:1134–1140, 2004
are associated with type 2 diabetes in the Amish. Diabetes 53:3337–3341, 2004
P2 promoter are associated with type 2 diabetes in the U.K. population. Diabetes 53:3002–3006, 2004This article has been cited by other articles:
![]() |
I. Barroso, J. Luan, E. Wheeler, P. Whittaker, J. Wasson, E. Zeggini, M. N. Weedon, S. Hunt, R. Venkatesh, T. M. Frayling, et al. Population-Specific Risk of Type 2 Diabetes Conferred by HNF4A P2 Promoter Variants: A Lesson for Replication Studies Diabetes, November 1, 2008; 57(11): 3161 - 3165. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Johansson, H. Raeder, S. A Eide, K. Midthjell, K. Hveem, O. Sovik, A. Molven, and P. R. Njolstad Studies in 3,523 Norwegians and Meta-Analysis in 11,571 Subjects Indicate That Variants in the Hepatocyte Nuclear Factor 4{alpha} (HNF4A) P2 Region Are Associated With Type 2 Diabetes in Scandinavians Diabetes, December 1, 2007; 56(12): 3112 - 3117. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Watanabe, M. H. Black, A. H. Xiang, H. Allayee, J. M. Lawrence, and T. A. Buchanan Genetics of Gestational Diabetes Mellitus and Type 2 Diabetes Diabetes Care, July 1, 2007; 30(Supplement_2): S134 - S140. [Full Text] [PDF] |
||||
![]() |
D. M. Lehman, D. K. Richardson, C. P. Jenkinson, K. J. Hunt, T. D. Dyer, R. J. Leach, R. Arya, H. E. Abboud, J. Blangero, R. Duggirala, et al. P2 Promoter Variants of the Hepatocyte Nuclear Factor 4{alpha} Gene Are Associated With Type 2 Diabetes in Mexican Americans Diabetes, February 1, 2007; 56(2): 513 - 517. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Benoit, A. Cooney, V. Giguere, H. Ingraham, M. Lazar, G. Muscat, T. Perlmann, J.-P. Renaud, J. Schwabe, F. Sladek, et al. International Union of Pharmacology. LXVI. Orphan Nuclear Receptors Pharmacol. Rev., December 1, 2006; 58(4): 798 - 836. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. L. Bonnycastle, C. J. Willer, K. N. Conneely, A. U. Jackson, C. P. Burrill, R. M. Watanabe, P. S. Chines, N. Narisu, L. J. Scott, S. T. Enloe, et al. Common Variants in Maturity-Onset Diabetes of the Young Genes Contribute to Risk of Type 2 Diabetes in Finns Diabetes, September 1, 2006; 55(9): 2534 - 2540. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. E. Pavkov, P. H. Bennett, W. C. Knowler, J. Krakoff, M. L. Sievers, and R. G. Nelson Effect of youth-onset type 2 diabetes mellitus on incidence of end-stage renal disease and mortality in young and middle-aged Pima Indians. JAMA, July 26, 2006; 296(4): 421 - 426. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Weissglas-Volkov, A. Huertas-Vazquez, E. Suviolahti, J. Lee, C. Plaisier, S. Canizales-Quinteros, T. Tusie-Luna, C. Aguilar-Salinas, M.-R. Taskinen, and P. Pajukanta Common Hepatic Nuclear Factor-4{alpha} Variants Are Associated With High Serum Lipid Levels and the Metabolic Syndrome. Diabetes, July 1, 2006; 55(7): 1970 - 1977. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Raeder, L. Bjorkhaug, S. Johansson, K. Mangseth, J. V. Sagen, A. Hunting, I. Folling, O. Johansen, M. Bjorgaas, P. N. Paus, et al. A Hepatocyte Nuclear Factor-4{alpha} Gene (HNF4A) P2 Promoter Haplotype Linked With Late-Onset Diabetes: Studies of HNF4A Variants in the Norwegian MODY Registry Diabetes, June 1, 2006; 55(6): 1899 - 1903. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Diabetes | Diabetes Care | Clinical Diabetes | Diabetes Spectrum |