DOI: 10.2337/db06-0189 © 2006 by the American Diabetes Association
Common Polymorphisms in the Promoter of the Visfatin Gene (PBEF1) Influence Plasma Insulin Levels in a French-Canadian Population
1 Department of Human Genetics, McGill University, Montréal, Québec, Canada Address correspondence and reprint requests to Dr. James C. Engert, McGill University, Division of Cardiology, Royal Victoria Hospital, H7.11, 687 Pine Ave. West, Montréal, Québec, Canada H3A 1A1. E-mail: jamie.engert{at}mcgill.ca
Abbreviations:
apoB, apolipoprotein B; CHD, coronary heart disease; Coup-TF, chick ovalbumin upstream promoter–transcription factor; DR, direct repeat; HRE, hormone response element; LD, linkage disequilibrium; PPAR, peroxisome poliferator–activated receptor; SNP, single nucleotide polymorphism; QFS, Quebec Family Study
The adipokine visfatin (PBEF1) exhibits insulin-mimetic effects and correlates strongly with visceral adiposity. We sequenced visfatin gene exons and 1,480 bp of the promoter in 23 individuals, including 18 individuals from the Quebec Family Study (QFS) with varying degrees of abdominal visceral fat, assessed by computed tomography, and 5 individuals from the Saguenay-Lac-Saint-Jean region of Québec. We identified a synonymous polymorphism in exon 7 (SER301SER) but no nonsynonymous mutations. We observed an additional 10 polymorphisms, including 5 intronic, 4 within the promoter, and 1 within the 3' untranslated region. Further promoter sequencing (816 bp) identified five additional single nucleotide polymorphisms (SNPs) in the QFS population. To investigate the role of visfatin gene variants in obesity-related phenotypes, we genotyped a total of 13 SNPs in the promoter region of the gene. From these, we analyzed the seven common SNPs in the QFS sample (918 participants from 208 families). A significant association was found between two SNPs (rs9770242 and rs1319501), in perfect linkage disequilibrium, and fasting insulin levels (P = 0.002). These SNPs were also associated with fasting glucose (P 0.02). In addition, a more distal SNP (rs7789066) was significantly associated with the apolipoprotein B component of VLDL (P = 0.012). One of the biggest threats to global health is the rising epidemic of obesity. Observed first in the wealthier Western world (1), it has now been shown to accompany the Westernization of many rising economies (2). Obesity frequently occurs with a constellation of clinical parameters that are all poor prognostic indicators for diabetes and coronary heart disease (CHD). Obesity can be defined as an excessive amount of body fat. However, the distribution of fat may also be important for determining risk. The existence of an association between visceral adipose tissue accumulation and an increased risk of type 2 diabetes and CHD is well established (3). Visceral obesity is one of several features of the metabolic syndrome, a clustering of risk factors associated with CHD, which includes insulin resistance, dyslipidemia, and hypertension (4). Insulin resistance plays a central role in the pathogenesis of the metabolic abnormalities associated with visceral adiposity. Chemical messengers secreted by adipocytes, termed adipokines, are thought to modulate insulin action (5). The newly identified adipokine visfatin, also known as pre–B-cell colony enhancing factor (PBEF1) (6), was discovered to be preferentially expressed in visceral as opposed to subcutaneous adipose tissue by use of a subtraction library (7). However, visfatin/PBEF1 is expressed in a variety of other tissues and is especially high in the liver (6).
Intriguingly, not only were circulating levels of visfatin shown to correlate strongly with visceral adiposity, but visfatin also exhibited insulin-mimetic effects that were mediated by the insulin receptor and resulted in lower blood glucose levels in mice (7). In addition, chronic administration of visfatin, in mice, lowered both plasma glucose and insulin levels (7). Visfatin stimulated triglyceride synthesis and accumulation in preadipocytes (7) and appears to be regulated by other adipokines, namely interleukin-6 (8) and tumor necrosis factor- These findings make visfatin a strong candidate for mediating the complex interplay between visceral obesity and associated metabolic complications. Thus, we hypothesized that genetic variants at the visfatin gene locus may contribute to interindividual variation in visceral obesity or insulin resistance.
We undertook a thorough examination of the visfatin gene locus for coding and noncoding variants. In the polymorphism discovery phase of this study, we sequenced all coding regions and 1,480 bp of the promoter of the visfatin gene in 18 individuals with varying degrees of visceral adipose tissue (assessed by computed tomography) from the Quebec Family Study (QFS) (10), as well as 5 individuals with various waist-to-hip ratios selected from ongoing CHD studies in the Saguenay-Lac-Saint-Jean region of Québec. Our initial sequencing identified a total of 11 polymorphisms (Figs. 1 and online appendix Table 1 [available at http://diabetes.diabetesjournals.org]). We observed one synonymous polymorphism in exon 7 (g.21327 A>C, SER301SER, rs2302559) but no nonsynonymous polymorphisms. This is in accordance with dbSNP, which lists the SER301SER variant but no other coding single nucleotide polymorphisms (SNPs). Two deletions were discovered: one in intron 6 (g.21196 delGAT) and the other in the 3' untranslated region (g.34283 delACATT). We also found four SNPs within the proximal promoter (g.-1001 T>G [rs9770242], g.-948 G>T, g.-423 A>G [rs1319501], and g.-295 G>C), while the remaining polymorphisms were intronic: one in intron 4 (g.12464 C>G), one in intron 6 (g.21179 T>C), and two in intron 8 (g.23516 A>C and g.30370 A>C). Minor allele frequencies varied from 0.02 to 0.50 (online appendix Table 1).
Because no nonsynonymous polymorphisms were discovered, we initiated a study of common (minor allele frequency >0.05) promoter polymorphisms (g.-1001 T>G [rs9770242], g.-948 G>T, and g.-423 A>G [rs1319501]), which may influence expression of the visfatin gene. Since these SNPs are in close proximity to one another, we genotyped them by sequencing. In addition, a comprehensive search of publicly available databases (HapMap and Perlegen) was performed, and four additional common SNPs, from the more distal promoter, were chosen for genotyping (rs2110385, rs1737358, rs11977021, and rs7789066). During the follow-up sequence-based genotyping, we discovered five rare SNPs (g.-882 G>C, g.-520 G>A, g.-404 C>A, g.-356 G>A, and g.-348 C>A) with allele frequencies between 0.002 and 0.009 (Table 1, Fig. 1). None of the common visfatin SNPs deviated significantly from Hardy-Weinberg equilibrium (these calculations were not performed for rare SNPs) (Table 1). To examine whether the common promoter SNPs of the visfatin gene explain some of the interindividual variation of clinical measures of obesity and related traits, we conducted a family-based association test on the QFS sample. Characteristics of the study subjects are shown in Table 2. Using a variance-components approach with a dominance parameter included in the model (adjusting for age and sex), a significant association was found between fasting plasma insulin levels and two SNPs, rs9770242 and rs1319501 (P = 0.002 for both SNPs) (Table 3). This association remained significant when corrected for the multiple testing of seven SNPs (empirical P value of 0.011 for both rs9770242 and rs1319501). In addition, an association was observed for both variants and fasting glucose levels (P = 0.020 and 0.017, respectively). Thus, individuals homozygous for the G-allele, of either the rs1319501 and rs9770242 variants, had lower fasting plasma insulin levels and lower fasting glucose levels. These two SNPs appear to be in perfect linkage disequilibrium (LD) with each genotype, predicting the other in all individuals for which genotype information was available. The discrepancy in the observed P values is a consequence of the distribution of missing genotypes.
We also observed a significant association between the rs7789066 variant and the apolipoprotein B (apoB) component of VLDL particles (P = 0.012) (Table 3). This result was comparable between the two models tested (Table 3). In addition, the rs11977021 variant was found to be significantly associated, under the additive model of inheritance, with total plasma cholesterol levels and LDL cholesterol (P = 0.03 and P = 0.049, respectively) (Table 3). Finally, the g.-948 G>T variant was found to be significantly associated, when a dominance parameter was included in the model, with fasting insulin levels (P = 0.011) and total apoB and LDL apoB levels (P = 0.031 and P = 0.033, respectively) (Table 3). (This variant is located between the rs9770242 and rs1319501 variants and is in significant LD with them [data not shown].)
Our sequencing of the coding regions of the visfatin gene in 23 unrelated individuals from Québec failed to identify any coding nonsynonymous variants. Human visfatin shares 95% sequence identity with both the rat and mouse protein sequences (data not shown). The absence of common missense mutations and the conservation of the visfatin protein suggests an important biological function, which is supported by the observation that homozygous knockout mice were embryonic lethals (7). The finding that certain 5' flanking variants of the visfatin gene (rs1319501 and rs9770242) are associated with plasma insulin levels (P = 0.002) and with plasma glucose levels (P < 0.02) is consistent with its demonstrated insulin-mimetic properties (7). Interestingly, Berndt et al. (11) found a significant correlation between visfatin gene expression in visceral fat tissue and plasma insulin levels, although these same authors observed no correlation between plasma visfatin and plasma insulin levels. In a recent study (12), plasma visfatin levels were found to be increased in individuals with type 2 diabetes and there was a significant association between plasma visfatin levels and plasma insulin levels. In our study, we did not observe a significant association between any of the visfatin promoter SNPs and type 2 diabetes (data not shown). However, the number of individuals with type 2 diabetes in the QFS sample is low (51 individuals). The change in plasma insulin levels associated with the rs1319501 and rs9770242 variants are consistent with the alteration of a transcription factor binding site at one of these SNPs that effects visfatin gene expression. Because of the perfect LD, it is not possible to say which is the functional variant. Of note, the G-allele of the rs9770242 (–1001 T>G) variant was recently associated with acute lung injury and sepsis, though a transfection assay did not reveal a significant change in visfatin/PBEF1 gene expression in human microvascular endothelial cells from the lung cells (13). However, these authors did not investigate the rs1319501 variant, for which perhaps a much stronger argument can be made. This SNP occurs within an imperfect nuclear hormone response element (HRE), which contains a direct repeat (DR) of the concensus hexamer (A/G)GGTCA separated by a range of nucleotide spacings (Pu-GGTCANiPu-GGTCA); in this case, a spacing of six nucleotides (DR6) (14,15). The receptor family that recognizes these HREs include the peroxisome poliferator–activated receptors (PPARs), the chick ovalbumin upstream promoter–transcription factors (Coup-TFs), and the vitamin D receptor (14,15).
In a recent study, PPAR-
The finding that the rs7899066 variant was suggestively associated with the apoB component of VLDL (P = 0.012) is intriguing because diabetic patients and insulin-resistant individuals have increased VLDL apoB and triglyceride levels (19). The rs7899066 SNP is located within a nuclear factor of activated T-cells (NFAT) binding sequence ([T/A]GGAAAA) (20). The NFAT transcription factors are involved in adipocyte differentiation (21,22) and can regulate the expression of fatty acid–binding protein aP2 (21), PPAR- Finally, the rs11977021 variant, which was suggestively associated with total cholesterol and LDL cholesterol levels (P < 0.049), is also located within an HRE sequence; in this case, a perfect DR separated by four nucleotides (DR4), which is a binding site for the thyroid hormone receptor, and is another member of the nuclear receptor family mentioned above (14). In addition, the g.-948 G>T variant was found to be associated with fasting plasma insulin levels, total serum apoB, and LDL apoB levels. These findings, while only suggestive, may indicate that visfatin is involved in other aspects of the dyslipidemia associated with insulin resistance, such as the presence of small dense LDL particles. It is worth noting that we obtained a greater number of significant results when our model included a dominance parameter (Table 3). This included our best result, being the two SNPs that were strongly associated with fasting insulin. Because we have used a family-based design with families from Québec (the site of a well-documented founder effect), this study may have increased power to detect genetic effects that exhibit a dominance component. Further, this could be due to the increased probability of sharing two alleles identical by descent. Our results demonstrate a significant association between common promoter polymorphisms in the visfatin gene and quantitative measures of insulin resistance among French Canadians in Québec. While we did not see a direct association with type 2 diabetes; our power to detect an association is low in the QFS population. The association between fasting insulin and glucose levels and the rs9770242 and rs1319501 variants, as well as the association between the rs7899066 and rs11977021 variants and components of diabetic dyslipidemia, warrant additional genetic studies in larger obesity and type 2 diabetes population samples, as well as functional characterization of these promoter variants.
The QFS is a long-term study of French-Canadian families from Québec City and its surrounding area and has been described in detail (10). Briefly, the current study involved 208 families and a total of 918 individuals for which phenotypic information and visfatin genotypes were available from QFS. Though QFS contains some extended families, the majority (165 families) are nuclear families. The five unrelated individuals from the Saguenay-Lac-Saint-Jean region of Québec were selected from ongoing CHD studies.
Phenotypes.
Sequencing.
Genotyping.
Restriction fragment–length polymorphisms.
TaqMan assays.
Statistical analysis.
SNP selection, transcription factor binding site identification, and protein sequence alignments. Potential transcription factor binding sites were identified on forward and backward strands of genomic DNA using the transcription element search system TESS program (available at http://cbil.upenn.edu/tess). Between-species protein sequence alignments were performed using the BLAST algorithm (available at http://www.ncbi.nlm.nih.gov/BLAST/).
After the current article was submitted, we became aware of similar work that also demonstrates a significant association between the visfatin promoter polymorphism (g.-948 G>T) and fasting plasma insulin levels (26).
J.C.E. and M.-C.V. are research scholars from the Fonds de la recherche en santé du Québec. D.G. is the chairholder of the Canadian Research Chair in preventive genetics and community genomics. T.J.H. is a recipient of an Investigator Award, from the Canadian Institutes of Health Research, and of the Clinician-Scientist Award in Translational Research, from the Burroughs Wellcome Fund. C.B. is partially supported by the George A. Bray Chair in Nutrition. We thank the individuals who volunteered to participate in the studies in the Saguenay-Lac St. Jean region and the Québec City area, the staff of the Lipid Research Centre in Québec City, and of the Lipid Research Group at Chicoutimi Hospital. Gratitude is expressed to G. Theriault and G. Fournier, L. Allard, M. Chagnon, and C. Leblanc for their contributions to the recruitment and data collection of the QFS. We thank C. Doré for technical assistance and L. Coderre, T. Pastinen, A. Sniderman, and K. Cianflone for helpful discussions.
Additional information for this article can be found in an online appendix at http://diabetes.diabetesjournals.org. 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. Received for publication February 9, 2006 and accepted in revised form June 27, 2006
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||