Abstract
Fibroblast growth factor 21 (FGF-21) is a novel regulator for metabolic syndrome (MetS), diabetes, and obesity. However, no study has been performed on the association of these diseases with FGF-21 gene polymorphism. The aim of the study was to investigate the association of 3′ untranslated region (UTR) single-nucleotide polymorphisms (SNPs) in the FGF-21 gene with MetS, obesity, and diabetes in the Han Chinese population. A total of 291 subjects were recruited from the Han Chinese population in Sichuan province. The genotypes of FGF-21 were determined by polymerase chain reaction–restriction fragment length polymorphism. The genotypes were confirmed by sequencing. No polymorphisms were found in rs11665841 (1 of 291) and rs3745706 (2 of 291). We did not find an association between genotype frequencies of SNP rs11665896 and lipid concentration, glucose concentration, or blood pressure. The TG/GG genotype relative to the TT genotype had an age- and sex-adjusted odds ratio of 1.41 for MetS (p=0.149), 1.84 (p=0.016) for obesity (body mass index ≥25 kg/m2), and 1.19 (p=0.492) for diabetes. Genetic variation of the 3′ UTR of the FGF-21 gene was associated with obesity, however, not with MetS or diabetes.
Introduction
Fibroblast growth factor 21 (FGF-21) is a novel metabolic regulator with multiple beneficial effects on glucose and lipid homeostasis in animal models and humans. It was suggested that FGF-21 could elevate glucose uptake in 3T3-L1 adipocytes and reduce blood glucose and triglyceride concentrations. Moreover, its overexpression in transgenic mice could prevent obesity (Kharitonenkov et al., 2005). FGF-21 was considered potentially therapeutic for obesity and fatty liver disease (Xu et al., 2009). Zhang et al. (2008) found that FGF-21 mRNA expression in both liver and adipose tissues was significantly higher in obese ob/ob mice than in normal-weight ob/ob mice, and FGF-21 expression levels in subcutaneous fat were in positive correlation with circulating FGF-21 concentrations.
MicroRNAs (miRNAs) constitute a growing class of noncoding RNAs that suppress target gene protein translation by combining with 3′ untranslated region (UTR) completely or incompletely (Bartel, 2004). miRNA binding may lead to phenotypic changes including disease susceptibility, cancer risk and drug-resistance of tumor cells and so on (Sethupahty et al., 2008). It is reported that ERα 3′UTR, which is the target binding site of miR-206 in vitro test, contains a C-T variant (rs 93410170) that enhanced binding with miR-206, and then reduced the level of ERa protein (Adams et al., 2007). In addition, SNPs near the miRNA target area also resulted in different phenotypes (Mishra et al., 2007). Therefore, the aim of our study is to investigate the association of 3′ UTR SNPs in FGF-21 gene with MetS, diabetes, and obesity in Chinese Han population.
Materials and Methods
Subjects
Subjects were recruited from Department of Physical Examination Center, West China Hospital, Sichuan University, China. Their age ranged from 50 to 90 years old. Initially, 550 subjects were recruited for this investigation. All subjects participated in a complete physical examination, including family history, history of past illness, and standard anthropometrical measurements. After exclusion of subjects with tumor, anorexia nervosa, serious liver disease, and serious kidney disease, 365 subjects were invited for further investigation, including 75 g oral glucose tolerance test (OGTT) and conventional biochemical tests. Finally, 291 subjects agreed to participate, including 130 MetS subjects, 103 diabetic subjects, 100 obesity subjects, and 141 healthy subjects. The diagnosis of MetS was according to International Diabetes Federation definition, which included central obesity (waist circumference ≥90 cm for men, ≥80 cm for women) plus any two of the following items: (1) triacylglycerol >1.7 mM; (2) high-density lipoprotein (HDL) cholesterol <1.03 mM for men and <1.29 mM for women; (3) systolic blood pressure ≥130 mmHg or diastolic blood pressure ≥85 mmHg or hypertension history; and (4) fasting plasma glucose ≥5.6 mM or diabetes history. The diagnosis of type 2 diabetes was based on 1999 World Health Organization criteria (fasting plasma glucose level ≥7.0 mM or 2-h OGTT plasma glucose level ≥11.1 mM). The World Health Organization's definition of obesity in Asians is body mass index (BMI; in kg/m2) ≥25.
All studies were approved by Ethics Committee at West China Hospital of Sichuan University. Written informed consent was obtained from all subjects before participation.
Measurement of biomarkers
Height, weight, waist circumferences, and blood pressure were measured in the fasting state by the same physician. The fasting and 2-h OGTT plasma glucose, serum triglycerides, total cholesterol, HDL and low-density lipoprotein were determined by Roche P800. Fasting plasma insulin was measured by electrochemiluminescence immunoassay (ECLIA) (Roche E170). Plasma FGF-21 levels were determined by radioimmunoassay (RIA). (Phoenix Pharmaceuticals, Inc.) using 125I-labeled FGF-21 as tracer. Blood samples were collected with aprotinin (0.6 TIU/mL of blood) and detection was centralized.
The FGF-21 antibody was rabbit antihuman FGF-21 IgG, which is specific for human FGF-21; there are no cross-reaction between FGF-21 antibody and human FGF-6, FGF-10, FGF-18, FGF-19, FGF-20, adiponectin, visfatin, leptin, and retinol-binding protein-4 resistin. The linear range is 0.5–8.5 ng/mL, intraassay coefficient of variation (CV) <5%, and inter-CV <14%.
Selection of SNPs and genotyping
There are a total of 17 SNPs in the genomic region from ∼200 bp upstream to ∼200 bp downstream FGF-21. Thirty subjects (containing 10 healthy subjects, 8 MetS subjects, 6 diabetic subjects, and 6 obesity subjects) were randomly selected for screening the 17 SNPs by high-resolution melting. We found that only 3 SNPs had a high frequency (>0.1) of the minor allele in the 30 subjects.
Genomic DNA was extracted from peripheral blood leukocytes using the QIAamp DNA mini kit (Qiagen GmbH) following the manufacturer's instructions. Genotyping of the SNPs was performed by using polymerase chain reaction (PCR)–restriction fragment length polymorphism method. Locus-specific PCR primers were designed using the primer 5.0 software; primer pairs and specific restriction endonuclease are shown in Table 1.
Data are mean±SD.
BMI, body mass index; SBP, systolic blood pressure; DBP, diastolic blood pressure; OGTT, oral glucose tolerance test; LDL, low-density lipoprotein; HDL, high-density lipoprotein; SD, standard deviation.
The rs11665896 PCR primers were 5′ TGTGTGAGTGTC TGAGGGAAG 3′ (forward) and 5′ GAAGTCAAGAGAT GGAGAGCA 3′ (reverse). The rs3745706 PCR primers were 5′ GGGGGTTGTCCGGGCCAAGG 3′ (forward) and 5′ GGC CGAAGCCCCAGCTACGCTT 3′ (reverse). For rs11665896 and rs3745706, the PCR conditions were as follows: initial denaturation at 95°C for 3 min, amplification for 35 cycles of 95°C for 40 s, 64°C for 40 s, 72°C for 40 s. A final extension step at 72°C for 5 min followed the last PCR cycle. The PCR was performed in a final volume of 25 μL that contained 100 ng genomic DNA, 25 pM of each primer, 200 μM dNTP, 0.5 units Amplitaq Gold DNA polymerase (Applied Biosystems), and 2 mM MgCl2. To analyze the rs11665896 and rs3745706 polymorphisms, PCR products were individually digested with 1 unit of BlnI and MboII restriction endonucleases at 37°C for 2 h and fragments obtained were visualized on a 2.5% agarose gel.
The rs11665841 PCR primers were 5′ CAGAGGAGGATA AGAAGAATGAGCA 3′ (forward) and 5′ TGGGCACAGT GACTCAAGC 3′ (reverse). For rs11665841, PCR conditions were as follows: initial denaturation at 95°C for 3 min, amplification for 35 cycles of 95°C for 40 s, 62°C for 40 s, and 72°C for 1 min. A final extension step at 72°C for 5 min followed the last PCR cycle. The rs11665896 PCR products were digested using 1 unit BsrI restriction endonuclease at 65°C for 4 h and fragments obtained were identified on a 2.5% agarose gel.
Statistical analysis
The normal distribution of data was expressed as mean±SD. Hardy–Weinberg equilibrium was computed for the expected genotype distribution. Genotypic and allelic distributions were compared using the Chi-square test. Odds ratios (ORs) were derived from multiple logistic regression analysis, which was carried out to evaluate the effect of the FGF-21 gene polymorphism. All statistical analyses were performed applying SPSS 11.5. A p<0.05 was considered statistically significant.
Results
Clinical and biochemical characteristics and 3′ UTR SNPs of FGF-21 in study subjects
The clinical and biochemical characteristics of all subjects are shown in Table 1. Frequencies of MetS, obesity, and diabetes were 44.67%, 34.36%, and 35.40%, respectively. No polymorphisms were found in rs11665841 (1 of 291) and rs3745706 (2 of 291). For rs11665896, the frequencies of genotypes TT, TG, and GG in our study subjects were 8.25%, 39.86%, and 51.89%, respectively. The genotypes of SNP rs11665896 were confirmed by sequencing. We did not find an association between rs11665896 and biochemical characteristics, as shown in Table 2.
Data are mean±SD.
Association of FGF-21 rs11665896 with MetS
The genotypes of rs11665896 were in Hardy–Weinberg equilibrium in both case and noncase groups (p>0.05 for all analyses). The genotype and allele frequencies in MetS, non-MetS group, and healthy groups are shown in Table 3. The genotype and allele frequencies were not significantly different between the MetS group and non-MetS group. No significant difference in genotype and allele frequencies was found between MetS group and healthy subjects. We also did not find an association between the genotype and allele frequencies of rs11665896 with MetS components (Table 4). As shown in Table 5, when the TG and GG genotypes were combined, the age- and sex-adjusted OR for MetS (vs. the TT genotype) was 1.41 (95% confidence interval [CI]: 0.88, 2.26; p=0.149). We estimated that the sample had a power of 83.40%.
Compared with noncase subjects.
Compared with healthy subjects.
Healthy subjects: non-MetS, non-diabetes, and non-obese.
MetS, metabolic syndrome.
GG or TG versus TT.
Control for age and sex.
OR, odds ratio; CI, confidence interval.
Association of FGF-21 rs11665896 with diabetes
There were 103 diabetic subjects. The genotypic and allelic frequencies in the diabetes, nondiabetes, and healthy groups are shown in Table 3. The genotype and allele frequencies were not significantly different between the diabetes and nondiabetes groups. No significant difference of genotype and allele frequencies was found between the diabetes group and healthy subjects. When the TG and GG genotypes were combined, the OR of diabetes versus the TT genotype was 1.19 (95% CI: 0.73, 1.94; p=0.492) after adjustment for age and sex, as shown in Table 5. We estimated that the sample had a power of 96.08%.
Association of FGF-21 rs11665896 with obesity
The genotype and allele frequencies in the obese group, nonobese group, and healthy subjects are shown in Table 3. The genotype and allele frequencies were not significantly different between the obese and nonobese groups. No significant difference of genotype and allele frequencies was found between the obese group and healthy subjects. When the TG and GG genotypes were combined, the OR of obesity versus the TT genotype was 1.84 (95% CI: 1.12, 3.01; p=0.016) after adjustment for age and sex, as shown in Table 5. We estimated that the sample had a power of 84.61%.
Discussion
FGF-21, which is secreted by hepatic cells, has been proposed as a potent metabolic regulator in glucose and lipid metabolism. Many studies have shown a strong association between FGF-21 and obesity, hyperlipidemia, and diabetes (Kharitonenkov et al., 2005, 2007; Chen et al., 2008). When ketogenic diet was given to adenovirus-mediated FGF-21 knockdown mice, the mice had a fatty liver and hyperglycemia (Badman et al., 2007). In our study, plasma FGF-21 was positively correlated with BMI (data not shown), which is consistent with a previous study by Mashili et al. (2011). A list of obesity candidate genes, including FGF-21, was selected to draft “The Human Obesity Gene Map” (Rankinen et al., 2006). Here, we have sought to address whether 3′ UTR SNPs in the FGF-21 gene are associated with MetS, obesity, and diabetes in the Han Chinese population.
Our efforts focused on the 3′ UTR SNPs of FGF-21 because other 3′ UTR SNPs, such as PYY, have been associated with MetS traits (Shih et al., 2009). 3′ UTR SNPs impact on the regulation of miRNA by generating new targets for miRNA regulation, which could lower protein expression (Chen and Rajewsky, 2006; Landi et al., 2008); or altering miRNA targets, which could elevate protein expression (Clop et al., 2006; Adams et al., 2007). FGF-21 gene SNPs that had a high frequency (>0.1) of the minor allele in Asians were mainly concentrated in the 3′ region of the gene. In addition, it is reported that many novel loci located in the 3′ region, such as rs10811661 of CDKN2B, were associated with type 2 diabetes by genome-wide association analysis (Saxena et al., 2007). We identified 3 SNPs by screening the 3′ UTR of the FGF-21 gene, but only rs11665896 showed polymorphism in our study subjects.
FGF-21–transgenic mice were viable and resistant to diet-induced obesity (Kharitonenkov et al., 2005). Wang et al. (2008) identified an amino acid (F372) within helix 7 of the ligand-binding domain that is required for the response of PPARγ to endogenous ligands and also required for the expression of FGF-21, which shows that the FGF-21 gene was a direct target of PPARγ. Notably, rs11665896 was associated with higher risk of obesity in our study (OR: 1.84; 95% CI: 1.12, 3.01; p=0.016). Kaess et al. (2010) examined the associations of genetic variation in the FGF-21 signaling pathway with metabolic risk in subjects of white European ancestry. They found that FGF receptor 2 polymorphism (rs2071616) showed association with low-density lipoprotein cholesterol in both the CoLaus cohort (p=0.009) and men from the German Myocardial Infarction Family Study (p=0.017), but did not find association between FGF-21 gene polymorphism and metabolic risk (Kaess et al., 2010). We did not find association of rs11665896 with concentrations of plasma FGF-21, lipid concentration, and blood pressure. No study reported the relationship of rs11665896 with FGF-21 mRNA levels in liver or fat tissue. rs11665896 is about 100 bp from the 3′ end of the FGF-21 gene. rs11665896 is not predicted to affect mRNA stability and protein expression by bioinformaticsmethods (
Moreover, our study population was >50 years old. For the younger control subjects, the probability that they will develop a metabolic disease in the future is relatively high; therefore, the use of elderly subjects may be better. The sample size and the number of patients in each group were moderate in our study. In addition, the interaction between rs11665896 and the stability of FGF-21 mRNAs and protein expression levels remains to be evaluated.
In summary, FGF-21 rs11665896 was associated with a higher risk of obesity. We lack the power to assess the relationship between FGF-21 variants and MetS and diabetes in the Han Chinese population. However, further study on genetic targets may potentially improve understanding of metabolic disorders and lay a foundation for the treatment of MetS, obesity, and diabetes.
Footnotes
Acknowledgments
This study was supported by grants from the National Natural Science Foundation of China (#30900658), Sichuan University Young Scientist Funds (#2008095), and Sichuan Science and Technology Support Project (#2008SG0013).
Authors Contributions
Mei Zhang, Zhen-Mei An, and BinWu Ying participated in the research design. Mei Zhang, You-Juan Wang, and Li Zeng contributed to the performance of laboratory measurements. Mei Zhang and BinWu Ying were involved in data collection and analysis. Mei Zhang, Binwu Ying, and Zhen-Mei An wrote this article.
Disclosure Statement
The authors declare that they have no competing financial interests.
