INTRODUCTION
Previous study observed that the zinc finger homeobox 3 (ZFHX3) gene was isolated and reported (Morinaga et al., 1991), as to be encoding a protein which works as enhancer to human alpha-fetoprotein (AFP) gene. It was proposed that ZFHX3 binds to AT-rich motif in human AFP gene region, and it has 17 zinc finger motifs with large size of 306 kDa (Morinaga et al., 1991). AFP is considered as fetal gene which related to albumin, and its level rapidly declines by first 2 years of life (Schieving et al., 2014). Additionally, AFP gene is a famous marker of cancer in adult (Schieving et al., 2014; Staples, 1986).
However, Dong et al. (2010) proposed that ZFHX3 (which was shown as ATBF1) has an ATPase A-motif, two DEAH box-like sequences, and 23 zinc finger motifs with lengths of 3703 amino acids. They suggested that ZFHX3 may have a role in inhibition of estrogen receptor (ER)-mediated gene transcription and proliferation of breast cancer cells (Dong et al., 2010). Moreover, somatic gene mutation of ZFHX3 may also affect prostate cancer risk (Xu et al., 2006), and it was also involved in cell proliferation of prostate cancer (Sun et al., 2005). Breast cancer and prostate cancer are closely related diseases with obesity (Engin, 2017).
A recent study performed in Chinese population shows that level of serum AFP in metabolic syndrome were significantly higher (Chen et al., 2016), although the AFP level was not directly associated with central obesity alone in their result (Chen et al., 2016). However, increased serum AFP level was coexisted with fatty liver (Xu et al., 2014). Plus, there have been many studies presenting the relationship between ER and obesity, or adipocyte responses (Luglio, 2014; Miao et al., 2016; Pedram et al., 2016; Strong et al., 2015; Zhu et al., 2016). Therefore, it might be observed that those previous studies suggest the possible relationship, between the level of ZFHX3 and obesity.
In the genetic control of expressions, it has been suggested that intronic single nucleotide polymorphisms may be associated with protein levels, because there are parts of nucleotides which may serve as regulatory elements (Barrett et al., 2012), and which may also affect transcription as intron-mediated enhancement (Shaul, 2017). Cis-acting regulatory elements may affect expression levels of transcripts (Martin et al., 2014; Sadee, 2009). Furthermore, recent studies report novel mechanisms of noncoding RNA regulating cellular functions, such as circular RNAs (Hsiao et al., 2017) and heterochromatins (Shimada et al., 2016).
However, up to now, there was no previous study performing direct analysis of whether association between SNPs of ZFHX3 gene and obesity. Therefore, we investigated the SNPs of ZFHX3 gene in the groups of normal and obese individuals in Korean population.
MATERIALS AND METHODS
Study subjects
In the present study, 209 overweight/obese subjects and 159 controls were recruited. These subjects were recruited among participants that examined a general health check-up program. Subjects with severe diseases such as stroke, psychiatric disorders, and cancers were excluded. The biochemical characteristics of individuals were measured such as fasting plasma glucose, fasted glycated hemoglobin, total cholesterol. Body mass index (BMI) is calculated as weight (kg) divided by the square of height (m). According to the classification of Korean Society for the Study of Obesity (underweight, BMI<18 kg/m2; normal, BMI 18 to <23 kg/m2; moderately obese, BMI 23 to <25 kg/m2; obesity I, BMI 25 to <30 kg/m2; obesity II, BMI≥30 kg/m2), subjects were divided into two subgroups, the abnormal (overweight/obese) group (BMI ≥23 kg/m2) and the normal group (18 kg/m2≤BMI<23 kg/m2).
SNP selection and genotyping
Peripheral bloods of all subjects were collected in ethylenediamine tetraacetic acid or heparin tube. Genomic DNAs were extracted by QIAamp DNA mini kit (QIAGEN, Valencia, CA, USA). We selected 54 tagging SNPs in ZFHX3 gene. Genotype of each SNP was performed by custom DNA chip.
Statistical analysis
SNPStats (http://bioinfo.iconcologia.net/index.php) and IBM SPSS Statistics ver. 23.0 (IBM Co., Armonk, NY, USA) were used to determine the odds ratio (OR), 95% confidence interval (CI), and P-value. Logistic regression models (dominant [A/A genotype vs. A/B genotype+B/B genotype], recessive [A/A genotype+A/B genotype vs. B/B genotype], and log-additive models [A/A genotype vs. A/B genotype vs. B/B genotype]) were applied. The P-value below 0.05 was considered significant.
RESULTS
In order to evaluate the association between ZFHX3 gene and susceptibility of obesity, we genotyped and analyzed the 54 SNPs of ZFHX3 gene. In genotypic analysis, we observed that the seven SNPs (rs4788480, rs8055870, rs1010852, rs16971447, rs 9930445, rs4788489, and rs879324) were associated with susceptibility of obesity (P<0.05 in each SNP). In the control group, the genotype distributions for seven SNPs (rs4788480, rs8055870, rs1010852, rs16971447, rs9930445, rs4788489, and rs879324) were in Hardy–Weinberg equilibrium (rs4788480, P=0.058; rs8055870, P=0.51; rs1010852, P=0.16; rs16971447, P=1.00; rs9930445, P=0.34; rs4788489, P=0.10; and rs879324, P= 0.86, data not shown).
The genotype frequencies of the polymorphisms were compared between the control group and the overweight/obese group by using logistic regression model (dominant, recessive, and log-additive models). Among significant SNPs (rs4788480, rs8055870, rs1010852, rs16971447, rs9930445, rs4788489, and rs879324), minor allele of five SNPs (rs4788480, rs8055870, rs16971447, rs9930445, and rs879324) showed protective factor in susceptibility of obesity and minor alleles of two SNPs (rs1010852 and rs4788489) showed risk factor in present results.
The five SNPs showed association with obesity in each model (rs4788480: OR, 0.62; 95% CI, 0.41–0.95, P=0.260 in dominant model; rs8055870: OR, 0.64; 95% CI, 0.42–0.98; P= 0.0410 in dominant model; OR, 0.68; 95% CI, 0.47–0.99; P=0.440 in log-additive model; rs16971447: OR, 0.58; 95% CI, 0.38–0.89; P=0.0110 in dominant model; OR, 0.68; 95% CI, 0.50–0.93; P=0.0130 in log-additive model; rs9930445: OR, 0.52; 95% CI, 0.33–0.82; P=0.0041 in dominant model; OR, 0.72; 95% CI, 0.54–0.97; P=0.030 in log-additive model; and rs879324: OR, 0.70; 95% CI, 0.50–0.97; P=0.0320 in log-additive model).
And two SNPs (rs1010852 and rs4788489) also showed association with obesity (rs1010852: OR, 2.19; 95% CI, 1.16–4.13; P=0.012 in recessive model and rs4788489: OR, 1.93; 95% CI, 1.17–3.18; P=0.0083 in recessive model; OR, 1.49; 95% CI, 1.49; 95% CI, 1.11–1.99; P=0.0068 in log-additive model).
However, other SNPs did show any significant association with susceptibility of obesity. In haplotype analysis, we analyzed haplotype in significant SNPs (rs4788480, rs8055870, rs1010852, rs16971447, rs9930445, rs4788489, and rs879324) using Haploview 4.2. There were four haplotypes in one LD block consisted of rs16971447, rs9930445, and rs4788489 (TTT, CCG, TCG, and TTG haplotype) in ZFHX3 gene. Among four haplotypes, frequencies TTT haplotype and CCG in ZFHX3 gene between the control group and the overweight/obese group showed significant difference (TTT haplotype, P=0.0119 and CCG haplotype, P=0.0173).
DISCUSSION
ZFHX3 gene is located on chromosome 16, q arm position 22.2–22.3 (http://www.ncbi.nlm.nih.gov/gene/463). The chromosomal position is frequently associated with deletion in various cancer patients.
ZFHX3 was regarded as AFP enhancer, however, recent studies of ZFHX3 report its association to heart diseases. Interestingly, chromosome 16q22 was also associated with cardiovascular phenotypes. Several studies demonstrated that ZFHX3 SNP were significantly associated with ZFHX3 expression, which is in transcript-specific manner, regarding each of ZFHX3 isoforms, and they did not find variant associated to the expression (Hsiao et al., 2017). CAA repeat polymorphisms in exon 9 of the ZFHX3 gene was related to coronary heart disease in Chinese population, that 10 repeat allele was significantly associated with hypertension, diabetes, or dyslipidemia. There were also 9 or 11 repeats of CAA, however, rare alleles (Sun et al., 2015).
Zhai et al. (2015) found that rs7193343 may affect atrial fibrillation (AF) in Caucasians population, however it was not associated in Asian population (Zhai et al., 2015). In another study in Chinese population, rs2106261 and rs6499600 were significantly associated, and rs16971436 was borderline associated with AF (Liu et al., 2014). Liu et al. (2014) hypothetically suggested that ZFHX3 expression change may vary STAT3 regulation to result susceptibility to AF.
In addition, in the National Center for Biotechnology Information phenotype-genotype integrator (PheGenI) of ZFHX3 SNPs (https://www.ncbi.nlm.nih.gov/gap/phegeni?tab=1&gene=463), rs1548374 was associated with waist circumference, obesity, and BMI.
In our study result, seven SNPs of ZFHX3 gene (rs4788480, rs8055870, rs1010852, rs16971447, rs9930445, rs4788489, and rs879324) were associated with obesity. However, PheGenI obesity SNP rs1548374 was not associated in our result,
Our study has some limitations, that only a small Korean population was analyzed, whole variants near the ZFHX3 gene were not assessed, and actual expression in the individuals were not measured. Furthermore, ZFHX3 expression study will be needed in the adipocytes in the future studies. However, our study may partly agree with previous results that ZFHX3 variants may significantly affect phenotypic occurrence, and such disease phenotypes are in closely related spectrum to obesity.
In conclusion, we suggest that polymorphisms of ZFHX3 gene (rs4788480, rs8055870, rs1010852, rs16971447, rs9930445, rs4788489, and rs879324) may be contributed to susceptibility of obesity in Korean population.